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.

Line-driven Disk Winds in Activ Galactic Nuclei: The Critical Importance of Ionization and Radiative Transfer

Accretion disk winds are thought to produce many of the characteristic features seen in the spectra of active galactic nuclei (AGNs) and quasi-stellar objects (QSOs). These outflows also represent a natural form of feedback between the central supermassive black hole and its host galaxy. The mechanism for driving this mass loss remains unknown, although radiation pressure mediated by spectral lines is a leading candidate. Here, we calculate the ionization state of, and emergent spectra for, the hydrodynamic simulation of a line-driven disk wind previously presented by Proga & Kallman. To achieve this, we carry out a comprehensive Monte Carlo simulation of the radiative transfer through, and energy exchange within, the predicted outflow. We find that the wind is much more ionized than originally estimated. This is in part because it is much more difficult to shield any wind regions effectively when the outflow itself is allowed to reprocess and redirect ionizing photons. As a result, the calculated spectrum that would be observed from this particular outflow solution would not contain the ultraviolet spectral lines that are observed in many AGN/QSOs. Furthermore, the wind is so highly ionized that line driving would not actually be efficient. This does not necessarily mean that line-driven winds are not viable. However, our work does illustrate that in order to arrive at a self-consistent model of line-driven disk winds in AGN/QSO, it will be critical to include a more detailed treatment of radiative transfer and ionization in the next generation of hydrodynamic simulations.

Tidally Distorted Exoplanets: Density Corrections for Short-period Hot-Jupiters Based Solely on Observable Parameters

The close proximity of short-period hot-Jupiters to their parent star means they are subject to extreme tidal forces. This has a profound effect on their structure and, as a result, density measurements that assume that the planet is spherical can be incorrect. We have simulated the tidally distorted surface for 34 known short-period hot-Jupiters, assuming surfaces of constant gravitational equipotential for the planet, and the resulting densities have been calculated based only on observed parameters of the exoplanet systems. Comparing these results to the density values, assuming the planets are spherical, shows that there is an appreciable change in the measured density for planets with very short periods (typically less than two days). For one of the shortest-period systems, WASP-19b, we determine a decrease in bulk density of 12% from the spherical case and, for the majority of systems in this study, this value is in the range of 1%-5%. On the other hand, we also find cases where the distortion is negligible (relative to the measurement errors on the planetary parameters) even in the cases of some very short period systems, depending on the mass ratio and planetary radius. For high-density gas planets requiring apparently anomalously large core masses, density corrections due to tidal deformation could become important for the shortest-period systems.

INTERPLAY OF CAF2 SUBSTRATE ROUGHENING AND FILM DEPOSITION RATE ON RADIATIVE SCATTERING OF SURFACE-PLASMON POLARITONS FROM AG FILMS

The radiative decay of surface plasmon polaritons has been investigated in an attempt to characterize the surface roughness of Ag films prepared under different conditions. The polaritons were excited by the method of attenuated total reflection of light. The films were deposited on the face of a 60-degrees BK-7 glass prism at a rate that was deliberately fixed in two different ranges (centred on 0.1 and 10 nm s-1) and in some cases a CaF2 underlayer was used to roughen the film surfaces. The intensity of the scattered light emitted from the opposite face of the films was measured as a function of direction for each using the same sensitivity scale and was correlated with the preparation of the film. It was found that on nominally smooth substrates fast-deposited thinner films give out more light and are deduced to have greater short wavelength (300-600 nm) roughness amplitude. There is also evidence for long wavelenth (7 mum) periodic roughness due to the prism substrate itself. On CaF2 roughened surfaces the light output from the films is further increased and the peak intensity is backward directed with respect to the exciting laser beam direction. Here roughness on a lateral scale of 350 nm is responsible. Also, elastic scattering of surface plasmon polaritons at grain boundaries reduces the light output from fast deposited, small grain, films on CaF2 roughened surfaces. Overall, a consistent picture of roughness induced radiative polariton decay emerges for all cases studied.

Scaling of magneto-quantum-radiative hydrodynamic equations: From laser-produced plasmas to astrophysics

We introduce the equations of magneto-quantum-radiative hydrodynamics. By rewriting them in a dimensionless form, we obtain a set of parameters that describe scale-dependent ratios of characteristic hydrodynamic quantities. We discuss how these dimensionless parameters relate to the scaling between astrophysical observations and laboratory experiments.

Line-driven radiative outflows in luminous quasars

