Source author record

Jörg Main

Jörg Main appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

61works
9topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

61 published item(s)

preprint2022arXiv

Mean first-passage times for solvated LiCN isomerization at intermediate to high temperatures

The behavior of a particle in a solvent has been framed using stochastic dynamics since the early theory of Kramers. A particle in a chemical reaction reacts slower in a diluted solvent because of the lack of energy transfer via collisions. The flux-over-population reaction rate constant rises with increasing density before falling again for very dense solvents. This Kramers turnover is observed in this paper at intermediate and high temperatures in the backward reaction of the LiNC $\rightleftharpoons$ LiCN isomerization via Langevin dynamics and mean first-passage times (MFPTs). It is in good agreement with the Pollak-Grabert-Hänggi (PGH) reaction rates at lower temperatures. Furthermore, we find a square root behavior of the reaction rate at high temperatures and have made direct comparisons of the methods in the intermediate- and high- temperature regimes; all suggesting increased ranges in accuracy of both the PGH and MFPT approaches.

preprint2022arXiv

Signatures of exciton orbits in quantum mechanical recurrence spectra of Cu$_2$O

The seminal work by T. Kazimierczuk et al. [Nature 514, 343 (2014)] has shown the existence of highly excited exciton states in a regime, where the correspondence principle is applicable and quantum mechanics turns into classical mechanics, however, any interpretation of exciton spectra based on a classical approach to excitons is still missing. Here, we close this gap by computing and comparing quantum mechanical and semiclassical recurrence spectra of cuprous oxide. We show that the quantum mechanical recurrence spectra exhibit peaks, which, by application of semiclassical theories and a scaling transformation, can be directly related to classical periodic exciton orbits. The application of semiclassical theories to exciton physics requires the detailed analysis of the classical exciton dynamics, including three-dimensional orbits, which strongly deviate from hydrogenlike Keplerian orbits. Our findings illuminate important aspects of excitons in semiconductors by directly relating the quantum mechanical band-structure splittings of excitons to the corresponding classical exciton dynamics.

preprint2022arXiv

Transition state theory characterizes thin film macrospin dynamics driven by an oscillatory magnetic field: Inertial effects

Understanding the magnetization switching process in ferromagnetic thin films is essential for many technological applications. We investigate the effects of periodic driving via magnetic fields on a macrospin system under explicit consideration of inertial dynamics. This is usually achieved by extending the Landau-Lifshitz-Gilbert equation with a term including the second time derivative of the magnetization. The dynamics of the magnetization switching can then be characterized by its switching rate. We apply methods from transition state theory for driven systems to resolve the rate of magnetization switching in this general case. In doing so, we find that magnetization exhibits resonance-like behavior under certain driving conditions, and it can be affected strongly by the system's relaxation rate.

preprint2021arXiv

Dynamics and decay rates of a time-dependent two-saddle system

The framework of transition state theory (TST) provides a powerful way for analyzing the dynamics of physical and chemical reactions. While TST has already been successfully used to obtain reaction rates for systems with a single time-dependent saddle point, multiple driven saddles have proven challenging because of their fractal-like phase space structure. This paper presents the construction of an approximately recrossing-free dividing surface based on the normally hyperbolic invariant manifold in a time-dependent two-saddle model system. Based on this, multiple methods for obtaining instantaneous (time-resolved) decay rates of the underlying activated complex are presented and their results discussed.

preprint2021arXiv

Exchange interaction in the yellow exciton series of cuprous oxide

We experimentally and numerically investigate the exchange interaction of the yellow excitons in cuprous oxide. By varying the material parameters in the numerical calculations, we can interpret experimental findings and understand their origin in the complex band structure and central-cell corrections. In particular, we experimentally observe the reversal of the ortho- and paraexciton for the $2S$ yellow exciton, and explain this phenomenon by an avoided crossing with the green $1S$ orthoexciton in a detailed numerical analysis. Furthermore, we discuss the exchange splitting as a function of the principal quantum number $n$ and its deviation from the $n^{-3}$ behavior expected from a hydrogenlike model. We also explain why the observed exchange splitting of the green $1S$ exciton is more than twice the splitting of the yellow $1S$ state.

preprint2021arXiv

Thermal decay rates of an activated complex in a driven model chemical reaction

Recent work has shown that in a non-thermal, multidimensional system, the trajectories in the activated complex possess different instantaneous and time-averaged reactant decay rates. Under dissipative dynamics, it is known that these trajectories, which are bound on the normally hyperbolic invariant manifold (NHIM), converge to a single trajectory over time. By subjecting these dissipative systems to thermal noise, we find fluctuations in the saddle-bound trajectories and their instantaneous decay rates. Averaging over these instantaneous rates results in the decay rate of the activated complex in a thermal system. We find, that the temperature dependence of the activated complex decay in a thermal system can be linked to the distribution of the phase space resolved decay rates on the NHIM in the non-dissipative case. By adjusting the external driving of the reaction, we show that it is possible to influence how the decay rate of the activated complex changes with rising temperature.

preprint2020arXiv

Balancing gain and loss in symmetrised multi-well potentials

Balanced gain and loss leads to stationary dynamics in open systems. This occurs naturally in PT-symmetric systems, where the imaginary part of the potential describing gain and loss is perfectly antisymmetric. While this case seems intuitive, stationary dynamics are also possible in asymmetric open systems. Open multi-well quantum systems can possess completely or partly real spectra if their Hamiltonian is symmetrised or semi-symmetrised, respectively. In contrast to similar concepts, symmetrisation allows for the description of physical multi-well potentials with gain and loss. A simple matrix model for the description of two and three-mode systems is used as an example, for which analytical symmetrised solutions are derived. It is explicitly shown how symmetrisation can be used to systematically find two-mode systems with a stable, stationary ground state and why only PT-symmetric two-mode systems can have stationary excited states.

preprint2020arXiv

Classical and semiclassical description of Rydberg excitons in cuprous oxide

Experimental and theoretical investigations of excitons in cuprous oxide have revealed a significant fine-structure splitting of the excitonic Rydberg states caused by a strong impact of the valence band structure. We provide a semiclassical interpretation of that splitting by investigating the classical dynamics of the excitonic electron-hole pair beyond the hydrogen-like model. Considering the slow motion of Rydberg excitons in coordinate space compared to the fast dynamics of quasispin and hole spin we use an adiabatic approach and energy surfaces in momentum space for the computation of the exciton dynamics. We observe quasi-periodic motion on near-integrable tori. Semiclassical torus quantization yields the energy regions of the fine-structure splitting of $n$-manifolds in agreement with quantum mechanical computations.

preprint2020arXiv

