The RMT method for many-electron atomic systems in intense short-pulse laser light

We describe a new ab initio method for solving the time-dependent Schrödinger equation for multi-electron atomic systems exposed to intense short-pulse laser light. We call the method the R-matrix with time-dependence (RMT) method. Our starting point is a finite-difference numerical integrator (HELIUM), which has proved successful at describing few-electron atoms and atomic ions in strong laser fields with high accuracy. By exploiting the R-matrix division-of-space concept, we bring together a numerical method most appropriate to the multi-electron finite inner region (R-matrix basis set) and a different numerical method most appropriate to the one-electron outer region (finite difference). In order to exploit massively parallel supercomputers efficiently, we time-propagate the wavefunction in both regions by employing Arnoldi methods, originally developed for HELIUM.

Liouville-space R-matrix-Floquet description of atomic radiative processes involving autoionizing states in the presence of intense electromagnetic fields

A reduced-density-operator description is developed for coherent optical phenomena in many-electron atomic systems, utilizing a Liouville-space, multiple-mode Floquet–Fourier representation. The Liouville-space formulation provides a natural generalization of the ordinary Hilbert-space (Hamiltonian) R-matrix-Floquet method, which has been developed for multi-photon transitions and laser-assisted electron–atom collision processes. In these applications, the R-matrix-Floquet method has been demonstrated to be capable of providing an accurate representation of the complex, multi-level structure of many-electron atomic systems in bound, continuum, and autoionizing states. The ordinary Hilbert-space (Hamiltonian) formulation of the R-matrix-Floquet method has been implemented in highly developed computer programs, which can provide a non-perturbative treatment of the interaction of a classical, multiple-mode electromagnetic field with a quantum system. This quantum system may correspond to a many-electron, bound atomic system and a single continuum electron. However, including pseudo-states in the expansion of the many-electron atomic wave function can provide a representation of multiple continuum electrons. The 'dressed' many-electron atomic states thereby obtained can be used in a realistic non-perturbative evaluation of the transition probabilities for an extensive class of atomic collision and radiation processes in the presence of intense electromagnetic fields. In order to incorporate environmental relaxation and decoherence phenomena, we propose to utilize the ordinary Hilbert-space (Hamiltonian) R-matrix-Floquet method as a starting-point for a Liouville-space (reduced-density-operator) formulation. To illustrate how the Liouville-space R-matrix-Floquet formulation can be implemented for coherent atomic radiative processes, we discuss applications to electromagnetically induced transparency, as well as to related pump–probe optical phenomena, and also to the unified description of radiative and dielectronic recombination in electron–ion beam interactions and high-temperature plasmas.

Many-body theory calculations of positron binding to negative ions

A many-body theory approach developed by the authors [Phys. Rev. A 70, 032720 (2004)] is applied to positron bound states and annihilation rates in atomic systems. Within the formalism, full account of virtual positronium (Ps) formation is made by summing the electron-positron ladder diagram series, thus enabling the theory to include all important many-body correlation effects in the positron problem. Numerical calculations have been performed for positron bound states with the hydrogen and halogen negative ions, also known as Ps hydride and Ps halides. The Ps binding energies of 1.118, 2.718, 2.245, 1.873 and 1.393 eV and annihilation rates of 2.544, 2.482, 1.984, 1.913 and 1.809 ns^{-1}, have been obtained for PsH, PsF, PsCl, PsBr and PsI, respectively.

Time delay between photoemission from the 2p and 2s subshells of neon

The R-matrix incorporating time (RMT) method is a method developed recently for solving the time-dependent Schrödinger equation for multielectron atomic systems exposed to intense short-pulse laser light. We have employed the RMT method to investigate the time delay in the photoemission of an electron liberated from a 2p orbital in a neon atom with respect to one released from a 2s orbital following absorption of an attosecond xuv pulse. Time delays due to xuv pulses in the range 76-105 eV are presented. For an xuv pulse at the experimentally relevant energy of 105.2 eV, we calculate the time delay to be 10.2±1.3 attoseconds (as), somewhat larger than estimated by other theoretical calculations, but still a factor of 2 smaller than experiment. We repeated the calculation for a photon energy of 89.8 eV with a larger basis set capable of modeling correlated-electron dynamics within the neon atom and the residual Ne ion. A time delay of 14.5±1.5 as was observed, compared to a 16.7±1.5 as result using a single-configuration representation of the residual Ne⁺ ion. 

Transferring entanglement to the steady state of flying qubits

The transfer of entanglement from optical fields to qubits provides a viable approach to entangling remote qubits in a quantum network. In cavity quantum electrodynamics, the scheme relies on the interaction between a photonic resource and two stationary intracavity atomic qubits. However, it might be hard in practice to trap two atoms simultaneously and synchronize their coupling to the cavities. To address this point, we propose and study entanglement transfer from cavities driven by an entangled external field to controlled flying qubits. We consider two exemplary non-Gaussian driving fields: NOON and entangled coherent states. We show that in the limit of long coherence time of the cavity fields, when the dynamics is approximately unitary, entanglement is transferred from the driving field to two atomic qubits that cross the cavities. On the other hand, a dissipation-dominated dynamics leads to very weakly quantum-correlated atomic systems, as witnessed by vanishing quantum discord.

