Excited Eigenstates Research Articles

ON THE SUPERSYMMETRY OF THE DIRAC–KEPLER PROBLEM PLUS A COULOMB-TYPE SCALAR POTENTIAL IN (D+1) DIMENSIONS AND THE GENERALIZED LIPPMANN–JOHNSON OPERATOR

We will study the Dirac–Kepler problem plus a Coulomb-type scalar potential by generalizing the Lippmann–Johnson operator to D spatial dimensions. From this operator, we construct the supersymmetric generators to obtain the energy spectrum for discrete excited eigenstates and the radial spinor for the SUSY ground state.

Communication: Transfer of more than half the population to a selected rovibrational state of H2 by Stark-induced adiabatic Raman passage

By using Stark-induced adiabatic Raman passage (SARP) with partially overlapping nanosecond pump (532 nm) and Stokes (683 nm) laser pulses, 73% ± 6% of the initial ground vibrational state population of H(2) (v = 0, J = 0) is transferred to the single vibrationally excited eigenstate (v = 1, J = 0). In contrast to other Stark chirped Raman adiabatic passage techniques, SARP transfers population from the initial ground state to a vibrationally excited target state of the ground electronic surface without using an intermediate vibronic resonance within an upper electronic state. Parallel linearly polarized, co-propagating pump and Stokes laser pulses of respective durations 6 ns and 4.5 ns, are combined with a relative delay of ~4 ns before orthogonally intersecting the molecular beam of H(2). The pump and Stokes laser pulses have fluences of ~10 J/mm(2) and ~1 J/mm(2), respectively. The intense pump pulse generates the necessary sweeping of the Raman resonance frequency by ac (second-order) Stark shifting the rovibrational levels. As the frequency of the v = 0 → v = 1 Raman transition is swept through resonance in the presence of the strong pump and the weaker delayed Stokes pulses, the population of (v = 0, J = 0) is coherently transferred via an adiabatic passage to (v = 1, J = 0). A quantitative measure of the population transferred to the target state is obtained from the depletion of the ground-state population using 2 + 1 resonance enhanced multiphoton ionization (REMPI) in a time-of-flight mass spectrometer. The depletion is measured by comparing the REMPI signal of (v = 0, J = 0) at Raman resonance with that obtained when the Stokes pulse is detuned from the Stark-shifted Raman resonance. No depletion is observed with either the pump or the Stokes pulses alone, confirming that the measured depletion is indeed caused by the SARP-induced population transfer from the ground to the target state and not by the loss of molecules from photoionization or photodissociation. The two-photon resonant UV pulse used for REMPI detection is delayed by 20 ns with respect to the pump pulse to avoid the ac Stark shift originating from the pump and Stokes laser pulses. This experiment demonstrates the feasibility of preparing a large ensemble of isolated molecules in a preselected single quantum state without requiring an intermediate vibronic resonance.

Adiabatic and nonadiabatic contributions to the energy of a system subject to a time-dependent perturbation: Complete separation and physical interpretation

When a time-dependent perturbation acts on a quantum system that is initially in the nondegenerate ground state ∣0> of an unperturbed Hamiltonian H(0), the wave function acquires excited-state components ∣k> with coefficients c(k)(t) exp(-iE(k)t/ℏ), where E(k) denotes the energy of the unperturbed state ∣k>. It is well known that each coefficient c(k)(t) separates into an adiabatic term a(k)(t) that reflects the adjustment of the ground state to the perturbation--without actual transitions--and a nonadiabatic term b(k)(t) that yields the probability amplitude for a transition to the excited state. In this work, we prove that the energy at any time t also separates completely into adiabatic and nonadiabatic components, after accounting for the secular and normalization terms that appear in the solution of the time-dependent Schrödinger equation via Dirac's method of variation of constants. This result is derived explicitly through third order in the perturbation. We prove that the cross-terms between the adiabatic and nonadiabatic parts of c(k)(t) vanish, when the energy at time t is determined as an expectation value. The adiabatic term in the energy is identical to the total energy obtained from static perturbation theory, for a system exposed to the instantaneous perturbation λH'(t). The nonadiabatic term is a sum over excited states ∣k> of the transition probability multiplied by the transition energy. By evaluating the probabilities of transition to the excited eigenstates ∣k'(t)> of the instantaneous Hamiltonian H(t), we provide a physically transparent explanation of the result for E(t). To lowest order in the perturbation parameter λ, the probability of finding the system in state ∣k'(t)> is given by λ(2) ∣b(k)(t)∣(2). At third order, the transition probability depends on a second-order transition coefficient, derived in this work. We indicate expected differences between the results for transition probabilities obtained from this work and from Fermi's golden rule.