Green exciton series in cuprous oxide

We numerically investigate the odd parity states of the green exciton series in cuprous oxide. Taking into account the coupling to the yellow series and especially to the yellow continuum, the green excitons are quasi-bound resonances with a finite lifetime which cannot be described with Hermitian operators. To calculate their positions and linewidths, we use the method of complex-coordinate rotation, leading to a non-Hermitian complex eigenvalue problem. We find that the behavior of the dominant P states is very well approximated by a modified Rydberg formula using a negative quantum defect. The corresponding linewidths induced by the coupling to the yellow continuum decrease with the third power of the principal quantum number.

preprint2020arXiv

Neural network approach for the dynamics on the normally hyperbolic invariant manifold of periodically driven systems

Chemical reactions in multidimensional systems are often described by a rank-1 saddle, whose stable and unstable manifolds intersect in the normally hyperbolic invariant manifold (NHIM). Trajectories started on the NHIM in principle never leave this manifold when propagated forward or backward in time. However, the numerical investigation of the dynamics on the NHIM is difficult because of the instability of the motion. We apply a neural network to describe time-dependent NHIMs and use this network to stabilize the motion on the NHIM for a periodically driven model system with two degrees of freedom. The method allows us to analyze the dynamics on the NHIM via Poincaré surfaces of section (PSOS) and to determine the transition state (TS) trajectory as a periodic orbit with the same periodicity as the driving saddle, viz. a fixed point of the PSOS surrounded by near-integrable tori. Based on Transition State Theory and a Floquet analysis of a periodic TS trajectory we compute the rate constant of the reaction with significantly reduced numerical effort compared to the propagation of a large trajectory ensemble.

preprint2020arXiv

Phase-space resolved rates in driven multidimensional chemical reactions

Chemical reactions in multidimensional driven systems are typically described by a time-dependent rank-1 saddle associated with one reaction and several orthogonal coordinates (including the solvent bath). To investigate reactions in such systems, we develop a fast and robust method ---viz., local manifold analysis (LMA)--- for computing the instantaneous decay rate of reactants. Specifically, it computes the instantaneous decay rates along saddle-bound trajectories near the activated complex by exploiting local properties of the stable and unstable manifold associated with the normally hyperbolic invariant manifold (NHIM). The LMA method offers substantial reduction of numerical effort and increased reliability in comparison to direct ensemble integration. It provides an instantaneous flux that can be assigned to every point on the NHIM and which is associated with a trajectory ---regardless of whether it is periodic, quasi-periodic, or chaotic--- that is bound on the NHIM. The time average of these fluxes in the driven system corresponds to the average rate through a given local section containing the corresponding point on the NHIM. We find good agreement between the results of the LMA and direct ensemble integration obtained using numerically constructed, recrossing-free dividing surfaces.

preprint2020arXiv

Second Harmonic Generation of cuprous oxide in magnetic fields

Recently Second Harmonic Generation (SHG) for the yellow exciton series in cuprous oxide has been demonstrated [J. Mund et al., Phys. Rev. B 98, 085203 (2018)]. Assuming perfect $O_{\mathrm{h}}$ symmetry, SHG is forbidden along certain high-symmetry axes. Perturbations can break this symmetry and forbidden transitions may become allowed. We investigate theoretically the effect of external magnetic fields on the yellow exciton lines of cuprous oxide. We identify two mechanisms by which an applied magnetic field can induce a second harmonic signal in a forbidden direction. First of all, a magnetic field by itself generally lifts the selection rules. In the Voigt configuration, an additional magneto-Stark electric field appears. This also induces certain SHG processes differing from those induced by the magnetic field alone. Complementary to the manuscript by A. Farenbruch et al. [Phys. Rev. B, submitted], we perform a full numerical diagonalization of the exciton Hamiltonian including the complex valence band structure. Numerical results are compared with experimental data.

preprint2019arXiv

Balanced gain and loss in spatially extended non-$\mathcal{PT}$-symmetric multi-well potentials

The experimental realization of balanced gain and loss in a quantum system has been a long standing goal in quantum mechanics since the introduction of the concept of $\mathcal{PT}$ symmetry and has only recently been achieved. In this paper we analyze balanced gain and loss in Gaussian multi-well potentials with either only gain or loss in each well. By means of symmetrization via matrix models we can construct asymmetric extended potentials with partially real or complex conjugate spectra. This will be demonstrated explicitly for double-well and triple-well systems. Such systems can be realized with Bose-Einstein condensates in optical trapping potentials in the presence of localized particle gain and loss. The usage of asymmetric potentials in the process is more versatile and is considered beneficial in real experimental implementations.

preprint2019arXiv

Rydberg excitons in electric and magnetic fields obtained with the complex-coordinate-rotation method

The complete theoretical description of experimentally observed magnetoexcitons in cuprous oxide has been achieved by F. Schweiner et al [Phys. Rev. B 95, 035202 (2017)], using a complete basis set and taking into account the valence band structure and the cubic symmetry of the solid. Here, we extend these calculations by investigating numerically the autoionising resonances of cuprous oxide in electric fields and in parallel electric and magnetic fields oriented in [001] direction. To this aim we apply the complex-coordinate-rotation method. Complex resonance energies are computed by solving a non-Hermitian generalised eigenvalue problem, and absorption spectra are simulated by using relative oscillator strengths. The method allows us to investigate the influence of different electric and magnetic field strengths on the position, the lifetime, and the shape of resonances.

preprint2016arXiv

$K$ dependent exchange interaction of the $1S$ ortho exciton in Cu$_2$O

When treating the exchange interaction of Wannier excitons, usually only the leading terms of the analytic and the nonanalytic exchange interaction are considered. However, higher order terms can lead to a splitting of exciton states, for which reason a splitting of the $1S$~exciton in cuprous oxide $\left(\mathrm{Cu_{2}O}\right)$ depending on its total momentum $\hbar K$ has been attributed to a $K$ dependent analytic exchange interaction by Dasbach et al [Phys. Rev. Lett. 91, 107401 (2003)]. Going beyond the common treatment of the exchange interaction, we derive the correct expressions for these $K$ dependent higher order terms using $\boldsymbol{k}\cdot\boldsymbol{p}$ perturbation theory. We prove that the appearance of a $K$ dependent exchange interaction is inseparably connected with a $K$ independent exchange interaction of $P$ and $D$ excitons. We estimate the magnitude of these terms for $\mathrm{Cu_{2}O}$ from microscopic calculations and show that they are far too small to explain the observed $K$ dependent splitting. Instead, this splitting has to be treated in terms of the dispersion of the excitons. Furthermore, we prove the occurence of a coupling between longitudinal and transverse excitons in $\mathrm{Cu_{2}O}$ due to the $K$ dependent nonanalytic exchange interaction.

