Source author record

Lorenz S. Cederbaum

Lorenz S. Cederbaum 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

53works
9topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

53 published item(s)

preprint2023arXiv

Non-Hermitian Rayleigh-Schroedinger Perturbation Theory

We devise a non-Hermitian Rayleigh-Schroedinger perturbation theory for the single- and the multireference case to tackle both the many-body problem and the decay problem encountered, for example, in the study of electronic resonances in molecules. A complex absorbing potential (CAP) is employed to facilitate a treatment of resonance states that is similar to the well-established bound-state techniques. For the perturbative approach, the full CAP-Schroedinger Hamiltonian, in suitable representation, is partitioned according to the Epstein-Nesbet scheme. The equations we derive in the framework of the single-reference perturbation theory turn out to be identical to those obtained by a time-dependent treatment in Wigner-Weisskopf theory. The multireference perturbation theory is studied for a model problem and is shown to be an efficient and accurate method. Algorithmic aspects of the integration of the perturbation theories into existing ab initio programs are discussed, and the simplicity of their implementation is elucidated.

preprint2021arXiv

Many-body effects in the excitations and dynamics of trapped Bose-Einstein condensates

This review explores the dynamics and the low-energy excitation spectra of Bose-Einstein condensates (BECs) of interacting bosons in external potential traps putting particular emphasis on the emerging many-body effects beyond mean-field descriptions. To do so, methods have to be used that, in principle, can provide numerically exact results for both the dynamics and the excitation spectra in a systematic manner. Numerically exact results for the dynamics are presented employing the well-established multicongurational time-dependent Hartree for bosons (MCTDHB) method. The respective excitation spectra are calculated utilizing the more recently introduced linear-response theory atop it (LR-MCTDHB). The latter theory gives rise to an, in general, non-hermitian eigenvalue problem. The theory and its newly developed implementation are described in detail and benchmarked towards the exactly-solvable harmonic-interaction model. Several applications to BECs in one- and two-dimensional potential traps are discussed. With respect to dynamics, it is shown that both the out-of-equilibrium tunneling dynamics and the dynamics of trapped vortices are of many-body nature. Furthermore, many-body effects in the excitation spectra are presented for BECs in different trap geometries. It is demonstrated that even for essentially-condensed systems, the spectrum of the lowest-in-energy excitations computed at the many-body level can differ substantially from the standard mean-field description. In general, it is shown that bosons carrying angular momentum are more sensitive to many-body effects than bosons without. These effects are present in both the dynamics and the excitation spectrum.

preprint2021arXiv

Polaritonic states of matter in rotating cavity

The interaction of quantum light with matter like that inside a cavity is known to give rise to mixed light-matter states called polaritons. We discuss the impact of rotation of the cavity on the polaritons. It is shown that the number of polaritons increases due to this rotation. The structure of the original polaritons is modified and new ones are induced by the rotation that strongly depend on the angular velocity and the choice of axis of rotation. In molecules the rotation can change the number of light-induced conical intersections and their dimensionality and hence strongly impact their quantum dynamics. General consequences are discussed.

preprint2021arXiv

Suppression of X-Ray-Induced Radiation Damage to Biomolecules in Aqueous Environments by Immediate Intermolecular Decay of Inner-Shell Vacancies

The predominant reason for the damaging power of high-energy radiation is multiple ionization of a molecule, either direct or via the decay of highly excited intermediates, as e.g., in the case of X-ray irradiation. Consequently, the molecule is irreparably damaged by the subsequent fragmentation in a Coulomb explosion. In an aqueous environment, however, it has been observed that irradiated molecules may be saved from fragmentation presumably by charge and energy dissipation mechanisms. Here, we show that the protective effect of the environment sets in even earlier than hitherto expected, namely immediately after single inner-shell ionization. By combining coincidence measurements of the fragmentation of X-ray-irradiated microsolvated pyrimidine molecules with theoretical calculations, we identify direct intermolecular electronic decay as the protective mechanism, outrunning the usually dominant Auger decay. Our results demonstrate that such processes play a key role in charge delocalization and have to be considered in investigations and models on high-energy radiation damage in realistic environments.

preprint2020arXiv

Born-Oppenheimer approximation in optical cavities: from success to breakdown

The coupling of a molecule and a cavity induces non-adiabaticity in the molecule which makes the description of its dynamics complicated. For polyatomic molecules, reduced-dimensional models and the use of the Born-Oppenheimer approximation (BOA) may remedy the situation. It is demonstrated that contrary to expectation, BOA may even fail in a one-dimensional model and generally expected to fail in two- or more-dimensional models due to the appearance of conical intersections induced by the cavity.

preprint2020arXiv

Fragmentation of molecules by virtual photons from remote neighbors

It is shown that a molecule can dissociate by the energy transferred from a remote neighbor. This neighbor can be an excited neutral or ionic atom or molecule. If it is an atom, the transferred energy is, of course, electronic and in the case of molecules it can also be vibrational. Explicit examples are given which demonstrate that the transfer can be highly efficient at distances where there is no bonding between the transmitter and the dissociating molecule.

preprint2020arXiv

Striking light-induced nonadiabatic fingerprint in the low-energy vibronic spectra of polyatomic molecules

Nonadiabaticity, i.e., the effect of mixing electronic states by nuclear motion, is a central phenomenon in molecular science. The strongest nonadiabatic effects arise due to the presence of conical intersections of electronic energy surfaces. These intersections are abundant in polyatomic molecules. Laser light can induce in a controlled manner new conical intersections, called light-induced conical intersections, which lead to strong nonadiabatic effects similar to those of the natural conical intersections. These effects are, however, controllable and may even compete with those of the natural intersections. In this work we show that the standard low-energy vibrational spectrum of the electronic ground state can change dramatically by inducing nonadiabaticity via a light-induced conical intersection. This generic effect is demonstrated for an explicit example by full-dimensional high-level quantum calculations using a pump-probe scheme with a moderate-intensity pump laser and a weak probe laser.

