Researcher profile

Jörg Main

Jörg Main contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
32works
0followers
9topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

32 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.

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.

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.