Source author record

Günter Wunner

Günter Wunner 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

60works
12topics
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

60 published item(s)

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.

preprint2016arXiv

$\mathcal{PT}$-symmetric currents of a Bose-Einstein condensate in a triple well

We study the case of $\mathcal{PT}$-symmetric perturbations of Hermitian Hamiltonians with degenerate eigenvalues using the example of a triple-well system. The degeneracy complicates the question, whether or not a stationary current through such a system can be established, i.e.\ whether or not the $\mathcal{PT}$-symmetric states are stable. It is shown that this is only the case for perturbations that do not couple to any of the degenerate states. The physical explanation for the inhibition of stable currents is discussed. However, introducing an on-site interaction restores the capability to support stable currents.

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

Purity oscillations in Bose-Einstein condensates with balanced gain and loss

In this work we present a new generic feature of PT-symmetric Bose-Einstein condensates by studying the many-particle description of a two-mode condensate with balanced gain and loss. This is achieved using a master equation in Lindblad form whose mean-field limit is a PT-symmetric Gross-Pitaevskii equation. It is shown that the purity of the condensate periodically drops to small values but then is nearly completely restored. This has a direct impact on the average contrast in interference experiments which cannot be covered by the mean-field approximation, in which a completely pure condensate is assumed.

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

Mean-field predictions for a dipolar Bose-Einstein condensate with $^{164}$Dy

We study dipolar Bose-Einstein condensates for a realistic set of parameters close to actual experimental setups with dysprosium. Our analysis is based on the extended Gross-Pitaevskii equation, which we solve numerically exact on a three-dimensional grid. We present stability and phase diagrams, and study the expansion dynamics of dipolar condensates. Our calculations show signatures of the dipole-dipole interaction in terms of structured states and a deviation of the well-known inversion of the aspect ratio of the cloud during a time of flight.

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

A Bose-Einstein Condensate with $\mathcal{PT}$-Symmetric Double-Delta Function Loss and Gain in a Harmonic Trap: A Test of Rigorous Estimates

We consider the linear and nonlinear Schrödinger equation for a Bose-Einstein condensate in a harmonic trap with $\cal {PT}$-symmetric double-delta function loss and gain terms. We verify that the conditions for the applicability of a recent proposition by Mityagin and Siegl on singular perturbations of harmonic oscillator type self-adjoint operators are fulfilled. In both the linear and nonlinear case we calculate numerically the shifts of the unperturbed levels with quantum numbers $n$ of up to 89 in dependence on the strength of the non-Hermiticity and compare with rigorous estimates derived by those authors. We confirm that the predicted $1/n^{1/2}$ estimate provides a valid upper bound on the the shrink rate of the numerical eigenvalues. Moreover, we find that a more recent estimate of $\log(n)/n^{3/2}$ is in excellent agreement with the numerical results. With nonlinearity the shrink rates are found to be smaller than without nonlinearity, and the rigorous estimates, derived only for the linear case, are no longer applicable.

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

Quantum master equation with balanced gain and loss

We present a quantum master equation describing a Bose-Einstein condensate with particle loss on one lattice site and particle gain on the other lattice site whose mean-field limit is a non-Hermitian PT-symmetric Gross-Pitaevskii equation. It is shown that the characteristic properties of PT-symmetric systems, such as the existence of stationary states and the phase shift of pulses between two lattice sites, are also found in the many-particle system. Visualizing the dynamics on a Bloch sphere allows us to compare the complete dynamics of the master equation with that of the Gross-Pitaevskii equation. We find that even for a relatively small number of particles the dynamics are in excellent agreement and the master equation with balanced gain and loss is indeed an appropriate many-particle description of a PT-symmetric Bose-Einstein condensate.

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

Supersymmetric model of a Bose-Einstein condensate in a $\mathcal{PT}$-symmetric double-delta trap

The most important properties of a Bose-Einstein condensate subject to balanced gain and loss can be modelled by a Gross-Pitaevskii equation with an external $\mathcal{PT}$-symmetric double-delta potential. We study its linear variant with a supersymmetric extension. It is shown that both in the $\mathcal{PT}$-symmetric as well as in the $\mathcal{PT}$-broken phase arbitrary stationary states can be removed in a supersymmetric partner potential without changing the energy eigenvalues of the other state. The characteristic structure of the singular delta potential in the supersymmetry formalism is discussed, and the applicability of the formalism to the nonlinear Gross-Pitaevskii equation is analysed. In the latter case the formalism could be used to remove $\mathcal{PT}$-broken states introducing an instability to the stationary $\mathcal{PT}$-symmetric states.

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

Circular orbits in the extreme Reissner-Nordstrøm dihole metric

We study the motion of neutral test particles in the gravitational field of two charged black holes described by the extreme Reissner-Nordstrøm dihole metric where the masses and charges of the black holes are chosen such that the gravitational attraction is compensated by the electrostatic repulsion. We investigate circular orbits in the equatorial plane between the two black holes with equal masses as well as the case of circular orbits outside this symmetry plane. We show that the first case reduces to an effective two-body problem with a behavior similar to a system described by the Reissner-Nordstrøm spacetime. The main focus is directed to the second case with circular orbits outside the equatorial plane.

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