preprint2016arXiv

Electron-correlation driven capture and release in double quantum dots

We recently predicted that the interatomic Coulombic electron capture (ICEC) process, a long-range electron correlation driven capture process, is achievable in gated double quantum dots (DQDs). In ICEC an incoming electron is captured by one QD and the excess energy is used to remove an electron from the neighboring QD. In this work we present systematic full three-dimensional electron dynamics calculations in quasi-one dimensional model potentials that allow for a detailed understanding of the connection between the DQD geometry and the reaction probability for the ICEC process. We derive an effective one-dimensional approach and show that its results compare very well with those obtained using the full three-dimensional calculations. This approach substantially reduces the computation times. The investigation of the electronic structure for various DQD geometries for which the ICEC process can take place clarify the origin of its remarkably high probability in the presence of two-electron resonances.

preprint2016arXiv

Towards controlling the dissociation probability by light-induced conical intersections

Light-induced conical intersections (LICIs) can be formed both by standing or by running laser waves. The position of a LICI is determined by the laser frequency while the laser intensity controls the strength of the nonadiabatic coupling. Recently, it was shown within the LICI framework that linearly chirped laser pulses have an impact on the dissociation dynamics of the $D_{2}^{+}$ molecule (J. Chem. Phys. 143, 014305, (2015); ibid 144, 074309, (2016)). In this work we exploit this finding and perform calculations using chirped laser pulses in which the time dependence of the laser frequency is designed so as to force the LICI to move together with the field-free vibrational wave packet as much as possible. Since nonadiabaticity is strongest in the vicinity of the conical intersection, this is the first step towards controlling the dissociation process via the LICI. Our showcase example is again the $D_{2}^{+}$ molecular ion. To demonstrate the impact of the LICIs on the dynamical properties of diatomics, the total dissociation probabilities and the population of the different vibrational levels after the dissociation process are studied and discussed.

preprint2016arXiv

Tracking the photodissociation probability of D$_2^+$ induced by linearly chirped laser pulses

In the presence of linearly varying frequency chirped laser pulses the photodissociation dynamics of D$_2^+$ is studied theoretically after ionization of D$_{2}$ . As a completion of our recent work (J. Chem. Phys. 143, 014305 (2015)) a comprehensive dependence on the pulse duration and delay time is presented in terms of total dissociation probabilities. Our numerical analysis carried out in the recently introduced light-induced conical intersection (LICI) framework clearly shows the effects of the changing position of the LICI which is induced by the frequency modulation of the chirped laser pulses. This impact is presented for positively, negatively and zero chirped short pulses.

preprint2015arXiv

Many-body tunneling dynamics of Bose-Einstein condensates and vortex states in two spatial dimensions

In this work, we study the out-of-equilibrium many-body tunneling dynamics of a Bose-Einstein condensate in a two-dimensional radial double well. We investigate the impact of interparticle repulsion and compare the influence of angular momentum on the many-body tunneling dynamics. Accurate many-body dynamics are obtained by solving the full many-body Schrödinger equation. We demonstrate that macroscopic vortex states of definite total angular momentum indeed tunnel and that, even in the regime of weak repulsions, a many-body treatment is necessary to capture the correct tunneling dynamics. As a general rule, many-body effects set in at weaker interactions when the tunneling system carries angular momentum.

preprint2015arXiv

The role of metal ions in X-ray induced photochemistry

Metal ions play numerous important roles in biological systems being central to the function of biomolecules. In this letter we show that the absorption of X-rays by these ions leads to a complicated chain of ultrafast relaxation steps resulting in the complete degradation of their nearest environment. We conducted high quality ab initio studies on microsolvated Mg clusters demonstrating that ionisation of an 1s-electron of Mg leads to a complicated electronic cascade comprising both intra- and intermolecular steps and lasting only a few hundreds femtoseconds. The metal cation reverts to its original charge state at the end of the cascade, while the nearest solvation shell becomes multiply ionised and large concentrations of radical and slow electron species build up in the metal vicinity. We conclude that such cascades involving metal ions are essential for understanding the radiation chemistry of solutions and radiation damage to metal containing biomolecules.

preprint2014arXiv

Breaking the resilience of a two-dimensional Bose-Einstein condensate to fragmentation

A two-dimensional Bose-Einstein condensate (BEC) split by a radial potential barrier is investigated. We determine on an accurate many-body level the system's ground-state phase diagram as well as a time-dependent phase diagram of the splitting process. Whereas the ground state is condensed for a wide range of parameters, the time-dependent splitting process leads to substantial fragmentation. We demonstrate for the first time the dynamical fragmentation of a BEC despite its ground state being condensed. The results are analyzed by a mean-field model and suggest that a large manifold of low-lying fragmented excited states can significantly impact the dynamics of trapped two-dimensional BECs.

preprint2014arXiv

Direct signature of light-induced conical intersections in diatomics

Nonadiabatic effects are ubiquitous in physics, chemistry and biology. They are strongly amplified by conical intersections (CIs) which are degeneracies between electronic states of triatomic or larger molecules. A few years ago it has been revealed that CIs in molecular systems can be formed by laser light even in diatomics. Due to the prevailing strong nonadiabatic couplings, the existence of such laser-induced conical intersections (LICIs) may considerably change the dynamical behavior of molecular systems. By analyzing the photodissociation process of the D2+ molecule carefully, we found a robust effect in the angular distribution of the photofragments which serves as a direct signature of the LICI providing undoubted evidence for its existence.

preprint2014arXiv

Universality of Fragmentation in the Schrödinger Dynamics of Bosonic Josephson Junctions