Mathematical modeling of the population dynamics of “dark” states of a three-level system

Mathematical modeling of the population dynamics is performed for states of a three-level system (atom) with a V-type configuration transforming a light pulse. It is assumed that the excited eigenstates of the atom are degenerate and coupled by coherent interaction, one of the states being radiating (radiative), while the other state is nonradiating (“dark”). The population dynamics of atomic states is described on the basis of numerical solutions of equations for the matrix elements of the density operator. The dependence of the efficiency of population of the atomic dark state from the values of the parameters of an irradiation pulse and from the ratio of the period of population oscillations of excited atomic states (caused by their coherent interaction) to the lifetime of the atomic radiating state is determined. Typical examples of the time dependence of the population of states of the atom considered are presented for the cases of irradiation by a short (as compared to the lifetime of the radiating state) sinusoidal light pulse and by a long rectangular light pulse with the resonance carrier frequency.

Conformal blocks in Virasoro and W theories: Duality and the Calogero–Sutherland model

We study the properties of the conformal blocks of the conformal field theories with Virasoro or W-extended symmetry. When the conformal blocks contain only second-order degenerate fields, the conformal blocks obey second order differential equations and they can be interpreted as ground-state wave functions of a trigonometric Calogero–Sutherland Hamiltonian with non-trivial braiding properties. A generalized duality property relates the two types of second order degenerate fields. By studying this duality we found that the excited states of the Calogero–Sutherland Hamiltonian are characterized by two partitions, or in the case of WAk−1 theories by k partitions. By extending the conformal field theories under consideration by a u(1) field, we find that we can put in correspondence the states in the Hilbert state of the extended CFT with the excited non-polynomial eigenstates of the Calogero–Sutherland Hamiltonian. When the action of the Calogero–Sutherland integrals of motion is translated on the Hilbert space, they become identical to the integrals of motion recently discovered by Alba, Fateev, Litvinov and Tarnopolsky in Liouville theory in the context of the AGT conjecture. Upon bosonization, these integrals of motion can be expressed as a sum of two, or in general k, bosonic Calogero–Sutherland Hamiltonian coupled by an interaction term with a triangular structure. For special values of the coupling constant, the conformal blocks can be expressed in terms of Jack polynomials with pairing properties, and they give electron wave functions for special Fractional Quantum Hall states.

Full-dimensional multi configuration time dependent Hartree calculations of the ground and vibrationally excited states of He2,3Br2 clusters

Quantum dynamics calculations are reported for the tetra-, and penta-atomic van der Waals He(N)Br(2) complexes using the multiconfiguration time-dependent Hartree (MCTDH) method. The computations are carried out in satellite coordinates, and the kinetic energy operator in this set of coordinates is given. A scheme for the representation of the potential energy surface based on the sum of the three-body HeBr(2) interactions at CSSD(T) level plus the He-He interaction is employed. The potential surfaces show multiple close lying minima, and a quantum description of such highly floppy multiminima systems is presented. Benchmark, full-dimensional converged results on ground vibrational/zero-point energies are reported and compared with recent experimental data available for all these complexes, as well as with previous variational quantum calculations for the smaller HeBr(2) and He(2)Br(2) complexes on the same surface. Some low-lying vibrationally excited eigenstates are also computed by block improved relaxation calculations. The binding energies and the corresponding vibrationally averaged structures are determined for different conformers of these complexes. Their relative stability is discussed, and contributes to evaluate the importance of the multiple-minima topology of the underlying potential surface.

Positive-parity excited states of the nucleon in quenched lattice QCD

Positive-parity spin-1/2 excitations of the nucleon are explored in lattice QCD. The variational method is used in this investigation and several correlation matrices are employed. As our focus is on the utility and methodology of the variational approach, we work in the quenched approximation to QCD. Various sweeps of Gaussian fermion-field smearing is applied at the source and at the sink of $\chi_{1}\bar\chi_{1}$ and $\chi_{1}\chi_{2}$ correlation functions to obtain a large basis of operators. Using several different approaches for constructing basis interpolators, we demonstrate how improving the basis can split what otherwise might be interpreted as a single state into multiple eigenstates. Consistency of the extracted excited energy states are explored over various dimensions of the correlation matrices. The use of large correlation matrices is emphasized for the reliable extraction of the excited eigenstates of QCD.

