The barrier height of the F+H2 reaction revisited: Coupled-cluster and multireference configuration-interaction benchmark calculations

Theoretical Study on Ethylamine Dissociation Reactions Using VRC-VTST and SS-QRRK Methods

Barrierless bond dissociation reactions play an important role in fuel combustion. In this work, the pressure-dependent dissociation rate constants of ethylamine (EA) are accurately determined using variable-reaction-coordinate variational transition-state theory combined with the system-specific quantum Rice-Ramsperger-Kassel method. Before the kinetics calculations, the performances of four density functional theory methods in describing the bond dissociation of EA are evaluated against the benchmark method, FIC-MRCISD(T)+Q/cc-pVTZ, and the MN15-L/cc-pVTZ method is the best choice. By comparison of the Gibbs free energies and the rate constants for the bond dissociation reactions of EA, ethanol, and propane, the influence of functional groups on the reaction kinetics is discussed. The kinetics calculations show that the dissociation rate constants of EA are sensitive to pressure at low pressures and high temperatures, and the dominant channel is the reaction that yields C2H5 and NH2 radicals. A literature combustion model of EA is updated with our calculations, and the satisfactory agreement between the model predictions and reported ignition delay times of EA suggests the reliability of our calculations.

Spectroscopic detection of the gallium methylene (GaCH2 and GaCD2) free radical in the gas phase by laser-induced fluorescence and emission spectroscopy

GaCH2, a free radical thought to play a role in the chemical vapor deposition of gallium-containing thin films and semiconductors, has been spectroscopically detected for the first time. The radical was produced in a pulsed discharge jet using a precursor mixture of trimethylgallium vapor in high pressure argon and studied by laser-induced fluorescence and wavelength resolved emission techniques. Partially rotationally resolved spectra of the hydrogenated and deuterated species were obtained, and they exhibit the nuclear statistical weight variations and subband structure expected for a 2A2-2B1 electronic transition. The measured spectroscopic quantities have been compared to our own ab initio calculations of the ground and excited state properties. The electronic spectrum of gallium methylene is similar to the corresponding spectrum of the aluminum methylene radical, which we reported in 2022.

Multi-d-Occupancy as an Alternative Definition for the Double d-Shell Effect

Despite the prevalence of first-row transition metal-containing compounds in virtually all areas of chemistry, the accurate modeling of these systems is a known challenge for the theoretical chemistry community. Such a challenge is shown in a myriad of facets; among them are difficulties in defining ground-state multiplicities, disagreement in the results from methods considered highly accurate, and convergence problems in calculations for excited states. These problems cause a scarcity of reliable theoretical data for transition metal-containing systems. In this work, we explore the double d-shell effect that plagues and makes the application of multireference methods to this type of system difficult. We propose an alternative definition for this effect based on the mixing among d-occupancy configurations or the multi-d-occupancy character of the wave function. Moreover, we present a protocol able to include this effect in multireference calculations using an active space smaller than that usually used in the literature. A molybdenum-copper model system and its copper subsystem are used as example study cases, in particular, the molybdenum-copper charge transfer of the former and the electron affinity of the latter. We have shown that our alternative definition can be used to analyze their reference wave functions qualitatively. Based on this qualitative description, it is possible to optimize an active space without a second d-shell able to obtain relative energies accurately. Seeing the double d-shell effect through the lens of a multi-d-occupancy character, it is possible to correctly describe the wave function, improve the accuracy of the relative energies, and reduce the computational cost of multireference calculations. That way, we believe that this alternative definition has the potential to improve the modeling of first-row transition metal-containing compounds both for their ground and excited electronic structures.

Ab Initio Rovibrational Spectroscopy of the Acetylide Anion

In this work the rovibrational spectrum of the acetylide anion HCC- is investigated using high-level electronic structure methods and variational rovibrational calculations. Using a composite approach the potential energy surface and dipole surface is constructed from explicitly correlated coupled-cluster accounting for corrections due to core-valence correlation, scalar relativistic effects and higher-order excitation effects. Previous approaches for approximating the latter are critically evaluated. Employing the composite potential, accurate spectroscopic parameters determined from variational calculations are presented. In comparison to the few available reference data the present results show excellent agreement with ground state rotational constants within 0.005% of the experimental value. Intensities determined from the variational calculations suggest the bending fundamental transition ν2 around 510 cm-1 to be the best target for detection. The rather weak CD stretching fundamental ν1 in deuterated isotopologues show a second-order resonance with the (0,20,1) state and the consequences are discussed in some detail. The spectroscopic parameters and band intensities provided for a number of vibrational bands in isotopologues of the acetylide anion should facilitate future spectroscopic investigations.

