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.

SOLVING THE INHOMOGENEOUS SCHRODINGER-EQUATION ON A BIDIRECTIONAL PIPELINE OF TRANSPUTERS

In this paper we describe the design of a parallel solution of the inhomogeneous Schrodinger equation, which arises in the construction of continuum orbitals in the R-matrix theory of atomic continuum processes. A prototype system is described which has been programmed in occam2 and implemented on a bi-directional pipeline of transputers. Some timing results for the prototype system are presented, and the development of a full production system is discussed.

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.

Oscillator strengths and transition probabilities for the W xlv ion

In this paper we present oscillator strengths and transition probabilities for W xlv transitions between levels arising from configurations 3d¹⁰4s²,4p²,4d², 3d¹⁰4k4l (k = s,p,d,f and l = p,d,f), 3d⁹4s²4l (l = p,d,f) and 3d⁹4s4p². The model used to calculate these contained all configurations which can be constructed from the available orbitals (up to n = 4), with either a 3d¹⁰ or 3d⁹ core. The calculations were performed with the configuration interaction CIV3 program with the inclusion of relativistic effects achieved through the use of the Breit-Pauli approximation. We compare our ab initio energy levels, oscillator strengths and transition rates with other experimental and theoretical values available in the literature. There is generally good agreement when only levels with 3d¹⁰ cores are considered. The literature is sparse for levels in which the 3d-subshell is opened: for the majority of the fine-structure lines considered, there is either no comparison data available or substantial differences are found. This paper also investigates how the inclusion of relativistic effects can result in a significant redistribution of the oscillator strength from the LS calculations.

Radiative charge transfer in cold and ultracold sulphur atoms colliding with protons

Radiative decay processes at cold and ultra cold temperatures for sulfur atoms colliding with protons are investigated. The MOLPRO quantum chemistry suite of codes was used to obtain accurate potential energies and transition dipole moments, as a function of internuclear distance, between low-lying states of the SH+ molecular cation. A multi-reference configuration-interaction approximation together with the Davidson correction is used to determine the potential energy curves and transition dipole moments, between the states of interest, where the molecular orbitals are obtained from state-averaged multi-configuration-self-consistent field calculations. The collision problem is solved approximately using an optical potential method to obtain radiative loss, and a fully two-channel quantum approach for radiative charge transfer. Cross sections and rate coefficients are determined for the first time for temperatures ranging from 10 μK up to 10 000 K. Results are obtained for all isotopes of sulfur, colliding with H+ and D+ ions and comparison is made to a number of other collision systems.

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

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.

Coherent and stochastic contributions of compound resonances in atomic processes: Electron recombination, photoionization, and scattering

In open-shell atoms and ions, processes such as photoionization, combination (Raman) scattering, electron scattering, and recombination are often mediated by many-electron compound resonances. We show that their interference (neglected in the independent-resonance approximation) leads to a coherent contribution, which determines the energy-averaged total cross sections of electron- and photon-induced reactions obtained using the optical theorem. In contrast, the partial cross sections (e.g., electron recombination or photon Raman scattering) are dominated by the stochastic contributions. Thus, the optical theorem provides a link between the stochastic and coherent contributions of the compound resonances. Similar conclusions are valid for reactions via compound states in molecules and nuclei.

Assessment of Atomic Data: Problems and Solutions

For the reliable analysis and modeling of astrophysical, laser-produced, and fusion plasmas, atomic data are required for a number of parameters, including energy levels, radiative rates, and electron impact excitation rates. Such data are desired for a range of elements (H to W) and their many ions. However, measurements of atomic data, mainly for radiative and excitation rates, are not feasible for many species, and therefore, calculations are needed. For some ions (such as of C, Fe, and Kr), there is a variety of calculations available in the literature, but often, they differ significantly from one another. Therefore, there is a great demand from the user community to have data "assessed" for accuracy so that they can be confidently applied to the modeling of plasmas. In this paper we highlight the difficulties in assessing atomic data and offer some solutions for improving the accuracy of calculated results.

Stout: Cloudy's Atomic and Molecular Database

We describe a new atomic and molecular database we developed for use in the spectral synthesis code Cloudy. The design of Stout is driven by the data needs of Cloudy, which simulates molecular, atomic, and ionized gas with kinetic temperatures and densities spanning the low-to high-density limits. The radiation field between photon energies 10−8 Ry and 100 MeV is considered, along with all atoms and ions of the lightest 30 elements, and ~102 molecules. For ease of maintenance, the data are stored in a format as close as possible to the original data sources. Few data sources include the full range of data we need. We describe how we fill in the gaps in the data or extrapolate rates beyond their tabulated range. We tabulate data sources both for the atomic spectroscopic parameters and for collision data for the next release of Cloudy. This is not intended as a review of the current status of atomic data, but rather a description of the features of the database which we will build upon.