Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ab initio data”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,171 records · Page 65
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Atlas of putative minima and low-lying energy networks of water clusters n = 3 - 25

We report a database consisting of the putative minima and ~3.2 million local minima lying within 5 kcal/mol from the putative minima for water clusters of sizes n = 3 – 25 using an improved version of the Monte Carlo Temperature Basin Paving (MCTBP) global optimization procedure in conjunction with the ab-initio based, flexible, polarizable Thole-Type Model (TTM2.1-F, version 2.1) interaction potential for water. Several of the low-lying structures, as well as lowlying penta-coordinated water networks obtained with the TTM2.1-F potential, were further refined at the Møller-Plesset second order perturbation MP2/aug-cc-pVTZ level of theory. In total, we have identified 3,138,303 networks corresponding to local minima of the clusters n = 3 – 25, whose Cartesian coordinates and relative energies can be obtained from the webpage https://sites.uw.edu/wdbase/. Networks containing penta-coordinated water molecules start to appear at n = 11 and, quite surprisingly, are energetically close (within 1 – 3 kcal/mol) to the putative minima, a fact that has been confirmed from the MP2 calculations. This large database of water cluster minima spanning quite dissimilar hydrogen bonding networks is expected to influence the development and assessment of the accuracy of interaction potentials for water as well as lower scaling electronic structure methods (such as different density functionals). Furthermore, it can also be used in conjunction with data science approaches (including but no limited to neural networks, machine and deep learning) to understand the properties of water, nature’s most important substance.

Rakshit, Avijit↗
DOE OSTI2021

Improved first-principles equation-of-state table of deuterium for high-energy-density applications

Here, we present a first-principles equation-of-state (EOS) table of deuterium aimed at improving the previously established first-principles equation-of-state table (FPEOS). The EOS table presented here, referred to as iFPEOS, introduces (1) a universal density functional theory (DFT) treatment of all density and temperature conditions, (2) a fully consistent treatment of exchange-correlation (XC) thermal effects across the entire range of temperatures covered, and (3) quantum treatment of ions. Based on ab initio molecular dynamics driven by thermal density functional theory, iFPEOS includes density points in the range 1×10 –3 $\leqslant$ ρ $\leqslant$ 1.6×10 3 g/cm 3 and temperature points in the range 800 K $\leqslant$ T $\leqslant$ 256 MK, thus covering the challenging warm dense matter (WDM) regime. For an improved description of the electronic structure, iFPEOS employs an advanced free-energy XC density functional with explicit temperature dependence, which is at the metageneralized gradient approximation level of DFT. We use the latest orbital-free free-energy density functional for the high-temperature regime where it shows excellent agreement with standard Mermin-Kohn-Sham DFT. For quantum treatment of ions we use path-integral molecular dynamics in order to take into account nuclear quantum effects. Results are compared to other EOS models and most recent experimental measurements of deuterium properties such as the molecular-toatomic fluid transition, the principal and reshock Hugoniot, and sound speed. We find that iFPEOS provides an improved agreement with experimental data compared to other first-principles EOS models in the WDM regime for pressures up to 200 GPa and temperatures up to 60 000 K. For higher pressures and temperatures, however, iFPEOS is in agreement with other models in predicting lower compressibility and higher sound speed along the Hugoniot, compared to experiment.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Electrons to Reactors Multiscale Modeling: Catalytic CO Oxidation over RuO 2

