Showing 1–50 of 53 results for author: Stancil, P C

PDRs4All IX. Sulfur elemental abundance in the Orion Bar

Authors: Asunción Fuente, Evelyne Roueff, Franck Le Petit, Jacques Le Bourlot, Emeric Bron, Mark G. Wolfire, James F. Babb, Pei-Gen Yan, Takashi Onaka, John H. Black, Ilane Schroetter, Dries Van De Putte, Ameek Sidhu, Amélie Canin, Boris Trahin, Felipe Alarcón, Ryan Chown, Olga Kannavou, Olivier Berné, Emilie Habart, Els Peeters, Javier R. Goicoechea, Marion Zannese, Raphael Meshaka, Yoko Okada , et al. (9 additional authors not shown)

Abstract: One of the main problems in astrochemistry is determining the amount of sulfur in volatiles and refractories in the interstellar medium. The detection of the main sulfur reservoirs (icy H$_2$S and atomic gas) has been challenging, and estimates are based on the reliability of models to account for the abundances of species containing less than 1% of the total sulfur. The high sensitivity of the Ja… ▽ More One of the main problems in astrochemistry is determining the amount of sulfur in volatiles and refractories in the interstellar medium. The detection of the main sulfur reservoirs (icy H$_2$S and atomic gas) has been challenging, and estimates are based on the reliability of models to account for the abundances of species containing less than 1% of the total sulfur. The high sensitivity of the James Webb Space Telescope provides an unprecedented opportunity to estimate the sulfur abundance through the observation of the [S I] 25.249 $μ$m line. We used the [S III] 18.7 $μ$m, [S IV] 10.5 $μ$m, and [S l] 25.249 $μ$m lines to estimate the amount of sulfur in the ionized and molecular gas along the Orion Bar. For the theoretical part, we used an upgraded version of the Meudon photodissociation region (PDR) code to model the observations. New inelastic collision rates of neutral atomic sulfur with ortho- and para- molecular hydrogen were calculated to predict the line intensities. The [S III] 18.7 $μ$m and [S IV] 10.5 $μ$m lines are detected over the imaged region with a shallow increase (by a factor of 4) toward the HII region. We estimate a moderate sulfur depletion, by a factor of $\sim$2, in the ionized gas. The corrugated interface between the molecular and atomic phases gives rise to several edge-on dissociation fronts we refer to as DF1, DF2, and DF3. The [S l] 25.249 $μ$m line is only detected toward DF2 and DF3, the dissociation fronts located farthest from the HII region. The detailed modeling of DF3 using the Meudon PDR code shows that the emission of the [S l] 25.249 $μ$m line is coming from warm ($>$ 40 K) molecular gas located at A$_{\rm V}$ $\sim$ 1$-$5 mag from the ionization front. Moreover, the intensity of the [S l] 25.249 $μ$m line is only accounted for if we assume the presence of undepleted sulfur. △ Less

Submitted 4 June, 2024; v1 submitted 14 April, 2024; originally announced April 2024.

Comments: 16 pages, 6 figures. Accepted for publication in Astronomy and Astrophysics

A New Definition of Exoplanet Habitability: Introducing the Photosynthetic Habitable Zone

Authors: C. Hall, P. C. Stancil, J. P. Terry, C. K. Ellison

Abstract: It may be possible to detect biosignatures of photosynthesis in an exoplanet's atmosphere. However, such a detection would likely require a dedicated study, occupying a large amount of telescope time. It is therefore prudent, while searching for signs of life that we may recognise, to pick the best target possible. In this work, we present a new region, the ``photosynthetic habitable zone'' \texte… ▽ More It may be possible to detect biosignatures of photosynthesis in an exoplanet's atmosphere. However, such a detection would likely require a dedicated study, occupying a large amount of telescope time. It is therefore prudent, while searching for signs of life that we may recognise, to pick the best target possible. In this work, we present a new region, the ``photosynthetic habitable zone'' \textemdash the distance from a star where both liquid water and oxygenic photosynthesis can occur. It is therefore the region where detectable biosignatures of oxygenic photosynthesis are most likely to occur. Our analysis indicates that in the most ideal conditions for life and no atmospheric effects, the photosynthetic habitable zone is almost as broad as the habitable zone. On the other hand, if conditions for life are anything less than excellent and atmospheric effects are even moderate, the photosynthetic habitable zone is concentrated at larger separations around more massive stars. Such cases are also not tidally locked to their host star, which could result in planetary rotation periods similar to the Earth's. We identify five planets, Kepler-452 b, Kepler-1638 b, Kepler-1544 b and Kepler-62 e and Kepler-62 f, that are consistently in the photosynthetic habitable zone for a variety of conditions, and we predict their day lengths to be between 9 and 11 hours. We conclude that the parameter space in which we should search for signs of life is much narrower than the standard habitable zone. △ Less

