Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userPath Integral Molecular Dynamics for Bosons
88
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
Trapped Bosons exhibit fundamental physical phenomena and are potentially useful for quantum technologies. We present a method for simulating Bosons using path integral molecular dynamics. A main challenge for simulations is including all permutations due to exchange symmetry. We show that evaluation of the potential can be done recursively, avoiding explicit enumeration of permutations, and scales cubically with system size. The method is applied to Bosons in a 2D trap and agrees with essentially exact results. An analysis of the role of exchange with decreasing temperature is also presented.
rate research
Read More
We present a method for performing path integral molecular dynamics (PIMD) simulations for fermions and address its sign problem. PIMD simulations are widely used for studying many-body quantum systems at thermal equilibrium. However, they assume that the particles are distinguishable and neglect bosonic and fermionic exchange effects. Interacting fermions play a key role in many chemical and physical systems, such as electrons in quantum dots and ultracold trapped atoms. A direct sampling of the fermionic partition function is impossible using PIMD since its integrand is not positive definite. We show that PIMD simulations for fermions are feasible by employing our recently developed method for bosonic PIMD and reweighting the results to obtain fermionic expectation values. The approach is tested against path integral Monte Carlo (PIMC) simulations for up to 7 electrons in a two-dimensional quantum dot for a range of interaction strengths. However, like PIMC, the method suffers from the sign problem at low temperatures. We propose a simple approach for alleviating it by simulating an auxiliary system with a larger average sign and obtaining an upper bound to the energy of the original system using the Bogoliubov inequality. This allows fermions to be studied at temperatures lower than would otherwise have been feasible using PIMD, as demonstrated in the case of a three-electron quantum dot. Our results extend the boundaries of PIMD simulations of fermions and will hopefully stimulate the development of new approaches for tackling the sign problem.
We investigate the continuum limit that the number of beads goes to infinity in the ring polymer representation of thermal averages. Studying the continuum limit of the trajectory sampling equation sheds light on possible preconditioning techniques for sampling ring polymer configurations with large number of beads. We propose two preconditioned Langevin sampling dynamics, which are shown to have improved stability and sampling accuracy. We present a careful mode analysis of the preconditioned dynamics and show their connections to the normal mode, the staging coordinate and the Matsubara mode representation for ring polymers. In the case where the potential is quadratic, we show that the continuum limit of the preconditioned mass modified Langevin dynamics converges to its equilibrium exponentially fast, which suggests that the finite-dimensional counterpart has a dimension-independent convergence rate. In addition, the preconditioning techniques can be naturally applied to the multi-level quantum systems in the nonadiabatic regime, which are compatible with various numerical approaches.
Expanding upon previous work, using the path-integral formalism we derive expressions for the one-particle reduced density matrix and the two-point correlation function for a quadratic system of bosons that interact through a general class of memory kernels. The results are applied to study the density, condensate fraction and pair correlation function of trapped bosons harmonically coupled to external distinguishable masses.
For a system of bosons that interact through a class of general memory kernels, a recurrence relation for the partition function is derived within the path-integral formalism. This approach provides a generalization to previously known treatments in the literature of harmonically coupled systems of identical particles. As an example the result is applied to the specific heat of a simplified model of an open quantum system of bosons, harmonically coupled to a reservoir of distinguishable fictitious masses.
We propose a generalization of the Feynman path integral using squeezed coherent states. We apply this approach to the dynamics of Bose-Einstein condensates, which gives an effective low energy description that contains both a coherent field and a squeezing field. We derive the classical trajectory of this action, which constitutes a generalization of the Gross Pitaevskii equation, at linear order. We derive the low energy excitations, which provides a description of second sound in weakly interacting condensates as a squeezing oscillation of the order parameter. This interpretation is also supported by a comparison to a numerical c-field method.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
