HeH+ Collisions with H2: Rotationally Inelastic Cross Sections and Rate Coefficients from Quantum Dynamics at Interstellar Temperatures.

We report for the first time an accurate ab initio potential energy surface for the HeH+-H2 system in four dimensions (4D) treating both diatomic species as rigid rotors. The computed ab initio potential energy point values are fitted using an artificial neural network method and used in quantum close coupling calculations for different initial states of both rotors, in their ground electronic states, over a range of collision energies. The state-to-state cross section results are used to compute the rate coefficients over a range of temperatures relevant to interstellar conditions. By comparing the four dimensional quantum results with those obtained by a reduced-dimensions approach that treats the H2 molecule as an averaged, nonrotating target, it is shown that the reduced dimensionality results are in good accord with the four dimensional results as long as the HeH+ molecule is not initially rotationally excited. By further comparing the present rate coefficients with those for HeH+-H and for HeH+-He, we demonstrate that H2 molecules are the most effective collision partners in inducing rotational excitation in HeH+ cation at interstellar temperatures. The rotationally inelastic rates involving o-H2 and p-H2 excitations are also obtained and they turn out to be, as in previous systems, orders of magnitude smaller than those involving the cation. The results for the H2 molecular partner clearly indicate its large energy-transfer efficiency to the HeH+ system, thereby confirming its expected importance within the kinetics networks involving HeH+ in interstellar environments.

Efficiency of rovibrational cooling of HeH+ by collisions with He: Cross sections and rate coefficients from quantum dynamics.

By extending an earlier study [Gianturco et al., J. Chem. Phys. 154, 054311 (2021)] on the purely rotational excitation of HeH+ by He atoms, we report in this paper integral cross sections and rate coefficients for rovibrational excitation and de-excitation processes in HeH+ due to collisions with He. The data were obtained using a new ab initio potential energy surface that includes the vibrational degree of freedom. The results are compared with those computed using the earlier potential energy surface by Panda and Sathyamurthy [J. Phys. Chem. A 107, 7125 (2003)] that additionally accounts for the proton-exchange reaction between HeH+ and He. It is shown that the exchange channel contributes nearly as much as the inelastic channel to the vibrational excitation and de-excitation processes and that the total rate constants pertaining to the purely inelastic processes are largely of the same magnitude as those obtained when both inelastic and reactive channels are included in the dynamics. The inelastic rovibrational rate coefficients involving this astrophysical cation are also found to be much larger than those obtained for anions present in similar interstellar environments.

Dynamics of HeHHe+ Rotational State Changes Induced by Collision with He: A Possible New Path in Early Universe Chemistry.

Ab initio calculations are employed to generate the rigid rotor (RR) potential energy surface (PES) describing the interaction of the linear molecular cation HeHHe+, at its equilibrium geometry, with the neutral He atom. The resulting interaction is employed to investigate the efficiency of rotational state-changing collisions at the temperatures relevant to the early universe conditions, where the latter molecule has been postulated to exist, albeit not yet observed. The inelastic rate coefficients are found to be fairly large and are compared with those found for another important cation just recently observed in the interstellar medium: the HeH+ polar molecule. The possibility for this cation to provide new options to energy dissipation routes under early universe conditions after the recombination era is briefly discussed.

Energy-transfer quantum dynamics of HeH+ with He atoms: Rotationally inelastic cross sections and rate coefficients.

Two different ab initio potential energy surfaces are employed to investigate the efficiency of the rotational excitation channels for the polar molecular ion HeH+ interacting with He atoms. We further use them to investigate the quantum dynamics of both the proton-exchange reaction and the purely rotational inelastic collisions over a broad range of temperatures. In current modeling studies, this cation is considered to be one of the possible cooling sources under early universe conditions after the recombination era and has recently been found to exist in the interstellar medium. The results from the present calculations are able to show the large efficiency of the state-changing channels involving rotational states of this cation. In fact, we find them to be similar in size and behavior to the inelastic and reaction rate coefficients obtained in previous studies, where H atoms were employed as projectiles. The same rotational excitation processes, occurring when free electrons are the collision partners of this cation, are also compared with the present findings. The relative importance of the reactive, proton-exchange channel and the purely inelastic channels is also analyzed and discussed. The rotational de-excitation processes are also investigated for the cooling kinetics of the present cation under cold trap conditions with He as the buffer gas. The implications of the present results for setting up more comprehensive numerical models to describe the chemical evolution networks in different environments are briefly discussed.

On the Formation of Interstellar CH- Anions: Exploring Mechanism and Rates for CH2 Reacting with H.

