An ignition key for atomic-scale engines

A current-carrying resonant nanoscale device, simulated by non-adiabatic molecular dynamics, exhibits sharp activation of non-conservative current-induced forces with bias. The result, above the critical bias, is generalized rotational atomic motion with a large gain in kinetic energy. The activation exploits sharp features in the electronic structure, and constitutes, in effect, an ignition key for atomic-scale motors. A controlling factor for the effect is the non-equilibrium dynamical response matrix for small-amplitude atomic motion under current. This matrix can be found from the steady-state electronic structure by a simpler static calculation, providing a way to detect the likely appearance, or otherwise, of non-conservative dynamics, in advance of real-time modelling.

Atomic harmonic generation in time-dependent R-matrix theory

We have developed the capability to determine accurate harmonic spectra for multielectron atoms within time-dependent R-matrix (TDRM) theory. Harmonic spectra can be calculated using the expectation value of the dipole length, velocity, or acceleration operator. We assess the calculation of the harmonic spectrum from He irradiated by 390-nm laser light with intensities up to 4 x 10(14) W cm(-2) using each form, including the influence of the multielectron basis used in the TDRM code. The spectra are consistent between the different forms, although the dipole acceleration calculation breaks down at lower harmonics. The results obtained from TDRM theory are compared with results from the HELIUM code, finding good quantitative agreement between the methods. We find that bases which include pseudostates give the best comparison with the HELIUM code, but models comprising only physical orbitals also produce accurate results.

Current-induced atomic dynamics, instabilities, and Raman signals: Quasiclassical Langevin equation approach

We derive and employ a semiclassical Langevin equation obtained from path integrals to describe the ionic dynamics of a molecular junction in the presence of electrical current. The electronic environment serves as an effective nonequilibrium bath. The bath results in random forces describing Joule heating, current-induced forces including the nonconservative wind force, dissipative frictional forces, and an effective Lorentz-type force due to the Berry phase of the nonequilibrium electrons. Using a generic two-level molecular model, we highlight the importance of both current-induced forces and Joule heating for the stability of the system. We compare the impact of the different forces, and the wide-band approximation for the electronic structure on our result. We examine the current-induced instabilities (excitation of runaway "waterwheel" modes) and investigate the signature of these in the Raman signals.

Deployment on GPUs of an application in computational atomic physics

This paper describes the deployment on GPUs of PROP, a program of the 2DRMP suite which models electron collisions with H-like atoms and ions. Because performance on GPUs is better in single precision than in double precision, the numerical stability of the PROP program in single precision has been studied. The numerical quality of PROP results computed in single precision and their impact on the next program of the 2DRMP suite has been analyzed. Successive versions of the PROP program on GPUs have been developed in order to improve its performance. Particular attention has been paid to the optimization of data transfers and of linear algebra operations. Performance obtained on several architectures (including NVIDIA Fermi) are presented.

Entanglement control via reservoir engineering in ultracold atomic gases

We study the entanglement of two impurity qubits immersed in a Bose-Einstein condensate (BEC) reservoir. This open quantum system model allows for interpolation between a common dephasing scenario and an independent dephasing scenario by modifying the wavelength of the superlattice superposed to the BEC, and how this influences the dynamical properties of the impurities. We demonstrate the existence of rich dynamics corresponding to different values of reservoir parameters, including phenomena such as entanglement trapping, revivals of entanglement, and entanglement generation. In the spirit of reservoir engineering, we present the optimal BEC parameters for entanglement generation and trapping, showing the key role of the ultracold-gas interactions. Copyright (C) EPLA, 2013

Effect of interparticle interaction in a free-oscillation atomic interferometer

We investigate the dynamics of two interacting bosons repeatedly scattering off a beam-splitter in a free oscillation atom interferometer. Using the interparticle scattering length and the beam-splitter probabilites as our control parameters, we show that even in a simple setup like this a wide range of strongly correlated quantum states can be created. This in particular includes the NOON state, which maximizes the quantum Fisher information and is a foremost state in quantum metrology. DOI: 10.1103/PhysRevA.87.043630

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.

Collective atomic recoil lasing including friction and diffusion effects

We extend the collective atomic recoil lasing (CARL) model including the effects of friction and diffusion forces acting on the atoms due to the presence of optical molasses fields. The results from this model are consistent with those from a recent experiment by Kruse [ Phys. Rev. Lett. 91, 183601 (2003) ]. In particular, we obtain a threshold condition above which collective backscattering occurs. Using a nonlinear analysis we show that the backscattered field and the bunching evolve to a steady state, in contrast to the nonstationary behavior of the standard CARL model. For a proper choice of the parameters, this steady state can be superfluorescent.

Full characterization of the quantum linear-zigzag transition in atomic chains

A string of repulsively interacting particles exhibits a phase transition to a zigzag structure, by reducing the transverse trap potential or the interparticle distance. Based on the emergent symmetry Z2 it has been argued that this instability is a quantum phase transition, which can be mapped to an Ising model in transverse field. An extensive Density Matrix Renormalization Group analysis is performed, resulting in an high-precision evaluation of the critical exponents and of the central charge of the system, confirming that the quantum linear-zigzag transition belongs to the critical Ising model universality class. Quantum corrections to the classical phase diagram are computed, and the range of experimental parameters where quantum effects play a role is provided. These results show that structural instabilities of one-dimensional interacting atomic arrays can simulate quantum critical phenomena typical of ferromagnetic systems.

Nonconservative dynamics in long atomic wires

