Ivana Vidanović

Geometric resonances in Bose-Einstein condensates with two- and three-body interactions

We investigate geometric resonances in Bose?Einstein condensates by solving the underlying time-dependent Gross?Pitaevskii equation for systems with two- and three-body interactions in an axially symmetric harmonic trap. To this end, we use a recently developed analytical method (Vidanovi? et al 2011 Phys. Rev. A 84 013618), based on both a perturbative expansion and a Poincar??Lindstedt analysis of a Gaussian variational approach, as well as a detailed numerical study of a set of ordinary differential equations for variational parameters. By changing the anisotropy of the confining potential, we numerically observe and analytically describe strong nonlinear effects: shifts in the frequencies and mode coupling of collective modes, as well as resonances. Furthermore, we discuss in detail the stability of a Bose?Einstein condensate in the presence of an attractive two-body interaction and a repulsive three-body interaction. In particular, we show that a small repulsive three-body interaction is able to significantly extend the stability region of the condensate.

Ultra-fast converging path-integral approach for rotating ideal Bose–Einstein condensates

A recently developed efficient recursive approach for analytically calculating the short-time evolution of the one-particle propagator to extremely high orders is applied here for numerically studying the thermodynamical and dynamical properties of a rotating ideal Bose gas of 87 Rb atoms in an anharmonic trap. At first, the one-particle energy spectrum of the system is obtained by diagonalizing the discretized short-time propagator. Using this, many-boson properties such as the condensation temperature, the ground-state occupancy, density profiles, and time-of-flight absorption pictures are calculated for varying rotation frequencies. The obtained results improve previous semiclassical calculations, in particular for smaller particle numbers. Furthermore, we find that typical time scales for a free expansion are increased by an order of magnitude for the delicate regime of both critical and overcritical rotation.

Recursive Schrödinger equation approach to faster converging path integrals.

By recursively solving the underlying Schrödinger equation, we set up an efficient systematic approach for deriving analytic expressions for discretized effective actions. With this, we obtain discrete short-time propagators for both one and many particles in arbitrary dimension to orders that have not been accessible before. They can be used to substantially speed up numerical Monte Carlo calculations of path integrals, as well as for setting up an alternative analytical approximation scheme for energy spectra, density of states, and other statistical properties of quantum systems.

Parametric and geometric resonances of collective oscillation modes in Bose-Einstein condensates

We analytically and numerically study nonlinear dynamics in Bose-Einstein condensates (BECs) induced either by a harmonic modulation of the interaction or by the geometry of the trapping potential. To analytically describe BEC dynamics, we use a perturbative expansion based on the Poincaranalysis of a Gaussian variational ansatz, whereas in the numerical approach we use numerical solutions of both a variational system of equations and the full time-dependent Gross-Pitaevskii equation. The harmonic modulation of the atomic s-wave scattering length of a BEC of 7 Li was achieved recently via Feshbach resonance, and such a modulation leads to a number of nonlinear effects, which we describe within our approach: mode coupling, higher harmonics generation and significant shifts in the frequencies of collective modes. In addition to the strength of atomic interactions, the geometry of the trapping potential is another key factor for the dynamics of the condensate, as well as for its collective modes. The asymmetry of the confining potential leads to important nonlinear effects, including resonances in the frequencies of collective modes of the condensate. We study in detail such geometric resonances and derive explicit analytic results for frequency shifts for the case of an axially symmetric condensate with two- and three-body interactions. Analytically obtained results are verified by extensive numerical simulations.

Fast converging path integrals for time-dependent potentials: I. Recursive calculation of short-time expansion of the propagator

We calculate the short-time expansion of the propagator for a general many-body quantum system in a time-dependent potential to orders that have not yet been accessible before. To this end the propagator is expressed in terms of a discretized effective potential, for which we derive and analytically solve a set of efficient recursion relations. Such a discretized effective potential can be used to substantially speed up numerical Monte-Carlo simulations for path integrals, or to set up various analytic approximation techniques to study dynamic properties of quantum systems in timedependent potentials. The analytically derived results are numerically verified by treating several simple one-dimensional models.