ηCondensate of Fermionic Atom Pairs via Adiabatic State Preparation

We discuss how an η condensate, corresponding to an exact excited eigenstate of the Fermi-Hubbard model, can be produced with cold atoms in an optical lattice. Using time-dependent density matrix renormalization group methods, we analyze a state preparation scheme beginning from a band insulator state in an optical superlattice. This state can act as an important test case, both for adiabatic preparation methods and the implementation of the many-body Hamiltonian, and measurements on the final state can be used to help detect associated errors.

Multiple atomic dark solitons in cigar-shaped Bose-Einstein condensates

We consider the stability and dynamics of multiple dark solitons in cigar-shaped Bose-Einstein condensates. Our study is motivated by the fact that multiple matter-wave dark solitons may naturally form in such settings as per our recent work [Phys. Rev. Lett. 101, 130401 (2008)]. First, we study the dark soliton interactions and show that the dynamics of well-separated solitons (i.e., ones that undergo a collision with relatively low velocities) can be analyzed by means of particle-like equations of motion. The latter take into regard the repulsion between solitons (via an effective repulsive potential) and the confinement and dimensionality of the system (via an effective parabolic trap for each soliton). Next, based on the fact that stationary, well-separated dark multisoliton states emerge as a nonlinear continuation of the appropriate excited eigenstates of the quantum harmonic oscillator, we use a Bogoliubov-de Gennes analysis to systematically study the stability of such structures. We find that for a sufficiently large number of atoms, multiple soliton states are dynamically stable, while for a small number of atoms, we predict a dynamical instability emerging from resonance effects between the eigenfrequencies of the soliton modes and the intrinsic excitation frequencies of the condensate. Finally, we present experimental realizations of multisoliton states including a three-soliton state consisting of two solitons oscillating around a stationary one and compare the relevant results to the predictions of the theoretical mean-field model.

Tracing molecular electronic excitation dynamics in real time and space

We present a method for studying the movement of electrons and energy within and between electronically excited molecules. The dynamically changing state is a many-electron wavepacket, for which we numerically integrate the Schrodinger equation using the ADC(2) effective Hamiltonian for the particle-hole propagator. We develop the tools necessary for following the separate motions of the particles and holes. Total particle and hole densities can be used to give an overview of the dynamics, which can be atomically decomposed in a Mulliken fashion, or individual particle and hole states give a more detailed look at the structure of an excitation. We apply our model to a neon chain, as an illustrative example, projecting an excited eigenstate of an isolated atom onto the coupled system as the initial state. In addition to demonstrating our propagation and analysis machinery, the results show a dramatic difference in excitation-energy transfer rates as a consequence of initial polarization. Furthermore, already in a system with three constituents, an important aspect of multiple coupled systems appears, in that one absorbing system essentially shields another, changing the effective sitewise coupling parameters.

Complete transfer of populations from a single state to a preselected superposition of states using piecewise adiabatic passage: Experiment

We demonstrate a method of adiabatic population transfer from a single quantum state into a coherent superposition of states. The transfer is executed with femtosecond pulses, spectrally shaped in a simple and intuitive manner, which does not require iterative feedback-controlled loops. In contrast to nonadiabatic methods of excitation, our approach is not sensitive to the exact value of laser intensity. We show that the population transfer is complete, and analyze the possibility of controlling the relative phases and amplitudes of the excited eigenstates. We discuss the limitations of the proposed control methods due to the dynamic level shifts and suggest ways of reducing their influence.

Full dimensional (15-dimensional) quantum-dynamical simulation of the protonated water-dimer III: Mixed Jacobi-valence parametrization and benchmark results for the zero point energy, vibrationally excited states, and infrared spectrum