We present accurate ab initio calculations on the structural properties of a gas-phase reaction of possible interest for Saturn's outer atmosphere chemistry, in which the CH2 molecule has been detected. In the present study, that molecule is made to react with the H- anion to form the CH- species, one considered as a possible intermediate in ionic processes networks. The results indicate that this reaction is markedly exothermic and proceeds with the formation of an intermediate, which occurs via only a shallow barrier below the reagents and progresses directly to the product region. The corresponding rate coefficients of reactions are also computed by making use of the variational transition state theory modeling and found to efficiently lead to the formation of the final anion even at the lower temperatures of interstellar medium conditions.

Modeling state-selective photodetachment in cold ion traps: Rotational state "crowding" in small anions.

Using accurate ab initio calculations of the interaction forces, we employ a quantum mechanical description of the collisional state-changing processes that occur in a cold ion trap with He as a buffer gas. We generate the corresponding inelastic rates for rotational transitions involving three simple molecular anions OH-(1Σ), MgH-(1Σ), and C2H-(1Σ) colliding with the helium atoms of the trap. We show that the rotational constants of these molecular anions are such that within the low-temperature regimes of a cold ion trap (up to about 50 K), a different proportion of molecular states are significantly populated when loading helium as a buffer gas in the trap. By varying the trap operating conditions, population equilibrium at the relevant range of temperatures is reached within different time scales. In the modeling of the photodetachment experiments, we analyze the effects of varying the chosen values for photodetachment rates as well as the laser photon fluxes. Additionally, the changing of the collision dynamics under different buffer gas densities is examined and the best operating conditions, for the different anions, for yielding higher populations of specific rotational states within the ion traps are extracted. The present modeling thus illustrates possible preparation of the trap conditions for carrying out more efficiently state-selected experiments with the trapped anions.

Modeling Ionic Reactions at Interstellar Temperatures: The Case of NH2- + H2 ⇔ NH3 + H.

Structural features and enthalpy details are presented for the title reactions, both for the exothermic (forward) path to NH3 formation and for the endothermic (reverse) reaction to NH2- formation. Both pathways have relevance for the nitrogen chemistry in the interstellar medium (ISM). They are also helpful to document the possible role of H- in molecular clouds at temperatures well below room temperature. The structural calculations are carried out using different ab initio methods and are further employed to obtain the reaction rates down to the interstellar temperatures detected in earlier experiments. The reaction rates are obtained from the computed minimum energy path (MEP) using the variational transition-state theory (VTST) approach. The results indicate very good accord with experiment results at room temperature, while measured low temperature data down to 8 K are well described using an appropriately modified VTST approach. This is done by employing a temperature-dependent scaling, from room temperature conditions down to the lower ISM temperatures, which acknowledges the noncanonical behavior of the fast, barrierless exothermic reaction. The reasons for this behavior and the need for improving on the VTST method when used away from room temperatures are discussed.

N2+(2Σg) and Rb(2S) in a hybrid trap: modeling ion losses from radiative association paths.

By employing ab initio computed intermolecular potential energy surfaces we calculate the radiative association probabilities and rates for two different associative mechanisms involving trapped molecular ions N2+(2Σg) interacting either directly with ultracold Rb atoms or undergoing charge-exchange (CE) processes leading to the formation of complexes of the strongly exothermic products N2(X1Σg) plus Rb+(1S0). The two processes are expected to provide possible paths to ion losses in the trap within the timescale of experiments. The present calculations suggest that the associative rates for the 'vibrational' direct process are too small to be of any significant importance at the millikelvin temperatures considered in the experiments, while the 'vibronic' path into radiatively associating the CE products has a probability of occurring which is several orders of magnitude larger. However the reaction rate constants attributed to non-adiabatic CE [F. H. J. Hall and S. Willist, Phys. Rev. Lett., 2012, 109, 233202] are in turn several orders of magnitude larger than the radiative ones calculated here, thereby making the primary experimental process substantially unaffected by the radiative losses channel.

Rotational 'cooling' and 'heating' of OH+(3Σ-) by collisions with He: quantum dynamics revealing propensity rules under ion trap conditions.

Multichannel scattering calculations are presented for the low-energy collisions of the OH+ cation and He atoms, using an ab initio evaluation of the interaction potential, which had been obtained in earlier work, and a time-independent, multichannel treatment of the quantum dynamics carried out in this study using our in-house scattering code ASPIN. Given the presence of spin-rotation coupling effects, within an essentially electrostatic formulation of the interaction forces with He atoms in the trap, the ensuing propensity rules which control the relative size of the state-changing cross sections and of the corresponding inelastic rates, also computed at the most likely temperatures in an ion trap, are presented and analysed in detail.

Modeling Quantum Kinetics in Ion Traps: State-changing Collisions for OH+ (3Σ- ) Ions with He as a Buffer Gas.

