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: in field theory , currents which are conserved classically may not be quantum mechanically  @xcite . for example , in massless @xmath2 the conservation of the axial current is violated by the axial anomaly . the adler - bardeen theorem states that with the proper regularization scheme , the coefficient of the anomaly , as computed at one loop order , is exact to all orders in perturbation theory . moreover , the axial anomaly does not change if the fermions propagate in either a thermal bath or a fermi sea  @xcite .    in vacuum , one of the most striking manifestations of the axial anomaly is the decay of a neutral pion into two photons , as the amplitude is directly proportional to the coefficient of the axial anomaly in qed  @xcite . a natural supposition is then that because the axial anomaly does not change with temperature , neither does the amplitude for @xmath3  @xcite .    in this paper we show that the story is more involved . we compute with a gauged nonlinear sigma model  @xcite which properly incorporates all anomalies by inclusion of the wess - zumino - witten ( wzw ) term  @xcite . the effective lagrangian for @xmath3 is @xmath4 where @xmath5 is the pion decay constant , @xmath6 is the number of colors , _ etc . _    in sec . [ sec : vac ] we start by computing the effects of pion loops on the amplitude of ( [ eq:1.1 ] ) , using the wzw action to one loop order in vacuum  @xcite . the form of the wzw action is constrained by topology  @xcite , so after the dust of calculation settles , in vacuum the result is trivial : the only effect of the pion loops is to change a bare pion decay constant into a renormalized @xmath7 . we then extend the calculations to soft , cool pions at low temperature  @xcite . in this paper we work exclusively with two flavors in the chiral limit . the restriction to two flavors is done for ease of calculation , and is otherwise inessential . the chiral limit , @xmath8 , is assumed because then the pion decay amplitude is directly related to the axial anomaly ; for calculations at @xmath9 at nonzero temperature , see  @xcite . we believe that our results are relevant for @xmath9 ( as in vacuum ) , but more detailed analysis is required to establish this .    at nonzero temperature , calculations in a background field formalism  @xcite show that to @xmath10 , the zero temperature pion decay constant is replaced by a temperature dependent form  @xcite , @xmath11 thus a second guess for the change of @xmath12 with temperature would be that the zero temperature @xmath7 is replaced by @xmath13 . since @xmath13 decreases to @xmath10 , if true the amplitude , @xmath14 , would increase to this order .    in sec [ sec : lowt ] we evaluate precisely the same diagrams as at zero temperature to @xmath10 . the result , ( [ eq:3.9 ] ) , is the sum of two terms : one has exactly the form of ( [ eq:1.1 ] ) , with @xmath7 replaced by @xmath13 , but there is also a second term , special to nonzero temperature . this type of term was derived recently in nonlinear sigma models in the absence of gauge fields  @xcite : it is nonlocal , analogous to the hard thermal loops of hot gauge theories  @xcite . for @xmath3 , to order @xmath10 we find that the sum of these two terms is such that instead of increasing , like @xmath14 , the amplitude decreases , like @xmath15 , ( [ eq:3.8a ] ) .    in sec . [ sec : ward ] we give a general analysis of the relationship between the chiral ward identities and the amplitude for @xmath3 . as is standard @xcite , we use the anomalous ward identity to relate a three point function of currents to the amplitude for @xmath3 . at zero temperature , this relationship is precise because of the sutherland - veltman theorem @xcite ( in a slight abuse of terminology ) . the proof of the sutherland - veltman theorem depends crucially upon lorentz invariance . a heat bath , however , provides a preferred rest frame ; extending an analysis of itoyama and mueller  @xcite , we show that consequently , the sutherland - veltman theorem does not apply at any nonzero temperature . this is why @xmath3 changes with temperature , even though the anomaly does nt : besides the contributions from @xmath3 , because there is no sutherland - veltman theorem at nonzero temperature , there are other terms which enter to ensure that the adler - bardeen theorem is satisfied .    in sec . [ sec : adlerbardeen ] we demonstrate these general arguments by computing the correlator between one axial and two vector currents in the nonlinear sigma model to one loop order , @xmath10 . at nonzero temperature , new tensor structures arise in this correlator ; these structures are why the sutherland - veltman theorem is inapplicable at @xmath16 . nevertheless , when all terms are added together , we find that the adler - bardeen theorem remains valid to one loop order . this is a useful and nontrivial check of our result for @xmath3 . technical details are relegated to several appendices . the wzw action is discussed in appendix  [ awzw ] . various formulas for hard thermal loops are collected in appendix  [ ahtl ] . we use the imaginary time formalism at nonzero temperature in this paper , but show in appendix  [ sec : realtime ] how the same results follow in the real time formalism . lastly , in appendix  [ gamma3pi ] we compute another anomalous amplitude , that for @xmath17  @xcite , at low temperature .    while the principal concern of our work are thermal field theories , we hope that some of our dicussion , especially that in sec . [ sec : ward ] , might be of more general interest . perhaps understanding why @xmath3 is not tied to the axial anomaly at nonzero temperature helps us better understand this relation at zero temperature . we start by computing the effects of pion loops on the amplitude for @xmath18 in vacuum . at one loop order , the relevant diagrams are figs . ( 1.a)-(1.d ) : fig . ( 1.a ) gives the pion field renormalization constant , @xmath19 ; fig . ( 1.b ) , the renormalized pion decay constant @xmath7 , while corrections to the amplitude itself are given by figs . ( 1.c ) and ( 1.d ) .    for the `` tadpole '' type diagrams of figs . ( 1.a ) , ( 1.b ) , and ( 1.c ) we use a trick . to compute fig . ( 1.a ) we expand the full lagrangian , ( [ eq : a.10 ] ) , to quartic order in the pion field , @xmath20 + \ldots \end{aligned}\ ] ] in all expressions from the appendix , we need to use the bare pion decay constant , @xmath21 , instead of the renormalized quantity @xmath7 , as one loop effects change @xmath21 into @xmath7 . for the quartic terms , contracting two out of the four pion fields in all possible ways gives @xmath22 where @xmath23 while this integral is quadratically divergent , we ignore regularization , since its actual value is irrelevant for our purposes . in ( [ eq:2.4 ] ) @xmath24 is the renormalized pion field , and so @xmath25 for the pion decay constant , instead of fig . ( 1.b ) we expand the axial current @xmath26 of ( [ eq : a.11 ] ) to cubic order in the pion field , @xmath27 contracting all pairs of pion fields , @xmath28 so that @xmath29 as is typical of nonlinear sigma models , unphysical , off - shell quantities such as @xmath19 , ( [ eq:2.5a ] ) , depend upon the parametrization of the coset space , while physical expressions , such as that for @xmath7 in ( [ eq:2.8 ] ) , do not .    turning to the amplitude for @xmath3 , from  ( [ eq:1.1 ] ) it is @xmath30 @xmath31 and @xmath32 , and @xmath33 and @xmath34 are the momenta and the polarization vectors of the two photons , both of which lie on the mass shell , @xmath35 . at tree level , the bare coupling @xmath36 satisfies @xmath37 the right hand side of ( [ eq:2.10 ] ) is precisely the coefficient of the axial anomaly in qed  @xcite .    to evaluate fig . ( 1.c ) we expand the anomalous current for the coupling to two photons , ( [ eq : a.7 ] ) , ( [ eq : a.11 ] ) , and ( [ eq : a.14 ] ) , to cubic order in the pion field , @xmath38 + \ldots \end{aligned}\ ] ] contracting two pion fields , @xmath39 the two distinct diagrams of fig . ( 1.d ) require more effort : @xmath40 where @xmath41 this peculiar separation of terms in @xmath42 is done in anticipation of the results at nonzero temperature . in vacuum , for @xmath43 , @xmath44 , and so because of the levi - civita symbol , @xmath45 does not contribute to @xmath46 . hence ( [ eq:2.13 ] ) reduces to a form proportional to the original term in ( [ eq:2.9 ] ) . putting everything together , the renormalized coupling @xmath47 equals @xmath48 the @xmath49 comes from a factor of @xmath50 for the renormalized pion field , fig . ( 1.a ) and ( [ eq:2.5a ] ) , the @xmath51 from fig . ( 1.c ) , ( [ eq:2.12 ] ) , and the @xmath52 from fig . ( 1.d ) , ( [ eq:2.13 ] ) and ( [ eq:2.14 ] ) . hence at one loop order , @xmath53 comparing ( [ eq:2.10 ] ) and ( [ eq:2.17 ] ) , we see that the anomaly is not renormalized to one loop order  @xcite : seperate divergences in @xmath7 and @xmath47 cancel in the product  @xcite . we now compute the decay for a cool pion , at a temperature @xmath54  @xcite . the diagrams are identical , the only difference is that we need to compute at @xmath55 . for the tadpole diagrams of figs . ( 1.a ) , ( 1.b ) , and ( 1.c ) , the integral is the analogy of ( [ eq:2.5 ] ) : @xmath56 we use the imaginary time formalism , @xmath57 , @xmath58 , and @xmath59 for integral @xmath60 ; after doing the sum over @xmath60 , the bose - einstein statistical distribution function , @xmath61 , appears : @xmath62 @xmath63 is the value of the integral at zero temperature . henceforth we drop the @xmath64 part of any integrals , assuming that they turn bare into renormalized quantities , such as @xmath21 into @xmath7 , throughout . the calculation of temperature dependent corrections to the pion decay constant proceeds as in the previous section . ignoring ultraviolet renormalization , in ( [ eq:2.8 ] ) we replace @xmath63 by @xmath65 , and @xmath21 by @xmath7 , to obtain @xmath66 which was quoted in the introduction , ( [ eq:1.2 ] ) . thus if anything unusual happens at nonzero temperature , it can only be from the diagram of fig . ( 1.d ) . unlike the tadpole diagrams , this diagram has nontrivial momentum dependence , and so we must be more precise in specifying the external momenta . to compute scattering in a thermal bath , we continue the euclidean momenta @xmath67 to a minkowski energy @xmath68 by @xmath69 . following  @xcite we further assume that each momentum is not only cool but soft , taking both @xmath70 . for scattering between soft , cool pions , in  @xcite we showed that the leading temperature corrections are directly analogous to the hard thermal loops of hot gauge theories  @xcite . we used the background field method , but only in the absence of external gauge fields . while the perturbative calculations which follow are thus less elegant , they illustrate the physics more directly . from the perspective of @xcite , there is nothing special about @xmath3 ; the connection to the axial anomaly will be clarified later . we introduce @xmath71 , @xmath72 with @xmath73 the @xmath74 sign denotes that only the hard thermal loops in the integral are retained , which we denote by @xmath75 and @xmath76 . the hard thermal loops are the terms @xmath77 , and are given explicitly in appendix [ ahtl ] . up to an overall constant , @xmath78 is the same hard thermal loop as appears in the polarization tensor for a photon in thermal equilibrium . for a thermal photon , the screening of time dependent electric and magnetic fields implies that the mass shell is at @xmath79 . for the sake of simplicity we assume that the photons do not thermalize ; then the only photons which propagate are transverse modes on the light cone , @xmath80 . the polarization vectors for these modes are purely spatial vectors which satisfy @xmath81 .    from ( [ eq:2.13 ] ) , the contribution to the amplitude from fig . ( 1.d ) , @xmath46 , involves @xmath82 for one of the two photons on their mass shell . using ( [ eq : b.5 ] ) of appendix  [ ahtl ] , @xmath83 where @xmath84 is the polarization vector for the photon with momentum @xmath85 . only the first term on the right hand side of ( [ eq : b.5 ] ) contributes , as terms in @xmath86 which are @xmath87 or @xmath88 drop out after contraction with @xmath89 . from the definition of @xmath90 , ( [ eq:3.3 ] ) , @xmath91 consequently , while at zero temperature fig . ( 1.d ) contributes to the amplitude for @xmath3 , to leading order at nonzero temperature its contribution vanishes identically , ( [ eq:3.7 ] ) .    knowing that fig . ( 1.d ) does nt contribute , it is then easy to read off the one loop corrections to the coupling @xmath47 to @xmath10 , @xmath92 . as in ( [ eq:3.2a ] ) , we start with ( [ eq:2.16 ] ) , and then replace @xmath63 by @xmath65 , and @xmath21 by @xmath7 . we keep the @xmath49 from fig . ( 1.a ) , the the @xmath51 from fig . ( 1.c ) , but replace the @xmath93 from fig . ( 1.d ) by @xmath94 , to obtain @xmath95 notice that while ( [ eq:2.8 ] ) is precisely analogous to ( [ eq:3.2a ] ) , because of the difference in fig . ( 1.d ) , ( [ eq:2.16 ] ) is not analogous to ( [ eq:3.8 ] )  there is a difference in sign . consequently , @xmath96 as discussed in the introduction , naively one might guess that ( [ eq:2.17 ] ) generalizes to nonzero temperature just by replacing @xmath7 and @xmath47 with @xmath13 and @xmath92 , respectively . this is wrong : to @xmath10 , instead of @xmath97 , as would be guessed from ( [ eq:2.17 ] ) , instead @xmath98 , ( [ eq:3.8a ] ) . we do not know why , to leading order about low temperature , @xmath92 decreases in exactly the same manner as @xmath13 . our result in ( [ eq:3.8 ] ) differs from that found by dobado , alvarez - estrada , and gomez  @xcite . these authors consider the same model , but find @xmath99 to @xmath10 . our results agree except for fig . ( 1.d ) , which we believe was treated incorrectly  @xcite .    before continuing , following @xcite we construct the effective lagrangian for @xmath3 to @xmath10 . in @xmath100 of ( [ eq:3.3 ] ) , the term @xmath101 is easy to include . at zero temperature , ( [ eq:2.14 ] ) , this term is the only part of fig . ( 1.d ) which contributes , @xmath93 in ( [ eq:2.16 ] ) , and turns @xmath102 in @xmath47 into @xmath103 in @xmath104 . thus to @xmath10 , the effect of figs . ( 1.a ) , ( 1.c ) , and the term @xmath101 in @xmath100 of fig . ( 1.d ) is just to change @xmath103 into @xmath105 in the original lagrangian , @xmath12 of ( [ eq:1.1 ] ) . including the term @xmath76 is less trivial . because of ( [ eq : b.1 ] ) , it must be constructed out of transverse quantities . using ( [ eq : b.7 ] ) , we find that to @xmath10 , the effective lagrangian for @xmath3 is @xmath106 @xmath107 in this expression @xmath108 , @xmath109 and @xmath110 the vector @xmath111 and the integration over the angle @xmath112 arediscussed following ( [ eq : b.2 ] ) . the non - local term in  ( [ eq:3.9 ] ) is specific to finite temperature . at zero temperature , there is no other term besides  ( [ eq:1.1 ] ) that contributes to @xmath113 for photons on the mass - shell . in the terminology of the nonlinear sigma model @xcite , the operator of  ( [ eq:1.1 ] ) is @xmath114 , while operators at next to leading order are @xmath115 . these operators , however , are either proportional to @xmath116 or @xmath117 , and so vanish on the photon(s ) mass shell , or @xmath118 , and so vanish in the chiral limit . thus in the vacuum , the only possible change in ( [ eq:1.1 ] ) is the transmutation of @xmath21 into @xmath7 , sec . [ sec : vac ] . at nonzero temperature , however , there are new nonlocal terms which arise ,  ( [ eq:3.9 ] ) . because they are nonlocal , these new terms are also @xmath114 , and so as important as ( [ eq:1.1])@xcite . this is the technical reason why the amplitude for @xmath119 depends nontrivially upon temperature . in the previous section we found that the amplitude for @xmath1 diminishes to leading order in an expansion about low temperature , ( [ eq:3.8a ] ) . the question we address in this section is why is this amplitude tied to the coefficient of the axial anomaly at zero temperature , ( [ eq:2.17 ] ) , but not at nonzero temperature ? we work in the chiral limit to leading order about zero temperature , @xmath10 , because then we can make certain technical assumptions which simplify the discussion . the general case is considered at the end of this section . define the vector current , @xmath120 , and the axial current in the isospin-@xmath121 direction , @xmath122 . the vector curent is conserved , @xmath123 while the axial current is anomalous , @xmath124 by the adler - bardeen theorem , the coefficient of the right hand side is exact to one loop order  @xcite , and is independent of temperature and density  @xcite . one quantity we can compute is the ( thermal ) three point green s function between two vector , and one axial vector , current : @xmath125 where @xmath126 is the hamiltonian . then @xmath127 satisfies current conservation , @xmath128 and the anomalous ward identity , @xmath129 @xmath130 . to relate the anomalous ward identity to the amplitude for pion decay we follow shore and veneziano  @xcite . at low temperature the pion couples to the axial current as @xmath131 ( this is not valid to @xmath132 ; then the relation is more complicated , @xcite . )    to obtain the amplitude for @xmath133 , we introduce @xmath134 times the matrix element between two @xmath2 currents and a pion , @xmath135 this is related to the pion decay amplitude , ( [ eq:2.9 ] ) , as @xmath136 subtracting the one pion pole term from ( [ eq:4.1 ] ) gives @xmath137 , which by construction is one pion irreducible , @xmath138 again , by current conservation @xmath139 while the anomalous ward identity , ( [ eq:4.2 ] ) , becomes @xmath140 the trick is now to try to deduce general relations using just ( [ eq:4.9 ] ) , ( [ eq:4.8 ] ) , and bose symmetry between the two photons , @xmath141 .    we first discuss zero temperature , where we can invoke euclidean invariance . the most general pseudo - tensor @xmath137 which satisfies all of our conditions involves three terms : @xmath142 current conservation , ( [ eq:4.9 ] ) , gives @xmath143 while from the anomalous ward identity , ( [ eq:4.8 ] ) , @xmath144 combining these two relations we obtain @xmath145 putting the photons on their mass shell @xmath146 , the left hand side in  ( [ eq:4.13 ] ) reduces to @xmath147 . this is zero on the pion mass shell , @xmath148 , since by definition @xmath149 is one pion irreducible , and so can not have a pole @xmath150 . hence the left hand side of ( [ eq:4.13 ] ) vanishes , and we obtain the desired relation between @xmath47 and the coefficient of the axial anomaly , @xmath151 which is  ( [ eq:2.14 ] ) . this analysis , and especially the tensor decomposition of ( [ eq:4.10 ] ) , is identical to the derivation of the sutherland - veltman theorem  @xcite . historically , this theorem predated the anomaly , and was used to conclude that @xmath152 . by adding the axial anomaly through the anomalous ward identity of ( [ eq:4.8 ] ) , however , we obtain @xmath153 , ( [ eq:4.14 ] ) . from a modern perspective , then , it is precisely the sutherland - veltman theorem which relates the axial anomaly to the amplitude for @xmath3 , and tells us at zero temperature the left hand side of ( [ eq:4.14 ] ) vanishes . this is no longer true at nonzero temperature . following itoyama and mueller  @xcite , we write the most general tensor decomposition of @xmath154 . in a thermal bath , however , euclidean symmetry is lost , and the rest frame of the thermal bath , which we take as @xmath155 , enters . some of the possible tensors include @xmath156 we have only included the terms in @xmath149 which contribute to the ward identities of ( [ eq:4.9 ] ) and ( [ eq:4.8 ] ) . current conservation gives @xmath157 while the anomalous ward identity fixes    @xmath158    from which follows @xmath159 @xmath160 putting all fields on their mass shell , @xmath80 , @xmath148 , and assuming as before that @xmath161 has no pole @xmath150 , we obtain @xmath162 since the terms on the left hand side involve only @xmath163 and @xmath164 , there is no reason for them to vanish even if @xmath165 ; thus the direct connection between @xmath92 , @xmath13 , and the anomaly is lost . we stress that , as always , the adler - bardeen theorem remains valid , and gives ( [ eq:4.18a ] ) . it is only the sutherland - veltman theorem which no longer applies at nonzero temperature . the above analysis only applies to leading order in low temperature , @xmath10 . this is because beyond leading order , the pion mass shell is no longer at @xmath166 @xcite ; also , as photons thermalize , their mass shell moves off the light cone this is incorporated by using ( [ eq:4.18a ] ) instead of ( [ eq:4.18 ] ) . for instance , we can understand how the anomalous ward identity is satisfied in a chirally symmetric phase . from explicit calculation in a constituent quark model @xcite , @xmath3 vanishes once chiral symmetry is restored . this does not conflict with the anomalous ward identity since even if @xmath167 , there are other terms which can ensure that ( [ eq:4.18a ] ) is satisfied . for photons which do not thermalize , so @xmath80 , even assuming that at the chiral critical point that the pion mass shell is @xmath166 , the tensors @xmath168 and/or @xmath169 will in general be nonzero . our analysis agrees with the results of contreras and loewe  @xcite , who computed the triangle diagram at @xmath55 with massive fermions . their result , ( 1.2 ) , is directly related to @xmath137 . there are two terms : the first is regular , temperature dependent , and gives the amplitude for @xmath119 , while the second is the anomaly , and is independent of temperature .    at first sight it might appear peculiar that the sutherland - veltman theorem applies at zero temperature , but fails at any nonzero temperature . even at zero temperature , however , the sutherland - veltman theorem only applies in the chiral limit when both photons are on their mass shell . without all of these conditions , the left hand side of ( 4.13 ) does not vanish , and @xmath47 is not given by ( 4.14 ) . an example of this occurs when one ( or both ) of the photons are off the mass shell @xcite , even if @xmath148 . in particular , in the limit of large @xmath170 , it is known that @xmath171 @xcite . thus on the right hand side of ( 4.13 ) , we can neglect @xmath172 relative to the anomaly term , @xmath173 , and ( 4.13 ) tells us that a combination of @xmath174 and @xmath161 are given entirely by the anomaly , @xmath175 . in this section we calculate the correlator @xmath176 to @xmath10 in the nonlinear sigma model . this allows us to check explicitly the temperature dependence of @xmath13 and @xmath92 found previously in sec . [ sec : lowt ] . it is also illuminating to see exactly which amplitudes enter at one loop order in the nonlinear sigma model , and how they conspire to satisfy the adler - bardeen theorem .    at tree level in the nonlinear sigma model , there is no direct coupling between two vector and one axial vector currents , so @xmath177 . thus the anomalous ward identity , ( [ eq:4.2 ] ) , is satisfied entirely by the one pion reducible term , ( [ eq:4.8 ] ) , with @xmath178 . this is illustrated in fig . contributions to @xmath179 are generated at one loop order . in order to compute these , it is necessary to compute corrections to the axial current . we have computed these in two ways . the most direct is to follow the original method of wess and zumino @xcite . besides the the photon field @xmath180 , which couples to the electromagnetic current @xmath120 , we also introduce an external field @xmath181 which couples to the axial current @xmath182 . one then differentiates the generating functional of wess and zumino with respect to the external fields @xmath180 and @xmath181 . at one loop order , the terms required are @xmath183 and @xmath184 ; @xmath185 are the charged pion fields . after lengthy calculation , we find @xmath186 @xmath187 to leading order in @xmath188 , this @xmath189 is invariant under an electromagnetic gauge transformation , @xmath190 the terms in @xmath189 can be viewed as corrections to the axial current . in appendix a we compute these corrections using the noether construction of the axial current . this is somewhat delicate , since it is necessary to start with a lagrangian density ( and not merely a lagrangian ) which is manifestly gauge invariant . the result , ( [ eq : a.19 ] ) , agrees with ( [ eq:5.1 ] ) .    through the diagrams of fig . ( 3 ) , the couplings in @xmath191 generate the one pion irreducible terms in @xmath154 . we begin with the results at zero temperature . given the couplings in @xmath192 , the only tensor structure which arises from fig . ( 3 ) is @xmath193 . comparing with the tensors in  ( [ eq:4.10 ] ) , then , at one loop order automatically @xmath194 . for @xmath195 , we find @xmath196 the contribution @xmath197 is from fig . ( 3.a ) , that @xmath198 from the two diagrams of fig . ( 3.b ) . because all of the @xmath199 s vanish , current conservation  ( [ eq:4.12 ] ) and the anomaly equation  ( [ eq:4.11 ] ) , are satisfied rather trivially . that the latter vanishes is equivalent to the sutherland - veltman theorem ,  ( [ eq:4.14 ] ) this explains our results in sec . [ sec : vac ] , where we found that at zero temperature , the renormalization of @xmath7 and @xmath47 exactly compensate each other .    at nonzero temperature we can evaluate @xmath200 using the results of sec . [ sec : lowt ] . ( 3.a ) is a tadpole diagram , and so as at zero temperature , just contributes to @xmath195 . thus the only possibility for a new tensor structure is from the hard thermal loop in fig . for one ordering of momenta , fig . ( 3.b ) gives @xmath201 where @xmath202 is given in  ( [ eq:3.3 ] ) . putting the photon momentum @xmath31 on its mass shell , we require the result for @xmath203 in  ( [ eq : b.5 ] ) . we drop all terms in @xmath204 , ( [ eq : b.5 ] ) , since they vanish upon contraction with the photon polarization tensor . the term @xmath205 in  ( [ eq : b.5 ] ) cancels against the same term in  ( [ eq:3.3 ] ) ; thus at nonzero temperature , there is no contribution to @xmath195 from fig . ( 3.b ) , and @xmath206 comparing to ( [ eq:5.1a ] ) , again fig . ( 3.a ) gives the term @xmath207 , and fig . ( 3.b ) the term @xmath208 . that fig . ( 3.b ) does nt contribute is just like the same result for @xmath92 , fig . ( 1.d ) and ( [ eq:3.8 ] ) . ( 3.b ) will contribute , however , through the term @xmath209 in @xmath203 in  ( [ eq : b.5 ] ) . comparing with the tensor decomposition of  ( [ eq:4.15 ] ) , we find : @xmath210 since these are the only diagrams at one loop order , @xmath211 consequently , @xmath212 thus we can see that our results satisfy current conservation , ( [ eq:4.17 ] ) and the anomaly equation of ( [ eq:4.18 ] ) . last but not least , the latter only holds given @xmath13 and @xmath213 in ( [ eq:1.1 ] ) and ( [ eq:3.8 ] ) , and so provides a nontrivial check of these calculations . notice that while there is a factor of @xmath214 in @xmath169 , it is essentially kinematic in origin , as envisioned by itoyama and mueller @xcite . one factor of @xmath215 arises from the definition of @xmath169 , ( [ eq:4.15 ] ) , while the other can be seen to arise from a directional singularity , @xmath216 , in the integral of  ( [ eq : b.5 ] ) . further , notice that there is a ( logarithmic ) collinear divergence in the integral of  ( [ eq : b.5 ] ) . this singularity drops out of the full amplitude after contracting with the polarization tensor of the photon . in this paper we have concentrated exclusively on the anomalous decay of @xmath3 . we have done so because it is the most familiar anomalous decay , and because the connection to the axial anomaly is especially close . for the collisions of heavy ions at ultrarelativistic energies , though , this is a purely academic point , since any high temperature plasma flies apart long before a @xmath0 has a chance to decay electromagnetically . our basic point is much more general , however . while the axial anomaly for fermions is independent of temperature , anomalous decays of mesonic fields change with temperature . processes of obvious interest for hadronic systems include @xmath217 into two gluons and the decays of the @xmath68 meson @xcite . electroweak processes include couplings of the axion . whether the changes in these couplings are physically relevant can only be determined after detailed calculation ; what is certain is that they do change . in a cosmological context , the mechanisms of electroweak baryogenesis involve the effective coupling of a pseudoscalar field to the pontryagin density for the @xmath218 gauge field . in vacuum , this coupling is directly related to the anomalous divergence of the baryon current . it was argued in  @xcite that this coupling must be suppressed above the electroweak phase transition ; presumably this can be understood from our analysis . lastly , we note that the t hooft anomaly matching conditions @xcite strongly constrain the appearance of massless bound states in confining theories . due to the failure of the sutherland - veltman theorem , sec . iv , these conditions are clearly much less restrictive at nonzero temperature .    if nothing else , perhaps this gives us a greater appreciation of the wonder of the fermion axial anomaly . while every other anomalous decays changes in complicated and detailed ways , that alone remains inviolate , always . we thank w. marciano for useful conversations . this work is supported by a doe grant at brookhaven national laboratory , de - ac02 - 76ch00016 . in this appendix we review the wess - zumino - witten ( wzw ) model coupled to an external photon field . along with establishing notation , this also enables us to discuss a novel form of the wzw model , mentioned recently  @xcite , and to comment about the two currents in the two - flavor version of the wzw model .    for @xmath219 flavors , the model is constructed from a @xmath220 unitary matrix @xmath221 , @xmath222 . in the absence of gauge fields , the action is the sum of two terms , @xmath223 , with @xmath224 the usual action for a nonlinear sigma model , @xmath225 the generators of @xmath226 are the matrices @xmath227 , normalized as @xmath228 . this lagrangian is invariant under global @xmath229 unitary transformations , @xmath230 . for example , under axial transformations , for which @xmath231 , the corresponding conserved current is @xmath232 the second piece of the action is the wess - zumino - witten term , @xmath233 where the integral is over a five - dimensional region whose boundary is four - dimensional spacetime . we wish to couple @xmath221 to a photon field @xmath234 in a gauge invariant manner . for three flavors the charge matrix is @xmath235 , and so we introduce the covariant derivative , @xmath236 $ ] . by construction , both @xmath221 and @xmath237 transform covariantly under local @xmath238 gauge transformations , @xmath239 , @xmath240 , @xmath241 , @xmath242 . as usual , to make @xmath243 gauge invariant we simply replace the ordinary by the covariant derivative , @xmath244 . we can also do this for @xmath245 by replacing @xmath246 with the gauge covariant @xmath247 . using @xmath248 in @xmath245 gives something which is manifestly gauge invariant , but incomplete , since it can be shown that the ensuing equations of motion depend on the fifth dimension . in five dimensions , however , there are several other gauge invariant terms which can be added , involving powers of the ( abelian , gauge invariant ) field strength , @xmath249 . the correct action then follows uniquely by requiring that the equations of motion are independent of the fifth dimension@xcite , @xmath250 this expression is manifestly gauge invariant but not obviously independent of the fifth dimension  @xcite . this is in contrast to the usual action , in which the terms which couple to the gauge field are manifestly four dimensional , but not evidently gauge invariant : @xmath251 , where @xmath252 where @xmath253 and @xmath254 the lagrangian  [ eq : a.5 ] is the form given by witten  @xcite : it is gauge invariant up to boundary terms .    henceforth we follow brihaye , pak , and rossi  @xcite and restrict ourselves to two flavors . introducing the pion field @xmath255 , @xmath256 then for @xmath257 @xmath258 @xmath259 . the original lagrangian @xmath243 becomes @xmath260 the axial current of ( [ eq : a.2 ] ) is @xmath261 in the wzw term , @xmath262 vanishes , while using the identity , @xmath263 we find that @xmath264 the charge @xmath265 ; in going from ( [ eq : a.6 ] ) to ( [ eq : a.13 ] ) , it turns out that the contribution from the piece @xmath266 drops out . this demonstrates that @xmath267 is directly proportional to the baryon current  @xcite . further , we find that @xmath268 is proportional to the axial current in the isospin-@xmath121 direction , ( [ eq : a.11 ] ) , @xmath269 that for two flavors @xmath267 and @xmath268 are proportional to the baryon and axial isospin currents does not seem to have been recognized previously . notice that @xmath268 is only proportional to the axial current for the original lagrangian , @xmath243 ; as is discussed in sec  [ sec : adlerbardeen ] , the complete axial current includes contributions from the wzw term . we do not know if @xmath267 and @xmath268 are equal to the analogous currents for three or more flavors .    to conclude , we discuss how to compute the axial current @xmath270 of sec . [ sec : adlerbardeen ] using the noether construction . the general form for the gauged axial current is rather involved . since we compute perturbatively , however , we just compute the current in the same way . under an infinitesimal axial rotation , @xmath271 , the pion field transforms nonlinearly , @xmath272 to @xmath273 and @xmath274 . after performing such a transformation , the requisite current is then the coefficient of @xmath275 .    in the present case , however , while the lagrangian is ( of course ) gauge invariant , the lagrangian density need not be . in particular , the lagragian density of ( [ eq : a.5 ] ) is not gauge invariant ; under a local gauge transformation , it transforms by a total derivative . to the order at which we compute , we add the following term to the lagrangian density : @xmath276 @xmath277 for example , the first piece contributes a chern - simons term to the current . with the addition of ( [ eq : a.20 ] ) , the lagrangian density is manifestly gauge invariant , and the axial current is @xmath278 when @xmath279 , @xmath182 reduces to the axial vector current in the absence of electromagnetism , @xmath280 , as can be verified by expanding ( [ eq : a.11 ] ) . the terms linear and quadratic in @xmath180 agree with the calculation from the wess - zumino consistency condition , ( [ eq:5.1 ] ) . this ambiguity in the construction of the axial current is familiar . instead of the gauge invariant , anomalous current used in sec . [ sec : ward ] , by subtracting off a chern - simons term , we can choose to work instead with a current which is conserved but not gauge invariant . we collect some results of hard thermal loops , with minor differences in notation from previous work  @xcite . the simplest hard thermal loop is the integral of ( [ eq:3.1 ] ) , @xmath281 . after that , there is @xmath86 of ( [ eq:3.3 ] ) and ( [ eq:3.4 ] ) . by definition the hard thermal loop includes only only the terms @xmath77 in the integral , in the limit of soft external momentum @xmath282 . in this approximation , @xmath45 is tranverse : @xmath283 one way of writing @xmath284 is @xmath285 where @xmath286 and @xmath111 ; @xmath112 is a three vector of unit norm , @xmath287 , so that @xmath288 is null , @xmath289 , and @xmath290 . it is a dummy variable , in that one integrates over all directions of @xmath112 , as @xmath291 . in component form , @xmath292 @xmath293 @xmath294 @xmath295 where @xmath296 , and the @xmath297 are legendre functions of the second kind , @xmath298 for @xmath3 , we need the value of @xmath86 near the light cone , @xmath299 , @xmath300 @xmath301 where @xmath302 . an equivalent but useful expression for ( [ eq : b.2 ] ) is @xmath303 \ ; ; \end{aligned}\ ] ] from which follows @xmath304 @xmath305 the abelian field strengths of the vector fields @xmath306 and @xmath307 enter because @xmath78 is transverse , ( [ eq : b.1 ] ) . we found in sec .  [ sec : lowt ] that the contribution from fig . ( 1.d ) vanishes to order @xmath309 . in this appendix we check this result using the real time approach . in this instance , the real time method is simple and not problematic . it also shows that the diagram has no temperature dependent terms whatsoever : the entire diagram is equal , identically , to its value at zero temperature .    after dropping irrelevant terms proportional to @xmath310 , fig . ( 1.d ) is proportional to the integral @xmath311 @xmath312 , @xmath313 . the physical amplitude is obtained by contraction with the polarization vectors for the two photons . if @xmath314 lies in the @xmath315-direction , say , then @xmath316 must lie along the @xmath317- or @xmath318-directions . the angular @xmath319 integration will then vanish unless @xmath320 . thus the integral reduces to @xmath321 @xmath322 this integral vanishes after integration over @xmath323 ; note that this physical amplitude is free of any collinear divergences . this confirms our results obtained with the imaginary time formalism . in this section we compute the one loop corrections to the amplitude @xmath325 for soft , cool pions . this provides another , less trivial , example of hard thermal loops . corrections to the five dimensional wess - zumino - witten term , ( [ eq : a.3 ] ) , have been computed in  @xcite using a background field method ; it would be interesting using this method to compute the one loop corrections to the gauged wzw model , ( [ eq : a.5 ] ) .        at tree level , the amplitude for @xmath327 scattering is @xmath328 where @xmath329 is the momentum of the pion with isospin @xmath330 , _ etc _ , and @xmath331 from ( [ eq : a.10 ] ) , the amplitude for the scattering between four pions is @xmath332 \ ; , \nonumber\end{aligned}\ ] ] where @xmath333 , _ etc_.      in the vacuum , @xmath336 , up to terms @xmath337 which drop out of the amplitude . to @xmath114 , @xmath338 the @xmath339 comes from @xmath340 , fig . ( 1.a ) and ( [ eq:2.5a ] ) , the @xmath341 from fig . ( 4.a ) and ( [ eq : d.1 ] ) , the @xmath342 from fig . ( 4.b ) and  ( [ eq : d.5 ] ) . using ( [ eq:2.8 ] ) , then , to one loop order in vacuum  ( [ eq : d.3 ] ) renormalizes with no change in form , @xmath343 analogous to ( [ eq:2.10 ] ) and ( [ eq:2.17 ] ) .    at low temperature , @xmath344 is replaced by @xmath345 of ( [ eq:3.3 ] ) and ( [ eq:3.4 ] ) . using  ( [ eq : b.7 ] ) , we find that the effective lagrangian for @xmath346 is similar to that for @xmath3 , ( [ eq:3.9 ] ) . one term is as at zero temperature , with @xmath7 replaced by @xmath13 , while the second is a hard thermal loop : @xmath347 @xmath348 where @xmath349    s. l. adler , phys . rev . * 177 * , 2426 ( 1969 ) ; j. s. bell and r. jackiw , nuovo cim . * 60 * , 147 ( 1969 ) ; s.l . adler and w.a . bardeen , phys . rev . * 182 * , 1517 ( 1969 ) ; b. zumino , y. s. wu , and a. zee , nucl . phys . * b239 * , 477 ( 1984 ) ; m. a. shifman , phys . rept . * 209 * , 341 ( 1991 ) . l. dolan and r. jackiw , phys . rev . * d9 * , 3320 ( 1974 ) ; m. reuter and w. dittrich , phys . rev . * d32 * , 513 ( 1985 ) ; f. ruiz ruiz and r. alvarez - estrada , phys . lett . * b180 * , 153 ( 1986 ) ; a. das and a. karev , phys . rev . * d36 * , 623 ( 1987 ) ; r. baier and e. pilon , zeit . * c52 * , 339 ( 1991 ) . h.  itoyama and a.  h. mueller , nucl . phys . * b218 * , 349 ( 1983 ) . l . liu and g ni , phys . rev . * d38 * , 3840 ( 1988 ) ; x .- g . he and g. c. joshi , _ ibid . _ * 41 * , 3796 ( 1990 ) ; r .- t . wang , jour . of phys . * a23 * , 5555 ( 1990 ) ; a. a. rawlinson , d. jackson , and r. j. crewther , zeit . * c56 * , 679 ( 1992 ) ; a.  v. smilga , phys . rev . * d45 * , 1378 ( 1992 ) ; a. gomez nicola and r. f. alvarez - estrada , intl . jour . of mod a9 * , 1423 ( 1994 ) . r.  jackiw , in _ current algebra and anomalies _ , ( princeton univ . press , princeton , 1985 ) ; c. itzykson and j .- b . zuber , _ quantum field theory _ , ( mc - graw - hill , new york , 1980 ) g. m. shore and g. veneziano , nucl * b381 * , 3 ( 1992 ) . a.  dobado , r. f. alvarez - estrada , and a.  gomez nicola , phys . lett . * b319 * , 238 ( 1993 ) . a.  dobado , r. f. alvarez - estrada , and a.  gomez nicola , phys . lett . * b324 * , 345 ( 1994 ) . s. weinberg , physica * 96a * , 327 ( 1979 ) ; j. gasser and h. leutwyler , ann . of phys . * 158 * , 142 ( 1984 ) ; j. f. donoghue , e. golowich , and b. r. holstein , _ dynamics of the standard model _ ( cambridge univ . press , ny , 1992 ) ; g. ecker , prog . part . . phys . * 35 * , 1 ( 1995 ) . j.  wess and b.  zumino , phys . lett * b37 * , 95 ( 1971 ) ; e.  witten , nucl . phys . * b223 * , 422 ( 1983 ) o. kaymakcalan and j. schechter , phys . rev . * d31 * , 1109 ( 1984 ) ; y. brihaye , n. pak , and p. rossi , nucl . phys . * b254 * , 71 ( 1985 ) ; u. g. meissner , phys . * 161 * , 213 ( 1988 ) . c. m. hull and b. spence , mod . . lett . * a6 * , 969 ( 1991 ) . m. creutz and m. tytgat , phys . lett . * 76 * , 4671 ( 1996 ) ; m. creutz and m. tytgat , proceedings of lattice 96 , bnl preprint bnl-63465 , hep - lat/9608007 b. r.  holstein and j. f.  donoghue and y. c. r. lin , phys . lett . * 55 * , 2766 ( 1985 ) ; j. f.  donoghue and d. wyler , nucl . phys . * b316 * , 289 ( 1989 ) ; j. bijnens , a. bramon and f. cornet z. phys . * c46 * , 599 ( 1990 ) ; r. akhoury and a. alfakih , ann . phys . * 210 * , 81 ( 1991 ) ; j.  bijnens , nucl . phys . * b367 * , 709 ( 1991 ) ; int . j. mod * a8 * , 3045 ( 1993 ) , and references therein . r. d. pisarski and m. tytgat , bnl preprint bnl - pt-961 , hep - ph/9611206 . c. contreras and m. loewe , * c40 * , 253 ( 1988 ) ; t. hashimoto , k. hirose , t. kanki , and o. miyamura , phys . rev . * d37 * , 3331 ( 1988 ) ; b. pin - zhen and j. rafelski , mod . a7 * , 2493 ( 1992 ) ; ( e ) _ ibid . _ * a7 * , 3439 ( 1992 ) ; a. gomez nicola and r. f. alvarez - estrada , zeit . * c60 * , 711 ( 1993 ) ; s. p. klevansky , nucl . * a575 * , 605 ( 1994 ) ; d. blaschke , m. jaminon , yu . l. kalinovsky , p. petrow , s. schmidt , and b. van den bossche , _ ibid . _ * a2415 * , 1 ( 1995 ) . p. binetruy and m. k. gaillard , phys * d32 * , 931 ( 1985 ) ; j. gasser and h. leutwyler , phys . lett . * b184 * , 83 ( 1987 ) ; m. dey , v. l. eletsky , and b. l. ioffe , _ ibid . _ * b252 * , 620 ( 1990 ) ; a. bochkarev and j. kapusta , phys . rev . * d54 * , 4066 ( 1996 ) . r. d. pisarski and m. tytgat , phys . rev . * d54 * , 2989 ( 1996 ) . r. d. pisarski , phys . * 63 * , 1129 ( 1989 ) ; e. braaten and r. d. pisarski , nucl . phys . * b337 * , 569 ( 1990 ) ; _ ibid . _ * b339 * , 310 ( 1990 ) ; phys . rev . * d45 * , r1827 ( 1992 ) ; j. frenkel and j. c. taylor , nucl . phys . * b334 * , 199 ( 1990 ) ; _ ibid . _ * b374 * , 156 ( 1992 ) ; j. c. taylor and s. m. h. wong , _ ibid . _ * b346 * , 115 ( 1990 ) ; j. p. blaizot and e. iancu , _ ibid . _ * b390 * , 589 ( 1993 ) ; _ ibid . _ * b417 * , 608 ( 1994 ) ; _ ibid . _ * b434 * , 662 ( 1995 ) ; r. efraty and v. p. nair , phys . rev . * d47 * , 5601 ( 1993 ) ; r. jackiw and v. p. nair , _ ibid . _ * d48 * , 4991 ( 1993 ) ; r. jackiw , q. liu , and c. lucchesi , _ ibid . _ * d49 * , 6787 ( 1994 ) ; p. f. kelly , q. liu , c. lecchesi , and c. manuel , _ ibid . _ * d50 * , 4209 ( 1994 ) . d. g. sutherland , nucl . b2 * , 433 ( 1967 ) ; m. veltman , proc . a301 * , 107 ( 1967 ) . pisarski , in _ from thermal field theory to neural networks : a day to remember tanguy altherr _ , edited by p. aurenche , p. sorba , and g. veneziano . ( world scientific publishing , singapore , 1996 ) ; phys . rev . lett . * 76 * , 3084 ( 1996 ) ; r. baier , m. dirks , and o. kober , phys . rev . * d54 * , 2222 ( 1996 ) . l. ametller , j. bijnens , a. bramon and f. cornet , phys . * b303 * , 140 ( 1993 ) ; j. bijnens , a. bramon and f. cornet , z. phys . * c46 * , 599 ( 1990 ) . @xmath350 , ( 2.13 ) , equals @xmath351 , ( 5 ) of ref . the terms @xmath77 in @xmath350 define @xmath75 , ( 3.4 ) . from fig . ( 1.b ) of ref . @xcite , which is our fig . ( 1.d ) , in ( 12 ) they obtain @xmath352 this does not agree with ( [ eq:3.7 ] ) , for which the coefficient of the term @xmath353 vanishes . s. j. brodsky and g. p. lepage , phys . lett . * 87b * , 359 ( 1979 ) ; phys . rev . * d22 * , 2157 ( 1980 ) ; _ ibid . _ , * d24 * , 1808 ( 1981 ) ; l. arnellos , w. j. marciano , and z. parsa , nucl . * b 196 * , 371 ( 1982 ) ; a. anselm , a. johansen , e. leader , and l. lukaszuk , hep - ph/9603444 . m. dine and s. thomas , phys . lett . * b328 * , 73 ( 1994 ) . g. t hooft , in _ recent developments in gauge theories _ , ed . g. t hooft _ et al . _ ( plenum press , n.y . , 1980 ) . this method also allows one to gauge an arbitary subgroup @xmath126 of @xmath354 , whether or not @xmath126 is anomaly free . if the @xmath126 anomalies do not cancel , however , there is a doubling of the effective number of degrees of freedom , similar to the doubling which occurs for fermions on a lattice  @xcite .
Abstract:

Answer: at zero temperature , in the chiral limit the amplitude for @xmath0 to decay into two photons is directly related to the coefficient of the axial anomaly . at any nonzero temperature , this direct relationship is lost : while the coefficient of the axial anomaly is independent of temperature , in a thermal bath the anomalous ward identities do not uniquely constrain the amplitude for @xmath1 . explicit calculation shows that to lowest order about zero temperature , this amplitude decreases .


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 study of the electronic structure of materials usually begins with noninteracting electrons due to the vast number of particles involved . many - body methods then provide a hierarchy of corrections to account for the coulomb interaction to various order . many - body approaches such as gw@xcite and the bethe - salpeter equation ( bse)@xcite are frequently and successfully employed in the calculation of the electronic structure and excitations of materials . though accurate and physically sound , these many - body methods can become cumbersome and impractical for large systems due to the steep scaling of the numerical cost versus the system size .    alternatively , density - functional theory ( dft ) and time - dependent density - functional theory ( tddft)@xcite are popular methods for calculating electronic ground states and excitations , respectively , and are widely used in chemistry , physics , materials science , and other areas . density - functional methods solve the many - body problem by constructing a noninteracting system which reproduces the electronic density of the interacting , physical system . the favorable balance between accuracy and efficiency makes the resulting dft and tddft schemes unrivaled for large but finite system sizes.@xcite    considerable effort has been spent to replicate this success of tddft for periodic solids.@xcite generally speaking , tddft works very well for simple metallic systems , where the excitation spectrum is dominated by collective plasmon modes . the reason is that common local and semilocal exchange - correlation ( xc ) functionals are based on the homogeneous electron liquid as reference system , which is an ideal starting point to describe electrons in metals . the situation is more complicated in insulators and semiconductors . the first problem that comes to mind is that of the band gap , which is typically strongly underestimated by most popular xc functionals of dft . in principle , tddft provides a mechanism to obtain the correct band gap,@xcite but this puts very strong demands on the xc kernel @xmath0 ( it has to simulate a discontinuity , which requires a strong frequency dependence ) . the second difficulty are excitonic effects . it is a well - known fact that standard local and semilocal xc functionals do not produce any excitonic binding;@xcite again , the proper choice of @xmath0 is crucial . there are many examples in the literature of successful tddft calculations of excitonic effects , using exact exchange , @xcite an effective xc kernel engineered from the bse,@xcite a meta - gga kernel,@xcite and a recent ` bootstrap ' xc kernel.@xcite these kernels all have in common that they have a long spatial range ; however , it has also been shown that certain excitonic features can be equally well reproduced by simple short - range kernels . @xcite this calls for further explanation .    due to the complexity of real solids , the question of the general requirements for excitonic binding in tddft has been difficult to analyze . as a first step towards a simplified tddft approach for excitons , a two - band model was recently developed , which was used to test the performance of simple xc kernels for calculating excitonic binding energies in several iii - v and ii - vi semiconductors.@xcite in this paper we will push this reductionist approach further and propose a minimal tddft model for excitons . our model is one - dimensional ( 1d ) and uses two simple kronig - penney - type bands as input . we show that the minimal model reproduces and reveals many aspects of excitons . the model is accessible and relatively easy to implement , and it can be used to identify important aspects of the xc functional for excitonic effects . it clearly shows that excitons are collective excitations of the many - body system : the appropriate phase - coherent mixing of the single - particle excitations is accomplished via a coupling matrix featuring @xmath0 . the properties of this coupling matrix are analyzed and compared with its bse counterpart .    in textbooks , excitons are usually introduced through the two - body wannier equation , which describes an electron and a hole which interact via a screened coulomb potential . while this arguably constitutes the simplest model for excitons , it is based on several drastic assumptions which are not fulfilled in general.@xcite we discuss how and under what circumstances a wannier - like equation emerges from our minimal tddft model . the paper is structured as follows . we give an introduction to wannier excitons in sec . [ sect : bkgnd ] and to tddft for solids in sec . [ sect : tddft ] , followed by a description of the minimal model in sec . [ sect : model ] . section [ sect : results ] then presents results for the minimal exciton model comparing tddft with the bse , and discusses various implications . in sec . [ sect : wannier ] we show how the wannier equation emerges from tddft . conclusions are given in in sec . [ sect : app : details ] . details of the bse method are provided in appendix [ sect : app : bse ] . atomic units ( @xmath1 ) are used throughout unless otherwise stated , and we will only consider spin - unpolarized systems . the electronic structure of crystalline solids is described by the bloch theory , where the electrons move in a periodic effective single - particle potential which reflects the crystal symmetry . as a result , the electronic states form energy bands . in insulators and semiconductors , electronic excitations take place between the occupied ( valence ) and unoccupied ( conduction ) bands . these interband transitions can be described within a simple independent - particle approach based on fermi s golden rule ; one thus obtains a reasonable qualitative account of the optical properties in these materials.@xcite        however , experiments reveal that there are important modifications to this picture , as illustrated schematically in fig . [ fig : bkgnd : illustrationexcitons ] . above the band gap , the spectrum appears strongly enhanced , and below the band gap one may find discrete absorption peaks known as bound excitons . the origin of these modifications are coulomb interactions : the simple picture of independent single - particle excitations is replaced by a more complex scenario where these excitations are dynamically coupled . excitonic effects are ubiquitous in nature , and occur in 3d , 2d and 1d systems alike.@xcite the details of the excitonic modifications to the noninteracting spectrum depend on the dimensionality of the system.@xcite        the standard textbook explanation of excitons is based on the simple picture of an electron - hole pair held together by coulomb interactions , see fig . [ fig : bkgnd : wannierexciton ] . one thus arrives at a two - body problem similar to the positronium atom , with a center - of - mass momentum @xmath2 and relative motion described by the wannier equation:@xcite @xmath3\psi_\nu({{\bf r}})=e_\nu\psi_\nu({{\bf r } } ) , \label{eqn : bkgnd : wannier}\ ] ] where @xmath4 is the reduced mass , defined as @xmath5 . @xmath6 and @xmath7 are the effective masses of the conduction band electrons and the valence band holes . @xmath8 is the coulomb interaction between the electron and the hole , divided by the static dielectric constant of the system . @xmath9 and @xmath10 are the excitonic wave function and binding energy , respectively .    in 1d systems , the coulomb interaction is ill - defined and requires , in general , some parametrized form;@xcite we will use the following soft - coulomb interaction : @xmath11 where @xmath12 and @xmath13 are parameters . in the following , we will use @xmath14 throughout . the wannier equation ( [ eqn : bkgnd : wannier ] ) has the form of a hydrogenic schrdinger equation , so it possesses a rydberg series ( even for 1d cases with soft - coulomb interaction ) with infinitely many eigenvalues below the band gap , and a continuum above the gap.@xcite the excitons in the wannier model for 3d and 2d cases both enhance the optical absorption near the band gap , while the presence of excitons in 1d systems suppresses the optical absorption just above the band gap . sham and rice@xcite showed how the wannier equation can be derived from first principles starting from the bse , under the assumption that the bohr radius of the exciton is much bigger than the lattice constant . the resulting wannier picture of excitons appears clear and intuitive , but this simplicity is somewhat deceptive . in reality , excitons are a dynamical many - body phenomenon and require a subtle coordination and cooperation of many single - particle transitions between two bands . in the following , we will develop a model based on tddft which will illustrate the true physical nature of excitons , but which will remain sufficiently transparent to allow a simple interpretation of the collective many - body effects that are responsible for the excitonic binding . tddft@xcite is an in principle exact approach for electron dynamics , based on the uniqueness of the mapping between the time - dependent electronic density and the external potential.@xcite the key equation of tddft is the time - dependent kohn - sham ( tdks ) equation : @xmath15\phi_i({{\bf r}},t ) , \label{eqn : bkgnd : tddft}\end{aligned}\ ] ] where @xmath16 are the tdks orbitals of the noninteracting kohn - sham system which reproduces the density of the real interacting system . @xmath17 and @xmath18 are the external potential of the physical system and the hartree potential , respectively , and the xc potential @xmath19 is the only piece that needs to be approximated in practice .    the excitation spectrum of a system can be calculated via time propagation of eq . ( [ eqn : bkgnd : tddft ] ) following a suitably chosen initial perturbation . alternatively , one can obtain excitation energies and optical spectra directly from linear - response tddft , using the so - called casida equation:@xcite @xmath20 equation ( [ eqn : bkgnd : casida ] ) is a generalized eigenvalue equation . one obtains the optical transition frequencies from the eigenvalues @xmath21 , and the corresponding eigenvectors tell us how the kohn - sham single - particle transitions are mixed to form the transitions of the interacting system . the matrices @xmath22 and @xmath23 in eq . ( [ eqn : bkgnd : casida ] ) are defined as @xmath24 where @xmath25,@xmath26,@xmath27,@xmath28 are labels for ground - state kohn - sham orbitals , and the @xmath29 s are the associated kohn - sham orbital energies . @xmath30 in eq . ( [ eqn : bkgnd : ab ] ) is defined as @xmath31 where the @xmath32 s are the ground - state kohn - sham orbitals , and @xmath33 is the hartree - exchange - correlation ( hxc ) kernel , defined as a fourier transform of @xmath34 the xc kernel @xmath35 has to be approximated in practice . @xmath36 and @xmath37 make up the eigenvector in eq . ( [ eqn : bkgnd : casida ] ) , and they describe excitations and de - excitations , respectively . a commonly used approximation to eq . ( [ eqn : bkgnd : casida ] ) , known as tamm - dancoff approximation ( tda),@xcite is to set @xmath38 , so that the casida equation reduces to @xmath39\rho_{mn}(\omega)=\omega\rho_{ij}(\omega ) . \label{eqn : bkgnd : tddftworking}\ ] ] this decouples excitations and de - excitations , and the computational cost is reduced . there are situations , for instance for molecular excitations of open - shell systems,@xcite in which the tda is preferred in practice over the full casida equation ( [ eqn : bkgnd : casida ] ) . we find that the tda can also be advantageous for excitons ( see sec . [ sect : results ] ) .    by considering only a single kohn - sham transition in eq . ( [ eqn : bkgnd : casida ] ) one arrives at the small - matrix approximation ( sma):@xcite @xmath40 where @xmath41 is the kohn - sham transition frequency to be corrected . one can further simplify this by making the tda , which yields the single - pole approximation ( spa ) : @xmath42 the sma and spa are valid when the considered excitation is far away from other transitions in the system . though they are not usually accurate enough for real calculations , their simplicity makes them very useful for theoretical analysis and development .    in periodic solids , the kohn - sham orbitals are labeled with the band index @xmath25 and the wavevector @xmath43 and have the bloch form : @xmath44 it is in principle possible to adapt the casida equation ( [ eqn : bkgnd : casida ] ) for the case of orbitals of the form ( [ eqn : bkgnd : blochform]),@xcite and use this to calculate the excitation spectrum . however , to obtain the optical spectrum in solids it is more convenient to calculate the macroscopic dielectric function:@xcite @xmath45 where the @xmath46 s are reciprocal lattice vectors , and @xmath47 is the linear response function . @xmath48 can be expressed as@xcite @xmath49 where @xmath50 labels the solutions of eq . ( [ eqn : bkgnd : tddftworking ] ) for an extended system .    beyond linear response , few real - time tddft calculations exist for periodic solids.@xcite instead of directly solving the tdks equation ( [ eqn : bkgnd : tddft ] ) , the tdks orbitals can be expanded in terms of ground - state kohn - sham bloch functions as @xmath51 the time - dependent density matrix is defined as @xmath52 the equation of motion of the density matrix is then @xmath53 , \label{eqn : bkgnd : tddmeom}\ ] ] with the tdks hamiltonian matrix @xmath54 defined by @xmath55 where @xmath56 is the volume of the unit cell , and @xmath57 is the tdks hamiltonian of eq . ( [ eqn : bkgnd : tddft ] ) . this density - matrix approach has been used to derive the tddft version of the semiconductor bloch equations.@xcite in this formalism we only consider vertical transitions , where the bloch wavevector @xmath2 does not change during the dynamics . nonvertical excitations are not considered , since they involve indirect ( i.e. , phonon - assisted ) processes which we ignore here . solids are formally described by the many - body schrdinger equation . since exact solutions are not possible , it is instructive to resort to model systems to eliminate undesired details of the many - body system and provide clear illustrations of the specific features one is interested in . the wannier equation for excitons presented in sec . [ sect : bkgnd ] is such a model .    while intuitive , the wannier model assumes the electron - hole interaction as given , so the excitonic effects are already built in by default however , this does not explain under what conditions one expects to see the formation of excitons , and the many - body nature of excitonic effects remains hidden . we therefore propose a minimal model for excitons which lowers the abstraction level of the wannier model , and where excitonic effects show up without any _ ad hoc _ assumptions .    for excitations near the band gap , a reasonable approximation is to use a two - band model , i.e. , only to consider the highest valence band @xmath58 and the lowest conduction band @xmath59 . this means that we only need to consider those elements of the time - dependent density matrix @xmath60 [ eq . ( [ tddm ] ) ] for which @xmath61 and @xmath62 ; the latter index will be dropped in the following .    for the case when a small perturbative electric field is applied to the system , it is sufficient , to lowest order in the perturbation , to consider only the time evolution of the off - diagonal part of the density matrix,@xcite @xmath63 . one then obtains from eq . ( [ eqn : bkgnd : tddmeom ] ) @xmath64 where @xmath65 , and @xmath66 . here , @xmath67 denotes the matrix elements of @xmath68 , defined similarly to eq . ( [ hmat ] ) . there is no external perturbation in eq . ( [ eqn : bkgnd : rhocvintime ] ) ; we assume free propagation of the system . our reasoning is that individual excitations can be viewed as the eigenmodes of the system and do not depend on the specific form of the perturbative field . fourier transformation of eq . ( [ eqn : bkgnd : rhocvintime ] ) gives @xmath69 where the factor 2 accounts for the spin , and @xmath70    equation ( [ eqn : bkgnd : rhocv ] ) is equivalent to the sma for finite systems . while the sma only refers to the transition between one individual occupied and one individual unoccupied orbital , eq . ( [ eqn : bkgnd : rhocv ] ) considers the transitions between the valence and the conduction band as a whole . ignoring the coupling between excitations and de - excitations by setting @xmath71 ( i.e. , making the tda ) , one arrives at the solid - state analog of the spa : @xmath72\rho_{{\mathbf{k}}'}^{cv}(\omega ) = \omega\rho_{\mathbf{k}}^{cv}(\omega ) . \label{eqn : bkgnd : rhocvspa}\ ] ] this is the central equation which we will use to describe excitonic effects . it requires as input the kohn - sham bloch functions for the valence and the conduction band of an insulator or semiconductor . to keep things as simple as possible , we will consider the kronig - penney ( kp ) model@xcite rather than a real material . the kp model is a 1d noninteracting system with a periodic potential of square wells . within the unit cell @xmath73 $ ] , the potential is @xmath74        a typical example for the band structure of the kp model is plotted in fig . [ fig : bkgnd : kppot ] . despite its simple appearance , the kp model is very versatile : by varying the values of the lattice constant @xmath75 , barrier width @xmath76 , and barrier height @xmath77 , a wide range of band gaps and band curvatures can be achieved . a square - well potential of finite depth does not support an infinite number of bound states as the coulomb potential does ; but this , in fact , closely reflects the reality of the effective potential felt by the valence electrons in a solid , which is relatively shallow due to the screening of the bare nuclear charges by the core electrons . in practice , many solid - state calculations account for this screening by using muffin - tin potentials or pseudopotentials . @xcite the kp model can be viewed as an elementary version of this approach .    in the following , we always choose the first two bands to be fully occupied and the higher bands to be empty . we then make the two - band approximation for band 2 and 3 , since the shapes of these bands resemble the highest valence band and lowest conduction band in direct - gap materials such as gaas . the bands in the kp model are sufficiently well separated , so the two - band approximation is justified .    to establish a connection with tddft , we assume the solution to the noninteracting kp model as our ground - state kohn - sham system . in other words , the potential in eq . ( [ eqn : bkgnd : kppot ] ) represents the _ exact _ kohn - sham potential @xmath78 which corresponds to a physical system whose external potential @xmath79 is uniquely determined thanks to the hohenberg - kohn theorem of dft . for our purpose , it is not necessary to know what this external potential looks like . the kohn - sham bloch functions can then be determined in an elementary fashion.@xcite    for comparison , we will also carry out bse calculations in our model system ( see appendix [ sect : app : bse ] for technical details ) . bse calculations are typically based on ground - state quasiparticle states obtained from the gw method.@xcite this is because the single - particle gap in gw is usually closer to experiment than the approximate kohn - sham gap . however , in our case this distinction is not important because we use the given kp band structure as input for both bse and tddft .    in summary , our minimal tddft model for excitons consists of the following two ingredients :    \(1 ) a two - band model for the vertical transitions between the highest valence band and the lowest conduction band , see eq . ( [ eqn : bkgnd : rhocvspa ] ) ;    \(2 ) the band structure from a 1d kp model . of course , the model is not complete without a choice for the xc kernel @xmath35 . this will be discussed below . the exact xc kernel @xmath80 is unknown and must be approximated ; we restrict ourselves to adiabatic kernels that have no frequency dependence . the adiabatic local - density approximation ( alda ) , as well as all semilocal , gradient - corrected xc kernels , are known to be unable to describe excitonic effects.@xcite the exact xc kernel has a long - range decay of @xmath81 , which is absent in all ( semi)local xc kernels derived from the uniform electron gas . this long - range part is thought to be essential for excitons.@xcite    the long - range behavior of @xmath35 depends on the dimensionality , and in our 1d model system we define the following long - ranged ( or ` soft - coulomb ' ) xc kernel : @xmath82 we also consider an extremely short - ranged contact xc kernel : @xmath83 these model xc kernels depend on the constants @xmath84 and @xmath85 , which we will treat as fitting parameters in the following . the idea is to tune the parameters in the model xc kernels so that bound excitons are produced , and to align the lowest bound exciton in the tddft spectrum with the lowest bound exciton in the bse spectrum . results for the imaginary part of the dielectric function are presented in figs . [ fig : results : tddftoneexciton ] and [ fig : results : tddfttwoexcitons ] . we find that both the long - ranged @xmath86 and the short - ranged @xmath87 produce bound excitons , and thus , strictly speaking , the long - range behavior of the xc kernel is not really required for excitonic effects . the bse results in figs . [ fig : results : tddftoneexciton ] and [ fig : results : tddfttwoexcitons ] show several identifiable bound excitons , in agreement with the rydberg series predicted by the wannier model ; the number of visible bound excitons somewhat depends on the numerical resolution in momentum space .    for the kp model parameters of fig . [ fig : results : tddftoneexciton ] , we find that the adiabatic tddft can only bind a single excitonic state . for other kp parameters ( specifically those in which the lowest conduction band is above the barrier ) , tddft produces two excitons , see fig . [ fig : results : tddfttwoexcitons ] , which agree well with the lowest two excitons in bse . there are additional , higher - lying bound excitons in bse which are very faint and difficult to resolve numerically . for all the kp systems we tested , we never found more than two bound excitons with tddft . this indicates the limitations of the adiabatic xc kernels used here .     +    as mentioned in sec . [ sect : model ] , the wannier model does not clearly demonstrate that excitons are collective excitations . since the wannier model assumes a single electron - hole pair picture , one can not immediately see that excitonic excitations are composed of a coherent superposition of many single - particle excitations . in our minimal model , we solve the eigenvalue equation ( [ eqn : bkgnd : rhocvspa ] ) , and the eigenvectors @xmath88 ( which depends on @xmath21 parametrically ) describe how the transitions between noninteracting orbitals form the transitions in the interacting system . @xmath89 is the percentage of a noninteracting transition in the transition of the interacting system . two typical cases are plotted in figs . [ fig : results : excitoniceigenvect ] and [ fig : results : nonexcitoniceigenvect ] . figure [ fig : results : excitoniceigenvect ] clearly shows that excitons are collective excitations which are formed by mixing a wide distribution of single - particle transitions . as expected , the lowest exciton eigenfunction is nodeless and the second excitonic eigenfunction has a single node . with purely parabolic bands , the results from the wannier model would be recovered , as we will show below . in contrast , the transitions in the continuum shown in fig . [ fig : results : nonexcitoniceigenvect ] have a strong single - particle character ( the two peaks arise from the @xmath90 degeneracy in the kp model ) .    equation ( [ eqn : bkgnd : rhocvspa ] ) is equivalent to the spa for finite systems , which ignores the coupling between excitations and de - excitations ( tda ) . we also investigated what happens when we do not make the tda , i.e. , when we work instead of eq . ( [ eqn : bkgnd : rhocvspa ] ) with the full equation for the two - band model , eq . ( [ eqn : bkgnd : rhocv ] ) . as long as we describe relatively weakly bound excitons that are not too far below the band gap , we find that the difference between the two methods is very minor . however , we also discovered that , under rare circumstances , eq . ( [ eqn : bkgnd : rhocv ] ) can lead to tddft excitonic binding energies that are purely imaginary . in our minimal model , such instabilities arise when the interaction strength @xmath12 in eq . ( [ eqn : results : fxcsc ] ) and ( [ eqn : results : fxccont ] ) increases so that the excitonic binding energy becomes greater than the band gap . this situation is comparable to the well - known triplet instability in tddft , for which the tda generally leads to an overall better behavior;@xcite for excitons binding energies in our minimal model , we draw similar conclusions . the bse scheme is commonly implemented within an adiabatic scheme ( see appendix [ sect : app : bse ] ) ; as we have seen , it produces a series of bound excitons . since we assumed that the kp model is the kohn - sham ground state in tddft and the gw quasiparticle ground state in bse , the difference between tddft and bse becomes easily comparable , since the central equation to be solved have the same form , eq . ( [ eqn : bkgnd : rhocvspa ] ) . the @xmath91 coupling matrices for tddft and bse are @xmath92 where @xmath93 is the hartree kernel ( the 1d soft - coulomb interaction ) , and @xmath94 is the screened interaction . aside from the change from @xmath35 to @xmath94 , the most prominent difference between bse and tddft is the order of the indices for @xmath94 in eq . ( [ eqn : results : fbse ] ) . since the noninteracting ground - state wave functions have the bloch form ( [ eqn : bkgnd : blochform ] ) , we can see from eq . ( [ eqn : bkgnd : fijmn ] ) that @xmath95 has a strong @xmath96 dependence ; in fig . [ fig : results : bsecoupling ] this shows up as a dominance along the diagonal . by contrast , this @xmath96 dependence is clearly absent in @xmath97 , as demonstrated in fig . [ fig : results : tddftcoupling ] . @xmath98 and @xmath99 with only vertical transitions can be expressed in momentum space as @xmath100 @xmath101 the xc matrix ( [ eqn : results : fxctrans ] ) only depends on the long - range ( @xmath102 ) behavior of its momentum space representation @xmath103 , while the @xmath94 matrix ( [ eqn : results : wtrans ] ) also depends on other @xmath104 values in its momentum space representation @xmath105 . it is impossible to find an adiabatic @xmath35 that reproduces the bse coupling matrix as in fig . [ fig : results : bsecoupling ] , since @xmath105 has an extra degree of freedom over @xmath106 . one can only hope to reproduce a portion of the bse coupling matrix with adiabatic tddft ( as pointed out in ref . ) , or make the xc kernel frequency dependent so that the information from the @xmath104-dependence in @xmath105 is mapped into the frequency dependence in @xmath107 .    considering the nature of the objects involved in this mapping , a highly nontrivial frequency dependence in @xmath35 is required to reproduce a series of bound excitons . for example , one can easily construct an @xmath35 which reproduces a given series of bound excitons by using a different contact kernel in the region @xmath108 surrounding each exciton at frequency @xmath109 : @xmath110\theta[(\omega_i+\eta_i^+)-\omega].\ ] ] here , the @xmath111 s are parameters which are adjusted so that a tddft calculation with the frequency - independent kernel @xmath112 would produce @xmath109 as the lowest excitonic binding energy . such an @xmath35 is of course completely ad hoc , but the fact that the excitonic series can be reproduced in this way demonstrates that the inclusion of the frequency dependence would greatly improve the flexibility of the tddft scheme .    on the other hand , within the adiabatic approximation the characteristics of the @xmath113 coupling matrix are important for excitonic effects . to emphasize this , we now show that in very special cases the number of discrete excitonic eigenvalues can be derived . consider a real matrix @xmath114 , where @xmath115 is a diagonal real matrix with @xmath116 , and @xmath117 has the symmetry @xmath118 . within second order perturbation theory , there is at most one discrete eigenvalue of @xmath114 in the limit where @xmath119 become continuous.@xcite though this case does not correspond to the matrices that would occur in real calculations , it indicates the close relationship between the properties of the coupling matrix and excitonic effects . it is also possible to derive properties of the discrete eigenvalues if @xmath120 and @xmath121 are completely decoupled in the xc kernel . owing to the symmetry @xmath122 implied by eq . ( [ eqn : bkgnd : fijmn ] ) , such separable kernels can only have the form @xmath123 for an excitation below the band gap with frequency @xmath21 , we can show@xcite that it must satisfy @xmath124 where the sum is carried out over the first brillouin zone ( fbz ) . equation ( [ eqn : results : eqnsepkernel ] ) shows that eq . ( [ eqn : results : sepkernel ] ) must have the negative sign in order to have bound excitons . the left - hand side of eq . ( [ eqn : results : eqnsepkernel ] ) is monotonically increasing with @xmath21 , so for separable kernels of the form of eq . ( [ eqn : results : sepkernel ] ) , there is only one bound excitonic solution .    as shown in fig . [ fig : results : tddftcoupling ] , tddft coupling matrices lack the strong dependence of @xmath125 as in bse coupling matrices . expanding the tddft coupling matrices into a power series of separable matrices and truncating at the first order would be a reasonable approximation , explaining why tddft produces fewer bound excitons ( if any at all ) than many - body methods such as bse . the contact xc kernel and the soft - coulomb xc kernel in eqs . ( [ eqn : results : fxcsc ] ) and ( [ eqn : results : fxccont ] ) have the following simple form in momentum space : @xmath126 where @xmath127 , @xmath128 and @xmath128 are reciprocal lattice vectors , and @xmath129 is the modified bessel function of the second kind . it is customary to refer to the matrix elements where @xmath130 as ` head ' , @xmath131 or @xmath132 as ` wings ' , and @xmath133 as ` body ' . the 3d coulomb potential has the form @xmath134 in momentum space . however , in 1d systems there is no real coulomb interaction which behaves as @xmath135 for @xmath136 , and one has to use the soft - coulomb interaction instead . though there are many flavors of the soft - coulomb interaction , they all have the same @xmath137 behavior for @xmath136 . however , the linear response function @xmath138 always behaves as @xmath135 for @xmath136 and does not depend on the dimensionality . this renders quantities like the macroscopic dielectric function ( [ eqn : bkgnd : macrodielec ] ) ill - defined for strictly 1d systems . furthermore , the bootstrap xc kernel@xcite and other xc kernels that depend on the cancellation of the 3d coulomb @xmath135 singularity will not work as designed in strictly 1d and 2d systems . therefore @xmath139 shown in fig . [ fig : results : tddftoneexciton ] and [ fig : results : tddfttwoexcitons ] are calculated at a small but finite @xmath140 . the coupling matrix @xmath141 can be written in momentum space as @xmath142 \nonumber\\ & & \times{\left\langle c , k \left| e^{igx } \right| v , k \right\rangle}{\left\langle v , k ' \left| e^{-ig'x } \right| c , k ' \right\rangle}. \label{eqn : app : detail : fvcvc}\end{aligned}\ ] ]    for any xc kernel that behaves as @xmath135 for @xmath136 , one can further simplify the calculation by ignoring the so - called local field effects,@xcite i.e. instead of summing over @xmath128 and @xmath143 in eq . ( [ eqn : app : detail : fvcvc ] ) , only the head is considered . in 3d systems , a prominent example is the long - range kernel @xmath144 , which is obtained as an effective xc kernel with only head matrix elements from inverting the bse of contact excitons.@xcite    on the other hand , any xc kernel that diverges more slowly than @xmath135 for @xmath136 changes the spectrum only through the local field effects , i.e. all @xmath128 and @xmath143 must be summed in eq . ( [ eqn : app : detail : fvcvc ] ) . in other words , effective xc kernels with only the head are not feasible in strictly 1d systems due to the asymptotic behavior of the soft - coulomb potential discussed above . for 1d systems with @xmath131 , we have @xmath145 and @xmath35 s with 1d long - range behavior such as the soft - coulomb kernel also behave as @xmath146 . considering eq . ( [ eqn : app : detail : fvcvc ] ) , these asymptotic properties imply that the head and wing contributions to @xmath91 always vanish in strictly 1d systems for physically meaningful xc kernels . due to these dimensionality restrictions , the xc kernel changes the strictly 1d system only through the local field effects .    in 3d the head contribution to the coupling matrix @xmath147 is much more important than the local field effects , which is the reason that long - range kernels ( with nonzero head ) work much better than local xc kernels ( with vanishing head ) such as alda . in our strictly 1d model system , the head contribution is zero even for the bse , and thus the long - range kernel does not outperform local kernels such as the contact kernel . these peculiarities only occur when one considers strictly 1d and 2d systems . in a more realistic picture , one encounters quasi-2d systems@xcite and quasi-1d systems ( such as quantum wires with finite radius or nanotubes@xcite ) , in which the movement of electrons is confined in certain directions such that the transverse motion can be averaged in comparison with the longitudinal motion . though these systems show low - dimensional characteristics in various properties due to confinement , in the limit of @xmath136 they eventually differ from strictly low - dimensional systems . our minimal model and the wannier exciton picture can be connected by considering the fourier transform of eq . ( [ eqn : bkgnd : rhocvspa ] ) . we define an effective two - body potential @xmath148 via the fourier transform of @xmath149 : @xmath150 where @xmath151 is a direct lattice vector . the fourier transform of the density matrix is @xmath152 since wannier excitons extends over many lattice constants , we approximate @xmath151 as a continuous variable @xmath153 . assuming the effective mass approximation , eq . ( [ eqn : bkgnd : rhocvspa ] ) becomes @xmath154 where @xmath155 is the excitonic binding energy , and the integration is carried out over all space . we call eq . ( [ eqn : bkgnd : tddftwannier ] ) the tddft wannier equation , since it has the same form as eq . ( [ eqn : bkgnd : wannier ] ) . with proper choice of the approximated xc kernel , the nonlocal effective electron - hole interaction potential @xmath148 supports bound excitonic states .    since the bse and the tddft are formally similar within the minimal model , eq . ( [ eqn : bkgnd : rhocvspa ] ) can also be applied to the bse results . [ fig : results : veh ] shows the effective interaction potential @xmath148 for tddft and bse . \(a ) + ( b )    the tddft wannier equation provides an intuitive way of describing the effective nonlocal electron - hole interaction , and of explaining why adiabatic tddft usually has fewer excitons than bse and the wannier model . however , in most cases the tddft wannier equation is not suitable for quantitative use due to the approximations involved . the approximation where we take the lattice vector @xmath151 as a continuous variable assumes that the exciton radius is much larger than the lattice constant ; this works fine in most cases we tested . but the effective mass approximation where @xmath156 is approximated by @xmath157 is only good for transitions near the band gap , thus requiring these transitions of the noninteracting system to dominate the exciton , which is equivalent to the exciton extending over many lattice constants . one obtains the @xmath158 term in eq . ( [ eqn : bkgnd : tddftwannier ] ) from @xmath157 in the limit where the lattice constant @xmath159 , and this approximation is not valid for most systems . although @xmath148 is a nonlocal potential , in most cases we find that the @xmath148 s for both tddft and bse are dominated by the diagonal part , so the exciton problem is in analogy to one - body systems . [ fig : results : vehdiag ] shows the diagonal part of @xmath148 , which can be taken as the effective one - body potential . the wannier model in 1d has the soft - coulomb potential , which supports an infinite number of bound excitons ( the soft - coulomb interaction is fitted so that the binding energy of the first exciton matches that of the bse ) . we find in general that the diagonal parts of @xmath148 for both bse and tddft are much more shallow than the soft - coulomb potential , and the tddft one is more narrow than the bse one . thus , bse and tddft are not able to produce a complete excitonic rydberg series , and tddft in general produces fewer bound excitons than the bse . the tddft wannier equation is not suitable for quantitative use for most of our model systems , despite the success of the wannier model in describing real semiconductors.@xcite since the approximations involved in eq . ( [ eqn : bkgnd : tddftwannier ] ) require that the exciton radius is large compared to the lattice constant , this suggests that this discrepancy is due to the special nature of 1d systems : namely , for similar effective masses the exciton radius in 1d is much smaller than in 3d and 2d.@xcite the purpose of this paper was to construct a transparent and accessible minimalist model system that produces excitonic effects in a non-_ad hoc _ fashion using tddft . the model , as presented here , is not intended to be a testing ground for xc kernels . thus despite the dimensionality restrictions for the strictly 1d system , our results in sec . [ sect : results ] carry over to 3d systems .    with our minimal model , we show that adiabatic tddft is capable of producing bound excitons through the local field effect even when the xc kernel is local in space , provided the strength of the kernel is strong enough . this statement is still true in 3d ; however , due to the non - vanishing head contribution of the exact @xmath35 , we expect that the deviation of the effective interaction strength of a local @xmath35 from the real , nonlocal @xmath35 becomes larger than the in our strictly 1d model . in this sense the long - range kernel , though very favorable , is not a necessary condition for excitonic effects . we show the connection between tddft and the wannier model for excitons by deriving the tddft wannier equation , which describes a real - space system featuring a nonlocal effective electron - hole interaction . such a connection intuitively demonstrates how adiabatic tddft generally produces fewer bound excitons than bse , and does not have a complete rydberg series . the eigenvectors of the excitonic excitations in the minimal model clearly show their collective nature , which is not obvious from the wannier model alone . excitonic instabilities may show up in tddft with approximate xc kernels , and this suggests that the tda tends to be more reliable for excitons than the formally exact method .    the frequency dependence of the exact xc kernel , @xmath160 , is usually ignored . despite the fact that adiabatic xc kernels have met with considerable recent success in producing optical spectra of insulators and semiconductors ( see the discussion in the introduction ) , they are incapable of producing excitonic rydberg series . our model system gives an explanation for why this is the case . this failure of the adiabatic approximation for @xmath35 is quite different from that which is responsible for the inability of adiabatic tddft to produce double excitations in finite systems or certain classes of charge - transfer excitations.@xcite this calls for continuing efforts in the search for nonadiabatic xc kernels for excitons . this work was supported by nsf grant no . dmr-1005651 . electrons and holes near the fermi surface are well described in the quasiparticle picture . the quasiparticle green s function @xmath128 is related to that of the noninteracting system , @xmath161 , through the use of the self - energy @xmath162 : @xmath163 where the arguments denote sets of space and time variables . a widely used approximation for the self - energy is the gw approximation : @xcite @xmath164 where @xmath94 is the screened interaction , @xmath165 @xmath166 is the bare coulomb interaction , and the inverse dielectric function @xmath167 is obtained as @xmath168 the linear response function @xmath138 can be calculated by its lehmann representation : @xmath169 where @xmath170 are quasiparticle states , @xmath171 are quasiparticle energies , and @xmath172 are occupation numbers . in practice , evaluating @xmath138 through eq . ( [ eqn : app : bse : chi ] ) can be quite time - consuming , and @xmath138 is often calculated with the plasmon - pole model.@xcite however , our minimal model is simple enough to allow us to use eq . ( [ eqn : app : bse : chi ] ) directly . the gw quasiparticle green s function obtained from eq . ( [ eqn : app : bse : gwsigma ] ) misses important dynamical many - body effects , such as the electron - hole ( excitonic ) interaction . the two - particle green s function includes these effects . the bse@xcite describes the relation between the four - point polarization function @xmath173 of an interacting system and the corresponding object of the quasiparticle system : @xmath174 in which @xmath175 is @xmath176 and assuming the gw approximation , the kernel @xmath129 is @xmath177 here , @xmath178 denotes the coulomb interaction with the long - range part removed.@xcite in practice the bse is often solved in the transition space , which is spanned by single - particle excitations . a four - point function such as @xmath179 then becomes @xmath180 where the @xmath32 s can be any complete basis set . ( [ eqn : app : bse : l ] ) in the transition space becomes @xmath181^{-1}_{(ij)(mn)}(f_m - f_n),\ ] ] where the excitonic hamiltonian matrix is @xmath182 we make the adiabatic approximation for @xmath183 and arrive at the following eigenvalue problem : @xmath184 and @xmath185 can be expressed in terms of these eigenvectors by @xmath186 only the transitions between the valence and conduction bands contribute . within our two - band model , the excitonic hamiltonian has the following block matrix form : @xmath187 ignoring the off - diagonal part in eq . ( [ eqn : app : bse : hexc ] ) is equivalent to the tda .    as shown in sec . [ sect : results ] , it is possible that instabilities show up in the full bse results when the underlying ground - state calculation is not exact . such instabilities in the minimal model are an artifact originating from the assumption that the solution of the kp model constitutes the ground state of the many - body system . however , this is not a matter of great concern in practice .    in principle , the transition space spans all possible combinations of valence and conduction orbitals , including nonvertical transitions connecting different bloch wavevectors . the kernel @xmath188 of the bse in momentum space , eq . ( [ bsekernel ] ) , has the following ingredients : @xmath189 and @xmath190 only the excitations with the same momentum transfer @xmath191 are coupled due to the @xmath192 functions in eq . ( [ eqn : app : bse : v ] ) and ( [ eqn : app : bse : w ] ) , so we only need to include vertical transitions in the calculations for optical properties .
Abstract:

Answer: the accurate description of the optical spectra of insulators and semiconductors remains an important challenge for time - dependent density - functional theory ( tddft ) . evidence has been given in the literature that tddft can produce bound as well as continuum excitons for specific systems , but there are still many unresolved basic questions concerning the role of dynamical exchange and correlation ( xc ) . in particular , the role of the long spatial range and the frequency dependence of the xc kernel @xmath0 for excitonic binding are still not very well explored . we present a minimal model for excitons in tddft , consisting of two bands from a one - dimensional kronig - penney model and simple approximate xc kernels , which allows us to address these questions in a transparent manner . depending on the system , it is found that adiabatic xc kernels can produce a single bound exciton , and sometimes two bound excitons , where the long spatial range of @xmath0 is not a necessary condition . it is shown how the wannier model , featuring an effective electron - hole interaction , emerges from tddft . the collective , many - body nature of excitons is explicitly demonstrated .