Submitted 12 April, 2023; v1 submitted 31 January, 2023; originally announced January 2023.

Comments: 12 pages, 3 figures, accepted to ApJL

Relativistic Atomic Structure of Au IV and the Os Isoelectronic Sequence: opacity data for kilonova ejecta

Authors: Z. S. Taghadomi, Y. Wan, A. Flowers, P. C. Stancil, B. M. McLaughlin, S. Bromley, J. P. Marler, C. E. Sosolik, S. Loch

Abstract: Direct detection of gravitational waves (GW) on Aug. 17, 2017, propagating from a binary neutron star merger, opened the era of multimessenger astronomy. The ejected material from neutron star mergers, or kilonova, is a good candidate for optical and near infrared follow-up observations after the detection of GWs. The kilonova from the ejecta of GW1780817 provided the first evidence for the astrop… ▽ More Direct detection of gravitational waves (GW) on Aug. 17, 2017, propagating from a binary neutron star merger, opened the era of multimessenger astronomy. The ejected material from neutron star mergers, or kilonova, is a good candidate for optical and near infrared follow-up observations after the detection of GWs. The kilonova from the ejecta of GW1780817 provided the first evidence for the astrophysical site of the synthesis of heavy nuclei through the rapid neutron capture process or rprocess. Since properties of the emission are largely affected by opacities of the ejected material, enhancements in the available, but limited rprocess atomic data have been motivated recently. However, given the complexity of the electronic structure of these heavy elements, considerable efforts are still needed to converge to a reliable set of atomic structure data. The aim of this work is to alleviate this situation for low charge state elements in the Os-like isoelectronic sequence. In this regard, the general purpose relativistic atomic structure packages (GRASP0 and GRASP2K) were used to obtain energy levels and transition probabilities (E1 and M1). We provide line lists and expansion opacities for a range of r-process elements. We focus here on the Os isoelectronic sequence (Os I, Ir II, Pt III, Au IV, Hg V). The results are benchmarked against existing experimental data and prior calculations, and predictions of emission spectra relevant to kilonovae are provided. Finestructure (M1) lines in the infrared potentially observable by the James Webb Space Telescope are highlighted. △ Less

Submitted 5 July, 2022; originally announced July 2022.

Revisiting the Formation of HeH$^+$ in the Planetary Nebula NGC 7027

Authors: R. C. Forrey, J. F. Babb, E. D. S. Courtney, R. McArdle, P. C. Stancil

Abstract: 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 li… ▽ More 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. △ Less

Submitted 15 June, 2020; originally announced June 2020.

Comments: etable of rate coefficients available

2018 NASA Laboratory Astrophysics Workshop: Scientific Organizing Committee Report

Authors: N. Brickhouse, G. J. Ferland, S. Milam, E. Sciamma-O'Brien, A. Smale, A. Spyrou, P. C. Stancil, L. Storrie-Lombardi, G. M. Wahlgren

Abstract: 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 Worksh… ▽ More 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. △ Less

Submitted 28 January, 2020; originally announced January 2020.

Perspectives on Astrophysics Based on Atomic, Molecular, and Optical (AMO) Techniques

Authors: Daniel Wolf Savin, James F. Babb, Paul M. Bellan, Crystal Brogan, Jan Cami, Paola Caselli, Lia Corrales, Gerardo Dominguez, Steven R. Federman, Chris J. Fontes, Richard Freedman, Brad Gibson, Leon Golub, Thomas W. Gorczyca, Michael Hahn, Sarah M. Hörst, Reggie L. Hudson, Jeffrey Kuhn, James E. Lawler, Maurice A. Leutenegger, Joan P. Marler, Michael C. McCarthy, Brett A. McGuire, Stefanie N. Milam, Nicholas A. Murphy , et al. (13 additional authors not shown)