First-principles kinetic Monte Carlo (1p-kMC) simulations for CO oxidation on two RuO 2 facets, RuO 2 (110) and RuO 2 (111), were coupled to the computational fluid dynamics (CFD) simulations package MFIX, and reactor-scale simulations were then performed. 1p-kMC coupled with CFD has recently been shown as a feasible method for translating molecular scale mechanistic knowledge to the reactor scale, enabling comparisons to in situ and online experimental measurements. Only a few studies with such coupling have been published. This work incorporates multiple catalytic surface facets into the scale-coupled simulation, and three possibilities were investigated: the two possibilities of each facet individually being the dominant phase in the reactor, and also the possibility that both facets were present on the catalyst particles in the ratio predicted by an ab initio thermodynamics-based Wulff construction. When lateral interactions between adsorbates were included in the 1p-kMC simulations, the two surfaces, RuO 2 (110) and RuO 2 (111), were found to be of similar order-of-magnitude in activity for the pressure range of 1 × 10 –4 bar to 1 bar, with the RuO 2 (110) surface-termination showing more simulated activity than the RuO 2 (111) surface-termination. Coupling between the 1p-kMC and CFD was achieved with a lookup table generated by the error-based modified Shepard interpolation scheme. Isothermal reactor scale simulations were performed and compared to two separate experimental studies, conducted with reactant partial pressures of ≤0.1 bar. Simulations without an isothermality restriction were also conducted and showed that the simulated temperature gradient across the catalytic reactor bed is <0.5 K, which validated the use of the isothermality restriction for investigating the reactor-scale phenomenological temperature dependences. The approach with the Wulff construction based reactor simulations reproduced a trend similar to one experimental data set relatively well, with the (110) surface being more active at higher temperaures; in contrast, for the other experimental data set, our reactor simulations achieve surprisingly and perhaps fortuitously good agreement with the activity and phenomenological pressure dependence when it is assumed that the (111) facet is the only active facet present. Lastly, the active phase of catalytic CO oxidation over RuO 2 remains unsettled, but the present study presents proof of principle (and progress) toward more accurate multiscale modeling from electrons to reactors and new simulation results.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2017

A density functional theory based approach for predicting melting points of ionic liquids

Accurate prediction of melting points of ILs is important both from the fundamental point of view and from the practical perspective for screening ILs with low melting points and broadening their utilization in a wider temperature range. In this work, we present an ab initio approach to calculating melting points of ILs with known crystal structures and illustrate its application for a series of 11 ILs containing imidazolium/pyrrolidinium cations and halide/polyatomic fluoro-containing anions. The melting point is determined as a temperature at which the Gibbs free energy of fusion is zero. The Gibbs free energy of fusion can be expressed through the use of the Born-Fajans-Haber cycle via the lattice free energy of forming a solid IL from gaseous phase ions and the sum of the solvation free energies of ions comprising IL. Dispersion-corrected density functional theory (DFT) involving (semi)local (PBE-D3) and hybrid exchange-correlation (HSE06-D3) functionals is applied to estimate the lattice enthalpy, entropy, and free energy. The ions solvation free energies are calculated with the SMD-generic-IL solvation model at the M06-2X/6-31+G(d) level of theory under standard conditions. The melting points of ILs computed with the HSE06-D3 functional are in good agreement with the experimental data, with a mean absolute error of 30.5 K and a mean relative error of 8.5%. The model is capable of accurately reproducing the trends in melting points upon variation of alkyl substituents in organic cations and replacement one anion by another. The results verify that the lattice energies of ILs containing polyatomic fluoro-containing anions can be approximated reasonably well using the volume-based thermodynamic approach. However, there is no correlation of the computed lattice energies with molecular volume for ILs containing halide anions. Moreover, entropies of solid ILs follow two different linear relationships with molecular volume for halides and polyatomic fluoro-containing anions. As a result, continuous progress in predicting crystal structures of organic salts with halide anions will be a key factor for successful prediction of melting points with no prior knowledge of the crystal structure.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2015

The benzene radical anion: A computationally demanding prototype for aromatic anions

The benzene radical anion is studied with ab initio coupled-cluster theory in large basis sets. Unlike the usual assumption, we find that, at the level of theory investigated, the minimum energy geometry is non-planar with tetrahedral distortion at two opposite carbon atoms. The anion is well known for its instability to auto-ionization which poses computational challenges to determine its properties. Despite the importance of the benzene radical anion, the considerable attention it has received in the literature so far has failed to address the details of its structure and shape-resonance character at a high level of theory. Here, we examine the dynamic Jahn-Teller effect and its impact on the anion potential energy surface. We find that a minimum energy geometry of C{sub 2} symmetry is located below one D{sub 2h} stationary point on a C{sub 2h} pseudo-rotation surface. The applicability of standard wave function methods to an unbound anion is assessed with the stabilization method. The isotropic hyperfine splitting constants (A{sub iso}) are computed and compared to data obtained from experimental electron spin resonance experiments. Satisfactory agreement with experiment is obtained with coupled-cluster theory and large basis sets such as cc-pCVQZ.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Nuclear properties with semilocal momentum-space regularized chiral interactions beyond N 2 LO

