ترغب بنشر مسار تعليمي؟ اضغط هنا

Fortran and C programs for the time-dependent dipolar Gross-Pitaevskii equation in an anisotropic trap

153   0   0.0 ( 0 )
 نشر من قبل Sadhan Adhikari K
 تاريخ النشر 2015
  مجال البحث فيزياء
والبحث باللغة English




اسأل ChatGPT حول البحث

Many of the static and dynamic properties of an atomic Bose-Einstein condensate (BEC) are usually studied by solving the mean-field Gross-Pitaevskii (GP) equation, which is a nonlinear partial differential equation for short-range atomic interaction. More recently, BEC of atoms with long-range dipolar atomic interaction are used in theoretical and experimental studies. For dipolar atomic interaction, the GP equation is a partial integro-differential equation, requiring complex algorithm for its numerical solution. Here we present numerical algorithms for both stationary and non-stationary solutions of the full three-dimensional (3D) GP equation for a dipolar BEC, including the contact interaction. We also consider the simplified one- (1D) and two-dimensional (2D) GP equations satisfied by cigar- and disk-shaped dipolar BECs. We employ the split-step Crank-Nicolson method with real- and imaginary-time propagations, respectively, for the numerical solution of the GP equation for dynamic and static properties of a dipolar BEC. The atoms are considered to be polarized along the z axis and we consider ten different cases, e.g., stationary and non-stationary solutions of the GP equation for a dipolar BEC in 1D (along x and z axes), 2D (in x-y and x-z planes), and 3D, and we provide working codes in Fortran 90/95 and C for these ten cases (twenty programs in all). We present numerical results for energy, chemical potential, root-mean-square sizes and density of the dipolar BECs and, where available, compare them with results of other authors and of variational and Thomas-Fermi approximations.



قيم البحث

اقرأ أيضاً

We present Open Multi-Processing (OpenMP) version of Fortran 90 programs for solving the Gross-Pitaevskii (GP) equation for a Bose-Einstein condensate in one, two, and three spatial dimensions, optimized for use with GNU and Intel compilers. We use t he split-step Crank-Nicolson algorithm for imaginary- and real-time propagation, which enables efficient calculation of stationary and non-stationary solutions, respectively. The present OpenMP programs are designed for computers with multi-core processors and optimized for compiling with both commercially-licensed Intel Fortran and popular free open-source GNU Fortran compiler. The programs are easy to use and are elaborated with helpful comments for the users. All input parameters are listed at the beginning of each program. Different output files provide physical quantities such as energy, chemical potential, root-mean-square sizes, densities, etc. We also present speedup test results for n
We describe a method for evolving the projected Gross-Pitaevskii equation (PGPE) for an interacting Bose gas in a harmonic oscillator potential, with the inclusion of a long-range dipolar interaction. The central difficulty in solving this equation i s the requirement that the field is restricted to a small set of prescribed modes that constitute the low energy c-field region of the system. We present a scheme, using a Hermite-polynomial based spectral representation, that precisely implements this mode restriction and allows an efficient and accurate solution of the dipolar PGPE. We introduce a set of auxiliary oscillator states to perform a Fourier transform necessary to evaluate the dipolar interaction in reciprocal space. We extensively characterize the accuracy of our approach, and derive Ehrenfest equations for the evolution of the angular momentum.
We report a theoretical description of the synthetic momentum-state lattices with a 3D Gross-Pitaevskii equation (GPE), where both the external trap potential and the mean-field spatial-density-dependent many-body interactions are naturally included and exactly treated. The GPE models exhibit better performance than the tight-binding model to depict the experimental observations. Since the trap modifies the dispersion relation for free particles and shapes the spatial density distribution that leads to inhomogeneous interactions, decoherences (damping oscillation) appear even for a short-time evolution. Our parametric calculations for the two-state oscillation suggest that we should work with a relatively shallow trap in the weakly interacting regime, especially when the long-term dynamics are concerned. The impact of the mean-field interaction, i.e., the self-trapping behavior, on the transport dynamics and the topological phase transition in a finite multiple-state lattice chain is also specifically investigated. Such an accurate treatment of the inhomogeneous interactions allows for further investigations on the interplay with disorder, the pair correlation dynamics, and the thermalization process in momentum space.
We consider an effective scaling approach for the free expansion of a one-dimensional quantum wave packet, consisting in a self-similar evolution to be satisfied on average, i.e. by integrating over the coordinates. A direct comparison with the solut ion of the Gross-Pitaevskii equation shows that the effective scaling reproduces with great accuracy the exact evolution - the actual wave function is reproduced with a fidelity close to unity - for arbitrary values of the interactions. This result represents a proof-of-concept of the effectiveness of the scaling ansatz, which has been used in different forms in the literature but never compared with the exact evolution.
التعليقات
جاري جلب التعليقات جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
mircosoft-partner

هل ترغب بارسال اشعارات عن اخر التحديثات في شمرا-اكاديميا