arxiv:cond-mat/ v1 31 Jan 2003

Transcription

1 Nosé-Hoover sampling of quantum entangle istribution functions D. Mentrup an J. Schnack 1 arxiv:con-mat/ v1 31 Jan 003 Abstract Universität Osnabrück, Fachbereich Physik Barbarastr. 7, D Osnabrück While thermostate time evolutions stan on firm grouns an are wiely use in classical molecular ynamics (MD) simulations [1], similar methos for quantum MD schemes are still lacking. In the special case of a quantum particle in a harmonic potential, it has been shown that the framework of coherent states permits to set up equations of motion for an isothermal quantum ynamics []. In the present article, these results are generalize to inistinguishable quantum particles. We investigate the consequences of the (anti-)symmetry of the many-particle wavefunction which leas to quantum entangle istribution functions. The resulting isothermal equations of motion for bosons an fermions contain new terms which cause Bose-attraction an Pauli-blocking. Questions of ergoicity are iscusse for ifferent coupling schemes. PACS: ; Ch; 0.70.Ns Keywors: Quantum statistics; Canonical ensemble; Ergoic behaviour; Thermostat; mixe quantum-classical system 1 Introuction an summary Classical MD simulations are performe by solving Hamilton s equations of motion numerically. Therefore, the internal energy of the system is conserve uring time evolution. If the ergoic hypothesis is satisfie, a time average of a macroscopic observable is, uner equilibrium conitions, the same as a microcanonical ensemble average. Hence, isoenergetic molecular ynamics entails a metho for the calculation of microcanonical ensemble averages. 1 jschnack@uos.e, WWW: Preprint submitte to PHYSICA A 31 January 003

2 In the canonical ensemble, however, the internal energy of the system is not constant, but free to fluctuate by thermal contact to an external heat bath. In orer to aapt classical MD simulations to the problem of the calculation of canonical ensemble averages, i. e. to pass over from an isoenergetic to an isothermal time evolution, powerful methos have been evelope since the 1980s, an they are commonly use nowaays [3,4]. Some of them, like the Nosé-Hoover chain technique [5] or the Kusnezov-Bulgac-Bauer thermostat [6], are base on eterministic time-reversal equations of motion leaing to trajectories that fill the phase space of the system accoring to the canonical thermal weight. This allows the calculation of canonical ensemble averages by means of molecular ynamics. In quantum mechanics, the problem is more involve, an the pioneering approaches of Grilli an Tosatti [7] an Kusnezov [8] were not applie very much. Moreover, the quantum mechanical time evolution itself is a har computational problem. However, a number of approximate quantum MD methos are available [9]. It is an open question whether they can be moifie in a manner appropriate to permit the calculation of quantum canonical averages. In a ifferent methoological approach to isothermal quantum ynamics, we have shown that in the special case of a quantum particle in a harmonic oscillator potential, the framework of coherent states permits to set up equations of motion for an isothermal quantum ynamics []. Following the approaches of Nosé an Hoover or the similar KBB-technique, time-epenent pseuofriction terms are ae to the equations of motion for the parameters of coherent states. The ynamics of the pseuofriction coefficients is esigne such that the esire thermal weight function is a stationary solution of a generalize Liouville equation on a mixe quantum-classical phase space. The principle of inistiguishability of ientical particles as an important quantum aspect to the problem of isothermal quantum ynamics. The present article focusses on the moifications of the metho presente in [] for a system of two non-interacting ientical particles. The (anti-)symmetry of the wavefunction leas to an entangle istribution function on quantum phase space. Surprisingly, we fin that the originally classical Nosé-Hoover approach succees to allow for this quantum feature. The entanglement leas to aitional terms in the equations of motion of the pseuofricional coefficients that act effectively like an attractive (for bosons) or repulsive force (for fermions). We examine whether the moifie equations of motion lea to an ergoic time evolution, an it turns out that the aitional terms improve the overall ergoicity of the various schemes. In aition, stanar techniques to generate ergoicity known from the classical methos are employe an teste. As a general result, ergoic time evolution is achieve even for low temperatures without serious ifficulties. However, it is inispensable to thermalize both particles since the effective forces o not suffice to yiel an ergoic time