We present a comprehensive investigation of few-nucleon systems as well as light and medium-mass nuclei up to A = 48 using the current Low Energy Nuclear Physics International Collaboration two-nucleon interactions in combination with the third-order (N 2 LO) three-nucleon forces. To address the systematic overbinding of nuclei starting from A ≈ 10 found in our earlier study utilizing the N 2 LO two- and three-nucleon forces, we take into account higher-order corrections to the two-nucleon potentials up through fifth order in chiral effective field theory. Here, the resulting Hamiltonian can be completely determined using the A = 3 binding energies and selected nucleon-deuteron cross sections as input. It is then shown to predict other nucleon-deuteron scattering observables and spectra of light p-shell nuclei, for which a detailed correlated truncation error analysis is performed, in agreement with experimental data. Moreover, the predicted ground state energies of nuclei in the oxygen isotopic chain from 14 O to 26 O as well as 40 Ca and 48 Ca show a remarkably good agreement with experimental values, given that the Hamiltonian is fixed completely from the A ≤ 3 data, once the fourth-order (N 3 LO) corrections to the two-nucleon interactions are taken into account. On the other hand, the charge radii are found to be underpredicted by ≈10% for the oxygen isotopes and by almost 20% for 40 Ca and 48 Ca.

20 ≤ A ≤ 38↗
DOE OSTI2023

Search for an s -wave resonance in Li 7 just above the proton-decay threshold

Near-threshold resonances play an important role in nucleosynthesis and applied nuclear science. The study of nuclei removed from stability has greatly extended the list of resonances close to decay thresholds. The no core shell model with continuum (NCSMC) recently predicted an s-wave resonance just above the proton-decay threshold of 7 Li at an excitation energy of 10.2 MeV. This potential case of a near-threshold resonance is dependent on the quantum mechanics of the p + 6 He fragments extended into the continuum. Further, the 6 He(d,n) 7 Li* reaction was employed in an attempt to populate this resonance and search for its proton decay via the invariant-mass method. No evidence of this resonance was found. However, the data collected in this search led to better constraints on the energy and width of the isobaric analog state (IAS) with (J π ,T) = (3/2 - , 3/2) and revealed a new weak resonance just above the IAS in energy, predicted by the NCSMC as a (J π , T ) = (3/2 – , 1/2) state and potentially a part of the antianalog structure.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Liquid–Vapor Phase Equilibrium in Molten Aluminum Chloride (AlCl 3 ) Enabled by Machine Learning Interatomic Potentials

Molten salts are promising candidates in numerous clean energy applications, where knowledge of thermophysical properties and vapor pressure across their operating temperature ranges is critical for safe operations. Due to challenges in evaluating these properties using experimental methods, fast and scalable molecular simulations are essential to complement the experimental data. In this study, we developed machine learning interatomic potentials (MLIP) to study the AlCl 3 molten salt across varied thermodynamic conditions (T = 473–613 K and P = 2.7–23.4 bar), which allowed us to predict temperature-surface tension correlations and liquid–vapor phase diagram from direct simulations of two-phase coexistence in this molten salt. Two MLIP architectures, a Kernel-based potential and neural network interatomic potential (NNIP), were considered to benchmark their performance for AlCl 3 molten salt using experimental structure and density values. The NNIP potential employed in two-phase equilibrium simulations yields the critical temperature and critical density of AlCl 3 that are within 10 K (∼3%) and 0.03 g/cm 3 (∼7%) of the reported experimental values. An accurate correlation between temperature and viscosities is obtained as well. In doing so, we report that the inclusion of low-density configurations in their training is critical to more accurately represent the AlCl 3 system across a wide phase-space. The MLIP trained using PBE-D3 functional in the ab initio molecular dynamics (AIMD) simulations (120 atoms) also showed close agreement with experimentally determined molten salt structure comprising Al 2 Cl 6 dimers, as validated using Raman spectra and neutron structure factor. Furthermore, the PBE-D3 as well as its trained MLIP showed better liquid density and temperature correlation for AlCl 3 system when compared to several other density functionals explored in this work. Overall, the demonstrated approach to predict temperature correlations for liquid and vapor densities in this study can be employed to screen nuclear reactors-relevant compositions, helping to mitigate safety concerns.