An analysis of ≃19 500 narrow (≲200 km s^-1) CIV λλ1548.2,1550.8 absorbers in ≃34 000 Sloan Digital Sky Survey quasar spectra is presented. The statistics of the number of absorbers as a function of outflow velocity shows that in approximately two-thirds of outflows, with multiple C IV absorbers present, absorbers are line-locked at the 500 km s^-1 velocity separation of the C IV absorber doublet; appearing as 'triplets' in the quasar spectra. Line-locking is an observational signature of radiative line-driving in outflowing material, where the successive shielding of 'clouds' of material in the outflow locks the clouds together in outflow velocity. Line-locked absorbers are seen in both broad absorption line (BAL) quasars and non-BAL quasars with comparable frequencies and with velocities out to at least 20 000 km s^-1. There are no detectable differences in the absorber properties and the dust content of single C IV doublets and line-locked C IV doublets. The gas associated with both single and line-locked CIV absorption systems includes material with a wide range of ionization potential (14-138 eV). Both single and line-locked CIV absorber systems show strong systematic trends in their ionization as a function of outflow velocity, with ionization decreasing rapidly with increasing outflow velocity. Initial simulations, employing CLOUDY, demonstrate that a rich spectrum of line-locked signals at various velocities may be expected due to significant opacities from resonance lines of Li-, He- and H-like ions of O, C and N, along with contributions from He II and HI resonance lines. The simulations confirm that line-driving can be the dominant acceleration mechanism for clouds with N(H I) ≃ 10¹⁹ cm^-2.

Ultracold, radiative charge transfer in hybrid Yb ion – Rb atom traps

Ultracold hybrid ion–atom traps offer the possibility of microscopic manipulation of quantum coherences in the gas using the ion as a probe. However, inelastic processes, particularly charge transfer can be a significant process of ion loss and has been measured experimentally for the ${\rm Y}{{{\rm b}}^{+}}$ ion immersed in a Rb vapour. We use first-principles quantum chemistry codes to obtain the potential energy curves and dipole moments for the lowest-lying energy states of this complex. Calculations for the radiative decay processes cross sections and rate coefficients are presented for the total decay processes; ${\rm Y}{{{\rm b}}^{+}}(6{\rm s}{{\;}^{2}}{\rm S})+{\rm Rb}(5{\rm s}{{\;}^{2}}{\rm S})\to {\rm Yb}(6{{{\rm s}}^{2}}{{\;}^{1}}{\rm S})+{\rm R}{{{\rm b}}^{+}}(4{{{\rm p}}^{6}}{{\;}^{1}}{\rm S})+h\nu $ and ${\rm Y}{{{\rm b}}^{+}}(6{\rm s}{{\;}^{2}}{\rm S})+{\rm Rb}(5{\rm s}{{\;}^{2}}{\rm S})\to {\rm YbR}{{{\rm b}}^{+}}({{X}^{1}}{{\Sigma }^{+}})+h\nu $. Comparing the semi-classical Langevin approximation with the quantum approach, we find it provides a very good estimate of the background at higher energies. The results demonstrate that radiative decay mechanisms are important over the energy and temperature region considered. In fact, the Langevin process of ion–atom collisions dominates cold ion–atom collisions. For spin-dependent processes [1] the anisotropic magnetic dipole–dipole interaction and the second-order spin–orbit coupling can play important roles, inducing coupling between the spin and the orbital motion. They measured the spin-relaxing collision rate to be approximately five orders of magnitude higher than the charge-exchange collision rate [1]. Regarding the measured radiative charge transfer collision rate, we find that our calculation is in very good agreement with experiment and with previous calculations. Nonetheless, we find no broad resonances features that might underly a strong isotope effect. In conclusion, we find, in agreement with previous theory that the isotope anomaly observed in experiment remains an open question.

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.

Marine or estuarine radiocarbon reservoir corrections for mollusks? A case study from a medieval site in the south of England

Mollusk shells are frequently radiocarbon dated and provide reliable calibrated age ranges when the regional marine reservoir correction is well-established. For mollusks from an estuarine environment the reservoir correction may be significantly different than the regional marine reservoir correction due to the input of bedrock or soil derived carbonates. Some mollusk species such as oysters are tolerant of a significant range of salinities which makes it difficult to determine which reservoir correction is appropriate. A case study is presented of an anomalous radiocarbon age for an oyster shell paint dish found in the fabric of the ruined nave walls of St Mary's Church, Shoreham-by-Sea, West Sussex, England. Stable isotopes (delta O-18 and delta C-13) were used to establish the type of environment in which the oyster had lived. Paired marine and terrestrial samples from a nearby medieval site were radiocarbon dated to provide an appropriate reservoir correction.

GRADIENT CORRECTIONS IN DENSITY-FUNCTIONAL THEORY CALCULATIONS FOR SURFACES - CO ON PD(110)

Ab initio total energy calculations have been performed for CO chemisorption on Pd(110). Local density approximation (LDA) calculations yield chemisorption energies which are significantly higher than experimental values but inclusion of the generalised gradient approximation (GGA) gives better agreement. In general, sites with higher coordination of the adsorbate to surface atoms lead to a larger degree of overbinding with LDA, and give larger corrections with GGA. The reason is discussed using a first-order perturbation approximation. It is concluded that this may be a general failure of LDA for chemisorption energy calculations. This conclusion may be extended to many surface calculations, such as potential energy surfaces for diffusion.

Radiative hydrodynamic modelling and observations of the X-class solar flare on 2011 March 9

