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: formulating strong - field transport in electron lattice has always been one of the most challenging theoretical goals in condensed matter physics  @xcite . this is more true with today s advanced nano - lithography techniques where we can now realize electron lattice driven far from equilibrium . even though this is an old problem , we are still in search of a firm theoretical paradigm to approach the problem in general . one of the central puzzles is dissipation . in equilibrium , the presence of an open environment in contact with a system introduces thermalization , and once the temperature is defined , we often use canonical or grand canonical ensemble , and do not consider the coupling to the environment as an explicit parameter . we naturally ask whether such simplifying ansatz may be possible in nonequilibrium , at least for steady - state description .    in nonequilibrium , we do not know such tremendously simplifying paradigms to take the role of the environment as an implicit parameter which can be hidden in thermalization process . on the contrary , the dissipation is considered as an integral part of the nonequilibrium process , and we need to include the dissipation mechanism explicitly for sound theoretical description . in some systems , however , the dissipation process can be simplified . in quantum dots ( qds ) under a finite dc bias , electron reservoirs coupled to the qd also act as energy source / drain and the energy relaxation is assumed to happen far away from qd and inside a battery . the electrical leads are then modeled as non - interacting reservoirs , and the electron transport is viewed as a transmission problem  @xcite . by taking the open limit , the excess energy can be taken infinitely far away from the qd , and the problem supports steady - state  @xcite . various quantum simulation methods have been proposed to study the transient behaviors of interacting models  @xcite . in the limit where a steady - state exists , the nonequilibrium state can also be described within the time - independent statistical mechanics framework  @xcite , from which hershfield  @xcite proposed the nonequilibrium density matrix @xmath6,\ ] ] with the reservoir energy @xmath7 for the source ( @xmath8 ) and drain ( @xmath9 ) with the continuum index @xmath10 . @xmath11 is the creation operator of the full scattering state as the solution of the _ whole _ system of quantum dot and the leads . @xmath12 is the chemical potential of the respective reservoirs . the scattering state formulation , although conceptually appealing , has initially been applied only to limited models  @xcite due to the difficulty of finding the scattering states . in recent years , several steady - state methods  @xcite have been developed and have been extended to general models . recently , the attention of the field has turned to lattice nonequilbrium problems . even at a very stage of the field , there have been important findings in the nonequilbrium processes , most notably that electrons under a dc electric field seem to build up internal energy quite quickly , reaching a quite different steady - state from equilibrium strongly correlated states  @xcite . one of the most popular technique of solving lattice many - body models has been the dynamical mean - field theory ( dmft )  @xcite . its success has been well documented in the description of the mott transition and strong - correlation physics in equilibrium . while there have been a fair amount of dmft works to electric - field driven lattice systems , the validity of the local approximation is still unconfirmed and subject to intense debate .    there have been numerous attempts to simulate nonequilibrium physics in lattice models , often through isolated hamiltonians  @xcite suited for quench dynamics of cold atom systems in optical lattice , periodically driven systems  @xcite , and some basic dissipation models  @xcite . however , in part due to the numerical difficulties in simulating long time - evolution , most of the efforts have focused on high - field phenomena such as the dielectric breakdown in mott insulators  @xcite . the main emerging picture of the calculations is that the external field drives the electronic lattice systems into hot temperature , generally regardless of the nature of many - body interaction . even though the picture is in agreement between many groups , the detailed understanding of the nature of the hot electron state and its eventual fate in more realistic setting is not available . to gain systematic understanding of such nonequilibrium state , the dissipation should be included in explicit models and their analytic behavior has to be studied with the damping as a controlled parameter . one of the goals of this paper is to introduce steady - state formulation via scattering state method for lattice with fermion baths and comprehensively analyze the model to show that the system possesses many properties which are expected of physical systems , for instance , consistent picture as the boltzmann transport theory . in the process , an argument will be made that the fermion bath model and the steady - state methods are a good minimal system to study nonequilibrium strong correlation physics . in the previous paper  @xcite by one of authors , the fermion bath model under a dc electric field has been shown to reproduce the key ingredients as predicted by the classical boltzmann transport theory , and to have a stead - state solution . the occupation number as a function of mechanical momentum has been shown to have the fermi sea shift by the drift velocity proportional to the scattering time and the electric field . furthermore , the dc current has been derived to be consistent with the boltzmann transport result applied to nanostructures  @xcite .    in this work , we further develop the solution to show that the scattering state formulation is applicable , and therefore a wide range of new techniques can be developed to solve the interacting lattice nonequilibrium phenomena . explicit calculations from temporal gauge and the coulomb gauge with scattering state formulation show that they are completely consistent with each other . the coulomb gauge enables the time - independent formalism , making physical interpretations more transparent . we calculate explicitly the local distribution function , from which we derive that the effective temperature scales as @xmath13 with a numerical constant @xmath1 , tight - binding parameter @xmath2 , the bloch oscillation frequency @xmath3 and the damping parameter @xmath4 . the effective temperature exhibits a singular limit of @xmath14 for @xmath5 . this proves that one should not take the damping as an implicit process , as treated in equilibrium theory . finally we derive , via the energy flux conservation with the joule heating , a general dc current relation as a functional of local green s functions as an extension of the meir - wingreen formula  @xcite to lattice models , and confirm the linear response theory . the main text of the paper is organized as follows . in section ii , the method introduced in ref . is further developed for green s functions in the temporal gauge . in section iii , we introduce the coulomb gauge and show that the green s functions in both gauges become identical in the long - time limit . in section iii we further discuss several important nonequilibrium quantities : the local distribution function and the effective temperature in a , time - evolution of wave - packet in b , dissipation and energy flux in c , and finally the derivation of the dc current in interacting models in d. appendices provide detailed analytic calculations . to demonstrate the equivalence of the time - dependent temporal gauge to the scattering - state formalism with time - independent coulomb gauge , we start with the one - dimensional non - interacting model considered earlier  @xcite . we study a quadratic model of a one - dimensional @xmath15-orbital tight - binding model connected to fermionic reservoirs ( see fig . [ fig1 ] ) under a uniform electric field @xmath16 . the effect of the electric field for time @xmath17 is absorbed in the temporal gauge as the peierls phase @xmath18 in the hopping integral  @xcite @xmath2 . here @xmath19 is the bloch oscillation frequency and @xmath20 is the step function . the time - dependent hamiltonian then reads @xmath21 with @xmath22 as the ( spinless ) electron operator on the tight - binding chain on site @xmath23 , @xmath24 with the reservoir fermion states connected to the site @xmath23 with the continuum index @xmath25 along each reservoir chain of length @xmath26 . the length @xmath26 is taken to infinity , and the time scale @xmath27 ( with fermi velocity of the chain @xmath28 ) for the wave to reach the end of the reservoir chain is considered larger than any other time scales in the problem . as discussed in ref . , the hamiltonian can be diagonalized in each @xmath10-sector as @xmath29 with @xmath30d^\dagger_kd_k+\sum_{\alpha}\epsilon_\alpha c^\dagger_{k\alpha}c_{k\alpha}\nonumber \\ & & -\frac{g}{\sqrt{l } } \sum_{\alpha}(c^\dagger_{k\alpha}d_k+{\rm h.c.}),\end{aligned}\ ] ] with the fermion operators fourier transformed to the wave - vector basis . here @xmath31 is the tight - binding dispersion at zero @xmath16-field . the reservoir states formed by @xmath32 acts as an open particle source with its chemical potential set at zero energy . the problem can be solved with @xmath33 as the unperturbed hamiltonian and @xmath34 as the time - dependent perturbation , @xmath35-\cos(k)\right\}d^\dagger_k d_k \equiv v(t)d^\dagger_k d_k.\label{eq : pert}\ ] ] this block - diagonal hamiltonian is nothing but a resonant level model coupled to a reservoir , with the level s energy oscillating in time  @xcite . with the perturbation of one - body terms of a finite degrees of freedom , one can write the dyson s equation for the retarded and lesser green s functions as  @xcite @xmath36{\bf g}_{k,0}^<[i+{\bf v}{\bf g}_k^a],\label{eq : glss}\end{aligned}\ ] ] where @xmath37 and @xmath38 are the lesser and retarded green s function matrices , respectively . @xmath39 and @xmath40 are for the non - interacting limit . the multiplication of green s function matrices denotes time integration .    following ref . , the retarded green s function is @xmath41 with the damping parameter @xmath42 for reservoirs of flat density of states of infinite bandwidth . the local retarded green s function @xmath43 becomes @xmath44 with the zero - th bessel function @xmath45 . here , the gauge - invariant local function becomes a function of only the relative time , @xmath46 . fourier transformation with respect to the relative time gives @xmath47 by using the bessel function relation  @xcite @xmath48 ^ 2e^{im\alpha}$ ] . the lesser green s function can be simplified in a straightforward calculation from eq . ( [ eq : glss ] ) following the similar procedures as in ref .   in the long - time limit ( @xmath49 ) as @xmath50 where @xmath51 with the self - energy @xmath52 from the damping taken as the perturbation . although the above equation has been derived with the time - dependent peierls term eq . ( [ eq : pert ] ) as the perturbation , the same result can be obtained when the damping is considered as perturbation in the limit that transient terms die out . the local lesser green s function can be computed as @xmath53 , which again renders the green s function only dependent on the relative time . after changing the dummy variables @xmath54 and @xmath55 , we have @xmath56 with @xmath57\nonumber \\ & & \times \sin\frac{\omega s}{2}\sin\frac{\omega s'}{2}.\nonumber\end{aligned}\ ] ] again by utilizing the bessel function relation  @xcite @xmath58 @xmath59 with @xmath60 , the fermi - dirac function at zero temperature . from the identity  @xcite , @xmath61 , @xmath62 although the two green s functions , eqs . ( [ eq : gtr ] ) and ( [ eq : gtlss ] ) , have been reduced to familiar forms resembling spectral representation , a clear relation between them is not available yet . in the following section , we discuss the scattering state formalism and find relations connecting the retarded and lesser green s functions . we have learned from ref .   and the above calculations that the fermion bath model has a well - defined time - independent limit for gauge - invariant quantities such as local green s function . this observation and the presence of infinite reservoir states prompt us to consider scattering state formalism  @xcite . as depicted in figs . [ fig2](a ) , for each site on the main chain there are infinite degrees of freedom coupled from each reservoirs . therefore , we can rewrite the quadratic hamiltonian in terms of the scattering states originating from the fermion reservoir chains  @xcite . to adopt the time - independent scattering theory , we use the coulomb gauge as shown in fig . [ fig2](a ) with the hamiltonian , @xmath63 where the static coulomb potential is applied as a potential slope to the chain and each reservoirs with their chemical potential are raised together with the corresponding tight - binding sites . since the hamiltonian is quadratic and the chain is coupled to open systems , we can decompose the hamiltonian in terms of the scattering state operators @xmath64 originating from the asymptotic state @xmath24 as given by the lippmann - schwinger equation  @xcite , @xmath65\nonumber \\ & = & c^\dagger_{\ell\alpha } -\frac{g}{\sqrt{l}}\frac{1}{\epsilon_\alpha-\ell\omega-{\cal   l}+i\eta}d^\dagger_{\ell},\end{aligned}\ ] ] with the liouvillian operator @xmath66 $ ] for any operator and @xmath67 . the scattering state embodies the openness of the fermion reservoirs . as implemented by the infinitesimal imaginary poles given by @xmath68 , once the electrons scatter into a reservoir they _ never _ come back to the tight - binding chain . with the scattering state basis , the hamiltonian is rewritten as @xmath69 with the fermi statistics applied separately within each @xmath23-sector @xmath70    since the hamiltonian is quadratic , the scattering state operator @xmath64 is a linear combination of the original basis of @xmath71 and @xmath72 , and @xmath73 with the expansion coefficient for a fermion annihilation operator @xmath74 given as @xmath75 it is important to realize that , for a quadratic hamiltonian , the anti - commutation of operators in @xmath76 is just a c - number . in fact , this c - number is nothing but the retarded green s function @xmath77 between @xmath78 and @xmath22 . from this argument , the retarded green s function in any time - independent quadratic hamiltonian is _ independent of statistics_. note that the scattering state @xmath64 originating from the @xmath23-th reservoir has admixture from any reservoir states @xmath72 . here , we have put an over - line to denote the green s function in the coulomb gauge . we then have @xmath79 this equation can be inverted to express @xmath80 in terms of the scattering state basis @xmath81 as @xmath82 similarly as above , @xmath83 since the scattering state operators also satisfy the anti - commutation relation  @xcite @xmath84 , and @xmath85 we will discuss below how @xmath86 is explicitly calculated . on - site retarded green s function at the central site @xmath87 can be expressed in terms of the scattering state basis as @xmath88^*}{\omega-\epsilon_\alpha+\ell\omega + i\eta}\nonumber \\ & & \times\langle\{\psi_{\ell\alpha},\psi^\dagger_{\ell'\alpha'}\}\rangle \nonumber \\ & = & \frac{g^2}{l}\sum_{\ell\alpha } \frac{|\overline{g}^r_{0\ell}(\epsilon_\alpha-\ell\omega)|^2 } { \omega-\epsilon_\alpha+\ell\omega+i\eta}. \label{eq : gwr}\end{aligned}\ ] ] here the retarded green s functions appear on both sides of the equation , and its self - consistency will be examined below . the lesser green s function can be calculated similarly as @xmath89^ * \nonumber \\ & & \times\delta(\omega-\epsilon_\alpha+m\omega ) \langle \psi^\dagger_{m\alpha}\psi_{m\alpha}\rangle \nonumber \\ & = & 2i\gamma\sum_m \overline{g}^r_{\ell m}(\omega)[\overline{g}^r_{\ell ' m}(\omega)]^ * f(\omega+m\omega ) , \label{eq : gwlss}\end{aligned}\ ] ] where each reservoir has its own fermi energy shifted by @xmath90 and @xmath91 . the above expression is quite appealing and physically transparent . with the dissipation provided by the particle reservoirs , all electron statistics are governed by the fermi statistics of the reservoirs and the effective tunneling between site @xmath23 and the reservoir attached at site @xmath92 is given by the retarded green s function @xmath93 . it is noted that we use the infinite - band approximation for each fermion reservoirs so that any reservoir can provide electrons to any other tight - binding lattice sites in principle , and all possible thermal factors mix throughout the lattice . now , we turn to calculation of retarded green s functions . with the time - independent hamiltonian , we only need to invert the matrix as @xmath94_{\ell\ell'}$ ] with @xmath95_{\ell\ell'}=(\omega+\ell\omega+i\gamma)\delta_{\ell\ell ' } + \gamma\delta_{|\ell-\ell'|,1},\ ] ] where the retarded self - energy @xmath96 is attached to each site @xmath23 of the tight - binding lattice with the potential slope . solution to the matrix inversion can be found  @xcite as @xmath97 which can be easily verified from @xmath98 . substituting this green s function into eq . ( [ eq : gwlss ] ) gives the identical result eq . ( [ eq : gtlss ] ) as derived from the time - dependent temporal gauge . coming back to the retarded green s function , we can easily confirm the identity eq . ( [ eq : gwr ] ) from a straightforward calculation after substituting eq . ( [ eq : grll ] ) into eq . ( [ eq : gwr ] ) and by using the contour integral and the completeness relation of bessel functions . from eq . ( [ eq : grll ] ) , it follows that @xmath99 from which eq . ( [ eq : gwlss ] ) satisfies @xmath100 in an interacting model , the self - energy is expressed in @xmath101 and inherits the same property , @xmath102      the discussion so far has demonstrated explicitly that the dissipative system with fermion baths can be described within the steady - state formalism using the scattering state basis . one of the central quantities to calculate is the effective local distribution function , @xmath103 where eq . ( [ eq : gwr ] ) has been used for @xmath104(a ) shows the numerical evaluation of the local distribution function @xmath105 . for @xmath106 , @xmath105 is a superposition of small steps coming from the thermal factors @xmath107 in eq . ( [ eq : floc ] ) with the envelope following a smooth profile similar to the fermi - dirac function . even though there is no reason to expect that the nonequilibrium distribution mimics the fermi - dirac function , we can nevertheless fit the result to the function with an effective temperature @xmath108 as shown , despite some deviation ( see appendix a for more details ) . as @xmath108 grows towards the finite tight - binding bandwidth @xmath109 , the fit becomes inaccurate . the envelope of the local distribution function plays a similar role of the fermi - dirac function which dictates the abundance of electron - hole pairs available for interaction . in the presence of additional many - body interactions such as the electron - phonon coupling to local optical phonons , the available electron - hole pairs for inelastic dissipation are given by the @xmath105 profile and the effective temperature @xmath108 in the fermi - dirac function will play the role of hot electron temperature effectively . it is remarkable that @xmath108 seems to approach infinity as the damping parameter @xmath4 becomes smaller . although it may look counter - intuitive at first , this is only the manifestation of the _ short - circuit _ behavior where a finite voltage applied across a low resistance conductor induces an extremely hot temperature . this is also consistent with the numerical calculations with the general conclusion that the electron temperature reaches an infinity in closed interacting models . numerical fit indicates that @xmath108 is an increasing function of @xmath110 for a wide range of @xmath111 , although the functional form eventually deviates from the form as shown in fig . [ fig3](b ) . for small @xmath110 , the effective temperature behaves as @xmath112 with a dimensionless numerical constant @xmath1 . this equation is one of the key results of this paper . in appendix a , we derive the above linear dependence of @xmath113 and approximately estimate that the constant @xmath114 by analyzing the @xmath115 step in @xmath116 . the effective temperature has been previously observed in the momentum distribution function and has been speculated  @xcite to behave as @xmath117 based on the dc conductivity analogy . more careful and quantitative analysis now shows that the correct dependence is the above relation , eq . ( [ eq : teff ] ) . the @xmath108 relation is even further corroborated with @xmath1 derived from the kubo formula ( see appendix b and discussions in section iii.d . ) as the blue line in fig . [ fig3](b ) . the kubo formula result @xmath118 should be exact for the limit @xmath119 . the two analytical estimates bracket the numerical @xmath108 [ see fig . [ fig3](b ) ] , which shows that the expressions eqs .  ( [ eq : teff ] ) and ( [ eq : teff2 ] ) are a reliable approximation for @xmath4 and @xmath3 up to @xmath2 .    the divergent effective temperature should be taken with a caution to interpret in finite bandwidth systems . unlike with the quadratic dispersion relation for continuum models  @xcite , the kinetic energy in the single - band tight - binding model is always bounded and thus can not give off extremely _ hot_-electrons to the environment .    as @xmath3 and @xmath4 become comparable to the bandwidth [ see fig . [ fig3](c ) ] , the signature of bloch oscillation steps become more obvious and the definition of the effective temperature as determined by the shape of the overall @xmath105 is not very robust . while the infinite temperature in a finite bandwidth system may be questionable , the trend observed in fig . [ fig3](a - b ) continues . in the small damping limit at @xmath120 , even a population inversion happened in the local distribution function so far , we have seen that the steady - state formalism provides a convenient theoretical framework for nonequilibrium lattice model of fermion baths . to better understand the nonequilibrium steady - state , we now look at time - dependent quantity indicative of wave - packet drift . a steady - state by definition is a time - independent reference state where a direct observation of time - evolution of a moving particle can not be made . to confirm that an electric charge moves down the potential slope as a function of time , we create a hole out of occupied states  @xcite at a central position and observe its movement to a different position @xmath121 as a function of time in @xmath122 the lesser green s function can be easily decomposed in terms of scattering states and fig . [ fig4 ] shows the wave - packet traveling in the direction of the applied field . amplitude of the observable decays as @xmath123 due to the dephasing of electrons from the fermion baths .    as depicted in fig . [ fig2](b ) , electrons travel down the potential slope in the coulomb gauge by creating a trail of electron - hole ( e - h ) pairs in the reservoirs . as long as the bandwidth of the reservoirs is greater than the potential drop @xmath3 between neighboring sites ( we assumed that the bandwidth is infinite in explicit calculations ) , each reservoirs can accommodate an e - h pair with its energy matching @xmath3 by particle exchange via tunneling . for narrower reservoir bandwidths , multiple e - h pairs should be created to establish a dc current . since the open reservoirs are of infinite length , the created e - h pairs travel indefinitely inside the reservoirs and therefore the fermion baths can produce effects similar to the inelastic processes in bosonic baths . we discuss this further in the subsequent sections . we turn to discussions of energy dissipation . the hamiltonian , eq . ( [ eq : hcoul ] ) , can be divided into three parts as @xmath124 with each term representing each line in eq . ( [ eq : hcoul ] ) , respectively . in the steady - state limit , the energies stored in @xmath125 and @xmath126 are stationary @xmath127 , as will be demonstrated below . unlike the case with @xmath125 and @xmath126 which are of finite spatial extent , the energy flux in @xmath128 can be non - zero . in the scattering theory  @xcite , the scattering states are formulated in the limit that the spatial extent of the scattered wave ( @xmath129 ) into the reservoirs is much shorter than the length of the reservoir chain ( @xmath130 ) . therefore the scattering state represents a solution that the scattered wave _ constantly propagating inside the reservoirs _ without being backscattered from the edge of the reservoirs , and quantities involving the extended states @xmath131 or @xmath132 do not have to be stationary in general  @xcite .    to be more concrete , we discuss explicit calculations . with the fermion baths , a dc electric - field establishes a dc current @xmath133 , as calculated in ref .  . to investigate the effect of the joule heating , we consider @xmath125 @xmath134 \right\rangle \\ & = & \omega\langle \hat{j}\rangle   + i\gamma g\sum_\ell\langle ( \bar{c}^\dagger_{\ell+1}+\bar{c}^\dagger_{\ell-1})d_\ell -{\rm h.c . } \rangle\nonumber\end{aligned}\ ] ] with the current operator @xmath135 within the main chain and @xmath136 . we used the steady - state condition for the occupation @xmath137 . the first term represents the joule heating and the second term the energy flux of electrons from the kinetic energy of the main chain into the coupling @xmath126 . denoting the energy flux per site as @xmath138 , we show that @xmath139 . we present the detailed calculations in appendix c and show the equality of the above equation based on the scattering - state formalism . it can be shown further that @xmath140 . as shown in appendix b , the energy influx to each of the reservoirs is nothing but the joule heating @xmath141d\omega\nonumber \\ & = & \omega\langle j\rangle , \label{joule}\end{aligned}\ ] ] with the local spectral function defined as @xmath142 the lowercase hamiltonian @xmath143 denotes the corresponding hamiltonian per tight - binding site . it might sound paradoxical that the energy in the electronic system is non stationary , @xmath144 ( here the last equality is from the steady - state current taken from ref .  . ) this is due to the fact that , although @xmath145 governs the electron dynamics , there is another part of hamiltonian which should be included for a closed system  the battery connected across the tight - binding chain . since the battery loses its stored charge @xmath146 with the rate of @xmath147 , the electrostatic energy decrease per unit cell of the tight - binding chain becomes @xmath148 , and the total energy @xmath149 is stationary in the steady - state . the discussion here again confirms the picture depicted in fig . [ fig2](b ) where the fermion baths act as energy reservoirs while the net electron number flux into the reservoirs is zero . despite their simplicity , the fermion baths through their particle - hole excitations play the role of bosonic baths , apart from the boson s explicit dispersion relation ( with the exception of the luttinger liquid bath ) and the physics that might occur from the nonlinear effect of the bosonic statistics .      from the energy dissipation relations above , we obtain the useful formula for the steady - state current , @xmath150d\omega , \label{eq : current}\ ] ] where only on - site green s functions are needed as in meir - wingreen formula in quantum dot transport  @xcite . to recover the ohm s law for small field that @xmath151 one should have that the integral goes as @xmath152 as the leading order . this is justified since applying a field of opposite direction @xmath153 should not change the local properties and the integral should be of order @xmath152 . this argument can be used to analyze the linear response limit , as described below . the above relation eq . ( [ eq : current ] ) , verified explicitly for the non - interacting model in appendix c , can be extended to interacting models . the key identities are steady - state conditions @xmath154 which we expect to hold generally for interacting systems as long as the interaction potential does not hold infinite amount of energy per site , as in hubbard model . here we used the spin index @xmath155 . with on - site interaction , @xmath156 ensures zero particle flux into the baths . the energy flux equation and the dysone equation hold as eq . ( [ app : bath ] ) and ( [ app : dyson ] ) , respectively . this immediately shows that the current ,  ( [ eq : current ] ) , holds for a wide range of interacting fermion bath models . steady - state current derived by meir - wingreen  @xcite has been widely used in quantum dot calculations . the formula  ( [ eq : current ] ) can be seen as its extension for lattice models with fermion baths . the equation is a functional of only local green s functions . however , it should be made clear that while the simplified meir - wingreen formula for a single quantum model only requires @xmath157 for the quantum dot , both of @xmath158 and @xmath159 are necessary in the lattice models . using the key equation ( [ eq : current ] ) , and eq . ( [ eq : teff ] ) in the non - interacting limit , a linear response limit can be analyzed . in the limit of @xmath119 , the effective temperature is expected to be small and , therefore , we can use the sommerfeld expansion  @xcite to derive the linear electrical current @xmath160 and we obtain the linear dc conductivity @xmath161 here @xmath162 is the equilibrium spectral function evaluated at the fermi energy . with @xmath163 , @xmath164 and the scattering time @xmath165 , we recover the drude conductivity  @xcite . comparing this with the linear response theory using the kubo formula , we obtain @xmath166 as detailed in appendix b. noting that @xmath167 , we have from eq . ( [ eq : cond ] ) that @xmath168 . this result is shown as the blue line in fig . [ fig3](b ) in comparison to the numerically obtained @xmath108 .    in the interacting limit , the effective temperature expression eq .  ( [ eq : teff ] ) should be modified . we expect that the same form @xmath169 holds with renormalized parameters @xmath170 , @xmath171 and @xmath172 . then the linear response equation becomes @xmath173 . in this work , we have reformulated the electron transport in tight - binding lattice driven by a dc electric field using both time - dependent and time - independent gauges . the time - independent coulomb gauge with fermion baths leads to the scattering state description for steady - state , which makes the calculation and interpretation more intuitive . nonequilibrium quantum statistics of quantum dot model , as proposed by hershfield  @xcite , can be extended to nonequilibrium lattice as summarized in the scattering state expressions , @xmath174\\ \hat{h } & = & \sum_{\ell\alpha}(\epsilon_\alpha-\ell\omega)\psi^\dagger_{\ell\alpha } \psi_{\ell\alpha},\end{aligned}\ ] ] with the inverse temperature @xmath175 of the baths . the reservoir scattering states ( represented by @xmath176 ) are shifted by the applied electrostatic potential @xmath90 , and the chemical potential is simultaneously shifted with the electrostatic potential . therefore , the energy spectra governing the dynamics and statistics are different in the above expression . the formalism provides a natural framework for approximations such as the dynamical mean - field theory ( dmft ) . it has been shown that the fermion bath model , although quite rudimentary , produces dissipation mechanism consistent with the boltzmann transport theory . in particular , the steady - state effective temperature induced by the external field depends quite strongly on the electric field and the damping . the effective temperature becomes divergent as @xmath177 for small damping @xmath4 versus the bloch frequency @xmath19 ( @xmath74 is the lattice constant , and @xmath16 the electric field ) . although this might look surprising at first , this phenomenon is simply the manifestation of the short - circuit effect . it also verifies various numerical calculations with the infinite electron temperature resulting in isolated lattice models . these findings have fundamental implications in nonequilibrium quantum statistics in that dissipation processes can not be implicitly included as thermalization as in the boltzmann factor of equilibrium gibbsian statistics . through the energy dissipation and the joule heating in the fermion reservoirs , a general dc current relation eq . ( [ eq : current ] ) has been derived for interacting models , as an extension of the meir - wingreen formula to nonequilibrium lattice systems . the linear response limit has been confirmed within this formalism .    despite the lack of momentum scattering and explicit inelastic processes , the generic features of the fermion bath model which are consistent with semi - classical theory are quite significant . furthermore , for its simplicity the fermion bath model can be used as an ideal building block for studying strong correlation effects in lattice driven out of equilibrium . particularly , with the time - independent coulomb gauge dmft can be readily formulated using the scattering state method  @xcite it is well - known in equilibrium strong correlation physics that electrons undergo collective state when a strong interaction is present , with some emergent energy scale @xmath178 . one may speculate that an electric field of order @xmath179 would significantly alter the strongly correlated state . however , our study suggests that the dissipation strongly interplays with the nonequilibrium condition and non - trivial physics may arise even at @xmath180 . further systematic studies are necessary to understand the interplay of nonequilibrium and strong correlation effects . we thank helpful discussions with kwon park , woo - ram lee , jainendra jain , anthony leggett , natan andrei and gabi kotliar . this work has been financially supported by the national science foundation through grant no . dmr- 0907150 . here we analytically justify the relation @xmath181 by considering the low frequency steps as shown in fig . [ fig5 ] . the first step @xmath182 at @xmath115 can be expressed as @xmath183 \right\}^{-1}.\ ] ] here we look at the limit of small @xmath3 and approximate @xmath184 by the equilibrium green function @xmath185^{1/2}.\ ] ] therefore , the analytic expression for the slope of the fit becomes @xmath186 by equating this to the slope of the effective fermi - dirac function @xmath187 , we obtain @xmath188 note that while the actual numerical fit overestimates @xmath108 from the analytic expression due to the high - frequency contribution , the overall functional dependence is quite reliable for @xmath189 . from the kubo formula  @xcite , the linear conductivity can be exactly calculated in the small @xmath3 limit . for convenience , we calculate the current - current correlation function in the imaginary - time formalism and then analytically continue to the real - frequency in the optical conductivity  @xcite . for the uniform ( @xmath190 ) response function in the matsubara frequency @xmath191 , the conductivity is expressed as @xmath192 here @xmath193 is the group velocity and the matsubara green s function for the electron is given as @xmath194 with @xmath195 . performing the matsubara summation and then the analytic continuation @xmath196 for finite @xmath197 , we have @xmath198 taking its real part and the static limit @xmath199 at zero temperature , we obtain the dc linear conductivity @xmath200 with the restored constants @xmath201 and @xmath74 . the current expectation value measured at the site @xmath87 is expressed as @xmath202 . using the scattering state basis , we have @xmath203^ * \overline{g}^r_{0n}(\epsilon_\alpha - n\omega)f(\epsilon_\alpha ) \nonumber \\ & = & \frac{\gamma}{\pi}\sum_{nkk'}\int_{-\infty}^0 d\omega \frac{j_{1-k - n}j_{-k}j_{-k'-n}j_{-k'}}{(\omega+k\omega - i\gamma ) ( \omega+k'\omega+i\gamma ) } \nonumber \\ & = & -\frac{\gamma}{\pi}\sum_{k}\int_{-\infty}^0   \frac{j_{k}j_{k-1}d\omega } { ( \omega+k\omega - i\gamma ) [ \omega+(k-1)\omega+i\gamma]},\nonumber\end{aligned}\ ] ] by using eq . ( [ eq : grll ] ) and the bessel function identities . we have suppressed the argument in the bessel functions . with integration of elementary functions we obtain @xmath204,\end{aligned}\ ] ] with @xmath205 . this immediately confirms that the current evaluated from the scattering state basis matches the result in ref .  .    for the operator @xmath206 we need to calculate @xmath207 . the @xmath208-operators are expressed with the scattering state basis as @xmath209^ * \psi^\dagger_{\ell'\alpha ' } } { \epsilon_{\alpha'}-\epsilon_\alpha-(\ell'-\ell)\omega - i\eta}\right].\ ] ] a lengthy but straightforward calculation gives @xmath210 which confirms the identity @xmath211 . for the energy flux into the fermion baths , we examine @xmath212\rangle\nonumber \\ & = & \frac{ig}{\sqrt{l}}\sum_{\ell\alpha}(\epsilon_\alpha-\ell\omega)\langle c^\dagger_{\ell\alpha}d_\ell - d^\dagger_\ell c_{\ell\alpha}\rangle . \label{app : bath}\end{aligned}\ ] ] first we show that @xmath213 . from the steady - state condition of @xmath214 , we have @xmath215\rangle = -\gamma\langle d^\dagger_{\ell + 1}d_\ell -d^\dagger_\ell d_{\ell + 1}+d^\dagger_{\ell -1}d_\ell   -d^\dagger_\ell d_{\ell -1 } \rangle -\frac{g}{\sqrt{l}}\sum_\alpha \langle c^\dagger_{\ell\alpha}d_\ell - d^\dagger_\ell c_{\ell\alpha}\rangle = 0 $ ] . the first term is the total flux into the @xmath23-th site due to the current along the tb chain , and in the steady - state it is zero . therefore we have zero particle - flux into the reservoir , @xmath213 . the remaining summation @xmath216 is the energy flux measured with respect to the @xmath23-th reservoir chemical potential level and it should be independent of @xmath23 . setting @xmath87 , we can rewrite the expression as follows . consider @xmath217 and @xmath218 . for the energy flux per reservoir , we can write @xmath219\frac{d\omega}{2\pi}.\nonumber\end{aligned}\ ] ] from the dyson s equation  @xcite , @xmath220\ } \nonumber \\ & & = \frac{4\pi^2 g}{\sqrt{l}}\delta(\omega-\epsilon_\alpha ) a_{\rm loc}(\omega)[f_{\rm loc}(\omega)-f(\omega ) ] . \label{app : dyson}\end{aligned}\ ] ] then @xmath221d\omega.\ ] ]    taking the time - derivative of the above green s functions at @xmath222 , we have from eq . ( [ app : dyson ] ) , @xmath223\right|_{t=0 }   \nonumber \\ & & = 2\gamma\int\omega a_{\rm loc}(\omega)[f_{\rm loc}(\omega)-f(\omega)]d\omega=\langle \dot{\hat{h}}_{\rm bath}\rangle.\nonumber\end{aligned}\ ] ] this can be equated to @xmath224+[\hat{h}_{\rm coul},d^\dagger_0]c_{0\alpha}\rangle \nonumber \\ & = & -ig\gamma\langle \bar{c}^\dagger_0(d_{1}+d_{-1})-(d^\dagger_1+d^\dagger_{-1})\bar{c}_0\rangle \nonumber \\ & = & -ig\gamma\langle ( \bar{c}^\dagger_1+\bar{c}^\dagger_{-1})d_0-d^\dagger_0(\bar{c}_1+\bar{c}_{-1})\rangle,\end{aligned}\ ] ] where the translational invariance of the steady - state has been used . using eq . ( [ app : cd0 ] ) , eq . ( [ joule ] ) is confirmed .
Abstract:

Answer: electron transport in electric - field - driven tight - binding lattice coupled to fermion baths is comprehensively studied . we reformulate the problem by using the scattering state method within the coulomb gauge . calculations show that the formulation justifies direct access to the steady - state bypassing the time - transient calculations , which then makes the steady - state methods developed for quantum dot theories applicable to lattice models . we show that the effective temperature of the hot - electron induced by a dc electric field behaves as @xmath0 with a numerical constant @xmath1 , tight - binding parameter @xmath2 , the bloch oscillation frequency @xmath3 and the damping parameter @xmath4 . in the small damping limit @xmath5 , the steady - state has a singular property with the electron becoming extremely hot in an analogy to the short - circuit effect . this leads to the conclusion that the dissipation mechanism can not be considered as an implicit process , as treated in equilibrium theories . finally , using the energy flux relation , we derive a steady - state current for interacting models where only on - site green s functions are necessary .


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 stellar initial mass function ( imf ) has been traditionally interpreted to be an invariant probability density distribution function ( e.g. kroupa 2001 , 2002 ; bastian et al . recent observational evidence has however ruled out this simple interpretation in that low - density regions have been found to be lacking the massive stars they ought to have produced given the large number of stars formed ( kirk & myers 2011 ; hsu et al . the mass function of young star clusters has also been shown to not be a probability density distribution function ( pflamm - altenburg et al . 2013 ; kroupa 2014 ) . furthermore , several independent observational evidences have been put forward in the last few years , showing that the galaxy - wide imf varies with star formation rate ( sfr ) ( e.g. hoversten & glazebrook 2008 ; gunawardhana et al . 2011 ; cappellari et al . 2012 ; conroy & van dokkum 2012 ) and probably varies also within a single galaxy ( dutton et al . 2013 ) .    a comprehensive theory , attempting to explain the evolution of the imf within galaxies is the so - called integrated galactic imf ( igimf , i.e. the galaxy - wide imf ) theory ( kroupa & weidner 2003 ; weidner & kroupa 2005 ; see the recent review by kroupa et al . 2013 ) . in two previous papers ( recchi et al . 2009 ; recchi et al . 2014 ) we have investigated the effect of the igimf on the abundance ratios ( global and of individual stars ) in galaxies . one paper ( calura et al . 2010 ) has been devoted to the study of the chemical evolution of the solar neighborhood by adopting the ( sfr - dependent ) igimf formulation . eventually , numerical hydrodynamical simulations of dwarf galaxies with a variable igimf have been pioneered , too ( ploeckinger et al . 2014 ) . it remains to be shown that the igimf theory is consistent with the observed mass - metallicity ( mz ) relation in galaxies ( see e.g. lee et al . 2006 ; kirby et al . 2013 ) and with the more general fundamental metallicity relation ( fmr , see mannucci et al . 2010 ; lara - lpez et al . 2010 ) . we build on a previous attempt ( kppen et al . 2007 ) and solve the equations of the so - called `` simple model '' but within the framework of the igimf theory , i.e. we assume that the galaxy - wide imf depends on the sfr and on the metallicity of the parent galaxy . analytical and semi - analytical calculations of the chemical evolution of galaxies based on simple models are nowadays quite popular , as they enlighten in a simple way complex correlations among galaxies ( spitoni et al . 2010 ; dayal et al . 2013 ; lilly et al . 2013 ; birrer et al . 2014 ; pipino et al . it is worth remarking that all these authors consider invariant imfs whereas , as explained above , there is significant observational evidence for imf variations among different galaxies . therefore , it is important and timely to understand how to modify the formalism of the simple model to take into account imf variations . this is indeed the main aim of the present paper . it is also worth stressing that , although this paper is based on the igimf theory , other theories of imf variation have been put forward in the literature ( e.g. larson 1998 ; padoan & nordlund 2002 ; hopkins 2013 . see also martinelli & matteucci 2000 ) . the formalism developed in this paper can be applied equally well to these imf theories . the paper is organized as follows : in sect . [ sec : igimf ] we recall the main assumptions of the igimf theory and describe the way in which this theory has been implemented in the present work . in sect . [ sec : simple ] we describe the simple models of chemical evolution  a simplified tool to describe the evolution of the metallicity of a galaxy . in sect . [ sec : cbrev ] we re - interpret the so - called closed - box model within the framework of the igimf theory , i.e. we relax the hypothesis that the imf is universal . similarly , in sect . [ sec : lbrev ] we re - interpret simple models of chemical evolution with outflows and infall and consider again variations of the imf between galaxies and within galaxies . finally , in sec . [ sec : disc ] the results of this work are discussed and some conclusions are drawn . in a companion paper ( recchi et al . 2014 ; hereafter paper ii ) , we will apply this methodology to the study of the mass - metallicity relation in tidal dwarf galaxies ( tdgs ) . the igimf theory is based on the assumption that stars in a galaxy are formed in space - time correlated entities called embedded star clusters ( lada & lada 2003 ) , thereby accounting naturally for the above mentioned results by kirk & myers ( 2011 ) , hsu et al . ( 2012 ) and pflamm - altenburg et al . the notion underlying the igimf theory lies in the conjecture that the freshly born galaxy - wide stellar population is the sum of all local star - formation activity and that this activity occurs in the density maxima of molecular cloud cores by being made up of embedded groups or clusters of stars on a spacial scale of less than a  pc ( marks & kroupa 2012 ) . the mass distribution of these is equivalent to the mass function of embedded clusters which do not need to be gravitationally bound upon expulsion of their residual gas ( lada et al . 1984 ; kroupa et al . 2001 ; banerjee & kroupa 2013 ) . an `` isolated '' or low - density mode of star formation is naturally accounted for in this picture by low - mass and/or the radially extended part of the embedded clusters .    in defining the igimf theory one posits the following six axioms derived from observations ( see kroupa et al . 2013 ; weidner et al . 2013 ) and @xmath0 are in units of @xmath1 . ]    1 . the stellar imf , @xmath2 , in an embedded star cluster is canonical ( a two - part power law @xmath3 with @xmath4 for 0.08 m@xmath5 @xmath6 0.5 m@xmath5 , and @xmath7 for stellar masses larger than 0.5 m@xmath5 ) for cloud core densities @xmath8 m@xmath5 pc@xmath9 . a star formation efficiency @xmath10 in molecular cloud cores of 0.33 has been assumed ( alves et al . 2007 ) . here @xmath0 is the stellar mass in the embedded cluster and @xmath11 is the residual gas mass within it . 2 .   the embedded clusters populate an embedded cluster mass function ( ecmf ) , which is assumed to be a power law of the form @xmath12 ( @xmath13 would be the salpeter index ) . the half - mass radii of embedded clusters follow @xmath14 ( pc ) = 0.1 @xmath15 ( marks & kroupa 2012 ) yielding for the total density ( gas plus stars ) log@xmath16(@xmath17 ) = 0.61 log@xmath16(@xmath0 ) + 2.85 , in units of m@xmath5 pc@xmath9 . 4 .   the most massive star in a cluster , @xmath18 , is a function of the stellar mass of the embedded cluster , @xmath0 ( weidner & kroupa 2004 , 2006 ; weidner et al . 2010 ) , @xmath19 @xmath20 ( e.g. eq . 10 in pflamm - altenburg 2007 ) . 5 . there exists a relation between the sfr of a galaxy and the most massive young ( age@xmath2110 myr ) star cluster , log@xmath16(@xmath0 ) = 0.764 log@xmath16(sfr)+4.93 ( weidner et al . 2004 ) , where the sfr is in units of m@xmath5 yr@xmath22 . the mass of the least - massive embedded `` cluster '' is assumed to be 5 m@xmath5 the dependence of the imf slope , @xmath23 , of stars above 1 m@xmath5 on the initial density and metallicity of the embedded cluster is given by ( marks et al . 2012 ) @xmath24 where @xmath25}\,+0.99 \log_{10}\left(\frac{\rho_{\rm cl}}{10 ^ 6   { \rm m}_{\odot } { \rm pc}^{-3}}\right).\ ] ] this dependency of the stellar imf on density and metallicity of a massive star - burst cluster has been derived by combining independent evidence from globular clusters ( marks et al . 2012 ) and ultra - compact dwarf galaxies ( dabringhausen et al . 2009 , 2012 ) .    for a value of @xmath26 the igimf built in this way nicely reproduces many observed properties of galaxies ( see kroupa et al . it also qualitatively reproduces the trend of decreasing @xmath23 with increasing sfr observed by gunawardhana et al . ( 2011 ) . in order to improve this fit , a seventh axiom has been introduced ( weidner et al . 2013 ) : the exponent @xmath27 of the ecmf is not a constant but it varies with the sfr according to : @xmath28 axiom  7 implies the mass function of embedded clusters to be top - heavy in star bursts . weidner et al . ( 2013 ) explored also the possibility that the mass of the smallest embedded cluster depends on the sfr , thus partially relaxing the fifth axiom above . they showed however that this assumption produces results almost indistinguishable from the results obtained with axiom 7 , thus we will not further consider this hypothesis . a class of galactic chemical evolution models , commonly dubbed as simple models , admit analytical solutions . they are thus quite useful and easy - to - apply diagnostics for the evolution of galaxies ( see spitoni et al . 2010 ; lilly et al . 2013 ; pipino et al . 2014 ; see also the introduction ) , thus it is useful to understand advantages and limitations of this formalism . the basic assumptions behind the simple models are ( see tinsley 1980 ; maeder 1992 ) :    1 . stars more massive than a fixed mass ( usually 1 m@xmath5 ) die instantaneously . stars less massive than this threshold live forever . this assumption ( instantaneous recycling approximation ) allows us to neglect the lifetime of single stars in what follows . 2 .   the gas is well mixed at any time . 3 . the flow rates ( infall and outflow ) are proportional to the sfr @xmath29 . the imf is universal : it is independent of time and it is the same in all galaxies . 5 .   the yields and the remnant masses are functions of the stellar mass alone . this assumption ( together with the universality of the imf ) implies that the quantities : @xmath30 are constant . here , @xmath31 is the adopted ( fixed , for the moment ) imf , @xmath32 is the remnant mass of the star with initial mass @xmath33 , and @xmath34 is the mass - fraction of newly synthesised metals by the star of initial mass @xmath33 . according to this terminology , the total amount of heavy elements released by a star of initial mass @xmath33 , initial metallicity @xmath35 and remnant mass @xmath32 is @xmath36 , i.e. it includes newly synthesised heavy elements but also metals present in the stars at birth and not processed . the quantity @xmath37 in eq . [ eq : ryz ] represents thus the fraction of a stellar population not locked into long - living ( dark ) remnants and @xmath38 represents the ratio between the mass of heavy elements ejected by a stellar generation and the mass locked up in remnants . it is well known that these assumptions lead to a set of differential equations for the time evolution of the total baryonic mass @xmath39 , total gas mass @xmath40 and total mass in metals @xmath41 within a galaxy ( see tinsley 1980 ; maeder 1992 ; pagel 1997 ; matteucci 2001 ) : @xmath42\end{cases } \label{eq : system}\ ] ] here , @xmath43 is the metallicity of the infalling material , whereas @xmath44 and @xmath45 are proportionality constants relating the sfr to the infall and outflow rate , respectively . in particular , it is assumed that : @xmath46 ( see matteucci & chiosi 1983 ) where @xmath47 and @xmath48 are the infall or accretion and the outflow rate , respectively . given the initial condition @xmath49 , @xmath50 , a solution as a function of the mass ratio @xmath51 can be found ( see also recchi et al . 2008 ) : @xmath52^ {        \frac{\lambda}{\lambda-\lambda-1}}\right\rbrace . \label{eq : simplechevsol}\ ] ] this solution expresses the gas - phase metallicity of a galaxy , as a function of the mass fraction @xmath53 . the time dependence of this relation is hidden in the time dependence of @xmath53 , which in turn depends on how quickly the gas is consumed to form stars ( i.e. it depends on the sfr ) and on the flows in and out of the galaxy . hence , to have a complete picture of the evolution of the galaxy , we need also some assumption concerning the sfr . in what follows we will assume that the sfr depends on the available gas reservoir through the simple law @xmath54 i.e. we assume here that the sfr is linearly proportional to the gas mass and that this proportionality constant is an invariant quantity . this equation also implies that the gas consumption timescale of a galaxy is always about 3 gyr . notice that many studies of the chemical evolution of galaxies suggest that the star formation efficiency increases with the mass of the galaxy ( e.g. matteucci 1994 ; pipino & matteucci 2004 ) . these studies are however based on invariant imfs and indeed matteucci ( 1994 ) showed that a variable imf or a variable star formation efficiency might produce similar results . in the framework of the igimf , pflamm - altenburg & kroupa ( 2009 ) showed that the gas consumption time - scale is similar in all galaxies , irrespective of their masses , and it is indeed of the order of 3 gyr . we consider also the average metallicity @xmath55 of the stars in a galaxy . a simple way of calculating this quantity is : @xmath56 this expression represents the mass - weighted average of the metallicities of all the stellar populations ever born in the galaxy ( see pagel 1997 ) . a more accurate definition of @xmath55 would involve a luminosity - weighted average , but it has been shown that these two averages do not differ much ( e.g. recchi et al . 2009 ) .    a solution can be found also in the presence of metal - enhanced galactic outflows ( see sect . [ sec : lbrev ] ) . a solution can also be obtained assuming generic infall and outflow laws ( see again recchi et al . 2008 ) . from the given solution [ eq : simplechevsol ] , other particular cases can be found . for instance , setting @xmath57 one obtains the well - known solution of the closed - box model @xmath58 a @xmath44 is in the denominator , hence @xmath59 diverges as @xmath60 . the previous sentence might thus look suspicious . the fact is that also the expression within the curly brackets tends to zero for @xmath61 . take @xmath62 ; the solution reduces to @xmath63        \simeq \frac{y_{\rm z}}{\lambda}\left [ 1 -          \left(\frac{1}{\mu}\right)^{-\lambda}\right],\label{eq : foot}\ ] ] where we have taken the limit for small @xmath44 in the second equality . write @xmath64 as @xmath65 and expand it in a taylor series , e.g. @xmath66 . by substituting this expansion in the expression for @xmath35 in eq . [ eq : foot ] , the closed box solution @xmath67 is recovered . ] it is generally accepted that this solution gives a first - order accurate description of the overall evolution of chemical elements ( such as the @xmath68-elements ) mainly released by massive stars on short timescales . the relevant equations for the closed - box model with a linear sfr such as in eq . [ eq : sflaw ] are : @xmath69 \cr      \psi ( t ) = s m_g ( t ) \end{cases } \label{eq : cbsystem}\ ] ] it is trivial to solve for @xmath70 , @xmath71 . assuming that there is only gas at the beginning and that @xmath72 , one obtains : @xmath73}. \label{eq : solcb}\ ] ] it is important to remark on this equation that , according to the assumptions of the simple model , all constants that enter in the calculation of @xmath70 are indeed constants , i.e. they do not depend on the mass of the model galaxy taken into consideration . the consequence is that the time evolution of the metallicity is always the same , independent of the mass of the galaxy . consequently , these models can not reproduce the observed mz relation . a viable solution is to assume that small galaxies can experience galactic winds , due to their shallow potential wells . observations show indeed that mass - loading factors in galactic winds are larger in small galaxies ( martin 2005 ) and theoretical investigations confirm this finding ( hopkins et al . 2012 ; puchwein & springel 2013 ) . if the parameter @xmath45 introduced in sect . [ sec : simple ] is large in small galaxies and small or zero in large ones , the mz relation can be nicely reproduced ( tremonti et al . 2004 ; spitoni et al . 2010 ) . although we are convinced that galactic winds can play an important role in the chemical evolution of dwarf galaxies ( recchi & hensler 2013 ; see also sect . [ sec : lbrev ] ) , we want to explore here the possibility that the imf is not constant but that it depends on time though its dependence on the sfr and on z , as per the igimf theory described in sect . [ sec : igimf ] . as explained in the introduction , this exploration is purely didactic and serves to illustrate in the simplest possible setting how a variable imf affects the formalism of the simple model described in sect . [ sec : simple ] .    clearly , the assumption of a variable imf introduces a differentiation in the metallicity evolution of model galaxies with different masses , as the sfr of different model galaxies are different and that implies variations in the igimf , which in turn imply variations in the chemical evolution .    to be more quantitative , the quantities @xmath37 and @xmath74 introduced in eq . [ eq : ryz ] become :    @xmath75dm,\notag\\    & y_{\rm z}(t)=\frac{1}{1-r}\int_1^\infty m \ , p_{{\rm z } m } \ , \xi_{\rm igimf }    \left[m,\psi(t),z(t)\right]dm ,    \label{eq : ryz2}\end{aligned}\ ] ]    i.e. they are not constant any more , but are functions of time , through the time dependence of @xmath29 and @xmath35 . the imf @xmath76 entering here is the igimf , i.e. the imf obtained adopting the axioms described in sect . [ sec : igimf ] . the formal solution of the closed - box equations [ eq : cbsystem ] is now : @xmath77y_z(\tau)d\tau , \notag\\    & m_g(t)=m_g(0)\exp{\left[-s\int_0^t \left[1-r(\tau)\right]d\tau\right ] } ,    \notag\\    & \psi(t)=sm_g(t)=sm_g(0)\exp{\left[-s\int_0^t      \left[1-r(\tau)\right]d\tau\right]}. \label{eq : sol1}\end{aligned}\ ] ] we have included also the solution for @xmath78 to show that this system of equations is indeed implicit ( the integrals in the right hand sides depend on the quantities on the left hand sides ) . to solve these equations thus , an implicit procedure must be employed . examples of the calculation of the yields @xmath74 using the igimf theory and the axioms described in sect . [ sec : igimf ] are shown in fig . [ fig : y_z ] . in particular , the left panel refers to a model ( model 1 ) in which the axioms 16 of sect . [ sec : igimf ] are adopted . the right panel instead refers to model 2 , i.e. a model where also the last axiom of sect . [ sec : igimf ] is implemented . as mentioned in the introduction , the methodology described in this paper is not limited to the igimf theory . other theories of variable imf could be used to calculate @xmath74 and @xmath37 as a function of galactic properties . once @xmath74 and @xmath37 have been calculated , the rest of the theory is general and not related to a specific imf formulation .    to construct these functions @xmath79 , we have adopted the yields of woosley & weaver ( 1995 ) for massive stars and the ones of van den hoek & groenewegen ( 1997 ) for intermediate - mass stars . both sets of yields are metallicity - dependent . since we are interested in the early phases of the evolution of galaxies , we have thus taken the set of yields with z=0.001 from van den hoek & groenewegen ( 1997 ) and the one with z=10@xmath80 z@xmath5 from woosley & weaver ( 1995 ) . we have checked that our results depend little on the chosen metallicity , although larger initial metallicities tend to produce slightly higher yields . the increase of @xmath74 with sfr is due to the fact that the igimf is flatter for high sfrs , thus more metals are produced . because of eq . [ eq : alpha3 ] ( and because of the dependence of @xmath81 on [ fe / h ] ) , the yield increases also with decreasing metallicity . the yields of model 2 are slightly larger than the yields of model 1 , because of the axiom 7 of sect . [ sec : igimf ] . the differences between model 1 and model 2 however disappear for sfrs smaller than 1 m@xmath5 yr@xmath22 because of eq . [ eq : betad ] . given the small differences between model 1 and model 2 , we will focus from now on on results pertaining model 1 . notice also that the yields tend to flatten out at very low sfrs . this is due to the fact that low- and intermediate - mass agb stars still produce some metals according to the yields of van den hoek & groenewegen and the fraction of agb stars ( in particular the fraction of low - mass agb stars ) changes very little with the sfr . the function @xmath82 is qualitatively similar to the shown function @xmath83 , although the magnitude of the variations is smaller . [ fig : traces ] depicts the time evolution of four representative model galaxies with different initial gaseous masses , ranging from 10@xmath84 m@xmath5 ( lower curves ) to 10@xmath85 m@xmath5 ( upper curves ) . the left panel refers to the model 1 ( axioms 16 of sect . [ sec : igimf ] ) and the right panel refers to model 2 ( where also the last axiom of sect . [ sec : igimf ] is implemented ) . in these models it is assumed that all the mass initially is in the form of gas . clearly , the metallicity evolution depends now , at variance with the invariant imf case , on the model galaxy , as larger galaxies have on average higher sfr , hence flatter imfs , with a higher production rate of heavy elements . as a comparison , a model with an invariant imf is shown in fig . [ fig : traces ] , too . typical parameters ( @xmath86 and @xmath87 ) have been used to plot these curves . it is true that low - metallicity models tend to have flatter imfs , too ( see eq . [ eq : alpha3 ] and the definition of the parameter @xmath81 ) , but since the proportionality constant in front of [ fe / h ] in the definition of @xmath81 is small , the metallicity effect on the imf is outweighted by the sfr effect . due to the larger yields , the final metallicities for model 2 are ( slightly ) larger than the ones for model 1 .    for all these model galaxies , we have also calculated the average stellar metallicities , according to eq . [ eq : zstar ] . in principle , these stellar metallicities could be compared with the best observational data available to us , i.e. the stellar metallicities of dwarf satellites of the milky way and in general of dwarf galaxies of the local group . the stellar metallicities of these galaxies are well determined and they build a clear and tight mz relation ( see e.g. kirby et al . unfortunately , the [ fe / h ] abundance is usually measured in stars ( see again kirby et al . this makes the comparison with simple model predictions much less straightforward . in fact , as it is well known , the majority of iron is produced by type ia supernovae , over long timescales . as recalled in sect . [ sec : simple ] , the simple models give acceptable results only for elements ( such as the @xmath68-elements ) produced on short timescales by massive stars . hence , simple models adopting the instantaneous recycling approximation can not predict the right fe abundance since they can not take into account properly the contribution from type ia sne . we are working on an extension of the simple model formalism that takes into account also iron produced by sneia . at the moment , results of simple models can be safely compared to gas - phase abundances of @xmath68-elements only , and this is what we do in fig . [ fig : mz ] . we use the data of lee et al . ( 2006 ) and we compare them with the results of model galaxies shown in fig . [ fig : traces ] . the initial metallicity considered here ( the @xmath88 in eq . [ eq : sol1 ] ) is zero . initial metallicities larger than zero are important when studying tdgs , therefore they will be considered extensively in paper ii . we anticipate however here that , if tdgs form at very large redshifts out of galaxies whose chemical enrichment is still limited , their initial metallicity is not large and @xmath89 is a reasonable assumption . clearly , we can see a mz in fig . [ fig : mz ] , in the sense that low - mass models attain smaller stellar metallicities than high - mass ones . however , the comparison with observations is not particularly good . in particular at low masses the theoretical predictions overestimate the observations by 0.30.4 dex . this is due in part to the inherent uncertainties in the lee et al . ( 2006 ) mz relation , in part to the fact that galactic flows in low - mass systems can not be neglected . we analyze in detain the effect of galactic flows in the following section . we revise now the system of equations [ eq : cbsystem ] to take into account possible gas inflows and outflows in galaxies . as already mentioned , this will be done by introducing parameters @xmath44 and @xmath45 that relate the flow rates to the sfr ( see eq . [ eq : flows ] ) . additionally , we consider the possibility that the outflows are metal - enriched ( _ differential winds _ ) , i.e. that the metallicity of the outflowing gas is larger than the average gas - phase metallicity in the galaxy . both observations ( martin et al . 2002 ; ott et al . 2005 ) and numerical models ( mac - low & ferrara 1999 ; recchi & hensler 2013 ) support this assumption . notice that probably also the efficiency with which the various chemical species are ejected vary from element to element ( _ selective winds _ ; see recchi et al . 2001 ; 2008 ) but for simplicity we will not consider further this hypothesis .    assuming that the metallicity of the outflowing gas is @xmath90 , with @xmath91 , we can write :    @xmath92 , \cr      \psi ( t ) = s m_g ( t ) . \end{cases } \label{eq : lbsystem}\ ] ]    the solution of this system of equations for @xmath71 and @xmath70 is given by : @xmath93d\tau\right ] } ,    \notag\\    & z(t)=z_0+\frac{\int_0^t[1-r(\tau)]s[\lambda z_a + y_z(\tau ) ]      i(\tau)d\tau}{i(t ) } , \notag\\    & i(t)=\exp{\left[-\int_0^t[1-r(\tau)]s[\lambda(\gamma-1)+\lambda ]        d\tau . \right ] }    \label{eq : sol2}\end{aligned}\ ] ] here , @xmath94 is an appropriate integration factor . notice that the expression for @xmath70 reduces to the analogous closed - box expression shown in eq . [ eq : sol1 ] in the case @xmath95 , @xmath96 ( i.e. no differential winds , no infall ) . as already mentioned , it is to be expected that low - mass galaxies are much more affected by galactic winds than are more - massive ones , i.e. we expect @xmath45 to decrease as the galactic mass increases . a simple functional form satisfying this condition is to decrease with increasing total ( baryons plus dark matter ) galactic mass . if we assume for simplicity a constant ratio between dark matter and baryonic mass , a correlation between @xmath45 and @xmath97 as in eq . [ eq : mglambda ] is recovered . ] : @xmath98 this formula is normalized in order to reproduce the outflow rate to sfr ratio in tdg simulations in which the igimf prescriptions are adopted ( ploeckinger et al . according to this formula , @xmath45 tends to be quite large for very low mass systems ( it is @xmath99 for @xmath100 m@xmath5 ) but it tends to zero for very high mass ones ( it is @xmath101 for @xmath102 m@xmath5 ) . on the other hand , we do not expect the infall to significantly depend on the mass of the parent galaxy , therefore we take , in compliance with spitoni et al . ( 2010 ) , @xmath103 for all model galaxies . notice that the values of @xmath45 are considerably smaller than the ones needed in order to reproduce the mz relation based on an invariant imf . the comparison between our adopted values of @xmath45 and the ones of spitoni et al . ( 2010 ) for @xmath103 are given in fig . [ fig : lambda_comp ] . notice that we are interested in this study in smaller systems compared to the ones considered by spitoni et al . ( 2010 ) , therefore the two curves in fig . [ fig : lambda_comp ] do not fully overlap . it is nevertheless clear that the outflow rates in spitoni et al . ( 2010 ) are much larger than the ones adopted in this study . this is due to the fact that the reduction of the metallicity at low masses is , according to the igimf theory , largely due to the steepening of the imf slope at large stellar masses in galaxies with small sfrs . moreover , we assume @xmath104 , based on the results of the chemo - dynamical simulations of recchi & hensler ( 2013 ) .    the mz relation obtained with leaky box models with this dependence of @xmath45 with @xmath97 is shown in fig . [ fig : mz_w ] . the fit with observations is very good . it could be better if we better fine - tune the relation @xmath105 but this is not the main aim of the paper . the igimf theory ( kroupa et al . 2013 ) predicts a coupling between some properties of a galaxy ( the sfr and the metallicity ) and the imf . since the imf in turn strongly affects the dynamical evolution of the galaxy , the feedback between galaxy evolution and imf is difficult and the fully complexity of a variable imf has not been yet included in hydrodynamical simulations ( but see bekki 2013 ; ploeckinger et al . 2014 ; recchi 2014 ) . even the here treated approach based on the so - called `` simple model of chemical evolution '' leads to implicit integral equations that must be solved iteratively . we note in passing that the dependence of the imf on the metallicity is well established theoretically , as low - metallicity gas cools less efficiently and self - gravitating clumps are more resilient to fragmentation . this leads to the formation of dense cores of higher mass . for this reason , the imf is supposed to be extremely biased towards massive stars if the metallicity is smaller than @xmath106 10@xmath80 z@xmath5 ( schneider et al . also the dependence of the imf on the sfr is nowadays observationally well established ( hoversten & glazebrook 2008 ; meurer et al . 2009 ; gunawardhana et al . 2011 ; kroupa 2014 ) , and the igimf theory is the only existing computable access accounting for these observations ( weidner et al . 2013 ) . thus , in spite of the inherent complexity , detailed simulations of galaxy evolution based on variable imfs need to be performed . with this paper , we aimed at showing in a simple setting how to take into account imf variations in models of the chemical evolution of galaxies ( see also martinelli & matteucci 2000 ; calura et al . the method outlined here can be used in order to extend analytical studies of the evolution of galaxies based on simple models of chemical evolution ( e.g. spitoni et al . 2010 ; lilly et al . 2013 ; pipino et al . we described how the solution of the simple model differs from the standard , textbook analytical solutions of the simple model . we also showed ( building on previous results of kppen et al . 2007 ) that the igimf theory naturally leads to a mass - metallicity relation . in fact , low - mass galaxies are characterized on average by smaller sfrs . according to the igimf theory , a low sfr leads to a steep , top - light imf , in which the production of heavy elements by massive stars is extremely limited . more massive galaxies instead produce many more massive stars because of the higher level of sfr , hence the attained present - day metallicity is larger . we outline here once again that the main aim of the paper was not to provide the ultimate explanation for the mass - metallicity relation . many detailed theoretical studies have been performed on this subject ( see e.g. matteucci & chiosi 1983 ; tremonti et al . 2004 ; finlator & dav 2008 ; calura et al . 2009 ; peeples & shankar 2011 ; yin et al . 2011 among many others ) . as shown by fig . [ fig : mz ] , the mz relation obtained with the igimf theory does not fit well the data and assumptions about flow rates are required in order to obtain a good fit . nevertheless it is important to show that @xmath107 the igimf theory ( as well as any other theory according to which the imf is top - heavier in more massive systems ) naturally accounts for the mz relation ; @xmath108 even within the igimf formalism , flow rates ( infall and outflows ) remain crucial ingredients for the evolution of galaxies ; @xmath109 the intensity of the galactic wind must be assumed to be inversely proportional to the initial gaseous mass of the model galaxy , in accordance with observations and theoretical studies . the outflow rates are however considerably smaller than the ones needed to reproduce the mz relation in models with invariant imfs such as the one of spitoni et al . ( 2010 ) .    in a companion paper this methodology will be applied to the study of the mass - metallicity relation in dark - matter - free tidal dwarf galaxies . at the same time , efforts are ongoing to model by means of hydrodynamical codes the full chemo - dynamical evolution of galaxies within the igimf theory , in order to assess the reliability of the results obtained in this paper . we thank an anonymous referee for suggestions and constructive criticisms which considerably improved the quality of the paper . we thank e. spitoni for providing us the data reported in fig . [ fig : lambda_comp ]
Abstract:

Answer: standard analytical chemical evolution modelling of galaxies has been assuming the stellar initial mass function ( imf ) to be invariant and fully sampled allowing fractions of massive stars to contribute even in dwarf galaxies with very low star formation rates ( sfrs ) . recent observations show the integrated galactic initial mass function ( igimf ) of stars , i.e. the galaxy - wide imf , to become systematically top - heavy with increasing sfr . this has been predicted by the igimf theory , which is here used to develop the analytical theory of the chemical evolution of galaxies . this theory is non - linear and requires the iterative solution of implicit integral equations due to the dependence of the igimf on the metallicity and on the sfr . it is shown that the mass  metallicity relation of galaxies emerges naturally , although at low masses the theoretical predictions overestimate the observations by 0.30.4 dex . a good agreement with the observation can be obtained only if gas flows are taken into account . in particular , we are able to reproduce the mass  metallicity relation observed by lee et al . ( 2006 ) with modest amounts of infall and with an outflow rate which decreases as a function of the galactic mass . the outflow rates required to fit the data are considerably smaller than required in models with invariant imfs .    stars : abundances  stars : luminosity function , mass function supernovae : general  galaxies : evolution  galaxies : dwarf  galaxies : star clusters : general