Rotational Quenching of H\u3csub\u3e2\u3c/sub\u3eO by He: Mixed Quantum/Classical Theory and Comparison with Quantum Results

The mixed quantum/classical theory (MQCT) formulated in the space-fixed reference frame is used to compute quenching cross sections of several rotationally excited states of water molecule by impact of He atom in a broad range of collision energies, and is tested against the full-quantum calculations on the same potential energy surface. In current implementation of MQCT method, there are two major sources of errors: one affects results at energies below 10 cm−1, while the other shows up at energies above 500 cm−1. Namely, when the collision energy E is below the state-to-state transition energy ΔE the MQCT method becomes less accurate due to its intrinsic classical approximation, although employment of the average-velocity principle (scaling of collision energy in order to satisfy microscopic reversibility) helps dramatically. At higher energies, MQCT is expected to be accurate but in current implementation, in order to make calculations computationally affordable, we had to cut off the basis set size. This can be avoided by using a more efficient body-fixed formulation of MQCT. Overall, the errors of MQCT method are within 20% of the full-quantum results almost everywhere through four-orders-of-magnitude range of collision energies, except near resonances, where the errors are somewhat larger

Mixed Quantum/Classical Theory for Inelastic Scattering of Asymmetric-top-rotor + Atom in the Body-fixed Reference Frame and Application to the H\u3csub\u3e2\u3c/sub\u3eO + He System

The mixed quantum/classical theory (MQCT) for inelastic molecule-atom scattering developed recently [A. Semenov and D. Babikov, J. Chem. Phys.139, 174108 (2013)] is extended to treat a general case of an asymmetric-top-rotor molecule in the body-fixed reference frame. This complements a similar theory formulated in the space-fixed reference-frame [M. Ivanov, M.-L. Dubernet, and D. Babikov, J. Chem. Phys.140, 134301 (2014)]. Here, the goal was to develop an approximate computationally affordable treatment of the rotationally inelastic scattering and apply it to H2O + He. We found that MQCT is somewhat less accurate at lower scattering energies. For example, below E = 1000 cm−1 the typical errors in the values of inelastic scattering cross sections are on the order of 10%. However, at higher scattering energies MQCT method appears to be rather accurate. Thus, at scattering energies above 2000 cm−1 the errors are consistently in the range of 1%–2%, which is basically our convergence criterion with respect to the number of trajectories. At these conditions our MQCT method remains computationally affordable. We found that computational cost of the fully-coupled MQCT calculations scales as n 2, where n is the number of channels. This is more favorable than the full-quantum inelastic scattering calculations that scale as n 3. Our conclusion is that for complex systems (heavy collision partners with many internal states) and at higher scattering energies MQCT may offer significant computational advantages

New model for datasets citation and extraction reproducibility in VAMDC

In this paper we present a new paradigm for the identification of datasets extracted from the Virtual Atomic and Molecular Data Centre (VAMDC) e-science infrastructure. Such identification includes information on the origin and version of the datasets, references associated to individual data in the datasets, as well as timestamps linked to the extraction procedure. This paradigm is described through the modifications of the language used to exchange data within the VAMDC and through the services that will implement those modifications. This new paradigm should enforce traceability of datasets, favour reproducibility of datasets extraction, and facilitate the systematic citation of the authors having originally measured and/or calculated the extracted atomic and molecular data.Comment: 48 page

Rate Coefficients for the Collisional Excitation of Molecules: Estimates from an Artificial Neural Network

An artificial neural network (ANN) is investigated as a tool for estimating rate coefficients for the collisional excitation of molecules. The performance of such a tool can be evaluated by testing it on a dataset of collisionally-induced transitions for which rate coefficients are already known: the network is trained on a subset of that dataset and tested on the remainder. Results obtained by this method are typically accurate to within a factor ~ 2.1 (median value) for transitions with low excitation rates and ~ 1.7 for those with medium or high excitation rates, although 4% of the ANN outputs are discrepant by a factor of 10 more. The results suggest that ANNs will be valuable in extrapolating a dataset of collisional rate coefficients to include high-lying transitions that have not yet been calculated. For the asymmetric top molecules considered in this paper, the favored architecture is a cascade-correlation network that creates 16 hidden neurons during the course of training, with 3 input neurons to characterize the nature of the transition and one output neuron to provide the logarithm of the rate coefficient.Comment: 23 pages including 9 figures. Accepted for publication in Ap