Aims. We investigated the response of the solar atmosphere to non-thermal electron beam heating using the radiative transfer and hydrodynamics modelling code RADYN. The temporal evolution of the parameters that describe the non-thermal electron energy distribution were derived from hard X-ray observations of a particular flare, and we compared the modelled and observed parameters.

Methods. The evolution of the non-thermal electron beam parameters during the X1.5 solar flare on 2011 March 9 were obtained from analysis of RHESSI X-ray spectra. The RADYN flare model was allowed to evolve for 110 s, after which the electron beam heating was ended, and was then allowed to continue evolving for a further 300 s. The modelled flare parameters were compared to the observed parameters determined from extreme-ultraviolet spectroscopy.

Results. The model produced a hotter and denser flare loop than that observed and also cooled more rapidly, suggesting that additional energy input in the decay phase of the flare is required. In the explosive evaporation phase a region of high-density cool material propagated upward through the corona. This material underwent a rapid increase in temperature as it was unable to radiate away all of the energy deposited across it by the non-thermal electron beam and via thermal conduction. A narrow and high-density (ne ≤ 1015 cm-3) region at the base of the flare transition region was the source of optical line emission in the model atmosphere. The collision-stopping depth of electrons was calculated throughout the evolution of the flare, and it was found that the compression of the lower atmosphere may permit electrons to penetrate farther into a flaring atmosphere compared to a quiet Sun atmosphere.

Radiative rates for E1, E2, M1, and M2 transitions in the Br-like ions Sr IV, Y V, Zr VI, Nb VII, and Mo VIII

Energies and lifetimes are reported for the lowest 375 levels of five Br-like ions, namely SrIV, YV, ZrVI, NbVII, and MoVIII, mostly belonging to the 4s²4p⁵, 4s²4p⁴4ℓ, 4s4p⁶, 4s²4p⁴5ℓ, 4s²4p³4d², 4s4p⁵4ℓ, and 4s4p⁵5ℓ configurations. Extensive configuration interaction has been included and the general-purpose relativistic atomic structure package (grasp) has been adopted for the calculations. Additionally, radiative rates are listed among these levels for all E1, E2, M1, and M2 transitions. From a comparison with the measurements, the majority of our energy levels are assessed to be accurate to better than 2%, although discrepancies between theory and experiment for a few are up to 6%. An accuracy assessment of the calculated radiative rates (and lifetimes) is more difficult, because no prior results exist for these ions.

Energy levels, radiative rates and lifetimes for transitions in Br-like ions with 38 ≤ Z ≤ 42

Energy levels and radiative rates for transitions in five Br-like ions (Sr IV, Y V, Zr VI, Nb VII and Mo VIII) are calculated with the general-purpose relativistic atomic structure package (GRASP). Extensive configuration interaction has been included and results are presented among the lowest 31 levels of the 4s²4p⁵, 4s²4p⁴4d and 4s4p⁶ configurations. Lifetimes for these levels have also been determined, although unfortunately no measurements are available with which to compare. However, recently theoretical results have been reported by Singh et al (2013 Phys. Scr. 88 035301) using the same GRASP code. But their reported data for radiative rates and lifetimes cannot be reproduced and show discrepancies of up to five orders of magnitude with the present calculations.

Energy levels, radiative rates and electron impact excitation rates for transitions in Be-like Cl XIV, K XVI and Ge XXIX

Results for energy levels, radiative rates and electron impact excitation (effective) collision strengths for transitions in Be-like Cl XIV, K XVI and Ge XXIX are reported. For the calculations of energy levels and radiative rates the general-purpose relativistic atomic structure package is adopted, while for determining the collision strengths and subsequently the excitation rates, the Dirac atomic R-matrix code is used. Oscillator strengths, radiative rates and line strengths are listed for all E1, E2, M1 and M2 transitions among the lowest 98 levels of the n ≤ 4 configurations. Furthermore, lifetimes are provided for all levels and comparisons made with available theoretical and experimental results. Resonances in the collision strengths are resolved in a fine energy mesh and averaged over a Maxwellian velocity distribution to obtain the effective collision strengths. Results obtained are listed over a wide temperature range up to 10^7.8 K, depending on the ion.

Energy levels, radiative rates, and lifetimes for transitions in W LVIII

Energy levels and radiative rates are reported for transitions in Cl-like W LVIII. Configuration interaction (CI) has been included among 44 configurations (generating 4978 levels) over a wide energy range up to 363 Ryd, and the general-purpose relativistic atomic structure package (grasp) adopted for the calculations. Since no other results of comparable complexity are available, calculations have also been performed with the flexible atomic code (fac), which help in assessing the accuracy of our results. Energies are listed for the lowest 400 levels (with energies up to ~98 Ryd), which mainly belong to the 3s²3p⁵, 3s3p⁶, 3s²3p⁴3d, 3s²3p³3d², 3s3p⁴3d², 3s²3p²3d³, and 3p⁶3d configurations, and radiative rates are provided for four types of transitions, i.e.E1, E2, M1, and M2. Our energy levels are assessed to be accurate to better than 0.5%, whereas radiative rates (and lifetimes) should be accurate to better than 20% for a majority of the strong transitions.