From four independent calculations using three different theoretical approaches, rate coefficients for the formation of HeH$^+$ via the radiative association of He$^+$ and H were computed. Good agreement is found between our new calculations and prior results obtained two decades ago for kinetic temperatures between $\sim$800 and 20,000 K. This finding is inconsistent with a recent claim in the literature of a wide variation in published values and establishes the robustness of our knowledge of this process for the formation of HeH$^+$. The implications of the current results to the first detection of HeH$^+$ and its modeled abundance in the planetary nebula NGC 7027 are discussed.

This report provides detailed findings on the critical laboratory astrophysics data needs that are required to maximize the scientific return for NASA's current and near-term planned astrophysics missions. It also provides prioritized rankings on said laboratory astrophysics data, generally by waveband. The Report is based on community input gathered at the 2018 NASA Laboratory Astrophysics Workshop (LAW) from presentations, from discussions during workshop breakout sessions, and from other solicited input deemed appropriate by the Scientific Organizing Committee (SOC) obtained prior to and after the meeting. Hence, the Report is a direct reflection of the spirit and participant make-up of LAW 2018. The Report also outlines specific opportunities and threats facing NASA's Laboratory Astrophysics Program, and articulates concrete actions by which the Agency can capitalize on the opportunities and mitigate the challenges. The Report was prepared by the SOC, with help from some invited speakers, and input and review from community members.

Thanks to high-resolution and non-dispersive spectrometers onboard future X-ray missions such as XRISM and Athena, we are finally poised to answer important questions about the formation and evolution of galaxies and large-scale structure. However, we currently lack an adequate understanding of many atomic processes behind the spectral features we will soon observe. Large error bars on parameters as critical as transition energies and atomic cross sections can lead to unacceptable uncertainties in the calculations of e.g., elemental abundance, velocity, and temperature. Unless we address these issues, we risk limiting the full scientific potential of these missions. Laboratory astrophysics, which comprises theoretical and experimental studies of the underlying physics behind observable astrophysical processes, is therefore central to the success of these missions.

About two generations ago, a large part of AMO science was dominated by experimental high energy collision studies and perturbative theoretical methods. Since then, AMO science has undergone a transition and is now dominated by quantum, ultracold, and ultrafast studies. But in the process, the field has passed over the complexity that lies between these two extremes. Most of the Universe resides in this intermediate region. We put forward that the next frontier for AMO science is to explore the AMO complexity that describes most of the Cosmos.

Accurate photodissociation cross sections have been computed for transitions from the X $^1\Sigma^+$ ground electronic state of CS to six low-lying excited electronic states. New ab initio potential curves and transition dipole moment functions have been obtained for these computations using the multi-reference configuration interaction approach with the Davidson correction (MRCI+Q) and aug-cc-pV6Z basis sets. State-resolved cross sections have been computed for transitions from nearly the full range of rovibrational levels of the X $^1\Sigma^+$ state and for photon wavelengths ranging from 500 $\text{\AA}$ to threshold. Destruction of CS via predissociation in highly excited electronic states originating from the rovibrational ground state is found to be unimportant. Photodissociation cross sections are presented for temperatures in the range between 1000 and 10,000 K, where a Boltzmann distribution of initial rovibrational levels is assumed. Applications of the current computations to various astrophysical environments are briefly discussed focusing on photodissociation rates due to the standard interstellar and blackbody radiation fields.

We report the first full-dimensional potential energy surface (PES) and quantum mechanical close-coupling calculations for scattering of SiO due to H$_2$. The full-dimensional interaction potential surface was computed using the explicitly correlated coupled-cluster (CCSD(T)-F12b) method and fitted using an invariant polynomial approach. Pure rotational quenching cross sections from initial states $v_1=0$, $j_1$=1-5 of SiO in collision with H$_2$ are calculated for collision energies between 1.0 and 5000 cm$^{-1}$. State-to-state rotational rate coefficients are calculated at temperatures between 5 and 1000 K. The rotational rate coefficients of SiO with para-H$_2$ are compared with previous approximate results which were obtained using SiO-He PESs or scaled from SiO-He rate coefficients. Rovibrational state-to-state and total quenching cross sections and rate coefficients for initially excited SiO($v_1=1, j_1$=0 and 1) in collisions with para-H$_2$($v_2=0,j_2=0$) and ortho-H$_2$($v_2=0,j_2=1$) were also obtained. The application of the current collisional rate coefficients to astrophysics is briefly discussed.

This editorial introduces the J. Phys. B: Atomic, Molecular and Optical Physics Special Issue "Atomic and Molecular Processes in the Ultracold Regime, the Chemical Regime and Astrophysics" dedicated to Professor Alexander Dalgarno (1928-2015). After a brief biographical review, short summaries of the contributed papers and their relations to some of Prof. Dalgarno's work are given.

Charge exchange (CX) has emerged in X-ray emission modeling as a significant process that must be considered in many astrophysical environment--particularly comets. Comets host an interaction between solar wind ions and cometary neutrals to promote solar wind charge exchange (SWCX). X-ray observatories provide astronomers and astrophysicists with data for many X-ray emitting comets that are impossible to accurately model without reliable CX data. Here, we utilize a streamlined set of computer programs that incorporate the multi-channel Landau-Zener theory and a cascade model for X-ray emission to generate cross sections and X-ray line ratios for a variety of bare and non-bare ion single electron capture (SEC) collisions. Namely, we consider collisions between the solar wind constituent bare and H-like ions of C, N, O, Ne, Na, Mg, Al, and Si and the cometary neutrals H2O, CO, CO2, OH, and O. To exemplify the application of this data, we model the X-ray emission of Comet C/2000 WM1 (linear) using the CX package in SPEX and find excellent agreement with observations made with the XMM-Newton RGS detector. Our analyses show that the X-ray intensity is dominated by SWCX with H, while H2O plays a secondary role. This is the first time, to our knowledge, that CX cross sections have been implemented into a X-ray spectral fitting package to determine the H to H2O ratio in cometary atmospheres. The CX data sets are incorporated into the modeling packages SPEX and Kronos.