A comprehensive benchmark ab initio survey of the stationary points and products of the OH· + CH3OH system

Reactions between methanol and the hydroxyl radical are of significant interest for combustion-, atmospheric-, and astrochemistry. While the two primary product channels (the formation of H₂O with either CH₃O· or ·CH₂OH) have been the subject of numerous studies, the possibility of other products has seen little attention. Here, we present a comprehensive thermochemical survey of the stationary points and plausible products of the reaction, featuring 29 geometries optimized at the UCCSD(T)-F12b/aug-cc-pVTZ level, followed by accurate composite ab initio computations for all stationary points (including ·CH₂OH dissociation and isomerization) and five product channels, with a detailed evaluation of basis set convergence and efficiency. The computations reveal that the formation of methanediol and the hydroxymethoxy radical is thermodynamically favorable and the endothermicity of formaldehyde formation is low enough to be a plausible product channel. We also observe unexpectedly large energy deviations between the partially-spin-adapted ROHF-RCCSD(T) method and ROHF-UCCSD(T) as well as between UHF-UCCSDT(Q) and ROHF-UCCSDT(Q) results.

Assessment of State-Averaged Driven Similarity Renormalization Group on Vertical Excitation Energies: Optimal Flow Parameters and Applications to Nucleobases

We present a comprehensive excited-state benchmark for the state-averaged (SA) driven similarity renormalization group (DSRG) [Li, C.; Evangelista, F. A. J. Chem. Phys. 2018, 148, 124106]. Following the QUEST database [Véril, M.; Scemama, A.; Caffarel, M.; Lipparini, F.; Boggio-Pasqua, M.; Jacquemin, D.; Loos, P.-F. Wiley Interdiscip. Rev. Comput. Mol. Sci. 2021, 11, e1517], 280 vertical transition energies of 35 medium-sized molecules are computed using the SA-DSRG derived second- and third-order perturbation theories (PT2/PT3) along with a nonperturbative approach [sq-LDSRG(2)]. Comparing to the theoretical best estimates, the optimal flow parameter is found to be 0.35 and 2.0 Eh-2 for SA-DSRG-PT2 and SA-DSRG-PT3, respectively. For SA-sq-LDSRG(2), a flow parameter of 1.5 Eh-2 provides converged equations without compromising the accuracy. We then assess the accuracy of the SA-DSRG hierarchy using these parameters. The SA-DSRG-PT2 scheme outperforms the level-shifted CASPT2 by 0.10 eV in mean absolute error (MAE), yet this accuracy is slightly inferior than that of CASPT2 with the ionization-potential-electron-affinity shift. Both SA-DSRG-PT3 and SA-sq-LDSRG(2) yield a MAE of 0.10 eV, which is comparable to that of CASPT3 (0.09 eV). Finally, we compute vertical excitation energies of several low-lying singlet states of nucleobases. The SA-sq-LDSRG(2) approach provides highly accurate results for π → π* excitations, while n → π* transitions are better described by SA-DSRG-PT3.

Rovibrational investigation of a new high-lying 0u + state of Cu2 by using two-color resonant four-wave-mixing spectroscopy

A highly excited electronic state of dicopper is observed and characterized for the first time. The [39.6]0_u ⁺-X¹Σ_g ⁺(0_g ⁺) system is measured at rotational resolution by using degenerate and two-color resonant four-wave-mixing, as well as laser induced fluorescence spectroscopy. Double-resonance experiments are performed by labeling selected rotational levels of the ground state by tuning the probe laser wavelength to transitions in the well-known (1-0) band of the B0_u ⁺-X¹Σ_g ⁺(0_g ⁺) electronic system. Spectra obtained by scans of the pump laser in the UV wavelength range were then assigned unambiguously by the stringent double-resonance selection rules. The absence of a Q-band suggests a parallel transition (ΔΩ = 0) and determines the term symbol of the state as 0_u ⁺ in Hund's case (c) notation. The equilibrium constants for ⁶³Cu₂ are T_e = 39 559.921(92) cm^-1, ω_e = 277.70(14) cm^-1, B_e = 0.104 942(66) cm^-1, and r_e = 2.2595(11) Å. These findings are supported by high-level ab initio calculations at the MRCI+Q level, which clearly identifies this state as resulting from a 4p ← 3d transition. In addition, three dark perturber states are found in the v = 1 and v = 2 vibrational levels of the new state. A deperturbation analysis characterizes the interaction and rationalizes the anomalous dips in the excitation spectrum of the [39.6]0_u ⁺-X¹Σ_g ⁺(0_g ⁺) system.