Quantum dynamical calculations are reported for the zero point energy, several low-lying vibrational states, and the infrared spectrum of the H(5)O(2)(+) cation. The calculations are performed by the multiconfiguration time-dependent Hartree (MCTDH) method. A new vector parametrization based on a mixed Jacobi-valence description of the system is presented. With this parametrization the potential energy surface coupling is reduced with respect to a full Jacobi description, providing a better convergence of the n-mode representation of the potential. However, new coupling terms appear in the kinetic energy operator. These terms are derived and discussed. A mode-combination scheme based on six combined coordinates is used, and the representation of the 15-dimensional potential in terms of a six-combined mode cluster expansion including up to some 7-dimensional grids is discussed. A statistical analysis of the accuracy of the n-mode representation of the potential at all orders is performed. Benchmark, fully converged results are reported for the zero point energy, which lie within the statistical uncertainty of the reference diffusion Monte Carlo result for this system. Some low-lying vibrationally excited eigenstates are computed by block improved relaxation, illustrating the applicability of the approach to large systems. Benchmark calculations of the linear infrared spectrum are provided, and convergence with increasing size of the time-dependent basis and as a function of the order of the n-mode representation is studied. The calculations presented here make use of recent developments in the parallel version of the MCTDH code, which are briefly discussed. We also show that the infrared spectrum can be computed, to a very good approximation, within D(2d) symmetry, instead of the G(16) symmetry used before, in which the complete rotation of one water molecule with respect to the other is allowed, thus simplifying the dynamical problem.

Decoding the Dynamical Information Embedded in Highly Excited Vibrational Eigenstates: State Space and Phase Space Viewpoints

We study the nature of highly excited eigenstates of strongly coupled multimode systems with three degrees of freedom. Attempts to dynamically assign the eigenstates using classical-quantum correspondence techniques poses a considerable challenge, due to both the number of degrees of freedom and the extensive chaos in the underlying classical phase space. Nevertheless, we show that sequences of localized states interspersed between delocalized states can be identified readily by using the parametric variation technique proposed earlier by us. In addition, we introduce a novel method to lift the quantum eigenstates onto the classical resonance web using a wavelet-based local time-frequency approach. It is shown that the lifting procedure provides clear information on the various dominant nonlinear resonances that influence the eigenstates. Analyzing the spectroscopic Hamiltonians for CDBrClF and CF(3)CHFI as examples, we illustrate our approach and demonstrate the consistency between state space and phase space perspectives of the eigenstates. Two exemplary cases of highly mixed quantum states are discussed to highlight the difficulties associated with their assignment. In particular, we provide arguments to distinguish between the states in terms of their predominantly classical or quantum nature of the mixing.

On the Absence of Excited Eigenstates of Atoms in QED

For the standard model of QED with static nuclei, nonrelativistic electrons and an ultraviolet cutoff, a new simple proof of absence of excited eigenstates with energies above the groundstate energy and below the ionization threshold of an atom is presented. Our proof is based on a multi-scale virial argument and exploits the fact that, in perturbation theory, excited atomic states decay by emission of one or two photons. Our arguments do not require an infrared cutoff (or regularization) and are applicable for all energies above the groundstate energy, except in a small (α-dependent) interval around the ionization threshold.

Theoretical approach for computing magnetic anisotropy in single molecule magnets

We present a theoretical approach to calculate the molecular magnetic anisotropy parameters, $D_{M}$ and $E_M$ for single molecule magnets in any eigenstate of the exchange Hamiltonian, treating the anisotropy Hamiltonian as a perturbation. Neglecting intersite dipolar interactions, we calculate molecular magnetic anisotropy in a given total spin state from the known single-ion anisotropies of the transition metal centers. The method is applied to $Mn_{12}Ac$ and $Fe_8$ in their ground and first few excited eigenstates, as an illustration. We have also studied the effect of orientation of local anisotropies on the molecular anisotropy in various eigenstates of the exchange Hamiltonian. We find that, in case of $Mn_{12}Ac$, the molecular anisotropy depends strongly on the orientation of the local anisotropies and the spin of the state. The $D_M$ value of $Mn_{12}Ac$ is almost independent of the orientation of the local anisotropy of the core Mn(IV) ions. In the case of $Fe_8$, the dependence of molecular anisotropy on the spin of the state in question is weaker. We have also calculated the anisotropy constants for several sets of exchange parameters and found that in $Mn_{12}Ac$ the anisotropy increases with spin excitation gap, while in $Fe_{8}$, the anisotropy is almost independent of the gap.

Towards Toroidal Hydrogen Bonds