We present quantum scattering calculations for rotational state-changing cross sections and rates, up to about 50 K of ion translational temperatures, for the OH+ molecular ion in collision with He atoms as the buffer gas in the trap. The results are obtained both by using the correct spin-rotation coupling of angular momenta and also within a recoupling scheme that treats the molecular target as a pseudo-(1Σ+ ) state, and then compares our findings with similar data for the OH- (1Σ+ ) molecular partner under the same conditions. This comparison intends to link the cation behaviour to the one found earlier for the molecular anion. The full calculations including the spin-rotation angular momenta coupling effects have been recently reported (L. González-Sánchez and R. Wester and F.A. Gianturco, Mol.Phys.2018, DOI 10.1080/00268976.2018.1442597) with the aim of extracting specific propensity rules controlling the size of the cross sections. The present study is instead directed to modelling trap cooling dynamics by further obtaining the solutions of the corresponding kinetics equations under different trap schemes so that, using the presently computed rates can allow us to indicate specific optimal conditions for the experimental setup of the collisional rotational cooling in an ion trap for the system under study.

Rotationally inelastic collisions of H2+ ions with He buffer gas: Computing cross sections and rates.

We present quantum calculations for the inelastic collisions between H2+ molecules, in rotationally excited internal states, and He atoms. This work is motivated by the possibility of experiments in which the molecular ions are stored and translationally cooled in an ion trap and a He buffer gas is added for deactivation of the internal rotational population, in particular at low (cryogenic) translational temperatures. We carry out an accurate representation of the forces at play from an ab initio description of the relevant potential energy surface, with the molecular ion in its ground vibrational state, and obtain the cross sections for state-changing rotationally inelastic collisions by solving the coupled channel quantum scattering equations. The presence of hyperfine and fine structure effects in both ortho- and para-H2+ molecules is investigated and compared to the results where such a contribution is disregarded. An analysis of possible propensity rules that may predict the relative probabilities of inelastic events involving rotational state-changing is also carried out, together with the corresponding elastic cross sections from several initial rotational states. Temperature-dependent rotationally inelastic rates are then computed and discussed in terms of relative state-changing collisional efficiency under trap conditions. The results provide the essential input data for modeling different aspects of the experimental setups which can finally produce internally cold molecular ions interacting with a buffer gas.

Forming NCO(-) in Dense Molecular Clouds: Possible Gas-Phase Chemical Paths From Quantum Calculations.

The existence of NCO(-) anions in the interstellar medium (ISM) has been suggested and searched for over the years but without any formal definitive sighting of that molecule. We discuss in this work the possible formation of either NCO(-) directly or of NCO neutral as a precursor to NCO(-) formation by electron attachment. We follow simple, gas-phase chemical reactions for which the general features are obtained from accurate quantum calculations. The results are shedding some additional light on the likely presence of this anion in the ISM environment, drawing further information from the specific features of the considered reactions on the additional chemical options that exist for its formation.

Nucleophilic substitution with two reactive centers: The CN(-) + CH3I case.

The nucleophilic substitution reaction CN(-) + CH3I allows for two possible reactive approaches of the reactant ion onto the methyl halide, which lead to two different product isomers. Stationary point calculations predict a similar shape of the potential and a dominant collinear approach for both attacks. In addition, an H-bonded pre-reaction complex is identified as a possible intermediate structure. Submerged potential energy barriers hint at a statistical formation process of both CNCH3 and NCCH3 isomers at the experimental collision energies. Experimental angle- and energy differential cross sections show dominant direct rebound dynamics and high internal excitation of the neutral product. No distinct bimodal distributions can be extracted from the velocity images, which impedes the indication of a specific preference towards any of the product isomers. A forward scattering simulation based on the experimental parameters describes accurately the experimental outcome and shows how the possibility to discriminate between the two isomers is mainly hindered by the large product internal excitation.

CH(+) Destruction by Reaction with H: Computing Quantum Rates To Model Different Molecular Regions in the Interstellar Medium.

A detailed analysis of an ionic reaction that plays a crucial role in the carbon chemistry of the interstellar medium (ISM) is carried out by computing ab initio reactive cross sections with a quantum method and by further obtaining the corresponding CH(+) destruction rates over a range of temperatures that shows good overall agreement with existing experiments. The differences found between all existing calculations and the very-low-T experiments are discussed and explored via a simple numerical model that links these cross section reductions to collinear approaches where nonadiabatic crossing is expected to dominate. The new rates are further linked to a complex chemical network that models the evolution of the CH(+) abundance in the photodissociation region (PDR) and molecular cloud (MC) environments of the ISM. The abundances of CH(+) are given by numerical solutions of a large set of coupled, first-order kinetics equations that employs our new chemical package krome. The analysis that we carry out reveals that the important region for CH(+) destruction is that above 100 K, hence showing that, at least for this reaction, the differences with the existing laboratory low-T experiments are of essentially no importance within the astrochemical environments discussed here because, at those temperatures, other chemical processes involving the title molecule are taking over. A detailed analysis of the chemical network involving CH(+) also shows that a slight decrease in the initial oxygen abundance might lead to higher CH(+) abundances because the main chemical carbon ion destruction channel is reduced in efficiency. This might provide an alternative chemical route to understand the reason why general astrochemical models fail when the observed CH(+) abundances are matched with the outcomes of their calculations.