Accurate molecular data for the low-lying states of SiO are computed and used to calculate rate constants for radiative association of Si and O. Einstein A-coefficients are also calculated for transitions between all of the bound and quasi-bound levels for each molecular state. The radiative widths are used together with elastic tunneling widths to define effective radiative association rate constants which include both direct and indirect (inverse predissociation) formation processes. The indirect process is evaluated for two kinetic models which represent limiting cases for astrophysical environments. The first case scenario assumes an equilibrium distribution of quasi-bound states and would be applicable whenever collisional and/or radiative excitation mechanisms are able to maintain the population. The second case scenario assumes that no excitation mechanisms are available which corresponds to the limit of zero radiation temperature and zero atomic density. Rate constants for SiO formation in realistic astrophysical environments would presumably lie between these two limiting cases.

We describe the 2017 release of the spectral synthesis code Cloudy. A major development since the previous release has been exporting the atomic data into external data files. This greatly simplifies updates and maintenance of the data. Many large datasets have been incorporated with the result that we can now predict well over an order of magnitude more emission lines when all databases are fully used. The use of such large datasets is not realistic for most calculations due to the time and memory needs, and we describe the limited subset of data we use by default. Despite the fact that we now predict significantly more lines than the previous Cloudy release, this version is faster because of optimization of memory access patterns and other tuning. The size and use of the databases can easily be adjusted in the command-line interface. We give examples of the accuracy limits using small models, and the performance requirements of large complete models. We summarize several advances in the H- and He-like iso-electronic sequences. We use our complete collisional-radiative models of the ionization of these one and two-electron ions to establish the highest density for which the coronal or interstellar medium (ISM) approximation works, and the lowest density where Saha or local thermodynamic equilibrium can be assumed. The coronal approximation fails at surprisingly low densities for collisional ionization equilibrium but is valid to higher densities for photoionized gas clouds. Many other improvements to the physics have been made and are described. These include the treatment of isotropic continuum sources such as the cosmic microwave background (CMB) in the reported output, and the ability to follow the evolution of cooling non-equilibrium clouds.

Charge exchange is an important process for the modeling of X-ray spectra obtained by the Chandra, XMM-Newton, and Suzaku X-ray observatories, as well as the anticipated Astro-H mission. The understanding of the observed X-ray spectra produced by many astrophysical environments is hindered by the current incompleteness of available atomic and molecular data -- especially for charge exchange. Here, we implement a streamlined program set that applies quantum defect methods and the Landau-Zener theory to generate total, n-resolved, and nlS-resolved cross sections for any given projectile ion/ target charge exchange collision. Using this data in a cascade model for X-ray emission, theoretical spectra for such systems can be predicted. With these techniques, Fe25+ and Fe26+ charge exchange collisions with H, He, H2, N2, H2O, and CO are studied for single electron capture. These systems have been selected as they illustrate computational difficulties for high projectile charges. Further, Fe XXV and Fe XXVI emission lines have been detected in the Galactic center and Galactic ridge. Theoretical X-ray spectra for these collision systems are compared to experimental data generated by an electron beam ion trap study. Several l- distribution models have been tested for Fe25+ and Fe26+ single electron capture. Such analysis suggests that commonly used l-distribution models struggle to accurately reflect the true distribution of electron capture as understood by more advanced theoretical methods.

A five-dimensional coupled states (5D-CS) approximation is used to compute cross sections and rate coefficients for CO+H2 collisions. The 5D-CS calculations are benchmarked against accurate six-dimensional close-coupling (6D-CC) calculations for transitions between low-lying rovibrational states. Good agreement between the two formulations is found for collision energies greater than 10 cm-1. The 5D-CS approximation is then used to compute two separate databases which include highly excited states of CO that are beyond the practical limitations of the 6D-CC method. The first database assumes an internally frozen H2 molecule and allows rovibrational transitions for v < 5 and j < 30, where v and j are the vibrational and rotational quantum numbers of the initial state of the CO molecule. The second database allows H2 rotational transitions for initial CO states with v < 5 and j < 10. The two databases are in good agreement with each other for transitions that are common to both basis sets. Together they provide data for astrophysical models which were previously unavailable.

Our computation effort is primarily concentrated on support of current and future measurements being carried out at various synchrotron radiation facilities around the globe, and photodissociation computations for astrophysical applications. In our work we solve the Schr\"odinger or Dirac equation for the appropriate collision problem using the R-matrix or R-matrix with pseudo-states approach from first principles. The time dependent close-coupling (TDCC) method is also used in our work. A brief summary of the methodology and ongoing developments implemented in the R-matrix suite of Breit-Pauli and Dirac-Atomic R-matrix codes (DARC) is presented.

Gas-grain and gas-phase reactions dominate the formation of molecules in the interstellar medium (ISM). Gas-grain reactions require a substrate (e.g. a dust or ice grain) on which the reaction is able to occur. The formation of molecular hydrogen (H$_2$) in the ISM is the prototypical example of a gas-grain reaction. In these reactions, an atom of hydrogen will strike a surface, stick to it, and diffuse across it. When it encounters another adsorbed hydrogen atom, the two can react to form molecular hydrogen and then be ejected from the surface by the energy released in the reaction. We perform in-depth classical molecular dynamics (MD) simulations of hydrogen atoms interacting with an amorphous water-ice surface. This study focuses on the first step in the formation process; the sticking of the hydrogen atom to the substrate. We find that careful attention must be paid in dealing with the ambiguities in defining a sticking event. The technical definition of a sticking event will affect the computed sticking probabilities and coefficients. Here, using our new definition of a sticking event, we report sticking probabilities and sticking coefficients for nine different incident kinetic energies of hydrogen atoms [5 K - 400 K] across seven different temperatures of dust grains [10 K - 70 K]. We find that probabilities and coefficients vary both as a function of grain temperature and incident kinetic energy over the range of 0.99 - 0.22.