Abstract: 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 i… ▽ More 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. △ Less

Submitted 14 November, 2018; originally announced November 2018.

Comments: White paper submission to the Decadal Assessment and Outlook Report on Atomic, Molecular, and Optical (AMO) Science (AMO 2020)

Photodissociation of CS from Excited Rovibrational Levels

Authors: Ryan J. Pattillo, René Cieszewski, Phillip C. Stancil, Robert C. Forrey, James F. Babb, James F. McCann, Brendan M. McLaughlin

Abstract: Accurate photodissociation cross sections have been computed for transitions from the X $^1Σ^+$ 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 bas… ▽ More Accurate photodissociation cross sections have been computed for transitions from the X $^1Σ^+$ 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Σ^+$ state and for photon wavelengths ranging from 500 $\textÅ$ 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. △ Less

Submitted 29 April, 2018; v1 submitted 3 April, 2018; originally announced April 2018.

Comments: 16 pages, 9 figures, 4 tables

Journal ref: ApJ 858, 10 (2018)

Full-dimensional Quantum Dynamics of SiO in Collision with H$_2$

Authors: Benhui Yang, P. Zhang, Chen Qu, X. H. Wang, P. C. Stancil, J. M. Bowman, N. Balakrishnan, B. M. McLaughlin, R. C. Forrey

Abstract: 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… ▽ More 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. △ Less

Submitted 24 January, 2018; originally announced February 2018.

Editorial: Special Issue on the Atomic and Molecular Processes in the Ultracold Regime, the Chemical Regime, and Astrophysics

Authors: J. F. Babb, R. Côté, H. R. Sadeghpour, P. C. Stancil

Abstract: 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. 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. △ Less

Submitted 2 August, 2017; originally announced August 2017.

Comments: 7 pp

Journal ref: J. Phys. B: Atomic, Molecular and Optical Physics 50 (2017) 140202

Line Ratios for Solar Wind Charge Exchange with Comets

Authors: P. D. Mullen, R. S. Cumbee, D. Lyons, L. Gu, J. Kaastra, R. L. Shelton, P. C. Stancil

Abstract: 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 imposs… ▽ More 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. △ Less

Submitted 27 July, 2017; originally announced July 2017.

Journal ref: ApJ, 844, 1 (2017)

Rate constants for the formation of SiO by radiative association

Authors: M. Cairnie, R. C. Forrey, J. F. Babb, P. C. Stancil, B. M. McLaughlin

Abstract: 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 const… ▽ More 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. △ Less

Submitted 6 July, 2017; originally announced July 2017.

Comments: To appear in MNRAS

The 2017 Release of Cloudy

Authors: G. J. Ferland, M. Chatzikos, F. Guzmán, M. L. Lykins, P. A. M. van Hoof, R. J. R. Williams, N. P. Abel, N. R. Badnell, F. P. Keenan, R. L. Porter, P. C. Stancil

Abstract: 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… ▽ More 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. △ Less

Submitted 5 October, 2023; v1 submitted 30 May, 2017; originally announced May 2017.

Comments: replace Fig 6 with the published version

Journal ref: 2017RMxAA..53..385F

Inelastic cross sections and rate coefficients for collisions between CO and H2

Authors: Christina Castro, Kyle Doan, Michael Klemka, Robert C. Forrey, B. H. Yang, Phillip C. Stancil, N. Balakrishnan

Abstract: 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… ▽ More 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. △ Less

Submitted 18 January, 2017; originally announced January 2017.

PAMOP project: computations in support of experiments and astrophysical applications

Authors: B M McLaughlin, C P Ballance, M S Pindzola, P C Stancil, S Schippers, A Müller

Abstract: 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ödinger or Dirac equation for the appropriate collision problem using the R-matrix or R-matrix with pseudo-states approach from f… ▽ More 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ödinger 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. △ Less

Submitted 14 January, 2017; originally announced January 2017.

Comments: 17 pages, 10 figures: chapter in the book, High Performance Computing in Science and Engineering'16, edited by W. E. Nagel, D. B. Kröner, and M. Reich (Springer, New York and Berlin, 2017)