preprint2016arXiv

Bifurcations and exceptional points in a PT-symmetric dipolar Bose-Einstein condensate

We investigate the bifurcation structure of stationary states in a dipolar Bose-Einstein condensate located in an external PT-symmetric potential. The imaginary part of this external potential allows for the effective description of in- and out-coupling of particles. To unveil the complete bifurcation structure and the properties of the exceptional points we perform an analytical continuation of the Gross-Pitaevskii equation, which is used to describe the system. We use an elegant and numerically efficient method for the analytical continuation of the Gross-Pitaevskii equation with dipolar interactions containing bicomplex numbers. The Bose-Einstein condensate with dipole interaction shows a much richer bifurcation scenario than a condensate without long-range interactions. The inclusion of analytically continued states can also explain property changes of bifurcation points which were hidden before, and allows for the examination of the properties of the exceptional points associated with the branch points. With the new analytically continued states we are able to prove the existence of an exceptional point of fifth order.

preprint2016arXiv

Bose-Einstein condensates with balanced gain and loss beyond mean-field theory

Most of the work done in the field of Bose-Einstein condensates with balanced gain and loss has been performed in the mean-field approximation using the PT-symmetric Gross-Pitaevskii equation. In this work we study the many-particle dynamics of a two-mode condensate with balanced gain and loss described by a master equation in Lindblad form whose purity periodically drops to small values but then is nearly completely restored. This effect cannot be covered by the mean-field approximation, in which a completely pure condensate is assumed. We present analytic solutions for the dynamics in the non-interacting limit and use the Bogoliubov backreaction method to discuss the influence of the on-site interaction. Our main result is that the strength of the purity revivals is almost exclusively determined by the strength of the gain and loss and is independent of the amount of particles in the system and the interaction strength. For larger particle numbers, however, strong revivals are shifted towards longer times, but by increasing the interaction strength these strong revivals again occur earlier.

preprint2016arXiv

Dynamics of Solitons in the One-Dimensional Nonlinear Schrödinger Equation

We investigate bright solitons in the one-dimensional Schrödinger equation in the framework of an extended variational approach. We apply the latter to the stationary ground state of the system as well as to coherent collisions between two or more solitons. Using coupled Gaussian trial wave functions, we demonstrate that the variational approach is a powerful method to calculate the soliton dynamics. This method has the advantage that it is computationally faster compared to numerically exact grid calculations. In addition, it goes far beyond the capability of analytical ground state solutions, because the variational approach provides the ability to treat excited solitons as well as dynamical interactions between different wave packets. To demonstrate the power of the variational approach, we calculate the stationary ground state of the soliton and compare it with the analytical solution showing the convergence to the exact solution. Furthermore, we extend our calculations to nonstationary solitons by investigating the breathing oscillations of excited solitons and the coherent collisions of several wave packets in both the low- and high-energy regime. Comparisons of the variational approach with numerically exact simulations on grids reveal excellent agreement in the high-energy regime while deviations can be observed for low energies.

preprint2016arXiv

Impact of the valence band structure of Cu$_{2}$O on excitonic spectra

We present a method to calculate the exciton spectra of all direct semiconductors with a complex valence band structure. The Schrödinger equation is solved using a complete basis set with Coulomb Sturmian functions. This method also allows for the computation of oscillator strengths. Here we apply this method to investigate the impact of the valence band structure of cuprous oxide $\left(\mathrm{Cu_{2}O}\right)$ on the yellow exciton spectrum. Results differ from those in [Phys. Rev. Lett. 115, 027402 (2015)]; the differences are discussed and explained. The difference between the second and third Luttinger parameter can be determined by comparisons with experiments, however, the evaluation of all three Luttinger parameters is not uniquely possible. Our results are consistent with band structure calculations. Considering also a finite momentum $\hbar K$ of the center of mass, we show that the large $K$-dependent line splitting observed for the $1S$ exciton state [Phys. Rev. Lett. 91, 107401 (2003)] is not related to an exchange interaction but rather to the complex valence band structure of $\mathrm{Cu_{2}O}$.

preprint2016arXiv

Linewidths in excitonic absorption spectra of cuprous oxide

We present a theoretical calculation of the absorption spectrum of cuprous oxide $\left(\mathrm{Cu_{2}O}\right)$ based on the general theory developed by Y. Toyozawa. An inclusion not only of acoustic phonons but also of optical phonons and of specific properties of the excitons in $\mathrm{Cu_{2}O}$ like the central-cell corrections for the $1S$-exciton allows us to calculate the experimentally observed line widths in experiments by T. Kazimierczuk et al [Nature 514, 343, (2014)] within the same order of magnitude, which demonstrates a clear improvement in comparison to earlier work on this topic. We also discuss a variety of further effects, which explain the still observable discrepancy between theory and experiment but can hardly be included in theoretical calculations.

preprint2016arXiv

Magnetoexcitons in cuprous oxide

Two of the most striking experimental findings when investigating exciton spectra in cuprous oxide using high-resolution spectroscopy are the observability and the fine structure splitting of $F$ excitons reported by J. Thewes et al. [Phys. Rev. Lett. 115, 027402 (2015)]. These findings show that it is indispensable to account for the complex valence band structure and the cubic symmetry of the solid in the theory of excitons. This is all the more important for magnetoexcitons, where the external magnetic field reduces the symmetry of the system even further. We present the theory of excitons in $\mathrm{Cu_{2}O}$ in an external magnetic field and especially discuss the dependence of the spectra on the direction of the external magnetic field, which cannot be understood from a simple hydrogen-like model. Using high-resolution spectroscopy, we also present the corresponding experimental spectra for cuprous oxide in Faraday configuration. The theoretical results and experimental spectra are in excellent agreement as regards not only the energies but also the relative oscillator strengths. Furthermore, this comparison allows for the determination of the fourth Luttinger parameter $κ$ of this semiconductor.

preprint2016arXiv

Rydberg systems in parallel electric and magnetic fields: an improved method for finding exceptional points

Exceptional points are special parameter points in spectra of open quantum systems, at which resonance energies degenerate and the associated eigenvectors coalesce. Typical examples are Rydberg systems in parallel electric and magnetic fields, for which we solve the Schrödinger equation in a complete basis to calculate the resonances and eigenvectors. Starting from an avoided crossing within the parameter-dependent spectra and using a two-dimensional matrix model, we develop an iterative algorithm to calculate the field strengths and resonance energies of exceptional points and to verify their basic properties. Additionally, we are able to visualise the wave functions of the degenerate states. We report the existence of various exceptional points. For the hydrogen atom these points are in an experimentally inaccessible regime of field strengths. However, excitons in cuprous oxide in parallel electric and magnetic fields, i. e., the corresponding hydrogen analogue in a solid state body, provide a suitable system, where the high-field regime can be reached at much smaller external fields and for which we propose an experiment to detect exceptional points.