Properties of quantum systems via diagonalization of transition amplitudes. II. Systematic improvements of short-time propagation.

In this paper, building on a previous analysis [I. Vidanović, A. Bogojević, and A. Belić, preceding paper, Phys. Rev. E 80, 066705 (2009)] of exact diagonalization of the space-discretized evolution operator for the study of properties of nonrelativistic quantum systems, we present a substantial improvement to this method. We apply recently introduced effective action approach for obtaining short-time expansion of the propagator up to very high orders to calculate matrix elements of space-discretized evolution operator. This improves by many orders of magnitude previously used approximations for discretized matrix elements and allows us to numerically obtain large numbers of accurate energy eigenvalues and eigenstates using numerical diagonalization. We illustrate this approach on several one- and two-dimensional models. The quality of numerically calculated higher-order eigenstates is assessed by comparison with semiclassical cumulative density of states.

Fast convergence of path integrals for many-body systems

Abstract We generalize a recently developed method for accelerated Monte Carlo calculation of path integrals to the physically relevant case of generic many-body systems. This is done by developing an analytic procedure for constructing a hierarchy of effective actions leading to improvements in convergence of N -fold discretized many-body path integral expressions from 1 / N to 1 / N p for generic p . In this Letter we present explicit solutions within this hierarchy up to level p = 5 . Using this we calculate the low lying energy levels of a two particle model with quartic interactions for several values of coupling and demonstrate agreement with analytical results governing the increase in efficiency of the new method. The applicability of the developed scheme is further extended to the calculation of energy expectation values through the construction of associated energy estimators exhibiting the same speedup in convergence.

Fast converging path integrals for time-dependent potentials: II. Generalization to many-body systems and real-time formalism

Based on a previously developed recursive approach for calculating the short-time expansion of the propagator for systems with time-independent potentials and its time-dependent generalization for simple single-particle systems, in this paper we present a full extension of this formalism to a general quantum system with many degrees of freedom in a time-dependent potential. Furthermore, we also present a recursive approach for the velocity-independent part of the effective potential, which is necessary for calculating diagonal amplitudes and partition functions, as well as an extension from the imaginary-time formalism to the real-time one, which enables to study the dynamical properties of quantum systems. The recursive approach developed here allows an analytic derivation of the short-time expansion to orders that have not been accessible before, using the implemented SPEEDUP symbolic calculation code. The analytically derived results are extensively numerically verified by treating several models in both imaginary and real time.

Properties of quantum systems via diagonalization of transition amplitudes. I. Discretization effects.

We analyze the method for calculation of properties of nonrelativistic quantum systems based on exact diagonalization of space-discretized short-time evolution operators. In this paper we present a detailed analysis of the errors associated with space discretization. Approaches using direct diagonalization of real-space discretized Hamiltonians lead to polynomial errors in discretization spacing Delta . Here we show that the method based on the diagonalization of the short-time evolution operators leads to substantially smaller discretization errors, vanishing exponentially with 1/Delta(2). As a result, the presented calculation scheme is particularly well suited for numerical studies of few-body quantum systems. The analytically derived discretization errors estimates are numerically shown to hold for several models. In the follow up paper [I. Vidanović, A. Bogojević, A. Balaz, and A. Belić, Phys. Rev. E 80, 066706 (2009)] we present and analyze substantial improvements that result from the merger of this approach with the recently introduced effective-action scheme for high-precision calculation of short-time propagation.

SPEEDUP Code for Calculation of Transition Amplitudes via the Effective Action Approach

We present Path Integral Monte Carlo C code for calculation of quantum mechanical transition amplitudes for 1D models. The SPEEDUP C code is based on the use of higher-order short-time effectiveactions and implemented to the maximal order p=18 in the time of propagation (Monte Carlo time step), which substantially improves the convergence of discretized amplitudes to their exact continuum values. Symbolic derivation of higher-order effective actions is implemented in SPEEDUP Mathematica codes, using the recursive Schrodinger equation approach. In addition to the general 1D quantum theory, developed Mathematica codes are capable of calculating effective actions for specific models, for general 2D and 3D potentials, as well as for a general many-body theory in arbitrary number of spatial dimensions.