Ab initio molecular dynamics↗
DOE OSTI2019

Pulse-duration dependence of the double-to-single ionization ratio of Ne by intense 780-nm and 800-nm laser fields: Comparison of simulations with experiments

Accurate ab initio calculations of the ratio of double-to-single ionization of Ne atoms in strong laser fields are difficult due to the many-electron nature of the target. Here, with accurate total cross sections carefully evaluated by using the state-of-the-art many-electron R-matrix theory for both electron-impact ionization and electron-impact excitation of Ne + , we simulate the total double-ionization yields of Ne 2+ in strong laser fields at 780 and 800 nm for pulse durations in the range from 7.5 to 200 fs based on the improved quantitative rescattering model. The corresponding single-ionization yields of Ne + are calculated within the nonadiabatic tunneling model of Perelomov, Popov, and Terent'ev. The ratio of double-to-single ionization of Ne is then obtained from the calculated double- and single-ionization yields. By normalizing the ratio to the one calculated from solving the time-dependent Schrödinger equation for a short few-cycle pulse, we make quantitative comparisons of our results with experimental data to show that our model predicts the experimental findings very well. Lastly, we analyze the pulse-duration dependence of the double-to-single ionization ratio.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2018

Theoretical and experimental investigation of the equation of state of boron plasmas