Collisional excitation of water by hydrogen atoms

We present quantum dynamical calculations that describe the rotational excitation of H$_2$O due to collisions with H atoms. We used a recent, high accuracy potential energy surface, and solved the collisional dynamics with the close-coupling formalism, for total energies up to 12 000 cm$^{-1}$. From these calculations, we obtained collisional rate coefficients for the first 45 energy levels of both ortho- and para-H$_2$O and for temperatures in the range T = 5-1500 K. These rate coefficients are subsequently compared to the values previously published for the H$_2$O / He and H$_2$O / H$_2$ collisional systems. It is shown that no simple relation exists between the three systems and that specific calculations are thus mandatory

Ro-vibrational relaxation of HCN in collisions with He: Rigid bender treatment of the bending-rotation interaction

We present a new theoretical method to treat atom-rigid bender inelastic collisions at the Close Coupling level (RBCC) in the space fixed frame. The coupling between rotation and bending is treated exactly within the rigid bender approximation and we obtain the cross section for the rotational transition between levels belonging to different bending levels. The results of this approach are compared with those obtained when using the rigid bender averaged approximation (RBAA) introduced in our previous work dedicated to this system. We discuss the validity of this approximation and of the previous studies based on rigid linear HCN

On the influence of collisional rate coefficients on the water vapour excitation

Water is a key molecule in many astrophysical studies. Its high dipole moment makes this molecule to be subthermally populated under the typical conditions of most astrophysical objects. This motivated the calculation of various sets of collisional rate coefficients (CRC) for H$_2$O (with He or H$_2$) which are necessary to model its rotational excitation and line emission. We performed accurate non--local non--LTE radiative transfer calculations using different sets of CRC in order to predict the line intensities from transitions that involve the lowest energy levels of H$_2$O (E $<$ 900 K). The results obtained from the different CRC sets are then compared using line intensity ratio statistics. For the whole range of physical conditions considered in this work, we obtain that the intensities based on the quantum and QCT CRC are in good agreement. However, at relatively low H$_2$ volume density ($n$(H$_2$) $<$ 10$^7$ cm$^{-3}$) and low water abundance ($\chi$(H$_2$O) $<$ 10$^{-6}$), these physical conditions being relevant to describe most molecular clouds, we find differences in the predicted line intensities of up to a factor of $\sim$ 3 for the bulk of the lines. Most of the recent studies interpreting early Herschel Space Observatory spectra used the QCT CRC. Our results show that although the global conclusions from those studies will not be drastically changed, each case has to be considered individually, since depending on the physical conditions, the use of the QCT CRC may lead to a mis--estimate of the water vapour abundance of up to a factor of $\sim$ 3

The Excitation of N2_2H+^+ in Interstellar Molecular Clouds. I - Models

We present LVG and non-local radiative transfer calculations involving the rotational and hyperfine structure of the spectrum of N$_2$H$^+$ with collisional rate coefficients recently derived by us. The goal of this study is to check the validity of the assumptions made to treat the hyperfine structure and to study the physical mechanisms leading to the observed hyperfine anomalies. We find that the usual hypothesis of identical excitation temperatures for all hyperfine components of the $J$=1-0 transition is not correct within the range of densities existing in cold dense cores, i.e., a few 10$^4$ \textless n(H$_2$) \textless a few 10$^6$ cm$^{-3}$. This is due to different radiative trapping effects in the hyperfine components. Moreover, within this range of densities and considering the typical abundance of N$_2$H$^+$, the total opacity of rotational lines has to be derived taking into account the hyperfine structure. The error made when only considering the rotational energy structure can be as large as 100%. Using non-local models we find that, due to saturation, hyperfine anomalies appear as soon as the total opacity of the $J$=1-0 transition becomes larger than $\simeq$ 20. Radiative scattering in less dense regions enhance these anomalies, and particularly, induce a differential increase of the excitation temperatures of the hyperfine components. This process is more effective for the transitions with the highest opacities for which emerging intensities are also reduced by self-absorption effects. These effects are not as critical as in HCO$^+$ or HCN, but should be taken into account when interpreting the spatial extent of the N$_2$H$^+$ emission in dark clouds.Comment: 13 pages, 12 figure