preprint2015arXiv

Cusp bifurcation in the eigenvalue spectrum of PT-symmetric Bose-Einstein condensates

A Bose-Einstein condensate in a double-well potential features stationary solutions even for attractive contact interaction as long as the particle number and therefore the interaction strength do not exceed a certain limit. Introducing balanced gain and loss into such a system drastically changes the bifurcation scenario at which these states are created. Instead of two tangent bifurcations at which the symmetric and antisymmetric states emerge, one tangent bifurcation between two formerly independent branches arises [D. Haag et al., Phys. Rev. A 89, 023601 (2014)]. We study this transition in detail using a bicomplex formulation of the time-dependent variational principle and find that in fact there are three tangent bifurcations for very small gain-loss contributions which coalesce in a cusp bifurcation.

preprint2015arXiv

Fano resonances in scattering: an alternative perspective

In a previous paper it has been shown that the interference of the first and second order pole of the Green's function at an exceptional point, as well as the interference of the first order poles in the vicinity of the exceptional point, gives rise to asymmetric scattering cross section profiles. In the present paper we demonstrate that these line profiles are indeed well described by the Beutler-Fano formula, and thus are genuine Fano resonances. Also further away from the exceptional points excellent agreement can be found by introducing energy dependent Fano parameters.

preprint2015arXiv

Normal form expansions and thermal decay rates of Bose-Einstein condensates with short- and long-range interaction

The thermally induced coherent collapse of Bose-Einstein condensates at finite temperature is the dominant decay mechanism near the critical scattering length in condensates with at least partially attractive interaction. The collapse dynamics out of the ground state is mediated by a transition state whose properties determine the corresponding decay rate or lifetime of the condensate. In this paper, we perform normal form expansions of the ground and the transition state of condensates with short-range scattering interaction as well as with anisotropic and long-range dipolar interaction in a variational framework. This method allows one to determine the local properties of these states, i.e. their mean-field energy, their normal modes, the coupling between different modes, and the structure of the reaction channel to any desired order. We discuss the physical interpretation of the transition state as a certain density distribution of the atomic cloud and the behavior of the single normal form contributions in dependence on the s-wave scattering length. Moreover, we investigate the convergence of the local normal form when using extended Gaussian variational approaches, and present the condensate's decay rate.

preprint2015arXiv

Realizing PT-symmetric BEC subsystems in closed hermitian systems

In open double-well Bose-Einstein condensate systems which balance in- and outfluxes of atoms and which are effectively described by a non-hermitian PT-symmetric Hamiltonian PT-symmetric states have been shown to exist. PT-symmetric states obey parity and time reversal symmetry. We tackle the question of how the in- and outfluxes can be realized and introduce a hermitian system in which two PT-symmetric subsystems are embedded. This system no longer requires an in- and outcoupling to and from the environment. We show that the subsystems still have PT-symmetric states. In addition we examine what degree of detail is necessary to correctly model the PT-symmetric properties and the bifurcation structure of such a system. We examine a four-mode matrix model and a system described by the full Gross-Pitaevskii equation in one dimension. We see that a simple matrix model correctly describes the qualitative properties of the system. For sufficiently isolated wells there is also quantitative agreement with the more advanced system descriptions. We also investigate which properties the wave functions of a system must fulfil to allow for PT-symmetric states. In particular the requirements for the phase difference between different parts of the system are examined.

preprint2015arXiv

State Flip at Exceptional Points in Atomic Spectra

We study the behavior of the non-adiabatic population transfer between resonances at an exceptional point in the spectrum of the hydrogen atom. It is known that, when the exceptional point is encircled, the system always ends up in the same state, independent of the initial occupation within the two-dimensional subspace spanned by the states coalescing at the exceptional point. We verify this behavior for a realistic quantum system, viz. the hydrogen atom in crossed electric and magnetic fields. It is also shown that the non-adiabatic hypothesis can be violated when resonances in the vicinity are taken into account. In addition, we study the non-adiabatic population transfer in the case of a third-order exceptional point, in which three resonances are involved.

preprint2015arXiv

Tilted optical lattices with defects as realizations of PT symmetry in Bose-Einstein condensates

A PT-symmetric Bose-Einstein condensate can theoretically be described using a complex optical potential, however, the experimental realization of such an optical potential describing the coherent in- and outcoupling of particles is a nontrivial task. We propose an experiment for a quantum mechanical realization of a PT-symmetric system, where the PT-symmetric currents are implemented by an accelerating Bose-Einstein condensate in a titled optical lattice. A defect consisting of two wells at the same energy level then acts as a PT-symmetric double-well if the tilt in the energy offsets of all further wells in the lattice is varied in time. We map the time-dependence of the amplitudes of a frozen Gaussian variational ansatz to a matrix model and increase the system size step by step starting with a six-well setup. In terms of this simple matrix model we derive conditions under which two wells of the Hermitian multi-well system behave exactly as the two wells of the PT-symmetric system.

preprint2014arXiv

Classical dynamics and localization of resonances in the high energy region of the hydrogen atom in crossed fields

When superimposing the potentials of external fields on the Coulomb potential of the hydrogen atom a saddle point appears, which is called the Stark saddle point. For energies slightly above the saddle point energy one can find classical orbits, which are located in the vicinity of this point. We follow those so-called quasi-Penning orbits to high energies and field strengths observing structural changes and uncovering their bifurcation behavior. By plotting the stability behavior of those orbits against energy and field strength the appearance of a stability apex is reported. A cusp bifurcation, located in the vicinity of the apex, will be investigated in detail. In this cusp bifurcation another orbit of similar shape is found, which becomes completely stable in the observed region of positive energy, i.e., in a region of parameter space, where the Kepler-like orbits located around the nucleus are already unstable. By quantum-mechanically exact calculations we prove the existence of signatures in quantum spectra belonging to those orbits. Husimi distributions are used to compare quantum-Poincaré sections with the extension of the classical torus structure around the orbits. Since periodic orbit theory predicts that each classical periodic orbit contributes an oscillating term to photoabsorption spectra, we finally give an estimation for future experiments, which could verify the existence of the stable orbits.

preprint2014arXiv

Coupling approach for the realization of a $\mathcal{PT}$-symmetric potential for a Bose-Einstein condensate in a double well