Abstract We present hydrogen bonds with toroidal densities of the protons, called toroidal hydrogen bonds, in systems with cylindrical symmetry e.g. triatomic molecules AHB or ions AHB+ or AHB–. These may be prepared by excitation of the degenerate bending vibrations and related pseudorotation such that the hydrogen atom (or the corresponding proton) is no longer located on the symmetry axis between atoms A and B, but – classically speaking – it rotates around that axis. Quantum mechanically, the toroidal hydrogen bond is represented by an excited nuclear eigenstate with nuclear wavefunction and corresponding nuclear density which have toroidal shapes around the central nodal line which conincides with the AB symmetry axis. The properties of these bonds are analyzed, including the pseudorotational angular momentum. Toroidal hydrogen bonds may be excited by means of circularly polarized infrared (IR) laser pulses. The results are demonstrated exemplarily for the oriented model system FHF–, by means of combined quantum chemistry calculations of the potential energy and dipole surfaces (adapted from L. González, G. Pérez-Hernández, J. González-Vázquez, (2008), submitted), calculations of the vibrational and pseudorotational states in the frame of Watson’s isomorphic Hamiltonian for linear molecules (J. K. G. Watson, Mol. Phys. 19 (1970) 465), and quantum dynamics simulations of the laser driven nuclear dynamics, analogous to recent applications to CdH2 (I. Barth, J. Manz, P. Sebald, Chem. Phys. 346 (2008) 89).

Single-particle density matrix and the momentum distribution of dark “solitons” in a Tonks-Girardeau gas

We study the reduced single-particle density matrix (RSPDM), the momentum distribution, and natural orbitals and their occupancies of dark ``soliton'' (DS) states in a Tonks-Girardeau gas. DS states are specially tailored excited many-body eigenstates, which have a dark solitonic notch in their single-particle density. The momentum distribution of DS states has a characteristic shape with two sharp spikes. We find that the two spikes arise due to the high degree of correlation observed within the RSPDM between the mirror points ($x$ and $\ensuremath{-}x$) with respect to the dark notch at $x=0$; the correlations oscillate rather than decay as the points $x$ and $\ensuremath{-}x$ are being separated.

Design of UV laser pulses for the preparation of matrix isolated homonuclear diatomic molecules in selective vibrational superposition states

The preparation of matrix isolated homonuclear diatomic molecules in a vibrational superposition state c0Phie=1,v=0+cjPhie=1,v=j, with large (|c0|2 approximately 1) plus small contributions (|cj|2<<1) of the ground v=0 and specific v=j low excited vibrational eigenstates, respectively, in the electronic ground (e=1) state, and without any net population transfer to electronic excited (e>1) states, is an important challenge; it serves as a prerequisite for coherent spin control. For this purpose, the authors investigate two scenarios of laser pulse control, involving sequential or intrapulse pump- and dump-type transitions via excited vibronic states Phiex,k with a dominant singlet or triplet character. The mechanisms are demonstrated by means of quantum simulations for representative nuclear wave packets on coupled potential energy surfaces, using as an example a one-dimensional model for Cl2 in an Ar matrix. A simple three-state model (including Phi1,0, Phi1,j and Phiex,k) allows illuminating analyses and efficient determinations of the parameters of the laser pulses based on the values of the transition energies and dipole couplings of the transient state which are derived from the absorption spectra.

Calculation of highly excited vibrational states using a Richardson-Leja-Davidson scheme

An efficient computational scheme for calculating highly excited vibrational eigenstates is proposed, combining a Richardson-Leja spectral filter with a novel version of the Davidson method [J. Comput. Phys. 17, 87 (1975)]. Highly excited eigenstates of the Rb2 and H2O molecules are computed to test and verify the method. On the average less than 2.5 outer recursions per eigenstate are needed. For each outer Davidson recursion, less than 20 inner filter recursions per eigenstate are needed on the average.

Theoretical investigation of highly excited vibrational states in DFCO: Calculation of the out-of-plane bending states and simulation of the intramolecular vibrational energy redistribution

A previously developed modified Davidson scheme [C. Iung and F. Ribeiro, J. Chem. Phys. 121, 174105 (2005)] is applied to compute and analyze highly excited (nu2,nu6) eigenstates in DFCO. The present paper is also devoted to the simulations of the intramolecular vibrational energy redistribution (IVR) initiated by an excitation of the out-of-plane bending vibration (nnu6, n=2,4,6, . . . ,18, and 20). The multiconfiguration time-dependent Hartree method is exploited to propagate the corresponding six-dimensional wave packets. A comprehensive comparison with experimental data as well as with previous simulations of IVR in HFCO [G. Pasin et al. J. Chem. Phys. 124, 194304 (2006)] is presented.