Potential Energy Surface Research Articles

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.

New Potential Energy Surface for the H + Cl2 Reaction and Quantum Dynamics Studies.

The reaction of H + Cl2 → HCl + Cl plays a crucial role in various fields. However, no previous study has investigated this reaction using accurate quantum mechanical methods. In this paper, we construct a global potential energy surface (PES) using the neural network method with more than 20,000 ab initio energies obtained by the MRCI-F12+Q method with the aug-cc-pV5Z basis and extrapolated to the complete basis set limit. The spin-orbit coupling of the Cl atom has been considered in the PES. With this new PES, product state-resolved quantum dynamics calculations for the H + Cl2 (v0 = 0, j0 = 0-2) → HCl + Cl reaction was carried out. Numerical results show that the initial rotational excitation of the Cl2 has negligible effects on the reactivity. Product state-resolved integral cross sections (ICS) and rate constants reveal that the HCl is most favorably produced in its v' = 2 vibrational state. The calculated product vibrational state-resolved and total reaction rate constants suggest that the new global PES is accurate enough, as compared with the available experimental measurements.

Exploring the photochemistry of OAlOH: Photodissociation pathways and electronic spectra.

This study was focused on the photochemistry of OAlOH and three possible pathways, which were studied with high-level multireference configuration interaction abinitio calculations. We computed cuts of the six-dimensional potential energy surfaces for the ground, the lowest singlet and triplet excited states, and probed the photodissociation mechanisms and the stabilities. The OAlOH electronic spectrum, with an energy reaching 7.15eV, contained four prominent peaks. Photodissociation to AlO, OH, and AlOH constituted a plausible mechanism within the deep-UV range (λ = 250.4nm). Our data indicated the photostability of OAlOH in the near-UV‒Vis region, so detection with laser-induced fluorescence is possible. Fluorescence and phosphorescence may occur upon excitation at 363.5nm. The roles of OAlOH in the photochemical reactions of Al-bearing molecules in the upper atmosphere and VY Canis Majoris are discussed.

Formation of Methanol via Fischer-Tropsch Catalysis by Cosmic Iron Sulphide.

Chemical reactions in the gas phase of the interstellar medium face significant challenges due to its extreme conditions (i.e., low gas densities and temperatures), necessitating the presence of dust grains to facilitate the synthesis of molecules inaccessible in the gas phase. While interstellar grains are known to enhance encounter rates and dissipate energy from exothermic reactions, their potential as chemical catalysts remain less explored. Here, we present mechanistic insights into the Fischer-Tropsch-type methanol (FTT-CH3OH) synthesis by reactivity of CO with H2 and using cosmic FeS surfaces as heterogeneous catalysts. Periodic quantum chemical calculations were employed to characterize the potential energy surface of the reactions on the (011) and (001) FeS surfaces, considering different Fe coordination environments and S vacancies. Kinetic calculations were also conducted to assess catalytic capacity and allocate reaction processes within the astrochemical framework. Findings demonstrate the feasibility of FeS-based astrocatalysis in the FTT-CH3OH synthesis. The reactions and their energetics were elucidated from a mechanistic standpoint. Kinetic analysis demonstrates the temperature dependency of the simulated processes, underscoring the compulsory need of energy sources considering the astrophysical scenario. Our results provide insights into the presence of CH3OH in diverse regions where current models struggle to explain its observational quantity.

Constructing and Evaluating Machine-Learned Interatomic Potentials for Li-Based Disordered Rocksalts.

