Potential Energy Surface Research Articles

Revisiting the Structure and Bonding in Li5H6- and the Exploration of Reactivity: Planar Pentacoordinate Hydrogen.

Recently, Guha and co-workers (Sarmah, K.; Kalita, A.; Purkayastha, S.; Guha, A. K. Pushing The Extreme of Multicentre Bonding: Planar Pentacoordinate Hydride. Angew. Chem. Int. Ed. 2024, e202318741) reported a highly intriguing bonding motif: planar pentacoordinate hydrogen (ppH) in Li5H6-, featuring C2v symmetry in the singlet state with two distinct H-Li (center-ring) bond distances. We herein revisited the potential energy surface of Li5H6- by using a target-oriented genetic algorithm. Our investigation revealed that the lowest-energy structure of Li5H6- exhibits a ppH configuration with very high D5h symmetry and a 1A1' electronic state. We did not find any electronic effect like Jahn-Teller distortion that could be responsible for lowering its symmetry. Moreover, our calculations demonstrated significant differences in the relative energies of other low-lying isomers. An energetically very competitive planar tetracoordinate hydrogen (ptH) isomer is also located, but it corresponds to a very shallow minimum on the potential energy surface depending on the used level of theory. Chemical bonding analyses, including AdNDP and EDA-NOCV, uncover that the optimal Lewis structure for Li5H6- involves H- ions stabilized by the Li5H5 crown. Surprisingly, despite the dominance of electrostatic interactions, the contribution from covalent bonding is also significant between ppH and the Li5H5 moiety, derived from H-(1s) → Li5H5 σ donation. Magnetically induced current density analysis revealed that due to minimal orbital overlap and the highly polar nature of the H-Li covalent interaction, the ppH exhibits local diatropic ring currents around the H centers, which fails to result in a global aromatic ring current. The coordination of Li5H6- with Lewis acids, BH3 and BMe3, instantly converts the ppH configuration to (quasi) ptH. These Lewis acid-bound ptH complexes show high electronic stability and high thermochemical stability against dissociation and, therefore, will be ideal candidates for the experimental realization.

To Pair or not to Pair? Machine-Learned Explicitly-Correlated Electronic Structure for NaCl in Water.

The extent of ion pairing in solution is an important phenomenon to rationalize transport and thermodynamic properties of electrolytes. A fundamental measure of this pairing is the potential of mean force (PMF) between solvated ions. The relative stabilities of the paired and solvent shared states in the PMF and the barrier between them are highly sensitive to the underlying potential energy surface. However, direct application of accurate electronic structure methods is challenging, since long simulations are required. We develop wave function based machine learning potentials with the random phase approximation (RPA) and second order Møller-Plesset (MP2) perturbation theory for the prototypical system of Na and Cl ions in water. We show both methods in agreement, predicting the paired and solvent shared states to have similar energies (within 0.2 kcal/mol). We also provide the same benchmarks for different DFT functionals as well as insight into the PMF based on simple analyses of the interactions in the system.

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.

The anion photoelectron spectrum and diabatization of tetrazolyl.

The potential energy surface of tetrazolyl [cyclic (N4CH)] has a conical intersection seam between the two lowest-energy electronic states near the ground state minimum geometry. This work treats that molecule. The potential energy surfaces used in this study are based on a least-squares fitting procedure that includes abinitio energies, energy gradients, and derivative couplings described using polynomials up to fourth-order and abinitio data obtained from multireference configuration interaction wave functions. A five-electronic-state description was generated with a root mean square absolute energy error of 9.6cm-1, compared to 326.8cm-1 when only second-order terms were used. The time-independent multimode vibronic coupling in the KDC approximation was used to simulate and analyze the anion ultraviolet photoelectron spectrum of tetrazolide.

Using a pruned basis and a sparse collocation grid with more points than basis functions to do efficient and accurate MCTDH calculations with general potential energy surfaces.

We propose a new collocation multi-configuration time-dependent Hartree (MCTDH) method. It reduces point-set error by using more points than basis functions. Collocation makes it possible to use MCTDH with a general potential energy surface without computing any integrals. The collocation points are associated with a basis larger than the basis used to represent wavefunctions. Both bases are obtained from a direct product basis built from single-particle functions by imposing a pruning condition. The collocation points are those on a sparse grid. Heretofore, collocation MCTDH calculations with more points than basis functions have only been possible if both the collocation grid and the basis set are direct products. In this paper, we exploit a new pseudo-inverse to use both more points than basis functions and a pruned basis and grid. We demonstrate that, for a calculation of the lowest 50 vibrational states (energy levels and wavefunctions) of CH2NH, errors can be reduced by two orders of magnitude by increasing the number of points, without increasing the basis size. This is true also when unrefined time-independent points are used.

A fast and smooth one-electron approach for investigating charge transfer states and D1-D0 crossings for systems with odd numbers of electrons.

We propose an efficient version of ensemble Hartree-Fock/density functional theory to calculate a set of two charge-transfer states for systems with odd-numbers of electrons. The approach does require definitions of donor/acceptor fragments; however, the approach is not very sensitive to such definitions-even in the limit of very strong electronic coupling. The key ansatz is that, by mandating that the vector space spanned by the active orbitals projects equally onto the donor and acceptor fragments, such a constraint eliminates all intra-molecular local excitations and makes it far easier to generate potential energy surfaces that are smooth over a wide region of configuration space. The method is fast, working with only two electron configurations, and should be useful for abinitio non-adiabatic dynamics in the near future.

