Ab Initio Potential Energy Surface Research Articles

Ortho-Para Conversion for H+ + H2 Collision in Low Temperature: A Fully Close-Coupled Time-Dependent Wave Packet Study.

The collision-induced rate coefficients of ortho-para conversion for the H+ + H2 reaction provide accurate information to probe the lifetime of cold environments in interstellar media. Rotationally resolved reaction probabilities are calculated at the low collision energy regime (0 < Ecol ≤ 0.3 eV) by employing the coupled three-dimensional (3D) time-dependent wave packet (TDWP) formalism in hyperspherical coordinates on a recently constructed ab initio ground adiabatic potential energy surface of H3+ [J. Chem. Phys. 2014, 141, 204306] for the process H+ + H2 (v = 0, j = 0-5) → H+ + H2 (v' = 0, j'). Cross-sections are then computed from the converged reaction probabilities as a function of total angular momentum (J) over the same energy regime and subsequently employed to obtain the rate constants for the ortho-to-para (O-P) and para-to-ortho (P-O) conversions and their ratio. The ratio of ortho-para conversion shows (a) appropriate convergence with the inclusion of a higher number of initial rotational states as well as a reasonable agreement with the results from a quantum statistical method and (b) a peak at lower temperature that could be due to the available collision energy for transitions involving lower rotational states (j = 0 → j' = 1).

Coarse-grained modeling of high-enthalpy air flows based on the updated vibrational state-to-state kinetics

The state-to-state (StS) model can accurately describe high-temperature thermochemical nonequilibrium flows. For the five-species air gas mixture, we develop a comprehensive database for the state-specific rate coefficients for temperatures 300–25 000 K in this paper. The database incorporates recent molecular dynamics simulations (based on the ab initio potential energy surfaces) in the literature, and theoretical methods, including the forced harmonic oscillator model and the Marrone–Treanor model, are employed to complement the rate coefficients that are unavailable from molecular dynamics calculations. The post-shock StS simulations using the present database agree with the experimental NO infrared radiation. Based on this updated StS kinetics database, we investigate the post-shock high-enthalpy air flows by employing both the StS and coarse-grained models (CGM). The CGM, which lumps molecular vibrational states into groups, shows results that align with the StS model, even utilizing only two groups for each molecule. However, the CGM-1G model, with only one group per molecule and belonging to the multi-temperature model (but uses StS kinetics), fails to reproduce the StS results. Analysis of vibrational energy source terms for different kinetic processes and fractions of vibrational groups reveals that the deficiency of the CGM-1G model stems from the overestimation of high-lying vibrational states, leading to higher dissociation rates and increased consumption of vibrational energy in dissociation. Furthermore, the presence of the Zeldovich-exchange processes indirectly facilitates energy transfer in N2 and O2, a phenomenon not observed in binary gas systems. These findings have important implications for developing the reduced-order model based on coarse-grained treatment.

Inelastic scattering of PO+ by H2 at interstellar temperatures

ABSTRACT Phosphorous species are of great interest in interstellar chemistry since they are the basic blocks for building life here on Earth. Modelling the abundance and environment of recently detected PO$^{+}$ under non-local thermodynamic equilibrium (LTE) requires rotational spectra of the molecule along with accurate collisional rates with the most abundant species, hydrogen and helium. A new 4D ab initio potential energy surface (PES) of PO$^{+}$ – H$_{2}$ collision is calculated using CCSD(T)/CBS(DTQ) methodology considering rigid rotor approximation. The region containing the minima of the PES is augmented using neural networks (NNs) model while very high potentials ($\gt 2500$ cm$^{-1}$) and asymptotic region have been approximated using Slater and R$^{-4}$ functions, respectively. The close coupling calculations have been performed using molscat software for both ortho and para-H$_{2}$. The rate coefficients have been reported for transitions $j-j^{\prime }=$$1-0$, $2-1$, $3-2$, and $5-4$ through which PO$^{+}$ has been experimentally detected in interstellar medium (ISM). The rate coefficients for even and odd transitions of PO$^{+}$ with para-H$_{2}$ are compared with that of helium and are found to be 1.1–2.0 times higher. For even transitions ($\Delta j = 2$), the ortho-H$_{2}$ rates are 10 per cent higher than para-H$_{2}$ rates. However, the trend reverses in the case of odd transitions ($\Delta j = 1$) when higher J transitions are considered at low temperatures. At higher temperatures, the ortho rates cross the para-H$_{2}$ rates and become larger than the latter. The new rate coefficients with both ortho and para-H$_{2}$ will enable accurate modelling of the PO$^{+}$ abundance in the ISM under non-LTE conditions.

First Ab Initio Line Lists for Triatomic Sulfur-Containing Molecules: S2O and S3.

We present in this paper the first room temperature rovibrational line lists for the main isotopologues of disulfur monoxide (S2O) and thiozone (S3) variationally computed from newly developed ab initio potential energy and dipole moment surfaces (PESs and DMSs). S2O and S3 are supposed to be potential candidates for astronomical detection, especially in the Venusian atmosphere where sulfur chemistry plays a major role. Contrary to other stable sulfur-containing species like SO2 or H2S, there is much less experimental data for the short-lived reactive S2O and S3 molecules. Indeed, the infrared spectra of S2O and S3 are quite complicated to analyze without consistent theoretical predictions because they contain both high J (∼100) and a lot of hot band transitions, even at T = 296 K. In this work, we have constructed PESs based on the single-reference coupled cluster approach [CCSD(T)] and including corrections due to the scalar relativistic effects, DBOC, and highly excited Slater determinants. The structure of the excited electronic states of S3 was also discussed. The nuclear-motion Eckart-Watson Hamiltonian expressed in terms of normal-mode irreducible tensor operators was used to compute the energy levels. For line intensity calculation, ab initio DMSs were computed at the CCSD(T)/aug-cc-pV5Z level of the theory. Rotation-vibration patterns in the fundamental bands of S2O and S3 have been simulated taking into account the recent Fourier-transform spectra of S2O recorded at the SOLEIL synchrotron. The present work aims at providing line intensities of S2O and S3 that are missing in the literature as well as new spectroscopic support for atmospheric applications.

Dynamics of the Cl + CH3CN reaction on an automatically-developed full-dimensional abinitio potential energy surface.

A full-dimensional analytical potential energy surface (PES) is developed for the Cl + CH3CN reaction following our previous work on the benchmark abinitio characterization of the stationary points. The spin-orbit-corrected PES is constructed using the Robosurfer program and a fifth-order permutationally invariant polynomial method for fitting the high-accuracy energy points determined by a ManyHF-based coupled-cluster/triple-zeta-quality composite method. Quasi-classical trajectory simulations are performed at six collision energies between 10 and 60kcal mol-1. Multiple low-probability product channels are found, including isomerization to isonitrile (CH3NC), but out of the eight possible channels, only the H-abstraction has significant reaction probability; thus, detailed dynamics studies are carried out only for this reaction. The cross sections and opacity functions show that the probability of the H-abstraction reaction increases with increasing collision energy (Ecoll). Scattering angle, initial attack angle, and product relative translational energy distributions indicate that the mechanism changes with the collision energy from indirect/rebound to direct stripping. The distribution of initial attack angles shows a clear preference for methyl group attack but with different angles at different Ecoll values. Post-reaction energy distributions show that the energy transfer is biased toward the products' relative translational energy instead of their internal energy. Rotational and vibrational energy have about the same amount of contribution to the internal energy in the case of both products (HCl and CH2CN), i.e., both of them are formed with high rotational excitations. HCl is produced mostly in the ground vibrational state, while a notable fraction of CH2CN is formed with vibrational excitation.

Rovibrational Line Lists of Triplet and Singlet Methylene.

The methylene molecule (CH2) is a short-lived radical with lacking data on its spectral line intensities. Although the lifetime of CH2 is extremely short under Earth's conditions, it exists in a free form in interstellar media. CH2 is an important intermediate species in chemical reactions associated with the formation and destruction of complex hydrocarbons. We present the first rovibrational line lists of CH2 in its ground triplet and first excited singlet electronic state. To this end, our previously developed accurate ab initio potential energy surface (PES) was used for the ground electronic triplet state [Egorov et al. J. Comp. Chem. 2024. V. 45. (2). P. 83] while a new PES for the singlet state was constructed in this work using the single-reference coupled cluster approach [CCSD(T)] combined with the extrapolation to the complete basis set (CBS) limit based on the correlation-consistent orbital basis sets with the core-valence electron correlation effects [aug-cc-pCVXZ, X = T, Q, 5, and 6]. In addition, the contributions to the correlation energy from highly excited Slater determinants [CC(n), n = 3-5] were included as well as the scalar relativistic effects and DBOC. The most accurate description of the infrared band origins of singlet CH2 was thus achieved for the energy range where the impact of the nonadiabatic coupling due to the Renner-Teller effect can be neglected. To obtain the probabilities of the rovibrational transitions, new ab initio DMSs were constructed both for the triplet and singlet CH2 using the CCSD(T)/aug-cc-pCVQZ approach. Finally, the absorption spectra of triplet and singlet methylene were predicted from the variationally computed line lists.

Ab Initio Potentials for the Ground S0 and the First Electronically Excited Singlet S1 States of Benzene-Helium with Application to Tunneling Intermolecular Vibrational States.

We present new ab initio intermolecular potential energy surfaces for the benzene-helium complex in its ground (S0) and first excited (S1) states. The coupled-cluster level of theory with single, double, and perturbative triple excitations, CCSD(T), was used to calculate the ground state potential. The excited state potential was obtained by adding the excitation energies S0 → S1 of the complex, calculated using the equation of motion approach EOM-CCSD, to the ground state potential interaction energies. Analytical potentials are constructed and applied to study the structural and vibrational dynamics of benzene-helium. The binding energies and equilibrium distances of the ground and excited states were found to be 89 cm-1, 3.14 Å and 77 cm-1, 3.20 Å, respectively. The calculated vibrational energy levels exhibit tunneling of He through the benzene plane and are in reasonable agreement with recently reported experimental values for both the ground and excited states [Hayashi, M.; Ohshima, Y. J. Phys. Chem. Lett. 2020, 11, 9745]. Prospects for the theoretical study of complexes with large aromatic molecules and He are also discussed.

Accurate reference spectra of HD in an H2–He bath for planetary applications

Context. The hydrogen deuteride (HD) molecule is an important deuterium tracer in astrophysical studies. The atmospheres of gas giants are dominated by molecular hydrogen, and the simultaneous observation of H2 and HD lines provides reliable information on the D/H ratios on these planets. The reference spectroscopic parameters play a crucial role in such studies. Under the thermodynamic conditions encountered in these atmospheres, spectroscopic studies of HD require not only the knowledge of line intensities and positions but also accurate reference data on pressure-induced line shapes and shifts. Aims. Our aim is to provide accurate collision-induced line-shape parameters for HD lines that cover any thermodynamic conditions relevant to the atmospheres of giant planets, namely any relevant temperature, pressure, and perturbing gas composition (the H2–He mixture). Methods. We performed quantum-scattering calculations on our new, highly accurate ab initio potential energy surface (PES), and we used scattering S matrices obtained in this way to determine the collision-induced line-shape parameters. We used cavity ring-down spectroscopy to validate our theoretical methodology. Results. We report accurate collision-induced line-shape parameters for the pure rotational R(0), R(1), and R(2) lines, the most relevant HD lines for investigations of the atmospheres of the giant planets. Besides the basic Voigt-profile collisional parameters (i.e., the broadening and shift parameters), we also report their speed dependences and the complex Dicke parameter, which can influence the effective width and height of the HD lines up to almost a factor of 2 for giant planet conditions. The sub-percent-level accuracy reached in this work is a considerable improvement over previously available data. All the reported parameters (and their temperature dependences) are consistent with the HITRAN database format, hence allowing for the use of the HITRAN Application Programming Interface (HAPI) for generating the beyond-Voigt spectra of HD.

Characteristic Vibrational and Rotational Relaxation Times for Air Species from First-Principles Calculations

We present molecular-scale computational rotational-vibrational relaxation studies for N2, O2, and NO. Characteristic relaxation times for diatom-diatom and diatom-atom interactions are calculated using direct molecular simulation (DMS), with ab initio potential energy surfaces (PESs) as the sole model input. Below approximately 8000 K our N2−N2, O2−O2, and O2−N2 vibrational relaxation times agree well with the Millikan–White (M&W) correlation, but gradually diverge at higher temperatures. Park’s high-temperature correction produces a relatively steeper temperature rise compared to our estimates. DMS further shows that, with increasing temperature, the gap between vibrational and rotational relaxation times shrinks for all species. At T>30,000K their magnitudes become comparable and a clear distinction between both energy modes becomes meaningless. For other interactions, our DMS results differ substantially from the M&W correlation, both in magnitude and temperature dependence. Our predicted N2−O2 vibrational relaxation times are noticeably shorter due to vibration-vibration transfer. For O2−O we observe minimal temperature dependence. Our O2−N and N2−N predictions follow the M&W temperature trend at values roughly one order of magnitude smaller. For NO−NO, N2−O, NO−N, and NO−O we generate partial data due to currently incomplete PES sets. These first-principles-derived relaxation times are useful for informing relaxation models in gas-kinetic and fluid-dynamics simulations of high-enthalpy flows.

Dynamics of the HCl + C2H5 Multichannel Reaction on a Full-Dimensional Ab Initio Potential Energy Surface.

We report a full-dimensional ab initio analytical potential energy surface (PES), which accurately describes the HCl + C2H5 multichannel reaction. The new PES is developed by iteratively adding selected configurations along HCl + C2H5 quasi-classical trajectories (QCTs), thereby improving our previous Cl(2P3/2) + C2H6 PES using the Robosurfer program package. QCT simulations for the H'Cl + C2H5 reaction reveal hydrogen-abstraction, chlorine-abstraction, and hydrogen-exchange channels leading to Cl + C2H5H', H' + C2H5Cl, and HCl + C2H4H', respectively. Hydrogen abstraction dominates in the collision energy (Ecoll) range of 1-80 kcal/mol and proceeds with indirect isotropic scattering at low Ecoll and forward-scattered direct stripping at high Ecoll. Chlorine abstraction opens around 40 kcal/mol collision energy and becomes competitive with hydrogen abstraction at Ecoll = 80 kcal/mol. A restricted opening of the cone of acceptance in the Cl-abstraction reaction is found to result in the preference for a backward-scattering direct-rebound mechanism at all energies studied. Initial attack-angle distributions show mainly side-on collision preference of C2H5 for both abstraction reactions, and in the case of the HCl reactant, H/Cl-side preference for the H/Cl abstraction. For hydrogen abstraction, the collision energy transfer into the product translational and internal energy is almost equally significant, whereas in the case of chlorine abstraction, most of the available energy goes into the internal degrees of freedom. Hydrogen exchange is a minor channel with nearly constant reactivity in the Ecoll range of 10-80 kcal/mol.

Theoretical Insights into the Vibrational Structure of Carbon Dioxide Rare-Gas Complexes.

Two new flexible-monomer two-body ab initio potential energy surfaces (PESs) for the neon and krypton van der Waals complexes with carbon dioxide were developed, extending our previous work on the Ar-CO2 molecule. The accuracy of the PESs was validated by their agreement with the vibrational spectrum of the rare-gas complexes. The intermolecular and intramolecular vibrational excitation energies were computed at the vibrational self-consistent field and vibrational configuration interaction levels of theory. Overall, the agreement between theory and experiment is excellent throughout the vibrational spectra. The observed slight splitting of the bending modes, resulting from their nondegeneracy in the complexes, is confirmed by our computations, and the results qualitatively agree with the experiment. The splitting increases with increasing polarizability of the rare-gas atom. Additionally, we explain a discrepancy in the mode assignment in the intermolecular region of the neon complex with our VCI character assignment.

Feature selection for high-dimensional neural network potentials with the adaptive group lasso

Neural network potentials are a powerful tool for atomistic simulations, allowing to accurately reproduce ab initio potential energy surfaces with computational performance approaching classical force fields. A central component of such potentials is the transformation of atomic positions into a set of atomic features in a most efficient and informative way. In this work, a feature selection method is introduced for high dimensional neural network potentials, based on the adaptive group lasso (AGL) approach. It is shown that the use of an embedded method, taking into account the interplay between features and their action in the estimator, is necessary to optimize the number of features. The method’s efficiency is tested on three different monoatomic systems, including Lennard–Jones as a simple test case, Aluminium as a system characterized by predominantly radial interactions, and Boron as representative of a system with strongly directional components in the interactions. The AGL is compared with unsupervised filter methods and found to perform consistently better in reducing the number of features needed to reproduce the reference simulation data at a similar level of accuracy as the starting feature set. In particular, our results show the importance of taking into account model predictions in feature selection for interatomic potentials.

Quantum Dynamics of O(3P)+HBr→OH+Br Reaction: Integral Cross Sections and Rate Constants and their Initial State Dependence.

The integral cross sections (ICSs) and rate constants (RCs) for the title reaction are calculated by means of accurate quantum wave packet method employing the recent ab initio potential energy surface developed by Peterson and co-workers. Hundreds of partial wave contributions (up to J=150) are calculated explicitly taking Coriolis coupling into account. The ICSs are found to increase monotonically with energy and their energy dependences are smooth except in the energy range below the potential energy barrier, where several resonance peaks are observed. The temperature dependences of the RCs obey in general the Arrhenius law. These behaviors manifest the dominance of abstract mechanism. Initial rotational state (ji) excitations are found to greatly enhance the reactivity starting from ji=4, while exhibiting a complicated dependence for ji<4. From Boltzmann average of the initial state specified RCs (for ji=0-15) we obtained thermal RCs which are in reasonable agreement with available experimental results.

Global Potential Energy Surfaces by Compressed-State Multistate Pair-Density Functional Theory for Hyperthermal Collisions in the O2+O2 System.

Interactions between oxygen molecules play an important role in atmospheric chemistry and hypersonic flow chemistry in atmospheric entries. Recently, high-quality ab initio potential energy surface (PES) of the quintet O4 was reported by Paukku et al. [J. Chem. Phys. 147, 034301 (2017)]. 10543 configurations were sampled and calculated at the level of MS-CASPT2/maug-cc-pVTZ with scaled external correlation. The PES was fitted to a many-body (MB) form with the many-body part described by the permutationally invariant polynomial approach (MB-PIP). In this work, the PIP-Neural Network (PIP-NN) and MB-PIP-NN methods were used to refit the PES based on the same data by Paukku et al. Three PESs were compared. It was found that the performances differ significantly in the O+O3 region as well as in the long-range region. Therefore, additional 1300 points were sampled, and the efficient compressed-state multistate pair-density functional theory (CMS-PDFT) was used to calculate the electronic structure of these 1300 points and 10543 points by Paukku et al. Then, a completely new quintet PES was fitted using the MB-PIP-NN method. Based on this PES, the quasi-classical trajectory (QCT) approach was used to reveal all possible reaction channels for hyperthermal O2-O2 collisions.

State-to-state study of non-equilibrium recombination of oxygen and nitrogen molecules.

Rapidly cooled mixtures are of interest for several applications, including hypersonic flows due to the presence of strong cooling temperature gradients in regions such as hypersonic boundary layers and expanding nozzles. There have been very few studies of rapidly cooled mixtures using the high-fidelity rovibrational databases afforded by abinitio potential energy surfaces. This work makes use of existing rovibrational state-specific databases to study rapidly cooled mixtures. In particular, we seek to understand the importance of thermal non-equilibrium in recombining mixtures using both rovibrational and vibrational state-to-state methods for oxygen and nitrogen molecules. We find that although there is significant non-equilibrium during recombination, it is well captured by the vibrational state-specific approach. Finally, we compare the global recombination rate computed based on the state-specific recombination rate coefficients and the global recombination rate computed based on the time local dissociation rate coefficient, which is reversed using the principle of detailed balance. The local dissociation rate coefficient is computed by weighting the state-specific dissociation rate coefficients with the state-specific distribution of energy states. We find a large difference between these rates, highlighting a potential source of errors in hypersonic flow predictions.

A new ab initio potential energy surface and rovibrational spectra for the N2–N2O complex

We propose a new potential energy surface for the N2–N2O complex, which encompasses the Q3 normal mode for the ν3 antisymmetric stretching coordinate of N2O. The interaction potential energies were obtained at the explicitly correlated coupled cluster with single, double, and perturbative triple excitations [CCSD(T)-F12a] level with aug-cc-pVTZ basis set plus mid-bond functions. Four-dimensional vibrationally averaged potentials with N2O in the ground and ν3 excited states were generated through integration over the Q3 coordinate. Both potentials exhibit two equivalent skew T-shaped global minima. The rovibrational energies were calculated by employing the radial discrete variable representation/angular finite basis representation approach. The determined fundamental band origin shift of the infrared spectra in the N2O ν3 region, transition frequencies and rotational constants are all consistent with the experimental values.

A New Ab Initio Potential Energy Surface and Rovibrational Spectra for the CO-N2O Complex.

We constructed a new ab initio potential energy surface (PES) for CO-N2O which includes the intramolecular Q3 normal coordinate for the N2O ν3 antisymmetric stretching vibration. The intermolecular potential was evaluated employing the supermolecular method at the [CCSD(T)]-F12a level, with the aug-cc-pVTZ basis set plus bond functions. By integral over the intramolecular Q3 coordinate, we obtained the vibrationally averaged PESs for the CO-N2O system in the ground and ν3 excited states of N2O. Each PES features one nearly T-shaped global minimum and one skewed T-shaped local minimum. Based on these obtained PESs of CO-N2O, the radial discrete variable representation/angle finite base representation method and the Lanczos algorithm were applied for the calculations of bound states and rovibrational energy levels. The calculated ν3 vibrational band origin shift of the N2O monomer in CO-N2O is 2.7570 cm-1, matching well with the observed value of 2.9048 cm-1. The computed microwave and infrared transition frequencies, as well as the rotational parameters, are consistent with the experimental observations.

Non-IRC Mechanism of Bimolecular Reactions with Submerged Barriers: A Case Study of Si+ + H2O Reaction.

Chemical reactions with submerged barriers may feature interesting dynamic behaviors that are distinct from those with substantial barriers or those entirely dominated by capture. The Si+ + H2O reaction is a prototypical example, involving even two submerged saddle points along the reaction path: one for the direct dissociation of H (H-dissociation SP) and another for H migration from the O-side to the Si-side (H-migration SP). We investigated the intricacies of this process by employing quasi-classical trajectory calculations on an accurate, full-dimensional ab initio potential energy surface. Through careful trajectory analysis, an interesting nonintrinsic reaction coordinate mechanism was found to play an important role in producing SiOH+ and H. This new pathway is featured as that the H atoms do not form HSiOH+ complexes along the minimum-energy path but directly dissociate into the products after passing through the H-migration SP. Furthermore, based on artificially scaled potential energy surfaces (PESs), the impact of barrier height on the reaction is also explored. This work provides new insights into the dynamics of the Si+ + H2O reaction and enriches our understanding of reactions with submerged barriers.

Translational eigenstates of He@C60 from four-dimensional abinitio potential energy surfaces interpolated using Gaussian process regression.

We investigate the endofullerene system 3He@C60 with a four-dimensional potential energy surface (PES) to include the three He translational degrees of freedom and C60 cage radius. We compare second order Møller-Plesset perturbation theory (MP2), spin component scaled-MP2, scaled opposite spin-MP2, random phase approximation (RPA)@Perdew, Burke, and Ernzerhof (PBE), and corrected Hartree-Fock-RPA to calibrate and gain confidence in the choice of electronic structure method. Due to the high cost of these calculations, the PES is interpolated using Gaussian Process Regression (GPR), owing to its effectiveness with sparse training data. The PES is split into a two-dimensional radial surface, to which corrections are applied to achieve an overall four-dimensional surface. The nuclear Hamiltonian is diagonalized to generate the in-cage translational/vibrational eigenstates. The degeneracy of the three-dimensional harmonic oscillator energies with principal quantum number n is lifted due to the anharmonicity in the radial potential. The (2l + 1)-fold degeneracy of the angular momentum states is also weakly lifted, due to the angular dependence in the potential. We calculate the fundamental frequency to range between 96 and 110 cm-1 depending on the electronic structure method used. Error bars of the eigenstate energies were calculated from the GPR and are on the order of ∼±1.5 cm-1. Wavefunctions are also compared by considering their overlap and Hellinger distance to the one-dimensional empirical potential. As with the energies, the two abinitio methods MP2 and RPA@PBE show the best agreement. While MP2 has better agreement than RPA@PBE, due to its higher computational efficiency and comparable performance, we recommend RPA as an alternative electronic structure method of choice to MP2 for these systems.

Hyperfine and Zeeman interactions in ultracold collisions of molecular hydrogen with atomic lithium.

We present a rigorous quantum scattering study of the effects of hyperfine and Zeeman interactions on cold Li-H2 collisions in the presence of an external magnetic field using a recent abinitio potential energy surface. We find that the low-field-seeking states of H2 predominantly undergo elastic collisions: the ratio of elastic-to-inelastic cross sections exceeds 100 for collision energies below 100 mK. Furthermore, we demonstrate that most inelastic collisions conserve the space-fixed projection of the nuclear spin. We show that the anisotropic hyperfine interaction between the nuclear spin of H2 and the electron spin of Li can have a significant effect on inelastic scattering in the ultracold regime, as it mediates two processes: the electron spin relaxation in lithium and the nuclear spin-electron spin exchange. Given the predominance of elastic collisions and the propensity of inelastic collisions to retain H2 in its low-field-seeking states, our results open up the possibility of sympathetic cooling of molecular hydrogen by atomic lithium, paving the way for future exploration of ultracold collisions and high-precision spectroscopy of H2 molecules.