Lithium-based disordered rocksalts (LDRs), which are an important class of positive electrode materials that can increase the energy density of current Li-ion batteries, represent a significantly complex chemical and configurational space for conventional density functional theory (DFT)-based high-throughput screening approaches. Notably, atom-centered machine-learned interatomic potentials (MLIPs) are a promising pathway to accurately model the potential energy surface of highly disordered chemical spaces, such as LDRs, where the performance of such MLIPs has not been rigorously explored yet. Here, we represent a comprehensive evaluation of the accuracy, transferability, and ease of training of five atom-centered MLIPs, including the artificial neural network potentials developed by the atomic energy network (AENET), the Gaussian approximation potential (GAP), the spectral neighbor analysis potential (SNAP) and its quadratic extension (qSNAP), and the moment tensor potential (MTP), in modeling a 11-component LDR chemical space. Specifically, we generate a DFT-calculated data set of 10,842 configurations of disordered LiTMO2 and TMO2 compositions, where TM = Sc, Ti, V, Cr, Mn, Fe, Co, Ni, and/or Cu. To provide a point-of-comparison on the performance of atom-centered MLIPs, we also trained the neural equivariant interatomic potential (NequIP) on a subset of our data. Importantly, we find AENET to be the best potential in terms of accuracy and transferability for energy predictions, while MTP is the best for atomic forces. While AENET is the fastest to train among the MLIPs considered at low number of epochs (300), the training time increases significantly as epochs increase (3300), with a corresponding reduction in training errors (∼60%). Note that AENET and GAP tend to overfit in small data sets, with the extent of overfitting reducing with larger data sets. Finally, we observe AENET to provide reasonable predictions of average Li-intercalation voltages in layered, single-TM LiTMO2 frameworks, compared to DFT (∼10% error on average). Our study should pave the way both for discovering novel disordered rocksalt electrodes and for modeling other configurationally complex systems, such as high-entropy ceramics and alloys.

New Diabatic Potential Energy Surfaces for the Li + H2 Reaction and Time-Dependent Quantum Wave Packet Studies.

New global diabatic potential energy surfaces (DPESs) for the ground (12A') and first excited (22A') states for the Li + H2 system were developed, with more than 30,000 energy points at the IC-MRCI+Q level of theory, utilizing the aug-cc-pV5Z basis set for the H atoms and the cc-pCV5Z basis set for the Li atom, fitted by a single neural network (NN) with symmetry. Product state-resolved quantum dynamics calculations of the nonadiabatic reaction Li (2P) + H2 (X 1 ∑g+, v0 = 0, j0 = 0) → LiH (X 1∑+) + H(2S) were carried out using these new DPESs and also the previous HYLC-DPESs. The numerical results suggested that our newly constructed DPESs provided an accurate description of the LiH2 system.

Comparing ANI-2x, ANI-1ccx neural networks, force field, and DFT methods for predicting conformational potential energy of organic molecules

In this study, the conformational potential energy surfaces of Amylmetacresol, Benzocaine, Dopamine, Betazole, and Betahistine molecules were scanned and analyzed using the neural network architecture ANI-2 × and ANI-1ccx, the force field method OPLS, and density functional theory with the exchange-correlation functional B3LYP and the basis set 6-31G(d). The ANI-1ccx and ANI-2 × methods demonstrated the highest accuracy in predicting torsional energy profiles, effectively capturing the minimum and maximum values of these profiles. Conformational potential energy values calculated by B3LYP and the OPLS force field method differ from those calculated by ANI-1ccx and ANI-2x, which account for non-bonded intramolecular interactions, since the B3LYP functional and OPLS force field weakly consider van der Waals and other intramolecular forces in torsional energy profiles. For a more comprehensive analysis, electronic parameters such as dipole moment, HOMO, and LUMO energies for different torsional angles were calculated at two levels of theory, B3LYP/6-31G(d) and ωB97X/6-31G(d). These calculations confirmed that ANI predictions are more accurate than density functional theory calculations with B3LYP functional and OPLS force field for determining potential energy surfaces. This research successfully addressed the challenges in determining conformational potential energy levels and shows how machine learning and deep neural networks offer a more accurate, cost-effective, and rapid alternative for predicting torsional energy profiles.

The N+ Formation Mechanism of Vibrationally Selected N2O+ Ions in the C2Σ+ State: A TPEPICO Imaging Study.

The N-NO bond fission of N2O+(C2Σ+) ions can produce two major fragment ions, NO+ or N+. In contrast to the dominant NO+ fragment ion, the N+ formation mechanism remains unclear to date. Here, dissociative photoionization of N2O via the C2Σ+ ionic state has been reinvestigated using a combined approach of threshold photoelectron-photoion coincidence (TPEPICO) velocity imaging and quantum chemical calculations. Accompanying the N+(3P) formation, the NO(X2Π) neutral fragment with low and high vi-rotational distributions was identified, based on the N+ speed and angular distributions derived from the TPEPICO images. In particular, the excitation of the symmetric stretching ν1+ mode promotes the formation of high rotational components, while the asymmetric stretching ν3+ mode shows the exact opposite effect. According to our calculated multistate potential energy surfaces, intersystem crossing from C2Σ+ to 14Π exclusively provides feasible decomposition pathways to produce the N+ fragment. In a slightly bent geometry, spin-orbit couplings between C2Σ+ and two substates of 14Π, 14A' or 14A″, play a crucial role in the N+ formation from vibrationally selected N2O+(C2Σ+) ions. The mechanism also provides new insights into the charge transfer reaction of N+ + NO → N + NO+.