The many-body Schrödinger dynamics of a one-dimensional bosonic Josephson junction is investigated for up to ten thousand bosons and long times. The initial states are fully condensed and the interaction strength is weak. We report on a universal fragmentation dynamics on the many-body level: systems consisting of different numbers of particles fragment to the same value at constant mean-field interaction strength. The phenomenon manifests itself in observables such as the correlation functions of the system. We explain this universal fragmentation dynamics analytically based on the Bose-Hubbard model. We thereby show that the extent to which many-body effects become important at later times depends crucially on the initial state. Even for arbitrarily large particle numbers and arbitrarily weak interaction strength the dynamics is many-body in nature and the fragmentation universal. There is no weakly interacting limit where the Gross-Piatevskii mean-field is valid for long times.

preprint2013arXiv

A new intermolecular mechanism to selectively drive photoinduced damages

Low-energy electrons (LEEs) are known to be effective in causing strand breaks in DNA. Recent experiments show that an important direct source of LEEs is the intermolecular Coulombic decay (ICD) process. Here we propose a new cascade mechanism initiated by core excitation and terminated by ICD and demonstrate its properties. Explicit calculations show that the energies of the emitted ICD-electrons can be controlled by selecting the initial atomic excitation. The properties of the cascade may have interesting applications in the fields of electron spectroscopy and radiation damage. Initiating such a cascade by resonant X-ray absorption from a high-Z element embedded in a cancerous cell nucleus, ICD will deliver genotoxic particles \textit{locally} at the absorption site, increasing in that way the controllability of the induced damage.

preprint2013arXiv

Controlled energy-selected electron capture and release in double quantum dots

Highly accurate quantum electron dynamics calculations demonstrate that energy can be efficiently transferred between quantum dots. Specifically, in a double quantum dot an incoming electron is captured by one dot and the excess energy is transferred to the neighboring dot and used to remove an electron from this dot. This process is due to long-range electron correlation and shown to be operative at rather large distances between the dots. The efficiency of the process is greatly enhanced by preparing the double quantum dot such that the incoming electron is initially captured by a two-electron resonance state of the system. In contrast to atoms and molecules in nature, double quantum dots can be manipulated to achieve this enhancement. This mechanism leads to a surprisingly narrow distribution of the energy of the electron removed in the process which is explained by resonance theory. We argue that the process could be exploited in practice.

preprint2013arXiv

Controlling the Velocities and Number of Emitted Particles in the Tunneling to Open Space Dynamics

A scheme to control the many-boson tunneling process to open space is derived and demonstrated. The number of ejected particles and their velocities can be controlled by two parameters, the threshold of the potential and the interparticle interaction. Since these parameters are fully under experimental control, this is also the case for the number of ejected particles and their emission spectrum. The process of tunneling to open space can hence be used, for example, for the quantum simulation of complicated tunneling ionization processes and atom lasers. To understand the many-body tunneling process, a generalization of the model introduced in [Proc. Natl. Acad. Sci. USA, 109, 13521 (2012)] for tunneling in the absence of a threshold is put forward and proven to apply for systems with a non-zero threshold value. It is demonstrated that the model is applicable for general interparticle interaction strengths, particle numbers and threshold values. The model constructs the many-body process from single-particle emission processes. The rates and emission momenta of the single-particle processes are determined by the chemical potentials and energy differences to the threshold value of the potential for systems with different particle numbers. The chemical potentials and these energy differences depend on the interparticle interaction. Both the number of confined particles and their rate of emission thus allow for a control by the manipulation of the interparticle interaction and the threshold. Numerically exact results for two, three and one hundred bosons are shown and discussed. The devised control scheme for the many-body tunneling process performs very well for the dynamics of the momentum density, the correlations, the coherence and of the final state, i.e., the number of particles that remain confined in the potential.

preprint2013arXiv

Elastic scattering of a Bose-Einstein condensate at a potential landscape

We investigate the elastic scattering of Bose-Einstein condensates at shallow periodic and disorder potentials. We show that the collective scattering of the macroscopic quantum object couples to internal degrees of freedom of the Bose-Einstein condensate such that the Bose-Einstein condensate gets depleted. As a precursor for the excitation of the Bose-Einstein condensate we observe wave chaos within a mean-field theory.

preprint2013arXiv

Excitation spectra of many-body systems by linear response: General theory and applications to trapped condensates

We derive a general linear-response many-body theory capable of computing excitation spectra of trapped interacting bosonic systems, e.g., depleted and fragmented Bose-Einstein condensates (BECs). To obtain the linear-response equations we linearize the multiconfigurational time-dependent Hartree for bosons (MCTDHB) method, which provides a self-consistent description of many-boson systems in terms of orbitals and a state vector (configurations), and is in principle numerically-exact. The derived linear-response many-body theory, which we term LR-MCTDHB, is applicable to systems with interaction potentials of general form. From the numerical implementation of the LR-MCTDHB equations and solution of the underlying eigenvalue problem, we obtain excitations beyond available theories of excitation spectra, such as the Bogoliubov-de Gennes (BdG) equations. The derived theory is first applied to study BECs in a one-dimensional harmonic potential. The LR-MCTDHB method contains the BdG excitations and, also, predicts a plethora of additional many-body excitations which are out of the realm of standard linear response. In particular, our theory describes the exact energy of the higher harmonic of the first (dipole) excitation not contained in the BdG theory. We next study a BEC in a very shallow one-dimensional double-well potential. We find with LR-MCTDHB low-lying excitations which are not accounted for by BdG, even though the BEC has only little fragmentation and, hence, the BdG theory is expected to be valid. The convergence of the LR-MCTDHB theory is assessed by systematically comparing the excitation spectra computed at several different levels of theory.

preprint2013arXiv

Generic regimes of quantum many-body dynamics of trapped bosonic systems with strong repulsive interactions