Conformational Stability and Vibrational Relaxation as a Function of Halogen Substitution in Benzene sulfonamides: A Theoretical Study

AbstractStudying the vibrational energy flow within molecules can help us better understand their reactivity and facilitate the design of novel functional groups for potential drug candidates. This study examined the conformation of benzene sulfonamide in solvents and vacuum. Potential energy surfaces were computed for the two possible conformations, showing that the conformers’ stability varied with the medium. Additionally, the vibrational energy redistribution of halogen‐substituted benzene sulfonamide conformers was studied in solvents and vacuums. The lifetimes of the halogen stretching mode C−X (where X=F, Cl, or Br) in the substituted benzene sulfonamides were calculated and compared. The solvent and vacuum media lifetimes followed the order: C−Cl>C−Br>C−F. However, the degree of delocalization of vibrational energy was found to vary between them.

Regularized relativistic corrections for polyelectronic and polyatomic systems with explicitly correlated Gaussians.

Drachmann's regularization approach is implemented for floating explicitly correlated Gaussians (fECGs) and molecular systems. Earlier applications of drachmannized relativistic corrections for molecular systems were hindered due to the unknown analytic matrix elements of 1/rix1/rjy-type operators with fECGs. In the present work, one of the 1/r factors is approximated by a linear combination of Gaussians, which results in calculable integrals. The numerical approach is found to be precise and robust over a range of molecular systems and nuclear configurations, and thus, it opens the route toward an automated evaluation of high-precision relativistic corrections over potential energy surfaces of polyatomic systems. Furthermore, the newly developed integration approach makes it possible to construct the matrix representation of the square of the electronic Hamiltonian relevant for energy lower-bound as well as time-dependent computations of molecular systems with a flexible and high-precision fECG basis representation.

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.

The Hydrogen Evolution Activity of BaZrS3, BaTiS3, and BaVS3 Chalcogenide Perovskites.

Chalcogenide perovskites are a class of materials with electronic and optoelectronic properties desirable for solar cells, infrared optics, and computing. The oxide counterparts of these chalcogenides have been studied extensively for their electrocatalytic and photoelectrochemical properties. As chalcogenide perovskites are more covalent, conductive, and stable, we hypothesize that they are more viable as electrocatalysts than oxide perovskites. The goal of this synthetic, experimental, and computational study is to examine the hydrogen evolution reaction (HER) activity of three Barium-based chalcogenides in perovskite and related structures: BaZrS3, BaTiS3, and BaVS3. Potential energy surfaces for hydrogen adsorption on surfaces of these materials are calculated using density functional theory and the computational hydrogen electrode model is used to contrast overpotentials with experiment. Although both experiments and computations agree that BaVS3 is the most active of the three materials, high overpotentials of these materials make them less viable than platinum for HER. Our work establishes a framework for future studies in the chemical and electrochemical properties of chalcogenide perovskites.

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+.

Multistate Dynamics and Kinetics of CO2 Activation by Ta+ in the Gas Phase: Insights into Single-Atom Catalysis.

The activation of carbon dioxide (CO2) by a transition-metal cation in the gas phase is a unique model system for understanding single-atom catalysis. The mechanism of such reactions is often attributed to a "two-state reactivity" model in which the high-energy barrier of a spin state correlating with ground-state reactants is avoided by intersystem crossing (ISC) to a different spin state with a lower barrier. However, such a "spin-forbidden" mechanism, along with the corresponding dynamics, has seldom been rigorously examined theoretically, due to the lack of global potential energy surfaces (PESs). In this work, we report full-dimensional PESs of the lowest-lying quintet, triplet, and singlet states of the TaCO2+ system, machine-learned from first-principles data. These PESs and the corresponding spin-orbit couplings enable us to provide an extensive theoretical characterization of the dynamics and kinetics of the reaction between the tantalum cation (Ta+) and CO2, which have recently been investigated experimentally at high collision energies using crossed beams and velocity map imaging, as well as at thermal energies using a selected-ion flow tube apparatus. The multistate quasi-classical trajectory simulations with surface hopping reproduce most of the measured product translational and angular distributions, shedding valuable light on the nonadiabatic reaction dynamics. The calculated rate coefficients from 200 to 600 K are also in good agreement with the latest experimental measurements. More importantly, these calculations revealed that the reaction is controlled by intersystem crossing, rather than potential barriers.

Breaking down walls: Continuous potential models for internal motions in NMR spin relaxation

Simple physical models for restricted diffusion in a potential, which provide important insights for NMR spin relaxation, usually are based on free diffusion within rigid boundaries or diffusion in relatively simple continuous potential energy surfaces. The diffusion-in-a-cone model is an example of the former and diffusion in an N-fold cosine potential is an example of the latter. The present work models restricted diffusion for arbitrary potential energy functions on the surface of a cone or a sphere, by expanding the potentials in Fourier or spherical harmonic series, respectively. The results exhibit simple relationships between generalized order parameters and effective correlation times, critical for analysis of experimental spin relaxation data, and illustrate the transition from diffusive-like to jump-like behavior in multi-well potentials.

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.