Vibrational relaxation of D2O induced by collision with He: A rigid bender close coupling study.

The vibrational relaxation of the first excited bending state ofD2O induced by collision with He is studied at the close coupling level and using the Rigid Bender approximation. A new 4D potential energy surface is calculated and reported for this system. It is then used to determine the low-lying bound states of the D2O-He van der Waals complex and to perform scattering calculations. Colli- sion rates are determined for pure rotational transitions as well as for rovibrational transitions within the first excited bending state. The results are compared with those obtained for the collision of D2O with other noble gases such as Ne and Ar. We also analyse the differences observed with respect to the H2O + He collisions and compare our results with experiment.

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.

Are Redox-active Centers Bridged by Saturated Flexible Linkers Systematically Electrochemically Independent?

The extent to which electrophores covalently bridged by a saturated linker are electrochemically independent was investigated considering the charge/spin duality of the electron and functionality of the electrophore as a spin carrier upon reduction. By combining computational modeling with electrochemical experiments, we investigated the mechanism by which tethered electrophores react together within 4,4'-oligo[n]methylene-bipyridinium assemblies (with n = 2 to 5). We show that native dicationic electrophores (redox state Z = +2) are folded prior to electron injection into the system, allowing the emergence of supra-molecular orbitals (supra-MOs) likely to support the process of the reductive s bond formation giving cyclomers. Indeed, for Z = +2, London Dispersion (LD) forces contribute to flatten the potential energy surface such that all-trans and folded conformers are approximately isoenergetic. Then, upon one-electron injection, for radical cations (Z = +1), LD forces significantly stabilize the folded conformers, except for the ethylene derivative deprived of supra-MOs. For radical cations equipped with supra-MOs, the unpaired electron is delocalized over both heterocycles through space. Cyclomer completion (Z = 0) upon the second electron transfer occurs according to the inversion of redox potentials. This mechanism explains why intramolecular reactivity is favored and why pyridinium electrophores are not independent.

Unifying O(3) equivariant neural networks design with tensor-network formalism

Many learning tasks, including learning potential energy surfaces from ab initio calculations, involve global spatial symmetries and permutational symmetry between atoms or general particles. Equivariant graph neural networks are a standard approach to such problems, with one of the most successful methods employing tensor products between various tensors that transform under the spatial group. However, as the number of different tensors and the complexity of relationships between them increase, maintaining parsimony and equivariance becomes increasingly challenging. In this paper, we propose using fusion diagrams, a technique widely employed in simulating SU(2)-symmetric quantum many-body problems, to design new spatial equivariant components for neural networks. This results in a diagrammatic approach to constructing novel neural network architectures. When applied to particles within a given local neighborhood, the resulting components, which we term ‘fusion blocks,’ serve as universal approximators of any continuous equivariant function defined on the neighborhood. We incorporate a fusion block into pre-existing equivariant architectures (Cormorant and MACE), leading to improved performance with fewer parameters on a range of challenging chemical problems. Furthermore, we apply group-equivariant neural networks to study non-adiabatic molecular dynamics of stilbene cis-trans isomerization. Our approach, which combines tensor networks with equivariant neural networks, suggests a potentially fruitful direction for designing more expressive equivariant neural networks.

Are Redox‐active Centers Bridged by Saturated Flexible Linkers Systematically Electrochemically Independent?

The extent to which electrophores covalently bridged by a saturated linker are electrochemically independent was investigated considering the charge/spin duality of the electron and functionality of the electrophore as a spin carrier upon reduction. By combining computational modeling with electrochemical experiments, we investigated the mechanism by which tethered electrophores react together within 4,4'‐oligo[n]methylene‐bipyridinium assemblies (with n = 2 to 5). We show that native dicationic electrophores (redox state Z = +2) are folded prior to electron injection into the system, allowing the emergence of supra‐molecular orbitals (supra‐MOs) likely to support the process of the reductive s bond formation giving cyclomers. Indeed, for Z = +2, London Dispersion (LD) forces contribute to flatten the potential energy surface such that all‐trans and folded conformers are approximately isoenergetic. Then, upon one‐electron injection, for radical cations (Z = +1), LD forces significantly stabilize the folded conformers, except for the ethylene derivative deprived of supra‐MOs. For radical cations equipped with supra‐MOs, the unpaired electron is delocalized over both heterocycles through space. Cyclomer completion (Z = 0) upon the second electron transfer occurs according to the inversion of redox potentials. This mechanism explains why intramolecular reactivity is favored and why pyridinium electrophores are not independent.