We show how non-Hermitian potentials used to describe probability gain and loss in effective theories of open quantum systems can be achieved by a coupling of the system to an environment. We do this by coupling a Bose-Einstein condensate (BEC) trapped in an attractive double-delta potential to a condensate fraction outside the double well. We investigate which requirements have to be imposed on possible environments with a linear coupling to the system. This information is used to determine an environment required for stationary states of the BEC. To investigate the stability of the system we use fully numerical simulations of the dynamics. It turns out that the approach is viable and possible setups for the realization of a $\mathcal{PT}$-symmetric potential for a BEC are accessible. Vulnerabilities of the whole system to small perturbations can be adhered to the singular character of the simplified delta-shaped potential used in our model.

preprint2014arXiv

Dipolar Bose-Einstein condensates in a PT-symmetric double-well potential

We investigate dipolar Bose-Einstein condensates in a complex external double-well potential that features a combined parity and time-reversal symmetry. On the basis of the Gross-Pitaevskii equation we study the effects of the long-ranged anisotropic dipole-dipole interaction on ground and excited states by the use of a time-dependent variational approach. We show that the property of a similar non-dipolar condensate to possess real energy eigenvalues in certain parameter ranges is preserved despite the inclusion of this nonlinear interaction. Furthermore, we present states that break the PT symmetry and investigate the stability of the distinct stationary solutions. In our dynamical simulations we reveal a complex stabilization mechanism for PT-symmetric, as well as for PT-broken states which are, in principle, unstable with respect to small perturbations.

preprint2014arXiv

Harmonic inversion analysis of exceptional points in resonance spectra

The spectra of, e.g. open quantum systems are typically given as the superposition of resonances with a Lorentzian line shape, where each resonance is related to a simple pole in the complex energy domain. However, at exceptional points two or more resonances are degenerate and the resulting non-Lorentzian line shapes are related to higher order poles in the complex energy domain. In the Fourier-transform time domain an $n$-th order exceptional point is characterised by a non-exponentially decaying time signal given as the product of an exponential function and a polynomial of degree $n-1$. The complex positions and amplitudes of the non-degenerate resonances can be determined with high accuracy by application of the nonlinear harmonic inversion method to the real-valued resonance spectra. We extend the harmonic inversion method to include the analysis of exceptional points. The technique yields, in the energy domain, the amplitudes of the higher order poles contributing to the spectra, and, in the time domain, the coefficients of the polynomial characterising the non-exponential decay of the time signal. The extended harmonic inversion method is demonstrated on two examples, viz. the analysis of exceptional points in resonance spectra of the hydrogen atom in crossed magnetic and electric fields, and an exceptional point occurring in the dynamics of a single particle in a time-dependent harmonic trap.

preprint2014arXiv

Nonlinear quantum dynamics in a PT-symmetric double well

We investigate the mean-field dynamics of a Bose-Einstein condensate (BEC) described by the Gross-Pitaevskii equation (GPE) in a double-well potential with particle gain and loss, rendering the system PT-symmetric. The stationary solutions of the system show a change from elliptically stable behavior to hyperbolically unstable behavior caused by the appearance of PT-broken solutions of the GPE and influenced by the nonlinear interaction. The dynamical behavior is visualized using the Bloch sphere formalism. However, the dynamics is not restricted to the surface of the sphere due to the nonlinear and non-Hermitian nature of the system.

preprint2014arXiv

Realizing PT-symmetric non-Hermiticity with ultra-cold atoms and Hermitian multi-well potentials

We discuss the possibility of realizing a non-Hermitian, i.e. an open two-well system of ultra-cold atoms by enclosing it with additional time-dependent wells that serve as particle reservoirs. With the appropriate design of the additional wells PT-symmetric currents can be induced to and from the inner wells, which support stable solutions. We show that interaction in the mean-field limit does not destroy this property. As a first method we use a simplified variational ansatz leading to a discrete nonlinear Schrödinger equation. A more accurate and more general variational ansatz is then used to confirm the results.

preprint2014arXiv

Stability of Bose-Einstein condensates in a $\mathcal{PT}$ symmetric double-$δ$ potential close to branch points

A Bose-Einstein condensate trapped in a double-well potential, where atoms are incoupled to one side and extracted from the other, can in the mean-field limit be described by the nonlinear Gross-Pitaevskii equation (GPE) with a $\mathcal{PT}$ symmetric external potential. If the strength of the in- and outcoupling is increased two $\mathcal{PT}$ broken states bifurcate from the $\mathcal{PT}$ symmetric ground state. At this bifurcation point a stability change of the ground state is expected. However, it is observed that this stability change does not occur exactly at the bifurcation but at a slightly different strength of the in-/outcoupling effect. We investigate a Bose-Einstein condensate in a $\mathcal{PT}$ symmetric double-$δ$ potential and calculate the stationary states. The ground state's stability is analysed by means of the Bogoliubov-de Gennes equations and it is shown that the difference in the strength of the in-/outcoupling between the bifurcation and the stability change can be completely explained by the norm-dependency of the nonlinear term in the Gross-Pitaevskii equation.

preprint2014arXiv

Verification of exceptional points in the collapse dynamics of Bose-Einstein condensates

In Bose-Einstein condensates with an attractive contact interaction the stable ground state and an unstable excited state emerge in a tangent bifurcation at a critical value of the scattering length. At the bifurcation point both the energies and the wave functions of the two states coalesce, which is the characteristic of an exceptional point. In numerical simulations signatures of the exceptional point can be observed by encircling the bifurcation point in the complex extended space of the scattering length, however, this method cannot be applied in an experiment. Here we show in which way the exceptional point effects the collapse dynamics of the Bose-Einstein condensate. The harmonic inversion analysis of the time signal given as the spatial extension of the collapsing condensate wave function can provide clear evidence for the existence of an exceptional point. This method can be used for an experimental verification of exceptional points in Bose-Einstein condensates.

preprint2013arXiv

Bifurcations and exceptional points in dipolar Bose-Einstein condensates

Bose-Einstein condensates are described in a mean-field approach by the nonlinear Gross-Pitaevskii equation and exhibit phenomena of nonlinear dynamics. The eigenstates can undergo bifurcations in such a way that two or more eigenvalues and the corresponding wave functions coalesce at critical values of external parameters. E.g. in condensates without long-range interactions a stable and an unstable state are created in a tangent bifurcation when the scattering length of the contact interaction is varied. At the critical point the coalescing states show the properties of an exceptional point. In dipolar condensates fingerprints of a pitchfork bifurcation have been discovered by Rau et al. [Phys. Rev. A, 81:031605(R), 2010]. We present a method to uncover all states participating in a pitchfork bifurcation, and investigate in detail the signatures of exceptional points related to bifurcations in dipolar condensates. For the perturbation by two parameters, viz. the scattering length and a parameter breaking the cylindrical symmetry of the harmonic trap, two cases leading to different characteristic eigenvalue and eigenvector patterns under cyclic variations of the parameters need to be distinguished. The observed structures resemble those of three coalescing eigenfunctions obtained by Demange and Graefe [J. Phys. A, 45:025303, 2012] using perturbation theory for non-Hermitian operators in a linear model. Furthermore, the splitting of the exceptional point under symmetry breaking in either two or three branching singularities is examined. Characteristic features are observed when one, two, or three exceptional points are encircled simultaneously.