Two generically different but universal dynamical quantum many-body behaviors are discovered by probing the stability of trapped fragmented bosonic systems with strong repulsive finite/long range inter-particle interactions. We use different time-dependent processes to destabilize the systems -- a sudden displacement of the trap is accompanied by a sudden quench of the strength of the inter-particle repulsion. A rather moderate non-violent evolution of the density in the first "topology-preserved" scenario is contrasted with a highly-non-equilibrium dynamics characterizing an explosive changes of the density profiles in the second scenario. The many-body physics behind is identified and interpreted in terms of self-induced time-dependent barriers governing the respective under- and over-a-barrier dynamical evolutions. The universality of the discovered scenarios is explicitly confirmed in 1D, 2D and 3D many-body computations in (a)symmetric traps and repulsive finite/long range inter-particle interaction potentials of different shapes. Implications are briefly discussed.

preprint2013arXiv

Nuclear wave packet quantum interference in the intense laser dissociation of the $\mathrm{D}_{2}^{+}$ molecule

Recently it has been recognized that electronic conical intersections in molecular systems can be induced by laser light even in diatomics. As is known a direct consequence of these accidental degeneracies is the appearence of nonadiabatic effects which has a strong impact on the nuclear quantum dynamics. Studying the photodissociation process of the $\mathrm{D}_{2}^{+}$ molecule, we report here some novel and observable quantum interference phenomena that arise from the topological singularity induced by a strong laser field.

preprint2013arXiv

Photoionization of hydrogen atom by coherent intense high-frequency short laser pulses: Direct propagation of electron wave packets on large spatial grids

The time-dependent Schrödinger equation for the hydrogen atom and its interaction with coherent intense high-frequency short laser pulses is solved numerically exactly by propagating the single-electron wave packets. Thereby, the wavefunction is followed in space and time for times longer than the pulse duration. Results are explicitly shown for 3 and 10 fs pulses. Particular attention is paid to identifying the effect of dynamic interference of photoelectrons emitted with the same kinetic energy at different times during the rising and falling sides of the pulse predicted in [\emph{Ph.V. Demekhin and L.S. Cederbaum}, Phys. Rev. Lett. \textbf{108}, 253001 (2012)]. In order to be able to see the dynamic interference pattern in the computed electron spectra, the photoelectron wave packet has to be propagated over long distances. Clearly, complex absorption potentials often employed to compute spectra of emitted particles cannot be used to detect dynamic interference. For the considered high-frequency pulses of 3 and 10 fs durations, this requires enormously large spatial grids. The presently computed photoionization and above-threshold ionization spectra are found to exhibit pronounced dynamic interference patterns. Where available, the patterns are in very good agreement with previously published results on the photoionization spectra which have been computed using a completely different method, thus supporting the previously made assumption that the above-threshold ionization processes are very weak for the considered pulse intensities and high carrier frequency. The quiver motion in space and time of a free electron in strong laser pulses is also investigated numerically. Finally, a discussion is presented of how fast the atom is ionized by an intense pulse.

preprint2013arXiv

The Exact Wavefunction Factorization of a Vibronic Coupling System

We investigate the exact wavefunction as a single product of electronic and nuclear wavefunction for a model conical intersection system. Exact factorized spiky potentials and nodeless nuclear wavefunctions are found. The exact factorized potential preserves the symmetry breaking effect when the coupling mode is present. Additionally the nodeless wavefunctions are found to be closely related to the adiabatic nuclear eigenfunctions. This phenomenon holds even for the regime where the non-adiabatic coupling is relevant, and sheds light on the relation between the exact wavefunction factorization and the adiabatic approximation.

preprint2013arXiv

Unified view on linear response of interacting identical and distinguishable particles from multiconfigurational time-dependent Hartree methods

A unified view on linear response of interacting systems utilizing multicongurational time-dependent Hartree methods is presented. The cases of one-particle and two-particle response operators for identical particles and up to all-system response operators for distinguishable degrees-of-freedom are considered. The working equations for systems of identical bosons (LR-MCTDHB) and identical fermions (LR-MCTDHF), as well for systems of distinguishable particles (LR-MCTDH) are explicitly derived. These linear-response theories provide numerically-exact excitation energies and system's properties, when numerical convergence is achieved in the calculations.

preprint2012arXiv

Coherent intense resonant laser pulses lead to interference in the time domain observable in the spectrum of the emitted particles

The dynamics of atomic levels resonantly coupled by a coherent and intense short high-frequency laser pulse is discussed and it is advocated that this dynamics is sensitively probed by measuring the spectra of the particles emitted. It is demonstrated that the time-envelope of this laser pulse gives rise to two waves emitted with a time delay with respect to each other at the rising and falling sides of the pulse, which interfere in the time domain. By computing numerically and analyzing explicitly analytically a show-case example of sequential two-photon ionization of an atom by resonant laser pulses, we argue that this dynamic interference should be a general phenomenon in the spectroscopy of strong laser fields. The emitted particles do not have to be photoelectrons. Our results allow also to interpret the already studied resonant Auger effect of an atom by intense free electron laser pulses, and also to envisage experiments in which photons are emitted.

preprint2012arXiv

Dynamics and symmetries of a repulsively bound atom pair in an infinite optical lattice

We investigate the dynamics of two bosons trapped in an infinite one-dimensional optical lattice potential within the framework of the Bose-Hubbard model and derive an exact expression for the wavefunction at finite time. As initial condition we chose localized atoms that are separated by a distance of $d$ lattice sites and carry a center of mass quasi-momentum. An initially localized pair ($d=0$) is found to be more stable as quantified by the pair probability (probability to find two atoms at the same lattice site) when the interaction and/or the center of mass quasi-momentum is increased. For initially separated atoms ($d \neq 0$) there exists an optimal interaction strength for pair formation. Simple expressions for the wavefunction, the pair probability and the optimal interaction strength for pair formation are computed in the limit of infinite time. Whereas the time-dependent wavefunction differs for values of the interaction strength that differ only by the sign, important observables like the density and the pair probability do not. With a symmetry analysis this behavior is shown to extend to the $N$-particle level and to fermionic systems. Our results provide a complementary understanding of the recently observed [Winkler \textit{et al.}, Nature (London) \textbf{441}, 853 (2006)] dynamical stability of atom pairs in a repulsively interacting lattice gas.