Recent X-ray observations of star-forming galaxies such as M82 have shown the Ly beta/Ly alpha line ratio of Ne X to be in excess of predictions for thermal electron impact excitation. Here we demonstrate that the observed line ratio may be due to charge exchange and can be used to constrain the ion kinetic energy to be <500 eV/u. This is accomplished by computing spectra and line ratios via a range of theoretical methods and comparing these to experiments with He over astrophysically relevant collision energies. The charge exchange emission spectra calculations were performed for Ne[10+] +H and Ne[10+] +He using widely applied approaches including the atomic orbital close coupling, classical trajectory Monte Carlo, and multichannel Landau- Zener (MCLZ) methods. A comparison of the results from these methods indicates that for the considered energy range and neutrals (H, He) the so-called "low-energy l-distribution" MCLZ method provides the most likely reliable predictions.

Photodissociation cross sections for the SH$^+$ radical are computed from all rovibrational (RV) levels of the ground electronic state X$~^3\Sigma^-$ for wavelengths from threshold to 500~\AA. The five electronic transitions, $2~ ^3\Sigma^- \leftarrow$ X$~^3\Sigma^-$, $3~ ^3\Sigma^- \leftarrow$ X$~^3\Sigma^-$, $A~ ^3\Pi \leftarrow$ X$~^3\Sigma^-$, $2~ ^3\Pi \leftarrow$ X$~^3\Sigma^-$, and $3~ ^3\Pi \leftarrow$ X$~^3\Sigma^-$, are treated with a fully quantum-mechanical two-state model, {i.e. no non-adiabatic coupling between excited states was included in our work.}. The photodissociation calculations incorporate adiabatic potentials and transition dipole moment functions computed in the multireference configuration interaction approach along with the Davidson correction (MRCI+Q), but adjusted to match available experimental molecular data and asymptotic atomic limits. Local thermodynamic equilibrium (LTE) photodissociation cross sections were computed which assume a Boltzmann distribution of RV levels in the X$~^3\Sigma^-$ molecular state of the SH$^+$ cation. The LTE cross sections are presented for temperatures in the range 1000-10,000~K. Applications of the current photodissociation cross sections to interstellar gas, photon-dominated regions, and stellar atmospheres are briefly discussed.

Motivated by recent claims of a compelling ~3.5 keV emission line from nearby galaxies and galaxy clusters, we investigate a novel plasma model incorporating a charge exchange component obtained from theoretical scattering calculations. Fitting this kind of component with a standard thermal model yields positive residuals around 3.5 keV, produced mostly by S XVI transitions from principal quantum numbers n > 8 to the ground. Such high-n states can only be populated by the charge exchange process. In this scenario, the observed 3.5 keV line flux in clusters can be naturally explained by an interaction in an effective volume of ~1 kpc^3 between a ~3 keV temperature plasma and cold dense clouds moving at a few hundred km/s. The S XVI lines at ~3.5 keV also provide a unique diagnostic of the charge exchange phenomenon in hot cosmic plasmas.

We present calculated rate coefficients for ro-vibrational transitions of CO in collisions with H atoms for a gas temperature range of 10 K $\leq T \leq$ 3000 K, based on the recent three-dimensional ab initio H-CO interaction potential of Song et al(2013). Rate coefficients for ro-vibrational $v=1,j=0-30 \rightarrow v'=0, j'$ transitions were obtained from scattering cross sections previously computed with the close-coupling method by Song et al(2015). Combining these with the rate coefficients for vibrational $v=1-5 \rightarrow v' < v$ quenching obtained with the infinite-order sudden approximation, we propose a new extrapolation scheme that yields the rate coefficients for ro-vibrational $v=2-5,j=0-30 \rightarrow v',j'$ de-excitation. Cross sections and rate coefficients for ro-vibrational $v=2, j=0-30 \rightarrow v'=1,j'$ transitions calculated with the close-coupling method confirm the effectiveness of this extrapolation scheme. Our calculated and extrapolated rates are very different from those that have been adopted in the modeling of many astrophysical environments. The current work provides the most comprehensive and accurate set of ro-vibrational de-excitation rate coefficients for the astrophysical modeling of the H-CO collision system. Application of the previously available and new data sets in astrophysical slab models shows that the line fluxes typically change by 20-70% in high temperature environments (800 K) with an H/H$_2$ ratio of 1; larger changes occur for lower temperatures.

Carbon monoxide is a simple molecule present in many astrophysical environments, and collisional excitation rate coefficients due to the dominant collision partners are necessary to accurately predict spectral line intensities and extract astrophysical parameters. We report new quantum scattering calculations for rotational deexcitation transitions of CO induced by H using the three-dimensional potential energy surface~(PES) of Song et al. (2015). State-to-state cross sections for collision energies from 10$^{-5}$ to 15,000~cm$^{-1}$ and rate coefficients for temperatures ranging from 1 to 3000~K are obtained for CO($v=0$, $j$) deexcitation from $j=1-45$ to all lower $j'$ levels, where $j$ is the rotational quantum number. Close-coupling and coupled-states calculations were performed in full-dimension for $j$=1-5, 10, 15, 20, 25, 30, 35, 40, and 45 while scaling approaches were used to estimate rate coefficients for all other intermediate rotational states. The current rate coefficients are compared with previous scattering results using earlier PESs. Astrophysical applications of the current results are briefly discussed.

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 2.8 K < T < 1e10 K 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 ~100 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.

