Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: the interaction between patchy colloids is valence limited , specific , and directional due to the presence of patches or interaction sites on the surface of the colloids . patchy colloids are excellent candidates to design new materials with properties on - demand by controlling the microscopic details of the interparticle interaction . the distribution , number , and types of patches can be tuned  @xcite such that patchy colloids self - assemble into new structures that are not found in isotropically interacting colloids . examples are the directed self - assembly of patchy colloids into a kagome lattice  @xcite , polyhedra  @xcite , and micelles  @xcite . in addition to complex regular structures , patchy colloids also form stable gels . the number of patches is a key ingredient determining the location of the liquid - gas critical point  @xcite . in systems of patchy colloids with low valence the liquid - gas phase separation region shrinks drastically , and it is possible to find stable liquids at very low densities and temperatures  @xcite . these , equilibrium , low density states are gels in which the patchy colloids assemble into amorphous percolated networks . mixtures of patchy particles provide an additional handle for tuning the morphology of the gel . in  @xcite we predicted theoretically the existence of new gel structures in binary mixtures with intriguing percolation properties , such a bicontinuous gel or bigel . here , two interpenetrating networks , each made of one type of colloids , span the system volume . other gel structures are possible , such as mixed gels where the mixture is percolated but neither species percolates independently  @xcite .    bigel - like percolated states are not specific of patchy colloidal systems . bicontinuous structures have been reported in e.g. , polymer blends  @xcite , mixtures of dipolar colloidal particles  @xcite , fumed silica - based systems  @xcite , mixtures of an aqueous gel and an oleogel  @xcite , dna coated colloids  @xcite , and mixtures of proteins  @xcite .    in this work , we investigate the structure of equilibrium gels in binary and ternary mixtures of patchy colloids . we test the validity of our extension of the flory - stockmayer theory of polymerization  @xcite to multicomponent systems with different types of patches  @xcite by comparing theoretical with monte carlo simulation results . the latter agree semi - quantitatively with the theoretical predictions and corroborate our previous findings  @xcite of bigels and mixed gels in binary systems of patchy colloids . additionally , we extend the analysis of the percolation of patchy colloids to ternary mixtures both theoretically and using computer simulations . we find eight fundamentally different types of percolated states that we classify according to their strength against breaking bonds . the strongest gel is a trigel or tricontinuous gel in which there are three interpenetrated networks , each one made of one type of colloids . the rich phenomenology of the ternary system can be understood as a competition between the entropy of mixing and the internal energy of bonding . at high temperature the entropy of mixing , which favors interspecies bonds , dominates and promotes the formation of different types of mixed gels . at low temperatures , however , the internal energy of bonding is the most important contribution to the free energy . here , we observe the occurrence of strong gels such as bigels and trigels provided that the intraspecies bonds are stronger than the interspecies ones . the stability of each type of gel can be controlled by changing the temperature and composition of the mixture . the model patchy colloidal particles , see fig . [ fig1 ] , consists of a hard core ( hard sphere of diameter @xmath0 ) with three interaction sites or patches on the surface . the patches are spheres of size @xmath1 centered at the surface of the hard core , uniformly distributed along the equator of the particles , i.e. the angle between two patches is @xmath2 . we set @xmath3 , which is the largest value that guarantees single bonding of each patch . in addition , the geometry of the particles ensures that two colloids can not form more than one bond between them . the interaction between patches @xmath4 is a square well potential , i.e. if two patches of type @xmath5 and @xmath6 overlap , the energy of the system decreases by a quantity @xmath7 : @xmath8 where @xmath9 indicates the position and orientations of colloids @xmath10 and @xmath11 . we consider binary and ternary mixtures . the only difference between the species is the type of their patches . in the case of binary mixtures the particles of species @xmath10 @xmath12 have patches of type a ( b ) . for ternary mixtures we add a third species with patches of type c.     with three bonding sites of size @xmath1 on the surface . the only difference between species is the type of their bonding sites : type a for species @xmath10 , type b for species @xmath11 , and type c for species @xmath13 ( not shown ) . a bond between patches of type @xmath5 and @xmath6 decreases the energy of the system by a quantity @xmath7 . ] we use a generalization of the flory - stockmayer theory of polymerization  @xcite to study the percolation transition . the theory neglects bond correlations , and therefore assumes that all clusters are tree - like . for details of the theory for multicomponent systems with different types of patches we refer the reader to refs . @xcite . here , we describe the theory for our patchy particle model : each particle has three identical patches , and there are three types of particles . the binary system may be obtained as a particular case . a schematic tree - like cluster is shown in fig . let @xmath14 be the number of patches of type @xmath5 that are bonded at level @xmath15 of the tree - like cluster , see fig . @xmath14 is related to the number of bonds that are bonded at the previous level @xmath16 via @xmath17 where @xmath18 is the probability of bonding a site @xmath5 to a site @xmath6 . note that , in general , @xmath19 since @xmath18 is proportional to the number of available patches of type @xmath6 and hence to the number of particles with @xmath6 patches , as we will see . the factor of two in all the terms of the above equation reflects the fact that one patch is bonded at level @xmath16 by construction , and therefore only two patches are available to form bonds at the next level @xmath15 . we can rewrite eq .   in matrix form @xmath20 where the components of the vector @xmath21 are @xmath22 , and the percolation matrix @xmath23 is @xmath24 the system percolates if the number of bonds diverges in the limit @xmath25 , or equivalently if @xmath26 with @xmath27 the largest eigenvalue of the percolation matrix , see @xcite for more details . in order to establish if subsets of the whole system are also percolated we have to consider a percolation matrix that includes only the desired bonds . for example , the pair formed by species @xmath10 ( patches of type a ) and species 2 ( patches of type b ) has percolated if the largest eigenvalue of the matrix : @xmath28 is larger than one . the same procedure applied to a single species leads to the known percolation threshold for tree - like clusters of identical particles with only one type of patch : @xmath29 , with @xmath30 the number of patches . is the probability of bonding a site @xmath5 to a site @xmath6 , with @xmath31 . ] we obtain the probabilities @xmath18 from wertheim s first order perturbation theory  @xcite , which also neglects correlations between bonds . the probability of bonding for a site @xmath5 is @xmath32 where @xmath33 is the probability that a site @xmath5 is not bonded . the probabilities @xmath33 are obtained from wertheim s theory . for our model , we find ( see @xcite for details ) @xmath34^{-1 } , \label{xa}\ ] ] where @xmath35 is the composition of species @xmath16 , the subscript @xmath6 denotes the type of patches of species @xmath16 ( i.e. , @xmath36 if @xmath37 , @xmath38 if @xmath39 , and @xmath40 if @xmath41 ) , @xmath42 is the total packing fraction of the system , and @xmath43d(12).\ ] ] here , @xmath44 is the volume of the hard - core of the patchy colloids , @xmath45 is the boltzmann constant , and @xmath46 is the radial distribution function for hard - spheres at a distance @xmath47 . the integral is calculated over all possible orientations and separations of the colloids . finally , we approximate @xmath48 by its contact value @xmath49 ( the patches are very small compared to the size of the colloids ) . hence , @xmath50,\ ] ] with @xmath51 and @xmath52 is the bonding volume , i.e. , the volume in which two colloids can form a bond averaged over all possible orientations of the colloids . for our model of the patches  @xcite @xmath53 .    analyzing term - by - term eqs . and   we obtain @xmath54 which has an intuitive explanation . the probability of bonding a site @xmath5 to a site @xmath6 , @xmath18 , is proportional to : ( i ) the density of available patches of type @xmath6 which is @xmath55 , and ( ii ) the probability that a patch @xmath5 is available ( not bonded ) which is @xmath33 . the parameter @xmath56 controls the strength of the bond . we have also obtained the probabilities directly from monte carlo simulations and used them to compute the percolation matrix . we compare both methods in the results section . we use monte carlo simulations in the canonical ensemble , that is , we fix the total number of particles of each species @xmath57 , the temperature @xmath58 , and the system volume @xmath59 . the simulation box is a cube with periodic boundary conditions . the total number of particles in the simulation box is @xmath60 and the packing fraction @xmath61 . we have also computed selected cases with @xmath62 and the same packing fraction and found no significant differences . therefore , we expect that finite size effects are irrelevant .    in each simulation we collect data from @xmath63 monte carlo steps ( mcs ) . each mcs is an attempt to move and rotate all particles in the system . in addition to translations and rotations of the particles we introduced a particle swap move  @xcite . the swap move randomly selects two particles and interchanges their coordinates ( spatial and orientational ) . the swap move is accepted according to the standard metropolis criterion and satisfies detailed balance . we perform @xmath64 swap attempts after every mcs . the particle swap is very useful to prevent trapping of the particles in local minima of the energy landscape , especially at low temperatures . before collecting the data , we equilibrate the system by running mcs until the energy fluctuates around a minimum value . typically @xmath65 mcs are used for equilibration , depending on the temperature .    for each composition we start by running a simulation at high temperature . next we decrease @xmath58 using the last configuration of the previous simulation as the initial state . we repeat this process for several @xmath58s . in order to check we have reached the equilibrium states we carried out the inverse process for selected cases . that is , we started at low @xmath58 and increased the temperature . we have found the same states both by increasing or decreasing @xmath58 . the main output of the simulation is the cluster size distribution , @xmath66 where @xmath67 is the number of clusters of size @xmath68 . we compute @xmath69 for the whole mixture , for each species @xmath70 , and in the case of ternary mixtures for pairs of species @xmath71 . using the cluster size distribution we can compute the probability of finding a particle in a cluster of size @xmath68 : @xmath72 and the moments of @xmath69 . we use the second moment of the cluster size distribution @xmath73 to characterize the percolation transition . we set the percolation threshold at @xmath74 , i.e. , we consider the whole system has percolated if @xmath75 . in order to characterize the type of percolated state we also analyze the second moments of the cluster size distributions of single species and pairs of species . e.g. , the species @xmath16 has percolated if @xmath76 . we start by describing the results for a binary mixture , a system that we have previously studied using wertheim s theory  @xcite . all possible gels are listed in fig . we order the gels according to their relative strength against breaking bonds . the weakest gel is a mixed gel ( mg ) . in a mg the mixture has percolated , but neither species @xmath10 nor species @xmath11 have percolated . the removal of one of the species will break the connectivity of the system . next we have a standard gel . here , one species percolates , and as a consequence the mixture also percolates ( in the mixture we add new bonds to a system that has percolated ) . the percolation of the mixture is driven ( induced ) by the percolation of one of the species . a standard gel is stronger than a mixed gel in the sense that we can remove one species , the one that has not percolated , and the system remains in a percolated state . the strongest gel is a bicontinuous ( bg ) gel , where both species percolate independently . we can remove one species , either @xmath10 or @xmath11 , and the other species remains in a percolated state . the relative stability of the different gels depends on many factors , such as temperature , composition , packing fraction and ratio between the different bonding energies . in fig . [ fig4 ] we show an example of a mixed and a bicontinuous gel . both mixtures are symmetric ( @xmath77 ) , and the packing fraction ( @xmath78 ) , and composition ( @xmath79 ) are also the same . a mixed gel ( see a snapshot in fig .  [ fig4]a1 ) can be easily stabilized by making the interspecies bonds stronger than the intraspecies bonds . in this way both the entropy of mixing and the internal energy of bonding favour the formation of interspecies bonds , which are dominant in a mixed gel ( see fig . [ fig4]a2 ) . the probabilities of finding a particle in a cluster of @xmath68 identical particles @xmath80 ( fig . [ fig4]a3 ) show that only small clusters of identical particles are present , whereas the probability of finding a particle in a cluster of @xmath68 particles ( independently of the particle type ) clearly indicates that the mixture has percolated . if interspecies bonds are weaker than intraspecies bonds , i.e. @xmath81 with @xmath82 , then the energy of bonding may overcome the entropy of mixing . in this case a bicontinuous gel appears ( see a snapshot in fig .  [ fig4]b1 ) . here , intraspecies bonds dominate ( see fig . [ fig4]b2 ) , and the cluster size distribution indicates that both species independently and the mixture have percolated ( see fig . [ fig4]b3 ) . the ratio between inter- and intra - species bonding energies plays a major role in determining the state of the sample but it is not the only relevant variable . temperature , packing fraction , and composition are also important . we illustrate this in fig . [ fig5 ] , where the second moment of the cluster size distributions as a function of the temperature is plotted for a symmetric mixture with @xmath83 and @xmath84 . at high temperatures there are only a few bonds and the system has not percolated . percolation occurs at @xmath85 . the system is always percolated below @xmath86 but the type of percolation changes . in a range of temperatures below @xmath86 there is a mixed gel ; only the two species together percolate . however , at low temperatures the state changes to a gel of species @xmath11 . that is , the particles of species @xmath11 , which are the majority , form a percolated network . the whole mixture also percolates as addition of another species only adds new bonds to the network .     and for the mixture @xmath87 as a function of the scaled temperature @xmath88 . symmetric mixture ( @xmath89 ) with @xmath90 , composition @xmath91 , and packing fraction @xmath61 . ] the full percolation behaviour is summarized in the composition - temperature percolation diagrams depicted in fig . the percolation diagrams have been obtained via monte carlo simulations analyzing the cluster size distribution ( left ) and theoretically , assuming that the network is a tree - like cluster and using the bonding probabilities measured in the simulations ( symbols ) and calculated using wertheim s theory ( lines ) . we first describe the global common features and then compare the two approaches .    * symmetric mixture with @xmath92*. in fig . [ fig6]a we show the results for a symmetric mixture with @xmath93 and @xmath94 . three percolated states are possible below the percolation threshold : two standard gels and a mixed gel . the mixed gel is stable in a wide range of compositions , and it is the first type of gel stabilized upon cooling the system . this is the expected behaviour as the whole mixture will always percolate before the individual species percolate . the temperature at the percolation threshold has a maximum for an equimolar mixture , @xmath95 , which is the composition where the number of possible interspecies bonds , with the highest bonding energy , is largest . at compositions @xmath96 ( @xmath97 ) we also find a gel of species @xmath11 ( @xmath10 ) at sufficiently low temperature . intraspecies bonds stabilize a standard gel , and they increase as the composition approaches that of a single component system . therefore , the temperature where the standard gel is stabilized increases as we approach the limit of pure systems @xmath98 and @xmath99 . in other words , the range in compositions where the mixed gel is stable decreases with the temperature ( in the range of temperatures analyzed here ) .     and scaled temperature @xmath100 for different bonding energies . @xmath101 is the composition of species @xmath10 . ( a ) @xmath93 , @xmath94 , ( b ) @xmath93 , @xmath102 , and ( c ) @xmath103 , @xmath104 . left column : results from the cluster size distribution obtained using monte carlo simulations . right column : results from the analysis of the percolation matrix with probabilities extracted from the monte carlo data ( symbols ) and computed using wertheim s theory ( lines ) . the mixture percolates at temperatures below the green line . species 1 ( 2 ) has percolated if @xmath58 is lower than the red ( blue ) line . ]    * symmetric mixture with @xmath105*. if the intraspecies bonds have the highest energy there is a competition between the entropy of mixing ( which favours a mixed gel ) and the internal energy of bonding ( which favours a bicontinuous gel ) . an example is shown in fig . [ fig6]b . at low temperatures the internal energy dominates and we observe a bigel ( bg ) in a range of compositions around @xmath95 . the lower the temperature the wider the composition range where a bg is stable . at high temperatures the entropy of mixing overcomes the internal energy of bonding and a mixed gel is stabilized . by contrast to the case with @xmath106 , the percolation temperature of the mixture exhibits a minimum at the equimolar composition @xmath95 rather than a maximum . this is simply because the number of interspecies bonds is maximal at @xmath95 , which are now the bonds with the lowest energy , and thus a lower temperature is required to stabilize the bonds . standard gels are also present in this mixture near the composition of single component systems .    * asymmetric mixture*. the diagrams discussed above are not present only in the very special case of symmetric mixtures . an example is shown in fig . [ fig6]c where we consider an asymmetric mixture with bonding energies @xmath107 and @xmath104 . the diagram is now asymmetric with respect to the equimolar composition . gel @xmath10 dominates over a gel of species 2 since the @xmath108 bonds are stronger than @xmath109 bonds . although @xmath110 there is a bigel at low temperatures . the reason is that the bonding energies satisfy @xmath111 , i.e. , it is energetically unfavorable to form two @xmath112 bonds rather than one @xmath108 and one @xmath109 . therefore a bicontinuous gel should be stable at sufficiently low temperature . we next compare the percolation diagrams obtained directly using the cluster size distribution ( left column of fig . [ fig6 ] ) and using the percolation matrix and the bonding probabilities ( right column of fig . [ fig6 ] ) . the derivation of the percolation matrix neglects the correlation between bonds and in particular neglects closed loops . the clusters are tree - like and their size is generally larger than that obtained from the simulations , where bond correlations and closed loops are present . consequently , the threshold temperatures predicted by the percolation matrix are higher than those obtained from the simulations . this is , however , the only significant difference . in all cases , the results are qualitatively the same . we note that the bonding probabilities were obtained in two different ways : ( i ) theoretically using wertheim s theory ( solid lines in fig . [ fig6 ] ) and ( ii ) directly from the simulations ( symbols in fig .  [ fig6 ] ) by measuring the number of bonds of each type . the resulting percolation diagrams are in perfect agreement and therefore we are confident that closed loops are the only relevant missing ingredient in our theoretical description . we emphasize that we have analyzed the percolation diagram only . some of the states considered above may be metastable against phase separation . for example , for the mixture with @xmath113 we observed signs of demixing at very low temperatures . we will return to this point in the discussion .      in this section , we consider ternary mixtures by adding a third component ( species @xmath13 ) with interaction sites of type @xmath114 . the addition of a new species increases significantly the number of percolated states . in fig . [ fig7 ] we depict all gels labelled according to their strength against breaking bonds :        * mixed gels . * the weakest type of percolated state is again a mixed gel . in a ternary mixture there are four different mixed gels @xmath115 with @xmath116 indicating the strength ( the lower the index the weaker the gel ) . in mixed gels the single species do not percolate . the difference between these gels relates to different percolating states of pairs of species . the weakest mixed gel is @xmath117 where only the ternary mixture percolates . a subset of species , e.g. , a mixture of species @xmath10 and @xmath11 does not percolate . next , we find gels @xmath118 where there is one and only one pair of species that percolates ( e.g. , species @xmath10 and @xmath11 , which in the following we denote by @xmath119 ) , and thus the ternary mixture also percolates . there are three @xmath118 gels related to the pair of species that percolates : @xmath119 , @xmath120 , and @xmath121 . the bonds of the species that is not involved in percolation may be removed without breaking the global connectivity of the system .    then we find @xmath122 gels where two pairs of species percolate . again , there are three different gels ( see fig . [ fig7 ] ) . we can remove only one species without breaking the global connectivity , as in @xmath118 , but now we have two choices . for example , if the pairs @xmath119 and @xmath120 are percolated , we can break the bonds of either species @xmath11 or @xmath13 without affecting the global connectivity . the strongest mixed gel is @xmath123 where any pair of species has percolated . as in the previous cases we can break the bonds of only one species within the percolated state , but we can choose any of the three species . * standard gels . * here , one and only one of the species has percolated . any pair of species where one component is the percolated species is also percolated , and the ternary mixture has percolated as well . the pair of species that does not involve the percolated component has not percolated . in a standard gel we can break the bonds of two species ( the non percolated components ) without affecting the connectivity of the remaining particles . * standard plus mixed gels . * in this case the @xmath16-th species ( @xmath124 or @xmath13 ) has percolated and the pair of species @xmath125 with @xmath126 has also percolated . as for a standard gel the connectivity is not lost if we break the bonds of species @xmath127 and @xmath128 . in addition , we have the possibility of breaking the bonds of species @xmath16 .    * bicontinuous gels . * as for binary mixtures bicontinuous gels may also occur . two species @xmath16 and @xmath127 percolate and any pair of species also percolates . global connectivity is not affected by breaking the bonds of two species provided that one is species @xmath129 . * trigel . * the trigel is the strongest gel . the three species percolate and all pairs of species also percolate . the bonds of any two species may be broken without affecting the global connectivity of the system .    in binary mixtures when the intraspecies bonds are stronger than interspecies ones the percolation diagram is very rich and exhibits all types of gels . the reason is a competition between entropy and internal energy . the entropy of mixing promotes interspecies bonds which stabilize the mixed gel . by contrast , the internal energy of bonding favour intraspecies bonds which are dominant in the bicontinuous gel . at high ( low ) temperature the entropy ( internal energy ) dominates and a mixed ( bicontinuous ) gel is formed . the same competition occurs in ternary mixtures . the richest percolation behaviour also arises when interspecies bonds are weaker than intraspecies ones . an example is shown in fig . [ fig8 ] where we plot the full percolation diagrams at different temperatures for a system with bonding energies @xmath130 ( intraspecies ) and @xmath131 ( interspecies ) with @xmath132 and @xmath133 . given the excellent agreement between theory and simulations for binary mixtures we have studied the percolation diagrams theoretically , and simulated only selected systems for each type of gel . decreasing the temperature leads to an increase in the number of bonds , strengthening the percolated phases . as we will see , decreasing the temperature decreases the regions of the percolation diagram occupied by the weakest gels and increases the regions occupied by the strongest ones . occasionally , two or more regions merge giving rise to a new stronger gel . this observation describes the most relevant features of the percolation behaviour of this type of ternary mixtures . the highest temperature considered is @xmath134 , see fig . although this is slightly below the percolation temperature of a single component system , the percolation diagram is dominated by a non percolated state . only at composition close to the pure systems does the ternary mixture percolate . in these regions the percolation behaviour is very rich ( see the close up view in fig .  [ fig8]a ) . in the limit @xmath135 a gel of species @xmath16 occurs . by moving away from the pure species we find the sequence of gels @xmath136 and finally a non percolated state . to understand this sequence consider the following . at @xmath137 the pure species has percolated and thus also the subsets @xmath138 and @xmath139 . as we decrease @xmath140 we replace particles of species @xmath16 by particles of species @xmath127 and @xmath128 and therefore reduce the number of intraspecies bonds @xmath16 . since the temperature is high ( slightly below the percolation temperature ) , the decrease of intraspecies bonds @xmath16 disrupts the global connectivity of this species . thus , the next gel state is either a @xmath122 where the subsets @xmath138 and @xmath139 percolate , or a @xmath118 where only one of the previous two subsets percolates . the compositions @xmath141 and @xmath142 determine which of the gels , @xmath122 or @xmath118 , is stable . further decreasing the number of intraspecies bonds @xmath16 ( i.e. moving away from the vertices of the percolation diagram ) the subsets @xmath138 or @xmath139 become non percolating and the state changes to a @xmath117 where only the ternary mixture percolates @xmath143 . finally , we find a transition between a @xmath117 and a non percolated state at compositions sufficiently far from the pure system . the reason is that approaching the center of the percolation diagrams increases the number of interspecies bonds , which have the lowest energy , effectively decreasing the percolation temperature .        the range of compositions where the mixture percolates increases by further decreasing the temperature . at @xmath144 see fig . [ fig8]b , only a small region close to the equimolar composition remains non percolated . there are no new percolated states with the previous ones occupying a larger area of the percolation diagram , especially the @xmath117 gel that has grown at the expense of the non percolated state .    at @xmath145 , see fig . [ fig8]c , the mixture has percolated over the whole range of composition . the standard gels , @xmath122 , and @xmath118 continued to grow ( occupying larger regions of the percolation diagram ) at the expense of the @xmath117 gel , which is the weakest and has been reduced to a small region around the equimolar composition .    at @xmath146 , see fig . [ fig8]d , the @xmath117 gel has disappeared . the three regions where the @xmath122 gel is stable have merged in the center of the percolation diagram , giving rise to a new state , the @xmath123 gel , where any pair of species percolates . the standard gels , which are the strongest at this temperature , continue to grow at the expense of the weakest gel , @xmath118 .    at the next temperature @xmath147 ( see fig .  [ fig8]e ) , the @xmath118 has disappeared and new states occur . two standard gels merge at compositions close to @xmath148 giving rise to bicontinuous gels of species @xmath16 and @xmath127 . next to the regions of stability of the bigels , a standard gel and a @xmath123 gel have merged and a gel+@xmath123 phase has appeared . the weakest gel is now @xmath122 and the strongest is no longer a standard gel but a bigel . bigels and gel+@xmath123 states grow by further decreasing @xmath58 , see fig . [ fig8]f , and completely replace the weak @xmath122 gel state . at this point , the central @xmath123 gel is the weakest , and starts to shrink . at @xmath149 , see fig . [ fig8]g , the @xmath123 is confined to a small area near the equimolar composition . the standard gels also shrink as the stronger gels grow . the trigel , the strongest of all gels , appears at the next temperature @xmath150 as a result of the merging of the three bigels close to the equimolar composition , see fig . if we continue to cool the system the trigel grows at the expense of the weaker gels , see fig . . no further topological change of the percolation diagram occurs at lower temperatures . in addition to the theoretical percolation diagrams , we have simulated selected systems at different compositions and temperatures . in all cases we have found perfect agreement of the sequences of percolated states as the temperature decreases . for example , by cooling an equimolar system ( @xmath151 ) the simulations exhibit the sequence non percolated-@xmath117- @xmath123-trigel , as predicted by the theory . as an example we show in fig . [ fig9 ] the simulation results for the gels @xmath117 ( top panels ) and @xmath123 ( bottom panels ) , obtained by decreasing the temperature from a non percolated state . the only relevant difference was the temperature at which each phase appears . we found , for example , that a trigel is stable at temperatures @xmath152 , whereas theoretically we observe stable trigels at @xmath153 . as for binary mixtures , the theory overestimates the percolation temperature of each type of gel . this discrepancy is attributed to the neglect of closed loops in the theoretical description . we have also computed the percolation diagrams for other sets of bonding energies , and the sequence of states and topology of the percolation diagrams is always similar ( the symmetry is obviously lost in cases with asymmetric interaction energies ) . the most significant differences occur if all interspecies bonds are energetically more favorable than the intraspecies bonds . in this case , the trigel and the bigels are not stable since they require a predominance of interspecies bonds . also , the percolation temperature of the weakest @xmath117 gel has a maximum at the equimolar composition , instead of a minimum ( like in fig . [ fig8]a - b ) . this is analogous to what happens in binary mixtures . we have considered the percolation transition in binary and ternary mixtures of patchy colloidal particles , at fixed packing fraction ( @xmath61 ) . the results of wertheim s theory and monte carlo simulations indicate that a pure system of patchy particles with three identical patches at this packing fraction is thermodynamically stable over the whole range of temperature  @xcite , i.e. , there is no liquid - gas or liquid - solid phase separation at @xmath61 . the addition of other species might change the stability of the fluid phase . for example , we expect demixing at sufficiently low temperatures if the intraspecies bonds are stronger than the interspecies ones , in line with wertheim s theory for binary mixtures  @xcite . nevertheless , in  @xcite we have shown that it is possible to design more complex patchy colloids ( with different types of patches ) in order to suppress the phase separation and stabilize the desired gels . in addition , arrested phase separation might occur , and effectively stabilize the gel . an example is the formation of bigels in mixtures of dna - coated colloids as a results of arrested demixing  @xcite . our theoretical description neglects the effect of closed loops . nevertheless , simulations ( closed loops included ) and theory ( no closed loops ) are in semi - quantitative agreement . the effect of loops in the present model simply shifts the percolation threshold to lower temperatures . the bonding probabilities are almost unaffected by the presence of closed loops . this is not a general feature , as closed loops could have a major impact in the percolation behaviour of the system . this will occur in systems where the formation of ring - like structures is favoured . for example , in a system of particles with two patches of type a near the poles and several patches of type b on the equator  @xcite the formation of rings through aa bonds leads to a lower critical point in the liquid - vapor binodal . see  @xcite for recent discussions of extensions of wertheim s theory to include the effect of loops . we have characterized the percolated states or gels according to the species that are percolated . for example , in a bigel there are two species that percolate independently . one can also characterize the states according to the types of bonds that are percolated . in this case one can for example distinguish between two types of bigels , one bigel where there are only two percolated networks formed by intraspecies bonds and another bigel where the interspecies bonds also form a percolated network . our model may help to understand network formation in other systems such as abc block terpolymers  @xcite , and ternary polymer blends  @xcite . in these systems one can find percolated states analogous to those described here . for example , the percolation behaviour of a ternary polymer blend made of high - density polyethylene , polystyrene , and pmma is topologically similar to the low temperature percolation behaviour of the ternary mixture analyzed here  @xcite ( fig . [ fig8]h ) . future studies will focus on the percolation behaviour under gravity . a gravitational field is , in general , unavoidable in sedimentation experiments , and might have a profound effect on colloidal mixtures  @xcite . for example , in a mixture of patchy colloids with two and three patches ( all of the same type ) gravity induces the formation of complex sequences in sedimentation - diffusion - equilibrium such as the presence of re - entrant percolation  @xcite . we thank jose m. tavares for useful discussions . this work was support in part by the portuguese foundation for science and technology ( fct ) through projects excl / fis - nan/0083/2012 and uid / fis/00618/2013 . 35ifxundefined [ 1 ] ifx#1 ifnum [ 1 ] # 1firstoftwo secondoftwo ifx [ 1 ] # 1firstoftwo secondoftwo `` `` # 1''''@noop [ 0]secondoftwosanitize@url [ 0 ]  + 12$12  & 12#1212_12%12@startlink[1]@endlink[0]@bib@innerbibempty @noop * * ,   ( ) link:\doibase 10.1038/nature11564 [ * * , ( ) ] http://stacks.iop.org/0953-8984/25/i=19/a=193101 [ * * ,   ( ) ] link:\doibase 10.1021/acsnano.5b07901 [ * * ,   ( ) ] link:\doibase 10.1038/nature09713 [ * * ,   ( ) ] link:\doibase 10.1021/nl0493500 [ * * , ( ) ] link:\doibase 10.1073/pnas.1116820109 [ * * ,   ( ) ] http://scitation.aip.org/content/aip/journal/jcp/128/14/10.1063/1.2888997 [ * * ,   ( ) ] link:\doibase    10.1103/physrevlett.97.168301 [ * * ,   ( ) ] link:\doibase 10.1039/c0sm01493a [ * * , ( ) ] link:\doibase 10.1039/c1sm06948a [ * * , ( ) ] link:\doibase 10.1021/ma010104 + [ * * ,   ( ) ] link:\doibase 10.1002/pen.760281205 [ * * ,   ( ) ] http://scitation.aip.org/content/aip/journal/jcp/133/6/10.1063/1.3477985 [ * * ,   ( ) ] link:\doibase 10.1039/b907873h [ * * ,   ( ) ] link:\doibase 10.1039/c4ra15437a [ * * ,   ( ) ] \doibase    http://dx.doi.org/10.1016/j.colsurfb.2014.09.056 [ * * ,   ( ) ] \doibase    http://dx.doi.org/10.1016/j.msec.2014.08.026 [ * * ,   ( ) ] link:\doibase    10.1073/pnas.1214971109 [ * * ,   ( ) ] link:\doibase    10.1039/c3sm52558a [ * * ,   ( ) ] link:\doibase 10.1039/c5tb00131e [ * * , ( ) ] link:\doibase 0.1021/ja01856a061 [ * * , ( ) ] link:\doibase 10.1063/1.1723803 [ * * , ( ) ] link:\doibase 10.1103/physreve.81.010501 [ * * ,   ( ) ] link:\doibase 10.1007/bf01127721 [ * * , ( ) ] \doibase http://dx.doi.org/10.1063/1.2730797 [ * * , ( ) ] link:\doibase 10.1103/physreve.63.045102 [ * * ,   ( ) ] link:\doibase 10.1103/physrevlett.111.168302 [ * * ,   ( ) ] link:\doibase 10.1103/physreve.87.052307 [ * * ,   ( ) ] link:\doibase 10.1039/c5sm00559k [ * * , ( ) ] link:\doibase 10.1021/ma9009593 [ * * , ( ) ] link:\doibase 10.1021/ma500302k [ * * ,   ( ) ] http://www.sciencedirect.com/science/article/pii/s0032386110006154 [ * * ,   ( ) ] link:\doibase 10.1039/c3sm51491a [ * * , ( ) ] link:\doibase 10.1103/physreve.93.030601 [ * * ,   ( ) ]
Abstract:

Answer: we investigate percolation in binary and ternary mixtures of patchy colloidal particles theoretically and using monte carlo simulations . each particle has three identical patches , with distinct species having different types of patch . theoretically we assume tree - like clusters and calculate the bonding probabilities using wertheim s first - order perturbation theory for association . for ternary mixtures we find up to eight fundamentally different percolated states . the states differ in terms of the species and pairs of species that have percolated . the strongest gel is a trigel or tricontinuous gel , in which each of the three species has percolated . the weakest gel is a mixed gel in which all of the particles have percolated , but none of the species percolates by itself . the competition between entropy of mixing and internal energy of bonding determines the stability of each state . theoretical and simulation results are in very good agreement . the only significant difference is the temperature at the percolation threshold , which is overestimated by the theory due to the absence of closed loops in the theoretical description .


Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: interactions between massive stars and the interstellar medium ( ism ) are among the dominant mechanisms involved in the chemical and dynamical evolution of galaxies . massive stars deposit thermal and kinetic energy into the ism through their energetic winds and eventual supernova explosions . one incarnation of these interactions is the `` superbubble , '' a large ( 100 - 200 pc diameter ) , roughly spherical interstellar shell around one or more ob associations . a superbubble shell consists of interstellar material swept up by the stellar winds and supernova ejecta . the shell interior contains hot ( 10@xmath2 k ) , shocked wind and ejecta with density of @xmath3 @xmath4 surrounded by the cooler , denser gas of the shell ( castor , mccray , & weaver 1975 ; weaver et al . 1977 ) .    a large number of superbubbles have been identified in local group galaxies . the large magellanic cloud ( lmc ) is a nearly ideal laboratory in which to study the detailed physics of these interstellar structures , as the lmc has a low foreground interstellar extinction and a small inclination angle of 40@xmath5  ( feast 1991 ) , minimizing confusion along the line - of - sight . at the lmc distance of 50 kpc , 1  corresponds to 15 pc , allowing us to study the gas portions of these structures in great detail with a variety of instruments , as well as to resolve the bulk of the individual massive stars with ground - based telescopes . dozens of superbubbles are known in the lmc , giving many examples for study . n44 is among the most luminous  complexes in the lmc , and consists of three ob associations and several regions of optical nebulosity ( henize 1956 ; lucke & hodge 1970 ) . a variety of -bright structures are seen in n44 , and identifications of the components , both bright knots and large filamentary arcs , have been made by henize ( 1956 ) with labels n44 a  n and davies , elliott , & meaburn ( 1976 ) with labels dem l140  170 . figure [ halpha]a shows an  ccd image of the region with the major emission regions marked . meaburn & laspias ( 1991 ) detected the expansion of the two of the large shells , dem l140 and dem l152 , confirming their description as superbubbles . diffuse x - ray emission from various portions of the n44 region was detected by the _ einstein observatory _ ( chu & mac low 1990 ; wang & helfand 1991 ) . observations of the region with rosat ( chu et al . 1993 , hereafter c93 ) allowed for a more detailed comparison between the x - ray and h@xmath6 morphologies , and showed several distinct features . c93 used the x - ray morphology to define five regions of particular interest : shells 1  3 , the south bar , and the north diffuse region . among the three brightest regions , shell 1 was identified as an x - ray - bright superbubble , the south bar was suggested to be a `` blowout structure '' from the superbubble , and shell 3 was identified as a previously unknown supernova remnant ( snr ) . these three bright x - ray features are marked in figure 2a , a rosat position sensitive proportional counter ( pspc ) image overlaid with h@xmath6 contours . the north diffuse region is outside of the bounds of this picture and shell 2 is not shown and will not be discussed further in this paper . shell 1 corresponds primarily with the interior of the optical shell dem l152 . the stellar content of dem l152 has recently been discussed by oey & massey ( 1995  hereafter om95 ) , who make model comparisons between the observed expansion velocities and the expected energy input from the stellar winds . two important results are drawn from their study . first , the stellar content interior to the shell is generally older than the population immediately exterior to the shell ( @xmath7 myr vs.  @xmath85 myr ) , suggestive of sequential star formation . second , their superbubble evolution models , adopting stellar wind and supernova energy inputs implied by realistic stellar contents , over - predict the shell diameter . several possible reasons for this size discrepancy are given by om95 ; however , we believe that part of the discrepancy may be caused by significant energy / pressure lost via the blowout shown in c93 .    in this paper , we have combined newly - obtained  ccd images , high - resolution echelle spectroscopy and asca x - ray observations , along with the existing rosat observations to test the hypothesis that the south bar is the result of a blowout . we have used optical spectroscopy to explore the velocity field in the vicinity of the proposed blowout region of the x - ray south bar . we obtained high - dispersion spectra at the cerro tololo inter - american observatory ( ctio ) using the echelle spectrograph with the long - focus red camera on the 4 m telescope . we used the instrument as a single - order , long - slit spectrograph by inserting a post - slit interference filter and replacing the cross - dispersing grating with a flat mirror so that a single echelle order around h@xmath6 and [ nii ] @xmath96548   , 6583    could be imaged over the entire slit . the 79 l / mm grating was used for these observations , and the data were recorded with a tek 2048@xmath102048 ccd using arcon 3.6 read - out . with a pixel size of 24@xmath11 m pixel@xmath0 , this ccd provided a sampling of 0.082    ( 3.75 km  s@xmath0 ) along the dispersion and 0.267 on the sky . spectral coverage was effectively limited by the bandwidth of the interference filter to 125   . coverage along the slit was 220 , including an unvignetted field of 2.5 the instrumental profile , as measured from th - ar calibration lamp lines , was 16.1 @xmath12 0.8 km  s@xmath0 fwhm . the ranges include variations in the mean focus and variations in the focus over a single spectrum . spatial resolution was determined by the seeing , which was 1 . the data were bias subtracted at the telescope . flat - fielding , wavelength calibration , distortion correction , and cosmic ray removal were performed using standard iraf routines . wavelength calibration and curvature correction were performed using th - ar lamp exposures taken in the beginning of each night . to eliminate velocity errors due to flexure in the spectrograph , the absolute wavelength scale in each spectrum was referenced to the telluric h@xmath6 airglow feature . figure [ halpha ] shows the echelle spectrum along with our recent  ccd image . the ccd image was taken at the curtis schmidt telescope at ctio . this image , along with several other filter images , was taken on february 18 , 1994 , using the thompson 1024@xmath101024 ccd . with this setup , the ccd has 1 pixels , giving a field of view of 31 . further details of the curtis schmidt observations can be found in smith et  al . the location of the echelle slit used for the spectrum is indicated on the ccd image . the star labeled `` star 1 '' shows up clearly at the center of the spectrum . several cuts in the velocity direction are plotted below the spectrum , and the heliocentric velocity scale is given to the left . the most important feature in this spectrum is the variation in the range of velocities seen . the regions to the east of star 1 have only a small range of velocities , roughly @xmath13 to @xmath14 km  sec@xmath0 , while to the west , there is clear evidence of significantly blue - shifted material , with velocity offsets up to @xmath15 km s@xmath0 . the region with this high - velocity feature coincides with the brightest part of the x - ray south bar . also in this area , optical filaments can be seen which extend away from the main optical shell of dem l152 . the high - velocity feature is not detected elsewhere in n44 ( meaburn & laspias 1991 ; hunter 1994 ) , except at another position within the x - ray south bar ( figure 15 of hunter 1994 ) . therefore , the high - velocity gas lends further support to the blowout nature of this region . if the south bar is indeed a blowout structure , with the gas escaping into a region of lower pressure beyond the superbubble shell , it is expected that the temperature of the hot escaping gas should drop as it undergoes adiabatic expansion . to test this hypothesis , we have used x - ray observations from both the rosat and asca x - ray observatories to compare the temperatures in shell 1 and the south bar . the rosat pspc observations have been discussed in some detail in c93 . we present here a brief summary : the observations were performed march 7  9 , 1992 , and are archived under observation number rp500093 . for these observations , n44 was centered in the pspc , where the angular resolution was 30 . the pspc is sensitive in the energy range of 0.1  2.4 kev and has an energy resolution of 43% at 1 kev . the total effective exposure time was 8871 seconds . our asca observations were performed march 1 , 1995 . the asca observatory has been discussed by tanaka et al . briefly , the satellite carries four imaging thin - foil grazing incidence x - ray telescopes . two of the telescopes are focused on solid - state ccd imaging spectrometers , called sis0 and sis1 . the other two telescopes use gas imaging spectrometers ( gis ) as the detectors . in this paper , we have only used the sis data because of the low sensitivity of the gis for energies below 1 kev , where our objects are the brightest , and the substantially lower energy resolution compared with the sis . the sis detectors have an energy resolution of 2% at 6 kev and cover the energy range 0.4  10 kev , with reduced throughput near the ends of the band . the spatial resolution is poor compared to the rosat pspc , with a narrow core of 1  diameter and a half - power radius of 3 . both of the sis detectors are composed of 4 separate ccds . not all ccds must be run for a given observation , and the presence of `` hot '' and `` flickering '' pixels , which fill the telemetry with false signals , have made it necessary to use only a subset of each detector for many observations . because of the telemetry limitations , our observations were performed in two - ccd mode , with the south bar placed at the standard point source location , near the center of the focal plane on sis0 , chip 1 ( sis1 chip 3 ; see figure [ xray]b ) . the observations were performed in faint mode when allowed by the telemetry , and converted to bright mode on the ground for analysis . using the standard processing software ( _ ascascreen _ & _ ftools _ ) , we removed bad time periods using the following selection criteria : aspect deviation @xmath80.01@xmath5 , angle to bright earth @xmath1620@xmath5 , elevation @xmath1610@xmath5 , minimum cutoff rigidity of 6 , pixl rejection of 75 . we also removed hot and flickering pixels . after screening , a total of 32,964 sec of usable data remained . we included high , medium , and low bit rate data . although low bit rate data have been seen to have difficulties in the past , we did not reject them as they constituted @xmath81% of our total integration time . we extracted spectra from the rosat and asca data using regions which were designed to be as similar as possible to those used in c93 . the background was chosen in a different location because the region used by c93 lies outside of the asca sis field of view . the background regions were chosen to exclude point sources and obvious diffuse emission which would contaminate the spectra . the background was taken from a source - free region on chip 0 of sis0 ( chip 2 of sis1 ) because there was no region on chip 1 of sis0 ( chip 3 of sis1 ) large enough to make a useful measurement . a comparison between a small source free region on chip 1 of sis0 ( chip 3 of sis1 ) showed no significant difference with the background we have used . it was necessary to use a local background region in the field because much of the background is contributed by diffuse emission from the lmc itself , in addition to the usual sources of background , such as charged particles , scattered solar light , and cosmic x - ray background . a comparison with the archive background fields shows significant excess emission in the local background field : 0.029 cts  sec@xmath0  for the local background vs 0.018 for the same region of the chip in the archive background field .    for the asca data , spectra were extracted from both sis0 and sis1 . the datasets from both sis telescopes and the pspc can not be simply combined directly ( i.e. , added together ) as the three telescopes have distinct response matrices . instead , the three spectra were fit jointly to the same models using xspec . the normalizations were allowed to fit separately to allow for differences in spatial and spectral resolution . spectral bins for all three instruments were add together in order to give sufficient counts in each spectral bin that the  statistics would be valid . optically thin thermal spectra from raymond & smith ( 1977 , hereafter rs ) and mewe - kaastra ( mewe et  al . 1985 ; kaastra 1992 , hereafter meka ) were used . the abundances of all elements except fe were set to 30% , consistent with the lmc abundances . because of some discrepancies in certain fe line ratios , we allowed the fe abundance to float to the best fit value , and fixed it at that value to determine the temperature confidence contours . in all fits , the fe abundances were in the range 1.5% to 7.9% solar , improving the fit substantially over models with the fe abundance set to 30% . we also experimented with letting other elemental abundances float , but no single element improved the fits substantially . we found that both rs and meka models with low fe abundances were in reasonable agreement with the observed spectra . derived parameters from the two models were in good agreement with each other . figure  [ s1sb ] presents the results of the spectral fits for shell 1 and the south bar using the meka model . the confidence contours and best - fit spectra are shown . figure  [ shell3 ] shows the confidence contours and best - fit spectra for shell 3 . the contours suggest that the expected temperature difference between shell 1 and the south bar is present : the blowout region ( south bar ) has a 99% confidence range for the temperature of 0.29  0.49 kev , while the main superbubble ( shell 1 ) has a 99% confidence range of 0.49  0.57 kev , with no overlap of the 99% confidence contours between the two fits . the supernova remnant proposed by c93 , shell 3 , does not show a significantly higher temperature as previously claimed  the 99% confidence interval is 0.24  0.46 kev . however , this does not refute the identification of this structure with an snr , and is perhaps more consistent with the observed size of the structure ( 30 pc ) and correspondingly large estimated age ( yr ) reported by c93 . the south bar and shell 1 have a low absorption , @xmath17 , consistent with the optical and radio measurements of the hydrogen column density ( see discussion in c93 ) . the best fit absorption for shell 3 is @xmath17 . although the error contour for shell 3 does not constrain the extinction strongly ( @xmath17 ) , it is worthwhile to note that the somewhat higher extinction of shell 3 , compared with shell 1 and the south bar is not unreasonable given the presence of faint dust lanes visible across shell 3 in the  images .    while these measurements appear to show quite a strong difference in the temperatures of the south bar and shell 1 , we caution the reader that the confidence contours do not completely describe the true errors involved . in particular the fact that , without allowing the fe abundance to float , the model spectra deviate systematically from the observed spectra is a point for concern . the problem of poorly - fitting fe lines has been noted already and blamed on errors in the atomic physics ( see e.g. , fabian et al . 1994 ; liedahl et al . several groups are currently working on updated model calculations incorporating newer atomic data in an effort to improve the situation . nonetheless , even without such improved models , the relative temperatures of the south bar and shell 1 determined with the meka and rs models are quite believable . one particular concern of the fits presented here is that there can be a small correlation between the fe abundance for a fit and the temperature . to show the effect of this correlation , which also plot confidence contours for the fit to the south bar , keeping the fe abundance fixed to the best fit value for shell 1 ( gray contours ) . the significance of the resulting temperature difference is reduced , but still present  now there is overlap between the 99% contours , but the 90% confidence contours do not touch . we have used optical echelle spectra along with rosat and asca x - ray spectra to test the hypothesis that the southern portion of the n44 x - ray bright region is the result of a blowout structure . three pieces of evidence now support this conclusion . first , the filamentary optical morphology corresponding with the location of the x - ray bright south bar suggests the blowout description ( c93 ) . second , optical echelle spectra show evidence of high velocity ( 90 km  sec@xmath0 ) clouds in the region of the blowout . third , x - ray spectral fits show a lower temperature for the south bar than the main superbubble region of shell 1 . if the south bar were the result of escaping gas expanding adiabatically into the surrounding region of lower pressure gas , the temperature would be expected to drop , as is observed . the combination of the optical and x - ray morphology , the high velocity gas seen in the echelle spectra , and these tentative x - ray temperature determinations give strong support to the interpretation of this structure of the south bar as a blowout from the main shell of the superbubble as suggested by c93 . such a blowout may explain the discrepancies discussed by om95 between the observed shell diameter and the diameter predicted on the basis of the stellar content and weaver et al.s ( 1977 ) pressure - driven bubble model . below we will discuss n44 s superbubble energetics and the effects of the blowout . first we use om95 s stellar content to estimate the total energy input to shell 1 of n44 . om95 have derived an age of 10 myr for the stars within the superbubble and 5 myr for the stars on the exterior . these ages will be used as references . we have integrated the stellar wind power implied by the stellar content ( figure 11 of om95 ) , and obtained a total stellar wind energy input of 3.2@xmath18 erg over the first 5 myr , and 3.5@xmath18 erg over the first 10 myr . apparently , the stellar wind input during the second 5 myr is only about 10% that of the first 5 myr . the initial mass function of the stars in n44 implies that 1  4 supernovae have occurred in the past 10 myr . it is likely that supernovae dominate the energy input in the second 5 myr as each explosion deposits erg . the total energy input from the stellar winds and supernovae is probably in the range of @xmath19 erg . we may use the model fits derived from the asca sis and rosat pspc data to determine more accurately the thermal energy in the superbubble interior . using a log n@xmath20 of 21.0 and a kt of 0.55 kev for the superbubble interior , the observed rosat flux gives a normalization factor log ( n@xmath21v/4@xmath22d@xmath23 ) = 11.05 in cgs units , where n@xmath24 is the electron density , v is the x - ray emitting volume , and d is the distance to the lmc ( 50 kpc ) . we use the rosat flux since the better angular resolution of the pspc implies that fewer photons are lost due to high - angle scattering . if we assume that the x - ray emitting volume is only 1/2 of the superbubble volume ( pc@xmath25 ) , we derive an n@xmath24 of 0.14 @xmath4 , a mass of 194 m@xmath26 , and a thermal energy of 2@xmath1010@xmath27 erg . if the x - ray emitting volume is only 1/100 of the superbubble volume , because of a small filling fraction , then n@xmath24 is increased to 1.0 @xmath4 , the mass is reduced to 28 m@xmath26 , and the thermal energy is reduced to 3@xmath1010@xmath28 erg . in weaver et al.s ( 1977 ) pressure - driven wind - blown bubble model , the thermal energy in the hot interior should be about 5/11 of the total stellar wind energy input . we see clearly that the thermal energy derived from our x - ray observations is at least an order of magnitude lower than expected . assuming the input energy from stellar winds and supernovae is correct , this energy discrepancy can be caused by two effects  energy leakage through the blowout and energy loss via radiation . to test the effect of energy loss via cooling , we note that the cooling timescale is roughly ( 1/n@xmath24 ) myr . for the possible range of n@xmath24 values , the cooling timescale would be a few myr ; therefore , we may expect a reduction of the thermal energy by a factor of 2 , or a few at most . another mechanism must account for the rest of the lost energy  one possibility is the blowout . the amount of energy lost in the blowout can be roughly estimated from the x - ray data . using the asca sis model fit of the blowout region , log n@xmath20 = 21.0 and kt = 0.41 kev , we obtain a normalization factor log ( n@xmath21 v / 4@xmath22d@xmath23 ) = 10.85 in cgs units . if we assume that the depth of the x - ray emission is the same as the width , the emitting volume is 4@xmath1010@xmath29 pc@xmath25 , n@xmath24 is 0.13 @xmath4 , the thermal energy is 1.0@xmath30 erg , and the mass is 133  m@xmath26 . note that the x - ray emitting volume could be overestimated by a factor of 2 or more , then the n@xmath24 would be a factor of @xmath31 higher , and the thermal energy and the mass would be a factor of @xmath31 lower . this amount of thermal energy , lost from the superbubble , is similar in magnitude to the thermal energy remaining in the superbubble . the largest uncertainty in this comparison is in the determinations of the emitting volume in the superbubble interior and the blowout . future high resolution imager images may help reduce the uncertainties . we now examine the thermal energy conversion problem with the blowout and radiative cooling taken into consideration . if we use the low estimates for the electron density , we obtain high estimates for the thermal energies in the superbubble interior and the blowout , 2@xmath1010@xmath27 erg and 1@xmath30 erg , but also high estimates for the radiative cooling timescale . applying the weaver et al . ( 1977 ) model , a total energy input of @xmath19 erg should have @xmath32 erg converted to thermal energy in hot gas . the observed thermal energy in the superbubble interior and the blowout together is only 3@xmath1010@xmath27 erg . even if we assume that 1/2 of the thermal energy has been radiated away ( thereby ignoring the high estimate for the cooling timescale ) , there is still a discrepancy by a factor of 3 - 5 . if we use the high estimates for the density , we would have much lower estimates for the thermal energies , and the discrepancy between observed and expected thermal energies would be even larger . therefore , we conclude that either the weaver et al . ( 1977 ) model is drastically wrong or the stellar wind energy input must have been over - estimated  we consider the latter to be more likely . a similar conclusion has been reached by garca - segura & mac low ( 1995 ) in their modeling of the single - star wind - blown bubble ngc 6888 ; they suggest that stellar wind strengths could have been over - estimated by a factor of @xmath33 because stellar winds are clumpy ( moffatt & robert 1994 ) .    besides the interior thermal energy , another important reservoir of energy in a superbubble is the expansion of the cool  shell . the radius of n44 superbubble is 30 pc , and the expansion velocity is @xmath3440 km s@xmath0 ( meaburn & laspias 1991 ) . the density in the shell is unfortunately in the low - density regime for the doublet , so it is necessary to follow the procedure outlined by chu et al . ( 1995 ) to derive a rms density from the peak emission measure along the shell rim . the shell thickness is calculated assuming an isothermal shock into a photo - ionized interstellar medium . for a peak emission measure of 1500  4000 @xmath35 pc along the rim of the shell ( c93 ) , we derive a shell thickness of 0.5 pc , a rms electron density of 12  20 @xmath4 in the  shell , an ambient density of 0.6  1.0 @xmath4 , a shell mass of @xmath36 m@xmath26 , and a shell kinetic energy of @xmath37 erg . this kinetic energy is at least an order of magnitude lower than expected from weaver et al.s model assuming a total energy input of @xmath19 erg . this discrepancy may also be caused by the energy leakage and radiation losses , or by an over - estimate of the stellar wind strengths as discussed above . the presence of a significant blowout in n44 , and the resulting effect on the evolution of the superbubble , has important implications for the hot phase of the ism and model calculations of superbubbles in galaxies . it is interesting to note that the blowout in n44 is probably within the gaseous disk , instead of into a halo . this assumption is based on the small size of n44 compared with the the scale height of the  layer , as well as the apparent interaction with the  regions to the south . if such blowouts are common in superbubbles in general , theoretical estimates of superbubble sizes will be systematically too large . such an error would affect model predictions of the ism distribution in galaxies . it is reasonable to conclude that blowouts of this nature may be common : the blowout discussed in this paper is not trivial to detect , and depended upon the x - ray morphology information . since n44 is the brightest superbubble in x - rays in the lmc , the blowout in this case is relatively easy to detect . for fainter superbubbles , such blowouts could be common and may simply have gone unnoticed so far . future work to search for evidence of other , similar blowout structures should be performed . we acknowledge the support of nasa grants nag 5 - 1900 , nag 5 - 2679 , nag 5 - 2973 , and nagw-4519 . sdp acknowledges support by a national science foundation graduate research fellowship . thanks to joel parker for helping to start our collaboration . = 1    castor , j. , mccray , r. , & weaver , r. 1975 , apj  200 , l107    = 1    chu y .- h . , chang h .- w . , su , & mac low m .- 1995 , apj  450 , 157    = 1    chu y .- h . & mac low m .- 1990 , apj  365 , 510    = 1    chu y .- h . , mac low m .- m . , garcia - segura g. , wakker b. , & kenicutt r.  c. 1993 , apj  414 , 213 ( c93 )    = 1    davies r.  d. , elliott k.  h. , & meaburn j. 1976 , mnras  81 , 89    = 1    fabian a.  c. , arnaud k.  a. , bautz m.  w. , & tawara y. 1994 , apj  436 , l63    = 1    feast m.  w. 1991 , in iau symp . 148 , the magellanic clouds , ed . r. haynes & d. milne ( dordrecht : kluwer ) , 1    = 1    garca - segura g. & mac low , m .- m . 1995 , apj  455 , 145    = 1    henize k.  g. 1956 , apjs  2 , 315    = 1    kaastra j.  s. 1992 , an x - ray spectral code for optically thin plasmas ( sron - leiden report ) . ( meka )    = 1    liedhal d.  a. , osterheld a.  l. , & goldstein w.  h. 1995 , apj  338 , l115    = 1    lucke p.  b. & hodge p.  w. 1970 , aj  75 , 171    = 1    meaburn j. & laspias n.  v. 1991 , a&a  245 , 635    = 1    mewe r. , gronenschild e.  h.  b.  m. , & van den oord g.  h.  j. 1985 , a&a  62 , 197 ( meka )    = 1    oey m.  s. & massey p. 1995 , apj  452 , 216    = 1    raymond j.  c. & smith b.  w. 1977 , apjs  35 , 419 ( rs )    = 1    smith r.  c. , chu y .- h . , mac low m .- m . , oey m.  s. , & klein u. 1994 , aj  108 , 1266    = 1    tanaka y. , inoue h. , & holt s.  s. 1994 , pasj  46 , l37    = 1    wang q. & helfand d.  j. 1991 , apj  373 , 497    = 1    weaver r. , mccray r. , castor j. , shapiro p. , & moore r. 1977 , apj  218 , 377 figure  [ halpha ] : top )  image of the n44 complex with rosat x - ray contours overlaid . the major regions of nebulosity are labeled ( henize 1956 ; dem ) . the location of the echelle slit is noted , along with the position of star 1 . bottom ) the echelle spectrum , including five representative cuts in the spectral direction . the velocity axis is at the left , showing the heliocentric velocity . in this figure , the contours are draw in pairs , with a white line at a slightly lower value than the black line . this allows the reader to distinguish positive and negative contours . figure  [ xray ] : top ) rosat x - ray image , with  contours from our ccd image overlaid . the image has been smoothed with a gaussian with a fwhm of 45 . four of the features described by chu et  al . , ( 1993 ) are labeled , and the regions from which spectra were extracted are marked . bottom ) asca x - ray image , with  contours overlaid . this image is at the same scale and orientation as ( a ) , and has also been smoothed with a gaussian with a fwhm of 60 . in both figures , the contours are draw in pairs , with a white line at a slightly lower value than the black line . this allows the reader to distinguish positive and negative contours . figure  [ s1sb ] : comparison of shell 1 and south bar . top ) confidence contours for the spectral fits for shell  1 ( thick , black line ) and the south bar ( thick , gray line ) in which the fe abundance is fixed at the best fit value . also plotted is the confidence contours for the spectral fit for the south bar , using the best - fit fe abundance value for shell  1 . the three sets of contours represent 68% , 90% , and 99% confidence levels . bottom ) observed spectra ( error bars ) and meka fits ( histograms ) for shell 1 ( left ) and the south bar ( right ) . in both of the figures , we have included data from asca sis 0 ( thick line ) , sis 1 ( thin line ) , and rosat pspc ( medium weight line ) . figure  [ shell3 ] : left ) confidence contours for the spectral fits for shell  3 in which the fe abundance is fixed at the best fit value . the three sets of contours represent 68% , 90% , and 99% confidence levels . right ) observed spectra ( error bars ) and meka fits ( histograms ) for shell 3 . in this figures , we have included data from asca sis 0 ( thick line ) , sis 1 ( thin line ) , and rosat pspc ( medium weight line ) .
Abstract:

Answer: we have used optical echelle spectra along with rosat and asca x - ray spectra to test the hypothesis that the southern portion of the n44 x - ray bright region is the result of a blowout structure . three pieces of evidence now support this conclusion . first , the filamentary optical morphology corresponding with the location of the x - ray bright south bar suggests the blowout description ( chu et  al.1993 ) . second , optical echelle spectra show evidence of high velocity ( 90  km  sec@xmath0 ) gas in the region of the blowout . third , x - ray spectral fits show a lower temperature for the south bar than the main superbubble region of shell  1 . such a blowout can affect the evolution of the superbubble and explain some of the discrepancy discussed by oey & massey ( 1995 ) between the observed shell diameter and the diameter predicted on the basis of the stellar content and weaver et al.s ( 1977 ) pressure - driven bubble model .    # 1#2#1@xmath1    -0.5 in 6.5 in