The Near-Sightedness of Many-Body Interactions in Anharmonic Vibrational Couplings.

Couplings between vibrational motions are driven by electronic interactions, and these couplings carry special significance in vibrational energy transfer, multidimensional spectroscopy experiments, and simulations of vibrational spectra. In this investigation, the many-body contributions to these couplings are analyzed computationally in the context of clathrate-like alkali metal cation hydrates, including Cs+(H2O)20, Rb+(H2O)20, and K+(H2O)20, using both analytic and quantum-chemistry potential energy surfaces. Although the harmonic spectra and one-dimensional anharmonic spectra depend strongly on these many-body interactions, the mode-pair couplings were, perhaps surprisingly, found to be dominated by one-body effects, even in cases of couplings to low-frequency modes that involved the motion of multiple water molecules. The origin of this effect was traced mainly to geometric distortion within water monomers and cancellation of many-body effects in differential couplings, and the effect was also shown to be agnostic to the identity of the ion. These outcomes provide new understanding of vibrational couplings and suggest the possibility of improved computational methods for the simulation of infrared and Raman spectra.

Competition between O-H and S-H Intermolecular Interactions in Conformationally Complex Systems: The 2-Phenylethanethiol and 2-Phenylethanol Dimers.

Noncovalent interactions involving sulfur centers play a relevant role in biological and chemical environments. Yet, detailed molecular descriptions are scarce and limited to very simple model systems. Here we explore the formation of the elusive S-H···S hydrogen bond and the competition between S-H···O and O-H···S interactions in pure and mixed dimers of the conformationally flexible molecules 2-phenylethanethiol (PET) and 2-phenylethanol (PEAL), using the isolated and size-controlled environment of a jet expansion. The structure of both PET-PET and PET-PEAL dimers was unraveled through a comprehensive methodology that combined rotationally resolved microwave spectroscopy, mass-resolved isomer-specific infrared laser spectroscopy, and quantum chemical calculations. This synergic experimental-computational approach offered unique insights into the potential energy surface, conformational equilibria, molecular structure, and intermolecular interactions of the dimers. The results show a preferential order for establishing hydrogen bonds following the sequence S-H···S < S-H···O ≲ O-H···S < O-H···O, despite the hydrogen bond only accounting for a fraction of the total interaction energy.

Resolving Photoinduced Femtosecond Three-Dimensional Solute-Solvent Dynamics through Surface Hopping Simulations.

Photoinduced dynamics in solution is governed by mutual solute-solvent interactions, which give rise to phenomena like solvatochromism, the Stokes shift, dual fluorescence, or charge transfer. Understanding these phenomena requires simulating the solute's photoinduced dynamics and simultaneously resolving the three-dimensional solvent distribution dynamics. If using trajectory surface hopping (TSH) to this aim, thousands of trajectories are required to adequately sample the time-dependent three-dimensional solvent distribution functions, and thus resolve the solvent dynamics with sub-Ångstrom and femtosecond accuracy and sufficiently low noise levels. Unfortunately, simulating thousands of trajectories with TSH in the framework of hybrid quantum mechanical/molecular mechanical (QM/MM) can be prohibitively expensive when employing ab initio electronic structure methods. To tackle this challenge, we recently introduced a computationally efficient approach that combines efficient linear vibronic coupling models with molecular mechanics (LVC/MM) via electrostatic embedding [Polonius et al., JCTC 2023, 19, 7171-7186]. This method provides solvent-embedded, nonadiabatically coupled potential energy surfaces while scaling similarly to MM force fields. Here, we employ TSH with LVC/MM to unravel the photoinduced dynamics of two small thiocarbonyl compounds solvated in water. We describe how to estimate the number of trajectories required to produce nearly noise-free three-dimensional solvent distribution functions and present an analysis based on approximately 10,000 trajectories propagated for 3 ps. In the electronic ground state, both molecules exhibit in-plane hydrogen bonds to the sulfur atom. Shortly after excitation, these bonds are broken and reform perpendicular to the molecular plane on timescales that differ by an order of magnitude due to steric effects. We also show that the solvent relaxation dynamics is coupled to the electronic dynamics, including intersystem crossing. These findings are relevant to advance the understanding of the coupled solute-solvent dynamics of solvated photoexcited molecules, e.g., biologically relevant thio-nucleobases.