Dipole Moment Calculations Using Multiconfiguration Pair-Density Functional Theory and Hybrid Multiconfiguration Pair-Density Functional Theory

The dipole moment is the molecular property that most directly indicates molecular polarity. The accuracy of computed dipole moments depends strongly on the quality of the calculated electron density, and the breakdown of single-reference methods for strongly correlated systems can lead to poor predictions of the dipole moments in those cases. Here, we derive the analytical expression for obtaining the electric dipole moment by multiconfiguration pair-density functional theory (MC-PDFT), and we assess the accuracy of MC-PDFT for predicting dipole moments at equilibrium and nonequilibrium geometries. We show that MC-PDFT dipole moment curves have reasonable behavior even for stretched geometries, and they significantly improve upon the CASSCF results by capturing more electron correlation. The analysis of a dataset consisting of 18 first-row transition-metal diatomics and 6 main-group polyatomic molecules with a multireference character suggests that MC-PDFT and its hybrid extension (HMC-PDFT) perform comparably to CASPT2 and MRCISD+Q methods and have a mean unsigned deviation of 0.2-0.3 D with respect to the best available dipole moment reference values. We explored the dependence of the predicted dipole moments upon the choice of the on-top density functional and active space, and we recommend the tPBE and hybrid tPBE0 on-top choices for the functionals combined with the moderate correlated-participating-orbitals scheme for selecting the active space. With these choices, the mean unsigned deviations (in debyes) of the calculated equilibrium dipole moments from the best estimates are 0.77 for CASSCF, 0.29 for MC-PDFT, 0.24 for HMC-PDFT, 0.28 for CASPT2, and 0.25 for MRCISD+Q. These results are encouraging because the computational cost of MC-PDFT or HMC-PDFT is largely reduced compared to the CASPT2 and MRCISD+Q methods.

Spin-free formulation of the multireference driven similarity renormalization group: A benchmark study of first-row diatomic molecules and spin-crossover energetics

We report a spin-free formulation of the multireference (MR) driven similarity renormalization group (DSRG) based on the ensemble normal ordering of Mukherjee and Kutzelnigg [J. Chem. Phys. 107, 432 (1997)]. This ensemble averages over all microstates of a given total spin quantum number, and therefore, it is invariant with respect to SU(2) transformations. As such, all equations may be reformulated in terms of spin-free quantities and they closely resemble those of spin-adapted closed-shell coupled cluster (CC) theory. The current implementation is used to assess the accuracy of various truncated MR-DSRG methods (perturbation theory up to third order and iterative methods with single and double excitations) in computing the constants of 33 first-row diatomic molecules. The accuracy trends for these first-row diatomics are consistent with our previous benchmark on a small subset of closed-shell diatomic molecules. We then present the first MR-DSRG application on transition-metal complexes by computing the spin splittings of the [Fe(H₂O)₆]²⁺ and [Fe(NH₃)₆]²⁺ molecules. A focal point analysis (FPA) shows that third-order perturbative corrections are essential to achieve reasonably converged energetics. The FPA based on the linearized MR-DSRG theory with one- and two-body operators and up to a quintuple-ζ basis set predicts the spin splittings of [Fe(H₂O)₆]²⁺ and [Fe(NH₃)₆]²⁺ to be -35.7 and -17.1 kcal mol^-1, respectively, showing good agreement with the results of local CC theory with singles, doubles, and perturbative triples.

Electronic Structure Methods for the Description of Nonadiabatic Effects and Conical Intersections

Nonadiabatic effects are ubiquitous in photophysics and photochemistry, and therefore, many theoretical developments have been made to properly describe them. Conical intersections are central in nonadiabatic processes, as they promote efficient and ultrafast nonadiabatic transitions between electronic states. A proper theoretical description requires developments in electronic structure and specifically in methods that describe conical intersections between states and nonadiabatic coupling terms. This review focuses on the electronic structure aspects of nonadiabatic processes. We discuss the requirements of electronic structure methods to describe conical intersections and nonadiabatic couplings, how the most common excited state methods perform in describing these effects, and what the recent developments are in expanding the methodology and implementing nonadiabatic couplings.