preprint2012arXiv

Exact quantum dynamics of bosons with finite-range time-dependent interactions of harmonic type

The exactly solvable quantum many-particle model with harmonic one- and two-particle interaction terms is extended to include time-dependency. We show that when the external trap potential and finite-range interparticle interaction have a time-dependency the exact solutions of the corresponding time-dependent many-boson Schrödinger equation are still available. We use these exact solutions to benchmark the recently developed multiconfigurational time-dependent Hartree method for bosons (MCTDHB) [Phys. Rev. Lett. {\bf 99}, 030402 (2007), Phys. Rev. A {\bf 77}, 033613 (2008)]. In particular, we benchmark the MCTDHB method for: (i) the ground state; (ii) the breathing many-body dynamics activated by a quench scenario where the interparticle interaction strength is suddenly turned on to a finite value; (iii) the non-equilibrium dynamic for driven scenarios where both the trap- and interparticle-interaction potentials are {\it time-dependent}. Excellent convergence of the ground state and dynamics is demonstrated. The great relevance of the self-consistency and time-adaptivity, which are the intrinsic features of the MCTDHB method, is demonstrated by contrasting the MCTDHB predictions and those obtained within the standard full configuration interaction method spanning the Fock space of the same size, but utilizing as one-particle basis set the fixed-shape eigenstates of the one-particle potential. Connections of the model's results to ultra-cold Bose-Einstein condensed systems are addressed.

preprint2012arXiv

Excitation spectra of fragmented condensates by linear response: General theory and application to a condensate in a double-well potential

Linear response of simple (i.e., condensed) Bose-Einstein condensates is known to lead to the Bogoliubov- de Gennes equations. Here, we derive linear response for fragmented Bose-Einstein condensates, i.e., for the case where the many-body wave function is not a product of one, but of several single-particle states (orbitals). Our approach is based on the number-conserving variational time-dependent mean field theory [O. E. Alon, A. I. Streltsov, and L. S. Cederbaum, Phys. Lett. A 362, 453 (2007)], which describes the time evolution of best-mean field states. Correspondingly, we call our linear response theory for fragmented states LR-BMF. In the derivation it follows naturally that excitations are orthogonal to the ground-state orbitals. As applications excitation spectra of Bose-Einstein condensates in double-well potentials are calculated. Both symmetric and asymmetric double-wells are studied for several interaction strengths and barrier heights. The cases of condensed and two-fold fragmented ground states are compared. Interestingly, even in such situations where the response frequencies of the two cases are computed to be close to each other, which is the situation for the excitations well below the barrier, striking differences in the density response in momentum space are found. For excitations with an energy of the order of the barrier height, both the energies and the density response of condensed and fragmented systems are very different. In fragmented systems there is a class of "swapped" excitations where an atom is transfered to the neighboring well. The mechanism of its origin is discussed. In asymmetric wells, the response of a fragmented system is purely local (i.e., finite in either one or the other well) with different frequencies for the left and right fragments.

preprint2012arXiv

Kinetic Energy Release in Fragmentation Processes following Electron Emission: A time dependent approach

A time-dependent approach for the kinetic energy release (KER) spectrum is developed for a fragmentation of a diatomic molecule after an electronic decay process, e.g. Auger process. It allows one to simulate the time-resolved spectra and provides more insight into the molecular dynamics than the time-independent approach. Detailed analysis of the time-resolved emitted electron and KER spectra sheds light on the interrelation between wave packet dynamics and spectra.

preprint2012arXiv

Two trapped particles interacting by a finite-ranged two-body potential in two spatial dimensions

We examine the problem of two particles confined in an isotropic harmonic trap, which interact via a finite-ranged Gaussian-shaped potential in two spatial dimensions. We derive an approximative transcendental equation for the energy and study the resulting spectrum as a function of the interparticle interaction strength. Both the attractive and repulsive systems are analyzed. We study the impact of the potential's range on the ground-state energy. Complementary, we also explicitly verify by a variational treatment that in the zero-range limit the positive delta potential in two dimensions only reproduces the non-interacting results, if the Hilbert space in not truncated. Finally, we establish and discuss the connection between our finite-range treatment and regularized zero-range results from the literature.

preprint2012arXiv

Wave chaos as signature for depletion of a Bose-Einstein condensate

We study the expansion of repulsively interacting Bose-Einstein condensates (BECs) in shallow one-dimensional potentials. We show for these systems that the onset of wave chaos in the Gross-Pitaevskii equation (GPE), i.e. the onset of exponential separation in Hilbert space of two nearby condensate wave functions, can be used as indication for the onset of depletion of the BEC and the occupation of excited modes within a many-body description. Comparison between the multiconfigurational time-dependent Hartree for bosons (MCTDHB) method and the GPE reveals a close correspondence between the many-body effect of depletion and the mean-field effect of wave chaos for a wide range of single-particle external potentials. In the regime of wave chaos the GPE fails to account for the fine-scale quantum fluctuations because many-body effects beyond the validity of the GPE are non-negligible. Surprisingly, despite the failure of the GPE to account for the depletion, coarse grained expectation values of the single-particle density such as the overall width of the atomic cloud agree very well with the many-body simulations. The time dependent depletion of the condensate could be investigated experimentally, e.g., via decay of coherence of the expanding atom cloud.

preprint2011arXiv