A Comprehensive Study of Hydrogen Adsorbing to Amorphous Water-Ice: Defining Adsorption in Classical Molecular Dynamics

Authors: John L. Dupuy, Steven P. Lewis, P. C. Stancil

Abstract: 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 d… ▽ More 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. △ Less

Submitted 17 August, 2016; originally announced August 2016.

Comments: 30 pages, 10 figures, 2 tables

Journal ref: Dupuy, J. L., Lewis, S. P., Stancil, P. C. 2016, ApJ, 831, 54

Ne X X-ray Emission due to Charge Exchange in M82

Authors: R. S. Cumbee, L. Liu, D. Lyons, D. R. Schultz, P. C. Stancil, J. G. Wang, R. Ali

Abstract: 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 ran… ▽ More 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. △ Less

Submitted 15 March, 2016; originally announced March 2016.

Comments: 8 pages, 9 figures

Rovibrationally-Resolved Photodissociaton of SH$^+$

Authors: E. C. McMillan, G. Shen, J. F. McCann, B. M. McLaughlin, P. C. Stancil

Abstract: Photodissociation cross sections for the SH$^+$ radical are computed from all rovibrational (RV) levels of the ground electronic state X$~^3Σ^-$ for wavelengths from threshold to 500~Å. The five electronic transitions, $2~ ^3Σ^- \leftarrow$ X$~^3Σ^-$, $3~ ^3Σ^- \leftarrow$ X$~^3Σ^-$, $A~ ^3Π\leftarrow$ X$~^3Σ^-$, $2~ ^3Π\leftarrow$ X$~^3Σ^-$, and $3~ ^3Π\leftarrow$ X$~^3Σ^-$, are treated with a fu… ▽ More Photodissociation cross sections for the SH$^+$ radical are computed from all rovibrational (RV) levels of the ground electronic state X$~^3Σ^-$ for wavelengths from threshold to 500~Å. The five electronic transitions, $2~ ^3Σ^- \leftarrow$ X$~^3Σ^-$, $3~ ^3Σ^- \leftarrow$ X$~^3Σ^-$, $A~ ^3Π\leftarrow$ X$~^3Σ^-$, $2~ ^3Π\leftarrow$ X$~^3Σ^-$, and $3~ ^3Π\leftarrow$ X$~^3Σ^-$, 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Σ^-$ 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. △ Less

Submitted 17 February, 2016; originally announced February 2016.

Comments: 21 pages, 13 figures, 2 tables. Accepted in J Phys. B: At. Mol. Opt. Phys. for the special issue on Atomic and molecular data for astrophysics

Charge Exchange Induced X-ray Emission of Fe XXV and Fe XXVI via a Streamlined Model

Authors: P. D. Mullen, R. S. Cumbee, D. Lyons, P. C. Stancil

Abstract: 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,… ▽ More 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. △ Less

Submitted 26 February, 2017; v1 submitted 7 February, 2016; originally announced February 2016.

Comments: -Keywords: Charge Exchange; Landau-Zener; Galactic center/Galactic ridge

Journal ref: 2016ApJS..224...31M

Quantum Calculation of Inelastic CO Collisions with H. III. Rate Coefficients for Ro-vibrational Transitions

Authors: L. Song, N. Balakrishnan, K. M. Walker, P. C. Stancil, W. F. Thi, I. Kamp, A. van der Avoird, G. C. Groenenboom

Abstract: 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 com… ▽ More 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. △ Less

Submitted 19 November, 2015; originally announced November 2015.

A novel scenario for the possible X-ray line feature at ~3.5 keV: Charge exchange with bare sulfur ions

Authors: L. Gu, J. Kaastra, A. J. J. Raassen, P. D. Mullen, R. S. Cumbee, D. Lyons, P. C. Stancil

Abstract: 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… ▽ More 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. △ Less

Submitted 20 November, 2015; originally announced November 2015.

Comments: 5 pages, 3 figures; Astronomy & Astrophysics, in press

Quantum Calculation of Inelastic CO Collisions with H. II. Pure Rotational Quenching of High Rotational Levels

