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 preparation of many - particle entangled states is becoming a fundamental task in modern quantum physics . entanglement lies at the heart of quantum communication and quantum information and is also a fundamental resource in precision spectroscopy . with recent advances in the manipulation of trapped ions and neutral atoms , there has been significant progress in the preparation of many - body entanglement states in atomic gases @xcite . such systems therefore offer rich opportunities to investigate new physics with both practical and fundamental applications .    here we propose an interferometric method that uses the entanglement that emerges during the many - body dynamics between interacting spinor atoms to perform precision measurements of their interatomic interaction strength , usually parameterized by the @xmath3-wave scattering length . our method is a many - body generalization of the technique reported in ref . @xcite , where information of the scattering properties was obtained by using the entanglement dynamics in pairs of atoms trapped in the ground state of a potential well . a precise determination of the scattering properties has broad and important applications beyond atomic physics . for example , ref . @xcite shows that such type of measurements might be useful for fundamental physics , as by monitoring the scattering length on the @xmath4 level one could detect variation in the electron to proton mass ratio on the level of @xmath5 . standard ramsey spectroscopy starts with a system of @xmath0 non - interacting spin @xmath1 particles initially prepared in the same internal state . subsequently , a ramsey pulse prepares it in a coherent superposition of the two internal states and the system is let to evolve freely for a time during which a relative phase accumulates , @xmath6 , due to the energy splitting between the two states . the latter is then decoded by using a second ramsey pulse which maps the phase onto a population difference , thus allowing to perform precision measurements of the atomic transition frequency . the statistical fluctuations associated with a finite sampling yields a lower limit in the phase accuracy @xmath7 called the shot noise limit @xcite . the heisenberg uncertainty principle , however , allows for phase accuracies consistent with the basic principles of quantum mechanics , as low as @xmath8 , called the heisenberg limit . the latter can be achieved by using initially entangled atoms instead of a polarized sample @xcite .    here we propose a spectroscopy technique which aims to measure the atomic interaction parameters arising form two - body collisions by using similar ramsey spectroscopic ideas , but applied to interacting atoms instead of non - interacting ones . the many - body interactions during the free evolution build in quantum mechanical correlations between the atoms and generates a large amount of entanglement even from an initially uncorrelated sample @xcite . this scheme provides a resolution which is fundamentally limited to @xmath2 , a factor of @xmath0 smaller than the heisenberg limit achievable with non - interacting atoms . moreover , for initially uncorrelated atoms the sensitivity scales as @xmath9 , implying a gain in resolution by a factor of @xmath10 with respect to the classical shot noise resolution . a more general analysis of beyond - heisenberg scaling for multi - body collisions or tensor - field interactions has been derived in ref . @xcite . many - particle entangled states , however , are difficult to prepare and maintain since they are extremely fragile : in practice , noise and decoherence rapidly collapse entangled states into classical statistical mixtures . for example in standard ramsey spectroscopy decoherence destroys the potential gain provided by entanglement , and , when it is taken into account both initially uncorrelated and maximally entangled states provide the same shot noise resolution @xcite . here we show that decoherence degrades the sensitivity of our method to @xmath11 for initially entangled atoms and to @xmath12 for initially uncorrelated ones . this scaling is interesting as it demonstrates that the entanglement build up by the many - body interactions helps , even with decoherence , to keep the sensitivity achievable with initially squeezed atoms above the maximally achievable with initially uncorrelated particles . we also discuss physical implementation of our scheme using a ) spinor bose- einstein condensates and b ) fermionic atoms loaded in optical lattices and compare and contrast the advantage and disadvantages of these two possible set - ups . the paper is organized as follows . in sec . ii we present our interferometry method based on the many - body collective dynamics . in sec . iii we apply it to initially uncorrelated states and derive the type of initially entangled states that lead to heisenberg - limited sensitivity . in this section we also propose a way to generate such states . in sec . iv we include the effect of decoherence . in sec . v we discuss the experimental implementation of the interferometry method using spinor condensates and cold fermionic atoms in optical lattices and discuss possible technical limitations . finally we conclude in sec . vi . let us consider a collection of @xmath0 spin-1/2 interacting particles described by the hamiltonian    @xmath13    we used @xmath14 to denote the collective spin operators of the @xmath0 atoms : @xmath15 , where @xmath16 and @xmath17 is a pauli operator acting on the @xmath18 atom . we set @xmath19 . our interferometric technique follows the same ideas used in standard ramsey spectroscopy , but replaces the free evolution of the atoms by evolution with @xmath20 and the goal is to estimate @xmath21 as accurately as possible . we start with the initial state @xmath22 and apply a @xmath23 rotation to all particles about the y axis , @xmath24 . after letting the system evolve for time @xmath25 under @xmath20 , @xmath26 , a second @xmath27 pulse is applied , @xmath28 . finally the collective spin is measured @xmath29 . if this scheme is repeated during a total time @xmath30 the achievable sensitivity is given by @xmath31 with variance @xmath32 . in contrast to standard ramsey spectroscopy , which is fundamentally limited by the heisenberg exclusion principle to @xmath33@xcite per realization , the fundamental limit of this scheme is @xmath34 . this follows from an application of the time - energy uncertainty principle : @xmath35 where @xmath36 is the variance of the hamiltonian and @xmath37 is the variance in estimating time from a measurement on the system . can be used to determine time with an uncertainty @xmath38 the time - energy uncertainty can be reexpressed as @xmath39 , where @xmath40 . we can establish an upper limit @xmath41 from @xmath42 is @xmath43 . this implies @xmath44 and therefore a fundamental uncertainty in determining @xmath45 after a time @xmath25 of : @xmath46 spectroscopy ------------    let us start by outlining the basic properties of @xmath20 . it commutes with @xmath47 and @xmath48 , so a good basis to describe the dynamics is the one spanned by the collective angular momentum eigenstates . in the @xmath49 dimensional hilbert space , there are @xmath50 orthogonal fully symmetric states which we denote by @xmath51 that satisfy : @xmath52 , and @xmath53 with @xmath54 . we refer to the @xmath55 states as the @xmath56 manifold . the quantum numbers @xmath57 and @xmath58 are conserved during the dynamics , and therefore , if the initial state @xmath59 belongs to @xmath56 , then the subsequent evolution of the system only takes place within @xmath56 .      first , consider the case when the initial state is fully polarized along @xmath60 : @xmath61 . after the first @xmath23 pulse it becomes @xmath62 . the effect of @xmath20 on the state is to imprint a phase @xmath63 to each @xmath64 component . in the limit of many atoms @xmath65 we may approximate the sums by integrals assuming @xmath66 is a continuous variable and to replace the binomial coefficient by a gaussian distribution with the same width . in this limit @xmath67 and @xmath68 become @xmath69 where we used the notation @xmath70 . the time evolution of the signal is a series of equally spaced gaussian pulses ( see fig . [ fig1 ] ) . as the system evolves , at first the different accumulated phases lead to a collapse of @xmath71 . for the initial coherent state in consideration . the collapse time , which depends on its variance , is @xmath72 . the evolution continues and at @xmath73 all the different components rephase inducing a perfect revival of the initial state , now with opposite polarization:@xmath74 . using eqs . ( [ jz ] ) and ( [ jz2 ] ) one can show that the optimal sensitivity is reached at the points with maximal signal intensity , @xmath75 . it scales as @xmath9 ( see fig . [ fig1 ] ) : @xmath76      entangled atomic states ( e.g. spin squeezed states ) potentially allow to significantly improve the sensitivity in precision measurements . for example , consider the initial state@xmath77 \right ) } { \sqrt{1+a^{2 } } } . \label{inisquee}\ ] ] here @xmath78 is a real number of order unity and @xmath0 is assumed to be even and large . the probability distribution of @xmath79 in the @xmath80 basis is plotted in the inset of fig.[fig2 ] . for even @xmath58 values , the amplitudes are almost constant ( except for those with @xmath81 close to @xmath82 ) , and for odd @xmath58 values they are proportional to @xmath83 . @xmath84 is a highly squeezed state with reduced variance in the @xmath60 direction : @xmath85 .    using properties of the wigner rotation functions , it is possible to show that :    @xmath86}{n ( \chi t-\pi/2 + k\pi ) } , \label{jzs } \\ \langle   \hat{j}_x^2\rangle_\psi^a & = & n^2a\sum_{k=0 } \left(\frac{% j_1[2n(\chi t- k \pi/2]}{(2n ( \chi t- k\pi/2))}\right ) + b , \label{jz2s}\end{aligned}\ ] ]    with @xmath87 and @xmath88 . as shown in fig . [ fig2 ] the system starts at @xmath89 with @xmath90 , and during the time evolution @xmath91 grows while oscillating . at @xmath92 the signal reaches a maximum and vanishes again at @xmath93 . then it evolves with opposite polarization until it returns to the initial state at @xmath94 . using eqs . ( [ jzs ] ) and ( [ jz2s ] ) , the optimal sensitivity can be shown to be reached at @xmath95 , with @xmath96 an odd integer ( see fig.[fig3 ] ) : @xmath97 these states lead to @xmath98 sensitivity . the best accuracy is obtained for @xmath99 , in which case @xmath100 . note , however , that estimating @xmath101 by measuring @xmath102 might be experimentally impractical due to the fact that not only the signal vanishes at @xmath95 ( regardless of the value of @xmath78 ) but also @xmath103 is very sharply peaked at @xmath95 and hence very sensitive to small variations of @xmath104 . to overcome this limitation we propose , following the ideas exposed in ref . @xcite , to measure @xmath105 instead of @xmath102 . in contrast to @xmath102 , the average square signal does not vanishes for finite @xmath78 ( it goes like @xmath106 ) and , as shown in fig . [ fig3 ] , measuring @xmath101 by means of @xmath107 provides the same @xmath2 accuracy with the advantage of a broader profile around the optimal value which becomes slightly shifted from @xmath108 . we calculated @xmath103 numerically according to eq . ( [ sig ] ) but replacing @xmath102 by @xmath109 . this idea , however , has the drawback that @xmath110 diverges exactly at @xmath108 . nevertheless as shown in fig . 3 this is a extremely narrow divergency and almost a one point behavior .    let us now discuss how to generate states of the form given by eq.([inisquee ] ) . these states are highly squeezed , and a robust method to generate them can have also applications in heisenberg limited phase detection @xcite . the idea is to adiabatically convert @xmath111 into @xmath112 by gradually increasing @xmath113 from @xmath114 to @xmath115 in the following time dependent hamiltonian@xcite : @xmath116the first term proportional to @xmath117 is required to keep the system within @xmath56 . as shown in fig . 4 , there is always a finite gap between the ( non - degenerate ) ground state . ] and the lowest excited state , so the adiabatic passage can be satisfied . the experimental implementation eq.([hadiabatic ] ) will be discussed in sec.vb    once @xmath118 is obtained , the @xmath119 components can be generated by evolving @xmath118 with the hamiltonian @xmath120 . as proposed in ref . @xcite , @xmath121 can be implemented by using two laser beams with frequency difference twice the hyperfine level splitting which coherently flips pairs of aligned spins . comparable to unity , excitations to states with @xmath122 , should also be considered . ] in realistic experiments decoherence effects are inevitably present . the main type of decoherence is dephasing due to processes that cause random changes in the relative phase of the quantum states , while preserving the total atomic population in the atomic levels . among the mechanisms that result in dephasing effects one can consider collisions , stray fields , and laser instabilities . we model the phase decoherence by adding a term to the hamiltonian given by @xmath123 where @xmath124 are independent random gaussian variables with zero mean and white noise correlation functions @xmath125 that represent couplings with the environment . here the bar denotes averaging over the different outcomes of the random variables . this model is equivalent to considering the time evolution of the reduced density operator for a single ion @xmath126 as given by the master equation @xmath127-\gamma /2(\sigma _ { z}\hat{\rho}\sigma _ { z}-\hat{\rho})$ ] .    to include the effect of decoherence it is simpler to go to the uncoupled basis which diagonalizes both @xmath128 and @xmath20 . each state in this basis can be label as @xmath129 , where @xmath130 and @xmath131 . using the uncorrelated nature of the @xmath132 variables one can show by straightforward calculations that @xcite    @xmath133    with @xmath134 and@xmath135 . note that these equations are valid for an arbitrary initial state . the noise induces an exponential decay of the average signal and its higher moments . however , unlike the average signal the higher moments are also affected by the presence of diffusion terms .    if we assume an initially coherent state with @xmath136 and @xmath137 given by eqs . ( [ jz ] ) and ( [ jz2 ] ) , one can calculate from eqs . ( [ dec1 ] ) and ( [ dec2])the sensitivity in the presence of decoherence . in contrast to the ideal dynamics where the longer the interrogation time @xmath25 the better is the phase accuracy , when decoherence is accounted for , one has to optimize with respect to both @xmath25 and @xmath104 , in order to obtain the best sensitivity . provided that @xmath138 , standard minimization yields :    @xmath139    at @xmath140{\frac{2}{n^3 } ( e^{\tau_{opt}}-1)}$ ] and @xmath141 . note that even with initially uncorrelated atoms decoherence reduces the sensitivity by a factor of @xmath142 . this has to be contrasted with standard ramsey spectroscopy done with non - interacting atoms , where the phase resolution with initially uncorrelated atoms has the same classical shot noise resolution so is not too much affected by decoherence . the reason of the different behavior is the fact that @xmath143 builds up particle correlations during the evolution and transforms uncorrelated states into highly squeezed ones which are vulnerable to decoherence . more specifically , as shown in ref . @xcite , @xmath143 ideally transforms initially coherent states into maximally correlated @xmath0-particle ghz states @xcite after a period of evolution of @xmath144 . these states are very fragile and decohere @xmath0 times faster than uncorrelated particles . let us now consider the case when the initial state is @xmath145 . using eqs . ( [ jzs ] ) and ( [ jz2s ] ) in eqs . ( [ dec1 ] ) and ( [ dec2 ] ) , one gets an optimal sensitivity , provided @xmath146 , given by : @xmath147 at @xmath148 , and @xmath149 . ( [ signdecma ] ) is very interesting as it demonstrates that by using many - body interactions , entangled states can outperform the signal to noise limit of classically unentangled states even in the presence of decoherence . the sensitivity of the squeezed states is a factor @xmath150 larger than the maximal achievable with uncorrelated atoms , see eq.([soptd ] ) . in fig . 3 we plot the optimal sensitivity in the presence of decoherence both by measuring @xmath102 and @xmath109 . for practical reasons , as the ones discussed above , measuring @xmath109 instead of @xmath102 might be experimentally convenient also with decoherence . moreover , note that while states with @xmath151 provide the best accuracy in the ideal dynamics , states with @xmath152 are optimal with decoherence . @xmath153 describes the spin dynamics of a @xmath154 bose einstein condensate , if one associates the internal @xmath155 hyperfine states of the atoms as the two states of an effective spin-1/2 particle ( the @xmath156 state can be decoupled by introducing a detuning much larger than the chemical potential ) @xcite . in this system the coupling strength @xmath21 is proportional to @xmath157 with @xmath158 the s wave scattering length between particles of type @xmath159 and @xmath160 , @xmath161 the atom mass and @xmath162 the volume .    due to the dependence of @xmath101 on the elastic scattering length difference , @xmath163 , using an inter - state feshbach resonance @xcite , such as the one predicted between the @xmath164 and @xmath154 states in @xmath165 atoms , one could widely tune @xmath101 by varying an external magnetic field . close to the feshbach resonance there is a singularity in the scattering length , @xmath166 with @xmath167 being the resonance width , @xmath168 the zero energy resonance and @xmath169 the background scattering length . this singularity is due to the near degeneracy of the collision energy of the atoms with the binding energy of an extremely loose diatomic molecular state . one can use our spectroscopy technique to measure @xmath170 in such system . however , in order for the method to work , it is important to reduce two and three body particle losses across the resonance , as they can considerably affect the sensitivity . we estimate for example the effect of two - body losses by solving the following master equation :    @xmath171-\gamma /2(\hat{c}^\dagger \hat{c } \hat{\rho}+ \hat{\rho } \hat{c}^\dagger \hat{c}- 2\hat{c}^\dagger \hat{\rho } \hat{c } ) , \ ] ]    with @xmath172 the operator that destroys a pair atoms in different spin states and which account for inelastic spin relaxation processes across the resonance . here the constant @xmath173 is related to the two body loss rate density , @xmath174 , which satisfies the relation @xmath175 @xmath176 depends on the imaginary part of @xmath177 as @xmath178 and @xmath179 varies with the applied magnetic field , the difference in magnetic moment between the feshbach resonance state and a pair of atoms in the entrance channel , @xmath180 and the inverse molecular lifetime , @xmath181 , as @xcite :    @xmath182 ^ 2}.\ ] ]    in fig.5 we summarize our numerical solutions by plotting the scaling exponent @xmath183 of the optimal sensitivity with the initial number of atoms , @xmath184 = - b \log[n(t=0 ) ] + { \rm const}$ ] , versus @xmath185 for both initially uncorrelated atoms and initially squeezed atoms . the plot demonstrates the drastic reduction of the sensitivity with @xmath173 . for both initially uncorrelated atoms and initially squeezed atoms one requires to limit @xmath186 to have at least shot noise resolution , i.e @xmath187 . experimentally both two and three body particle losses contribute . as three body losses scale as the four power of the scattering length @xcite they decrease faster than two - body losses away from the resonance . as such , the sensitivity of the method can be optimized by , on one hand , tuning the magnetic field far enough from the resonance so that mainly two - body losses are relevant and , on the other , by reducing @xmath176 by using a very narrow feshbach resonance , and by considering species of atoms that belong to the lowest magnetic manifold so they posses long relaxation times . recently the dynamics of bosonic atoms with two relevant internal states loaded in a deep optical lattice has been used to perform precision measurements of atomic scattering properties @xcite . explicitly , the experiment detected the modification in the ramsey fringes caused by frequency shifts induced by interatomic interactions in wells occupied with two atoms , and used it to measure the elastic scattering length difference . as the various wells in the optical lattice behave as independent wells due to the suppressed tunneling between them , the experimental phase sensitivity was limited to the shot noise resolution associated to statistical fluctuations in a finite number of wells . in practice however , the sensitivity was much lower due to the two - body losses close to the feshbach resonance in the double occupied wells . we now discuss how to generalize these ideas by using the collective spin dynamics ( instead of the single well one ) of fermionic atoms in an optical lattice loaded with one atom per site . specifically we demonstrate how to engineer an effective @xmath188 hamiltonian in such systems . the latter can be used to perform measurements of the scattering properties via the entanglement dynamics but without the limitation of particle losses at the resonance . the effective dynamics of a system of spinor ultra - cold atoms confined in a unit filled lattice , deep in the mott insulator regime @xcite can be described in terms of effective spin operators , @xmath189 , @xmath190 , where @xmath191 are the annihilation operator of a particle of type @xmath192 at site @xmath193 , and @xmath194 are number operators . in such spin variables the effective hamiltonian maps to a spin xxz hamiltonian @xcite : @xmath195    here @xmath196 means that the sum is over nearest neighbors , @xmath197 are spin - dependent tunneling energies and @xmath198 are on - site interaction energies proportional to the s - wave scattering lengths between the various components . for fermions @xmath199 due to the pauli exclusion principle . the heisenberg , @xmath200 , and ising @xmath201 , coupling constants are given by    @xmath202    where the upper and lower signs are for fermionic and bosonic atoms respectively .    using spin dependent optical potentials one can tune the various coupling constants and in particular one can engineer the condition @xmath203 . in this limit , we may treat the effect of the ising term , @xmath201 , by means of perturbation theory . assuming that at @xmath89 the initial state is prepared within the @xmath55 manifold , @xmath56 , a perturbative analysis predicts that for times @xmath25 such that @xmath204 , @xmath205 confines the dynamics to @xmath56 and transitions outside it can be neglected . as a consequence , only the projection of @xmath206 on it , which corresponds to @xmath207 is effective and @xmath208 acts as a long range hamiltonian @xcite . in eq . ( [ impl ] ) @xmath60 is the number of nearest neighbors . with this effective hamiltonian the spectroscopic method described above can be applied in lattice systems to measure interactions .    moreover , due to the fact that the large heisenberg term restricts the dynamics to the @xmath56 manifold , @xmath200 acts as an effective @xmath209 . therefore , by adding a @xmath210 term via a raman transition with effective two - photon rabi frequency @xmath211 @xcite one can also implement in such lattice systems eq . ( [ hadiabatic ] ) . this hamiltonian allows , as described above , for the squeezed state preparation . the realization of @xmath212 in unit filled optical lattices has the strong advantage that two and three body losses are suppressed and thus they are not longer a limitation in these systems . it is important to note , however , that only fermionic atoms can be used if @xmath101 is going to be tuned by means of a feshbach resonance . the reason is that for bosons , as @xmath213 grows , @xmath101 also grows and at some point the condition @xmath203 is violated . on the other hand , for fermions , as @xmath213 grows , the ratio @xmath214 remains constant . the drawback of the lattice implementation is the factor of @xmath0 in the denominator of the effective @xmath101 , because the slower dynamics limits the sensitivity of our method . typical experiments that use lattices deep enough for the effective spin model to be valid , work in a parameter regime where @xmath215 . using spin dependent optical lattices @xcite one can control the tunneling rates and set @xmath216 . additionally close to the resonance @xmath217 can be enhanced up to about @xmath218 times its off resonance value @xcite . therefore the revival time in these set - ups varies with the magnetic field between @xmath219 ( @xmath218 ms - @xmath220ms ) . thus , in order to keep the measurement time in a reasonable experimental time scale , one should limit the experiment to 1d lattice systems which have of the order of @xmath221 atoms per tube @xcite . the scalability problem certainly limits the achievable sensitivity . nevertheless , even with these reduced number of atoms , if initially squeezed atoms are used , the phase accuracy @xmath222 is within the @xmath223 sensitivity required to test the proton - electron mass ratio variation on the level of @xmath224 , in a narrow feshbach resonance ( @xmath225 mg)@xcite . in summary , here we proposed an inteferometric method that relays on the quantum dynamics of interacting spins to perform precision measurements of the scattering properties with a resolution fundamentally limited to @xmath2 . we discussed the class of squeezed states that are required to achieve heisenberg sensitivity and proposed a method for its generation . we studied the effect of decoherence and showed that it affects the achievable sensitivity even with initially uncorrelated atoms as they become vulnerable due the many - particle entanglement builded during the many - body dynamics . we also showed that the dynamically induced quantum correlations help to keep the sensitivity of initially squeezed states above the uncorrelated atom threshold even with decoherence . this has to be contrasted with standard ramsey spectroscopy with non - interacting atoms where all the potential gain in sensitivity due to initial quantum correlation is lost in the presence of decoherence . finally we discussed possible physical implementation of the proposed spectroscopy using spinor condensates and fermions in unit filled optical lattices . we showed that particle losses close to the feshbach resonance limits the sensitivity in the bec implementation , and thus the use of very narrow feshbach resonances and atomic species with a large spin relaxation time are required to optimize the sensitivity of the method . our implementation of the scheme in optical lattice set - ups is insensitive to atomic losses but on the other is limited by the slow dynamics . regardless of such problems , the proposed spectroscopy is a new scheme that takes advantages of many - particle entanglement to perform high precision measurements of the interaction parameters in spin mixtures . we are optimistic that emerging technology in optical and magnetic feshbach resonances might overcome the current experimental limitations in such systems or that new developments in the field might provide alternative set - ups for the experimental application of the presented method and make of it an useful spectroscopic technique . this work was supported by itamp , nsf ( career program ) , afosr and the david and lucille packard foundation .
Abstract:

Answer: we propose a method to perform precision measurements of the interaction parameters in systems of @xmath0 ultra - cold spin @xmath1 atoms . the spectroscopy is realized by first creating a coherent spin superposition of the two relevant internal states of each atom and then letting the atoms evolve under a squeezing hamiltonian . the non - linear nature of the hamiltonian decreases the fundamental limit imposed by the heisenberg uncertainty principle to @xmath2 , a factor of @xmath0 smaller than the fundamental limit achievable with non - interacting atoms . we study the effect of decoherence and show that even with decoherence , entangled states can outperform the signal to noise limit of non - entangled states . we present two possible experimental implementations of the method using bose - einstein spinor condensates and fermionic atoms loaded in optical lattices and discuss their advantages and disadvantages .


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 transition metals and their oxides are well known for exhibiting a variety of magnetic ordering phenomena @xcite . the nature of the low - temperature magnetically ordered states depends on the type of dominant magnetic exchange interactions , which in turn depend on the details of the crystal structure and electronic band structure . the theoretical models that describe magnetism in solids can be divided into three broad categories : ( i ) models that do not _ apriori _ contain a local moment but allow for moment formation via electron - electron interactions , such as the hubbard model @xcite , ( ii ) models consisting of localized spins only , such as the heisenberg model , and ( iii ) models that have itinerant electrons coupled to localized moments , e.g. , the kondo lattice model @xcite . starting with a general multi - orbital hubbard model , the other two types can be obtained in appropriate limits @xcite . hence , the hubbard model is the elementary model for describing magnetism . indeed , two of the most common magnetically ordered ground states , the ferromagnet and the staggered antiferromagnet , are present in the mean - field phase diagram of the hubbard model @xcite . search for non - collinear and non - coplanar magnetic ordering in models and materials has emerged as an important research topic in recent years . the reason being the fundamental connection between the nature of magnetic order and electrical properties of various magnetic materials . a non - coplanar magnetic order is known to give rise to anomalous hall response in transport @xcite . on the other hand , a planar spiral order can allow for a ferroelectric response via the spin - current mechanism @xcite . recent theoretical studies on kondo - lattice models have shown that a variety of unconventional magnetic ground states can be stabilized depending on the underlying lattice geometry and electronic filling fraction . in particular , triangular and checkerboard lattices allow for non - coplanar magnetism at quarter filling of the band @xcite . more recently , existence of spin - charge ordered phases has also been reported at the average filling fraction of two electrons per three sites @xcite . given that the hubbard model is the fundamental model for magnetism , it is important to know if the unconventional magnetic phases found in the kondo - lattice hamiltonian are also present in the hubbard model . moreover , for some of the magnetic phases the shape of the fermi surface plays a crucial role and the ordered phases should therefore be independent of the nature of interactions @xcite . to the best of our knowledge a systematic search across the electronic density range for non - collinear and non - coplanar magnetic phases in the hubbard model on triangular lattice has not been reported .    in this work , we map out the magnetic phase diagram of the hubbard model on anisotropic triangular lattice within a variational mean - field approach that captures the non - collinear and non - coplanar magnetism on same footing as the collinear magnetism . the phase diagram is rich and consists of many unconventional phases that are found in the corresponding kondo - lattice model . these include , ( i ) the four - sublattice 3q order at quarter filling , ( ii ) the flux state at quarter filling , ( iii ) the non - collinear charge ordered phase at two - third filling . in addition , the 120@xmath0state is stable near the isotropic triangular limit and a collinear stripe state is stable in the anisotropic regime . the influence of nearest neighbor coulomb repulsion on these phases is also investigated . we find that a pinball - liquid like charge ordering phase , which has so far been reported in spinless fermion models or multi - orbital models on triangular lattice @xcite , is stabilized by the coulomb interactions at a filling fraction of @xmath2 and occurs concomitantly with either collinear or non - collinear magnetic order . given the fundamental nature of the hubbard hamiltonian , our results are of general interest . in addition , the reported phases are of specific relevance to various triangular lattice systems , such as the layered agnio@xmath3 @xcite , the nacoo@xmath3 @xcite and the organic charge - transfer salts @xcite . the remainder of the paper is organized as follows . in section ii we describe the model and briefly discuss the previous investigations that focus on non - collinear magnetism . in section iii we describe the method used in this work . results and discussions follow in section iv , where we first discuss the density versus interaction phase diagram , and then discuss the phase diagrams corresponding to specific electronic filling fractions . in section v we analyze the effect of nearest neighbor coulomb repulsion on the magnetic phase diagram . conclusions are presented in section vi . the single - orbital hubbard hamiltonian on an anisotropic triangular lattice is given by ,    @xmath4    where @xmath5 is the creation ( annihilation ) operator , and @xmath6 is the number operator for electrons with spin @xmath7 at site @xmath8 . the total number operator for electrons at site @xmath8 is @xmath9 , and @xmath10 is the chemical potential . the electronic hopping amplitudes between nearest neighbor sites on triangular lattice are @xmath11 and @xmath12 as shown in fig . 1 , and @xmath13 denotes the on - site hubbard repulsion . the hubbard model is one of the most successful and well studied hamiltonians in condensed matter physics @xcite . however , due to its relevance to various materials , certain filling fractions have been of more interest than others . for instance , the model on square lattice ( @xmath14 ) has been studied in great detail near half filling . very few studies have explored the possibility for non - collinear and non - coplanar phases @xcite . the half - filled case on triangular lattice has also been explored in search of spin liquid phases , as the large @xmath13 limit of the hubbard model leads to the spin-@xmath15 heisenberg model which is frustrated on the triangular geometry . there have been studies for the quarter filled case , which describes the physics of organic charge transfer salts @xcite . the filling fraction of @xmath2 has also been of interest due to its relevance to layered agnio@xmath3 and nacoo@xmath3 . a systematic search for non - collinear and , in particular , non - coplanar phases over the full range of electronic densities has not been reported even at the mean field level . we determine the ground state phase diagrams by performing a systematic search in the phase space of various ordered magnetic configurations . this is achieved by making use of the rotational invariance of the hubbard hamiltonian @xcite . before describing the details of the rotationally invariant scheme , let us recall the commonly used unrestricted hartree fock ( uhf ) method @xcite . in this method , the interaction term is decoupled into charge and spin sectors . ignoring the 2nd order term in fluctuations , the hamiltonian reduces to the well known form ,    @xmath16 \nonumber   \\     & & - \lambda u \sum\limits_{i } [ \langle s^{+}_{i } \rangle s^{-}_i+     \langle s^{-}_{i } \rangle s^{+}_i-\langle s^{+}_{i}\rangle \langle s^{-}_{i } \rangle],\end{aligned}\ ] ]    where , @xmath17 and @xmath18 are the spin operators . the variational parameter @xmath19 determines the relative contribution of the hartree and the fock terms in the mean - field decoupling .    in the uhf approach , the site dependence of the quantum averages is retained . starting from a random guess for the various mean - field parameters , the hamiltonian is diagonalized iteratively until self - consistency is achieved . in principle , the self - consistent solution depends on the starting configuration of mean - field parameters . therefore , it is required to use a large number of starting configurations . the reliability of the final solution obtained in this way depends to a large extent on the complexity of the energy landscape . in addition , @xmath19 should be determined via an energy minimization over the various self - consistent solutions . although , this scheme can capture the conventional ordered magnetic phases very well , it can easily miss out on non - collinear and non - coplanar magnetic orderings . motivated by the appearance of such unusual magnetically ordered phases in a kondo lattice model , we search for similar magnetic phases in the hubbard model . for this purpose , we recast the mean - field decouplings in a way that allows for a systematical search of a variety of ordered phases . we begin by re - writing the hamiltonian in a reference frame with site - dependent spin - quantization axes @xcite . this is done by performing local su(2 ) rotations of the quantization axes , given by ,    @xmath20 & = & \left [ \begin{array}{cc }   cos(\frac{\theta_{i}}{2 } ) e^{i\frac{\phi_{i}}{2 } } &   -sin(\frac{\theta_{i}}{2 } ) e^{i\frac{\phi_{i}}{2}}\\   sin(\frac{\theta_{i}}{2 } ) e^{-i\frac{\phi_{i}}{2 } } & cos(\frac{\theta_{i}}{2 } ) e^{-i\frac{\phi_{i}}{2}}\\   \end{array } \right ] \left [ \begin{array}{c } d_{ip}\\ d_{ia}\\ \end{array } \right ] \nonumber \\ \vspace{0.5 cm }   & \equiv & { \cal r}(\theta_{i } , \phi_{i } ) \left [ \begin{array}{c } d_{ip}\\ d_{ia}\\ \end{array } \right],\end{aligned}\ ] ]    where @xmath21 are annihilation operators for electron at site @xmath8 with spin parallel ( antiparallel ) to the quantization axis . a simple hartree decoupling of the transformed hamiltonian leads to , @xmath22 , \end{aligned}\ ] ]    where , @xmath23 and @xmath24 can take two values . the coefficients @xmath25 are explicitly given by ,    @xmath26 = { \cal r}^{\dag}(\theta_{i } , \phi_{i } ) . { \cal r}(\theta_{j } , \phi_{j})\end{aligned}\ ] ]    we consider two broad categories of variational states . first one is the set of general non - coplanar spiral states which can be parameterized by a cone angle @xmath27 and a spiral wave vector @xmath28 @xcite . the polar and azimuthal angles of the spins are then given by , @xmath29 and @xmath30 . the other set of states are the block - periodic states , where the spin configuration of a @xmath31 block is repeated periodically to generate the configuration over the entire lattice . the choice of these variational states is motivated by the presence of various long - range ordered phases in the corresponding kondo - lattice model . the flux state and the non - coplanar states shown in fig . [ fig1 ] are examples of the @xmath32 block - periodic states . in the flux state the azimuthal angles of the neighboring spins differ by @xmath33 while the polar angle remains equal to @xmath33 . the non - coplanar state can be obtained from the flux state by changing the polar angles of the four spins as shown in fig . similarly we take a @xmath34 block to capture the 120@xmath0state , which can be connected to a ferromagnetic state by varying the polar angle from @xmath33 to 0 . in order to allow for the six - sublattice non - collinear charge ordered ( nc - co ) state at @xmath2 , we take a @xmath35 cluster as the building block @xcite . the calculations have been performed in both canonical and grand canonical approaches . in the grand canonical ( canonical ) approach , we seek a self - consistent solution for a given fixed value of chemical potential ( density ) for all the variational states discussed above . in cases where more than one self - consistent solutions exists , the lower energy solution is selected and the corresponding state is taken as the mean - field ground state at that value of the chemical potential ( density ) . for the grand canonical approach , the density corresponding to the ground state solution is calculated at the given value of the chemical potential and the phase diagrams are presented with density as a parameter . in this section we present the results obtained by using the mean - field decoupling scheme described in the previous section . for most of the calculations we used @xmath36 @xmath37-points in the first brillouin zone . results have also been checked for @xmath38 in some cases . in the following subsection we discuss the @xmath1-@xmath13 phase diagrams for different values of @xmath39 . the @xmath1-@xmath13 phase diagrams are obtained within the grand canonical approach . the energy minimization over the two sets of variational states is performed for a fixed value of chemical potential @xmath10 . the chemical potential is then systematically varied in order to obtain different average electronic densities . the total energy is computed at zero temperature , where the fermi function simply becomes a step function . the grand canonical approach has an advantage that it allows for the phase separation regions , which commonly arise in electronic systems , to be easily captured .    for @xmath40 , the hopping connectivity is that of a square lattice . therefore , in this case we display the density range @xmath41 in the phase diagram as the system is particle - hole symmetric ( see fig . [ fig2 ] ( @xmath42 ) ) . the square lattice phase diagram is dominated by three simple phases : a fermi liquid state for weak to intermediate @xmath13 and away from @xmath43 , an antiferromagnetic ( afm ) state near @xmath43 and a ferromagnetic ( fm ) state for large @xmath13 and away from @xmath43 . in addition , in the intermediate @xmath13 range we find a flux phase and a narrow window of phase separation close to @xmath44 . the flux phase is known to exits in the kondo - lattice model with a classical approximation for local moments @xcite . it has also been reported recently in he uhf study of a five - orbital hubbard model for iron pnictides @xcite . the planar spiral states are stable over a wide range of parameters . next , we present the phase diagram for the anisotropic case where the particle - hole symmetry does not hold . therefore , the density range displayed in fig . [ fig2 ] ( @xmath45)-(@xmath46 ) is @xmath47 . the fm state is strongly suppressed by the triangular anisotropy in the low - density regime , whereas it still dominates the large density regime . this can be explained by invoking the stoner picture for ferromagnetism : the triangular lattice density of states is higher if the fermi level is located close to @xmath48 and therefore the stoner criterion for ferromagnetism can be satisfied for relatively small values of @xmath13 . the flux state remains stable in a narrow window around intermediate @xmath13 and @xmath44 . a new magnetically ordered state enters the phase diagram near @xmath43 . this is the well known @xmath49-state which is the classical ground state of a triangular lattice heisenberg model . a wider window of phase separation appears in @xmath50 regime @xcite , where a conical spiral state also becomes stable . the phase diagram remains qualitatively similar between @xmath51 and @xmath52 . however , the isotropic triangular case ( @xmath53 ) gives rise to new and interesting phases . the first of these new phases is the 4-sublattice non - coplanar order near @xmath44 . this unusual magnetic state , also known as 3q state , is of wide interest as it supports a non - vanishing scalar spin chirality and an associated quantized hall conductivity @xcite . the quantized hall conductivity is related to the nontrivial topology of the electronic band structure that the 3q state induces in the kondo hamiltonians @xcite . interestingly , the 3q state was predicted to be stable in the hubbard model at @xmath48 based on fermi surface nesting arguments @xcite . we find that the 3q state is stable at @xmath48 only for @xmath54 and loses to fm order for larger values of @xmath13 . the second new state that appears exactly at @xmath2 is the recently reported @xmath55-sublattice non - collinear charge ordered ( nc - co ) state . in addition , the @xmath56 ordered afm state completely disappears , and instead the @xmath49-state dominates the phase diagram near @xmath43 , as expected @xcite . having discussed the overall structure of the @xmath1-@xmath13 phase diagrams , we now examine more closely the effect of anisotropic hopping connectivity on different magnetic phases at specific filling fractions .      in this subsection we describe the @xmath12-@xmath13 phase diagrams obtained within canonical ensemble approach , where the minimization is performed over the variational parameters for a fixed value of the average electronic density @xmath1 . we discuss the evolution of various phases as the lattice connectivity changes from square type to triangular type . the non - collinear flux phase is lowest energy state for unfrustrated square lattice ( @xmath60 ) at @xmath44 for intermediate coupling strength . this state remains stable in a wide parameter regime , @xmath61 . in fact , the stability window of the flux state widens in terms of the @xmath13 values upon increasing @xmath39 ( see fig . [ fig3](@xmath42 ) ) . for @xmath62 , a non - coplanar state becomes stable which evolves into the 3q state at @xmath63 . the evolution from the flux state to the 3q non - coplanar state is depicted by the change in variational angle @xmath64 ( see fig . [ fig4 ] ( @xmath42 ) ) . for the flux state @xmath65 , whereas for the perfect 3q state @xmath66 . the other effect of the triangular geometry for @xmath44 is to suppress the fm phase . the fm phase give way to planar spiral phases even for small values of the parameter @xmath39 . the transition from a non - magnetic fermi liquid state to a magnetic state also occurs at a slightly larger value of @xmath13 for the triangular lattice in comparison to the square lattice . another unusual six - sublattice magnetic order was recently reported in the kondo - lattice model for average density of two electrons per three sites @xcite . in order to explore the possibility of this phase in the hubbard model , we focus on the filling fraction of @xmath2 ( see fig . [ fig3](@xmath45 ) ) . we again find that the fm state is suppressed in favor of the planar spiral states . additionally the @xmath49 order is stable for @xmath67 . it is interesting to note that the stability window of the @xmath49 state at @xmath2 is even wider than that at @xmath43 . importantly , we find that the six - sublattice nc - co order is the ground state for @xmath68 and for @xmath69 . the half - filling , @xmath43 , is the most commonly studied filling fraction in the hubbard models on triangular and square lattice @xcite . within the mean - field , we find that the afm order is robust and is stable in the regime @xmath70 . beyond @xmath52 , the afm order continuously evolves towards the @xmath49 state . in fig . [ fig4](@xmath42 ) , we show the change in the spiral wave - vector @xmath71 from @xmath56 to @xmath72 . the half - filling case is also of interest for the possible existence of a spin - liquid state @xcite . however , our mean - field method does not capture the spin - liquid states and therefore we can not comment on this competition in the present paper . finally we discuss the filling fraction @xmath48 , which is interesting due to the presence of fermi - surface nesting feature in the non - interacting hamiltonian . in fact , the non - coplanar 3q state was first predicted based on the nesting property of the fermi surface . clearly , for small @xmath13 the 3q state is the ground state due to the presence of a peculiar fermi surface nesting in the non - interacting hamiltonian @xcite . interestingly , the 3q order loses to fm order for @xmath73 . the relative stability of the 3q state at @xmath44 and @xmath48 was also explored in the kondo - lattice hamiltonian with a conclusion that the window of stability of the 3q state is wider for @xmath44 compared to @xmath48 @xcite . our calculations on the hubbard model support the results of the corresponding kondo - lattice model . in sharp contrast to the quarter filling ( @xmath44 ) , the three - quarter filling favors fm state for increasing @xmath39 . as mentioned earlier , this can be understood from the particle - hole asymmetry and large dos in the @xmath50 regime for the triangular lattice . the electronic spectrum corresponding to different phases is displayed in the plots for the density of states in fig . [ fig4](@xmath45)-(@xmath74 ) . the dos is defined as , @xmath75},\end{aligned}\ ] ]    where @xmath76 are the eigenvalues corresponding to the lowest - energy self - consistent solution , @xmath77 is the dirac delta function and @xmath78 is the broadening parameter which we take to be @xmath79 for calculations . all the three possible varieties of the spectra are present : the spiral states are metallic , the 3q non - coplanar , the @xmath56 afm and the nc - co states are insulating , and the flux state is semimettalic with a graphene - like dos . a robust gap in the dos is suggestive of the stability of these phases against higher - order quantum corrections . in this section we analyze the influence of nearest neighbor coulomb repulsion on various magnetic phases discussed so far . we begin by extending the hubbard model eq . ( 1 ) by including a nn repulsion . the extended hubbard hamiltonian is given by ,    @xmath80    the parameters @xmath81 and @xmath82 denote the strength of the coulomb repulsion along inequivalent directions on the anisotropic triangular lattice . it is well known that nn coulomb repulsion in bipartite lattices favors charge ordering @xcite . however , the triangular geometry frustrates the checkerboard type charge ordering and leads to exotic charge ordered phases , such as the pinball liquid and inverse pinball liquid @xcite . here our primary focus is to study the effect of nn coulomb repulsions on the various magnetic phases found in the hubbard model . we present results for the quarter - filled case where the interesting mean - field ground states are the flux state , for @xmath83 , and non - coplanar states for @xmath62 . we use the variational scheme corresponding to the block - periodic states discussed in sec . this allows for four inequivalent sites in terms of the spin and charge structures . the planar structure is fixed to be of the flux type and the polar angle is varied to obtain magnetic phases that interpolate between the flux and the ( @xmath84 ) afm state via the 3q state ( see fig . [ fig1 ] ) . results for @xmath85 are summarized in fig . we find that the flux state ( @xmath40 ) is compatible with stripe charge ordering . indeed , the charge disproportionation @xmath86 becomes finite as soon as @xmath87 ( see fig . [ fig5 ] ( @xmath42 ) ) . the magnetic structure retains the flux character , however , the magnetic moments become unequal on sites @xmath88 and @xmath89 ( see fig . [ fig7 ] ( @xmath42)-(@xmath90 ) ) . for very large values of @xmath81 the magnetic moment becomes vanishingly small at one of the sites , leading to a strongly charge ordered afm state . in fact , a similar state has been reported in the quarter - filled hubbard - holstein , and kondo lattice models on square lattice @xcite . the nature of the magnetic order is depicted by the value of the variational polar angle @xmath64 corresponding to the minimum energy self - consistent solution . @xmath65 remains constant for @xmath40 ( see inset in fig [ fig5 ] ( @xmath46 ) ) . in contrast to the flux state , the 3q state competes with charge ordering . it is easy to understand this difference from the energetics : charge ordering leads to lowering of energy via opening an energy gap in the spectrum . since the flux state is gapless , a gap opening lowers the energy and therefore the ground state develops a charge ordering even for small values of @xmath81 . however , the 3q state is already gapped with the origin of the gap tied to peculiar magnetic structure . therefore , in order to further lower the energy the charge ordering must open a gap that off - sets the 3q gap . starting with @xmath63 , the charge disproportionation remains zero until @xmath91 . a weak reduction of the magnetic moment occurs near @xmath92 together with a deviation from the ideal 3q structure . near @xmath52 , @xmath86 becomes finite , and the magnetic moment at site @xmath88 vanishes , leading to the same charge ordered afm state as for @xmath14 . therefore , in this case the 3q state is destabilized at the onset of charge ordering . the density of states ( fig . [ fig5](@xmath90)-(@xmath46 ) ) show that a gap opens in the electronic spectrum owing to the charge ordering . thus , the semimettalic flux state is turned into a fully gapped flux state , and eventually into a gapped collinear afm for large @xmath81 . the results for the isotropic repulsive ( @xmath93 ) case are qualitatively similar to those discussed above . the main difference is that the charge disproportionation is smaller for @xmath93 compared to @xmath94 . finally , we study the effect of coulomb repulsion on the 120@xmath0phase . the results are presented in fig . the 120@xmath0state is destabilized in favor of a collinear state containing three inequivalent sites . two of the sites have equal charge density which is larger than that on the third ( @xmath95 ) . the magnetic moments on the low - charge sites vanish , while those on the other two sites are antialligned ( see fig . [ fig6 ] ( @xmath42)-(@xmath45 ) and [ fig7](@xmath46 ) ) . the resulting spin - charge order can be visualized as two interpenetrating lattices : an afm ordered honeycomb lattice and a nonmagnetic triangular lattice . this spontaneous separation into two sublattices is similar to what happens in the pinball - liquid phase . the crucial difference being the insulating character of the dos in the present case . for @xmath96 , the charge densities on the three sites are given by @xmath97 , and all three sites have finite magnetic moments . the resulting state with up - up - down structure have also been found as the ground state of the ising spin kondo lattice model . for @xmath98 the charge densities at three sites are given by @xmath99 , and the magnetic state again becomes non - collinear ( see inset in fig . [ fig6](@xmath45 ) ) . upon a further increase of @xmath81 the charge density as well as the magnetic moment on one of the sites keep reducing , giving rise to an effective honeycomb lattice for the electrons in the large @xmath81 limit a systematic search is carried out for non - collinear and non - coplanar magnetic ground states in the hubbard model . using a mean - field decoupling scheme that treats collinear , non - collinear and non - coplanar order on equal footing , we uncover a rich phase diagram for the hamiltonian on anisotropic triangular lattice . the most notable of the unconventional magnetic ground states are , ( i ) the four - sublattice flux order , ( ii ) the 3q non - coplanar state , and ( iii ) the six - sublattice non - collinear charge ordered state . these states have been found as ground states of the kondo lattice model on triangular lattice in previous studies @xcite . the effect of nearest neighbor coulomb interactions on these magnetic states is also investigated , and leads to some fascinating spin - charge ordered phases . the flux state is found to be compatible with a stripe - like charge ordered arrangement , whereas , the non - coplanar 3q state competes with charge ordering and becomes unstable upon increasing the coulomb interaction strength beyond a critical value . the spin - charge ordered state found at @xmath44 is a triangular lattice version of the state found in quarter - filled hubbard - holstein and kondo lattice models @xcite . the 120@xmath0state gives way to a sequence of spin - charge ordered states which lead to spontaneous decoupling of the triangular lattice into a honeycomb lattice and a triangular lattice . for a certain range of parameter values , one of the sublattice is magnetic while the other is nonmagnetic . this highly resembles the pinball - liquid state reported in theoretical studies of multiband extended hubbard model on triangular lattice @xcite . the phase diagrams presented here set the stage for further studies to analyze the stability of these phases . a number of states reported here possess a robust gap in the electronic spectrum , and therefore are likely to remain stable against higher order quantum effects . given the elementary status of the hubbard model and the recent interest in unconventional magnetic ordering , our results should be of general interest . in particular , the reported phases are of relevance to various triangular lattice systems , such as , the layered agnio@xmath3 @xcite , the nacoo@xmath3 @xcite and the organic charge - transfer salts @xcite . many of the phases presented in this work are similar to the spin - charge ordered phases reported experimentally in these materials . while material - specific models have been proposed to understand these phases , it is interesting to see that the most elementary model for magnetism  the single - band hubbard model  supports many of the unconventional spin - charge ordered states . the calculations were performed using the high performance computing facility at iiser mohali . k.p . acknowledges support via csir / ugc fellowship . s.k . acknowledges support from dst , india . 58ifxundefined [ 1 ] ifx#1 ifnum [ 1 ] # 1firstoftwo secondoftwo ifx [ 1 ] # 1firstoftwo secondoftwo `` `` # 1''''@noop [ 0]secondoftwosanitize@url [ 0 ]  + 12$12  & 12#1212_12%12@startlink[1]@endlink[0]@bib@innerbibempty link:\doibase 10.1103/physrev.82.403 [ * * , ( ) ] https://books.google.com/books?hl=en{&}lr={&}id=iynzfufnkbgc{&}pgis=1[__ ]  ( ,  ) p.   link:\doibase 10.1098/rspa.1963.0204 [ * * , ( ) ] link:\doibase 10.1103/revmodphys.56.755 [ * * ,   ( ) ] https://books.google.com/books?hl=en{&}lr={&}id=d-shcaaaqbaj{&}pgis=1[__ ] ( ,  )  p. link:\doibase 10.1103/physrevb.31.4403 [ * * , ( ) ] link:\doibase    10.1103/revmodphys.82.1539 [ * * ,   ( ) ] link:\doibase 10.1103/physrevb.72.094427 [ * * ,   ( ) ] link:\doibase 10.1103/physrevb.91.075104 [ * * ,   ( ) ] link:\doibase 10.1103/physrevlett.95.057205 [ * * ,   ( ) ] link:\doibase 10.1103/physrevlett.96.067601 [ * * ,   ( ) ] link:\doibase 10.1103/physrevlett.101.156402 [ * * ,   ( ) ] link:\doibase 10.1103/physrevlett.105.216405 [ * * ,   ( ) ] link:\doibase 10.1143/jpsj.79.083711 [ * * ,   ( ) ] link:\doibase 10.1103/physrevlett.108.096401 [ * * , ( ) ] link:\doibase 10.1103/physrevlett.109.237207 [ * * ,   ( ) ] link:\doibase 10.1103/physrevlett.108.257205 [ * * ,   ( ) ] link:\doibase 10.1103/physrevlett.109.166405 [ * * ,   ( ) ] link:\doibase    10.1103/physrevb.91.140403 [ * * ,   ( ) ] link:\doibase 10.1103/physrevlett.105.226403 [ * * ,   ( ) ] link:\doibase 10.1103/physrevb.74.193107 [ * * ,   ( ) ] link:\doibase    10.1143/jpsj.75.123704 [ * * , ( ) ] link:\doibase    10.1143/jpsj.78.014707 [ * * , ( ) ] link:\doibase 10.1103/physrevb.91.165139 [ * * ,   ( ) ] link:\doibase 10.1103/physrevb.77.094439 [ * * ,   ( ) ] link:\doibase    10.1103/physrevlett.99.157204 [ * * ,   ( ) ] link:\doibase 10.1134/s0021364014080116 [ * * , ( ) ] link:\doibase 10.1016/j.physb.2014.11.040 [ ( ) ,  10.1016/j.physb.2014.11.040 ] link:\doibase 10.1103/physrevlett.93.167601 [ * * , ( ) ] link:\doibase 10.1103/physrevlett.93.167003 [ * * , ( ) ] link:\doibase 10.1209/0295 - 5075/85/47006 [ * * ,   ( ) ] link:\doibase 10.3390/cryst2031155 [ * * ,   ( ) ] link:\doibase 10.1038/nphys2642 [ * * ,   ( ) ] http://arxiv.org/abs/1505.02290 [ ( ) ] ,   link:\doibase 10.1103/physrevb.51.12605 [ `` , ''  ] ( ) link:\doibase 10.1103/physrevb.60.6121 [ * * ,   ( ) ] link:\doibase 10.1103/physrevb.82.155102 [ * * ,   ( ) ] link:\doibase 10.1103/physrevlett.113.246405 [ * * , ( ) ] link:\doibase 10.1103/physrevlett.65.2462 [ * * ,   ( ) ] link:\doibase 10.1088/0953 - 8984/23/50/505601 [ * * ,   ( ) ] @noop * * ,   ( ) link:\doibase 10.1103/physrevb.62.13816 [ * * , ( ) ] link:\doibase 10.1016/j.physb.2005.01.233 [ * * ,   ( ) ] link:\doibase 10.1103/physrevb.81.064420 [ * * ,   ( ) ] link:\doibase 10.1103/physrevb.89.045115 [ * * ,   ( ) ] link:\doibase 10.1103/physrevb.39.9749 [ * * , ( ) ] link:\doibase 10.1103/physrevb.46.11110 [ * * , ( ) ] link:\doibase 10.1103/physrevb.81.094407 [ * * ,   ( ) ] link:\doibase 10.1103/physrevb.86.235118 [ * * ,   ( ) ] link:\doibase 10.1103/physrevlett.64.950 [ * * , ( ) ] link:\doibase 10.1209/0295 - 5075/15/6/011 [ * * ,   ( ) ] link:\doibase    10.1103/physrevb.91.245125 [ * * ,   ( ) ] link:\doibase    10.1103/physrevb.87.035143 [ * * ,   ( ) ] link:\doibase 10.1103/physrevlett.111.126403 [ * * , ( ) ] link:\doibase    10.1103/physrevb.77.214505 [ * * ,   ( ) ] link:\doibase 10.1143/jpsj.74.2901 [ * * ,   ( ) ] @noop * * ,   ( ) link:\doibase 10.1103/physrevlett.110.246401 [ * * , ( ) ]
Abstract:

Answer: motivated by the importance of non - collinear and non - coplanar magnetic phases in determining various electrical properties of magnetic materials , we investigate the phase diagrams of the extended hubbard model on anisotropic triangular lattice . we make use of a mean - field scheme that treats collinear , non - collinear and non - coplanar phases on equal footing . in addition to the ferromagnetic and 120@xmath0antiferromagnetic phases , we find the four - sublattice flux , the 3q non - coplanar and the non - collinear charge - ordered states to be stable at specific values of filling fraction @xmath1 . inter - site coulomb repulsion leads to intriguing spin - charge ordered phases . most notable of these are the collinear and non - collinear magnetic states at @xmath2 , which occur together with a pinball - liquid - like charge order . our results demonstrate that the elementary single - orbital extended hubbard model on a triangular lattice hosts unconventional spin - charge ordered phases , which have been reported in more complex and material - specific electronic hamiltonians relevant to layered triangular systems .


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 statistics and general structure of passive scalars in turbulent flows differ significantly from those of the turbulent velocity field . in particular , conclusive experimental ( e.g. stewart 1969 , sreenivasan , antonia & britz 1979 , mestayer 1982 , and mydlarski & warhaft 1998a ) and numerical ( holzer & siggia 1994 ) evidence shows that structure functions and the derivative skewness of the scalar field do not follow predictions from isotropy at inertial and dissipative scales , in the presence of a mean scalar gradient . in particular , the deviations are thought to be related to `` ramp and cliff structures '' and to imply a direct effect of large - scale structures on small - scale structures . the data relevant to this question has been reviewed by sreenivasan ( 1991 ) and more recently by warhaft ( 2000 ) . moreover , scalar spectra display a distinctly less universal structure than velocity spectra . this is manifested both in terms of spectral exponents , as well as in terms of the dimensionless spectral coefficient @xmath0 . the latter varies from values near 0.4 in many experiments in the atmospheric surface layer and grid turbulence ( see sreenivasan 1991 , sreenivasan 1996 and warhaft 2000 ) to @xmath1 for other atmospheric measurements ( antonia , ould - rouis , anselmet & zhu 1997 , and antonia , xu & zhou 1999 ) . from a fundamental point of view , these observations challenge the kolmogorov cascade phenomenology for the transfer of scalar - variance from large to small scales . in the classical phenomenology , the multiplicity of separate `` eddy breakdown '' events is assumed to gradually uncouple the small from the large scales , allowing the former to tend to a universal and isotropic structure more or less independent of the large scales . on the contrary , the scalar - field measurements suggest a direct linkage among largest and smallest scales . the coupling among scales is an important ingredient in large - eddy simulations ( les ) . in les , the turbulent fields ( velocity and scalar ) are decomposed into large and small ( subgrid - scale , sgs ) scale contributions by means of a spatial low - pass filter of characteristic width @xmath2 . the resulting equations which can be numerically discretized with mesh - spacing of the order of @xmath2 require closure of the unresolved momentum fluxes ( sgs stress tensor , @xmath3 , where a _ tilde _ represents filtering at scale @xmath2 ) and scalar fluxes ( e.g. the sgs heat flux , @xmath4 where @xmath5 is the passive scalar field ) . the promise of les is often predicated upon small - scale universality and isotropy , and the absense of a strong coupling across disparate length scales . hence , the observed anisotropy of the scalar field seems to pose a challenge to the very foundation of les . while these deviations from classical phenomenology are now quite well established , little is known about their impact on the closure problem for les . the present work quantifies the implications of small - scale scalar anisotropy on quantities that describe subgrid - scale physics and directly affect modeling for les . as reviewed in meneveau & katz ( 2000 ) , the most important statistical property of the fluxes @xmath6 and @xmath7 is how they affect the mean kinetic energy and scalar - variance budgets of the resolved fields . specifically , their dominant effect is through the kinetic energy and scalar - variance dissipations that arise from interactions between subgrid and resolved scales . therefore , in the present study , we mainly focus on the so - called sgs kinetic energy dissipation @xmath8 ( piomelli , 1991 ) and scalar - variance dissipation @xmath9 ( port - agel , meneveau & parlange 1998 ) . here @xmath10 and @xmath11 are the resolved strain - rate tensor and scalar gradient , respectively . the sgs dissipation rates represent the flux ( cascade ) of kinetic energy or scalar - variance from resolved towards subgrid scales ( when positive ) . when @xmath2 pertains to the inertial range , and when the flow is in equilibrium , one expects the mean sgs dissipation to equal the molecular dissipation rate . deviations from isotropy in the context of sgs dissipation can be probed by measuring the isotropy level of the tensors @xmath12 and @xmath13 , as a function of scale . the main question to be addressed in this work is whether the approach to isotropy ( if it exists ) is the same for kinetic energy and scalar - variance dissipation tensors . another goal is to test the ability of two popular model classes ( eddy diffusivity and nonlinear models , see meneveau & katz 2000 for a review ) to reproduce the observations . in order to span a sizeable range of inertial - range filter scales , a sufficiently high reynolds number must be considered . hence , this study is based on experimental data ( as opposed to dns which is limited to small reynolds numbers ) . the study is performed in a canonical shear flow , the heated cylinder wake . experiments were performed in the return type corrsin wind tunnel ( comte - bellot & corrsin 1966 ) . a heated smooth cylinder of diameter @xmath14 = 4.83 cm was located horizontally at the centerline of the test section . the measurement location in the streamwise direction @xmath15 was fixed at @xmath16 . to obtain the filtered and sgs quantities , an array of four custom - made miniature probes was used . each probe was composed of one x - type hot - wire and one i - type cold - wire for the velocities in the @xmath17 plane and the temperatures , respectively . here @xmath18 is the `` cross - wake '' direction , i.e. perpendicular to @xmath19 and to the cylinder axis . the separation distance @xmath20 between the probes in the cross - wake direction @xmath18 could be adjusted manually between 5 and 20 mm . three configurations ( @xmath21 , with @xmath22 , 10 and 20 mm ) were used in the present study . a 2.5 @xmath23 m platinum - coated tungsten wire which had been copper - plated was soldered on to the x - wire prong ends and etched , yielding an active length - to - diameter ratio of about 200 . the wire spacing between the hot - wires was 0.5 mm . a 0.625 @xmath23 m silver - coated pure platinum wire for a cold - wire sensor was soldered on the i - wire prong ends and subsequently etched . to minimize the low frequency amplitude attenuation , the active length - to - diameter ratio was about 1000 as suggested by @xcite . the separation distance between the cold - wire and its nearest hot - wire was 0.9 mm so that the thermal effect from the hot - wire on the cold - wire was negligible . the signals were low - pass filtered at a frequency of 20 khz and sampled at @xmath24 khz . sampling time was 60 seconds , so the total number of data points per channel for each measurement location was @xmath25 . the array was traversed across the wake , and data were recorded at 17 discrete cross - wake locations from the centerline to the wake edge at increments of 14.4 mm .    at the measurement location of @xmath16 , the mean centerline velocity @xmath26 was @xmath27 m / s , the defect velocity @xmath28 was @xmath29 m / s , the defect temperature @xmath30 was @xmath31c , and the half - width of the wake was @xmath32 m. to get the spatial quantities along the streamwise direction from the temporal data , taylor s hypothesis was invoked . the turbulence intensity of the streamwise velocity at the centerline was about @xmath33 . the molecular kinetic energy dissipation at the centerline @xmath34 was @xmath35 m@xmath36s@xmath37 , and the molecular scalar - variance dissipation at the centerline @xmath38 was @xmath39c@xmath40s@xmath41 . the latter two variables were obtained from ( corrected ) third - order structure functions as in cerutti ( 2000 ) and @xcite . it follows that the kolmogorov length scale @xmath42 was @xmath43 mm , the taylor micro - scale @xmath44 was @xmath45 mm , and the reynolds number based on taylor micro - scale @xmath46 was @xmath47 . the longitudinal integral scale obtained by integrating up to the first zero crossing of the @xmath48 correlation function was @xmath49 m. profiles of mean velocity , rms velocities , reynolds shear stress , rms temperature and heat flux distributions across the wake agreed quite well with results in the literature ( e.g. , matsumura & antonia 1993 and kiya & matsumura 1988 ) .    in the present study , to separate between large and small scales , the box filter is applied to the streamwise and cross - wake directions , and the trapezoidal rule is used for the spatial integrations . the filtering process consists in a discrete approximation to a two - dimensional box filter . in the @xmath18 direction , a four point discretization is used for evaluating the sgs fluxes while a three - point approximation is used for the filtered derivatives . filtered velocity and scalar gradients in the @xmath18 direction are evaluated using first - order finite differences over a distance @xmath20 . in the streamwise direction , the box filter is approximated using @xmath50 sampling points and the @xmath19 derivatives are evaluated using finite differences over a distance @xmath20 . the filtering and error analysis is documented in cerutti & meneveau ( 2000 ) and cerutti _ et al . _ ( 2000 ) .    figure 1(a ) shows a comparison between the longitudinal spectrum @xmath51 of the @xmath48 component and the longitudinal spectrum @xmath52 of the @xmath53 component multiplied by 3/4 at the centerline . here , @xmath54 is the longitudinal wave number . the three vertical dashed lines correspond to the filter sizes of @xmath55125 ( 10 mm ) , 250 ( 20 mm ) , 500 ( 40 mm ) . all the filter sizes are in the inertial range . since the noise peak in the longitudinal spectrum of @xmath48 is in the far dissipation region , quite removed from any of the filter frequencies and scales of interest in this study , no effort is made to remove the noise by additional filtering ( various attempts such as notch filtering showed no effect on the results ) . it can be clearly observed that @xmath56 as required by isotropy in the inertial range , over about one decade of wave numbers . the peak in @xmath52 at @xmath57 is due to the periodic von krmn vortex street behind the cylinder . the frequency of the vortices is 78.1 hz , and this gives strouhal number @xmath58 . the kolmogorov constant @xmath59 in @xmath60 is obtained as @xmath61 in the insert . this value is quite close to the standard value of 1.6 ( see review by sreenivasan 1995 ) and a similar value of 1.7 was observed by oneil & meneveau ( 1997 ) . figure 1(b ) shows the longitudinal spectrum of the temperature @xmath62 . as seen in the premultiplied spectrum shown in the insert , there is a fairly clear inertial range with a @xmath63 slope . this slope is steeper than those found in the round jet by tong & warhaft ( 1995 ) and miller & dimotakis ( 1996 ) , and in several other shear flows reviewed in sreenivasan ( 1996 ) . it is closer to results quoted in antonia & pearson ( 1997 ) , who report a scaling exponent of 0.65 - 0.66 for the @xmath64-order temperature structure function in the heated cylinder wake at @xmath65 , or to the spectra for grid turbulence of mydlarski & warhaft ( 1998b ) . the coefficient @xmath0 in @xmath66 deduced from this spectrum is about 1.4 . it is significantly higher than results quoted in sreenivasan ( 1996 ) for high reynolds numbers , but is within the range of results from various shear flow measurements reported in antonia 1999 . as discussed in sreenivasan ( 1996 ) , the dissipation measurements are based on assuming scalar isotropy and are thus subject to considerable uncertainty . on the other hand , at the centerline we find good isotropy ( see section 3 ) which would seem to support our current estimates of the dissipation and @xmath0 . while the universality of @xmath0 and spectral exponent is not the main subject of this paper , the scatter of results certainly supports the view that the scalar spectrum for shear flows at moderate reynolds numbers ( e.g. @xmath67 ) depends upon details of the generation of the flow ( sreenivasan , 1996 ) . specifically , it seems that the dependence of the spectral exponent and prefactor on @xmath68 is not universal . the longitudinal spectrum of the filtered temperature is shown as the dashed curve in figure 1(b ) . the lobes at scales below the filter size are due to the streamwise box filter used in the present study . spectra of filtered velocity ( not shown ) have similar shape . in discussing isotropy of various tensors , we distinguish second - rank and fourth - rank tensors . examples of second - rank tensors are the mean sgs stress @xmath69 , the filtered scalar - gradient @xmath70 , and the sgs scalar - variance dissipation @xmath71 . their isotropic form can be written as @xmath72 , @xmath73 , and @xmath74 respectively . on the other hand , the strain - rate - product tensor @xmath75 and sgs dissipation tensor @xmath76 are fourth - rank tensors . their isotropic form can be written as @xmath77,\ ] ] and a similar expression holds for @xmath75 by replacing @xmath6 with @xmath78 . to derive the above expression , we have used the tensor symmetry ( in @xmath79 or @xmath80 ) and the divergence - free condition ( @xmath81 ) . it follows that @xmath82 and @xmath83 . figure 2 shows isotropy ratios of the filtered gradient fields and of the mean sgs stresses , for @xmath84 . in isotropic conditions , all ratios have to be on the unit line except @xmath85 which should tend to zero with decreasing filter scale since the fraction of mean shear stress carried by the sgs scales is expected to vanish when @xmath86 . as can be seen , both the sgs stress and strain - rate fields are quite isotropic across the heated wake flow at this scale . for the larger filter sizes of @xmath87 and @xmath88 ( not shown ) , the trends are similar to figure 2 but slightly less isotropic . next , results are given for the sgs dissipation tensors across the wake . we define the following ` isotropy ratios ' : @xmath89    figures 3(a ) and ( b ) show the spatial distributions of the sgs kinetic energy dissipation isotropy ratios @xmath90 and the sgs scalar - variance dissipation isotropy ratios @xmath91 , respectively , across the wake flow with the different filter sizes . all moments in this study are statistically well converged . for example , running averages such as @xmath92 did not differ from their global average by more than 6% over the last 80% of data ( i.e. for @xmath93 $ ] where @xmath94 ) . for different filter sizes considered , the sgs kinetic energy and scalar - variance dissipations are isotropic only towards the centerline of the wake . towards the edges of the wake , large anisotropies persist ( even though there the mean shear and mean temperature gradients vanish also ) . the isotropy ratio of @xmath90 at a fixed @xmath95 in figure 3(a ) clearly increases with the filter size . however , the isotropy ratio of the scalar - variance dissipation @xmath91 shows almost no variation with filter size for @xmath96 , where the scalar gradient in the cross - wake direction , @xmath97 , is large . the increased anisotropy in the outer regions of the wake ( see figures 3(a ) and ( b ) ) is another striking result . the anisotropy may be associated with the highly intermittent character of the flow there . one may wonder whether the anisotropy arises from a superposition of distinct behaviors in the turbulent and the non - turbulent ( outer ) regions . using conditional averaging , oneil & meneveau ( 1997 ) already showed that the conditional averages of the sgs dissipations in the non - turbulent regions were negligible compared to those in the turbulent regions . the implication was that the global averages of dissipation could be explained entirely by their conditional mean value inside the turbulent part : @xmath98 , where @xmath99 is the intermittency function ( i.e. the fraction of time the signal is turbulent at any given @xmath18 ) , and the subscript `` @xmath100 '' stands for averaging conditioned on ` turbulence ' ( see oneil & meneveau 1997 for details ) . when replacing the conditional averages multiplied by @xmath101 in the expression for isotropy ratios , @xmath101 cancels from both numerator and denominator . this behavior suggests that the anisotropy ratios shown in figure 3 are equal to those inside the turbulent regions alone , and that their rise in the outer parts of the wake can not be explained by contributions coming from the non - turbulent parts . however , it is still possible that non - trivial contributions to the scalar - variance dissipations could originate at the interface separating turbulent from non - turbulent regions .    to highlight the variations with filter scale more directly , in figure 4 we plot the isotropy ratios as a function of scale , for the transverse location @xmath102 close to the peak mean temperature gradient . as is evident , the isotropy ratios of the kinetic energy dissipation decrease towards unity as the filter size decreases , whereas the isotropy ratio of the scalar - variance dissipation remains unchanged near @xmath103 as the filter size is decreased . as seen in figure 3 , near the centerline , where the mean shear and scalar gradient vanish , the sgs isotropy ratios are all near unity . therefore , we conclude that in terms of the most important of the resolved - sgs scale interactions ( the mean sgs dissipation ) , the scalar field maintains strong anisotropy in the presence of a mean scalar gradient . conversely , the velocity field has the trends expected from approach to isotropy at small scales . next , we quantify the isotropy level of model predictions , when the data are analyzed in an ` a - priori ' sense , i.e. by replacing @xmath6 and @xmath7 above by model expressions . first , the standard eddy - diffusion model is considered , i.e. , @xmath104 where @xmath105 is the modulus of the resolved strain rate , @xmath106 is the smagorinsky coefficient , and @xmath107 is the sgs prandtl number . consequently , and independent of the model coefficients , the isotropy ratios from the eddy - diffusion model are : @xmath108 in computing the filtered strain - rate magnitude from the data , the following approximation is used ( this is a 2d extension of the 1d approach of oneil & meneveau 1997 , and was also used in liu , katz & meneveau 1999 ) : @xmath109^\frac{1}{2}. $ ] this approximation is not expected to affect the accuracy of our measured isotropy ratios significantly . the reason is that @xmath110 and its 2d approximation are scalars which multiply equally all terms of the squared velocity and scalar gradients and is further supported by the isotropic behavior of second - order moments of @xmath111 shown in figure 2 . the isotropy ratios of the sgs dissipations from the eddy - diffusion model are shown in figure 5(a ) . results are near unity almost independently of scale , including the passive scalar dissipation . profiles across the wake ( not shown ) also are near unity . this result is consistent with the observed isotropy of the square gradient tensors shown in figure 2 ( the only difference here is the additional @xmath112 factor ) . hence , using the eddy - diffusion model one would ( incorrectly ) predict sgs isotropy since the resolved second - order moments are isotropic . the second modeling class to be considered here is the ` nonlinear model ' or ` tensor eddy - diffusion model ' ( leonard 1974 , clark , ferziger & reynolds 1979 , liu , meneveau & katz 1994 , leonard 1997 , borue & orszag 1998 and meneveau & katz 2000 ) . the nonlinear model reads as follows : @xmath113 where @xmath114 and @xmath115 are the nonlinear model coefficients . therefore , the corresponding isotropy ratios from the nonlinear model are defined as follows : @xmath116 where @xmath117 varies from 1 to 2 in the analysis of our 2d data .    the isotropy ratios of the sgs dissipations from the nonlinear model are shown in figure 5(b ) , as a function of filter scale . as is apparent in comparing with figure 4 , the main features of anisotropy and filter - size dependence are correctly reproduced qualitatively . quantitatively , the levels of anisotropy are overestimated . for instance , the level of anisotropy for the modeled sgs scalar - variance dissipation appears to stay near @xmath118 as opposed to @xmath103 for the real sgs scalar dissipation . in studying passive scalar statistics in a turbulent shear flow with a mean temperature gradient we focus on statistics of interest to subgrid modeling and large - eddy simulation . in order to obtain the filtered and subgrid velocities and temperatures , a probe array composed of four x - wire and four cold - wire sensors is used and two - dimensional box filtering in the streamwise and cross - wake directions is applied to the data . the isotropy ratios of the sgs kinetic energy and scalar - variance dissipations are investigated as function of position in the flow and as a function of filter scale . both dissipations are isotropic independent of filter size near the centerline where there is no mean shear or scalar gradient . however , at locations with high gradient ( and also further towards the outer wake regions ) , we find that the scalar - variance dissipation remains highly anisotropic , independent of filter size . conversely , the kinetic energy dissipation tensor approaches isotropy as the filter size is decreased . a mechanistic explanation of the observed trends in terms of possible orientations of ramp and cliff structures is not evident to us at this time . the persistence of scalar anisotropy even at small scales is consistent with prior results for structure functions and gradient statistics of unfiltered turbulence ( mestayer 1982 , sreenivasan 1991 , warhaft 2000 ) . present results quantify the impact of this anisotropy on the interactions among large and small scales in the context of sgs modeling and les . we find that the predictions of eddy - diffusion models are much more isotropic than the real phenomenon . this may at first glance seem obvious since eddy - diffusion models are often invoked under the banner of isotropy . however , the present result does not arise from an explicit isotropy assumption built into the model , but because the resolved gradients have second - order statistics that are isotropic . for instance , if the scalar gradient tensor @xmath119 had been found to be anisotropic , it would have implied anisotropic behavior of the eddy - diffusion model s predictions . instead , the isotropy that exists in the filtered scalar gradients is incorrectly applied to model the anisotropic statistics of the sgs heat flux . the main problem for the eddy - diffusion model appears to be that it uses second - order statistics to model third - order statistics . conversely , the anisotropy exists in the third - order moments that arise from the velocity and scalar product ( @xmath120 ) multiplied by the scalar gradient ( @xmath121 ) , but is not discernible in the second - order statistics of @xmath121 alone ( even when modulated by the strain - rate magnitude ) . the anisotropy is clearly discernible , however , in the third - order moments consisting of the filtered velocity gradients times scalar gradients squared that arise in the expression for modeled sgs dissipation of scalar - variance using the nonlinear model . these expressions are able to reproduce the detailed phase relationships among the velocity and scalar field that govern the sgs dissipation ( cascade ) of scalar variance . we remark that the overprediction of anisotropy by the nonlinear model ( and its underprediction by the eddy - diffusion model ) is reminiscent of the opposing trends of these two models documented in liu ( 1999 ) for rapidly strained turbulence in cold - flow . the opposing trends suggest that a linear combination of the two models , i.e. the ` mixed model ' , can be tuned to reproduce the correct amount of anisotropy ( for a discussion of the application of mixed models in les , see meneveau & katz 2000 ) . however , note that the data for @xmath90 and @xmath122 at @xmath87 and 125 do not provide clear justification for choosing the mixed over the non - linear model . finally , it is stressed that current results are obtained in a single flow for a single moderate reynolds number . even if the present evidence in figure 4 of an essentially scale - independent anisotropy for the scalar dissipation appears to be quite strong , the results could change in another flow , or at higher reynolds numbers ( we recall that according to sreenivasan 1996 , universal behavior for the scalar requires @xmath68 above 1000 or so ) . these considerations serve as motivation for further work in this area . we thank professors z. warhaft and l. mydlarski for useful comments and suggestions about the cold - wire calibration . this work was financially supported by the national science foundation ( grant cts-9803385 ) .
Abstract:

Answer: the effects of passive scalar anisotropy on subgrid - scale ( sgs ) physics and modeling for large - eddy simulations are studied experimentally . measurements are performed across a moderate reynolds number wake flow generated by a heated cylinder , using an array of four x - wire and four cold - wire probes . by varying the separation distance among probes in the array , we obtain filtered and subgrid quantities at three different filter sizes . we compute several terms that comprise the subgrid dissipation tensor of kinetic energy and scalar - variance and test for isotropic behavior , as a function of filter scale . we find that whereas the kinetic energy dissipation tensor tends towards isotropy at small scales , the sgs scalar - variance dissipation remains anisotropic independent of filter scale . the eddy - diffusion model predicts isotropic behavior , whereas the nonlinear ( or tensor eddy diffusivity ) model reproduces the correct trends , but overestimates the level of scalar dissipation anisotropy . these results provide some support for so - called mixed models but raise new questions about the causes of the observed anisotropy .