Hartree-Fock calculations for the photoionisation of helium and helium-like ions in neutron star magnetic fields

We derive the photoionisation cross section in dipole approximation for many-electron atoms and ions for neutron star magnetic field strengths in the range of 10^7 to 10^9 T. Both bound and continuum states are treated in adiabatic approximation in a self-consistent way. Bound states are calculated by solving the Hartree-Fock-Roothaan equations using finite-element and B-spline techniques while the continuum orbital is calculated by direct integration of the Hartree-Fock equations in the mean-field potential of the remaining bound orbitals. We take into account mass and photon density in the neutron star's atmosphere, finite nuclear mass as well as thermal occupation of the levels. The data may be of importance for the quantitative interpretation of observed x-ray spectra that originate from the thermal emission of isolated neutron stars. They can serve as input for modeling neutron star atmospheres as regards chemical composition, magnetic field strength, temperature, and redshift. Our main focus in this paper lies on helium and helium-like oxygen. These two-electron systems are simple enough to calculate all possible transitions when limiting the quantum numbers and should show all the basic structures and behaviour of other two-electron systems up to iron.

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

Creating two-dimensional bright solitons in dipolar Bose-Einstein condensates

We propose a realistic experimental setup for creating quasi-two-dimensional (2D) bright solitons in dipolar Bose-Einstein condensates (BECs), the existence of which was proposed in Phys. Rev. Lett. 100, 090406 (2008). A challenging feature of the expected solitons is their strong inherent anisotropy, due to the necessary in-plane orientation of the local moments in the dipolar gas. This may be the first chance of making multidimensional matter-wave solitons, as well as solitons featuring the anistropy due to their intrinsic dynamics. Our analysis is based on the extended Gross-Pitaevskii equation, which includes three-body losses and noise in the scattering length, induced by fluctuations of currents inducing the necessary magnetic fields, which are factors crucial to the adequate description of experimental conditions. By means of systematic 3D simulations, we find a ramping scenario for the change of the scattering length and trap frequencies which results in the creation of robust solitons, that readily withstand the concomitant excitation of the condensate.

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

Model of a PT symmetric Bose-Einstein condensate in a delta-functions double well

The observation of PT symmetry in a coupled optical wave guide system that involves a complex refractive index has been demonstrated impressively in the experiment by Rüter el al. (Nat. Phys. 6, 192, 2010). This is, however, only an optical analogue of a quantum system, and it would be highly desirable to observe the manifestation of PT symmetry and the resulting properties also in a real, experimentally accessible, quantum system. Following a suggestion by Klaiman et al. (Phys. Rev. Lett. 101, 080402, 2008), we investigate a PT symmetric arrangement of a Bose-Einstein condensate in a double well potential, where in one well cold atoms are injected while in the other particles are extracted from the condensate. We investigate, in particular, the effects of the nonlinearity in the Gross-Pitaevskii equation on the PT properties of the condensate. To study these effects we analyze a simple one-dimensional model system in which the condensate is placed into two PT symmetric delta-function traps. The analysis will serve as a useful guide for studies of the behaviour of Bose-Einstein condensates in realistic PT symmetric double wells.

preprint2012arXiv

Nonlinear Schrödinger equation for a PT symmetric delta-functions double well

The time-independent nonlinear Schrödinger equation is solved for two attractive delta-function shaped potential wells where an imaginary loss term is added in one well, and a gain term of the same size but with opposite sign in the other. We show that for vanishing nonlinearity the model captures all the features known from studies of PT symmetric optical wave guides, e.g., the coalescence of modes in an exceptional point at a critical value of the loss/gain parameter, and the breaking of PT symmetry beyond. With the nonlinearity present, the equation is a model for a Bose-Einstein condensate with loss and gain in a double well potential. We find that the nonlinear Hamiltonian picks as stationary eigenstates exactly such solutions which render the nonlinear Hamiltonian itself PT symmetric, but observe coalescence and bifurcation scenarios different from those known from linear PT symmetric Hamiltonians.

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

Twin Paradox in de Sitter Spacetime

The "twin paradox" of special relativity offers the possibility to make interstellar flights within a lifetime. For very long journeys with velocities close to the speed of light, however, we have to take into account the expansion of the universe. Inspired by the work of Rindler on hyperbolic motion in curved spacetime, we study the worldline of a uniformly accelerated observer in de Sitter spacetime and the communication between the traveling observer and an observer at rest.

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.

preprint2009arXiv

Phonon instability and self-organized structures in multi-layer stacks of confined dipolar Bose-Einstein condensates in optical lattices

In calculations to date [1,2] of multi-layer stacks of dipolar condensates, created in one-dimensional optical lattices, the condensates have been assumed to be two-dimensional. In a real experiment, however, the condensates do not extend to infinity in the oblate direction, but have to be confined by a trap potential, too. By three-dimensional numerical simulations of this realistic experimental situation we find a crucial dependence of the phonon instability boundary on the number of layers. Moreover, near the boundary of the phonon instability, a variety of structured ground-state wave functions emerges, which may indicate the onset of a roton instability [3,4].

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.