Authors: Kyle M. Walker, L. Song, B. H. Yang, G. C. Groenenboom, A. van der Avoird, N. Balakrishnan, R. C. Forrey, P. C. Stancil

Abstract: 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 po… ▽ More 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. △ Less

Submitted 16 November, 2015; originally announced November 2015.

Comments: 8 figures, 1 table

Journal ref: 2015, ApJ, 811, 27

Stout: Cloudy's Atomic and Molecular Database

Authors: M. L. Lykins, G. J. Ferland, R. Kisielius, M. Chatzikos, R. L. Porter, P. A. M. van Hoof, R. J. R. Williams, F. P. Keenan, P. C. Stancil

Abstract: 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… ▽ More 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. △ Less

Submitted 4 June, 2015; originally announced June 2015.

Comments: ApJ in press

Collisional quenching of highly rotationally excited HF

Authors: Benhui Yang, K. M. Walker, R. C. Forrey, P. C. Stancil, N. Balakrishnan

Abstract: 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 r… ▽ More 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. △ Less

Submitted 28 April, 2015; originally announced April 2015.

Comments: 26 pages, 14 figures, and 3 tables in A&A 2015

Journal ref: A&A 578, A65 (2015)

Radiative charge transfer in cold and ultracold Sulfur atoms colliding with Protons

Authors: G Shen, P C Stancil, J G Wang, J F McCann, B M McLaughlin

Abstract: 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) approxim… ▽ More 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 $μ$ 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. △ Less

Submitted 25 February, 2015; originally announced February 2015.

Comments: 16 pages, 9 figures, accepted for publication in J Phys B: At. Mol. Opt. Phys

Quantum dynamics of CO-H$_2$ in full dimensionality

Authors: Benhui Yang, P. Zhang, X. Wang, P. C. Stancil, J. M. Bowman, N. Balakrishnan, R. C. Forrey

Abstract: 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… ▽ More 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. △ Less

Submitted 9 February, 2015; originally announced February 2015.

XMM-Newton Measurement of the Galactic Halo X-ray Emission using a Compact Shadowing Cloud

Authors: David B. Henley, Robin L. Shelton, Renata S. Cumbee, Phillip C. Stancil

Abstract: 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… ▽ More 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. △ Less

Submitted 18 November, 2014; originally announced November 2014.

Comments: 10 pages, 6 figures. Accepted for publication in the Astrophysical Journal

On the Validity of Collider-Mass Scaling for Molecular Rotational Excitation

Authors: Kyle M. Walker, B. H. Yang, P. C. Stancil, N. Balakrishnan, R. C. Forrey

Abstract: 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 r… ▽ More 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. △ Less

Submitted 22 June, 2014; originally announced June 2014.

Comments: 23 pages, 8 figures. Accepted for publication in ApJ

Can Charge Exchange Explain Anomalous Soft X-ray Emission in the Cygnus Loop?

Authors: Renata S. Cumbee, David B. Henley, Phillip C. Stancil, Robin L. Shelton, Jeff L. Nolte, Yong Wu, David R. Schultz

Abstract: 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+}$,… ▽ More 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. △ Less

Submitted 5 May, 2014; originally announced May 2014.

Comments: Accepted for publication in ApJL, 5 pages, 4 figures, 1 table

Quenching of Highly-Rotationally-Excited HCl in Collisions with HE

Authors: Benhui Yang, Phillip C. Stancil

Abstract: 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… ▽ More 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. △ Less

Submitted 22 January, 2014; originally announced January 2014.

Quantum Calculation of Inelastic CO Collisions with H: I. rotational quenching of low-lying rotational levels

Authors: Benhui Yang, P. C. Stancil, N. Balakrishnan, R. C. Forrey, J. M. Bowman

Abstract: 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 numbe… ▽ More 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. △ Less

Submitted 10 May, 2013; originally announced May 2013.

Comments: 21 pages, 5 figures

The 2013 Release of Cloudy

Authors: G. J. Ferland, R. L. Porter, P. A. M. van Hoof, R. J. R. Williams, N. P. Abel, M. L. Lykins, Gargi Shaw, W. J. Henney, P. C. Stancil

Abstract: 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.… ▽ More 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. △ Less

Submitted 18 February, 2013; originally announced February 2013.

Comments: Major review article to be published in RevMexAA. 27 pages, 18 figures