Exploring Interatomic Coulombic Decay by Free Electron Lasers

To exploit the high intensity of laser radiation, we propose to select frequencies at which single-photon absorption is of too low energy and two or more photons are needed to produce states of an atom that can undergo interatomic Coulombic decay (ICD) with its neighbors. For Ne dimer it is explicitly demonstrated that the proposed scheme to investigate interatomic processes by multiphoton absorption is much more efficient than with single-photon absorption of sufficiently large frequency as used until now. Extensive calculations on Ne dimer including all the involved nuclear dynamics and the losses by ionization of the participating states show how the low-energy ICD electrons and Ne$^+$ pairs are produced for different laser intensities and pulse durations. At higher intensities the production of Ne$^+$ pairs by successive ionization of the two atoms becomes competitive and the respective emitted electrons interfere coherently with the ICD electrons. It is also demonstrated that a measurement after a time delay can be used to determine the contribution of ICD even at high laser intensity.

preprint2011arXiv

Interrelation between the distributions of kinetic energy release and emitted electron energy following a decay of electronic states

In an electronic decay process followed by fragmentation the kinetic energy release and electron spectra can be measured. Classically they are the mirror image of each other, a fact which is often used in practice. Quantum expressions are derived for both spectra and analyzed. It is demonstrated that these spectra carry complementary quantum information and are related to the nuclear dynamics in different participating electronic states. Illustrative examples show that the classical picture of mirror image can break down and shed light on the underlying physics.

preprint2011arXiv

Number fluctuations of cold spatially split bosonic objects

We investigate the number fluctuations of spatially split many-boson systems employing a theorem about the maximally and minimally attainable variances of an observable. The number fluctuations of many-boson systems are given for different numbers of lattice sites and both mean-field and many-body wave functions. It is shown which states maximize the particle number fluctuations, both in lattices and double-wells. The fragmentation of the states is discussed, and it is shown that the number fluctuations of some fragmented states are identical to those of fully condensed states.

preprint2011arXiv

Optimal time-dependent lattice models for nonequilibrium dynamics

Lattice models are abundant in theoretical and condensed-matter physics. Generally, lattice models contain time-independent hopping and interaction parameters that are derived from the Wannier functions of the noninteracting problem. Here, we present a new concept based on time-dependent Wannier functions and the variational principle that leads to optimal time-dependent lattice models. As an application, we use the Bose-Hubbard model with time-dependent Wannier functions to study a quench scenario involving higher bands. We find a separation of times scales in the dynamics and show that under some circumstances the multi-band nonequilibrium dynamics of a quantum system can be obtained essentially at the cost of a single-band model.

preprint2011arXiv

Recursive formulation of the multiconfigurational time-dependent Hartree method for fermions, bosons and mixtures thereof in terms of one-body density operators

The multiconfigurational time-dependent Hartree method (MCTDH) [Chem. Phys. Lett. {\bf 165}, 73 (1990); J. Chem. Phys. {\bf 97}, 3199 (1992)] is celebrating nowadays entering its third decade of tackling numerically-exactly a broad range of correlated multi-dimensional non-equilibrium quantum dynamical systems. Taking in recent years particles' statistics explicitly into account, within the MCTDH for fermions (MCTDHF) and for bosons (MCTDHB), has opened up further opportunities to treat larger systems of interacting identical particles, primarily in laser-atom and cold-atom physics. With the increase of experimental capabilities to simultaneously trap mixtures of two, three, and possibly even multiple kinds of interacting composite identical particles together, we set up the stage in the present work and specify the MCTDH method for such cases. Explicitly, the MCTDH method for systems with three kinds of identical particles interacting via all combinations of two- and three-body forces is presented, and the resulting equations-of-motion are briefly discussed. All four possible mixtures of fermions and bosons are presented in a unified manner. Particular attention is paid to represent the coefficients' part of the equations-of-motion in a compact recursive form in terms of one-body density operators only. The recursion utilizes the recently proposed Combinadic-based mapping for fermionic and bosonic operators in Fock space [Phys. Rev. A {\bf 81}, 022124 (2010)] and successfully applied and implemented within MCTDHB. Our work sheds new light on the representation of the coefficients' part in MCTDHF and MCTDHB without resorting to the matrix elements of the many-body Hamiltonian with respect to the time-dependent configurations. It suggests a recipe for efficient implementation of the schemes derived here for mixtures which is suitable for parallelization.

preprint2011arXiv

Resonant Auger decay of the core-excited C$^\ast$O molecule in intense X-ray laser fields