preprint2013arXiv

Dipolar Bose-Einstein condensates in triple-well potentials

Dipolar Bose-Einstein condensates in triple-well potentials are well-suited model systems for periodic optical potentials with important contributions of the non-local and anisotropic dipole-dipole interaction, which show a variety of effects such as self-organisation and formation of patterns. We address here a macroscopic sample of dipolar bosons in the mean-field limit. This work is based on the Gross-Pitaevskii description of dipolar condensates in triple-well potentials by Peter et al 2012 J. Phys. B: At. Mol. Opt. Phys. 45 225302. Our analysis goes beyond the calculation of ground states presented there and clarifies the role of excited and metastable states in such systems. In particular, we find the formation of phases originating from the interplay of several states with distinct stability properties. As some of the phases are formed by metastable states special attention is paid to the characteristics of phase transitions in real-time and the dynamical stabilisation of the condensate.

preprint2013arXiv

Eigenvalue structure of a Bose-Einstein condensate in a PT-symmetric double well

We study a Bose-Einstein condensate in a PT-symmetric double-well potential where particles are coherently injected in one well and removed from the other well. In mean-field approximation the condensate is described by the Gross-Pitaevskii equation thus falling into the category of nonlinear non-Hermitian quantum systems. After extending the concept of PT symmetry to such systems, we apply an analytic continuation to the Gross-Pitaevskii equation from complex to bicomplex numbers and show a thorough numerical investigation of the four-dimensional bicomplex eigenvalue spectrum. The continuation introduces additional symmetries to the system which are confirmed by the numerical calculations and furthermore allows us to analyze the bifurcation scenarios and exceptional points of the system. We present a linear matrix model and show the excellent agreement with our numerical results. The matrix model includes both exceptional points found in the double-well potential, namely an EP2 at the tangent bifurcation and an EP3 at the pitchfork bifurcation. When the two bifurcation points coincide the matrix model possesses four degenerate eigenvectors. Close to that point we observe the characteristic features of four interacting modes in both the matrix model and the numerical calculations, which provides clear evidence for the existence of an EP4.

preprint2013arXiv

Hermitian four-well potential as a realization of a PT-symmetric system

A PT-symmetric Bose-Einstein condensate can be theoretically described using a complex optical potential, however, the experimental realization of such an optical potential describing the coherent in- and outcoupling of particles is a nontrivial task. We propose an experiment for a quantum mechanical realization of a PT-symmetric system, where the PT-symmetric currents of a two-well system are implemented by coupling two additional wells to the system, which act as particle reservoirs. In terms of a simple four-mode model we derive conditions under which the two middle wells of the Hermitian four-well system behave exactly as the two wells of the PT-symmetric system. We apply these conditions to calculate stationary solutions and oscillatory dynamics. By means of frozen Gaussian wave packets we relate the Gross-Pitaevskii equation to the four-mode model and give parameters required for the external potential, which provides approximate conditions for a realistic experimental setup.

preprint2013arXiv

Transition states and thermal collapse of dipolar Bose-Einstein condensates

We investigate thermally excited, dipolar Bose-Einstein condensates. Quasi-particle excitations of the atomic cloud cause density fluctuations which can induce the collapse of the condensate if the inter-particle interaction is attractive. Within a variational approach, we identify the collectively excited stationary states of the gas which form transition states on the way to the BEC's collapse. We analyze transition states with different $m$-fold rotational symmetry and identify the one which mediates the collapse. The latter's symmetry depends on the trap aspect ratio of the external trapping potential which determines the shape of the BEC. Moreover, we present the collapse dynamics of the BEC and calculate the corresponding decay rate using transition state theory. We observe that the thermally induced collapse mechanism is important near the critical scattering length, where the lifetime of the condensate can be significantly reduced. Our results are valid for an arbitrary strength of the dipole-dipole interaction. Specific applications are discussed for the elements $^{52}$Cr, $^{164}$Dy and $^{168}$Er with which dipolar BECs have been experimentally realized.

preprint2012arXiv

A variational approach to Bogoliubov excitations and dynamics of dipolar Bose-Einstein condensates

We investigate the stability properties and the dynamics of Bose-Einstein condensates with axial symmetry, especially with dipolar long-range interaction, using both simulations on grids and variational calculations. We present an extended variational ansatz which is applicable for axial symmetry and show that this ansatz can reproduce the lowest eigenfrequencies of the Bogoliubov spectrum, and also the corresponding eigenfunctions. Our variational ansatz is capable of describing the roton instability of pancake-shaped dipolar condensates for arbitrary angular momenta. After investigating the linear regime we apply the ansatz to determine the dynamics and show how the angular collapse is correctly described within the variational framework.

preprint2012arXiv

Collisions of anisotropic two-dimensional bright solitons in dipolar Bose-Einstein condensates

We investigate the coherent collision of anisotropic quasi-two-dimensional bright solitons in dipolar Bose-Einstein condensates. Our analysis is based on the extended Gross-Pitaevskii equation, and we use the split-operator method for the grid calculations and the time-dependent variational principle with an ansatz of coupled Gaussian functions to calculate the time evolution of the ground state. We compare the results of both approaches for collisions where initially the solitons are in the repelling side-by-side configuration and move towards each other with a specific momentum. We change the relative phases of the condensates, and introduce a total angular momentum by shifting the solitons in opposite direction along the polarization axis. Our calculations show that collisions result in breathing-mode-like excitations of the solitons.

preprint2012arXiv

Macroscopic quantum tunneling of Bose-Einstein condensates with long-range interaction

The ground state of Bose-Einstein condensates with attractive particle interaction is metastable. One of the decay mechanisms of the condensate is a collapse by macroscopic quantum tunneling, which can be described by the bounce trajectory as solution of the time-dependent Gross-Pitaevskii equation in imaginary time. For condensates with an electromagnetically induced gravity-like interaction the bounce trajectory is computed with an extended variational approach using coupled Gaussian functions and simulated numerically exact within the mean-field approach on a space-time lattice. It is shown that the variational computations converge very rapidly to the numerically exact result with increasing number of Gaussians. The tunneling rate of the condensate is obtained from the classical action and additional parameters of the bounce trajectory. The converged variational and numerically exact results drastically improve by several orders of magnitude the decay rates obtained previously with a simple variational approach using a single Gaussian-type orbital for the condensate wave function.