Collisional excitation rate coefficients play an important role in the dynamics of energy transfer in the interstellar medium. In particular, accurate rotational excitation rates are needed to interpret microwave and infrared observations of the interstellar gas for nonlocal thermodynamic equilibrium line formation. Theoretical cross sections and rate coefficients for collisional deexcitation of rotationally excited HF in the vibrational ground state are reported. The quantum-mechanical close-coupling approach implemented in the nonreactive scattering code MOLSCAT was applied in the cross section and rate coefficient calculations on an accurate 2D HF-He potential energy surface. Estimates of rate coefficients for H and H$_2$ colliders were obtained from the HF-He collisional data with a reduced-potential scaling approach. The calculation of state-to-state rotational quenching cross sections for HF due to He with initial rotational levels up to $j=20$ were performed for kinetic energies from 10$^{-5}$ to 15000 cm$^{-1}$. State-to-state rate coefficients for temperatures between 0.1 and 3000 K are also presented. The comparison of the present results with previous work for lowly-excited rotational levels reveals significant differences. In estimating HF-H$_2$ rate coefficients, the reduced-potential method is found to be more reliable than the standard reduced-mass approach. The current state-to-state rate coefficient calculations are the most comprehensive to date for HF-He collisions. We attribute the differences between previously reported and our results to differences in the adopted interaction potential energy surfaces. The new He rate coefficients can be used in a variety of applications. The estimated H$_2$ and H collision rates can also augment the smaller datasets previously developed for H$_2$ and electrons.

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 (MRCI) 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 (MO's) are obtained from state-averaged multi configuration-self-consistent field (MCSCF) 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 $\mu$ 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.

Accurate rate coefficients for molecular vibrational transitions due to collisions with H$_2$, critical for interpreting infrared astronomical observations, are lacking for most molecules. Quantum calculations are the primary source of such data, but reliable values that consider all internal degrees of freedom of the collision complex have only been reported for H$_2$-H$_2$ due to the difficulty of the computations. Here we present essentially exact full-dimensional dynamics computations for rovibrational quenching of CO due to H$_2$ impact. Using a high-level six-dimensional potential surface, time-independent scattering calculations, within a full angular-momentum-coupling formulation, were performed for the deexcitation of vibrationally excited CO. Agreement with experimentally-determined results confirms the accuracy of the potential and scattering computations, representing the largest of such calculations performed to date. This investigation advances computational quantum dynamics studies representing initial steps toward obtaining CO-H$_2$ rovibrational quenching data needed for astrophysical modeling.

Observations of interstellar clouds that cast shadows in the soft X-ray background can be used to separate the background Galactic halo emission from the local emission due to solar wind charge exchange (SWCX) and/or the Local Bubble (LB). We present an XMM-Newton observation of a shadowing cloud, G225.60-66.40, that is sufficiently compact that the on- and off-shadow spectra can be extracted from a single field of view (unlike previous shadowing observations of the halo with CCD-resolution spectrometers, which consisted of separate on- and off-shadow pointings). We analyzed the spectra using a variety of foreground models: one representing LB emission, and two representing SWCX emission. We found that the resulting halo model parameters (temperature $T_h \approx 2 \times 10^6$ K, emission measure $E_h \approx 4 \times 10^{-3}$ cm$^{-6}$ pc) were not sensitive to the foreground model used. This is likely due to the relative faintness of the foreground emission in this observation. However, the data do favor the existence of a foreground. The halo parameters derived from this observation are in good agreement with those from previous shadowing observations, and from an XMM-Newton survey of the Galactic halo emission. This supports the conclusion that the latter results are not subject to systematic errors, and can confidently be used to test models of the halo emission.

Rate coefficients for collisional processes such as rotational and vibrational excitation are essential inputs in many astrophysical models. When rate coefficients are unknown, they are often estimated using known values from other systems. The most common example is to use He-collider rate coefficients to estimate values for other colliders, typically H$_2$, using scaling arguments based on the reduced mass of the collision system. This procedure is often justified by the assumption that the inelastic cross section is independent of the collider. Here we explore the validity of this approach focusing on rotational inelastic transitions for collisions of H, para-H$_2$, $^3$He, and $^4$He with CO in its vibrational ground state. We compare rate coefficients obtained via explicit calculations to those deduced by standard reduced-mass scaling. Not surprisingly, inelastic cross sections and rate coefficients are found to depend sensitively on both the reduced mass and the interaction potential energy surface. We demonstrate that standard reduced-mass scaling is not valid on physical and mathematical grounds, and as a consequence, the common approach of multiplying a rate coefficient for a molecule-He collision system by the constant factor of ~1.4 to estimate the rate coefficient for para-H$_2$ collisions is deemed unreliable. Furthermore, we test an alternative analytic scaling approach based on the strength of the interaction potential and the reduced mass of the collision systems. Any scaling approach, however, may be problematic when low-energy resonances are present; explicit calculations or measurements of rate coefficients are to be preferred.

Recent X-ray studies have shown that supernova shock models are unable to satisfactorily explain X-ray emission in the rim of the Cygnus Loop. In an attempt to account for this anomalously enhanced X-ray flux, we fit the region with a model including theoretical charge exchange (CX) data along with shock and background X-ray models. The model includes the CX collisions of $O^{8+}$, $O{7+}$, $N^{7+}$, $N^{6+}$, $C^{6+}$, and $C^{5+}$ with H with an energy of 1 keV/u (438 km/s). The observations reveal a strong emission feature near 0.7 keV that cannot fully be accounted for by a shock model, nor the current CX data. Inclusion of CX, specifically $O^{7+} + H$, does provide for a statistically significant improvement over a pure shock model.

We report rotational quenching cross sections and rate coefficients of HCl due to collisions with He. The close-coupling method and the coupled-states approximation are applied in quantum-mechanical scattering calculations of state-to-state cross sections for HCl with initial rotational levels up to $j=20$ for kinetic energies from 10$^{-5}$ to 15000 cm$^{-1}$. State-to-state rate coefficients for temperature between 0.1 and 3000 K are also presented. Comparison of the present rate coefficients with previous results reported in the literature for lowly-excited rotational levels shows reasonable agreement. Small differences are attributed to the differences in the interaction potential energy surfaces. The uncertainty in the computed cross sections and rate coefficients is estimated by varying the potential well depth. Applications of current results to astrophysical systems are also briefly discussed.

New quantum scattering calculations for rotational deexcitation transitions of CO induced by H collisions using two CO-H potential energy surfaces (PESs) from Shepler et al. (2007) are reported. State-to-state rate coefficients are computed for temperatures ranging from 1 to 3000 K for CO($v=0,j$) deexcitation from $j=1-5$ to all lower $j^\prime$ levels, with $j$ being the rotational quantum number. Different resonance structures in the cross sections are attributed to the differences in the anisotropy and the long-range van der Waals well depths of the two PESs. These differences affect rate coefficients at low temperatures and give an indication of the uncertainty of the results. Significant discrepancies are found between the current rate coefficients and previous results computed using earlier potentials, while the current results satisfy expected propensity rules. Astrophysical applications to modeling far infrared and submillimeter observations are briefly discussed.

This is a summary of the 2013 release of the plasma simulation code Cloudy. Cloudy models the ionization, chemical, and thermal state of material that may be exposed to an external radiation field or other source of heating, and predicts observables such as emission and absorption spectra. It works in terms of elementary processes, so is not limited to any particular temperature or density regime. This paper summarizes advances made since the last major review in 1998. Much of the recent development has emphasized dusty molecular environments, improvements to the ionization / chemistry solvers, and how atomic and molecular data are used. We present two types of simulations to demonstrate the capability of the code. We consider a molecular cloud irradiated by an X-ray source such as an Active Nucleus and show how treating EUV recombination lines and the full SED affects the observed spectrum. A second example illustrates the very wide range of particle and radiation density that can be considered.

Theoretical rotational quenching cross sections and rate coefficients of ortho- and para-H$_2$O due to collisions with He atoms are presented. The complete angular momentum close-coupling approach as well as the coupled-states approximation for angular momentum decoupling were applied to solve the scattering problem for a large range of rotationally-excited states of water. Results are obtained for quenching from initial levels 1$_{1,0}$, 2$_{1,2}$, 2$_{2,1}$, 3$_{0,3}$, 3$_{1,2}$, 3$_{2,1}$, 4$_{1,4}$, 3$_{3,0}$, and 4$_{2,3}$ of ortho-H$_2$O and from initial levels 1$_{1,1}$, 2$_{0,2}$, 2$_{1,1}$, 2$_{2,0}$, 3$_{1,3}$, 3$_{2,2}$, 4$_{0,4}$, 4$_{1,3}$, and 3$_{3,1}$ of para-H$_2$O for kinetic energies from 10$^{-5}$ to 10$^4$ cm$^{-1}$. State-to-state and total deexcitation cross sections and rate coefficients for temperatures between 0.1 and 3000 K are reported. The present state-to-state rate coefficients are found to be in good agreement with previous results obtained by Green and coworkers at high temperatures, but significant discrepancies are obtained at lower temperatures likely due to differences in the adopted potential energy surfaces. Astrophysical applications of the current rate coefficients are briefly discussed.

A group of 72 members of the astrochemistry and spectroscopy communities assembled at Emory University in October, 2012, to discuss the state-of-the art in molecular laboratory spectroscopic techniques used to study molecules of astrophysical interest, and highlight recent successes in laboratory efforts that complement new astronomical observations. The format of the meeting facilitated detailed discussion of the needs for collecting and interpreting molecular spectroscopic information to support observational astrophysics in the era of Far-IR astronomy. We present the findings and recommendations from this workshop in this report. The participants in this workshop identified the development of new laboratory capabilities that offer rapid, broadband, high-resolution, high-sensitivity measurements of molecular spectra in the frequency ranges that overlap with new observational facilities as the top priority for far-infrared laboratory astrophysics initiatives over the next 5-10 years. An additional major need is the development of suitable analysis tools for broadband molecular spectra from both the laboratory and space; though such efforts should be conducted within the framework of larger collaborative funding structures, rather than through single-PI efforts, so as to provide tools that are generally applicable rather than highly specialized. Additional efforts in theory, computation, and modeling that support these goals are also important for advancing this field.

Neutron(n)-capture elements (atomic number Z>30), which can be produced in planetary nebula (PN) progenitor stars via s-process nucleosynthesis, have been detected in nearly 100 PNe. This demonstrates that nebular spectroscopy is a potentially powerful tool for studying the production and chemical evolution of trans-iron elements. However, significant challenges must be addressed before this goal can be achieved. One of the most substantial hurdles is the lack of atomic data for n-capture elements, particularly that needed to solve for their ionization equilibrium (and hence to convert ionic abundances to elemental abundances). To address this need, we have computed photoionization cross sections and radiative and dielectronic recombination rate coefficients for the first six ions of Se and Kr. The calculations were benchmarked against experimental photoionization cross section measurements. In addition, we computed charge transfer (CT) rate coefficients for ions of six n-capture elements. These efforts will enable the accurate determination of nebular Se and Kr abundances, allowing robust investigations of s-process enrichments in PNe.

We present total and final-state resolved charge transfer (CT) rate coefficients for low-charge Ge, Se, Br, Kr, Rb, and Xe ions reacting with neutral hydrogen over the temperature range 10^2--10^6 K. Each of these elements has been detected in ionized astrophysical nebulae, particularly planetary nebulae. CT rate coefficients are a key ingredient for the ionization equilibrium solutions needed to determine total elemental abundances from those of the observed ions. A multi-channel Landau Zener approach was used to compute rate coefficients for projectile ions with charges q=2-5, and for singly-charged ions the Demkov approximation was utilized. Our results for five-times ionized species are lower limits, due to the incompleteness of level energies in the NIST database. In addition, we computed rate coefficients for charge transfer ionization reactions between the neutral species of the above six elements and ionized hydrogen. The resulting total and state-resolved CT rate coefficients are tabulated and available at the CDS. In tandem with our concurrent investigations of other important atomic processes in photoionized nebulae, this work will enable robust investigations of neutron-capture element abundances and nucleosynthesis via nebular spectroscopy.

Accurate photodissociation cross sections have been obtained for the A-X electronic transition of HeH+ using ab initio potential curves and dipole transition moments. Partial cross sections have been evaluated for all rotational transitions from the vibrational levels v"=0-11 and over the entire accessible wavelength range 100-1129 Angstrom. Assuming a Boltzmann distribution of the rovibrational levels of the X state, photodissociation cross sections are presented for temperatures between 500 and 12,000 K. A similar set of calculations was performed for the pure rovibrational photodissociation in the X-X electronic ground state, but covering photon wavelengths into the far infrared. Applications of the cross sections to the destruction of HeH+in the early Universe and in UV-irradiated environments such as primordial halos and protoplanetary disks are briefly discussed.

Rate coefficients for state-to-state rotational transitions in CO induced by both para- and ortho-H$_2$ collisions are presented. The results were obtained using the close-coupling method and the coupled-states approximation, with the CO-H$_2$ interaction potential of Jankowski & Szalewicz (2005). Rate coefficients are presented for temperatures between 1 and 3000 K, and for CO($v=0,j$) quenching from $j=1-40$ to all lower $j^\prime$ levels. Comparisons with previous calculations using an earlier potential show some discrepancies, especially at low temperatures and for rotational transitions involving large $|\Delta j|$. The differences in the well depths of the van der Waals interactions in the two potential surfaces lead to different resonance structures in the energy dependence of the cross sections which influence the low temperature rate coefficients. Applications to far infrared observations of astrophysical environments are briefly discussed.

Experimental and theoretical state-selective X-ray spectra resulting from single-electron capture in charge exchange (CX) collisions of Ne^10+ with He, Ne, and Ar are presented for a collision velocity of 933 km s^-1 (4.54 keV nucleon^-1), comparable to the highest velocity components of the fast solar wind. The experimental spectra were obtained by detecting scattered projectiles, target recoil ions, and X-rays in coincidence; with simultaneous determination of the recoil ion momenta. Use and interpretation of these spectra are free from the complications of non-coincident total X-ray measurements that do not differentiate between the primary reaction channels. The spectra offer the opportunity to test critically the ability of CX theories to describe such interactions at the quantum orbital angular momentum level of the final projectile ion. To this end, new classical trajectory Monte Carlo calculations are compared here with the measurements. The current work demonstrates that modeling of cometary, heliospheric, planetary, and laboratory X-ray emission based on approximate state-selective CX models may result in erroneous conclusions and deductions of relevant parameters.

A new and more comprehensive model of charge-exchange induced X-ray emission, due to ions precipitating into the Jovian atmosphere near the poles, has been used to analyze spectral observations made by the Chandra X-ray Observatory. The model includes for the first time carbon ions, in addition to the oxygen and sulfur ions previously considered, in order to account for possible ion origins from both the solar wind and the Jovian magnetosphere. By comparing the model spectra with newly reprocessed Chandra observations, we conclude that carbon ion emission provides a negligible contribution, suggesting that solar wind ions are not responsible for the observed polar X-rays. In addition, results of the model fits to observations support the previously estimated seeding kinetic energies of the precipitating ions (~0.7-2 MeV/u), but infer a different relative sulfur to oxygen abundance ratio for these Chandra observations.

The Orion Bar is one of the nearest and best-studied photodissociation or photon-dominated regions (PDRs). Observations reveal the presence of H2 lines from vibrationally or rotationally excited upper levels that suggest warm gas temperatures (400 to 700 K). However, standard models of PDRs are unable to reproduce such warm rotational temperatures. In this paper we attempt to explain these observations with new comprehensive models which extend from the H+ region through the Bar and include the magnetic field in the equation of state. We adopt the model parameters from our previous paper which successfully reproduced a wide variety of spectral observations across the Bar. In this model the local cosmic-ray density is enhanced above the galactic background, as is the magnetic field, and which increases the cosmic-ray heating elevating the temperature in the molecular region. The pressure is further enhanced above the gas pressure in the H+ region by the momentum transferred from the absorbed starlight. Here we investigate whether the observed H2 lines can be reproduced with standard assumptions concerning the grain photoelectric emission. We also explore the effects due to the inclusion of recently computed H2 + H2, H2 + H and H2 + He collisional rate coefficients.

Rate coefficients for rotational transitions in H_2 induced by H_2 impact are presented. Extensive quantum mechanical coupled-channel calculations based on a recently published (H_2)_2 potential energy surface were performed. The potential energy surface used here is presumed to be more reliable than surfaces used in previous work. Rotational transition cross sections with initial levels J <= 8 were computed for collision energies ranging between 0.0001 and 2.5 eV, and the corresponding rate coefficients were calculated for the temperature range 2 < T <10,000 K. In general, agreement with earlier calculations, which were limited to 100-6000 K, is good though discrepancies are found at the lowest and highest temperatures. Low-density-limit cooling functions due to para- and ortho-H_2 collisions are obtained from the collisional rate coefficients. Implications of the new results for non-thermal H_2 rotational distributions in molecular regions are also investigated.

Recent laboratory measurements for the S^+2 + H2 reaction find a total rate coefficient significantly larger than previously used in theoretical models of X-ray dominated regions (XDRs). While the branching ratio of the products is unknown, one energetically possible route leads to the SH+ molecule, a known XDR diagnostic. In this work, we study the effects of S^+2 on the formation of SH+ and the destruction of S^+2 in XDRs. We find the predicted SH+ column density for molecular gas surrounding an Active Galactic Nucleus (AGN) increases by as much as 2 dex. As long as the branching ratio for S^+2 + H2 -> SH+ + H+ exceeds a few percent, doubly ionized chemistry will be the dominant pathway to SH+, which then initiates the formation of other sulfur-bearing molecules. We also find that the high rate of S^+2 + H2 efficiently destroys S^+2 once H2 forms, while the S^+2 abundance remains high in the atomic hydrogen region. We discuss the possible consequences of S^+2 in the atomic hydrogen region on mid-infrared diagnostics. The enhanced SH+ abundance has important implications in the study of XDRs, while our conclusions for S^+2 could potentially impact the interpretation of Spitzer and SOFIA observations.

The spatial distribution of the cosmic-ray flux is important in understanding the Interstellar Medium (ISM) of the Galaxy. This distribution can be analyzed by studying different molecular species along different sight lines whose abundances are sensitive to the cosmic-ray ionization rate. Recently several groups have reported an enhanced cosmic-ray ionization rate in diffuse clouds compared to the standard value, zeta(average)=2.5e-17 s^-1, measured toward dense molecular clouds. In an earlier work we reported an enhancement in the cosmic ray rate of 20 towards HD185418. McCall et al. have reported enhancements of 48 towards zeta Persei based on the observed abundance of H$_{3}^+ while Le Petit et al. found a cosmic ray enhancement of ~10 to be consistent with their models for this same sight line. Here we revisit zeta Persei and perform a detailed calculation using a self-consistent treatment of the hydrogen chemistry, grain physics, energy and ionization balance, and excitation physics. We show that the value of zeta deduced from the H$_{3}^+$ column density in the diffuse region of the sightline depends strongly on the properties of the grains because they remove free electrons and change the hydrogen chemistry. The observations are largely consistent with a cosmic ray enhancement of 40, with several diagnostics indicating higher values. This underscores the importance of a full treatment of grain physics in studies of interstellar chemistry.

A recent modeling study of brightness ratios for CO rotational transitions in gas typical of the diffuse ISM by Liszt found the role of H collisions to be more important than previously assumed. This conclusion was based on quantum scattering calculations using the so-called WKS potential energy surface (PES) which reported a large cross section for the important 0->1 rotational transition. New close-coupling (CC) rigid-rotor calculations for CO(v=0,J=0) excitation by H are performed on four different PESs. Two of the PESs are obtained in this work using state-of-the-art quantum chemistry techniques at the CCSD(T) and MRCI levels of theory. Cross sections for the J=0->1, as well as other odd Delta J, transitions are significantly suppressed compared to even Delta J transitions in thermal energy CC calculations using the CCSD(T) and MRCI surfaces. This is consistent with a expected even Delta J propensity and in contrast to CC calculations using the earlier WKS PES which predict a dominating 0->1 transition. The current results suggest that the original astrophysical assumption that excitation of CO by H_2 dominates the kinetics of CO in diffuse ISM gas is likely to remain valid.

We present quantum mechanical close-coupling calculations of collisions between two hydrogen molecules over a wide range of energies, extending from the ultracold limit to the super-thermal region. The two most recently published potential energy surfaces for the H$_2$-H$_2$ complex, the so-called DJ (Diep and Johnson, 2000) and BMKP (Boothroyd et al., 2002) surfaces, are quantitatively evaluated and compared through the investigation of rotational transitions in H$_2$+H$_2$ collisions within rigid rotor approximation. The BMKP surface is expected to be an improvement, approaching chemical accuracy, over all conformations of the potential energy surface compared to previous calculations of H$_2$-H$_2$ interaction. We found significant differences in rotational excitation/de-excitation cross sections computed on the two surfaces in collisions between two para-H$_2$ molecules. The discrepancy persists over a large range of energies from the ultracold regime to thermal energies and occurs for several low-lying initial rotational levels. Good agreement is found with experiment (Mat\'e et al., 2005) for the lowest rotational excitation process, but only with the use of the DJ potential. Rate coefficients computed with the BMKP potential are an order of magnitude smaller.

Much of the baryonic matter in the Universe is in the form of H2 which includes most of the gas in galactic and extragalactic interstellar clouds. Molecular hydrogen plays a significant role in establishing the thermal balance in many astrophysical environments and can be important as a spectral diagnostic of the gas. Modeling and interpretation of observations of such environments requires a quantitatively complete and accurate treatment of H2. Using this micro-physical model of H2, illustrative calculations of prototypical astrophysical environments are presented. This work forms the foundation for future investigations of these and other environments where H2 is an important constituent.

Two peculiar stars showing an apparent extremely broadened and strong NaI D absorption have been discovered in surveys for cool white dwarfs by Oppenheimer et al. (2001) and Harris et al. (SDSS, 2003). We discuss the nature of these objects using PHOENIX atmosphere models for metal-poor brown dwarfs/very low mass stars, and new white dwarf LTE and NLTE models for hydrogen- and helium-dominated atmospheres with metals. These include complete molecular formation in chemical equilibrium and a model for the alkali resonance line broadening based on the damping profiles of Allard et al. (2003), as well as new molecular line opacities for metal hydrides. First results of our calculations indicate good agreement with a hydrogen-dominated WD atmosphere with a Na abundance roughly consistent with a state of high accretion. We analyse deviations of the abundances of Na, K, Mg and Ca from the cosmic pattern and comment on implications of these results for standard accretion scenarios.

The two most recently published potential energy surfaces (PESs) for the HeH$_2$ complex, the so-called MR (Muchnick and Russek) and BMP (Boothroyd, Martin, and Peterson) surfaces, are quantitatively evaluated and compared through the investigation of atom-diatom collision processes. The BMP surface is expected to be an improvement, approaching chemical accuracy, over all conformations of the PES compared to that of the MR surface. We found significant differences in inelastic rovibrational cross sections computed on the two surfaces for processes dominated by large changes in target rotational angular momentum. In particular, the H$_2$($\nu=1, j=0$) total quenching cross section computed using the BMP potential was found to be a factor of $~$1000 larger than that obtained with the MR surface.

We present a complete line list for the X 1Sigma+ electronic ground state of LiCl computed using fully quantum-mechanical techniques. This list includes transition energies and oscillator strengths in the spectral region 0.3-39,640.7 cm-1 for all allowed rovibrational transitions in absorption within the electronic ground state. The calculations were performed using an accurate hybrid potential constructed from a spectral inversion fit of experimental data and from recent multi-reference single- and double-excitation configuration interaction calculations. The line list was incorporated into the stellar atmosphere code PHOENIX to compute spectra for a range of young to old T dwarf models. The possibility of observing a signature of LiCl in absorption near 15.8 microns is addressed and the proposal to use this feature to estimate the total lithium elemental abundance for these cool objects is discussed.

Krstic has carried out the first quantum mechanical calculations near threshold for the charge transfer (CT) process H^+ + H_2(X ^1Sigma_g^+, nu=0, J=0) --> H(1s) + H_2^+. These results are relevant for models of primordial galaxy and first star formation that require reliable atomic and molecular data for obtaining the early universe hydrogen chemistry. Using the results of Krstic, we calculate the relevant CT rate coefficient for temperatures between 100 and 30,000 K. We also present a simple fit which can be readily implemented into early universe chemical models. Additionally, we explore how the range of previously published data for this reaction translates into uncertainties in the predicted gas temperature and H_2 relative abundance in a collapsing primordial gas cloud. Our new data significantly reduce these cosmological uncertainties that are due to the uncertainties in the previously published CT rate coefficients.

The cosmological recombination history of lithium, produced during Big--Bang nucleosynthesis, is presented using updated chemistry and cosmological parameters consistent with recent cosmic microwave background (CMB) measurements. For the popular set of cosmological parameters, about a fifth of the lithium ions recombine into neutral atoms by a redshift $z\sim 400$. The neutral lithium atoms scatter resonantly the CMB at 6708 \AA and distort its intensity and polarization anisotropies at observed wavelengths around $\sim 300 \mu$m, as originally suggested by Loeb (2001). The modified anistropies resulting from the lithium recombination history are calculated for a variety of cosmological models and found to result primarily in a suppression of the power spectrum amplitude. Significant modification of the power spectrum occurs for models which assume a large primordial abundance of lithium. While detection of the lithium signal might prove difficult, if offers the possibility of inferring the lithium primordial abundance and is the only probe proposed to date of the large-scale structure of the Universe for $z\sim 500-100$.

The opacity due to photodissociation of 24MgH is investigated in the atmospheres of cool stars. The lowest two electronic transitions A 2Pi -- X 2Sigma+ and B' 2Sigma+ -- X 2 Sigma+ are considered where the cross sections for the latter were published previously (Weck, Stancil, & Kirby 2002) while the former are presented in this work. Model atmospheres calculated with the PHOENIX code are used to investigate the effect of the photodissociation opacity on spectra of cool stars. The A -- X photodissociation cross sections are obtained using a combination of ab initio and experimentally derived potential curves and dipole transition moments. Partial cross sections have been evaluated over the accessible wavelength range 1770-4560 Angstrom for all rotational transitions from the vibrational levels v''=0-11. Assuming a Boltzmann distribution of the rovibrational levels of the X 2Sigma+ state, LTE photodissociation cross sections are presented for temperatures between 1000 and 5000 K. Shape resonances, arising from rotational predissociation of quasi-bound levels of the A 2Pi state near threshold, characterize the LTE photodissociation cross sections. A sum rule is proposed as a check on the accuracy of the photodissociation calculations.

A new, complete, theoretical rotational and vibrational line list for the A-X electronic transition in MgH is presented. The list includes transition energies and oscillator strengths for all possible allowed transitions and was computed using the best available theoretical potential energies and dipole transition moment function with the former adjusted to account for experimental data. The A-X line list, as well as new line lists for the B'-X and the X-X (pure rovibrational) transitions, were included in comprehensive stellar atmosphere models for M, L, and T dwarfs and solar-type stars. The resulting spectra, when compared to models lacking MgH, show that MgH provides significant opacity in the visible between 4400 and 5600 Angstrom. Further, comparison of the spectra obtained with the current line list to spectra obtained using the line list constructed by Kurucz (1993) show that the Kurucz list significantly overestimates the opacity due to MgH particularly for the bands near 5150 and 4800 Angstrom with the discrepancy increasing with decreasing effective temperature.

The recent announcement that sodium absorption has been observed in the atmosphere of HD209458b, the only EGP observed to transit its parent star, is the first direct detection of an EGP atmosphere. We explore the possibility that neutral sodium is {\em not} in local thermodynamic equilibrium (LTE) in the outer atmosphere of irradiated EGPs and that the sodium concentration may be underestimated by models that make the LTE assumption. Our results indicate that it may not be necessary to invoke excessive photoionization, low metallicity, or even high altitude clouds to explain the observations.