The dynamics of the resonant Auger (RA) process of the core-excited C$^\ast$O(1s$^{-1}π^\ast,v_r=0$) molecule in an intense X-ray laser field is studied theoretically. The theoretical approach includes the analogue of the conical intersections of the complex potential energy surfaces of the ground and `dressed' resonant states due to intense X-ray pulses, taking into account the decay of the resonance and the direct photoionization of the ground state, both populating the same final ionic states coherently, as well as the direct photoionization of the resonance state itself. The light-induced non-adiabatic effect of the analogue of the conical intersections of the resulting complex potential energy surfaces gives rise to strong coupling between the electronic, vibrational and rotational degrees of freedom of the diatomic CO molecule. The interplay of the direct photoionization of the ground state and of the decay of the resonance increases dramatically with the field intensity. The coherent population of a final ionic state via both the direct photoionization and the resonant Auger decay channels induces strong interference effects with distinct patterns in the RA electron spectra. The individual impact of these physical processes on the total electron yield and on the CO$^+(A^2Π)$ electron spectrum are demonstrated.

preprint2010arXiv

Death of soliton trains in attractive Bose-Einstein condensates

Experiments on ultra-cold attractive Bose-Einstein Condensates (BECs) have demonstrated that at low dimensions atomic clouds can form localized objects, propagating for long times without significant changes in their shapes and attributed to bright matter-wave solitons, which are coherent objects. We consider the dynamics of bright soliton trains from the perspective of many-boson physics. The fate of matter-wave soliton trains is actually to quickly loose their coherence and become macroscopically fragmented BECs. The death of the coherent matter-wave soliton trains gives birth to fragmented objects, whose quantum properties and experimental signatures differ substantially from what is currently assumed.

preprint2010arXiv

Fragmented Many-Body states of definite angular momentum and stability of attractive 3D Condensates

A three dimensional attractive Bose-Einstein Condensate (BEC) is expected to collapse, when the number of the particles $N$ in the ground state or the interaction strength $λ_0$ exceeds a critical value. We study systems of different particle numbers and interaction strength and find that even if the overall ground state is collapsed there is a plethora of fragmented excited states that are still in the metastable region. Utilizing the \emph{configuration interaction} expansion we determine the spectrum of the ground (`yrast') and excited many-body states with definite total angular momentum quantum numbers $0\leqslant L\leqslant N$ and $-L\leqslant M_L\leqslant L$, and we find and examine states that survive the collapse. This opens up the possibility of realizing a metastable system with overcritical numbers of bosons in a ground state with angular momentum $L\neq0$. The multi-orbital mean-field theory predictions about the existence of fragmented metastable states with overcritical numbers of bosons are verified and elucidated at the many-body level. The descriptions of the total angular momentum within the mean-field and the many-body approaches are compared.

preprint2010arXiv

Mechanism of Tunneling in Interacting Open Ultracold Few-Boson Systems

We investigate the mechanism in the tunneling dynamics of open ultracold few-boson systems by numerically solving the time-dependent few-boson Schrödinger equation exactly. By starting from a weakly repulsive, initially coherent two-boson system we demonstrate that the decay dynamics incorporate fragmentation. The wavefunction of the tunneling state exhibits a pronounced dynamically-stable pattern which we explain by an analytical model. By studying more bosons and stronger interactions we arrive to the conclusion that the decay by tunneling is not a coherent process and exhibits a wealth of phenomena depending on the interaction between the particles.

preprint2010arXiv

Molecular double core-hole electron spectroscopy for chemical analysis

We explore the potential of double core hole electron spectroscopy for chemical analysis in terms of x-ray two-photon photoelectron spectroscopy (XTPPS). The creation of deep single and double core vacancies induces significant reorganization of valence electrons. The corresponding relaxation energies and the interatomic relaxation energies are evaluated by CASSCF calculations. We propose a method how to experimentally extract these quantities by the measurement of single and double core-hole ionization potentials (IPs and DIPs). The influence of the chemical environment on these DIPs is also discussed for states with two holes at the same atomic site and states with two holes at two different atomic sites. Electron density difference between the ground and double core-hole states clearly shows the relaxations accompanying the double core-hole ionization. The effect is also compared with the sensitivity of single core hole ionization potentials (IPs) arising in single core hole electron spectroscopy. We have demonstrated the method for a representative set of small molecules LiF, BeO, BF, CO, N2, C2H2, C2H4, C2H6, CO2 and N2O. The scalar relativistic effect on IPs and on DIPs are briefly addressed.

preprint2010arXiv

On the Cholesky Decomposition for electron propagator methods: General aspects and application on C60

To treat the electronic structure of large molecules by electron propagator methods we developed a parallel computer program called P-RICD$Σ$. The program exploits the sparsity of the two-electron integral matrix by using Cholesky decomposition techniques. The advantage of these techniques is that the error introduced is controlled only by one parameter which can be chosen as small as needed. We verify the tolerance of electron propagator methods to the Cholesky decomposition threshold and demonstrate the power of the P-RICD$Σ$ program for a representative example (C60). All decomposition schemes addressed in the literature are investigated. Even with moderate thresholds the maximal error encountered in the calculated electron affinities and ionization potentials amount to a few meV only, and the error becomes negligible for small thresholds.

preprint2010arXiv

Radiation Generated by Charge Migration Following Ionization

Electronic many-body effects alone can be the driving force for an ultrafast migration of a positive charge created upon ionization of molecular systems. Here we show that this purely electronic phenomenon generates a characteristic IR radiation. The situation when the initial ionic wave packet is produced by a sudden removal of an electron is also studied. It is shown that in this case a much stronger UV emission is generated. This emission appears as an ultrafast response of the remaining electrons to the perturbation caused by the sudden ionization and as such is a universal phenomenon to be expected in every multielectron system.

preprint2010arXiv

Strong impact of light induced conical intersections on the spectrum of diatomic molecules

We show that dressing of diatomic molecules by running laser waves gives rise to conical intersections (CIs). Due to presence of such CIs, the rovibronic molecular motions are strongly coupled. A pronounced impact of the CI on the spectrum of $Na_2$ molecule is demonstrated via numerical calculation for weak and moderate laser intensity, and an experiment is suggested on this basis. The position of the light induced CI and the strength of its non-adiabatic couplings can be chosen by changing the frequency and intensity of the used running laser wave. This offers new possibilities to control the photo-induced rovibronic molecular dynamics.

preprint2010arXiv

Strong interference effects in the resonant Auger decay of atoms induced by intense X-Ray fields

The theory of resonant Auger decay of atoms in a high intensity coherent X-ray pulse is presented. The theory includes the coupling between the ground state and the resonance due to an intense X-ray pulse, taking into account the decay of the resonance and the direct photoionization of the ground state, both populating the final ionic states coherently. The theory also considers the impact of the direct photoionization of the resonance state itself which typically populates highly-excited ionic states. The combined action of the resonant decay and of the direct ionization of the ground state in the field induces a non-hermitian time-dependent coupling between the ground and the 'dressed' resonance stats. The impact of these competing processes on the total electron yield and on the 2s$^2$2p$^{4}(^1\mathrm{D})$3p $^2$P spectator and 2s$^1$2p$^{6}$ $^2$S participator Auger decay spectra of the Ne 1s$\to$3p resonance is investigated. The role of the direct photoionization of the ground state and of the resonance increases dramatically with the field intensity. This results in strong interference effects with distinct patterns in the electron spectra, different for the participator and spectator final states.