A Full-Dimensional Potential Energy Surface and Dynamics of the Multichannel Reaction between H and HO2

In addition to its vital significance in combustion and atmospheric chemistry, the reaction between H' and HO₂ on the ground triplet state represents a prototype with multiple product channels, including H₂ + O₂, OH + OH, O + H₂O, and H + H'O₂. In this work, a full-dimensional accurate potential energy surface (PES) for the title reaction was developed to provide reliable descriptions for all dynamically relevant regions. Using this PES, we adopted the quasi-classical trajectory approach to study the corresponding reaction dynamics, including the reactivity of each product channel and the associated product branching ratio, the product energy distributions, product angular distributions, and associated microscopic mechanisms. For representing distributions of the product energies, such as product translational energy as well as product rotational and vibrational energies, both the traditional histogram and the kernel density estimation (KDE) methods were used and compared. It seems that the features of the resulting distributions in this work are very similar to each other among different methods. The KDE method is suggested for statistics, particularly for those populations with small oscillations in the histogram plot.

A neural network potential energy surface for the F + H2O ↔ HF + OH reaction and quantum dynamics study of the isotopic effect

We report here a global and full dimensional neural network potential energy surface for the F + CH₄ reaction and investigate the isotopic effect on the total reaction probabilities using the time-dependent wave packet method.

Exploring the electronic states of the hydroxyl dication OH2+: thermodynamic (meta)stability, bound-free emission spectra, and charge transfer processes

Accurate potential energy curves were constructed for a manifold of electronic states of the hydroxyl dication using a highly correlated electronic structure approach (SA-CASSCF/MRCI+Q/aug-cc-pV5Z). The existence of a bound (meta)stable ground state and bound low-lying states for OH2+ are ruled out, but do not exclude the possibility of its transient formation and dissociation along the repulsive ground state potential energy curve. Our results do not support the conclusion reported for the observation of OH2+ by electron ionization from ground state OH+. Despite the repulsive character of the low-lying states, thermodynamic stability was indeed verified for the states 2 4Π and 3 4Σ- along with a series of metastable high-lying doublet states. For the (quasi)bound states, we obtained vibrational levels, spectroscopic parameters, and dipole moment functions. Using accurate transition dipole moment functions, we also evaluated bound-free emission transition probabilities and radiative lifetimes. For transitions from v'= 0, our estimates of 92.8 ns (4Π) and 9.3 ns (4Σ-) indicate that the ones obtained by a multichannel theory of predissociating states are too short (2-60 ps). Landau-Zener cross sections averaged over the Maxwellian distribution of relative velocities, and rate coefficients for the reaction O2+ + H → O+ + H+ were obtained using the potential energy curves of the states 4Π and 4Σ- associated with the channel O2+ + H and the repulsive ones dissociating into O+ + H+ leading to good results for the rate constant thus supporting its importance to explain the distribution of O+ in astrophysical plasmas.

Theoretical investigation of the electronic structure and spectra of sulfur monoiodide cation, SI+

A manifold of singlet, triplet, and quintet electronic states of the sulfur monoiodide cation (SI+) correlating with the two lowest-lying dissociation channels is characterized theoretically at a high level of theoretical treatment (SA-CASSCF/MRCI+Q/aug-cc-pV5Z) for the first time. Potential energy curves, also including the effect of spin-orbit couplings, are constructed and the associated spectroscopic parameters and dissociation energies determined. As to the molecular polarity, we computed the dipole moment as a function of the internuclear distance and the associated vibrationally averaged dipole moments. Transition dipole moment functions were also constructed, and transition probabilities, as expressed by the Einstein coefficients for spontaneous emission, were evaluated for selected pairs of states that we identify as more easily accessible to experimental investigation. An analysis of the bonding in this system is also presented. Together with previous studies on neutral and cationic sulfur-monohalides, one has a comprehensive view of this series of molecules.

Ab initio spectroscopy of the aluminum methylene (AlCH2) free radical