Formation of cyanopolyyne anions in the interstellar medium: the possible role of permanent dipoles.

The possibility of attaching near-threshold electrons to N-terminated carbon chains, like those observed in the outer envelopes of carbon-rich stars, is examined via accurate quantum chemistry orbital structures evaluation and quantum scattering analysis of the corresponding extra-electron wavefunctions at meV energies. It is shown that the differences in the signs and sizes of the permanent dipole moments which exist for both the neutral and anionic species of the C(n)N series of molecules play a significant role in suggesting or excluding possible energy paths to permanent anion formations of cyanopolyynes, for which the cases with n from 1 to 7 are examined in more detail.

Electron scattering cross sections from anthracene over a broad energy range (0.00001-10,000 eV).

We report a computational investigation of electron scattering by anthracene (C14H10) in the gas phase. Integral and differential cross sections have been calculated by employing two distinct ab-initio quantum scattering methods: the symmetry adapted-single centre expansion method (ePOLYSCAT) and a screening corrected form of the independent atom model (IAM-SCAR) at low and high energies, respectively. After a detailed evaluation of the current results, we present a complete set of integral scattering cross sections from 0.00001 to 10,000 eV.

Electron scattering cross section calculations for polar molecules over a broad energy range.

We report computational integral and differential cross sections for electron scattering by two different polar molecules, HCN and pyrimidine, over a broad energy range. We employ, for low energies, either the single-centre expansion (ePOLYSCAT) or the R-matrix method, while for the higher energies we select a corrected form of the independent-atom representation (IAM-SCAR). We provide complete sets of integral electron scattering cross sections from low energies up to 10,000 eV. Our present calculated data agree well with prior experimental results.

Low-energy positron scattering from gas-phase tetrahydrofuran: a quantum treatment of the dynamics and a comparison with experiments.

In this paper we report new quantum calculations of the dynamics for low-energy positrons interacting with gaseous molecules of tetrahydrofuran. The new quantum scattering cross sections are differential and integral cross sections at collision energies between 1.0 and 25.0 eV and include a careful treatment of the additional effects on the scattering process brought about by the permanent dipole moment of the target molecule. The present results are compared with an extensive range of measured data, both for the angular distributions and for the elastic integral cross sections and agree remarkably well with all findings. The new calculated quantities reported here also show the importance of correcting the experimental integral cross sections for the angular discrimination in the forward direction.

Electron scattering from gas phase cis-diamminedichloroplatinum(II): quantum analysis of resonance dynamics.

We present scattering calculations of electron collisions with the platinum-containing compound cis-diamminedichloroplatinum (CDDP), commonly known as cisplatin, between 0.5 eV and 6 eV, and the corresponding isolated Pt atom from 0.1 eV to 10 eV. We find evidence of resonances in e(-)-CDDP scattering, using an ab initio description of the target. We computed scattering matrix elements from equations incorporating exchange and polarization effects through the use of the static-exchange plus density functional correlation potential. Additionally, we made use of a purely local adiabatic model potential that allows Siegert eigenstates to be calculated, thereby allowing inspection of the possible resonant scattering wave functions. The total cross section for electron scattering from (5d(10)) (1)S Pt displays a large magnitude, monotonic decay from the initial collision energies, with no apparent resonance scattering features in any scattering symmetry. By contrast, the e(-)-CDDP scattering cross section shows a small feature near 3.8 eV, which results from a narrow, well localized resonance of b2 symmetry. These findings are then related to the possible electron-mediated mechanism of the action of CDDP on DNA replication as suggested by recent experiments.

Electron scattering cross sections from HCN over a broad energy range (0.1-10,000 eV): Influence of the permanent dipole moment on the scattering process.

We report theoretical integral and differential cross sections for electron scattering from hydrogen cyanide derived from two ab initio scattering potential methods. For low energies (0.1-100 eV), we have used the symmetry adapted-single centre expansion method using a multichannel scattering formulation of the problem. For intermediate and high energies (10-10,000 eV), we have applied an optical potential method based on a screening corrected independent atom representation. Since HCN is a strong polar molecule, further dipole-induced excitations have been calculated in the framework of the first Born approximation and employing a transformation to a space-fixed reference frame of the calculated K-matrix elements. Results are compared with experimental data available in the literature and a complete set of recommended integral elastic, inelastic, and total scattering cross sections is provided from 0.1 to 10,000 eV.