Rotational quenching of rotationally-excited H$_2$O in collisions with He

Authors: Benhui Yang, M. Nagao, W. Satomi, M. Kimura, P. C. Stancil

Abstract: 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 fo… ▽ More 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. △ Less

Submitted 15 February, 2013; originally announced February 2013.

Comments: 29 pages, 9 figures, 1 table

Advances in Atomic Data for Neutron-Capture Elements

Authors: N. C. Sterling, M. C. Witthoeft, D. A. Esteves, P. C. Stancil, A. L. D. Kilcoyne, R. C. Bilodeau, A. Aguilar

Abstract: 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… ▽ More 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. △ Less

Submitted 9 September, 2011; originally announced September 2011.

Comments: To be published in IAU Symp. 283: Planetary Nebulae, an Eye to the Future; 2 pages

Atomic data for neutron-capture elements III. Charge transfer rate coefficients for low-charge ions of Ge, Se, Br, Kr, Rb, and Xe

Authors: N. C. Sterling, P. C. Stancil

Abstract: 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… ▽ More 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. △ Less

Submitted 18 August, 2011; originally announced August 2011.

Comments: 11 pages, 4 figures, accepted for publication in Astronomy & Astrophysics

Rovibrationally resolved photodissociation of HeH+

Authors: S. Miyake, C. D. Gay, P. C. Stancil

Abstract: 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 l… ▽ More 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. △ Less

Submitted 29 April, 2011; originally announced April 2011.

Critical Test of Simulations of Charge-Exchange-Induced X-Ray Emission in the Solar System

Authors: R. Ali, P. A. Neill, P. Beiersdorfer, C. L. Harris, D. R. Schultz, P. C. Stancil

Abstract: 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… ▽ More 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. △ Less

Submitted 14 May, 2010; originally announced May 2010.

Comments: 4 figures

Rotational quenching of CO due to H$_2$ collisions

Authors: Benhui Yang, P. C. Stancil, N. Balakrishnan, R. C. Forrey

Abstract: 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… ▽ More 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 $|Δ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. △ Less

Submitted 1 June, 2010; v1 submitted 22 April, 2010; originally announced April 2010.

Comments: 28 pages, 10 figures

The ion-induced charge-exchange X-ray emission of the Jovian Auroras: Magnetospheric or solar wind origin?

Authors: Yawei Hui, David R. Schultz, Vasili A. Kharchenko, Phillip C. Stancil, Thomas E. Cravens, Carey M. Lisse, Alexander Dalgarno

Abstract: 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… ▽ More 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. △ Less

Submitted 9 July, 2009; originally announced July 2009.

Comments: 11 pages, 2 figures, 2 tables, submitted to ApJ Letter

Journal ref: Astrophys.J.702:L158-L162,2009

Rotationally Warm Molecular Hydrogen in the Orion Bar

Authors: Gargi Shaw, G. J. Ferland, W. J. Henney, P. C. Stancil, N. P. Abel, E. W. Pellegrini, J. A. Baldwin, P. A. M. van Hoof

Abstract: 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 obser… ▽ More 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. △ Less

Submitted 12 June, 2009; originally announced June 2009.

Comments: Accepted for publication in ApJ (34 pages, including 16 figures)

Journal ref: Astrophys.J.701:677-685,2009

Rotational quenching rate coefficients for H_2 in collisions with H_2 from 2 to 10,000 K

Authors: T. -G. Lee, N. Balakrishnan, R. C. Forrey, P. C. Stancil, G. Shaw, D. R. Schultz, G. J. Ferland

Abstract: 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… ▽ More 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. △ Less

Submitted 12 May, 2008; originally announced May 2008.

The Effects of Doubly Ionized Chemistry on SH+ and S^+2 Abundances in X-ray Dominated Regions

Authors: N. P. Abel, S. R. Federman, P. C. Stancil

Abstract: 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 th… ▽ More 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. △ Less

Submitted 24 January, 2008; originally announced January 2008.

Comments: 19 pages, 3 figures, Accepted for Publication in ApJ Letters

On the enhanced cosmic-ray ionization rate in the diffuse cloud towards Zeta Persei