3 evolution of the whole system if only one particle is couple to a thermostat. Metho an setup.1 Two-particle istribution function in a harmonic oscillator potential The iea of the single-particle quantum Nosé-Hoover metho [] is to moify the equations of motion of the coherent states [10] parameters r an p (which are combine to the complex parameter α = mω r + i m ω p) such that the quantum weight function of a single particle in a harmonic oscillator potential, enote by w qm (r, p), is sample in time. In precise analogy to this case, we etermine a thermal weight function w ε ()(α 1, α ) that permits to etermine canonical ensemble averages of two ientical particles [11]. The values of ε, + an, refer to bosons an fermions, respectively. We write S α 1, α A ε = S + α 1, α (1) for the (anti-)symmetrize two-particle wavefunctions. S enotes the antisymmetrizing projector (i. e., A is a two-particle Slater eterminant), S + the symmetrizing projector. The starting point for the calculation of w ε () is the following expression for the calculation of a thermal expectation value B as a phase space integral, B = 1 Z () ε = 1 Z () ε tr ( B e βh ) α 1 π α π α 1, α S ε B e βh α1, α, () with Z () ε = tr ( e βh ) being the respective partition function. Note that we have roppe the projector S ε acting upon the ket using the iempotency of projectors. Decomposition of e βh = e βh / e βh / an a cyclic shift of the operators uner the trace enables further simplification, taking avantage of 3