preprint2010arXiv

Ultrafast interatomic electronic decay in multiply excited clusters

An ultrafast mechanism belonging to the family of interatomic Coulombic decay (ICD) phenomena is proposed. When two excited species are present, an ultrafast energy transfer can take place bringing one of them to its ground state and ionizing the other one. It is shown that if large homoatomic clusters are exposed to an ultrashort and intense laser pulse whose photon energy is in resonance with an excitation transition of the cluster constituents, the large majority of ions will be produced by this ICD mechanism rather than by two-photon ionization. A related collective-ICD process that is operative in heteroatomic systems is also discussed.

preprint2009arXiv

A novel general mapping for bosonic and fermionic operators in Fock space

In this paper we provide a novel and general way to construct the result of the action of any bosonic or fermionic operator represented in second quantized form on a state vector, without resorting to the matrix representation of operators and even to its elements. The new approach is based on our proposal to compactly enumerate the configurations (i.e., determinants for fermions, permanents for bosons) which are the elements of the state vector. This extremely simplifies the calculation of the action of an operator on a state vector. The computations of statical properties and of the evolution dynamics of a system become much more efficient and applications to systems made of more particles become feasible. Explicit formulations are given for spin-polarized fermionic systems and spinless bosonic systems, as well as to general (two-component) fermionic systems, two-component bosonic systems, and mixtures thereof.

preprint2009arXiv

Accurate multi-boson long-time dynamics in triple-well periodic traps

To solve the many-boson Schrödinger equation we utilize the Multiconfigurational time-dependent Hartree method for bosons (MCTDHB). To be able to attack larger systems and/or to propagate the solution for longer times, we implement a parallel version of the MCTDHB method thereby realizing the recently proposed [Streltsov {\it et al.} arXiv:0910.2577v1] novel idea how to construct efficiently the result of the action of the Hamiltonian on a bosonic state vector. We study the real-space dynamics of repulsive bosonic systems made of N=12, 51 and 3003 bosons in triple-well periodic potentials. The ground state of this system is three-fold fragmented. By suddenly strongly distorting the trap potential, the system performs complex many-body quantum dynamics. At long times it reveals a tendency to an oscillatory behavior around a threefold fragmented state. These oscillations are strongly suppressed and damped by quantum depletions. In spite of the richness of the observed dynamics, the three time-adaptive orbitals of MCTDHB(M=3) are capable to describe the many-boson quantum dynamics of the system for short and intermediate times. For longer times, however, more self-consistent time-adaptive orbitals are needed to correctly describe the non-equilibrium many-body physics. The convergence of the MCTDHB($M$) method with the number $M$ of self-consistent time-dependent orbitals used is demonstrated.

preprint2009arXiv

Environment assisted electron capture

Electron capture by {\it isolated} atoms and ions proceeds by photorecombination. In this process a species captures a free electron by emitting a photon which carries away the excess energy. It is shown here that in the presence of an {\it environment} a competing non-radiative electron capture process can take place due to long range electron correlation. In this interatomic (intermolecular) process the excess energy is transferred to neighboring species. The asymptotic expression for the cross section of this process is derived. We demonstrate by explicit examples that under realizable conditions the cross section of this interatomic process can clearly dominate that of photorecombination.

preprint2009arXiv

Quantum dynamics of attractive versus repulsive bosonic Josephson junctions: Bose-Hubbard and full-Hamiltonian results

The quantum dynamics of one-dimensional bosonic Josephson junctions with attractive and repulsive interparticle interactions is studied using the Bose-Hubbard model and by numerically-exact computations of the full many-body Hamiltonian. A symmetry present in the Bose-Hubbard Hamiltonian dictates an equivalence between the evolution in time of attractive and repulsive Josephson junctions with attractive and repulsive interactions of equal magnitude. The full many-body Hamiltonian does not possess this symmetry and consequently the dynamics of the attractive and repulsive junctions are different.

preprint2008arXiv

Ultrafast electron dynamics following outer-valence ionization: The impact of low-lying relaxation satellite states

Low-lying relaxation satellites give rise to ultrafast electron dynamics following outer-valence ionization of a molecular system is studied. To demonstrate the impact of such satellites, the evolution of the electronic cloud after sudden removal of an electron from the highest occupied molecular orbital (HOMO) of the organic unsaturated nitroso compound 2-Nitroso[1,3]oxazolo[5,4-d][1,3]oxazole is traced in real time and space using ab initio methods only. Our results show that the initially created hole charge remains stationary but on top of it the system reacts by an ultrafast pi-pi* excitation followed by a cyclic excitation-deexcitation process which leads to a redistribution of the charge. The pi-pi* excitation following the removal of the HOMO electron takes place on a sub-femtosecond time scale and the period of the excitation-deexcitation alternations is about 1.4 fs. In real space the processes of excitation and de-excitation represent ultrafast delocalization and localization of the charge. The results are analyzed by a simple two- and three-state model.

preprint2005arXiv

Pathway from condensation via fragmentation to fermionization of cold bosonic systems

For small scattering lengths, cold bosonic atoms form a condensate the density profile of which is smooth. With increasing scattering length, the density {\it gradually} acquires more and more oscillations. Finally, the number of oscillations equals the number of bosons and the system becomes {\it fermionized}. On this pathway from condensation to fermionization intriguing phenomena occur, depending on the shape of the trap. These include macroscopic fragmentation and {\it coexistence} of condensed and fermionized parts that are separated in space.