We report a theoretical equation of state (EOS) table for boron across a wide range of temperatures (5.1 10 4 {5.2 10 8 K) and densities (0.25{49 g/cm 3 ), and experimental shock Hugoniot data at unprecedented high pressures (5608 118 GPa). The calculations are performed with first-principles methods combining path integral Monte Carlo (PIMC) at high temperatures and density functional theory molecular dynamics (DFT-MD) methods at lower temperatures. PIMC and DFT-MD cross-validate each other by providing coherent EOS (difference <1.5 Hartree/boron in energy and <5% in pressure) at 5.1 10 5 K. The Hugoniot measurement is conducted at the National Ignition Facility using a planar shock platform. The pressure-density relation found in our shock experiment is on top of the shock Hugoniot profile predicted with our first-principles EOS and a semi-empirical EOS table (LEOS 50). We investigate the self diffusivity and the effect of thermal and pressure-driven ionization on the EOS and shock compression behavior in high pressure and temperature conditions. We also study the sensitivity of a polar direct-drive exploding pusher platform to pressure variations based on applying pressure multipliers to LEOS 50 and by utilizing a new EOS model based on our ab initio simulations via 1D radiation-hydrodynamic calculations. The results are valuable for future theoretical and experimental studies and engineering design in high energy density research.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Effect of Oxygen on the Quantum, Magnetic, and Thermodynamic Properties of Co Nanowires on the Reconstructed Anisotropic (1 × 2)/Au(110) and (1 × 2)/Pt(110) Surfaces: Ab Initio Approach

Ab initio theoretical study of the quantum magnetic properties of Co nanowires on the pure and oxygen-reconstructed (1 × 2)/Au(110) and (1 × 2)/Pt(110) surfaces is performed. Their structures and electronic configurations are calculated using the electron density functional theory. High values of magnetic moment and magnetic anisotropy energies of Co atoms are found on both pure and oxygen-reconstructed (1 × 2)/Au(110) and (1 × 2)/Pt(110) surfaces. The adsorption of oxygen atoms on the (1 × 2)/Au(110) substrate is shown to affect the structural arrangement of Co nanowire atoms on this substrate and to increase the magnetic anisotropy energy (by 1.91 meV per nanowire atom). The adsorption of oxygen on the Pt(110) substrate substantially decreases the magnetic anisotropy energy of the Co nanowire on it (by 5.98 meV per atom). The origin of these changes is revealed by analyzing the local densities of states of the d electrons of nanowire atoms. The temperature ranges of the states with the lowest free surface energy are determined using the atomistic thermodynamics methods. These data and the available experimental data are used to predict the possibility of observing the structures under study in experiments.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Direct-methods structure determination of a trypanosome RNA-editing substrate fragment with translational pseudosymmetry

Using direct methods starting from random phases, the crystal structure of a 32-base-pair RNA (675 non-H RNA atoms in the asymmetric unit) was determined using only the native diffraction data (resolution limit 1.05 Å) and the computer program SIR 2014. The almost three helical turns of the RNA in the asymmetric unit introduced partial or imperfect translational pseudosymmetry (TPS) that modulated the intensities when averaged by the l Miller indices but still escaped automated detection. Almost six times as many random phase sets had to be tested on average to reach a correct structure compared with a similar-sized RNA hairpin (27 nucleotides, 580 non-H RNA atoms) without TPS. Lastly, more sensitive methods are needed for the automated detection of partial TPS.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2016

Rate constants of chemical reactions from semiclassical transition state theory in full and one dimension

Semiclassical Transition State Theory (SCTST), a method for calculating rate constants of chemical reactions, offers gains in computational efficiency relative to more accurate quantum scattering methods. In full-dimensional (FD) SCTST, reaction probabilities are calculated from third and fourth potential derivatives along all vibrational degrees of freedom. However, the computational cost of FD SCTST scales unfavorably with system size, which prohibits its application to larger systems. In this study, the accuracy and efficiency of 1-D SCTST, in which only third and fourth derivatives along the reaction mode are used, are investigated in comparison to those of FD SCTST. Potential derivatives are obtained from numerical ab initio Hessian matrix calculations at the MP2/cc-pVTZ level of theory, and Richardson extrapolation is applied to improve the accuracy of these derivatives. Reaction barriers are calculated at the CCSD(T)/cc-pVTZ level. Results from FD SCTST agree with results from previous theoretical and experimental studies when Richardson extrapolation is applied. Results from our implementation of 1-D SCTST, which uses only 4 single-point MP2/cc-pVTZ energy calculations in addition to those for conventional TST, agree with FD results to within a factor of 5 at 250 K. This degree of agreement and the efficiency of the 1-D method suggest its potential as a means of approximating rate constants for systems too large for existing quantum scattering methods.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1997

On the Role of Resonances in Nonradial Pulsators

Resonances or near resonances are ubiquitous among the excited nonradial pulsation modes of variable stars and they must play an important role in determining their pulsational behavior. Here in a first step at nonlinear asteroseismology, we explore some of the consequences of resonances by means of the amplitude equation formalism. We show how parity and angular momentum constraints can be used to eliminate many of the possible nonlinear resonant couplings between modes (and multiplets of modes), and how the amplitude equations can thus be simplified. Even when we may not be able, nor wish, to make an ab initio computation of the values of the coupling coefficients, it is still possible to obtain constraints on the nature of the excited modes if a resonance between observed frequencies can be identified. Resonances can cause nonlinear frequency locking of modes. This means that the observed frequencies appear in exact resonance even though the linear frequencies are only approximately in resonance. The nonlinear frequency lock, when it occurs, it does so over a range of departures from linear resonance, and it is accompanied by constant pulsation amplitudes. The locked, nonlinear frequencies can differ noticeably from their nonresonant counterparts which are usually used in seismology. This is particularly true for multiplets of modes split by rotation. Beyond the regime of the frequency lock, amplitude and frequency modulations can appear in the pulsations. Far from the resonance condition one recovers the regime of steady pulsations with nonresonant frequencies for which the seismological studies, as they are presently carried out, are justified (provided furthermore, of course, that nonlinear frequency shifts are negligible). Success in identifying a resonance in an observed power spectrum depends on the quality of the data. While keeping this limitation in mind, ew discuss the possible existence of peculiar resonances the pulsations specific variable white dwarfs and delta Scuti stars.

Buchler, J. R.↗
DOE OSTI2016

Structural, electronic and magnetic properties of carbon doped boron nitride nanowire: Ab initio study

Using spin-polarized density functional theory calculations, we demonstrated that carbon doped boron nitride nanowire (C-doped BNNW) has diverse electronic and magnetic properties depending on position of carbon atoms and their percentages. Our results show that only when one carbon atom is situated on the edge of the nanowire, C-doped BNNW is transformed into half-metal. The calculated electronic structure of the C-doped BNNW suggests that doping carbon can induce localized edge states around the Fermi level, and the interaction among localized edge states leads to semiconductor to half-metal transition. Overall, the bond reconstruction causes of appearance of different electronic behavior such as semiconducting, half-metallicity, nonmagnetic metallic, and ferromagnetic metallic characters. The formation energy of the system shows that when a C atom is doped on surface boron site, system is more stable than the other positions of carbon impurity. Our calculations show that C-doped BNNW may offer unique opportunities for developing nanoscale spintronic materials.

36 MATERIALS SCIENCE↗
DOE OSTI2019

First principle study of magnetism and vacancy energetics in a near equimolar NiFeMnCr high entropy alloy

We report the results of ab initio calculations of a novel NiFeMnCr high entropy alloy (HEA) with potential applications as a high performance structural material. The bulk and defect property variations due to chemical disordering and magnetic frustration have been studied using both supercell and coherent potential approximation-based techniques. While magnetic frustration due to the presence of multiple 3d transition metals can severely affect the accuracy of vacancy formation energy in first-principles calculations, this effect should be suppressed at intermediate and high temperatures. An efficient approach to evaluate the chemical potential in HEA is constructed and implemented. Vacancy formation energies are computed based on the chemical potential. The statistical distribution of formation energies is weakly dependent upon the chemical identity of the vacancy. On the other hand, the calculated vacancy migration energies show that Fe is more likely to have a large migration barrier than Cr, Mn, or Ni. Finally, atomic-level stresses are computed. A qualitative model to explain the elemental segregation trend in HEA is built upon the atomic-level stress calculation results and provides a reasonable qualitative agreement with ion irradiation experimental data of a NiFeMnCr HEA.

36 MATERIALS SCIENCE↗
DOE OSTI2013

A new ab initio potential energy surface for the collisional excitation of HCN by para- and ortho-H{sub 2}

We present a new four-dimensional potential energy surface for the collisional excitation of HCN by H{sub 2}. Ab initio calculations of the HCN–H{sub 2} van der Waals complex, considering both molecules as rigid rotors, were carried out at the explicitly correlated coupled cluster with single, double, and perturbative triple excitations [CCSD(T)-F12a] level of theory using an augmented correlation-consistent triple zeta (aVTZ) basis set. The equilibrium structure is linear HCN–H{sub 2} with the nitrogen pointing towards H{sub 2} at an intermolecular separation of 7.20 a{sub 0}. The corresponding well depth is −195.20 cm{sup −1}. A secondary minimum of −183.59 cm{sup −1} was found for a T-shape configuration with the H of HCN pointing to the center of mass of H{sub 2}. We also determine the rovibrational energy levels of the HCN–para-H{sub 2} and HCN–ortho-H{sub 2} complexes. The calculated dissociation energies for the para and ortho complexes are 37.79 cm{sup −1} and 60.26 cm{sup −1}, respectively. The calculated ro-vibrational transitions in the HCN–H{sub 2} complex are found to agree by more than 0.5% with the available experimental data, confirming the accuracy of the potential energy surface.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2015

Quantum calculation of inelastic co collisions with H. III. Rate coefficients for RO-vibrational transitions

We present calculated rate coefficients for ro-vibrational transitions of CO in collisions with H atoms for a gas temperature range of 10 K ≤ T ≤ 3000 K, based on the recent three-dimensional ab initio H–CO interaction potential of Song et al. Rate coefficients for ro-vibrational v=1,j=0−30→v{sup ′}=0,j{sup ′} transitions were obtained from scattering cross sections previously computed with the close-coupling (CC) method by Song et al. Combining these with the rate coefficients for vibrational v=1−5→v{sup ′}<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→v{sup ′},j{sup ′} de-excitation. Cross sections and rate coefficients for ro-vibrational v=2,j=0−30→v{sup ′}=1,j{sup ′} transitions calculated with the CC 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. The 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{sub 2} ratio of 1; larger changes occur for lower temperatures.

79 ASTRONOMY AND ASTROPHYSICS↗