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: there have been a number of attempts to study hypernuclei containing not only nucleons but also hyperons with strangeness @xcite . it has been of high interest to investigate how the presence of strangeness broadens and/or modifies the knowledge achieved in the field of ordinary nuclear physics . it is a fundamental problem to understand properties of a quantum many - body system in terms of basic interactions between constituent particles . in hypernuclear physics a unified treatment of all the baryon - baryon interactions , especially nucleon - nucleon ( @xmath4 ) and hyperon - nucleon ( @xmath5 ) interactions , will be necessary if one wish to understand general properties of hypernuclei . however , our knowledge on the @xmath5 interaction is quite inadequate in contrast to that of the @xmath6 interaction due to experimental difficulties . in such a situation many - body theoretical studies of hypernuclei could provide alternative information on the @xmath5 interaction . it has been discussed that physics of hypernuclei is different from conventional nuclear physics in some aspects . first , the channel couplings such as the @xmath7-@xmath8 coupling play a significant role in the structure of hypernuclei in contrast to the very limited role of the @xmath9-@xmath10 coupling appeared in nuclear many - body systems @xcite . second , the anti - symmetric spin - orbit force , which is forbidden in ordinary nuclei , emerges in hypernuclei in addition to the symmetric spin - orbit force . the spin - orbit splitting of @xmath11 single - particle levels in @xmath11 hypernuclei was considered experimentally to be much smaller than their nucleonic counterpart @xcite . however , some recent studies as in ref . @xcite suggest that the splitting would not be so small as discussed earlier . it will be highly desirable for the spin - orbit splitting to be established in experiment . the magnitude of the spin - orbit splitting will put further constraints on the @xmath12 interaction . there have been several possibilities of solving the structure of hypernuclei . one possibility is to treat few - body or light hypernuclei by solving the faddeev equation @xcite or by applying the cluster model @xcite . another possibility is microscopic studies of hypernuclei with a larger mass number . while a few lower partial waves play a decisive role in few - body or light hypernuclei , higher partial waves could give some contributions in heavier hypernuclei . these microscopic studies of hypernuclei could make us understand interesting physics of their structure itself , and also help us in determining the @xmath5 interaction . microscopic structure calculation of heavier hypernuclei is usually performed by introducing a @xmath12 effective interaction . for this purpose the @xmath13-matrix theory has been applied in two ways : in some works an effective @xmath5 interaction has been derived by simulating a nuclear matter @xmath13-matrix . in this approach the non - locality and energy dependence of the @xmath13-matrix are neglected , and the effective @xmath5 interaction is represented in the two- or three - range gaussian form with the adjustable fermi momentum @xmath14 . this type of the effective interaction has been referred to as the yng interaction @xcite . the yng interaction has been successfully used in many of hypernuclear structure calculations @xcite . it seems to us that the yng interaction has made an important contribution to the considerable progress in overall description of hypernuclei . the other approach , in the framework of the @xmath13-matrix theory , is to calculate the @xmath13-matrix in finite hypernuclei by treating rigorously the pauli exclusion operator and making a self - consistent calculation of single - particle potential of hyperons @xcite . in contrast to studies of using the yng interaction this approach possibly gives a better description of the state dependence of the effective @xmath5 interaction . it has been established that the @xmath13-matrix can be used as a basic ingredient in constructing a microscopic many - body theory starting with bare interactions .    in spite of a marked success in the applications of the @xmath13-matrix theory , we may say from a formal point of view that the @xmath13-matrix is energy - dependent and non - hermitian , and does not have the property of decoupling between a model space of low - lying two - particle states and its complement . the @xmath13-matrix is not considered to be an effective interaction in a formal sense of the effective interaction theory . in order to derive an effective interaction we should add some higher - order corrections such as folded diagrams @xcite . it would be desirable if we could have a theory describing many - body systems in terms of an energy - independent and hermitian effective interaction with the property of decoupling . two of the present authors , k.  s. and r.  o. , proposed a many - body theory , the unitary - model - operator approach ( umoa ) that was formulated on the basis of such an effective interaction . the umoa was applied to finite nuclei , @xmath15o and @xmath16ca @xcite , and some improvements have been attained . the method may be viewed as an alternative way of deriving an effective interaction in finite many - body systems . a unified study of hypernuclei and ordinary nuclei would be helpful in understanding many - body systems of baryons in a microscopic way . therefore , we extended the formulation of the umoa , in the previous work @xcite , to the description of hypernuclei and made a calculation of the properties in @xmath17o with a realistic @xmath12 interaction . the purposes of this work are to present a general formulation of the extended umoa and apply it to @xmath17o by employing some realistic @xmath12 interactions . this paper is organized as follows : in section 2 we present a general formulation of the umoa for deriving the @xmath18 effective interaction . in section 3 the approximation procedure is given for making actual calculations . in section 4 we apply the umoa to @xmath17o by employing the nijmegen soft - core ( nsc ) @xcite , jlich model -  ( j ) and model - b ( jb ) @xcite @xmath12 potentials . in section 5 we make some concluding remarks . we present a formulation of the umoa for applying to a calculation of the properties in @xmath11 hypernuclei , especially , closed - shell nucleus plus one @xmath11 systems . a main purpose is to give a method of calculating an effective interaction between @xmath11 and @xmath19 . in the derivation of the @xmath8 effective interaction we should note that , because of the small mass difference of about 77mev between @xmath11 and @xmath20 , the coupling of the @xmath8 and @xmath7 channels plays an important role . we first consider a two - body hamiltonian of a @xmath7-@xmath8 coupled system given by @xmath21 where @xmath22 @xmath23 and @xmath24 the terms @xmath25 and @xmath26 for @xmath27 and @xmath19 are the kinetic and single - particle potential energies of @xmath28 and @xmath19 , respectively , and @xmath29 is the bare two - body @xmath12 interaction including the @xmath8 and @xmath7 channels . the terms @xmath30 and @xmath31 are the rest masses of @xmath32 and @xmath33 , respectively . the hamiltonian @xmath34 describes a subsystem of interacting two particles @xmath8 or @xmath7 in the one - body potentials @xmath35 , @xmath36 and @xmath37 .    in order to derive a @xmath8 effective interaction we introduce a unitarily transformed @xmath5 interaction given by @xmath38 where @xmath39 is the correlation operator defined in a space of two - particle states of @xmath8 and @xmath7 . the @xmath39 is an anti - hermitian operator satisfying @xmath40 in the formulation of the umoa the correlation operator @xmath39 is determined by the equation of decoupling for @xmath41 between a certain model space consisting of low - momentum two - particle states and its complement . this determination of the correlation operator @xmath39 is one of the characteristics of the present approach .    for the determination of a @xmath8 effective interaction , we introduce a model space consisting of low - momentum @xmath8 states . let @xmath42 be the projection operator onto the model space . the complement of the @xmath43 space includes a space of high - momentum @xmath8 states referred to as the @xmath44 space and a space of all the @xmath7 states referred to as the @xmath45 space . the condition of decoupling for the transformed interaction @xmath41 is given by @xmath46 hereafter , we refer to the above equation as the decoupling equation . we assume that the one - body part @xmath47 in eq . ( [ eq : tvyn ] ) is decoupled as @xmath48 on the above assumption the decoupling equation ( [ eq : dec ] ) becomes @xmath49 in general , the decoupling equation ( [ eq : decr ] ) does not determine @xmath39 uniquely . the usual restrictive conditions are @xmath50 these conditions are called the minimal effect requirements @xcite , which mean that the transformation exp@xmath51 does not induce unnecessary transformation within each of the @xmath52 and @xmath53 spaces .    in the general theory of constructing a hermitian effective interaction , the solution for @xmath54 in eq . ( [ eq : decr ] ) has been known @xcite and is given by @xmath55 with latexmath:[\ ] ] the @xmath59 in eq . ( [ eq : oyn ] ) is the biorthogonal state of @xmath60 satisfying latexmath:[\[\begin{aligned } \label{eq : bio } \langle \overline{\psi_{{\sst { \it \lambda n}}k } }    where @xmath60 is the @xmath52-space component of @xmath57 defined by latexmath:[\[\begin{aligned } \label{eq : psc }     that the solution @xmath39 in eq . ( [ eq : syn ] ) is determined dependently on the choice of a set of @xmath63 eigenstates @xmath64 as given in eq . ( [ eq : oyn ] ) . we choose a set of @xmath65 so that they have the largest @xmath52-space overlaps among all the eigenstates in eq . ( [ eq : eige ] ) . this is a usual choice in the derivation of the effective interaction .    with the solution @xmath39 in eq . ( [ eq : syn ] ) , the @xmath8 effective interaction in the @xmath52 space is given by @xmath66 it has been known that this @xmath8 effective interaction can be written explicitly @xcite as @xmath67 where @xmath68 ( @xmath69 ) and @xmath70 ( @xmath71 ) are defined through the eigenvalue equation @xmath72 and @xmath73 is given by @xmath74 the @xmath75 is a @xmath42-space operator that agrees with the effective interaction of non - hermitian type . equation ( [ eq : evln ] ) is a formula for converting the non - hermitian effective interaction to the hermitian one . the term @xmath76 in eq . ( [ eq : shyn ] ) is introduced as a self - consistent potential of @xmath32 calculated with the @xmath8 effective interaction @xmath77 as @xmath78 ( @xmath79 ) and @xmath80 are the single - particle states of @xmath32 and @xmath19 , respectively .    from eqs . ( [ eq : evln])-([eq : nhevln ] ) we see that the @xmath8 effective interaction @xmath81 can be calculated once the operator @xmath82 is given . the operator @xmath82 is determined completely , as given in eq . ( [ eq : oyn ] ) , by solving the eigenvalue equation ( [ eq : eige ] ) and selecting @xmath63 eigenstates @xmath65 with the largest @xmath52-space overlaps . it is noted here that eqs . ( [ eq : syn])-([eq : ul ] ) give a set of equations that determines the @xmath8 effective interaction @xmath77 self - consistently with the single - particle potential @xmath76 . here we discuss a calculation procedure for obtaining the @xmath8 effective interaction acting in a small model space of low - lying @xmath8 states . we use the harmonic oscillator ( h.o . ) wave functions as basis states . it is desirable that the single - particle potentials of @xmath83 , @xmath33 and @xmath19 in intermediate states are taken into account in a self - consistent way . however , as our main concern is to derive a @xmath8 effective interaction , the self - consistent treatment of @xmath35 would have predominant importance in the present calculation . for this reason we make a self - consistent calculation only for @xmath35 , and as for @xmath37 we employ a fixed potential calculated previously for closed shell nuclei in the umoa @xcite . the effect of the single - particle potential @xmath36 on the @xmath8 effective interaction has not necessarily been clarified , although some calculations with @xmath36 have been made for the @xmath13-matrices in hyperonic nuclear matter @xcite . in the actual calculations made so far the single - particle potential @xmath36 has been considered to be rather shallow @xcite . although we have not known clearly how much the effect of @xmath36 on the @xmath8 effective interaction , we here assume @xmath84 .    in solving the self - consistent equation for @xmath35 with @xmath81 we should be careful in the following two respects : ( a ) in the calculation of energies of intermediate states in propagators we should use @xmath35 that is given for @xmath83 states including high - momentum states . in order to do this we need a @xmath8 effective interaction acting in a large model space , because in general the @xmath44 space contains high - momentum @xmath8 states and the @xmath85-space effective interaction @xmath86 is not always well behaved . we thus make a calculation by introducing two model spaces , namely , large and small model spaces . the @xmath8 effective interaction in each of the model spaces is calculated by following a two - step procedure . in the first - step procedure , we define a large model space and calculate the @xmath8 effective interaction and single - particle potential of @xmath83 in high - lying states . the @xmath8 effective interaction acting in this large model space will not be suitable for the effective interaction in a usual sense . by introducing a small model space and using the @xmath8 effective interaction determined in the first - step calculation , we proceed to the second - step calculation where we calculate a @xmath8 effective interaction acting in this small model space of low - lying @xmath87 states . ( b ) it is important in actual calculations how to determine the @xmath88 and @xmath85 spaces because an effective interaction is determined dependently on the choice of the @xmath88 and @xmath85 spaces . if a state in the @xmath88 space mixes strongly with @xmath85-space states in the eigenstate in eq . ( [ eq : eige ] ) , the matrix element of @xmath89 or @xmath54 defined in sec . 2 becomes large and as a result higher - order many - body correlations @xcite give rise to non - negligible contributions . in general , the strong mixing takes place when some of the @xmath88- and @xmath85-space states are quasi - degenerate in energy . for this reason , we should choose the @xmath88 and @xmath85 spaces which are well separated in energy .    with due regard to points ( a ) and ( b ) we perform a calculation of the effective interaction through the two - step procedure . in the following subsections we present the two - step method in detail . let us write a two - particle state of @xmath83 and @xmath19 as the product of the h.o . single - particle wave functions in the usual notation as @xmath90 we define a model space denoted by @xmath91 and its complement @xmath92 with a boundary number @xmath93 as @xmath94 we also define a space of all the @xmath7 states denoted by @xmath95 as @xmath96 note that the model space @xmath91 does not contain the @xmath7 states . if we make a graph of two - particle states , the @xmath97 , @xmath92 and @xmath98 spaces are represented in fig .  1 . the @xmath93 means a number such that @xmath99 gives the maximum kinetic energy of the @xmath8 two - particle states in the @xmath97 space . the @xmath100 denotes the unit of the h.o . energy which we use commonly for @xmath32 , @xmath33 and @xmath19 . as we require to introduce a sufficiently large model space in the first - step calculation , we should choose @xmath93 as large as possible . another requirement is that the states in the @xmath91 space should be separated in energy from the states in the @xmath101 space . for this reason , we adopt a model space @xmath97 of a triangular shape . this choice of the @xmath97 space guarantees the separation in energy . [ fig:1 ] = 13 cm    in fig . 1 the spaces @xmath102 for @xmath103 and @xmath33 , specified by the numbers @xmath104 and @xmath105 , stand for the excluded spaces due to the pauli principle for nucleons . the number @xmath106 is defined as @xmath107 with the quantum numbers @xmath108 of the highest occupied orbit in the closed - shell core nucleus . the @xmath109 means a number such that @xmath110 gives the maximum kinetic energy of the @xmath12 two - particle states in the @xmath111 and @xmath112 spaces . we then write the @xmath111 and @xmath112 spaces as @xmath113 and @xmath114 we also define the sum of the excluded spaces as @xmath115 furthermore , by using the projection operators , @xmath101 , @xmath116 , @xmath95 and @xmath117 , we introduce projection operators onto the spaces of allowed @xmath12 states as @xmath118 @xmath119 and we also define @xmath120 considering the pauli - exclusion principle for nucleons the terms @xmath121 and @xmath122 should be removed in treating the problem of two - body correlations . however , if we remove all the matrix elements of these terms , this leads to overmuch counting of the pauli principle . the matrix elements which are diagonal in nucleon states should be restored even though nucleons are in occupied states . therefore , instead of @xmath123 , we should use a @xmath12 interaction @xmath124 defined as @xmath125 where @xmath126 denotes the set of the h.o . quantum numbers @xmath127 of an occupied state of a nucleon , and the operator @xmath128 is defined as @xmath129    the two - body @xmath12 equation to be solved in the first - step calculation is now written as @xmath130 where we define @xmath131 as @xmath132 the @xmath133 is the single - particle potential of @xmath32 to be calculated with the @xmath8 effective interaction self - consistently . the number @xmath134 stands for a set of quantum numbers to specify a two - body @xmath12 eigenstate . the unitarily transformed @xmath12 interaction in the first - step calculation is given by @xmath135 the decoupling equation for determining the correlation operator @xmath136 becomes @xmath137 the @xmath8 effective interaction @xmath138 in the @xmath97 space is then given by @xmath139 in the first - step procedure we solve the decoupling equation ( [ eq : dece1 ] ) by introducing some assumptions . as we take a sufficiently large number as @xmath140 , we may assume that the one - body potential of @xmath83 for a state in the @xmath101 space is negligible in comparison with the sum of the energies , @xmath141 . therefore , @xmath142 takes nonzero values only for states in the @xmath91 space , and we have @xmath143 the definitions of @xmath97 and @xmath92 in eq . ( [ eq : dtps ] ) allow us to divide a space of @xmath8 states written in relative and center - of - mass ( c.m . ) coordinates into the same subspaces @xmath97 and @xmath92 as @xmath144 where @xmath145 and @xmath146 are the h.o . quantum numbers of relative and c.m . motions , respectively . the @xmath147 and @xmath148 denote , respectively , spin and angular momentum defined through @xmath149 . the isospin @xmath150 is omitted because @xmath150 takes a constant @xmath151 in the @xmath8 system . a @xmath7 state in the @xmath152 space in fig .  1 is also written by the product of the h.o . relative and c.m . states . the @xmath5 interaction @xmath153 in eq . ( [ eq : vbyn ] ) is not diagonal in the c.m . quantum numbers @xmath154 because the projection operators @xmath97 , @xmath155 and @xmath156 are not diagonal in @xmath154 . thus , the @xmath8 effective interaction @xmath157 is not diagonal in @xmath154 . it has been shown , however , in the calculation of ordinary nuclei that @xmath154 non - diagonal elements of the @xmath13-matrix are generally small @xcite . therefore , we make an assumption that the @xmath154 non - diagonal matrix elements of @xmath157 can be neglected . the two - body eigenvalue equation ( [ eq : eige1 ] ) is solved for the relative motion with a fixed set of @xmath158 written by @xmath159 where we rewrite @xmath160 as @xmath161 the @xmath162 and @xmath163 are the kinetic energies of the @xmath5 relative and c.m . motions , respectively . the state @xmath164 expresses a superposition of the @xmath8 and @xmath7 h.o . states and gives a general form of the eigenstate when we consider the @xmath7-@xmath8 coupling and the @xmath12 interaction with the tensor and antisymmetric spin - orbit couplings . the @xmath165 is an additional quantum number specifying a two - body @xmath5 eigenstate .    in general , the one - body part @xmath166 does not satisfies the assumption of decoupling in eq . ( [ eq : assume ] ) . however , it has been shown in the calculation in nuclei that if we take a larger number as @xmath167 , the effect of violating the assumption of decoupling becomes considerably smaller @xcite . we have confirmed that this non - decoupling effect becomes negligible in calculating nuclear effective interactions . we here assume this to be the case also in the determination of the @xmath8 effective interaction . furthermore , we introduce another approximation that the sum of @xmath168 and @xmath169 be diagonal in the c.m . quantum numbers @xmath170 and @xmath171 . this assumption has also been made sure to be acceptable in the previous works @xcite .    in order to solve the eigenvalue equation ( [ eq : eige1r ] ) we need to calculate the matrix elements of @xmath172 in the relative and c.m . coordinates representation . however , the calculation of these matrix elements is difficult because the term @xmath172 can not be separated into a sum of operators written in relative and c.m . coordinates . therefore , as have made in the previous work @xcite , we employ the angle - average approximation and have @xmath173 where , as was assumed previously , @xmath174 for @xmath175 . the @xmath176 is the h.o . transformation bracket @xcite for the @xmath12 two - body state which depends on the rest masses for @xmath177 and @xmath33 . here the angle - average approximation means that , when a function @xmath178 is given with an angular momentum @xmath179 and its @xmath180-component @xmath181 , we approximate @xmath178 as @xmath182 . in the present case the matrix element on the left - hand side in eq . ( [ eq : uyn ] ) is dependent on the @xmath180-components @xmath183 and @xmath184 of @xmath148 and @xmath171 , respectively . the average over @xmath183 and @xmath184 leads to the expression given in eq . ( [ eq : uyn ] ) . this approximation is essentially the same as has been made to derive an angle - averaged pauli operator by wong and sauer @xcite . the use of the basis states @xmath185 makes it difficult to represent the operator @xmath155 , but it has been confirmed that the angle - average technique works well . we thus make an approximation to the operator @xmath155 as @xmath186 where the summation should be made on the conditions @xmath187 here the coefficient @xmath188 is defined as @xmath189 where in the summation the following conditions should be satisfied @xmath190 furthermore , we need an expression of @xmath191 given in eq . ( [ eq : qynd ] ) in the angle - average approximation @xmath192 where the summation should be made on the same conditions as in eq . ( [ eq : sumc1 ] ) . the coefficient @xmath193 is given by @xmath194 where @xmath195 and @xmath196 should satisfy the conditions given in eq . ( [ eq : sumc2 ] ) .    once the set of the eigenstates @xmath197 in eq . ( [ eq : eige1r ] ) is solved , the @xmath8 effective interaction @xmath198 acting in the @xmath97 space is given according to eqs . ( [ eq : evln])-([eq : nhevln ] ) by substituting @xmath199 , where @xmath200 is defined as latexmath:[\[\begin{aligned } \label{eq : oyn1 } \hspace{-5mm}\omega _ { \sst { \it yn}}^{\sst ( 1 ) } = \sum_{k_{1},ll'ss'j_{r},n_{c}l_{c } } \overline{q}_{\sst { \it yn}}^{\sst ( 1 ) }    \langle \overline{{\it yn};k_{1},ll'ss'j_{r},n_{c}l_{c } }    the @xmath8 effective interaction @xmath157 thus determined , we calculate a new single - particle potential @xmath202 as an averaged matrix element of @xmath203 . the @xmath202 can be given explicitly in terms of the effective interaction @xmath157 as @xmath204 in the above treatment of the single - particle potential of @xmath32 in the first - step calculation , we do not consider the one - body spin - orbit splitting . the @xmath8 effective interaction @xmath205 in eq . ( [ eq : tvln1 ] ) is dependent on the single - particle potential @xmath133 , and @xmath133 is determined with @xmath205 as in eq . ( [ eq : ul1 ] ) . therefore , the equation for determining @xmath205 and @xmath133 should be solved simultaneously . the model space @xmath91 in the first - step procedure is chosen as a _ low - momentum _ space , but the @xmath91 space is still large when we wish to use @xmath205 in the calculation of bound or low excited states of @xmath32 hypernuclei . with @xmath205 , we further proceed to the calculation of the @xmath8 effective interaction acting in a smaller model space and correspondingly the single - particle potential of @xmath32 in low - lying states . we separate the @xmath91 space into the @xmath206 , @xmath207 , @xmath208 and @xmath209 spaces specified by the numbers @xmath140 , @xmath106 , and @xmath210 as shown in fig .  2 . the number @xmath211 is introduced to specify the uppermost bound state of @xmath32 . here these spaces are defined as @xmath212 the @xmath213 space is the excluded space due to the pauli principle for nucleons . we are principally interested in the @xmath8 effective interaction that gives a description of the bound states of @xmath83 . it would be desirable that such a @xmath8 effective interaction satisfies the condition of decoupling between two @xmath8 states consisting of the bound and unbound states of @xmath83 . this requirement is equivalent to the situation that the @xmath8 effective interactions should be decoupled between two @xmath8 states in the @xmath214 and @xmath209 spaces . therefore , we solve the decoupling equation for @xmath215 between the @xmath8 states in the @xmath214 and @xmath209 spaces . the space @xmath208 is not a space to be excluded due to the pauli principle . however , some states in the @xmath208 and @xmath214 spaces overlap each other in energy . therefore , we treat separately the effect of the coupling between the @xmath208 and @xmath214 spaces . actually the interaction between states in the @xmath208 and @xmath214 spaces induces the core polarization in the core nucleus , as discussed in refs . @xcite and @xcite . we take into account the effect of the core polarization in the usual perturbative calculation . = 6 cm    in the second - step calculation the unitarily transformed interaction is given by @xmath216 where @xmath217 the term @xmath218 is the self - consistent potential of @xmath32 to be determined in the second - step procedure . the decoupling equation for determining the correlation operator @xmath219 in eq . ( [ eq : tvyn2 ] ) becomes @xmath220 we note here that the assumption of decoupling @xmath221 is satisfied exactly because the two two - body states in the @xmath214 and @xmath209 spaces do not contain a common single - particle state of @xmath83 or @xmath19 . the two - body eigenvalue equation to be solved in the second - step calculation is written as @xmath222 where @xmath223 denotes the total angular momentum , and @xmath224 a set of quantum numbers of a @xmath8 two - body eigenstate . it should be noted that this eigenvalue equation can be solved exactly in the @xmath225 space by the diagonalization of the hamiltonian written with the basis states @xmath226 if the @xmath227 space is chosen suitably as a small space . once a set of the eigenstates @xmath228 in eq . ( [ eq : eige2 ] ) is solved , the @xmath8 effective interaction @xmath229 in the @xmath214 space is given from eqs . ( [ eq : evln])-([eq : nhevln ] ) with @xmath230 defined by @xmath231 is calculated with the @xmath8 effective interactions , @xmath229 and @xmath157 , as @xmath232 we solve the set of equations for determining @xmath229 and @xmath233 iteratively until the result converges . we remark here that we employ @xmath157 as the effective interaction in the @xmath213 space in the calculation of @xmath233 . since the effective interaction in the @xmath213 space is unchanged in the transformation made in the second - step procedure , we may write the @xmath8 effective interaction in the @xmath213 space as @xmath234 . we , therefore , have used partly the @xmath87 effective interaction @xmath157 in the calculation of @xmath233 as in eq . ( [ eq : ul2 ] ) . first we note that the one - body hamiltonian @xmath235 in eq . ( [ eq : hln2 ] ) has non - diagonal terms because we did not make the hartree - fock calculation . these terms could bring about sizable corrections to the single - particle energies of @xmath236 . we adopt the same h.o . energy both for @xmath83 and @xmath19 , though it has been considered that the wave function of @xmath83 is basically different from that of @xmath237 . therefore , we must evaluate the corrections arising from the non - diagonal terms of @xmath235 . by taking the diagonal part of @xmath235 as the unperturbed hamiltonian we calculate these corrections perturbatively up to second order as shown in diagram ( a ) of fig .  3 .    next , the interaction between the @xmath8 states in the @xmath214 and @xmath208 spaces induces the core polarization in the core nucleus , and is important in the calculation of the corrections . we here assume that these corrections can be calculated approximately , using the interaction @xmath238 it will be desirable to take into account the effect of intermediate nucleon states with high momentum . in this perturbative calculation we include only the nucleon states belonging to the @xmath208 space , that is , the nucleon h.o . state should satisfy @xmath239 , where @xmath240 are the h.o . quantum numbers of a @xmath83 state . the validity of this truncation of nucleon states will be checked in the numerical calculation in the next section by increasing the number @xmath241 . the diagrams corresponding these corrections are given in diagram ( b ) and ( c ) of fig . we performed calculations employing the jlich model -  ( j ) , model - b ( jb ) @xcite potentials in the momentum - space representation made by b.  holzenkamp and a.  reuber . ] and nijmegen soft - core ( nsc ) @xcite potentials for the _ yn _ interaction . we used the h.o . basis states with a common @xmath242mev for @xmath32 , @xmath33 and @xmath19 . the dependence of the effective interaction and single - particle energies on the value @xmath100 has been already discussed in @xmath15o @xcite and we confirmed that the dependence was quite small around @xmath243mev if perturbation corrections were included . as for the nucleon single - particle potential @xmath244 , we used the fixed data @xcite calculated using the paris potential @xcite , and we have assumed that @xmath245 . the boundary numbers @xmath246 , @xmath247 and @xmath109 are taken as @xmath248 , @xmath249 and @xmath250 for @xmath0o . = 13.8 cm    we present , in fig .  4 , the @xmath93 dependence of the single - particle energies of @xmath32 including the perturbation corrections shown in fig .  3 . the parameter @xmath93 specifies the @xmath91 space in which the @xmath8 effective interaction @xmath229 is determined rigorously without any approximation as discussed in sec . 3 . therefore , as @xmath140 becomes larger , the calculated result does more accurate . if the result converges at a certain number @xmath251 , we may say that the method of the two - step calculation of the @xmath8 effective interaction would be trustworthy . as shown in fig .  4 , the results are convergent and stable for @xmath252 for each @xmath12 potential . if we take @xmath93 larger than 6 , our model space is considered to be sufficiently large . we note that the model space with @xmath253 corresponds to that used in the @xmath13-matrix calculations made by halderson @xcite and by hao and kuo @xcite . .dependence of typical diagonal matrix elements @xmath254 on @xmath255 in units of mev , where @xmath256 and @xmath257 denote the @xmath32 and @xmath19 single - particle orbits , respectively , labeled as @xmath258 , @xmath259 and @xmath260 , and @xmath223 is the total angular momentum . [ cols="^,^,^,^,^,^,^,^ " , ]     next , the partial - wave contributions to the first - order potential energy of @xmath32 in the 0@xmath1 state are given for the different treatments of the pauli - blocking effect . we see from table  4 that the contributions of the @xmath261 and @xmath262-@xmath263 channels are very different dependently on the @xmath5 interactions employed . in the @xmath261 channel the largest pauli - blocking effect is brought about in the result for the j potential . on the other hand , in the @xmath262-@xmath263 channel we have an extremely large pauli - blocking effect for the nsc potential . this result shows that the nsc potential has the strong @xmath7-@xmath8 coupling in the @xmath262-@xmath263 channel as discussed in ref . the absolute values contributing to the first - order potential energy of @xmath32 in the channels with angular momenta @xmath264 are much smaller than those in the @xmath261 and @xmath262-@xmath263 channels .    in table  5 we display the single - particle energies of @xmath32 with the contributions of the perturbation corrections given in fig .  3 , where the pauli - blocking effect is treated accordingly to case ( 3 ) . from this table we see that the j and nsc potentials yield almost the same values of the first - order potential energy pe and the perturbation corrections in spite of having different partial - wave contributions as shown in table  4 . the single - particle energies obtained for the jb potential are most attractive among the @xmath265 interactions used . it should be noted here that the values of the ke plus pe in the 0@xmath266 states for three kinds of @xmath5 interactions are positive ( unbound ) . therefore , the perturbation corrections play a decisive role when we argue that the @xmath32 in the 0@xmath266 states is bound or not . correction ( a ) in fig . 3 is induced by the presence of non - diagonal terms in the one - body hamiltonian @xmath267 . this effect should be treated in the hartree - fock ( hf ) calculation . however , in the present calculation , we treat this effect perturbatively . these non - diagonal terms work to change the single - particle wave function of @xmath32 . therefore , this procedure is indispensable in this approach since we have taken a fixed value of @xmath100 commonly for @xmath32 and @xmath19 . this correction  ( a ) contributes most attractively to the single - particle energies among all the corrections in fig . the magnitude of every correction is considerably small in comparison with the first - order potential energies for all the @xmath5 interactions used . all the perturbation corrections contribute , however , attractively to the single - particle energies of @xmath32 , and do not change the mutual spacings of the single - particle states significantly . this trend is very different from the situation observed in ordinary nuclei such as @xmath15o and @xmath16ca , in which the correction terms contribute to about 40@xmath268 of the calculated spin - orbit splitting of the 0@xmath266 or 0@xmath269 states @xcite . the role of many - body correlations , such as the core polarization , does not seem to be important in the spin - orbit splitting of @xmath32 . we have presented a formulation of the umoa to apply it to a structure calculation of hypernuclei . we have introduced a unitary transformation exp(@xmath270 ) , with the correlation operator @xmath270 , that describes two - body correlation of hyperon and nucleon . an equation has been given for determining the correlation operator @xmath270 as the equation of decoupling so that the unitarily transformed interaction does not have non - zero matrix elements between the model space @xmath88 of low - lying @xmath8 states and its complement @xmath85 . the effective interaction has been given by the projected interaction onto the @xmath88 space . the effective interaction thus determined has properties of being hermitian , energy - independent and decoupled between the @xmath88 and @xmath85 spaces . in comparison with the usual @xmath13-matrix approach , the use of such an effective interaction would have some advantages in describing many - body systems . the umoa has been applied to the calculation of the structure of @xmath0o . the @xmath8 effective interaction has been calculated by employing three @xmath12 potentials , namely , the j , jb and nsc potentials . with the @xmath8 effective interaction , we have calculated the single - particle energies of @xmath32 in the 0@xmath1 , 0@xmath2 and 0@xmath3 states . we have seen that how to treat the pauli - blocking effect is very important , especially when we solve the problem of the @xmath7-@xmath8 coupling . in the usual treatment one excludes all the interactions acting in a space of @xmath271 and @xmath272 states in which a nucleon is in occupied states . however , this treatment of the pauli - blocking effect would lead to overmuch counting . we have introduced a new treatment in which , while we exclude the interactions acting among forbidden @xmath5 states as in the usual treatment , we restore the interactions that are diagonal in occupied nucleon states . we have observed that this treatment of the pauli - blocking effect has given rise to a significant effect as compared with the usual treatment . the present results have shown that @xmath32 in the 0@xmath3 state is unbound for the j and nsc potentials and , on the other hand , bound for the jb potential . the spin - orbit splitting of the single - particle levels of @xmath83 in the 0@xmath266 states has been given by 1.41 , 1.57 and 1.66 for the j , jb and nsc potentials , respectively . these values are very large compared with the results in ref . it will be of high interest for the single - particle nature of @xmath32 in @xmath0o to be established experimentally .    in the application of the umoa some problems have still remained . the unitarily transformed hamiltonian @xcite contains originally three - or - more - body interactions . the effects of these many - body interactions have been evaluated in the previous works of calculating nuclear properties , and we have concluded that the contributions of these many - body terms are much smaller than those of two - body interactions but sizable in some cases . we therefore should estimate the effects of these many - body interactions also in hypernuclei . another refinement to be made is that , in the present calculation , the single - particle energies of @xmath33 in intermediate states are neglected . however , the single - particle spectrum of @xmath33 has a possibility of giving rise to some contributions to the @xmath8 effective interaction . although the properties of the single - particle energies of @xmath33 have not always been made clear , we should investigate how much the spectrum of @xmath33 affects in the determination of the @xmath8 effective interaction . this problem is also a remaining task in the umoa . the authors would like to thank y.  yamamoto and k.  miyagawa for useful discussions and their help for using the nijmegen and jlich @xmath265 potential codes . one of the authors ( s.  f. ) would like to express his thank to k.  takada , m.  kamimura and y.  r.  shimizu for their continuous encouragement and instructive discussions . a.  gal , advances in nuclear physics , vol .  8 , eds . m.  baranger and e.  vogt ( plenum , new york , 1975 ) p.  1 . b.  povh , ann . nuc .  part . 28 ( 1978 ) 1 . h.  band et al . , prog . 81 ( 1985 ) . h.  band , t.  motoba and j.  @xmath273ofka , int .  j.  mod . a  5 ( 1990 ) 4021 , and references therein . t.  motoba , y.  akaishi and k.  ikeda , prog . 117 ( 1994 ) . b.  f.  gibson and e.  v.  hungerford iii , phys . 257 ( 1995 ) 349 . b.  f.  gibson , i.  r.  afnan , j.  a.  carlson and d.  r.  lehman , prog . 117 ( 1994 ) 339 . w.  brckner et al . , phys . . 79  b ( 1976 ) 157 . r.  h.  dalitz , d.  h.  davis , t.  motoba and d.  n.  tovee , nucl . phys . a  625 ( 1997 ) 71 . k.  miyagawa , and w.  glckle , phys . rev . c  48 ( 1993 ) 2576 ; + k.  miyagawa , h.  kamada , w.  glckle and v.  stoks , phys . rev . c  51 ( 1995 ) 2905 . t.  yamada , h.  ohkura , t.  motoba and k.  ikeda , prog . 117 ( 1994 ) 65 . e.  hiyama , m.  kamimura , t.  motoba , t.  yamada and y.  yamamoto , prog . 97 ( 1997 ) 881 . h.  band and y.  yamamoto , prog . theor .  phys . 81 ( 1985 ) 9 . y.  yamamoto , t.  motoba , h.  himeno , k.  ikeda and s.  nagata , prog . 117 ( 1994 ) 361 . h.  band , prog . 66 ( 1981 ) 1349 . d.  halderson , phys . rev . c  48 ( 1993 ) 581 . j.  hao , t.  t.  s.  kuo , a.  reuber , k.  holinde , j.  speth and d.  j.  millener , phys . 71 ( 1993 ) 1498 ; + t.  t.  s.  kuo and j.  hao , prog . 117 ( 1994 ) 351 . m.  hjorth - jensen , a.  polls , a.  ramos and h.  mther , nucl . phys . a  605 ( 1996 ) 458 . i.  vidaa , a.  polls , a.  ramos and m.  hjorth - jensen , nucl . phys . a  644 ( 1998 ) 201 . t.  t.  s.  kuo and e.  osnes , in : folded - diagram theory of the effective interaction in nuclei , atoms and molecules , lecture notes in physics , vol . 264 , ( springer , berlin , 1990 ) . j.  hao and t.  t.  s.  kuo , phys . 264 ( 1996 ) 233 . k.  suzuki , r.  okamoto and h.  kumagai , phys . c  36 ( 1987 ) 804 ; + k.  suzuki and r.  okamoto , prog .  theor . ( 1994 ) 1045 . h.  kumagai , k.  suzuki and r.  okamoto , prog . 97 ( 1997 ) 1023 . s.  fujii , r.  okamoto and k.  suzuki , prog . 99 ( 1998 ) 151 . m.  m.  nagels , t.  a.  rijken and j.  j.  de  swart , phys . rev . d  15 ( 1977 ) 2547 ; + p. m.  m.  maessen , t.  a.  rijken and j.  j.  de  swart , phys . c  40 ( 1989 ) 2226 . b.  holzenkamp , k.  holinde and j.  speth , nucl . phys . a  500 ( 1989 ) 485 ; + a.  reuber , k.  holinde and j.  speth , nucl . phys . a  570 ( 1994 ) 543 . m.  lacombe , b.  loiseau , j.  m.  richard , r.  vinh  mau , j.  ct , p.  pirs and r.  de  tourreil , phys . c  21 ( 1980 ) 861 . b.  h.  brandow , in : effective interactions and operators in nuclei , lecture notes in physics , vol .  40 , ed . b.  r.  barrett ( springer , berlin , 1975 ) p.  1 . i.  shavitt and l.  t.  redmon , j. chem . 73 ( 1980 ) 5711 . p.  westhaus , int . j. quantum chem . 20 ( 1981 ) 1243 . k.  suzuki , prog . 68 ( 1982 ) 246 . t.  t.  s.  kuo , p.  j.  ellis , j.  hao , z.  li , k.  suzuki , r.  okamoto and h.  kumagai , nucl . phys . a  560 ( 1993 ) 621 . y.  yamamoto and h.  band , prog . 83 ( 1990 ) 254 . schulze , a.  lejeune , j.  cugnon , m.  baldo and u.  lombardo , phys . b  355 ( 1995 ) 21 ; phys . c  57 ( 1998 ) 704 . k.  suzuki and r.  okamoto , prog . 76 ( 1986 ) 127 . e.  m.  krenciglowa , c.  l.  kung , t.  t.  s.  kuo and e.  osnes , ann .  of  phys . 101 ( 1976 ) 154 . b.  buck and a.  c.  merchant , nucl . a  600 ( 1996 ) 387 . c.  w.  wong , nucl . phys . a  91 ( 1967 ) 399 . p.  u.  sauer , nucl . phys . a  150 ( 1970 ) 467 .
Abstract:

Answer: a method is formulated for the description of lambda hypernuclei in the framework of the unitary - model - operator approach ( umoa ) . the method is applied to @xmath0o . a lambda - nucleon effective interaction is derived , taking the coupling of the sigma - nucleon channel into account . the lambda single - particle energies are calculated for the 0@xmath1 , 0@xmath2 and 0@xmath3 states employing the nijmegen soft - core ( nsc ) , jlich model -  ( j ) and model - b ( jb ) hyperon - nucleon potentials .    unitary model operator approach ; hypernuclei ; lambda single particle energy ; effective interaction ; nijmegen potential ; jlich potential


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: gas cooling plays a key role in galaxy formation . according to the hierarchical structure formation model , galaxies form in gravitationally collapsed dark matter halos . the dissipative baryonic matter cools by emitting radiation and condenses into halos to form galaxies @xcite . subsequently stars form through the radiative cooling of interstellar medium ( ism ; e.g. , * ? ? ? the cooling rate of gas varies with its chemical abundance , because different atoms and molecules have different cooling rates . the chemical abundance of cosmic gas evolves as a function of time and environment . in the early universe , the primordial gas consists of mainly hydrogen ( mass fraction of @xmath3 ) and helium ( mass fraction of @xmath4 ) with a negligible amount of light metals ( e.g. * ? ? ? * ) . throughout this paper , all the elements heavier than helium are collectively called ` metals ' . most metals are produced inside stars and spread out to the space by galactic winds driven by supernova ( sn ) explosions . in the metal - enriched gas , the total cooling rate is significantly enhanced relative to that of the primordial gas by the atomic emission lines from recombination of ionised metals . metals may influence galaxy formation in two different ways . firstly , the sf efficiency is increased owing to the shorter gas cooling time . star formation generally takes place in dusty metal - enriched ism , therefore , the enhanced cooling rate by metals would boost up the sf efficiency . however , the energy / momentum feedback by sne may suppress the subsequent star formation after the initial starburst by heating up the ambient gas . this complicates the situation , and it is not clear whether the net effect of metal cooling and sn feedback would be negative or positive . the feedback process is highly nonlinear , therefore , a direct numerical simulation would be a useful tool to explore the effects of feedback and metal cooling . secondly , the metals dispersed into the igm by sn feedback also enhance the cooling of igm . this can lead to the increase of igm accretion onto galaxies , because colder gas can sink into galaxies more easily . both of the above two effects may affect the galaxy growth and cosmic sf history significantly . therefore it is essential to include the effects of metal cooling and chemical enrichment by feedback in the studies of galaxy formation and evolution , and capture the complex two - way interactions between galaxies and igm . in addition , metal enrichment also affects the equation of state by altering the mean molecular weight of gas , which influences the hydrodynamic calculation of gas thermal state . cosmological hydrodynamic simulations are widely used in the studies of galaxy formation and cosmic sf history ( e.g. , * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? most of these works included the treatments of radiative cooling by h & he , star formation , and sn feedback . however , some of the simulations did not include the effects of metal cooling and/or chemical enrichment by galactic wind . furthermore , the effects of metal cooling on galaxy growth and cosmic sf history has not been explored in detail using cosmological hydrodynamic simulations and presented in the literature , as it is costly to run large cosmological simulations with and without the effect of metal cooling .    in this paper , we investigate the effects of metal cooling and metal enrichment on galaxy growth and cosmic sf history using a series of cosmological hydrodynamic simulations with and without metal cooling . the aim of this paper is to single out the effects of metal cooling among our simulations . this work can be regarded as our initial step towards the long - term goal of developing more complete cosmological hydrodynamic code with physically motivated models of star formation and feedback . the paper is organised as follows . in  [ sec : method ] , we describe our simulation method focusing on how we implement the metal cooling . in   [ sec : global ] , we study the metal cooling effects on the global properties , such as cosmic sfr , phase space ( @xmath5 vs. @xmath6 ) distribution of gas , and the evolution of four phases ( hot , warm - hot , diffuse , and condensed ) of baryons . we then study the galaxy mass functions (   [ sec : mf ] ) and gas mass fractions (   [ sec : gasfrac ] ) . finally , we discuss and summarise our findings in   [ sec : summary ] . we use the updated version of the tree - particle - mesh ( treepm ) smoothed particle hydrodynamics ( sph ) code gadget-2 @xcite for our cosmological simulations . the gravitational dynamics is computed by a treepm algorithm , which uses a particle - mesh method @xcite for the long - range gravitational force and a tree method for the short - range gravitational force @xcite . this hybrid algorithm makes the gravitational force calculation faster than a tree method and allows better force resolution than a pm method in dense regions . the gas dynamics is computed by an sph method . the sph is particularly useful if the simulation needs to resolve large dynamical range , which is an inevitable requirement for the study of galaxy formation in a cosmological context . therefore a treepm - sph simulation can provide a fast and high resolution calculation for both gravitational dynamics and hydrodynamics . the gadget-2 code adopts the entropy - conservative formulation @xcite , which alleviates the overcooling problem that previous sph codes suffered from . our basic simulations include radiative cooling and heating processes for hydrogen and helium using a method similar to @xcite . an external uv background radiation is treated as a spatially uniform photoionising radiation @xcite , and modified to match the ly@xmath7 forest observations @xcite . implementing star formation and feedback from first principles is not feasible in current cosmological simulations , because the spatial & mass - scales of molecular clouds are not resolved . star formation and sn feedback are represented by the subgrid multiphase ism model developed by @xcite . in the multiphase scheme , a single gas particle represents both hot and cold gas . the stars are formed in the cold portion when the density exceeds a given threshold , @xmath8 , which is derived self - consistently within the multiphase ism model . the related parameters , such as the normalisation of gas consumption time and evaporation efficiency of cold gas , are set to satisfy the empirical kennicutt - schmidt law @xcite . the sn feedback returns some fraction of the cold gas to the hot phase , and increases the thermal energy of the hot gas .    as an extension to the multiphase ism model , the simulation includes a phenomenological model for sn - driven galactic wind @xcite . the galactic wind is particularly important for distributing the metals produced by sne into the igm . we use the strong kinematic wind with a velocity of 484kms@xmath9 . it has been shown that this model produces favourable results for the luminosity function of lyman - break galaxies at the bright - end @xcite and the hi column density distribution function @xcite at @xmath1 , when compared to the runs without the wind . however , @xcite pointed out the problems of this galactic wind model by comparing with the observations of civ absorption lines in quasar spectra , and suggested that the momentum - driven wind model @xcite is a more viable model , which can carry more metals with lower wind velocities . in this paper we choose not to modify our galactic wind model in order to single out the effects of metal cooling , and to allow direct comparisons to the previous works @xcite . here we focus on the effects of metal cooling on galaxy growth , while @xcite focused on the civ statistics of the igm . the metallicity of gas particles are also tracked by the code , assuming a closed box model for each gas particle . the yield ( @xmath10 ) of 0.02 is assumed . in principle , there could be a time delay between sn explosions and chemical enrichment of the ambient gas . unfortunately , current cosmological simulations do not have sufficient resolution to track the detailed mixing process of metals . in our simulations , we ignore this time delay and assume an instantaneous mixing within each gas particle . this assumption is reasonable , because the overall time step of the simulations is longer than the mixing time - scale on small scales . ideally we would like to track individual metals and compute the cooling rate of each element at each time - step of the simulation . this approach has been used to study the evolution of individual galaxies @xcite . for cosmological simulations , it requires a large memory to track individual metals and the computation becomes exceedingly expensive ( but see * ? ? ? * ; * ? ? ? * for such an attempts ) . our simulations collectively track the total mass fraction of metals ( @xmath11 ) as the measure of metallicity . we obtain the chemical abundance and cooling rates for a given @xmath12 and temperature by interpolating the table given in ( * ? ? ? * hereafter sd93 ) for the standard collisional ionisation equilibrium model , and create a lookup table . metal cooling is implemented in the full range of metallicity in the table from sd93 , from [ fe / h]=-3 to [ fe / h]=1 . in addition , the metal cooling outside of this range is computed by extrapolating the table . we require that the extrapolation does not go lower than the primordial cooling rate . the primordial abundance pattern is used where [ fe / h ] @xmath13 @xcite , and the solar abundance pattern is used where [ fe / h ] @xmath14 @xcite . for @xmath15 [ fe / h ] @xmath16 , the abundance pattern is computed by interpolating between the primordial and solar abundance patterns . note that , [ fe / h ] = @xmath17 - @xmath18 , where @xmath19 is the number density . figure [ fig : xyz_table ] shows the relationships between the metal parameters ( [ fe / h ] , @xmath20 , @xmath21 , and @xmath22 ) and the total metallicity @xmath12 . because the abundance pattern is fixed for a given @xmath12 , [ fe / h ] is a monotonically increasing function of @xmath12 . the main difference between the solar and the primordial abundance pattern is the oxygen or @xmath7-particle enhancement . in the early phase of metal enrichment , type ii sne are the main source of metals and they tend to generate more @xmath7-particles than sn ia . therefore , it is expected that a low metallicity gas to have more @xmath7-particles than the high metallicity one . as mentioned above , we use two different chemical abundance pattern ( at low- and high - end of the @xmath12 values ) to model this evolution in chemical composition . we compute the oxygen abundance based on the table 4 of sd93 . in our code , we incorporate the metal cooling effect by adding the additional cooling contribution from metals on top of the primordial cooling rate , which is computed in a similar fashion as described in @xcite . metal enrichment changes the mean molecular weight @xmath23 , which is needed to compute the temperature and internal energy of gas . in order to compute @xmath23 , we need to know the electron number density @xmath24 , which varies with metallicity . we obtain the values of @xmath24 from sd93 in the same way as we obtained the metal cooling rates . our electron number density estimate ignores the electrons from photoionised metals by the uv background radiation , but we expect that the number of electrons from ionised metals will be small .    varying the mean molecular weight @xmath23 also influences the sf threshold density @xmath8 in the multiphase ism model of @xcite , as well as the cold and hot gas fractions of multiphase gas particles . in this paper , we consider two different treatments of @xmath8 . in the ` constant @xmath8 ' model , we take the value of @xmath8 as in the original formulation by @xcite and do not modulate @xmath8 by the varying @xmath23 . in the ` varying @xmath8 ' model , we take into account of the modulation of @xmath8 by the change in @xmath23 . we use three series of simulations with varying box size and resolution . the effects of metal enrichment and metal cooling are tested by comparing the runs listed in table  [ tbl : simulation ] . the effects of galactic wind are also tested by turning on and off the wind . the adopted cosmological parameters of all simulations are @xmath25 , which mostly concur with the results of the wilkinson microwave anisotropy probe ( wmap1 ) @xcite .    [ cols="<,^,^,^,^,^,^,<,^",options="header " , ]     our naming conventions for the simulations are as follows . the first part of the run name denotes the particle number and box size used for the simulation series . the run with the extension of ` mc ' ( for ` metal cooling ' ) adopts the ` constant @xmath8 ' model for metal cooling , the run with the extension ` mv ' ( for ` metal , varying ' ) adopts the ` varying @xmath8 ' model , the run with the extension ` nw ' has no galactic wind , and the run without any extension uses the cooling rates for the primordial chemical composition ( i.e. , only h and he ) . the ` n216l10 ' series is the fiducial simulation set in our study . we implement one simulation without galactic wind ( ` n216l10nw ' run ) . the ` n216l10 ' run has the same physical models as the ` q4 ' run used in @xcite and @xcite with no metal cooling .    to examine the resolution effect , we implement lower resolution simulations , the ` n144l10 ' series . owing to the comparably small box size and missing long wavelength perturbations , we evolve the n144l10 and n216l10 series only down to @xmath26 . to evolve the simulation toward lower redshifts , we need a simulation with a larger box size . therefore we implement the ` n288l34 ' series , which has a larger box size than other series , but a lower resolution . this series is evolved down to @xmath0 . the runs within the same series use the identical initial condition , which makes the comparison more robust and enables halo - by - halo comparison if needed .      in cosmological simulations , dark matter , gas , and stars are represented by particles , and these particles are monte carlo representations of matter distribution . in order to study galaxy formation using a cosmological simulation , we need a working definition for simulated galaxies . in this paper the simulated galaxies are defined as isolated groups of star and gas particles , identified by a simplified variant of the algorithm developed by @xcite . in more detail , the code first computes a smoothed baryonic density field to identify candidate galaxies with high density peaks . the full extent of these groups are found by adding gas and star particles to the groups in the order of declining density . if all @xmath27 is the minimum number of gas and star particles that constitute one isolated group . in this paper we set @xmath28 . ] nearest neighbour particles have lower densities , this group of particles is considered as a new group . if there is a denser neighbour , the particle is attached to the group to which its nearest denser neighbour already belongs to . if two nearest neighbours belong to different groups and one of them has less than @xmath27 particles , these the two groups are merged . if two nearest neighbours belong to different groups and both of them has more than @xmath27 particles , the particles are attached to the larger group , leaving the other group intact . in addition , the gas particles in groups should be denser than @xmath29 . in this paper we are not concerned with groups without star particles . the fraction of such groups are very small , therefore they do not affect our conclusions . % by metal cooling throughout the entire redshift range . ]    , but for the n288l34 series , which goes down to @xmath0 . the n216l10mc run is also shown for a resolution comparison . the green dot - dashed line ( hs03 ) is for the analytic model of ( * ? ? ? * eq.(2 ) ) . ] figure  [ fig : sfr_n216l10 ] shows the effect of metal cooling and galactic wind on the cosmic sf history for the n216l10 series . the comparison of n216l10mc and n216l10 runs clearly shows that metal cooling enhances the cosmic sfr throughout all redshifts by @xmath30% . previously , ( * ? ? ? 13 ) argued that the sfr density is enhanced by metal cooling mostly at lower redshifts and hardly at high-@xmath31 , however , our simulations suggest otherwise . ( see   [ sec : summary ] for more discussion on this point . ) metallicity needs to be higher than @xmath32 to noticeably increase the gas cooling rate , because the metal cooling rate is much smaller than that by h and he at @xmath33 . the metallicity of diffuse igm is generally lower than @xmath34 at @xmath35 , and the metal cooling effect on igm is small . however the local sfr could be enhanced significantly even at high-@xmath31 when the ism is enriched to @xmath36 , which can be easily achieved after a few events of sn ii . our results suggest that star formation at high-@xmath31 can also be enhanced by metal cooling , and not just at low-@xmath31 . we also find that the peak of the cosmic sf history hardly shifts by the introduction of metal cooling , consistently with the estimate by @xcite . the n216l10mv run has a lower sfr at @xmath37 and a higher sfr at @xmath38 than the n216l10mc run . this is because at high-@xmath31 , the value of @xmath8 is higher in the ` mv ' run than in the ` mc ' run due to lower metallicity and mean molecular weight , leading to a larger fraction of gas not being able to satisfy the sf criteria in the n216l10mv run than in the n216l10mc . the opposite is true at @xmath39 . figure  [ fig : sfr_n216l10 ] also shows that the sfr in the n216l10nw run is higher than other runs by an order of magnitude . this is because the gas can cool very efficiently without being ejected by the galactic wind . although the current galactic wind models in cosmological simulations are not well established yet and need to be improved , both theoretical and observational studies evidently suggest the important role of galactic wind feedback @xcite . as expected , the n144l10mc run has lower sfr densities at @xmath40 compared to the n216l10 series . this is because a lower resolution run is unable to resolve low - mass halos that would otherwise forms stars in them . when the star formation at high-@xmath31 is underestimated , some of the gas remain unconverted into stars until lower redshift , resulting in a higher sfr at @xmath41 in the n144l10mc run than in the n216l10mc run . figure  [ fig : sfr_n288l34 ] shows the cosmic sf history down to @xmath0 for the n288l34 series . again , the sfr density in the n288l34mc run is systematically higher than that in the n288l34 run at all redshifts . moreover , the difference in sfr density between the two runs increases at @xmath39 . by @xmath42 the igm metallicity becomes high enough to noticeably increase the cooling rate , and the accretion of igm onto galaxies becomes more efficient , which later fuels the star formation . our results suggest that the cosmic sfr density is enhanced by metal cooling at @xmath43 by both of the following two effects : 1 ) more efficient igm accretion onto galaxies due to higher metallicity of igm , and 2 ) increased sf efficiency due to higher cooling rate with the contribution from metal cooling .    in figure [ fig : sfr_n288l34 ] the analytic model of ( * ? ? ? * eq.(2 ) , hereafter h&s model ) is also shown in the green dot - dashed line . this model is a result of many cosmological sph simulations with increasing resolution , and the authors gave a theoretical interpretation to the empirical fitting formula for the cosmic sfr density . since our runs adopt a lower value of @xmath44 than the original value used to derive the h&s model ( @xmath45 ) , we modify the h&s model by rescaling the @xmath46 parameter in their eq.(2 ) and eq.(47 ) . this reduces the sfr density at @xmath47 by @xmath48% , but not much at @xmath49 . the results of the n288l34 and n216l10mc runs agree well with the h&s model at @xmath41 and @xmath50 , respectively . this is expected , because the n288l34 run adopts the same physical models as in the original simulations that were used to derive the h&s model . at high-@xmath31 , both runs fall short compared to the h&s model owing to insufficient resolution . the metal enrichment and metal cooling also change the phase space distribution ( @xmath51 diagram ) of cosmic gas . figure  [ fig : phase.n216l10 ] shows the state of cosmic gas at @xmath1 for the n216l10 series . as shown in panel ( @xmath52 ) , the baryons in the universe can be broadly categorised into four different phases according to their overdensity and temperature : ` hot ' , ` warm - hot ' , ` diffuse ' , and ` condensed ' @xcite . the first two phases are mostly the shock - heated gas in clusters and groups of galaxies . the ` diffuse ' phase is mostly the photoionised igm with lower temperature , which can be observed as the ly@xmath7 forest . the tight power - law relationship between @xmath5 and @xmath6 at @xmath53 and @xmath54k is governed by the ionisation equilibrium @xcite , where @xmath55 is the mean density of baryons at @xmath56 . the ` condensed ' phase is the high density , cold gas in galaxies . because the primordial cooling curve has a sharp cutoff at @xmath57k , the temperature of condensed phase decreases only down to this temperature in figure  [ fig : phase.n216l10]a , b . the characteristic finger - like feature at @xmath58 is due to the pressurisation of the star - forming gas by sn feedback in the multiphase ism model of @xcite . a fitting formula for this effective equation of state was derived by @xcite , which is shown by the blue dashed lines . above the sf threshold density @xmath8 , the gas particles in the simulation go into the multiphase mode and are allowed to form stars . the most noticeable and interesting change in the @xmath51 plot by the metal cooling is in the distribution of ` condensed ' gas . star formation and sn explosions occur in the multiphase star - forming gas , therefore the immediate effect of metal cooling appears in the condensed phase . the presence of gas at @xmath59k with metal cooling results from the effective temperature calculation for the multiphase medium , and not from the direct gas cooling down to @xmath59k . according to the multiphase scheme in @xcite , the star forming gas consists of two phases : hot gas ( @xmath60k @xmath61 @xmath62k ) and cold gas ( @xmath63k ) . the effective temperature for the multiphase gas is computed from the internal energy of hot and cold gases , which are derived using the cold gas fraction . the cold gas fraction is computed from the cooling rate of hot gas and the gas density . the original gadget imposes the effective temperature of @xmath64k at the sf threshold density . in our metal cooling scheme , however , the cooling rate of hot gas increases , and the cold gas fraction increases . therefore , the multiphase gas above the sf threshold density with higher metallicity could have lower effective temperature of @xmath59k in our implementation .    in the n216l10mc run , the value of @xmath8 was kept fixed as the original value , which causes a sharp edge at @xmath65 for the multiphase gas . in the n216l10mv run , the value of @xmath8 was varied according to the change in mean molecular weight , therefore the condensed gas has a spread in both temperature and density near @xmath65 . galactic wind returns the metal - enriched , star - forming gas to the hot and low density igm , as evidenced in the @xmath66 diagrams . this metal - enriched , low - density ( @xmath67 ) gas is absent in the n216l10nw run , which clearly shows the necessity of galactic wind to explain the metallicity observed in ly@xmath7 forest . the phase distribution of igm in the n216l10 , n216l10mc , and n216l10mv runs are very similar . this is presumably because the effects of metal cooling is not so significant yet in the igm at @xmath68 . as we will discuss in the next section , this is not be the case at @xmath2 . the connection between the four phases of cosmic gas is important in galaxy formation , and it is beneficial for us to study their evolution for a better understanding of metal cooling effects . a conventional view is that the diffuse gas is shock - heated to hot or warm - hot phase , which later cools down to the condensed phase @xcite . @xcite also suggested that some diffuse gas can migrate to the condensed phase without being shock - heated , which they called the ` cold mode ' accretion . figure  [ fig : evolv4 ] shows the evolution of mass fractions of different phases . all of our simulations show that the diffuse phase continuously decreases from high-@xmath31 to low-@xmath31 , while the ` warm - hot + hot ' and ` condensed ' phases continuously increase , consistently with the general expectations . the igm ( i.e. , hot + warm - hot + diffuse ) is accreted onto galaxies , to become the condensed phase . the diffuse component ( panel @xmath52 ) decreases from @xmath69% at @xmath70 to @xmath71% at @xmath0 . the ` warm - hot + hot ' component ( panel @xmath72 ) increases from almost zero at @xmath70 to @xmath73% at @xmath0 . the mass fraction of the total igm ( panel @xmath74 ) decreases from @xmath75% at @xmath70 to @xmath76% at @xmath1 . this reflects the increase of the condensed phase ( panel @xmath77 ) from a few percent at @xmath70 to @xmath78% at @xmath1 . these mass fractions are in general agreement with the results of @xcite . the comparison between n288l34 and n288l34mc runs shows interesting differences . at @xmath68 , the mass fractions of all four phases are similar in the two runs . however , at @xmath2 , the igm mass fraction continues to decrease in the n288l10mc run , while it becomes almost constant in the n288l10 run . a similar behaviour is seen for the condensed phase , and the n288l10mc continues to increase its condensed mass . this suggests that the metal cooling enhances the igm accretion onto galaxies mostly at @xmath39 , which supplies the fuel for star formation . the cooling of warm+hot phase into diffuse phase is also enhanced , leading to a larger mass fraction of diffuse component in the n288l34mc run than in the n288l34 run . ( the same is true for n216l10mc and n216l10 runs . ) the n216l10nw run shows stark differences from all other runs . it has a significantly lower igm mass fraction , and a very large mass fraction in the condensed phase . this shows that the igm is mostly enriched by galactic wind , and that most of the masses are returned into the ` warm - hot + hot ' phase , but not to the diffuse phase . galaxy stellar mass function ( gsmf ) tells us how the stellar mass is distributed in different galaxies , and its redshift evolution reflects the growth of structure in the hierarchical universe . figure  [ fig : mf_star ] compares the gsmf in different simulations . we only show the range of @xmath80 , which corresponds to the limiting mass of 32 star particles in n216l10 runs . we chose 32 star particles as our resolution limits because all mass function start to turn over around this mass range . panel  ( @xmath52 ) shows that the run with no galactic wind ( n216l10nw ) completely overestimates the gsmf , and the run with lower resolution ( n144l10 ) underestimates the number of lower mass galaxies at @xmath81 relative to the n216l10 run . panel  ( @xmath72 ) shows that the metal cooling increases the masses of galaxies by @xmath82% at @xmath1 . as a result , both the n216l10mc and n216l10mv runs show a slight increase in the number of galaxies at the massive end , with a slightly stronger enhancement in the n216l10mv run .    figures  [ fig : mf_star]c , d demonstrate that the metal cooling also changes the redshift evolution of gsmf . in the n216l10mc run , the peak of mass function is shifting more toward the massive side than in the n216l10 run at @xmath1 . from @xmath83 to @xmath1 , the low - mass galaxies merge to form more massive galaxies , increasing the number of galaxies with @xmath84 furthermore . it also seems that the formation of low - mass galaxies ( @xmath81 ) is enhanced by metal cooling at early times ( @xmath85 ) . however we find that the behaviour of mf at the low - mass end is somewhat dependent on the threshold density for the grouping . as mentioned in  [ sec : groups ] , our grouping code imposes a threshold gas density of @xmath29 for a gas particle to be part of a galaxy . when we lower this threshold density , we find that the two mfs for the n216l10 and n216l10mc converged , as more particles are incorporated into galaxies at the outskirts of galaxies . metal cooling increases the gas density in galaxies , enabling more low - mass galaxies to satisfy the gas density threshold and to be identified as simulated galaxies . therefore the behaviour of mf at @xmath86 is somewhat dependent on the choice of the outer threshold density of galaxies . figure  [ fig : mf_comp ] compares the gas and baryonic ( star + gas ) mfs for the n216l10 series at @xmath1 ( panels @xmath52 & @xmath72 ) , and the redshift evolution of the baryonic mf from @xmath87 to @xmath1 in the n216l10 ( panel @xmath74 ) and n216l10mc ( panel @xmath77 ) runs . here we show only the range of @xmath88 , which corresponds to the limiting mass of 32 gas particles . panel ( @xmath52 ) shows that more gas is converted into stars in the runs with metal cooling ( n216l10mc and n216l10mv ) compared to the n216l10 run , while panel ( @xmath72 ) shows that the baryonic mf is similar in all the runs . as we discussed in  [ sec : evolv4 ] , the igm accretion onto galaxies is not so much enhanced yet before @xmath89 due to relatively low igm metallicity , which explains the similarity in the total baryonic mf . in the runs with metal cooling , the peak of gsmf is shifted toward higher mass , whereas the peak of gas mf is shifted toward lower mass . this suggests that the enhanced gsmf at @xmath1 is due to the increased sf efficiency by metal cooling within the galaxies . again , the apparent enhancement in the number of low - mass galaxies in the n216l10mc run compared to the n216l10 run is somewhat dependent on the threshold density of grouping , therefore it should be interpreted with caution .          as the accretion of igm onto galaxies become more efficient at @xmath2 due to the increased igm cooling by the metals , we expect that the total baryonic mass of galaxies would be more enhanced at @xmath0 than at @xmath1 . figure  [ fig : mf_comp_n288l34 ] shows the mfs at @xmath1 and @xmath0 in the n288l34 series , and it clearly demonstrates that this expectation is true . here we show only the range of @xmath91 , which corresponds to the limiting mass of 32 gas particles . the general trend in the three ( star , gas , and total baryon ) mfs is similar to that we saw in figure  [ fig : mf_comp ] , although the difference at @xmath1 between the n288l34 and n288l34mc runs is slightly smaller than in the n216l10 series due to poorer resolution . figure  [ fig : mf_comp_n288l34]f shows most prominently the enhancement of the total baryonic mf at @xmath0 in the n288l34mc run owing to the metal cooling . in panels ( @xmath92 ) & ( @xmath93 ) , the n288l34 run actually has a longer tail at the most massive - end than the n288l34mc run . the reason for this feature is not fully clear , but it may be related to the balance between igm accretion and feedback . owing to metal cooling , igm accretion rate increases in the mc run , and the sfr is also enhanced , leading to a stronger feedback . the amount of mass loss is greater in low mass galaxies , but the net heating of igm is more significant for massive galaxies . the mc run has stronger feedback , therefore its feedback heating may become more significant than igm accretion for very massive galaxies . the significant igm heating suppresses the growth of massive - end of mass function from @xmath1 to @xmath0 and results in shorter tail for the mc run at @xmath0 as shown in figure  [ fig : mf_comp_n288l34]e , f . for a fixed galaxy baryonic mass , the enhancement of star formation by metal cooling would decrease the gas mass fraction , @xmath94 . figure  [ fig : gfrac_evol_n216l10 ] shows @xmath95 as a function of galaxy stellar mass for the n216l10 series . panel ( @xmath52 ) shows the data only from n216l10mc run at @xmath1 , and each data point corresponds to a simulated galaxy . to characterise the distribution , we compute the following two quantities in each logarithmic stellar mass bin : ` average ' and ` median ' . the ` average ' is the ratio of total gas mass to total baryonic mass for all the galaxies in each mass bin , i.e. , @xmath96 . the ` median ' case is simply the median of @xmath95 values in each mass bin . both quantities show a similar trend , however , there is a sharper drop - off at @xmath97 for the ` median ' case in figure  [ fig : gfrac_evol_n216l10]a . this mass - scale corresponds to 32 star particles in the n216l10 series . we find that there are many galaxies with @xmath98 above this mass - scale , which causes the sharp drop - off in the ` median ' line . below this limiting mass , galaxies are not resolved well , which results in an underestimate of star formation and an overestimate of @xmath95 . if we had a higher resolution simulation with finer particle masses , this limiting mass - scale would shift to a lower mass . therefore the location of this sharp drop - off is currently determined by the resolution of our simulation . however , dark matter halos would stop forming stars at some lower limiting halo mass , if we had an infinitely high - resolution simulation , this lower limit to the galaxy mass is presumably determined by the photoevaporation of gas by the uv background radiation @xcite . recent works suggest that star formation could be suppressed by the uv background in halos with @xmath99 at @xmath89 . galaxies with @xmath100 would reside in halos with @xmath101 . the n216l10 series would resolve such a halo with @xmath102 dark matter particles , and its mass resolution is actually close to the astrophysical limit for dwarf galaxy formation at @xmath89 . therefore the sharp drop - off in @xmath95 at @xmath100 may not be so far from the true answer . we find that @xmath103 increases with decreasing @xmath104 at all redshifts , regardless of metal cooling and wind effects . this trend is qualitatively consistent with current observations . @xcite estimate the gas fraction as a function of stellar mass using the rest - frame uv - selected star - forming galaxies at @xmath105 , and show that the gas fraction decreases with increasing stellar mass . and using the mean gas and stellar mass , they find the average @xmath106 , which agrees with the predicted gas fraction for massive galaxies in our simulation . in addition , @xcite reported that the average neutral gas fraction is @xmath107 for the local dwarf galaxies selected from the sloan digital sky survey . in our n216l10 series with metal cooling , the ` average ' @xmath95 reaches 0.6 for galaxies with @xmath108 at @xmath1 & 4 . figures[fig : gfrac_evol_n216l10]b , c , d show the redshift evolution of @xmath95 for the n216l10 series . in the runs with metal cooling ( n216l10mc and n216l10mv ) , @xmath95 is lower than in the n216l10 run by @xmath30% at all redshifts for galaxies with @xmath109 . this result suggests that the metal cooling reduces @xmath95 owing to more efficient star formation . the values of @xmath95 seem to be more convergent at the massive - end ( @xmath110 ) . in addition , we find that galaxies in the n216l10nw run are the most gas - rich , because almost no galactic gas is returned to the igm .    figures  [ fig : gfrac_evol_n216l10]b , c , d also show that @xmath95 is higher at higher redshifts . for example , the well - resolved galaxies with @xmath111 have @xmath112 at @xmath87 , but @xmath113 at @xmath1 . as we saw in figure  [ fig : mf_star]d , there is a considerable growth in gsmf from @xmath87 to @xmath1 . together with the decrease in the gas fraction , these results suggest that most of the gas accreted during @xmath114 has been converted into stars .    figure  [ fig : mgfrac_resolution]a highlights the redshift evolution of @xmath95 from @xmath87 to @xmath0 in the n288l34mc run . the gas mass fraction clearly decreases with decreasing redshift as the gas is converted into stars . the rate of decrease is greater at @xmath115 than at @xmath116 , with @xmath117 & 0.2 for @xmath118 , & 1 , respectively , for galaxies with @xmath119 . figure  [ fig : mgfrac_resolution]b shows the effect of resolution on @xmath95 . the location of the drop - off in @xmath95 shifts to lower masses as the resolution is increased from n144l10mc to n216l10mc run . the location of the drop - off is at around our resolution limit  the vertical lines in _ panel ( b)_. therefore the values of @xmath95 are overestimated in unresolved galaxies . using cosmological hydrodynamic simulations with metal enrichment and metal cooling , we studied their effects on galaxy growth and cosmic sfr . owing to metal cooling , the sfr density increases about 20% at @xmath1 and about 50% at @xmath0 . our results suggest that metal cooling enhances the star formation through two different processes : 1 ) more efficient conversion of local ism into stars ( i.e. , increase of local sf efficiency ) , and 2 ) the increase of igm accretion onto galaxies . the former process is in effect essentially at all times at @xmath120 , because the local ism can be instantaneously enriched by sn explosions . this process enhances the sfr , but does not noticeably change the total baryonic mass of galaxies . the latter process , on the other hand , can increase the total baryonic mass of galaxies , as well as enhancing the overall sfr density by supplying more gas for star formation . this process becomes effective only at lower redshifts ( @xmath39 ) , because it takes some time to enrich the igm up to @xmath121 .    in this paper , we used two different prescriptions for star formation in the runs with metal cooling : 1 ) the ` constant @xmath8 ' scheme and 2 ) the ` varying @xmath8 ' scheme . the value of @xmath8 was fixed in the former scheme , while it was modulated according to the metallicity in the latter scheme . both schemes show similar increase of cosmic sfr density ( see figure  [ fig : sfr_n216l10 ] ) , which implies that the overall sf enhancement by metal cooling is not so sensitive to the choice of @xmath8 ( but with some differences as we discussed in  [ sec : cos_sfr ] ) . the multiphase ism model for star formation by @xcite contains two free parameters : @xmath8 and the normalisation of gas consumption time - scale , @xmath122 . the values of these two parameters were originally chosen to match the empirical kennicutt - schmidt law @xcite using simulations of isolated disk galaxies . since we did not change the value of @xmath122 in the runs with metal cooling ( as well as @xmath8 in the ` mc ' runs ) , it is possible that our simulations may now violate the kennicutt law . however , as we show in figure  [ fig : kenni ] , the plots of @xmath123 vs. @xmath124 for the n216l10 and n216l10mc runs are not so different from each other . this can be understood as follows . as figure  [ fig : evolv4]c showed , metal cooling increases the density and lowers the temperature of star - forming gas . the fundamental scaling relationship between sfr and cold gas density does not need to change when the metal cooling is introduced ; the star - forming gas particles would simply slide upward along the kennicutt law , @xmath125 .    in the case of n216l10mv run , we varied @xmath8 according to the gas metallicity , primarily scattering it to lower values because higher metallicity increases the cooling rate and lowers @xmath8 . therefore more gas particles become eligible to form stars and @xmath124 is scattered upward above the kennicutt law for a given value of @xmath123 , resulting in a broader deviation from the power - law relationship of the kennicutt law near @xmath126 . at higher values of @xmath123 , @xmath124 is enhanced , but still within the range of the kennicutt law . this difference can be alleviated by adjusting the sf timescale ( see * ? ? ? * ) , however , we try to emphasise the effect of metal cooing on sfr by keeping this parameter the same in all the runs . further investigations of the hi  aspect of our simulations is beyond the scope of the present paper , and we will present the results elsewhere .    upon evaluating the effect of metal cooling , @xcite considered two different gas phases where the cooling become important . one is the diffuse gas in galactic halos , which must radiate its thermal energy in order to collapse onto the high - density ism , and the other is the multiphase star - forming gas . they argued in their  5.2 that the cosmic sfr density at @xmath127 would not be affected by the metal cooling very much , because at high-@xmath31 cooling is so efficient that the gas in diffuse atmospheres of halos cools nearly instantly , even without any metal cooling . they expected their model results to be largely independent of metal enrichment , because in this regime the evolution of sfr density is driven by the fast gravitational growth of the halo mass function . for the star - forming multiphase ism , they argued that the parameters of their model ( e.g. cold gas evaporation efficiency and gas consumption time - scale ) can be adjusted such that the normalisation of the kennicutt law can be maintained , yielding to first order unaltered model predictions . however , our simulations , in particular the comparison between the n216l10 and n216l10mc run , suggest that the metal cooling can enhance the sfr density even at @xmath127 by @xmath128% . this is because the gas density becomes higher and temperature becomes lower ( see figure  [ fig : evolv4]c ) by responding to the enhanced cooling rate by metals .    in summary , we consider that our study serves to demonstrate the importance of metal enrichment and metal cooling in galaxy formation , beyond our original motivation to perform a simple numerical study . our study highlights the role of metals at low-@xmath31 and high-@xmath31 to enhance the star formation in the universe , and demonstrates that the metals increase the sf efficiency in the star - forming ism , as well as enabling more igm to accrete onto galaxies and fuel star formation . we expect that metal cooling would also alter some of the galaxy properties such as the mass - metallicity relationship ( e.g. , * ? ? ? * ) and specific sfr ( e.g. , * ? ? ? we plan to investigate these issues in the future , as well as alternative models for star formation and feedback . we thank v. springel for allowing us to use the updated version of gadget-2 code for our study , and for useful comments on the manuscript . kn is grateful for the hospitality of institute for the physics and mathematics of the universe , university of tokyo , where part of this work was done . we are also grateful to the anonymous referee for constructive comments which improved this paper . this research was supported in part by the national aeronautics and space administration under grant / cooperative agreement no . nnx08ae57a issued by the nevada nasa epscor program , by the national science foundation through teragrid resources provided by the san diego supercomputer center ( sdsc ) , and by the president s infrastructure award at unlv . the simulations were performed at the unlv cosmology computing cluster and the datastar at sdsc . d.  n. , verde l. , peiris h.  v. , komatsu e. , nolta m.  r. , bennett c.  l. , halpern m. , hinshaw g. , jarosik n. , kogut a. , limon m. , meyer s.  s. , page l. , tucker g.  s. , weiland j.  l. , wollack e. , wright e.  l. , 2003 , apjs , 148 , 175
Abstract:

Answer: we present the results of a numerical study on the effects of metal enrichment and metal cooling on galaxy formation and cosmic star formation ( sf ) history using cosmological hydrodynamic simulations . we find following differences in the simulation with metal cooling when compared to the run without it : ( 1 ) the cosmic star formation rate ( sfr ) is enhanced by about 50 & 20% at @xmath0 & 3 , respectively ; ( 2 ) the gas mass fraction in galaxies is lower ; ( 3 ) the total baryonic mass function ( gas + star ) at @xmath1 does not differ significantly , but shows an increase in the number of relatively massive galaxies at @xmath0 ; ( 4 ) the baryonic mass fraction of intergalactic medium ( igm ) is reduced at @xmath2 due to more efficient cooling and gas accretion onto galaxies . our results suggest that the metal cooling enhances the galaxy growth by two different mechanisms : ( 1 ) increase of sf efficiency in the local interstellar medium ( ism ) , and ( 2 ) increase of igm accretion onto galaxies . the former process is effective throughout most of the cosmic history , while the latter is effective only at @xmath2 when the igm is sufficiently enriched by metals owing to feedback . method : numerical  galaxies : evolution  galaxies : formation  galaxies : high redshift  galaxies : mass function  cosmology : theory