The effect of nonconservative current-induced forces on the ions in a defect-free metallic nanowire is investigated using both steady-state calculations and dynamical simulations. Nonconservative forces were found to have a major influence on the ion dynamics in these systems, but their role in increasing the kinetic energy of the ions decreases with increasing system length. The results illustrate the importance of nonconservative effects in short nanowires and the scaling of these effects with system size. The dependence on bias and ion mass can be understood with the help of a simple pen and paper model. This material highlights the benefit of simple preliminary steady-state calculations in anticipating aspects of brute-force dynamical simulations, and provides rule of thumb criteria for the design of stable quantum wires.

Comet C/2012 S1 (ISON): Observations of the Dust Grains from SOFIA and of the Atomic Gas from NSO Dunn and McMath-Pierce Solar Telescopes (Invited)

Comet C/2012 S1 (ISON) is unique in that it is a dynamically new comet derived from the Oort cloud reservoir of comets with a sun-grazing orbit. Infrared (IR) and visible wavelength observing campaigns were planned on NASA's Stratospheric Observatory For Infrared Astronomy (SOFIA) and on National Solar Observatory Dunn (DST) and McMath-Pierce Solar Telescopes, respectively. We highlight our early results. SOFIA (+FORCAST [1]) mid- to far-IR images and spectroscopy (~5-35 μm) of the dust in the coma of ISON are to be obtained by the ISON-SOFIA Team during a flight window 2013 Oct 21-23 UT (r_h≈1.18 AU). Dust characteristics, identified through the 10 μm silicate emission feature and its strength [2], as well as spectral features from cometary crystalline silicates (Forsterite) at 11.05-11.2 μm, and near 16, 19, 23.5, 27.5, and 33 μm are compared with other Oort cloud comets that span the range of small and/or highly porous grains (e.g., C/1995 O1 (Hale-Bopp) [3,4,5] and C/2001 Q4 (NEAT) [6]) to large and/or compact grains (e.g., C/2007 N4 (Lulin) [7] and C/2006 P1 (McNaught) [8]). Measurement of the crystalline peaks in contrast to the broad 10 and 20 μm amorphous silicate features yields the cometary silicate crystalline mass fraction [9], which is a benchmark for radial transport in our protoplanetary disk [10]. The central wavelength positions, relative intensities, and feature asymmetries for the crystalline peaks may constrain the shapes of the crystals [11]. Only SOFIA can look for cometary organics in the 5-8 μm region. Spatially resolved measurements of atoms and simple molecules from when comet ISON is near the Sun (r_h<0.4 AU, near Nov-20--Dec-03 UT) were proposed for by the ISON-DST Team. Comet ISON is the first comet since comet Ikeya-Seki (1965f) [12,13] suitable for studying the alkalai metals Na and K and the atoms specifically attributed to dust grains including Mg, Si, Fe, as well as Ca. DST's Horizontal Grating Spectrometer (HGS) measures 4 settings: Na I, K, C2 to sample cometary organics (along with Mg I), and [O I] as a proxy for activity from water [14] (along with Si I and Fe I). State-of-the-art instruments that will also be employed include IBIS [15], which is a Fabry-Perot spectral imaging system that concurrently measures lines of Na, K, Ca II, or Fe, and ROSA (CSUN/QUB) [16], which is a rapid imager that simultaneously monitors Ca II or CN. From McMath-Pierce, the Solar-Stellar Spectrograph also will target ISON (320-900 nm, R~21,000, r_h

A density functional theory study of CO and atomic oxygen chemisorption on Pt(111)

Ab initio total energy calculations within a density functional theory framework have been performed for CO and atomic oxygen chemisorbed on the Pt(111) surface. Optimised geometries and chemisorption energies for CO and O on four high-symmetry sites, namely the top, bridge, fee hollow and hcp hollow sites, are presented, the coverage in all cases being 0.25 ML. The differences in CO adsorption energies between these sites are found to be small, suggesting that the potential energy surface for CO diffusion across Pt(111) is relatively flat. The 5 sigma and 2 pi molecular orbitals of CO are found to contribute to bonding with the metal. Some mixing of the 4 sigma and 1 pi molecular orbitals with metal states is also observed. For atomic oxygen, the most stable adsorption site is found to be the fee hollow site, followed in decreasing order of stability by the hcp hollow and bridge sites, with the top site being the least stable. The differences in chemisorption energies between sites for oxygen are larger than in the case of CO, suggesting a higher barrier to diffusion for atomic oxygen. The co-adsorption of CO and O has also been investigated. Calculated chemisorption energies for CO on an O/fcc-precovered surface show that of the available chemisorption sites, the top site at the oxygen atom's next-nearest neighbour surface metal atom is the most stable, with the other four sites calculated bring at least 0.29 eV less stable. The trend of CO site stability in the coadsorption system is explained in terms of a 'bonding competition' model. (C) 2000 Elsevier Science B.V. All rights reserved.

BUCKLED LAYER STRUCTURE FOR ATOMIC ADSORPTION ON W(100) - THE (ROOT-2X-ROOT-2)R45-DEGREES NITROGEN STRUCTURE FROM ATLEED

Adsorption of 0.5 monolayer of N adatoms on W{100} results in a sharp (root 2 X root 2)R45 degrees LEED pattern. The only previous quantitative LEED study of this system gave a simple overlayer model with a Pendry R-factor of 0.55. An exhaustive search has been made of possible structures, including a novel vacancy reconstruction, displacive reconstructions and underlayer adsorption. From this work a new overlayer structure is derived with an R(p) value of 0.22, displaying a considerable buckling of 0.27 +/- 0.05 Angstrom within the second W layer and consequently involving large changes in the interlayer spacings of the surface. The N adatom is pseudo-five-fold coordinated to the W surface, bonding to a second-layer W atom with a nearest-neighbour bond length of 2.13 Angstrom and with the four next-nearest-neighbour W atoms in the surface plane at 2.27 Angstrom. The structure does not resolve the work function anomaly observed on this surface.