4 the specific properties of coherent states [1]. The following expression B = 1 α 1 α Z ε () π π A e α 1 (e β ω 1) e α (e β ω 1) ε B A ε A ε A ε }{{} A ε A ε w ε () (α 1, α ) (3) finally efines w () ε as the thermal weight of the expectation value Aε B Aε A ε A ε. w ε () (α 1, α ) is not merely the prouct of two one-particle thermal weight functions (as it woul be the case for two istinguishable particles), but contains in aition the factor A ε A ε that accounts for the quantum effects of inistinguishability. Moreover, w ε ()(α 1, α ) cannot be written as a prouct of two functions epening only on α 1 an α, respectively, since the term A ε A ε cannot be separate in this way. This is a result of the quantum mechanical principle of inistinguishability. Therefore, we say that w ε () is entangle. In the case of fermions, we easily fin A A = 1 (1 e α 1 α ). This expression vanishes if α 1 = α. A quantum state with two ientical fermions in the same one-particle-state is forbien by the Pauli exclusion principle, an therefore oes not contribute to a thermal average. In contrast, for bosons we have A + A + = 1 (1 + e α 1 α ), which contains the opposite sign that enhances the thermal weight of the quantum state with two bosons in the same one-particle state. The thermal istribution function w ε () yiels the correct partition function for the thermal average 1. To show this, we calculate tre βh α 1 α = π π e α 1 (eβ ω 1) e α (eβ ω 1) 1 (1 + εe α 1 α ) (4) = 1 ( ) α 1 π e α 1 (eβ ω 1) + ε 1 ( α 1 ) α π π e α 1 (e β ω 1) e α (e β ω 1) e α 1 α = 1 ( 1 + ε e βω 1) 1 ( ) 1 e βω 1 = 1 [ Z (1) (β) + εz (1) (β) ], which is a correct, well-known recursion relation for the two-particle partition function. We inicate that the secon integral is solve most easily by a change of variables from α 1, α to α + = 1 (α 1 + α ), α = 1 (α 1 α ), which leas to a separation of the ouble integral. 4

5 . Moification of the equations of motion The Nosé-Hoover equations of motion for the coherent states parameters that we will investigate rea t r 1 = p 1 m, t r = p m, (5) t p 1 = mω p η1 r 1 p 1, Q 1 t p = mω p η r p. Q The time epenence of the pseuofriction coefficients has to be etermine in a proceure analog to the case of a single particle, i. e. we require that the esire istribution function f () ε (α 1, α, p η1, p η ) w () ε (α 1, α ) exp ( β ( p η1 + p ) ) η (6) Q 1 Q = e ( α 1 + α )(e β ω 1) 1 (1 + εe α 1 α ) ( ( p η1 exp β + p ) ) η Q 1 Q = e U (1 ± e V ) is a stationary solution of a generalize Liouville equation on the phase space with elements x = (α 1, α, p η1, p η ) (r 1, p 1, r, p, p η1, p η ). The abbreviations U an V are efine as U = ( α 1 + α )(e β ω 1) + β = ( p η1 Q 1 + p η Q ( 1 mω (r1 + r ) + 1 ) e β ω m (p 1 + p ) 1 + β ω V = α 1 α = mω (r 1 r ) + 1 m ω (p 1 p ). ) ( p η1 Q 1 + p η Q As in the one-particle case, p η1 an p η are regare as classical pseuofriction coefficients. For the Liouville equation, we calculate t f() ε (α 1, α, p η1, p η ) = ( U + V )f ± + V e U (8) ), (7) 5

6 an ( ) ( x ẋ = ṙ 1 + ṙ + ṗ 1 + ṗ + ṗ η1 + ) ṗ η r 1 r p 1 p p η1 p η = p η 1 Q 1 + p η Q, using the equations of motion (5) an imposing, as common in the present context [6], the constraint ṗ ηi / p ηi = 0. After further transformations, we obtain the following equations of motion for the pseuofriction coefficients from a comparison of the coefficients of the terms p η1 /Q 1 an p η /Q on both sies of the Liouville equation: ṗ η1 = 1 ( p 1 β m ṗ η = 1 β e β ω 1 ω 1 + εp 1 p 1 p m ω ( p e β ω 1 p 1 p 1 εp m ω m ω (9) ) 1, (10) e V + ε1 ) 1. e V + ε1 These equations, along with the set of equations (5), form a genuine quantum Nosé-Hoover thermostat for two ientical quantum particles. The part p i e β ω 1 1, i = 1,, in the equations of motion is familiar from the ynamics of a single thermostatte particle. However, in the present case of two m ω particles we fin aitional terms that reflect the effects of Bose-attraction an Pauli-blocking irectly in the thermostatte ynamics, see next section. The set of equations of motion (5), (10) conserves the quantity H = 1 β lnw() ε (α 1, α ) + p η 1 + p η + 1 t t Q 1 Q β ( pη1 (t ) + p η (t ) ) Q 1 Q. (11) The analogous approach using a KBB-scheme starts with the set of equations t r 1 = p 1 m g r 1 (ξ 1 )F r1 (r 1, p 1 ) (1) t p 1 = mω r 1 g p 1 (ζ 1 )G p1 (r 1, p 1 ) t r = p m g r (ξ )F r (r, p ) t p = mω r g p (ζ )G p (r, p ), 6

7 with the conserve quantity H = 1 β ln w() ε (α 1, α ) + g r 1 (ξ 1 ) + g r (ξ ) + g p 1 (ζ 1 ) + g p (ζ ) κ r1 κ r κ p1 κ p + 1 ( t t Fr1 g r β r 1 + F r g r 1 r + G p 1 g p p 1 + G ) p g p 1 p. (13) The functions g r1, g p1, g r, g p are chosen such that the canonical istribution of the emons can be normalize. The functions F r1, F r, G p1, G p are arbitrary. This scheme has the obvious avantage that positions an momenta are treate symmetrically, i. e. pseuofriction coefficients are present in all equations of motion. The time epenence of the pseuofriction coefficients that is obtaine from the Liouville equation in the phase space with elements x = {r 1, p 1, r, p, ξ 1, ξ, ζ 1, ζ } reas t ζ 1 = κ ( p 1 p1 β m G e β ω 1 p 1 G ) p 1 p 1 p 1 + εg p1 (14) ω p 1 m ω e V + ε1 t ζ = κ ( p p β m G e β ω 1 p G ) p p 1 p 1 εg p ω p m ω e V + ε1 t ξ 1 = κ ( r 1 mω e β ω 1 r 1 F r1 F ) r εf r1 m ω(r 1 r ) β ω r 1 e V + ε1 t ξ = κ ( r mω e β ω 1 r F r F ) r 1 εf r m ω(r 1 r ). β ω r e V + ε1 The effects of Pauli-blocking an Bose-attraction are now present in the equations of motion of all pseuofriction coefficients. 3 Results 3.1 Bose-attraction an Pauli-blocking The ifferent signs for bosons an fermions in the equations of motion of the pseuofriction coefficients stem from the respective two-particle wavefunction. We investigate the consequences for the movements of the particles. Obviously, the effects of Bose-attraction an Pauli-blocking will be most pronounce at low temperatures, when both particles ten to occupy the one-particle groun state an thereby get close to one another in phase space. We examine the most elementary example, the scheme given by the set of equations of motion (5) an (10). When two fermions are close in phase space, V = α 1 α becomes very small an the factor 1/(e V 1) gets very large, causing a strong acceleration of p η1 an p η, an thereby of p 1 an p into 7

8 Fig. 1. Left panel: isothermal ynamics of two bosons, right panel: two fermions. The initial values of {r 1,p 1,r,p } are {0, 0.1,0,0.1} in both cases, the temperature is T = 0.1 ω, an Q 1 = Q = 0.5. The integration time is 8.1 perios of the harmonic oscillator for the bosons, an 6.5 for the fermions. The time istance between the symbols is The figure illustrates the effect of the Bose attraction an the Pauli-blocking in the ynamics: While the bosons stay close to each other for a longer time, the fermions are immeiately riven away from each other ue to the exclusion principle. opposite irections ue to the ifferent signs. Effectively, a close approach of the particles in phase space, corresponing to V 0, is avoie by the ynamics. Thus, in the case of fermions, the aitional terms in the equations of motion (10) act like a repulsive force. In the case of bosons, the opposite signs in equations (5), (10) cause an acceleration of the parameters in the irection of one another, favoring a meeting of the particles in phase space. So the case V = 0 is not exclue at all; on the contrary, it is aime at with the maximum value of 1/(e V + 1). Fig. 1 illustrates the typical behaviour of two ientical particles neighboring in phase space. In essence, the resulting effect of inistinguishability on the thermostatte ynamics of ientical quantum particles looks like an attractive or repulsive interaction, although we treat a system of non-interacting particles. The interaction which is of purely statistical origin is meiate by the influence of the pseuofrictional forces. Nevertheless, one can hope that this statistical interaction influences the ergoicity of the system. As an extreme example, one coul think of thermalizing only one particle an hope that the secon one thermalizes ue to the interaction, see next section. 3. Ergoicity The classical Nosé-Hoover metho applie to a single particle in a harmonic oscillator features ergoicity problems that are encountere in the quantum 8

9 case, too, since the overall characteristics of the equations of motion are similar in both cases. However, there is no classical counterpart of the quantum statistical interaction between ientical particles. Therefore, the question of whether the more complex equations of motion presente in the preceent section lea to improve ergoic behavior eserves attention. For the ifferent methos, we look at marginal istributions of the thermal weight function f ε (). As an example, we give the analytical expression of the marginal istribution of r 1 : f ε (r 1 ) = 1 Z () ε = 1 Z () ε p1 r p π π mω 1 8π (e β ω 1) 3 e ε e β ω 1 w() ε (r 1, r, p 1, p ) (15) mω (eβ ω 1)r 1 mω e eβ ω sinh(β ω)r Nosé-Hoover an Nosé-Hoover chain metho Firstly, we present results that are obtaine using the equations of motion (5), (10), i. e. the original Nosé-Hoover technique generalize to the case of two ientical fermions. A etaile investigation of these equations of motion for a wie temperature range shows that the system is not ergoic in general. We fin non-ergoic motion at low an intermeiate temperatures T 1.. However, above that value, we fin ergoic motion, an the respective marginal istributions are well matche by the histograms obtaine by time averaging, see Fig.. This fact is to be contraste sharply to the finings in the singleparticle case, where the simple Nosé-Hoover scheme prouces non-ergoic motion even at very high temperature values [,3]. The ifferent behavior of the two-particle ynamics may be attribute to the statistical interaction appearing in this case. The more complicate form of the equations of motion improves the ergoicity of the scheme substantially compare to the case of a single particle. As a result, this technique turns out to be not recommenable in the lowtemperature regime. Therefore, we have investigate whether the problems of non-ergoicity can be overcome by stanar methos that are known from classical molecular ynamics, e. g. with a chain of thermostats [5]. Using one aitional thermostatting pseuofriction coefficient for each p ηi, i = 1,, the problems of ergoicity are immeiately an reliably resolve. The marginal istributions are exactly matche by the respective histograms, an the system is ergoic within the whole large temperature range that has been 9

10 Fig.. Results of time averaging with the simple Nosé-Hoover scheme for two fermions, left panel: T = 0.1, right panel: T = 1.. In both cases, we use ientical initial conitions, r 1 (0) = r (0) = 0.01, p 1 (0) = p (0) = 0.01, p η1 (0) = p η (0) = 0.01, Q 1 = Q = 1.0. We present ensity plots an histograms for r 1, p 1, an p η1 ; the respective results for r, p, an p η look similar. investigate [13]. As in the classical case, it is remarkable that a solution to the seemingly inaccessible problem of ergoicity is so easily obtaine. Moreover, from the point of view of computational time, it is relatively inexpensive, since only two more egrees of freeom are ae, which is negligible especially if 10

11 one eals with larger systems. We note that it is not sufficient to couple a secon thermostating pseuofriction coefficient to only one parameter, say p η1. In this case, the marginals of p η are not sample correctly at temperature values T 0.3. Even worse, if only p 1 is couple to a thermostat, the ynamics of the parameters r an p is not affecte at all. Hence, the energy of particle is an aitional conserve quantity, leaing to strongly non-ergoic motion. This illustrates the limits of the improve ergoicity ue to the statistical interaction. 3.. KBB metho with ifferent coupling schemes In the original paper [6], Kusnezov, Bulgac, an Bauer propose the following choice of functions g r1 (ξ 1 ) = 1 ξ 1, g r (ξ ) = 1 ξ, g p1 (ζ 1 ) = 1 4 ζ4 1, g p (ζ ) = 1 4 ζ4, (16) as a reliable scheme that provies ergoic behavior. Using their rules of thumb for the aaptation of the respective coupling constants, we fin that this cubic coupling scheme leas to ergoic motion also in the present case. However, consiering the low-temperature results for fermions presente in Fig. 3, it can be inferre that the marginal istributions of the pseuofriction coefficients ξ 1, an equivalently, ξ (not shown), that are linearly couple to the system are not sample properly in the range T 0.1. The histograms presente in Fig. 3 clearly inicate that ξ 1 oes not epart substantially from its initial value uring the time evolution. Nevertheless, the thermal weight function w ε () is sample correctly, since all marginal istributions coincie with the exact result. However, it is obvious that the cubic coupling scheme is not ergoic in a part of the phase space at low temperatures. This applies to bosons as well. We propose to resolve this problem by coupling the coefficents ξ 1 an ξ in precisely the same manner as ζ 1 an ζ, i. e. cubically. This appears to be a sensible improvement since the istributions of ζ 1 an ζ are sample correctly even at very low temperatures. The resulting equations of motion yiel the outcome shown in the right han panel of Fig. 3. All marginal istributions are now sample correctly. Recall that the erivatives of the functions given in (16) appear in the equations of motion. 11

12 Fig. 3. Marginal istributions obtaine at T = 0.05 with a cubic coupling scheme (left panel), an with the alternate scheme escribe in the text (right panel). The initial value of ξ 1 on the left panel is ξ 1 (0) = 5.1, an the freezing of the emon variable is striking. Surprisingly, all other marginals are well reprouce. The right panel illustrates that the problem is resolve by the moifie scheme. 1

13 3.3 Two-particle ensity As an example for the etermination of the canonical average of a two-particle observable, we present results of time averaging for the two-particle ensity, ρ () ε (x 1, x ) = x 1, x A ε A ε A ε, (17) which can also be calculate analytically. In terms of the relative an centerof-mass variables x = 1 (x 1 x ), x + = 1 (x 1 + x ), the result is given by ρ () ε (x, x + ) = 1 ( mω 1 Z ε () π e β ω 1 exp mω ( ) ) tanh 1 β ω (x + + x ) ( ( 1 + ε exp 4 mω e β ω ) ) e β ω 1 x. (18) In orer to give a useful representation of our finings, we present results of time averages along with the respective analytical results for fixe values of x + that have been chosen arbitrarily. Fig. 4 shows an excellent agreement. Fig. 4. Results of time averages for the bosonic (left figure) an fermionic (right figure) two-particle ensity. The soli curves correspon to the respective analytical results as given by equation (18). In the fermionic case, the crosses are obtaine from the simulations that gave the right han panels of figure 3. The upper line correspons to the value x + = 0, the lower line to the value x + = Summary an outlook This article presents an extension of the powerful methos of heat bath coupling in classical MD simulations to a quantum system of ientical particles. Surprisingly, the classical ansatz turns out to be suitable for the sampling of quantum entangle istribution functions. The resulting two-particle ynamics contains aitional terms that act like an attractive (bosons) or repulsive 13

14 (fermions) force meiate by the pseuofriction coefficients. Ergoicity problems are reliably resolve by Nosé-Hoover chains or a moification of the cubic coupling scheme, even at very low temperature values. An analytical treatment of the case of an N-particle Fermi system is possible, an the resulting equations are given in Ref. [13]. The most promising prospect that our metho offers lies in a combination with approximate quantum ynamics schemes. A variety of such schemes is available, some of which are base on the time-epenent quantum variational principle [14] which allows to erive approximations to the time-epenent Schröinger equation. How can we combine the thermostating metho evelope in this work with, e. g., Fermionic Molecular Dynamics (FMD) [9] in orer to obtain an isothermal ynamical scheme for a complex interacting fermion system? Does our generalization of the classical methos to quantum ynamics permit to make the power of approximate quantum MD schemes available for the calculation of quantum canonical averages without iagonalising the full many-boy Hamiltonian? This question is very timely in view of recent experiments investigating the behavior of trappe Fermi gases. In FMD, the variational trial state is a Slater eterminant of single-particle Gaussian wave packets parametrize by mean position, mean momentum, an complex with. These wave packets are frequently referre to as squeeze states an may be regare as generalisations of coherent states. Therefore, an application of the thermostats evelope in the present work to FMD appears feasible. Another iea of combining FMD with a thermostat is to cool the system of interest sympathetically, i. e. via an interaction between particles that are kept at a constant temperature by a quantum Nosé-Hoover chain an the physical system uner investigation. This correspons precisely to the experimental technique of sympathetic cooling currently employe to investigate ultracol fermionic gases [15,16]. Since the thermalizing of the particles couple to a Nosé-Hoover chain woul correspon to a thermalizing of non-interacting particles, such a metho is only approximately correct since we nee to employ an interaction to enable the sympathetic cooling. Despite these ifficulties, an important asset of an isothermal MD scheme is that it can possibly provie temporal information, in particular time correlation functions. Although it is not clear to which egree the extene system methos realistically mimic the heat bath interaction, the unerlying equations of motion are physically reasonable. Therefore, in principle, this metho is tailor-mae to moel the particle ynamics at constant temperature in a magnetic trap. Acknowlegments 14

15 The authors woul like to thank the Deutsche Forschungsgemeinschaft (DFG) for financial support of the project Isothermal ynamics of small quantum systems. References [1] M. E. Tuckerman, G. J. Martyna, J. Phys. Chem. B 104 (000) 159 [] D. Mentrup, J. Schnack, Physica A 97 (001) 337 [3] W. G. Hoover, Phys. Rev. A31 (1985) 1685 [4] S. Nosé, Prog. Theor. Phys. Suppl. 103 (1991) 1 [5] G. J. Martyna, M. L. Klein, M. Tuckerman, J. Chem. Phys. 97 (199) 635 [6] D. Kusnezov, A. Bulgac, W. Bauer, Ann. of Phys. 04 (1990) 155 [7] M. Grilli, E. Tosatti, Phys. Rev. Lett. 6 (1989) 889 [8] D. Kusnezov, Phys. Lett. A184 (1993) 50 [9] H. Felmeier, J. Schnack, Rev. Mo. Phys. 7 (000) 655 [10] J. R. Klauer, B.-S. Skagerstam, Coherent states, Worl Scientific Publishing, Singapore 1985 [11] J. Schnack, PhD thesis, TH Darmstat, 1996 [1] J. Schnack, Europhys. Lett. 45 (1999) 647 [13] D. Mentrup, PhD thesis, Universität Osnabrück, 003 [14] P. Kramer, M. Saraceno, Geometry of the Time-Depenent Variational Principle in Quantum Mechanics, Lecture Notes in Physics 140, Springer, Berlin (1981) [15] C. J. Myatt, E. A. Burt, R. W. Ghrist, E. A. Cornell, an C. E. Wieman, Phys. Rev. Lett 78 (1997) 586 [16] B. DeMarco, D. S. Jin, Phys. Rev. A 58 (1998) R467 15