preprint2012arXiv

Quantum-classical lifetimes of Rydberg molecules

A remarkable property of Rydberg atoms is the possibility to create molecules formed by one highly excited atom and another atom in the ground state. The first realisation of such a Rydberg molecule has opened an active field of physical investigations, and showed that its basic properties can be described within a simple model regarding the ground state atom as a small perturber that is bound by a low-energy scattering process with the Rydberg electron [C. H. Greene et al, Phys. Rev. Lett., 85: 2458 (2000)]. Besides the good agreement between theory and the experiment concerning the vibrational states of the molecule, the experimental observations yield the astonishing feature that the lifetime of the molecule is clearly reduced as compared to the bare Rydberg atom [B. Butscher et al, J. Phys. B: At. Mol. Opt. Phys. 44, 184004 (2011)]. With focus on this yet unexplained observation, we investigate in this paper the vibrational ground state of the molecule in a quantum-classical framework. We show that the Rydberg wave function is continuously detuned by the presence of the moving ground state atom and that the timescale on which the detuning significantly exceeds the natural linewidth is in good agreement with the observed reduced lifetimes of the Rydberg molecule.

preprint2012arXiv

Quantum-classical model for the formation of Rydberg molecules

A fascinating aspect of Rydberg atoms is their ability to form huge but very weakly bound molecules with a ground state atom, only held together by a scattering process between the latter and the Rydberg electron. Beyond the usual way of creating such molecules by laser excitation from two ground state atoms with a distance of less than the Rydberg radius, we demonstrate that Rydberg molecules can also be formed by capturing a ground state atom which is initially located outside the range of the Rydberg atom when it comes in contact with it. To demonstrate this effect, we investigate the scattering process between the Rydberg electron and the ground state atom within a quantum-classical framework. In this picture capturing results from a dissipative finite-mass correction term in the classical equations of motion. We show that and under which conditions the capturing takes place.

preprint2012arXiv

Relation between the eigenfrequencies of Bogoliubov excitations of Bose-Einstein condensates and the eigenvalues of the Jacobian in a time-dependent variational approach

We study the relation between the eigenfrequencies of the Bogoliubov excitations of Bose-Einstein condensates, and the eigenvalues of the Jacobian stability matrix in a variational approach which maps the Gross-Pitaevskii equation to a system of equations of motion for the variational parameters. We do this for Bose-Einstein condensates with attractive contact interaction in an external trap, and for a simple model of a self-trapped Bose-Einstein condensate with attractive 1/r interaction. The stationary solutions of the Gross-Pitaevskii equation and Bogoliubov excitations are calculated using a finite-difference scheme. The Bogoliubov spectra of the ground and excited state of the self-trapped monopolar condensate exhibits a Rydberg-like structure, which can be explained by means of a quantum defect theory. On the variational side, we treat the problem using an ansatz of time-dependent coupled Gaussians combined with spherical harmonics. We first apply this ansatz to a condensate in an external trap without long-range interaction, and calculate the excitation spectrum with the help of the time-dependent variational principle. Comparing with the full-numerical results, we find a good agreement for the eigenfrequencies of the lowest excitation modes with arbitrary angular momenta. The variational method is then applied to calculate the excitations of the self-trapped monopolar condensates, and the eigenfrequencies of the excitation modes are compared.

preprint2012arXiv

Symmetry-breaking thermally induced collapse of dipolar Bose-Einstein condensates

We investigate a Bose-Einstein condensate with additional long-range dipolar interaction in a cylindrically symmetric trap within a variational framework. Compared to the ground state of this system, little attention has as yet been payed to its unstable excited states. For thermal excitations, however, the latter is of great interest, because it forms the "activated complex" that mediates the collapse of the condensate. For a certain value of the s-wave scatting length our investigations reveal a bifurcation in the transition state, leading to the emergence of two additional and symmetry-breaking excited states. Because these are of lower energy than their symmetric counterpart, we predict the occurrence of a symmetry-breaking thermally induced collapse of dipolar condensates. We show that its occurrence crucially depends on the trap geometry and calculate the thermal decay rates of the system within leading order transition state theory with the help of a uniform rate formula near the rank-2 saddle which allows to smoothly pass the bifurcation.

preprint2012arXiv

Transition state theory for wave packet dynamics. I. Thermal decay in metastable Schrödinger systems

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 wave function 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 wave function. The decay rate obtained with a narrow trial wave function agrees perfectly with the results using the classical normal form of the corresponding point particle. The results with a broader trial wave function 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 [A. Junginger, M. Dorwarth, J. Main, and G. Wunner, submitted to J. Phys. A].

preprint2012arXiv

Transition state theory for wave packet dynamics. II. Thermal decay of Bose-Einstein condensates with long-range interaction

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 [A. Junginger, J. Main, and G. Wunner, submitted to J. Phys. A]. 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 wave functions as well as different normal form orders.

preprint2011arXiv

Variational calculations for anisotropic solitons in dipolar Bose-Einstein condensates

We present variational calculations using a Gaussian trial function to calculate the ground state of the Gross-Pitaevskii equation and to describe the dynamics of the quasi-two-dimensional solitons in dipolar Bose-Einstein condensates. Furthermore we extend the ansatz to a linear superposition of Gaussians improving the results for the ground state to exact agreement with numerical grid calculations using imaginary time and split-operator method. We are able to give boundaries for the scattering length at which stable solitons may be observed in an experiment. By dynamical calculations with coupled Gaussians we are able to describe the rather complex behavior of the thermally excited solitons. The discovery of dynamically stabilized solitons indicates the existence of such BECs at experimentally accessible temperatures.

preprint2010arXiv

Fractal Weyl law for three-dimensional chaotic hard-sphere scattering systems

The fractal Weyl law connects the asymptotic level number with the fractal dimension of the chaotic repeller. We provide the first test for the fractal Weyl law for a three-dimensional open scattering system. For the four-sphere billiard, we investigate the chaotic repeller and discuss the semiclassical quantization of the system by the method of cycle expansion with symmetry decomposition. We test the fractal Weyl law for various symmetry subspaces and sphere-to-sphere separations.

preprint2010arXiv

Pitchfork bifurcations in blood-cell shaped dipolar Bose-Einstein condensates