Bimolecular collision outcomes on multidimensional potential energy surfaces: infrared spectroscopy and activation of NO-alkane collision complexes.

In bimolecular collisions between open-shell radicals and increasingly-larger alkanes, the relative impact configurations open the possibility of reactive and nonreactive outcomes that are isomer specific. To model the interaction potential between molecular scattering partners, observables are needed from experiments that can quantify both the initial molecular orientations and internal energies on multidimensional potential energy surfaces. Recent work by our group demonstrated that upon infrared (IR) excitation, the dynamics of the nitric oxide-methane collision complex (NO-CH4) are dependent on the initial monomer geometries, as small changes in configuration substantially affect the energies, electronic couplings, and predissociation pathways due to the Jahn-Teller effect. This study focuses on the isomer-specific scattering mechanisms between NO and ethane (C2H6), encoded in the spectroscopic and dynamical signatures of the NO-C2H6 collision complex. IR action spectroscopy with 1 + 1 resonance-enhanced multiphoton ionization of NO products was employed to characterize the fundamental CH stretch transitions of NO-C2H6, as well as to initiate the nonreactive decay mechanisms of the complex. Furthermore, velocity map imaging (VMI) was utilized to explore the dynamics prior to and following IR excitation of NO-C2H6, imprinted on the NO photoproducts. This work compares the dynamics from NO-C2H6 and NO-CH4 vibrational predissociation, in which substantial differences are observed in the energy exchange mechanisms during the evolution of the collision complexes to products.

Computational exploration of the copper(I)-catalyzed conversion of hydrazones to dihalogenated vinyldiazene derivatives.

This computational study explores the copper (I) chloride catalyzed synthesis of (E)-1-(2,2-dichloro-1-phenylvinyl)-2-phenyldiazene (2Cl-VD) from readily available hydrazone derivative and carbon tetrachloride (CCl4). 2Cl-VD has been extensively utilized to synthesize variety of heterocyclic organic compounds in mild conditions. The present computational investigations primarily focus on understanding the role of copper (I) and N1,N1,N2,N2-tetramethylethane-1,2-diamine (TMEDA) in this reaction, TMEDA often being considered a proton scavenger by experimentalists. Considering TMEDA as a ligand significantly alters the energy barrier. In fact, it is only 8.3 kcal/mol higher compared to the ligand-free (LF) route for the removal of a chlorine atom to form the radical ·CCl3 but the following steps are almost barrierless. This intermediate then participates in attacking the electrophilic carbon in the hydrazone. Crucially, the study reveals that the overall potential energy surface is thermodynamically favorable, and the theoretical turnover frequency (TOF) value is higher in the case of Cu(I)-TMEDA complex catalyzed pathway.

Six-dimensional quantum dynamics of an Eley-Rideal reaction between gaseous and adsorbed hydrogen atoms on Cu(111).

In the form of direct abstraction of a surface adsorbate by a gaseous projectile, the Eley-Rideal (ER) reaction at the gas-surface interface manifests interesting dynamics. Unfortunately, high-dimensional quantum dynamical (QD) studies for ER reactions remain very challenging, which demands a large configuration space and the coordinate transformation of wavefunctions. Here, we report the first six-dimensional (6D) fully coupled quantum scattering method for studying the ER reaction between gas phase H(D) atoms and adsorbed D(H) atoms on a rigid Cu(111) surface. Reaction probabilities and product rovibrational state distributions obtained by this 6D model are found to be quite different from those obtained by reduced-dimensional QD models, demonstrating the high-dimensional nature of the ER reaction. Using two distinct potential energy surfaces (PESs), we further discuss the influence of the PES on the calculated product vibrational and rotational state distributions, in comparison with experimental results. The lateral corrugation and the exothermicity of the PES are found to play a critical role in controlling the energy disposal in the ER reaction.

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.