We demonstrate the application of transition state theory to wave packet dynamics in metastable Schrödinger systems which are approached by means of a variational ansatz for the wavefunction and whose dynamics is described within the framework of a time-dependent variational principle. The application of classical transition state theory, which requires knowledge of a classical Hamilton function, is made possible by mapping the variational parameters to classical phase space coordinates and constructing an appropriate Hamiltonian in action variables. This mapping, which is performed by a normal form expansion of the equations of motion and an additional adaptation to the energy functional, as well as the requirements to the variational ansatz are discussed in detail. The applicability of the procedure is demonstrated for a cubic model potential for which we calculate thermal decay rates of a frozen Gaussian wavefunction. The decay rate obtained with a narrow trial wavefunction agrees perfectly with the results using the classical normal form of the corresponding point particle. The results with a broader trial wavefunction go even beyond the classical approach, i.e. they agree with those using the quantum normal form. The method presented here will be applied to Bose–Einstein condensates in the following paper (Junginger et al 2012 J. Phys. A: Math. Theor. 45 155202).
We demonstrate the application of transition state theory to wave packet dynamics in metastable Schrödinger systems which are approached by means of a variational ansatz for the wavefunction and whose dynamics is described within the framework of a time-dependent variational principle. The application of classical transition state theory, which requires knowledge of a classical Hamilton function, is made possible by mapping the variational parameters to classical phase space coordinates and constructing an appropriate Hamiltonian in action variables. This mapping, which is performed by a normal form expansion of the equations of motion and an additional adaptation to the energy functional, as well as the requirements to the variational ansatz are discussed in detail. The applicability of the procedure is demonstrated for a cubic model potential for which we calculate thermal decay rates of a frozen Gaussian wavefunction. The decay rate obtained with a narrow trial wavefunction agrees perfectly with the results using the classical normal form of the corresponding point particle. The results with a broader trial wavefunction go even beyond the classical approach, i.e. they agree with those using the quantum normal form. The method presented here will be applied to Bose–Einstein condensates in the following paper (Junginger et al 2012 J. Phys. A: Math. Theor. 45 155202).
We apply transition state theory to coupled Gaussian wave packets and calculate thermal decay rates of Bose–Einstein condensates with additional long-range interaction. The ground state of such a condensate is metastable if the contact interaction is attractive and a sufficient thermal excitation may lead to its collapse. The use of transition state theory is made possible by describing the condensate within a variational framework and locally mapping the variational parameters to classical phase space as has been demonstrated in the preceding paper (Junginger et al 2012 J. Phys. A: Math. Theor. 45 155201). We apply this procedure to Gaussian wave packets and present results for condensates with monopolar 1/ r -interaction comparing decay rates obtained by using different numbers of coupled Gaussian trial wavefunctions as well as different normal form orders.
Bose-Einstein condensates are bosonic many-body systems in which the quantum mechanical ground state is macroscopically occupied. If the interparticle interaction between the single bosons is attractive, this ground state is metastable and separated from the collapse by a barrier. A sufficient thermal excitation of the condensate can enable the crossing of this barrier and induce its collapse. In this thesis, the corresponding decay mechanism of the thermally induced coherent collapse of Bose-Einstein condensates with additional long-range interaction is investigated for the first time. For the description of the condensate’s dynamics a variational approach to the quantum wave function is used, which allows for the precise determination of the system’s transition state that governs the reaction dynamics. In order to quantitatively calculate the decay rate and lifetime of the condensate, respectively, a novel transition state theory for quantum wave packets is developed. This is based on a variational approach as well as normal form expansions of the dynamical equations and the energy functional in the vicinity of the transition state. The dynamical equations obtained from the variational approach induce a symplectic structure onto the space of variational parameters. Furthermore, the normal form expansions serve as a tool to locally construct canonical coordinates on this symplectic manifold. For this purpose, the dynamical equations as well as the energy functional are expanded and diagonalized. The eigenvectors of the linearized equations of motion thereby define a symplectic basis of the system. Afterwards, the transformations of the higher-order terms are carried out in two steps via nonlinear near-identity transformations: In the first step, the nonresonant terms of the expansions are eliminated and in the second step the remaining resonant terms are transformed in such a way that they fulfill canonical equations. The combination of these methods finally allows for the application of the well-established classical transition state theory. The procedure is valid for a wide variety of trial wave functions, for arbitrary dimensions of the variational space, in arbitrary order of the expansions, and it is independent of the precise form of the Hamilton operator. In addition, a procedure to handle bifurcations in the transition state is developed using a uniform rate formula in connection with an appropriate normal form of the potential energy surface. In contrast to the usual treatment which breaks down at the point of the bifurcation, the new method allows to cross the bifurcation smoothly, because it is valid on both sides of the bifurcation as well as arbitrarily close to it. At the same time, one can calculate the flux over the whole configuration of saddles which emerge in the bifurcation. The applicability of the newly developed transition state theory for wave packets is demonstrated for a cubic potential well. The thermal decay rates of a point particle and a frozen Gaussian wave packet, respectively, are calculated within the classical and quantum normal form approaches as well as the variational one. It is shown that the variational approach is capable to perfectly reproduce both the classical and quantum normal form results of the point particle in different limits. In the limit of a narrow trial wave function it agrees perfectly with the results using the classical normal form. The results with a broader trial wave function go even beyond the classical approach, i.e. they coincide with those of the quantum normal form. In the field of Bose-Einstein condensates with long-range interaction the application of transition state theory for quantum wave packets is extended to the field of nonlinear Schrödinger systems. For condensates with attractive 1/r-interaction, the thermal decay rates are calculated for different numbers of coupled Gaussian wave functions as well as different normal form orders. Significant corrections of the decay rate are obtained with increasing order of the normal forms. In the low-temperature regime convergence is observed in eighth order of the expansions. Finally, the thermally induced coherent collapse of dipolar condensates is investigated. These exhibit the feature of a biconcave ground state density distribution for certain values of the physical parameters. It is shown that such condensates reveal the astonishing effect of a symmetry-breaking collapse dynamics. This effect is directly related to the occurrence of the structured ground state, and it depends on the geometry of the external trapping potential. In case of a “conventional” density distribution, the collapse possesses an s-wave symmetry, while it has a d-wave symmetry in case of a blood-cell shaped ground state. The thermally induced coherent collapse plays a major role close to the stability threshold where the energy barrier is small and the decay rates are high. In case of Cr atoms, dipolar condensate have been first realized with, the region of interest is on the order of a Bohr radius above the critical scattering length. There, the lifetime of the condensate can be reduced to a millisecond. For the elements Dy and Er which have also been condensed in the laboratory, the effects become even stronger, i.e. the parameter region of importance becomes larger, the energy barriers become smaller, and the decay rates are higher.
The decay of Bose-Einstein condensates from a metastable ground state into collapse due to macroscopic quantum tunneling is investigated using a semiclassical approximation to Feynman’s path integral formalism. Applying a variational ansatz of a single Gaussian to the wave function determined by the Gross-Pitaevskii equation, a special choice of the Gaussian width parameters yields a mean-field energy functional in Hamiltonian form. The temporal spatial extension of the condensate then is described in the picture of a particle moving in an external potential. In this picture the decay of the condensate wave function is investigated using the bounce trajectory method, accounting for the action of the tunneling orbit in imaginary time and quadratic fluctuations around it, the latter described by the Gelfand-Yaglom differential equation. Additionally, the tunneling formalism is used for describing the condensate wave function by a trial function of a superposition of several Gaussians. Using this ansatz a Hamiltonian form of the mean-field energy does not exist and the formalism of describing the tunneling process has to be extended. The bounce trajectory is described in the parameter space of the Gaussian width parameters and the equations of motion for the imaginary time evolution are obtained by applying a time-dependent variational principle. The action of the bounce trajectory is investigated as a function of the number of Gaussians taken into account and the contribution of the fluctuations is evaluated using the monodromy matrix.
The wave packet dynamics in atomic systems and in Bose-Einstein condensates is investigated by means of a time-dependent variational principle. The wave packets are assumed to be parametrized by a set of time-dependent parameters. The time evolution of the parameters of the coupled Gaussian wave packets can be calculated from a set of ordinary differential equations, as obtained from the time-dependent variational principle. Unfortunately, the set of equations is ill-behaved in most practical applications, depending on the number of propagated Gaussian wave packets, and methods for regularization are needed. A general method for regularization based on applying adequate nonholonomic inequality constraints to the evolution of the parameters, keeping the equations of motion well-behaved is presented. The power of the method is demonstrated for a non-integrable system with two degrees of freedom. The Gaussian wave packet (GWP) method is applied to the three-dimensional hydrogen atom. The regularization based on the introduction of Kustaanheimo-Stiefel coordinates and a fictitious time is performed. The regularization implies a transition from the three-dimensional physical position space to a four-dimensional parameter space, spanned by the Kustaanheimo-Stiefel coordinates. The regularization is accompanied by a restriction on physically allowed wave functions. The Coulomb potential is transformed to a harmonic potential and GWPs are the exact solutions, provided they fulfill the restriction. The effect of the restriction on the four-dimensional GWP is discussed and it is shown that the GWPs can satisfy the restriction if the Gaussian parameter space is reduced in a certain way. The exact analytic evolution of the restricted GWP is presented, and the expansion of a localized initial wave function in the restricted Gaussian basis set and its analytic propagation in the fictitious time are shown. Symmetry subspaces with conserved magnetic quantum number m and with conserved angular momentum quantum numbers l,m are treated separately. The method is also applied to the non-integrable H atom in in a homogeneous magnetic field and in perpendicular external electric and magnetic fields. The evolution of the wave packets is determined by the constrained time-dependent variational principle. The numerical results are compared to numerically exact values and show excellent agreement. Another class of systems where the variational Gaussian wave packet method yields good results are cold gases. In this thesis Bose-Einstein condensates are investigated, where in addition to the common short-range contact interaction two different long-range interactions, viz. a laser induced gravity-like 1/r-interaction or a magnetic dipole-dipole interaction are present. The dynamics as resulting from the time-dependent extended Gross-Pitaevskii equation for Bose-Einstein condensates with attractive 1/r-interaction is investigated with both the GWP method and numerically exact calculations. It is shown that these condensates exhibit signatures known from the nonlinear dynamics of autonomous Hamiltonian systems. The two stationary solutions created in a tangent bifurcation at a critical value of the scattering length are identified as elliptical and hyperbolical fixed points, corresponding to stable and unstable stationary states of the condensate. The stable stationary state is surrounded by elliptical islands, corresponding to condensates periodically oscillating in time, whereas condensate wave functions in the unstable region undergo a collapse within finite time. For negative scattering lengths below the tangent bifurcation no stationary solutions exist, i.e., the condensate is always unstable and collapses. The dynamics of condensates with inter-atomic magnetic dipole-dipole interaction is investigated in the mean field limit using a Gaussian trial function. The anisotropy of the magnetic dipole-dipole interaction breaks the spherical symmetry and for ordered dipoles an effectively two-dimensional system with cylindrical symmetry is obtained. Special attention is payed to the regularity of the dynamics and it is shown that a transition from regular to chaotic motion takes place with increasing energy where regions of regular and chaotic motion coexist. It is shown that stable modes exist at energies high above the saddle point energy, i.e. at energies where a collapse of the condensate is expected.
This work is about a quantum-mechanical generalisation of transition rate theory that accounts for temperature-dependent quantum-dynamical effects. The starting point of this investigation is a formulation of the classical transmission coefficient where the activated dynamics is needed to obtain the time evolution of the classical state function and the rate is expressed as an average flux in momentum space. Traditional approaches to quantum transition rate theory are briefly discussed, with emphasis on the fact that these do not generally require knowledge of the open quantum dynamics. The latter can be studied by means of the Lindblad equation for mixed states, which is shown to be analogous to an average over realisations of a certain type of stochastic Schrödinger dynamics. The stochastic evolution in Hilbert space is used in combination with an ansatz for the wave function as a Gaussian wave packet to show that under appropriate assumptions the stochastic Schrödinger equation converges to the classical Langevin equation. A novel formulation of the quantum transmission coefficient is presented which features a quantum-mechanical version of the state function obtained via the density operator, so that the activated quantum dynamics can be taken into account by solving the Lindblad equation. A numerical investigation of the quantum transmission coefficient thus obtained for different temperatures is carried out, along with a comparison to the classical results. The quantum-mechanical transmission coefficient is found to display some expected features, such as vanishing in the limit as t → 0, displaying a “plateau” value for certain values of the temperature, and decaying to zero exponentially fast when the fundamental assumptions of TRT do not hold. Additionally the plateau value of the quantum transmission coefficient is found to be smaller than the classical counterpart in all regimes, and to approach the corresponding classical value in the limit of large temperatures.
We propose a variational framework for the efficient simulation of stochastic quantum dynamics in interacting bosonic systems. Our method generalizes the time-dependent variational principle to accommodate quantum state diffusion processes and is tailored for describing quantum trajectories unraveled from a master equation. We adopt a wavefunction ansatz composed of a coherent superposition of Gaussian wavepackets, which allows the fully analytical assembly of the variational equations of motion while capturing both the quantum trajectory dynamics and the Lindblad evolution beyond semiclassical approximations. The method is carefully benchmarked on a Bose-Hubbard dimer, and applied to extended Bose-Hubbard lattices in both 1D and 2D. In particular, we show the emergence of a symmetry-breaking phase transition in 2D lattices, which is absent in 1D chains.