We demonstrate that the method of coupled Gaussian wave packets is a full-fledged alternative to direct numerical solutions of the Gross-Pitaevskii equation of condensates with electromagnetically induced attractive 1/r interaction, or with dipole-dipole interaction. Moreover, Gaussian wave packets are superior in that they are capable of producing both stable and unstable stationary solutions, and thus of giving access to yet unexplored regions of the space of solutions of the Gross-Pitaevskii equation. We apply the method to clarify the theoretical nature of the collapse mechanism of blood-cell shaped dipolar condensates: On the route to collapse the condensate passes through a pitchfork bifurcation, where the ground state itself turns unstable, before it finally vanishes in a tangent bifurcation.

preprint2010arXiv

Resonance wave functions located at the Stark saddle point

We calculate quantum mechanically exact wave functions of resonances in spectra of the hydrogen atom in crossed external fields and prove the existence of long-lived decaying quantum states localized at the Stark saddle point. A spectrum of ground and excited states reproducing the nodal patterns expected from simple quadratic and cubic expansions of the potential in the vicinity of the saddle point can be identified. The results demonstrate the presence of resonances in the vicinity of the saddle predicted by simple approximations.

preprint2010arXiv

Semiclassical quantization of the diamagnetic hydrogen atom with near action-degenerate periodic-orbit bunches

The existence of periodic orbit bunches is proven for the diamagnetic Kepler problem. Members of each bunch are reconnected differently at self-encounters in phase space but have nearly equal classical action and stability parameters. Orbits can be grouped already on the level of the symbolic dynamics by application of appropriate reconnection rules to the symbolic code in the ternary alphabet. The periodic orbit bunches can significantly improve the efficiency of semiclassical quantization methods for classically chaotic systems, which suffer from the exponential proliferation of orbits. For the diamagnetic hydrogen atom the use of one or few representatives of a periodic orbit bunch in Gutzwiller's trace formula allows for the computation of semiclassical spectra with a classical data set reduced by up to a factor of 20.

preprint2010arXiv

Variational calculations on multilayer stacks of dipolar Bose-Einstein condensates

We investigate a multilayer stack of dipolar Bose-Einstein condensates in terms of a simple Gaussian variational ansatz and demonstrate that this arrangement is characterized by the existence several stationary states. Using a Hamiltonian picture we show that in an excited stack there is a coupled motion of the individual condensates by which they exchange energy. We find that for high excitations the interaction between the single condensates can induce the collapse of one of them. We furthermore demonstrate that one collapse in the stack can force other collapses, too. We discuss the possibility of experimentally observing the coupled motion and the relevance of the variational results found to full numerical investigations.

preprint2010arXiv

Variational methods with coupled Gaussian functions for Bose-Einstein condensates with long-range interactions. I. General concept

The variational method of coupled Gaussian functions is applied to Bose-Einstein condensates with long-range interactions. The time-dependence of the condensate is described by dynamical equations for the variational parameters. We present the method and analytically derive the dynamical equations from the time-dependent Gross-Pitaevskii equation. The stability of the solutions is investigated using methods of nonlinear dynamics. The concept presented in this paper will be applied to Bose-Einstein condensates with monopolar 1/r and dipolar 1/r^3 interaction in the subsequent paper [S. Rau et al., Phys. Rev. A, submitted], where we will present a wealth of new phenomena obtained by using the ansatz with coupled Gaussian functions.

preprint2010arXiv

Variational methods with coupled Gaussian functions for Bose-Einstein condensates with long-range interactions. II. Applications

Bose-Einstein condensates with an attractive 1/r interaction and with dipole-dipole interaction are investigated in the framework of the Gaussian variational ansatz introduced by S. Rau, J. Main, and G. Wunner [Phys. Rev. A, submitted]. We demonstrate that the method of coupled Gaussian wave packets is a full-fledged alternative to direct numerical solutions of the Gross-Pitaevskii equation, or even superior in that coupled Gaussians are capable of producing both, stable and unstable states of the Gross-Pitaevskii equation, and thus of giving access to yet unexplored regions of the space of solutions of the Gross-Pitaevskii equation. As an alternative to numerical solutions of the Bogoliubov-de Gennes equations, the stability of the stationary condensate wave functions is investigated by analyzing the stability properties of the dynamical equations of motion for the Gaussian variational parameters in the local vicinity of the stationary fixed points. For blood-cell-shaped dipolar condensates it is shown that on the route to collapse the condensate passes through a pitchfork bifurcation, where the ground state itself turns unstable, before it finally vanishes in a tangent bifurcation.

preprint2009arXiv

Fictitious-time wave-packet dynamics in atomic systems

Gaussian wave packets (GWPs) are well suited as basis functions to describe the time evolution of arbitrary wave functions in systems with nonsingular smooth potentials. They are less so in atomic systems on account of the singular behavior of the Coulomb potential. We present a time-dependent variational method that makes the use of GWPs possible in the description of propagation of quantum states also in these systems. We do so by a regularization of the Coulomb potential, and by introduction of a fictitious time coordinate in which the evolution of an initial state can be calculated exactly and analytically for a pure Coulomb potential. Therefore in perturbed atomic systems variational approximations only arise from those parts of the potentials which deviate from the Coulomb potential. The method is applied to the hydrogen atom in external magnetic and electric fields. It can be adapted to systems with definite symmetries, and thus allows for a wide range of applications.

preprint2007arXiv

Discovery of exceptional points in the Bose-Einstein condensation of gases with attractive 1/r-interaction

The extended Gross-Pitaevskii equation for the Bose-Einstein condensation of gases with attractive 1/r-interaction has a second solution which is born together with the ground state in a tangent bifurcation. At the bifurcation point both states coalesce, i.e., the energies and the wave functions are identical. We investigate the bifurcation point in the context of exceptional points, a phenomenon known for linear non-Hermitian Hamiltonians. We point out that the mean field energy, the chemical potential, and the wave functions show the same behavior as an exceptional point in a linear, non-symmetric system. The analysis of the analytically continued Gross-Pitaevskii equation reveals complex waves at negative scattering lengths below the tangent bifurcation. These solutions are interpreted as a decay of the condensate caused by an absorbing potential.

preprint2006arXiv

Hydrogen atom in crossed electric and magnetic fields: Phase space topology and torus quantization via periodic orbits

A hierarchical ordering is demonstrated for the periodic orbits in a strongly coupled multidimensional Hamiltonian system, namely the hydrogen atom in crossed electric and magnetic fields. It mirrors the hierarchy of broken resonant tori and thereby allows one to characterize the periodic orbits by a set of winding numbers. With this knowledge, we construct the action variables as functions of the frequency ratios and carry out a semiclassical torus quantization. The semiclassical energy levels thus obtained agree well with exact quantum calculations.