A simplifie BBGKY hierarchy for correlate fermionic systems from a Stochastic Mean-Fiel approach arxiv:1511.03815v1 [nucl-th 1 Nov 015 Denis Lacroix, 1, Yusue Tanimura, 1 Sair Ayi, an Bulent Yilmaz 3 1 Institut e Physique Nucléaire, INP3-CNRS, Université Paris-Su, F-9106 Orsay Ceex, France Physics Department, Tennessee Technological University, Cooeville, TN 38505, USA 3 Physics Department, Faculty of Sciences, Anara University, 06100, Anara, Turey Date: August, 018 The stochastic mean-fiel SMF approach allows to treat correlations beyon mean-fiel using a set of inepenent mean-fiel trajectories with appropriate choice of fluctuating initial conitions. We show here, that this approach is equivalent to a simplifie version of the Bogolyubov-Born- Green-Kirwoo-Yvon BBGKY hierarchy between one-, two-,..., N-boy egrees of freeom. In this simplifie version, one-boy egrees of freeom are couple to fluctuations to all orers while retaining only specific terms of the general BBGKY hierarchy. The use of the simplifie BBGKY is illustrate with the Lipin-Meshov-Glic LMG moel. We show that a truncate version of this hierarchy can be useful, as an alternative to the SMF, especially in the wea coupling regime to get physical insight in the effect beyon mean-fiel. In particular, it leas to approximate analytical expressions for the quantum fluctuations both in the wea an strong coupling regime. In the strong coupling regime, it can only be use for short time evolution. In that case, it gives information on the evolution time-scale close to a sale point associate to a quantum phase-transition. For long time evolution an strong coupling, we observe that the simplifie BBGKY hierarchy cannot be truncate an only the full SMF with initial sampling leas to reasonable results. PACS numbers:.75.i, 1.60.Jz,7.90.b I. INTRODUCTION The BBGKY hierarchy [1 3 is an exact reformulation of the problem of interacting fermions where oneboy egrees of freeom DOFs are couple to two-boy DOFs that are themselves couple to three-boy DOFs an so one an so forth. One of the avantages of this hierarchy is that it illustrates how correlations can affect the evolution of the one-boy ensity matrix. For this reason, the set of equations between one-boy evolution an many-boy correlations are often use as a starting point to evelop approximations beyon mean-fiel see for instance [. However, ue to the increasing complexity occurring when complex correlations are inclue, a truncation scheme of the hierarchy is necessary. The proposal of an appropriate truncation of the BBGKY hierarchy has been the subject of extensive wor in the past [5 7 an leas to the so-calle Time- Depenent Density Matrix TDDM approach to interacting systems where not only the one-boy ensity is followe in time but also eventually the two-boy an eventually three-boy ensity matrices. The flexibility of choosing the proper equations of motion is still the subject of intensive wor [8. In particular the absence of a efinite strategy for truncating the BBGKY can lea to uncontrolle results with varying quality as illustrate in Ref. [9 It has been shown recently that the stochastic meanfiel approach can provie rather accurate escription of lacroix@ipno.inp3.fr many-boy effects beyon the Time-Depenent Hartree- Foc TDHF or Time-Depenent Hartree-Foc Bogolyubov TDHFB approach in several test cases [10 1 where the approximate treatment can be confronte to the exact solution. In this approach, the system is escribe by a set of initial conitions with fluctuating oneboy ensity, followe by eterministic TDHF or TD- HFB trajectories [13. Each trajectory evolves through its own self-consistent mean-fiel for a given initial conition an is inepenent from the others. Among the interesting aspects of the SMF approach one can mention that beyon mean-fiel effects are incorporate although only mean-fiel type evolution is neee. One of the surprising results turns out to be that the SMF technique can provie a better approximation compare to an approach where two-boy egrees of freeom are explicitly introuce in the escription [1. One of the main objectives of the present wor is to clarify how a theory that only involves mean-fiel evolution can incorporate properly many-boy effects. Below, we mae explicit connection between the SMF approach an the BBGKY hierarchy. We show that the SMF theory is equivalent to solve a set of couple equations between one-boy DOFs an higher-orer correlations without truncation. This set of equations can be seen as a simplifie version of the stanar BBGKY hierarchy. We finally explore if this simplifie BBGKY hierarchy can be use either as a numerical tool or to get physical insight in correlate systems if a truncation is assume.
II. THE SMF APPROACH AND ITS CONNECTION TO A NON-TRUNCATED BBGKY HIERARCHY A. Basic aspects of the stochastic mean-fiel approach In the SMF theory, the N-boy problem is replace by a set of eterministic time-epenent mean-fiel trajectories [13, i ρn t = [ hρ n, ρ n, 1 where n labels a given trajectory, ρ n is the one-boy ensity along this trajectory an hρ n is the associate self-consistent mean-fiel. For a system interacting through a two-boy interaction v 1, h is given by: h[ρ = t Tr {ṽ 1 ρ }, where t is the inetic term an ṽ enotes antisymmetric effective interaction. Here, we use the convention of Refs. [, 1, i.e. the inex i in ρ i means that the ensity acts on the i th particle, similarly Tr i. means that the partial trace is mae on the i th particle. Starting from the above expression an using this convention, we recover the stanar mean-fiel expression h ij [ρ = t ij m ṽ i,jm ρ m. In the SMF approach, each trajectory is eterministic an fluctuations stem only from the initial conitions ρ n t 0. The statistical properties of the ensity matrices are usually chosen in such a way that they reprouce at least in an approximate way the initial quantum phase-space of the problem. This is usually one by imposing specific properties of ifferent moments of the initial ensity fluctuations δρ n t 0 = ρ n t 0 ρ n t 0. X n enotes here the statistical average over the initial conitions. Most often, the moments are chosen to reprouce quantum fluctuations of the initial many-boy states. In its original form [13, the SMF theory has been formulate assuming i that the initial state is either a pure Slater eterminant or a statistical ensemble of inepenent particles. In both cases, the one-boy ensity is associate with occupation numbers enote by n i in the canonical basis. ii The initial quantum phase-space has been approximate by a Gaussian statistical ensemble. Then, initial phase-space sampling is completely specifie by the first an secon moments of ρ n t 0 that are conveniently chosen as: ρ n ij t 0 = δ ij n i, δρ n ij t 0δρ n l t 0 = 1 δ ilδ j [n i 1 n j n i 1 n j. Starting from this original formulation, it has been realize that the SMF approach can be extene to superflui systems by replacing the N-boy problem by a set of Time-Depenent Hartree-Foc Bogolyubov TD- HFB with initial fluctuations on both the normal an anomalous ensity matrices [10. In aition, we have shown recently that an improve treatment of the initial phase-space [15 without the Gaussian approximation can further ameliorate the many-boy ynamics beyon the mean-fiel an allows to escribe initially correlate many-boy states. More surprisingly, a recent stuy on the Hubbar moel [1 has shown that the SMF approach can lea to a better than state of the art Green function techniques incluing explicitly two-boy effects [16. This actually might appear as a surprising result since in the SMF approach only mean-fiel trajectories are require. It is one of the goals of the present article to show how correlations are incorporate in our theory. In particular, we show that the SMF approach is equivalent to a simplifie BBGKY hierarchy where -boy, 3-boy,..., N-boy correlations are approximately propagate in time. B. Many-boy evolution in SMF In the SMF theory, the quantum expectation value of any -boy operator is replace by the a statistical average over the ifferent trajectories. Consiering a set of one-boy operators, enote by A1,, A, their expectation values within SMF are given by A1 n A n = α iβ i A α1β 1 1 A α β ρ n β 1α 1 ρ n β α. Therefore, the nowlege of any many-boy observable is equivalent to the nowlege of the time evolution of the set of moments M 1, M 1,, M 1 efine through: M 1 = ρ n β 1α 1 ρ n β α. We show in appenix A, starting from the equation 1 that the ifferent moments evolve accoring to the following couple equations: [ i t M 1,, = t, M 1 Tr 1 [ṽ α1, M 1 1. 3 These equations are formally very close to the BBGKY hierarchy of ensity matrices in many-boy systems [, 1 except that here, the many-boy ensities are replace by ifferent moments of the ensity ρ n. Another important ifference is that part of the fermionic aspects are lost in this escription. Inee, the moments M 1,, here are symmetric with respect to the exchange of two inices. Therefore, part of the quantum correlations inuce by Fermionic statistic are lost in the theory.
3 This is one of the ifferences with the TDDM approach where many-boy ensities truly associate to fermions are solve in time. One can equivalently rewrite the above set of equations as couple equations between the average one-boy ensity enote for simplicity below as ρ an the centere moments C 1 efine through: C 1 = δρ n 1 δρ n. The resulting couple equation are see appenix A: i t ρt = [hρt, ρt Tr [ṽ 1, C 1, 5 together with for : i t C 1 = t α, C 1 α Tr 1 [ṽ α1, C 1 ρ 1 Tr 1 [ṽα1, C 1 α 1α1 1 ρ α Tr 1 [ṽα1, C 1 α 1α1 C α1 Tr 1 [ṽα1, C 1 1. 6 Equation 5 is similar to the first equation of the BBGKY hierarchy where usually C 1 stans for the twoboy correlation matrix. This equation clearly points out that effects beyon the stanar mean-fiel are accounte for in the SMF approach. A more precise iscussion on what many-boy effects are inclue is given below. C. Beyon mean-fiel effects in SMF theory The fact that the SMF approach is equivalent to solve an unrestricte set of couple equations between the oneboy ensity an higher-orer moments is a clear avantage of this technique. As a irect conclusion, it inclues not only two-boy effects but also higher orer correlation effects. This might explain why, in previous applications, it has given better results compare to the calculation using BBGKY hierarchy truncate at secon orer. Still, because of the replacement of quantum average by classical average, it is not anticipate that all manyboy effects are properly inclue with the SMF framewor. Let us write explicitly the equation on C 1. Starting from Eq. 6, after some straightforwar manipulation, we obtain: i t C 1 = [h 1 [ρ h [ρ, C 1 Tr 3 [ṽ 13 ṽ 3, C 13 ρ C 3 ρ 1 Tr 3 [ṽ 13 ṽ 3, C 13, 7 where h i is the mean-fiel Hamiltonian acting on the particle i. This equation correspons to a simplifie version of the secon BBGKY equation for correlation see Eq. 55-58 of Ref. [1. On the negative sie, we first observe that the term B 1 Eq. 57 of [1 that is responsible for both irect in-meium collisions an pairing effects is missing. Regaring in-meium collisions, it is anticipate that the SMF approach is only vali at low internal excitation of the system where the in-meium collisions are strongly hinere ue to the Pauli exclusion principle. It is worth mentioning that the inclusion of in-meium collisions can be mae eventually by consiering a Langevin process where fluctuations are introuce continuously in time [17. Regaring the pairing term, although the complete proof is out of the scope of the present wor, it is anticipate that this term is partially inclue when the superflui version of SMF is use [10. On the positive sie of SMF we can remar: The first line of Eq. 7 shows that the SMF approach properly propagate the initial fluctuations in the average self-consistent mean-fiel. The secon line of Eq. 7 shows that the term P 1 usually appearing in the secon BBGKY equation Eq. 58 of [1, is approximately accounte for. We note in particular, that the Pauli principle is not fully respecte in the SMF approach. This is inee not surprising ue to the simple classical assumption mae on the average moments. The last line of Eq. 7 shows that a great avantage of the resulting equation is that 3-boy an higher correlation effects are inclue. In orer to better unerstan what type of correlations are retaine in Eq. 7, it is interesting to note that the average ensity evolution can be rewritten as i t ρt = [hρt, ρt [ δh[ρ n, δρ n. 8 This alternative form eviences that the effect beyon mean-fiel stems from the correlation between the ensity fluctuation an the resulting mean-fiel fluctuation. Such correlations an their effects on one-boy ynamics have been stuie in Ref. [18 for small amplitue vibrations aroun the groun state of atomic nuclei. In that case, they lea to a coupling of the collective motion to surface vibrations the so-calle particle-phonon coupling. This coupling is nown as a ominant source of fragmentation of the nuclear collective response at low temperature [.
D. Discussion The recent success of the SMF approach an more particularly the goo agreement between the approximate treatment an exact solutions in some moel cases [10 1 inicate that this approach retains some of the most important correlation effects an properly incorporate their effects on one-boy egrees of freeom. The equivalent formulation given above in terms of a simplifie BBGKY hierarchy is helpful to clearly pinpoint the retaine terms. As a sie prouct of the present wor, we note that the present evelopment leas to a new approximate hierarchy of equations of motion that is much less involve than the stanar BBGKY one. In particular, it is commonly accepte that the main ifficulty in manyboy theories base on the truncation of the BBGKY hierarchy is the absence of a systematic prescription to truncate the couple equations at a given orer [9. The couple equations erive above can be seen as an alternative way to incorporate many-boy effects. Let us maes few remars on the avantages of the new equations of motion: i first, the equation of motion are much simpler than the stanar BBGKY hierarchy, as we will see below, in some cases, thir-orer or even fourth orer equation can be obtaine without ifficulty. ii When possible, it is clear that it is preferable to irectly perform the full SMF theory, i.e. sample a set of mean-fiel trajectories. However, in some cases, this might become rather cumbersome ue to number of trajectories require to get small statistical error-bars. Then, the corresponing BBGKY hierarchy can be an alternative approximate metho to solve SMF. iii As we pointe out in Ref. [15, a proper account for the initial quantum phase-space might require to go beyon the Gaussian assumption of the initial sampling. In simple moels, it is possible to irectly get the initial phasespace with some semi-classical techniques as in [15. In more complex systems, this seems more complicate an the sampling beyon the Gaussian approximation might not be possible. Alternatively, one can use the simplifie BBGKY hierarchy imposing that ifferent moments at initial time exactly equals the initial quantum moments. iv Last, as we will see below, the simplifie form of the BBGKY equation can lea to interesting physical information. This is illustrate below, where approximate analytical expressions are obtaine for the evolution incluing beyon mean-fiel effects. III. APPLICATION The goal of the present section is to illustrate the interesting aspects associate to the hierarchy of couple equations motivate by the SMF approach. For this purpose, we too as an example the Lipin-Meshov- Glic LMG Moel [0 3 that has alreay been use to benchmar the SMF theory using phase-space initial sampling both in the Gaussian [11 an non-gaussian [15 assumptions. In the two previous applications, a set of mean-fiel trajectories have been explicitly followe in time using the equations of motion: t jn x t jn y t jn z = j n y = j n x χj y n j z n, χj x n j z n, = χj y n j x n ; 9 with fluctuating initial conitions. Here, j x,y,z n are the three reuce quasi-spin components an χ is a constant that measures the strength of the two-boy interaction. All technical etails associate to the LMG moel, its SMF solution, as well as the possibility to get an exact solution have been extensively iscusse in our previous wor [11, 15. To get the close set of equations, we follow the general strategy epicte above an separate the average quasispin from its fluctuation as j n i = j i δj n i. Starting from Eq. 9, we obtain t j x = 1 χj z j y χ yz, t j y = 1 χj z j x χ xz, t j z = χj x j y χ xy, where ij δjn i j are the average fluctuations of the reuce quasi-spin. This evolution correspons to the first couple equation of the simplifie BBGKY hierarchy. Higher orer equation can be euce from the equations of motion of the fluctuating part, that are given by: δj n t δjn x = 1 χj z δj y n χj y δj z n χ[δj y n δj z n yz, t δjn y = 1 χj z δj x n χj x δj z n χ[δj x n δj z n xz, t δjn z = χj x δj y n j y δj n x χ[δj x n δj y n xy. From these equations, any equation of motion of the average quasi-spin centere moments can be erive rather easily. An illustration is given in appenix B for the evolution of secon, thir an fourth centere moments, enote respectively by ij, 3 ij an ijl. The average equation of motion greatly simplifies ue to the fact that the Hamiltonian is invariant with respect to the rotation R z = e iπĵz. Accoringly, if the initial conition is also invariant with respect to this symmetry, all moments containing an o number of components of δj x n an δj y n are equal to zero at all times. Here, we consier initial
5 states or phase-space that fulfills this conition. Then, we have j x t = j y t = 0 as well as xzt = yzt = 0. Accoringly, the average quasi-spin evolution reuces to: t j z = χ xy, 10 while the relevant secon moments evolutions are given by: t xx = 1 χj z xy χ 3 xyz, t yy = 1 χj z xy χ 3 xyz, t xy = 1 χj z yy 1 χj z 11 xx χ 3 yyz χ 3 xxz, = 8χ 3 xyz. t zz We see that the secon orer moments are couple to the thir moments. Similarly, the thir moments evolution inuces a coupling to the th moments an so on an so forth. The explicit form of the thir an fourth moments evolution are given in the appenix B. It shoul be note that these equations are obtaine in a rather straightforwar way starting from the SMF framewor. In the following, the evolution obtaine by truncating the set of equation to the secon moments, thir moments or fourth moments will be respectively referre to as Quasi- Classical TDDM QC-TDDM, QC-TDDM3 an QC- TDDM, while the original technique that samples the initial conitions an oes not assume any truncation is simply referre to SMF. We use as a benchmar the same initial conition as in Ref. [11. A system of N = 0 particles initially escribe by the Slater eterminant j, j is consiere. This initial conition correspons to the situation where all particles occupy the lowest level in the set of two-level system in the LMG moel an all single-particle spins are 1/. The initial non-zero moments corresponing to this state are given up to fourth orer by: j z = 1, xx = yy = 1 N, 3 yyz = 3 xxz = 1 1N, xxxx= yyyy = 1 3N, 16N 3 xxyy= 1 8N 3 3N, xxzz = yyzz = 1 N 3. As stresse in Ref. [11, this initial conition is particularly critical for the TDHF theory since it correspons exactly to the sale point of a spontaneous symmetry breaing. Accoringly, this initial state is just a stationary solution of TDHF while it is not stationary in the exact evolution. Below, we consier the ynamical evolution in the wea an strong coupling regime. The frontier between these two regimes is etermine by the interaction strength an by the presence χ > 1 or absence χ < 1 of a spontaneous symmetry breaing at the mean-fiel level. In practice, we assume that the wea coupling regime correspons to χ < 1, see for instance [11, 3. A. Wea coupling regime In many physical situations, we are mainly intereste in the average evolution of one-boy observables as well as their quantum fluctuations. Stanar time-epenent Hartree-Foc is usually convenient for the former one but provies a rather poor approximation for the fluctuations. We alreay now that the SMF approach can greatly improve the escription of two-boy egrees of freeom from previous applications [11, 15. We show here, that the QC-TDDM approach where only the average quasi-spin an their fluctuations are followe in time, alreay provies a great improvement compare to the original TDHF. In Fig. 1, the exact results obtaine for χ = 0.5 are compare with the SMF an QC-TDDM calculations. We see that the QC-TDDM approach oes also provie Fluctuations 30 5 0 15 10 5 0 0 8 time FIG. 1. color online Exact evolution of ispersions of quasispin operators obtaine when the initial state is j, j for χ = 0.5 blac soli line. The results of SMF an QC- TDDM are respectively shown with blue fille circles an re ashe line. Note that here the fluctuations are given by ii = N ii. Green open squares show the results of the analytical expression given in Eqs. 1 for xx an yy. the correct evolution of the secon moments xx an yy even if higher orer effects are neglecte. Only at large time, eviation with the exact solutions starts to be visible in Fig. 1. We see however that zz remains equal to zero in QC-TDDM. This coul have been anticipate from Eq. 11 where we see that the only way to have zz evolving in time is to inclue the effects of thir orer moments. The goo agreement between QC-TDDM an the exact evolution in the wea coupling regime, espite the fact that it correspons to a simplifie BBGKY hierarchy truncate at orer, is alreay an interesting result. This inirectly proves that SMF is able to grasp some important physical effects. Actually, because of the simplicity of the QC-TDDM approach, we can even obtain
6 approximate analytical solutions for xx an yy. Inee, let us assume that the average evolution of j z can be approximate by its mean fiel solution in the evolution of secon moments. Then j z j 0 = 1/ for all time. The equation of motion can be written as: t xx = Ω xy, t yy = Ω xy, t xy = Ω yy Ω xx, where Ω ± = χj 0 ± 1. Introucing the frequency ω = Ω Ω = ω 0, we immeiately euce that: t xy = ω xy, that coul be easily integrate to finally give using the fact that xy0 = 0: xxt = Ω xx0 sin ω 0 t, yyt = Ω yy0 yy 0 ω 0 xx 0 ω 0 xx 0 ω 0 yy 0 ω 0 1 sin ω 0 t. In Fig. 1, we see that the results of these analytical expressions are very close to the QC-TDDM case an henceforth to the exact case. This gives an example where the simple formula obtaine with QC-TDDM an the approximate treatment of the average evolution of j z provies interesting information in the many-boy correlate ynamics. In particular, it irectly gives the expression of the relevant frequency ω 0 that rives the two-boy evolution. In orer to escribe also the fluctuation of the z component of the quasi-spin, one nees to inclue higher orer moments effects. The equations of motion for thir an fourth moments are given in appenix B. In Fig., we show the evolution of zz for ifferent orer of truncation of the hierarchy up to QC-TDDM. In this figure, we see that incluing thir orer moments an fourth orer moments improves graually the escription of zz. In the latter case, we see that the result is almost on top of the exact result. In the bottom part of, we also show the result of SMF assuming a Gaussian sampling of the initial phasespace. We see in particular that the gaussian approximation has an impact on the quality of the result. For comparison, we also show the results of QC-TDDM3 an QC-TDDM imposing that all thir moments are initially zero as it shoul be for a Gaussian istribution. We observe again that the QC-TDDM results ten to the complete SMF case assuming initially Gaussian phasespace. These comparisons illustrate one of the point raise above. In more complicate situations of interacting many-boy systems, initial sampling beyon the Gaussian approximation might be impossible. Then, using a truncate version of the simplifie BBGKY hierarchy with non-zero thir an/or fourth orer moments might be a useful alternative. As a sie remar, we note zz zz 1.0 0.5 0.0 1.0 0.5 0.0 0 3 6 9 time 0 3 6 9 time FIG.. color online Top: Exact evolution of ispersions of zz obtaine when the initial state is j, j for χ = 0.5 blac soli line. This evolution is compare with QC-TDDM re ashe line,qc-tddm3 pin otte line an QC- TDDM green otte-ashe line results. Bottom: Same as top figure except that the QC-TDDM, QC-TDDM3 an QC-TDDM are performe assuming that all thir orer moments are zero initially. In both panels, we also show the result obtaine using the full SMF calculation with the Gaussian phase-space sampling of the initial conition from [11. that the truncation of the hierarchy an/or a wrong assumption on the initial conitions can lea to unphysical values of the observables, here negative values of the fluctuations. This is a phenomenon that also happens when the stanar BBGKY hierarchy is truncate. This problem oes not occur when the full SMF is use. B. Strong coupling regime 1. Secon orer truncation In Fig. 3, the results of QC-TDDM are compare to the exact an SMF case for χ = 1.8. It is nown
7 that the SMF result is expecte to be worse as the twoboy interaction increases. Still, we see that SMF results remain comparable with the exact results over the perio isplaye in Fig. 3. We observe that the QC-TDDM results are able to escribe the very short time evolution but then rather fast eviate from the correct evolution. In that case, the truncation to secon orer oes not seem to be appropriate. Similarly to the previous case, to unerstan the very short time ynamics, one can assume that the average evolution of j z ientifies with the mean-fiel one. Then, for χ > 1, we en-up with the set of equations: t xx = Ω xy, t yy = Ω xy, t xy = Ω yy Ω xx, that coul be solve to give: xxt = Ω xx0 yyt = yy0 yy 0ω 0 xx 0 ω 0 Ω xx 0ω 0 yy 0 ω0 sinh ω 0 t, 13. sinh ω 0 t. This ivergent behavior reflects the fact that the initial state is locate at a sale point associate to a spontaneous symmetry breaing. The above formula can grasp the short time behavior an give access to the typical instability time-scale associate to this symmetry breaing see Fig. 3. However, non-linear effects an higher orer effects are neee to escribe longer time ynamics. The complete SMF approach that inclues correlations to all orer not gives access not only to the instability time but also to longer time evolution. Fluctuations 350 300 50 00 150 100 50 0 0 1 time FIG. 3. color online Same as Fig. 1 for χ = 1.8. The results of the analytical formulas 13 are shown with green open squares. Blac soli lines an blue fille circles correspon respectively to exact an SMF results.. Thir an fourth orer truncation In Fig., the evolution obtaine with higher orer truncations QC-TDDM3 an QC-TDDM are compare with the exact solution. We see in this figure that the QC-TDDM follows the exact result for longer time compare to QC-TDDM. However, in both QC- TDDM3 an QC-TDDM, the trajectory iverges after some rather short time. We note that truncate ynam- Fluctuations 500 00 300 00 100 0 0 1 time FIG.. color online Evolution of xxt obtaine using the QC-TDDM re ashe line, QC-TDDM3 green otte line an QC-TDDM blue ot-ashe line. The exact result is shown in blac soli line. ics iverges while the full SMF calculation with initial fluctuations oesn t. This is a clear inication that either a simplifie BBGKY truncate at some orer or the sharp truncation strategy that is use here, i.e. assuming irectly fluctuations to be zero above a certain orer are not appropriate in the strong coupling regime an/or close to a spontaneous symmetry breaing case. In that case, in the absence of a clear prescription for truncation, there is no alternative than performing explicitly the initial sampling. IV. CONCLUSION In the present wor, we show that the recently propose Stochastic Mean-Fiel approach that consists in sampling initial quantum fluctuations associate to a system of interacting fermions an then follow each iniviual trajectory using the associate time-epenent mean-fiel is equivalent to a simplifie BBGKY hierarchy of equations of motion. Similarly to the stanar BBGKY approach, in this hierarchy, average evolution of one-boy egrees of freeom are couple to secon orer fluctuations that are themselves couple to thir orer fluctuation, an so on. We show that the simplifie hierarchy is retaining spe-
8 cific terms of the stanar BBGKY one. The quality of the SMF approach in many recent applications [10 1, 15, 5 inicates that the retaine terms are among the important ones beyon the mean-fiel approximation at least for short time evolution. This fining is helpful in particular to unerstan why the SMF technique can be more accurate than an approach base on the stanar BBGKY assuming a truncation at some orer, lie the TDDM approach. Inee, while some correlations are lost in SMF, it propagates correlations to any orers. Using the LMG moel as an illustrative example, we show that the simplifie BBGKY approach can be truncate an use as an alternative to the complete SMF theory in the wea coupling regime. In this regime, it is shown that it can even give better results that the complete SMF theory where a Gaussian assumption for the initial phase-space is mae. It aition, in some cases, it gives interesting physical insight beyon the inepenent particle approximation. In the strong coupling regime, a sharp truncation seems to be aequate only for very short time evolution, unstable behavior is observe in the LMG moel for longer times. Still it provies the timescale associate to the spontaneous symmetry breaing occurring in the LGM moel. For strong coupling, the complete SMF implementation remains more aequate an leas to better escription of correlations over longer time, without unstable evolution. ACKNOWLEDGMENT Y.T. acnowleges the grant receive by the label PIO. This wor is supporte in part by US DOE Grants No. DE-FG05-89ER0530. Appenix A: connection between SMF approach an BBGKY hierarchy In the present appenix, the intermeiate steps to obtain the set of couple equations 3 an 6 are given. The main assumption of the SMF theory is that all trajectories ientify with a mean-fiel evolution given by Eq. 1. Using the expression of the mean-fiel hamiltonian, we obtain: i ρn αβ t = γ t αγ ρ n γβ ρn αγ t γβ γλλ ṽ αλ,γλ ρ n λ λ ρn γβ γλλ ρ n αγ ρ n λ λ ṽ γλ,βλ. This equation coul be rewritten formally as: i ρn 1 t [ = t 1, ρ n 1 [ Tr ṽ 1, ρ n 1 ρn. A1 Note that here we use the fact that Tr ṽ 1 ρ n = Tr ρ n ṽ 1. Taing the average, we irectly obtain the first equation of the hierarchy 3. Higher orer equations are immeiately obtaine by using the property: i t [ρn 1 ρ n = = [ ρ n 1 i t ρn α ρ n ρ n 1 [t α, ρ α ρ n [ ρ n 1 Tr 1 ṽ α1 ρ n 1, ρn α ρ n Introucing the notation M n 1 = [ρn 1 ρ n n, we then en-up with the fact that the equation of motion of M couple to M n 1 1. [ i [ t M n 1 = t, M n 1 Tr 1 ṽ α1, M n 1 1. The set of equations 3 then correspon to the average version of the above couple equations. 1 is 1. Average evolution of centere moments In orer to get the equations of motion for the centere moments efine in Eq., it is first convenient to obtain the evolution of the fluctuations δρ n with respect to the average. Subtracting the average evolution of ρ n obtaine above to the Eq. A1 gives: i t δρn 1 = [t 1, δρ n 1 Tr [ [ ṽ 1, δρ n 1 ρ Tr ṽ 1, ρ 1 δρ n [ Tr ṽ 1, δρ n 1 δρn δρ n 1 δρn.
9 Similarly as in the previous section, we can then use the fact that: i t δρ n 1 δρ n = i to obtain the general equation of motion vali for : i δρ n 1 δρ n = t α, δρ n 1 δρ n t α δρ n 1 Tr 1 [ṽ α1, δρ n 1 δρ n ρ 1 Tr 1 [ṽ α1, δρ n 1 δρ n 1 δρ n α δρ n t, Tr 1 [ṽ α1, δρ n 1 ρ α δρ n δρ n 1 Tr 1 [ṽ α1, δρ n 1 C α1 δρ n. A Taing the average, we euce Eq. 6. Appenix B: Simplifie BBGKY hierarchy for the LGM moel Starting from the ynamical evolution of the fluctuations given in the main text, one can euce the evolution of the secon moments through the use of ij t = δj n j t δjn i δj n i The resulting equation are given by: t δjn t xx = 1 χj z xy χj y xz, t yy = 1 χj z xy χj x yz, t zz = 8χ[j x yz j y xz, j. B1 t xy = 1 χj z yy 1 χj z xx χj y yz χj x xz, t xz = 1 χj z yz χj y zz χj x xy χj y xx, t yz = 1 χj z xz χj x zz χj x yy χj y xy. It is worth mentioning that in the exact case, since the Hamiltonian is invariant uner the rotation R z = e iπjz, most of the moments appearing above will be zero uring the evolution. The same property hols if the initial phase-space is invariant uner this symmetry. Then, any moments that contains an o number of x or y are equal to zero. The resulting equation of motion for the average fluctuations is given in Eq. 11. Using a similar strategy, the evolution of the thir orer moments is given by: t 3 xxz = 1 χj z 3 xyz χ xx xy χ xyzz xxxy, t 3 yyz = 1 χj z 3 xyz χ yy xy χ xyzz xyyy, t 3 xyz = 1 χj z 3 yyz 1 χj z 3 xxz χ xy xy χ [ yyzz xxzz xxyy, t 3 zzz = 1χ zz xy 1χ xyzz, while the th moments evolution reas as neglecting higher orers terms: t xxxx = 1 χj z xxxy, t yyyy = 1 χj z yyyx, t zzzz = 16χ xy 3 zzz, t xxxy = 3 1 χj z xxyy 1 χj z xxxx, t yyyx = 3 1 χj z xxyy 1 χj z yyyy, t xxyy = 1 χj z xxxy 1 χj z yyyx, t yyzz = 1 χj z xyzz 8χ xy 3 yyz, t zzxx = 1 χj z xyzz 8χ xy 3 xxz, t xyzz = 1 χj z xxzz 1 χj z yyzz 8χ xy 3 xyz.
10 It shoul be note that the above equations are rather straightforwar to obtain from the SMF approach. In comparison, using stanar BBGKY approach woul be more cumbersome an woul lea to more complicate expressions. [1 N.N. Bogolyubov, J. Phys. URSS 10, 56 196. [ H. Born, H.S. Green, Proc. Roy. Soc. A 188, 10 196. [3 J.G. Kirwoo, J. Chem. Phys. 1, 180 196. [ D. Lacroix, S. Ayi, an Ph. Chomaz, Prog. Part. Nucl. Phys. 5, 97 00. [5 W. Cassing, S.J.Wang, Z. Phys. A 337, 1 1990. [6 M. Gong, M. Tohyama, Z. Phys. A 335, 153 1990. [7 K.-J. Schmitt, P.-G. Reinhar, an C. Toepffer, Z. Phys. A 336, 13 1990. [8 Mitsuru Tohyama an Peter Schuc, Eur. Phys. J. A 50, 77 01. [9 A. Abari, M.J. Hashemi, A. Rubio, R.M. Nieminen, an R. van Leeuwen, Phys. Rev. B 85, 3511 01. [10 D. Lacroix, D. Gambacurta an S. Ayi, Phys. Rev. C 87, 06130R 013. [11 D. Lacroix, S. Ayi, an B. Yilmaz, Phys. Rev. C 85, 0160R 01. [1 Denis Lacroix, S. Hermanns, C. M. Hinz, an M. Bonitz Phys. Rev. B 90, 1511 01. [13 S. Ayi, Phys. Lett. B 658, 17 008. [1 D. Lacroix an S. Ayi, Eur. Phys. J. A50, 95 01. [15 Bulent Yilmaz, Denis Lacroix, an Resul Curebal, Phys. Rev. C 90, 05617 01. [16 Sebastian Hermanns, Niclas Schlünzen, an Michael Bonitz Phys. Rev. B 90, 15111 01. [17 D. Lacroix, Phys. Rev. C 73, 0311 006. [18 S. Ayi an Y. Abe, Phys. Rev. C 6, 0609 1996. [19 G. F. Bertsch, P. F. Bortignon, an R. A. Broglia, Rev. Mo. Phys. 55, 87 1983 [0 H. J. Lipin, N. Meshov, an A. J. Glic, Nucl. Phys. A 6, 188 1965. [1 N. M. D. Agassi an H. J. Lipin, Nucl. Phys. A86, 31 1966. [ A. P. Severyuhin, M. Bener, an P.-H. Heenen, Phys. Rev. C7, 0311 006. [3 P. Ring an P. Schuc, The Nuclear Many-Boy Problem Springer-Verlag, New Yor, 1980. [ P. Bonche an H. Flocar, Nucl. Phys. A 37, 189 1985. [5 B. Yilmaz, S. Ayi, D. Lacroix, an O. Yilmaz, Phys. Rev. C 90, 0613 01. [6 Wei-Min Zhang, Da Hsuan Feng, an Robert Gilmore Rev. Mo. Phys. 6, 867 1990. [7 K.-K. Kan, P. C. Lichtner, M. Dworzeca, an J. J. Griffin, Phys. Rev. C 1, 1098 1980.