Extensive ab initio investigations of the ground and electronic excited states of the AlCH₂ free radical have been carried out in order to predict the spectroscopic properties of this, as yet, undetected species. Difficulties with erratic predictions of the ground state vibrational frequencies, both in the literature and in the present work, have been traced to serious broken-symmetry instabilities in the unrestricted Hartree-Fock orbitals at the ground state equilibrium geometry. The use of restricted open-shell Hartree-Fock or complete active space self consistent field orbitals avoids these problems and leads to consistent and realistic sets of vibrational frequencies for the ground state. Using the internally contracted multireference configuration interaction method with aug-cc-pV(T+d)Z basis sets, we have calculated the geometries, energies, dipole moments, and vibrational frequencies of eight electronic states of AlCH₂ and AlCD₂. In addition, we have generated Franck-Condon simulations of the expected vibronic structure of the Ã-X̃, B̃-X̃, C̃-X̃, and C̃-Ã band systems, which will be useful in searches for the electronic spectra of the radical. We have also simulated the expected rotational structure of the 0-0 absorption bands of these transitions at modest resolution under supersonic expansion cooled conditions. Our conclusion is that if AlCH₂ can be generated in sufficient concentrations in the gas phase, it is most likely detectable through the B̃²A₂-X̃²B₁ or C̃²A₁-X̃²B₁ electronic transitions at 515 nm and 372 nm, respectively. Both band systems have vibrational and rotational signatures, even at modest resolution, that are diagnostic of the aluminum methylene free radical.

Characterizing structures, energetics, and spectra of species on the 1,3[H, C, As] potential energy surfaces: A high-level theoretical contribution

The ground and the low lying electronic states of structures on the ^1,3[H, C, As] potential energy surfaces were investigated with the highly correlated theoretical approaches CCSD(T), CCSD(T)-F12b, and CASSCF/MRCI along with the series of correlation consistent (aug-cc-pVnZ, n = D, T, Q, 5) basis sets. Energetic and spectroscopic parameters were obtained at the complete basis set limit, and the effect of core-valence correlation on these properties evaluated. Fundamental frequencies were also computed with the variational configuration interaction (VCI) approach. Heats of formation at 0 and 298.15 K were estimated for HCAs and CH, AsH, CAs, and HCAs, as well as the calculation of ionization potentials for HCAs. Comparisons of the present results with literature ones for the systems HCN/HNC, HCP/HPC highlight similarities and differences among these systems. Altogether, this investigation provides a very reliable characterization of the species on the surfaces and should guide future experimental studies on these systems.

Two-State, Three-Mode Parametrization of the Force Field of a Retinal Chromophore Model

In recent years, the potential energy surfaces of the penta-2,4-dieniminium cation have been investigated using several electronic structure methods. The resulting pool of geometrical, electronic, and energy data provides a suitable basis for the construction of a topographically correct analytical model of the molecule force field and, therefore, for a better understanding of this class of molecules, which includes the chromophore of visual pigments. In the present contribution, we report the construction of such a model for regions of the force field that drive the photochemical and thermal isomerization of the central double bound of the cation. While previous models included only two modes, it is here shown that the proposed three-mode model and corresponding set of parameters are able to reproduce the complex topographical and electronic structure features seen in electronically correlated data obtained at the XMCQDPT2//CASSCF/6-31G* level of theory.

Ab initio interaction potentials of the Ba, Ba+ complexes with Ar, Kr, and Xe in the lowest excited states

The complexes of the Ba atom and Ba⁺ cation with the rare gas atoms Ar, Kr, and Xe in the states associated with the 6s → 5d, 6p excitations are investigated by means of the multireference configuration interaction techniques. Scalar relativistic potentials are obtained by the complete basis limit extrapolation through the sequence of aug-cc-pwCVnZ basis sets with the cardinal numbers n = Q, T, 5, combined with the suitable effective core potentials and benchmarked against the coupled cluster with singles, doubles, and non-iterative triples calculations and the literature data available for selected electronic states. Spin-orbit coupling is taken into account by means of the state-interacting multireference configuration interaction calculations performed for the Breit-Pauli spin-orbit Hamiltonian. The results show weak spin-orbit coupling between the states belonging to distinct atomic multiplets. General trends in the interaction strength and long-range anisotropy along the rare gas series are discussed. Vibronic spectra of the Ba and Ba⁺ complexes in the vicinity of the ¹S → ¹P° and ²S → ²P° atomic transitions and diffusion cross sections of the Ba(¹S₀, ³D_J) atom in high-temperature rare gases are calculated. Comparison with available experimental data shows that multireference calculations tend to underestimate the interaction strength for excited complexes.

A neural network potential energy surface for the F + CH4reaction including multiple channels based on coupled cluster theory

We report here a new global and full dimensional potential energy surface (PES) for the F + CH₄reaction.