Source author record

Holger Cartarius

Holger Cartarius 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

35works
10topics
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

35 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

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

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

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

Structural information about the $\text{Ar}_6$ cluster with the frozen Gaussian imaginary time propagator

A numerically cheap way to obtain structural information about clusters of rare gas atoms at low temperatures is developed. The semiclassical frozen Gaussian imaginary time propagator is extended such that it can account for the mean values of all inter-atomic distances in the cluster and their variances. To reduce the required numerical effort an approximation for the mean values is developed which preserves the quality of the results offered by the semiclassical ansatz. The method is applied to the $\text{Ar}_6$ cluster. It is found that the cluster dissociates almost in one step to six free atoms when the temperature is increased. Precursors of the dissociation are only observable in the distances of the atoms via the appearance of a second isomer. The process is almost classical. However, the method is able to resolve small differences in the temperatures at which the dissociation takes place and in the mean distances of the bound configuration.

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

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

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

preprint2013arXiv

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

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

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.

preprint2011arXiv

Fingerprints of exceptional points in the survival probability of resonances in atomic spectra

The unique time signature of the survival probability exactly at the exceptional point parameters is studied here for the hydrogen atom in strong static magnetic and electric fields. We show that indeed the survival probability S(t)=|<psi(0)|psi(t)>|^2 decays exactly as |1-a*t|^2 e^(-Gamma_EP*t/hbar) where Gamma_EP is associated with the decay rate at the exceptional point and a is a complex constant depending solely on the initial wave packet that populates exclusively the two almost degenerate states of the non-Hermitian Hamiltonian. This may open the possibility for a first experimental detection of exceptional points in a quantum system.

preprint2011arXiv

First-order corrections to semiclassical Gaussian partition functions for clusters of atoms

Gaussian approximations to the Boltzmann operator have proven themselves in recent years as useful tools for the study of the thermodynamic properties of rare gas clusters. They are, however, not necessarily correct at very low temperatures. In this article we introduce a first-order correction term to the frozen Gaussian imaginary time propagator and apply it to the argon trimer. Our findings show that the correction term provides objective access to the quality of the propagator's results and clearly defines the "best" Gaussian width parameter. The strength of the correction monitored as a function of the temperature indicates that the results of the Gaussian propagator become questionable below a certain temperature. The interesting thermodynamic transition from a bounded trimer to three body dissociation lies in the temperature range for which the Gaussian approximation is predicted to be accurate.

preprint2010arXiv

Imaginary time Gaussian dynamics of the Ar_3 cluster

Semiclassical Gaussian approximations to the Boltzmann operator have become an important tool for the investigation of thermodynamic properties of clusters of atoms at low temperatures. Usually, numerically expensive thawed Gaussian variants are applied. In this article, we introduce a numerically much cheaper frozen Gaussian approximation to the imaginary time propagator with a width matrix especially suited for the dynamics of clusters. The quality of the results is comparable to that of thawed Gaussian methods based on the single-particle ansatz. We apply the method to the argon trimer and investigate the dissociation process of the cluster. The results clearly show a classical-like transition from a bounded moiety to three free particles at a temperature T ~ 20 K, whereas previous studies of the system were not able to resolve this transition. Quantum effects, i.e., differences with the purely classical case manifest themselves in the low-temperature behavior of the mean energy and specific heat as well as in a slight shift of the transition temperature. We also discuss the influence of an artificial confinement of the atoms usually introduced to converge numerical computations. The results show that restrictive confinements often implemented in studies of clusters can influence the thermodynamic properties drastically. This finding may have implications on other studies of atomic clusters.

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

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.