Authors: Gargi Shaw, G. J. Ferland, R. Srianand, N. P. Abel, P. A. M. van Hoof, P. C. Stancil

Abstract: 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… ▽ More 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. △ Less

Submitted 2 November, 2007; originally announced November 2007.

Comments: 25 pages, 4 figures, accepted for publication in ApJ

Low energy H+CO scattering revisited: CO rotational excitation with new potential surfaces

Authors: B. C. Shepler, B. H. Yang, T. J. Dhilip Kumar, P. C. Stancil, J. M. Bowman, N. Balakrishnan, P. Zhang, E. Bodo, A. Dalgarno

Abstract: 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… ▽ More 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. △ Less

Submitted 17 September, 2007; originally announced September 2007.

Comments: 3 figures

State-to-state rotational transitions in H$_2$+H$_2$ collisions at low temperatures

Authors: Teck-Ghee Lee, N. Balakrishnan, R. C. Forrey, P. C. Stancil, D. R. Schultz, Gary J. Ferland

Abstract: 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 evaluate… ▽ More 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é 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. △ Less

Submitted 18 July, 2006; originally announced July 2006.

Comments: Accepted by J. Chem. Phys

Molecular Hydrogen in Star-forming regions: implementation of its micro-physics in Cloudy

Authors: G. Shaw, G. J. Ferland, N. P. Abel, P. C. Stancil, P. A. M. van Hoof

Abstract: 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 qu… ▽ More 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. △ Less

Submitted 22 January, 2005; originally announced January 2005.

Comments: Accepted for publication in ApJ (58 pages). Version with full resolution is available at http://www.pa.uky.edu/~gshaw/paper/manuscript.pdf

Journal ref: Astrophys.J. 624 (2005) 794-807

Modelling Alkali Line Absorption and Molecular Bands in Cool DAZs

Authors: Derek Homeier, Nicole F. Allard, France Allard, Peter H. Hauschildt, Andreas Schweitzer, Phillip C. Stancil, Philippe F. Weck

Abstract: 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 a… ▽ More 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. △ Less

Submitted 5 January, 2005; originally announced January 2005.

Comments: to appear in: Proceedings of the 14th European Workshop on White Dwarfs, http://www.astrophysik.uni-kiel.de/~eurowd/ -- ASP Conference Series, uses epsf.sty

Close-coupling calculations of low-energy inelastic and elastic processes in $^4$He collisions with H$_2$: A comparative study of two potential energy surfaces

Authors: T. G. Lee, C. Rochow, R. Martin, T. K. Clark, R. C. Forrey, N. Balakrishnan, P. C. Stancil, D. R. Schultz, A. Dalgarno, G. J. Ferland

Abstract: 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 PE… ▽ More 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$($ν=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. △ Less

Submitted 9 December, 2004; originally announced December 2004.

Comments: 26 pages, 11 figures

Molecular line opacity of LiCl in the mid-infrared spectra of brown dwarfs

Authors: P. F. Weck, A. Schweitzer, K. Kirby, P. H. Hauschildt, P. C. Stancil

Abstract: 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… ▽ More 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. △ Less

Submitted 6 June, 2004; originally announced June 2004.

Comments: 8 pages, 2 figures, 1 table. Accepted for publication in ApJ 613, Sept. 20 2004

Journal ref: Astrophys.J. 613 (2004) 567-571

Rate Coefficient for H^+ + H_2(X ^1Sigma_g^+, nu=0, J=0) --> H(1s) + H_2^+ Charge Transfer and Some Cosmological Implications

Authors: Daniel Wolf Savin, Predrag S. Krstic, Zoltan Haiman, Phillip C. Stancil

Abstract: 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… ▽ More 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. △ Less

Submitted 14 April, 2004; originally announced April 2004.

Comments: This version corrects a typo in Eq. 4 of our ApJL. 11 pages, 2 figures

The Molecular Continuum Opacity of MgH in Cool Stellar Atmospheres

Authors: P. F. Weck, A. Schweitzer, P. C. Stancil, P. H. Hauschildt, K. Kirby

Abstract: 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… ▽ More 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. △ Less

Submitted 1 July, 2002; originally announced July 2002.

Comments: 9 pages, 8 figures

Journal ref: Astrophys.J. 584 (2003) 459-464