Robust nonadiabatic molecular dynamics for metals and insulators

We present a new formulation of the correlated electron-ion dynamics (CEID) scheme, which systematically improves Ehrenfest dynamics by including quantum fluctuations around the mean-field atomic trajectories. We show that the method can simulate models of nonadiabatic electronic transitions and test it against exact integration of the time-dependent Schrodinger equation. Unlike previous formulations of CEID, the accuracy of this scheme depends on a single tunable parameter which sets the level of atomic fluctuations included. The convergence to the exact dynamics by increasing the tunable parameter is demonstrated for a model two level system. This algorithm provides a smooth description of the nonadiabatic electronic transitions which satisfies the kinematic constraints (energy and momentum conservation) and preserves quantum coherence. The applicability of this algorithm to more complex atomic systems is discussed.

Correspondence between continuous-variable and discrete quantum systems of arbitrary dimensions

We establish a mapping between a continuous-variable (CV) quantum system and a discrete quantum system of arbitrary dimension. This opens up the general possibility to perform any quantum information task with a CV system as if it were a discrete system. The Einstein-Podolsky-Rosen state is mapped onto the maximally entangled state in any finite-dimensional Hilbert space and thus can be considered as a universal resource of entanglement. An explicit example of the map and a proposal for its experimental realization are discussed.

Computational challenges in atomic, molecular and optical physics

Six challenges are discussed. These are the laser-driven helium atom; the laser-driven hydrogen molecule and hydrogen molecular ion: electron scattering (with ionization) from one-electron atoms; the vibrational and rotational structure of molecules such as H-3(+) and water at their dissociation limits; laser- heated clusters; and quantum degeneracy and Bose-Einstein condensation. The first four concern fundamental few-body systems where use of high-performance computing (HPC) is currently making possible accurate modelling from first principles. This leads to reliable predictions and support for laboratory experiment as well as true understanding of the dynamics. Important aspects of these challenges addressable only via a terascale facility are set out. Such a facility makes the last two challenges in the above list meaningfully accessible for the first time, and the scientific interest together with the prospective role for HPC in these is emphasized.

Multiple proton translocation in biomolecular systems: concerted to stepwise transition in a simple model

In this paper we study a simple model potential energy surface (PES) useful for describing multiple proton translocation mechanisms. The approach presented is relevant to the study of more complex biomolecular systems like enzymes. In this model, at low temperatures, proton tunnelling favours a concerted proton transport mechanism, while at higher temperatures there is a crossover from concerted to stepwise mechanisms; the crossover temperature depends on the energetic features of the PES. We illustrate these ideas by calculating the crossover temperature using energies taken from ab initio calculations on specific systems. Interestingly, typical crossover temperatures lie around room temperature; thus both concerted and stepwise reaction mechanisms should play an important role in biological systems, and one can be easily turned into another by external means such as modifying the temperature or the pH, thus establishing a general mechanism for modulation of the biomolecular function by external effectors.

Concentration and purification of entanglement for qubit systems with ancillary cavity fields

We propose schemes for entanglement concentration and purification for qubit systems encoded in flying atomic pairs. We use cavity-quantum electrodynamics as an illustrative setting within which our proposals can be implemented. Maximally entangled pure states of qubits can be produced as a result of our protocols. In particular, the concentration protocol yields Bell states with the largest achievable theoretical probability while the purification scheme produces arbitrarily pure Bell states. The requirements for the implementation of these protocols are modest, within the state of the art, and we address all necessary steps in two specific setups based on experimentally mature microwave technology.

Extended electrochemical windows made accessible by room temperature ionic liquid/organic solvent electrolyte systems

The electrochemical windows of acetonitrile solutions doped with 0.1 m concentrations of several ionic liquids were examined by cyclic voltammetry at gold and platinum microelectrodes. These results were compared with those observed in the commonly used 0.1 m tetrabutylammonium perchlorate/acetonitrile system as well as with neat ionic liquids. The use of a trifluorotris(pentofluoroethyl)phosphate-based ionic liquid, specifically, as supporting electrolyte in acetonitrile solutions affords a wider anodic window, which is attributed to the high stability of the anionic component of these intrinsically conductive and thermally robust compounds.

Newtonian origin of the spin motive force in ferromagnetic atomic wires

We demonstrate numerically the existence of a spin-motive force acting on spin carriers when moving in a time and space dependent internal ?eld. This is the case for electrons in a one-dimensional wire with a precessing domain wall. The effect can be explained solely by adiabatic dynamics and is shown to exist for both classical and quantum systems.

Structural changes in quasi-one-dimensional many-electron systems: From linear to zigzag and beyond

Many-electron systems confined to a quasi-one-dimensional geometry by a cylindrical distribution of positive charge have been investigated by density functional computations in the unrestricted local spin density approximation. Our investigations have been focused on the low-density regime, in which electrons are localized. The results reveal a wide variety of different charge and spin configurations, including linear and zig-zag chains, single-and double-strand helices, and twisted chains of dimers. The spin-spin coupling turns from weakly antiferromagnetic at relatively high density, to weakly ferromagnetic at the lowest densities considered in our computations. The stability of linear chains of localized charge has been investigated by analyzing the radial dependence of the self-consistent potential and by computing the dispersion relation of low-energy harmonic excitations.