Potential energy surfaces and dynamic properties via ab initio composite and density functional approaches

Vibrational spectroscopy enables critical insight into the structural and dynamic properties of molecules. Presently, the majority of theoretical approaches to spectroscopy employ wavefunction-based ab initio or density functional methods that rely on the harmonic approximation. This approximation breaks down for large molecules with strongly anharmonic bonds or for molecules with large internuclear separations. An alternative to these methods involves generating molecular anharmonic potential energy surfaces (potentials) and using them to extrapolate the vibrational frequencies. This study examines the efficacy of density functional theory (DFT) and the correlation consistent Composite Approach (ccCA) in generating anharmonic frequencies from potentials of small main group molecules. Vibrational self-consistent field Theory (VSCF) and post-VSCF methods were used to calculate the fundamental frequencies of these molecules from their potentials. Functional choice, basis set selection, and mode-coupling are also examined as factors in influencing accuracy. The absolute deviations for the calculated frequencies using potentials at the ccCA level of theory were lower than the potentials at the DFT level. With DFT resulting in bending modes that are better described than those of ccCA, a multilevel DFT:ccCA approach where DFT potentials are used for single vibrational mode potentials and ccCA is used for vibrational mode-mode couplings can be utilized for larger polyatomic systems. The frequencies obtained with this multilevel approach using VCIPSI-PT2 were closer to experimental frequencies than the scaled harmonic frequencies, indicating the success of utilizing post-VSCF methods to generate more accurate representations of computed infrared spectra.

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.

Ab initio molecular dynamics study of intersystem crossing dynamics for MH2 (M = Si, Ge, Sn, Pb) on spin-pure and spin-mixed potential energy surfaces

Recently, surface-hopping ab initio molecular dynamics (SH-AIMD) simulations have come to be used to discuss the mechanisms and dynamics of excited-state chemical reactions, including internal conversion and intersystem crossing. In dynamics simulations involving intersystem crossing, there are two potential energy surfaces (PESs) governing the motion of nuclei: PES in a spin-pure state and PES in a spin-mixed state. The former gives wrong results for molecular systems with large spin-orbit coupling (SOC), while the latter requires a potential gradient that includes a change in SOC at each point, making the computational cost very high. In this study, we systematically investigate the extent to which the magnitude of SOC affects the results of the spin-pure state-based dynamics simulations for the hydride MH2 (M = Si, Ge, Sn, Pb) by performing SH-AIMD simulations based on spin-pure and spin-mixed states. It is clearly shown that spin-mixed state PESs are indispensable for the dynamics simulation of intersystem crossing in systems containing elements Sn and Pb from the fifth period onward. Furthermore, in addition to the widely used Tully's fewest switches (TFS) algorithm, the Zhu-Nakamura (ZN) global switching algorithm, which is computationally less expensive, is applied to SH for comparison. The results from TFS- and ZN-SH-AIMD methods are in qualitative agreement, suggesting that the less expensive ZN-SH-AIMD can be successfully utilized to investigate the dynamics of photochemical reactions based on quantum chemical calculations.

A ReaxFF Potential for Modeling Organic Matter Degradation with Oxybromine Oxidants

Oxidation of organic matter with oxybromine oxidants is ushering in a new era of enhanced hydrocarbon recovery. While these potent reagents are being tested in laboratory and field experiments, there is a pressing demand to delineate the molecular processes governing oxidation reactions at geological depth. Here, we parameterize a ReaxFF potential to model the oxidative decompositions of aliphatic and aromatic hydrocarbons in the presence of water-NaBr solutions that contain oxybromine (BrOn)- oxidizers. Our parameterization results in a reliable empirical bond-order potential that accurately calculates bond energies, exhibiting an RMSE of ∼1.18 eV, corresponding to 1.36 % average error. Reproducing bond dissociation and binding energies from Density Functional Theory (DFT), our parameterization proves transferable to aqueous environments. This H/C/O/Na/Br ReaxFF potential accurately reproduces the oxidation pathways of small hydrocarbons with oxybromine oxidizers. This force field captures proton and oxygen transfer, C-C bond tautomerization, and cleavage, leading to ring-opening and chain fragmentation. Molecular dynamic simulations demonstrate the oxidative degradation of aromatic and aliphatic kerogen-like moieties in bulk solutions. We envision that such reactive force fields will be useful to understand better the oxidation reactions of organic matter formed in geological reservoirs for enhanced shale gas recovery and improved carbon dioxide treatments.

A Q-learning method based on coarse-to-fine potential energy surface for locating transition state and reaction pathway

Transition state (TS) on the potential energy surface (PES) plays a key role in determining the kinetics and thermodynamics of chemical reactions. Inspired by the fact that the dynamics of complex systems are always driven by rare but significant transition events, we herein propose a TS search method in accordance with the Q-learning algorithm. Appropriate reward functions are set for a given PES to optimize the reaction pathway through continuous trial and error, and then the TS can be obtained from the optimized reaction pathway. The validity of this Q-learning method with reasonable settings of Q-value table including actions, states, learning rate, greedy rate, discount rate, and so on, is exemplified in 2 two-dimensional potential functions. In the applications of the Q-learning method to two chemical reactions, it is demonstrated that the Q-learning method can predict consistent TS and reaction pathway with those by ab initio calculations. Notably, the PES must be well prepared before using the Q-learning method, and a coarse-to-fine PES scanning scheme is thus introduced to save the computational time while maintaining the accuracy of the Q-learning prediction. This work offers a simple and reliable Q-learning method to search for all possible TS and reaction pathway of a chemical reaction, which may be a new option for effectively exploring the PES in an extensive search manner.

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.

Advances in the Dynamics of Adsorbate Diffusion on Metal Surfaces: Focus on Hydrogen and Oxygen

Adsorbates on metal surfaces are typically formed from the dissociative chemisorption of molecules occurring at gas-solid interfaces. These adsorbed species exhibit unique diffusion behaviors on metal surfaces, which are influenced by their translational energy. They play crucial roles in various fields, including heterogeneous catalysis and corrosion. This review examines recent theoretical advancements in understanding the diffusion dynamics of adsorbates on metal surfaces, with a specific emphasis on hydrogen and oxygen atoms. The diffusion processes of adsorbates on metal surfaces involve two energy transfer mechanisms: surface phonons and electron-hole pair excitations. This review also surveys new theoretical methods, including the characterization of the electron-hole pair excitation within electronic friction models, the acceleration of quantum chemistry calculations through machine learning, and the treatment of atomic nuclear motion from both quantum mechanical and classical perspectives. Furthermore, this review offers valuable insights into how energy transfer, nuclear quantum effects, supercell sizes, and the topography of potential energy surfaces impact the diffusion behavior of hydrogen and oxygen species on metal surfaces. Lastly, some preliminary research proposals are presented.

Probing the dynamics and bottleneck of the key atmospheric SO2 oxidation reaction by the hydroxyl radical

SO2 (Sulfur dioxide) is the major precursor to the production of sulfuric acid (H2SO4), contributing to acid rain and atmospheric aerosols. Sulfuric acid formed from SO2 generates light-reflecting sulfate aerosol particles in the atmosphere. This property has prompted recent geoengineering proposals to inject sulfuric acid or its precursors into the Earth's atmosphere to increase the planetary albedo to counteract global warming. SO2 oxidation in the atmosphere by the hydroxyl radical HO to form HOSO2 is a key rate-limiting step in the mechanism for forming acid rain. However, the dynamics of the HO + SO2 → HOSO2 reaction and its slow rate in the atmosphere are poorly understood to date. Herein, we use photoelectron spectroscopy of cryogenically cooled HOSO2- anion to access the neutral HOSO2 radical near the transition state of the HO + SO2 reaction. Spectroscopic and dynamic calculations are conducted on the first ab initio-based full-dimensional potential energy surface to interpret the photoelectron spectra of HOSO2- and to probe the dynamics of the HO + SO2 reaction. In addition to the finding of a unique pre-reaction complex (HO⋯SO2) directly connected to the transition state, dynamic calculations reveal that the accessible phase space for the HO + SO2 → HOSO2 reaction is extremely narrow, forming a key reaction bottleneck and slowing the reaction rate in the atmosphere, despite the low reaction barrier. This study underlines the importance of understanding the full multidimensional potential energy surface to elucidate the dynamics of complex bimolecular reactions involving polyatomic reactants.

Effects of Chalcogen Atoms on Excited-State Double-Proton Transfer Behavior for 3,6-bis(4,5-Dihydroxyoxazo-2-yl)benzene-1,2-diol Derivatives: A Computational Investigation

The impact of the chalcogen atomic electronegativity (O, S, and Se atoms) of new organic molecules on excited-state dynamical reactions is self-evident. Inspired by this kind of distinguished photochemical characteristic, in this work, we performed a computational investigation of chalcogen-substituted 3,6-bis(4,5-dihydroxyoxazo-2-yl)benzene-1,2-diol (BDYBD) derivatives (i.e., BDYBD-O, BDYBD-S, and BDYBD-Se). In this paper, we pay close attention to characteristic BDYBD derivatives that contain intramolecular double hydrogen bonds (O1-H2···N3 and O4-H5···N6). The main goal of this study was to explore how changes in atomic electronegativity affect the way hydrogen bonds interact and how excited molecules affect transfer protons. We go into further detail in the main text of the paper. By fixing our attention to geometrical variations and infrared (IR) vibrational spectra between the S0 and S1 states, exploring hydrogen bonding behaviors using the core-valence bifurcation (CVB) index, and simulating hydrogen bonding energy (EHB) via the atom in molecule (AIM) method, we clarified the photo-induced strengthened dual hydrogen bonding interactions that facilitate the excited-state dual-proton transfer (ESDPT) behavior of BDYBD derivatives. The reorganization of charge stemming from photoexcitation further verifies the tendencies of ESDPT reactions. We relied on constructing potential energy surfaces (PESs) by adopting a restrictive optimization approach, and herein, we finally clarify the gradual ESDPT mechanism of BDYBD derivatives. Particularly, we confirm that the variation in chalcogen atomic electronegativity has a regulatory effect on the ESDPT behavior of BDYBD derivatives; that is, the lower the atomic electronegativity, the more favorable it is for the gradual ESDPT reaction.

Evolution Mechanism of Solid-Phase Catalysts During Catalytic Growth of Single-Walled Carbon Nanotubes

Using solid nanoparticles (NPs) as catalysts is the most effective method to achieve catalytic growth of single-walled carbon nanotubes (SWCNTs) with ultrapure chirality. Until now, SWCNTs with a suitable chirality purity have not been prepared in experiments. That is, the evolution of solid NPs during the catalytic growth of SWCNTs is in contradiction with the original concept of a changeless structure. Hence, in this work, the evolution mechanism of solid cobalt NPs during the nucleation process of SWCNTs is analyzed through molecular dynamics. Similar to the experimental observations, the results show that a drastic structural fluctuation of the NPs occurs during the nucleation of SWCNTs. This structural fluctuation is caused by the fact that the elastic strain energy and surface energy of the NPs can be tuned when a carbon gradient exists between the subsurface and interior of the NP. Furthermore, such a carbon gradient can be reduced by changing the carbon feeding rate. This work not only reveals the evolution mechanism of solid catalysts during the nucleation of SWCNTs but also provides prospects for realizing solid catalysts with a changeless structure by tuning the experimental parameters.

Multiple Dissociation Pathways in HNCO Decomposition Governed by Potential Energy Surface Topography

The exquisite features of molecular photochemistry are key to any complete understanding of the chemical processes governed by potential energy surfaces (PESs). It is well established that multiple dissociation pathways relate to nonadiabatic transitions between multiple coupled PESs. However, little detail is known about how the single PES determines reaction outcomes. Here we perform detailed experiments on HNCO photodissociation, acquiring the state-specific correlations of the NH (a1Δ) and CO (X1Σ+) products. The experiments reveal a trimodal CO rotational distribution. Dynamics simulations based on a full-dimensional machine-learning-based PES of HNCO unveil three dissociation pathways exclusively occurring on the S1 excited electronic state. One pathway, following the minimum energy path (MEP) via the transition state, contributes to mild rotational excitation in CO, while the other two pathways deviating substantially from the MEP account for relatively cold and hot CO rotational state populations. These peculiar dynamics are unambiguously governed by the S1 state PES topography, i.e., a narrow acceptance cone in the vicinity of the transition state region. The dynamical picture shown in this work will serve as a textbook example illustrating the importance of the PES topography in molecular photochemistry.

Alternating Stereospecificity upon Central-Atom Change: Dynamics of the F- +PH2 Cl SN 2 Reaction Compared to its C- and N-Centered Analogues

Central-atom effects on bimolecular nucleophilic substitution (SN 2) reactions are well-known in chemistry, however, the atomic-level SN 2 dynamics at phosphorous (P) centers has never been studied. We investigate the dynamics of the F- +PH2 Cl reaction with the quasi-classical trajectory method on a novel full-dimensional analytical potential energy surface fitted on high-level ab initio data. Our computations reveal intermediate dynamics compared to the F- +CH3 Cl and the F- +NH2 Cl SN 2 reactions: phosphorus as central atom leads to a more indirect SN 2 reaction with extensive complex-formation with respect to the carbon-centered one, however, the title reaction is more direct than its N-centered pair. Stereospecificity, characteristic at C-center, does not appear here either, due to the submerged front-side-attack retention path and the repeated entrance-channel inversional motion, whereas the multi-inversion mechanism discovered at nitrogen center is also undermined by the deep Walden-well. At low collision energies, 6 % of the PH2 F products form with retained configuration, mostly through complex-mediated mechanisms, while this ratio reaches 24 % at the highest energy due to the increasing dominance of the direct front-side mechanism and the smaller chance for hitting the deep Walden-inversion minimum. Our results suggest pronounced central-atom effects in SN 2 reactions, which can fundamentally change their (stereo)dynamics.

Theoretical Study of Structure and Photophysics of Homologous Series of Bis(arylydene)cycloalkanones

Photophysical properties of a series of bis(arylydene)cycloalkanone dyes with various donor substituents are studied using quantum chemistry. Their capacity for luminescence and nonradiative relaxation through trans-cis isomerization is related to their structure, in particular, to the donor capacity of the substituents and the degree of conjugation due to the central cycloalkanone moiety. It is shown that cyclohexanone central moiety introduces distortions and disrupts the conjugation, thus leading to a nonmonotonic change in their properties. The increasing donor capacity of the substituents causes increase in the HOMO energy (rise in the oxidation potential) and decrease in the HOMO-LUMO gap, which results in the red shift of the absorption spectra. The ability of the excited dye to relax through fluorescence or through trans-cis isomerization is governed by the height of the barrier between the Franck-Condon and S1-S0 conical intersection regions on the potential energy surface of the lowest π-π* excited state. This barrier also correlates with the donor capacity of the substituents and the degree of conjugation between the central and donor moieties. The calculated fluorescence and trans-cis isomerization rates are in good agreement with the observed fluorescence quantum yields.

Computational Insights into Excited State Intramolecular Double Proton Transfer Behavior Associated with Atomic Electronegativity for Bis(2'-benzothiazolyl)hydroquinone

Inspired by the distinguished regulated photochemical and photophysical properties of 2-(2'-hydroxyphenyl)benzazole derivatives, in this work, the novel bis(2'-benzothiazolyl)hydroquinone (BBTHQ) fluorophore is explored, looking at its photo-induced behaviors associated with different substituted atomic electronegativities, i.e., BBTHQ-SO, BBTHQ-SS and BBTHQ-Se compounds. From the structural changes, infrared (IR) vibrational variations and simulated core-valence bifurcation (CVB) indexes for the dual hydrogen bonds for the three BBTHQ derivatives, we see that low atomic electronegativity could be conducive to enhancing hydrogen bonding effects in the S1 state. Particularly, the O4-H5⋯N6 of BBTHQ-SO and the O1-H2⋯N3 of BBTHQ-SSe could be strengthened to be more intensive in the S1 state, respectively. Looking into the charge recombination induced by photoexcitation, we confirm a favorable ESDPT trend deriving from the charge reorganization of the dual hydrogen bonding regions. By constructing the potential energy surfaces (PESs) along with the ESDPT paths for the BBTHQ-SO, BBTHQ-SS and BBTHQ-Se compounds, we not only unveil stepwise ESDPT behaviors, but also present an atomic electronegativity-regulated ESDPT mechanism.

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.

Theoretical investigation of steric effects on the S1 potential energy surface of o-carborane-anthracene derivatives

TDDFT scan calculations were performed for s-carborane-anthracene derivatives (o-CB-X-Ant where X=-H, -CH₃, -C₂H₅ and tert-butyl or -tBu) in order to understand the interplay between the steric effects, S₁ potential energy surface (PES) and photophysical properties. The results show that all systems exhibit three local minima on the S₁ PES, which correspond to the emissive LE and TICT state, along with the nonemissive CT state respectively. In the case of the unsubstituted system (o-CB-H-Ant), and -CH₃ and -C₂H₅ substituted cases, S₁ PES is predicted to be quite flat for certain conformations indicating that it is possible for these systems to reach the nonemissive CT state without a large energy penalty. In comparison, conformational pathways for the nonemissive CT state are predicted to be energetically unfavorable for o-CB-tBu-Ant as a result of both steric and electronic effects. These results provide a mechanism for the enhanced emission of σ-CB-fluorophore molecules with bulky ligands.

Sequential Multistep Excited-State Structural Transformations in N,N'-Diphenyl-dihydrodibenzo[a,c]phenazine Fluorophores

We demonstrate that a single polycyclic π-scaffold can undergo sequential multistep excited-state structural evolution along the bent, planar, and twisted conformers, which coexist to produce intrinsic multiple fluorescence emissions in room-temperature solution. By installing a methyl or trifluoromethyl group on the ortho-site of N,N'-diphenyl-dihydrodibenzo[a,c]phenazine (DPAC), the enhanced steric effects change the fluorescence emission of DPAC from a dominant red band to well-resolved triple bands. The ultra-broadband triple emissions of ortho-substituted DPACs range from ~350 to ~850 nm, which is unprecedented for small fluorophores with molecular weight of < 500. Ultrafast spectroscopy and theoretical calculations clearly reveal that the above dramatic changes originate from the influence of steric hindrance on the shape of excited state potential energy surface (S1 PES). Compared to the steep S1 PES of parental DPAC, introduction of ortho-substituent is shown to make the path of structural evolution in S1 wider and flatter, so the ortho-substituted derivatives exhibit slower structural transformations from bent to planar and then to twisted forms, yielding intrinsic triple emission. The results provide the proof of concept that the bent, planar, and twisted emissive states can coexist in the same S1 PES, which greatly expand the fundamental understanding of the excited-state structural relaxation.

Planar Elongated B12 Structure in M3B12 Clusters (M = Cu-Au)

Here, it is shown that the M₃B₁₂ (M = Cu-Au) clusters' global minima consist of an elongated planar B₁₂ fragment connected by an in-plane linear M₃ fragment. This result is striking since this B₁₂ planar structure is not favored in the bare cluster, nor when one or two metals are added. The minimum energy structures were revealed by screening the potential energy surface using genetic algorithms and density functional theory calculations. Chemical bonding analysis shows that the strong electrostatic interactions with the metal compensate for the high energy spent in the M₃ and B₁₂ fragment distortion. Furthermore, metals participate in the delocalized π-bonds, which infers an aromatic character to these species.

On the Dynamics of the Carbon-Bromine Bond Dissociation in the 1-Bromo-2-Methylnaphthalene Radical Anion

This paper studies the mechanism of electrochemically induced carbon-bromine dissociation in 1-Br-2-methylnaphalene in the reduction regime. In particular, the bond dissociation of the relevant radical anion is disassembled at a molecular level, exploiting quantum mechanical calculations including steady-state, equilibrium and dissociation dynamics via dynamic reaction coordinate (DRC) calculations. DRC is a molecular-dynamic-based calculation relying on an ab initio potential surface. This is to achieve a detailed picture of the dissociation process in an elementary molecular detail. From a thermodynamic point of view, all the reaction paths examined are energetically feasible. The obtained results suggest that the carbon halogen bond dissociates following the first electron uptake follow a stepwise mechanism. Indeed, the formation of the bromide anion and an organic radical occurs. The latter reacts to form a binaphthalene intrinsically chiral dimer. This paper is respectfully dedicated to Professors Anny Jutand and Christian Amatore for their outstanding contribution in the field of electrochemical catalysis and electrosynthesis.

Cuts through the manifold of molecular H2O potential energy surfaces in liquid water at ambient conditions

Liquid water at ambient conditions is ubiquitous in chemistry and biology as well as in technology, energy, and atmospheric processes. Since parts of the phase diagram of water are unsettled—most notably the supercooled liquid homogeneous nucleation region—repercussions thereof on our molecular-level understanding for even the common ambient conditions remain. Breathtaking advances in X-ray–based approaches over the last decade give us now the tools to derive molecular potential energy surfaces as a quantitative view on the molecular manifold within the fluctuating hydrogen bonding network. With selective cuts along the local asymmetric O–H bond coordinate and the symmetric normal mode excitations an experimental foundation to benchmark competing molecular-level models of water has been achieved.

The fluctuating hydrogen bridge bonded network of liquid water at ambient conditions entails a varied ensemble of the underlying constituting H₂O molecular moieties. This is mirrored in a manifold of the H₂O molecular potentials. Subnatural line width resonant inelastic X-ray scattering allowed us to quantify the manifold of molecular potential energy surfaces along the H₂O symmetric normal mode and the local asymmetric O–H bond coordinate up to 1 and 1.5 Å, respectively. The comparison of the single H₂O molecular potentials and spectroscopic signatures with the ambient conditions liquid phase H₂O molecular potentials is done on various levels. In the gas phase, first principles, Morse potentials, and stepwise harmonic potential reconstruction have been employed and benchmarked. In the liquid phase the determination of the potential energy manifold along the local asymmetric O–H bond coordinate from resonant inelastic X-ray scattering via the bound state oxygen 1s to 4a1 resonance is treated within these frameworks. The potential energy surface manifold along the symmetric stretch from resonant inelastic X-ray scattering via the oxygen 1s to 2b2 resonance is based on stepwise harmonic reconstruction. We find in liquid water at ambient conditions H₂O molecular potentials ranging from the weak interaction limit to strongly distorted potentials which are put into perspective to established parameters, i.e., intermolecular O–H, H–H, and O–O correlation lengths from neutron scattering.

Current Status of the X + C2H6 [X ≡ H, F(2P), Cl(2P), O(3P), OH] Hydrogen Abstraction Reactions: A Theoretical Review

This paper is a detailed review of the chemistry of medium-size reactive systems using the following hydrogen abstraction reactions with ethane, X + C₂H₆ → HX + C₂H₅; X ≡ H, F(²P), Cl(²P), O(³P) and OH, and focusing attention mainly on the theoretical developments. These bimolecular reactions range from exothermic to endothermic systems and from barrierless to high classical barriers of activation. Thus, the topography of the reactive systems changes from reaction to reaction with the presence or not of stabilized intermediate complexes in the entrance and exit channels. The review begins with some reflections on the inherent problems in the theory/experiment comparison. When one compares kinetics or dynamics theoretical results with experimental measures, one is testing both the potential energy surface describing the nuclei motion and the kinetics or dynamics method used. Discrepancies in the comparison may be due to inaccuracies of the surface, limitations of the kinetics or dynamics methods, and experimental uncertainties that also cannot be ruled out. The paper continues with a detailed review of some bimolecular reactions with ethane, beginning with the reactions with hydrogen atoms. The reactions with halogens present a challenge owing to the presence of stabilized intermediate complexes in the entrance and exit channels and the influence of the spin-orbit states on reactivity. Reactions with O(³P) atoms lead to three surfaces, which is an additional difficulty in the theoretical study. Finally, the reactions with the hydroxyl radical correspond to a reactive system with ten atoms and twenty-four degrees of freedom. Throughout this review, different strategies in the development of analytical potential energy surfaces describing these bimolecular reactions have been critically analyzed, showing their advantages and limitations. These surfaces are fitted to a large number of ab initio calculations, and we found that a huge number of calculations leads to accurate surfaces, but this information does not guarantee that the kinetics and dynamics results match the experimental measurements.

Quantum Chemical Calculations to Trace Back Reaction Paths for the Prediction of Reactants

The long-due development of a computational method for the ab initio prediction of chemical reactants that provide a target compound has been hampered by the combinatorial explosion that occurs when reactions consist of multiple elementary reaction processes. To address this challenge, we have developed a quantum chemical calculation method that can enumerate the reactant candidates from a given target compound by combining an exhaustive automated reaction path search method with a kinetics method for narrowing down the possibilities. Two conventional name reactions were then assessed by tracing back the reaction paths using this new method to determine whether the known reactants could be identified. Our method is expected to be a powerful tool for the prediction of reactants and the discovery of new reactions.

Design and Nonadiabatic Photoisomerization Dynamics Study of a Three-Stroke Light-Driven Molecular Rotary Motor

Working cycle of conventional light-driven molecular rotary motors (LDMRMs), especially Feringa-type motors, usually have four steps, two photoisomerization steps, and two thermal helix inversion (THI) steps. THI steps hinder the ability of the motor to operate at lower temperatures and limit the rotation speed of LDMRMs. A three-stroke LDMRM, 2-(2,7-dimethyl-2,3-dihydro-1H-inden-1-ylidene)-1,2-dihydro-3H-pyrrol-3-one (DDIY), is proposed, which is capable of completing an unidirectional rotation by two photoisomerization steps and one thermal helix inversion step at room temperature. On the basis of trajectory surface-hopping simulation at the semi-empirical OM2/MRCI level, the EP→ZP and ZP→EM nonadiabatic photoisomerization dynamics of DDIY were systematically analyzed. Quantum yields of EP→ZP and ZP→EM photoisomerization of DDIY are ca. 34% and 18%, respectively. Both EP→ZP and ZP→EM photoisomerization processes occur on an ultrafast time scale (ca. 100-300 fs). This three-stroke LDMRM may stimulate further research for the development of new families of more efficient LDMRMs.

Search for transition states with external forces

It is much more difficult to find on the potential energy surface (PES) a transition state (TS) than a local minimum (LM). We propose a new methodology which makes this task much easier. Applying external forces to nuclei in a molecule we can locate on PES not just one particular TS but a number of consecutive transition states. With our approach it is possible to move over PES from one transition state to another without involving any local minima. The latter can be located in a separate step through the reaction path calculations performed for every transition state found before. Preliminary results for the 2-fluorofuran molecule illustrate the usefulness of the proposed method.

In-Depth Molecular Dynamics Study of All Possible Chondroitin Sulfate Disaccharides Reveals Key Insight into Structural Heterogeneity and Dynamism

GAGs exhibit a high level of conformational and configurational diversity, which remains untapped in terms of the recognition and modulation of proteins. Although GAGs are suggested to bind to more than 800 biologically important proteins, very few therapeutics have been designed or discovered so far. A key challenge is the inability to identify, understand and predict distinct topologies accessed by GAGs, which may help design novel protein-binding GAG sequences. Recent studies on chondroitin sulfate (CS), a key member of the GAG family, pinpointing its role in multiple biological functions led us to study the conformational dynamism of CS building blocks using molecular dynamics (MD). In the present study, we used the all-atom GLYCAM06 force field for the first time to explore the conformational space of all possible CS building blocks. Each of the 16 disaccharides was solvated in a TIP3P water box with an appropriate number of counter ions followed by equilibration and a production run. We analyzed the MD trajectories for torsional space, inter- and intra-molecular H-bonding, bridging water, conformational spread and energy landscapes. An in-house phi and psi probability density analysis showed that 1→3-linked sequences were more flexible than 1→4-linked sequences. More specifically, phi and psi regions for 1→4-linked sequences were held within a narrower range because of intra-molecular H-bonding between the GalNAc O5 atom and GlcA O3 atom, irrespective of sulfation pattern. In contrast, no such intra-molecular interaction arose for 1→3-linked sequences. Further, the stability of 1→4-linked sequences also arose from inter-molecular interactions involving bridged water molecules. The energy landscape for both classes of CS disaccharides demonstrated increased ruggedness as the level of sulfation increased. The results show that CS building blocks present distinct conformational dynamism that offers the high possibility of unique electrostatic surfaces for protein recognition. The fundamental results presented here will support the development of algorithms that help to design longer CS chains for protein recognition.

Benchmark study on deep neural network potentials for small organic molecules

There has been tremendous advancement in machine learning (ML) applications in computational chemistry, particularly in neural network potentials (NNP). NNPs can approximate potential energy surface (PES) as a high dimensional function by learning from existing reference data, thereby circumventing the need to solve the electronic Schrödinger equation explicitly. As a result, ML accelerates chemical space exploration and property prediction compared to quantum mechanical methods. Novel ML methods have the potential to provide efficient means for predicting the properties of molecules. However, this potential has been limited by the lack of standard comparative evaluations. In this work, we compare four selected models, that is, ANI, PhysNet, SchNet, and BAND-NN, developed to represent the PES of small organic molecules. We evaluate these models for their accuracy and transferability on two different test sets (i) Small organic molecules of up to eight-heavy atoms on which ANI and SchNet achieve root mean square error (RMSE) of 0.55 and 0.60 kcal/mol, respectively. (ii) On random selection of molecules from the GDB-11 database with 10-heavy atoms, ANI achieves RMSE of 1.17 kcal/mol and SchNet achieves RMSE of 1.89 kcal/mol. We examine their ability to produce smooth meaningful surface by performing PES scans for bond stretch, angle bend, and dihedral rotations on relatively large molecules to assess their possible application in molecular dynamics simulations. We also evaluate their performance for yielding minimum energy structures via geometry optimization using various minimization algorithms. All these models were also able to accurately differentiate different isomers of the same empirical formula C 10 H 20 . ANI and PhysNet achieve an RMSE of 0.29 and 0.52 kcal/mol, respectively, on C 10 H 20 isomers.

Vibrational Energy Transfer in CO+N2 Collisions: A Database for V-V and V-T/R Quantum-Classical Rate Coefficients

Knowledge of energy exchange rate constants in inelastic collisions is critically required for accurate characterization and simulation of several processes in gaseous environments, including planetary atmospheres, plasma, combustion, etc. Determination of these rate constants requires accurate potential energy surfaces (PESs) that describe in detail the full interaction region space and the use of collision dynamics methods capable of including the most relevant quantum effects. In this work, we produce an extensive collection of vibration-to-vibration (V-V) and vibration-to-translation/rotation (V-T/R) energy transfer rate coefficients for collisions between CO and N2 molecules using a mixed quantum-classical method and a recently introduced (A. Lombardi, F. Pirani, M. Bartolomei, C. Coletti, and A. Laganà, Frontiers in chemistry, 7, 309 (2019)) analytical PES, critically revised to improve its performance against ab initio and experimental data of different sources. The present database gives a good agreement with available experimental values of V-V rate coefficients and covers an unprecedented number of transitions and a wide range of temperatures. Furthermore, this is the first database of V-T/R rate coefficients for the title collisions. These processes are shown to often be the most probable ones at high temperatures and/or for highly excited molecules, such conditions being relevant in the modeling of hypersonic flows, plasma, and aerospace applications.

The Reaction of o-Benzyne with Vinylacetylene: An Unexplored Way to Produce Naphthalene

The mechanism and kinetics of the reaction of ortho-benzyne with vinylacetylene have been studied by ab initio and density functional CCSD(T)-F12/cc-pVTZ-f12//B3LYP/6-311G(d,p) calculations of the pertinent potential energy surface combined with Rice-Ramsperger-Kassel-Marcus - Master Equation calculations of reaction rate constants at various temperatures and pressures. Under prevailing combustion conditions, the reaction has been shown to predominantly proceed by the biradical acetylenic mechanism initiated by the addition of C₄ H₄ to one of the C atoms of the triple bond in ortho-benzyne by the acetylenic end, with a significant contribution of the concerted addition mechanism. Following the initial reaction steps, an extra six-membered ring is produced and the rearrangement of H atoms in this new ring leads to the formation of naphthalene, which can further dissociate to 1- or 2-naphthyl radicals. The o-C₆ H₄ +C₄ H₄ reaction is highly exothermic, by ∼143 kcal/mol to form naphthalene and by 31-32 kcal mol^-1 to produce naphthyl radicals plus H, but features relatively high entrance barriers of 9-11 kcal mol^-1 . Although the reaction is rather slow, much slower than the reaction of phenyl radical with vinylacetylene, it forms naphthalene and 1- and 2-naphthyl radicals directly, with their relative yields controlled by the temperature and pressure, and thus represents a viable source of the naphthalene core under conditions where ortho-benzyne and vinylacetylene are available.

Potential Energy Surfaces of Vinylogous Wolff Rearrangement: NBO Analysis and Molecular Dynamic Simulation

The potential energy surfaces of the Vinylogous Wolff Rearrangement, an alternative process for the Wolff Rearrangement that takes place β,γ-unsaturated diazoketones have been fully explored employing M062X model chemistry and in a complementary task by means of the CASSCF method. The NBO analysis has been invoked to reveal the alternations of orbital occupancies and their stabilization energies through some of the critical structures located on the pathways. The calculations establish a two steps process for each pathway involving a second higher energy transition state that seems likely to be the rate determining step in agreement with experimental results. Energy analysis indicates more feasibility of the pathway including diradical intermediate. The most striking feature of the study concerns the finding out a transition state providing a channel in between two pathways through interconversion of the two keys bicyclo[2.1.0]pentanone and diradical intermediates.

The Se-S Bond Formation in the Covalent Inhibition Mechanism of SARS-CoV-2 Main Protease by Ebselen-like Inhibitors: A Computational Study

The inhibition mechanism of the main protease (Mpro) of SARS-CoV-2 by ebselen (EBS) and its analog with a hydroxyl group at position 2 of the benzisoselenazol-3(2H)-one ring (EBS-OH) was studied by using a density functional level of theory. Preliminary molecular dynamics simulations on the apo form of Mpro were performed taking into account both the hydrogen donor and acceptor natures of the Nδ and Nε of His41, a member of the catalytic dyad. The potential energy surfaces for the formation of the Se-S covalent bond mediated by EBS and EBS-OH on Mpro are discussed in detail. The EBS-OH shows a distinctive behavior with respect to EBS in the formation of the noncovalent complex. Due to the presence of canonical H-bonds and noncanonical ones involving less electronegative atoms, such as sulfur and selenium, the influence on the energy barriers and reaction energy of the Minnesota hybrid meta-GGA functionals M06, M06-2X and M08HX, and the more recent range-separated hybrid functional wB97X were also considered. The knowledge of the inhibition mechanism of Mpro by the small protease inhibitors EBS or EBS-OH can enlarge the possibilities for designing more potent and selective inhibitor-based drugs to be used in combination with other antiviral therapies.

Energy contour exploration with potentiostatic kinematics

We introduce a method of exploring potential energy contours (PECs) in complex dynamical systems based on potentiostatic kinematics wherein the systems are evolved with minimal changes to their potential energy. We construct a simple iterative algorithm for performing potentiostatic kinematics, which uses an estimate curvature to predict new configuration-space coordinates on the PEC and a potentiostat term component to correct for errors in prediction. Our methods are then applied to atomic structure models using an interatomic potential for energy and force evaluations as would commonly be invoked in a molecular dynamics simulation. Using several model systems, we assess the stability and accuracy of the method on different hyperparameters in the implementation of the potentiostatic kinematics. Our implementation is open source and available within the atomic simulation environment package.

Spectroscopic and Theoretical Study of the Intramolecular π-Type Hydrogen Bonding and Conformations of 2-Cyclopenten-1-ol

The conformations of 2-cyclopenten-1-ol (2CPOL) have been investigated by high-level theoretical computations and infrared spectroscopy. The six conformational minima correspond to specific values of the ring-puckering and OH internal rotation coordinates. The conformation with the lowest energy possesses intramolecular π-type hydrogen bonding. A second conformer with weaker hydrogen bonding has somewhat higher energy. Ab initio coupled-cluster theory with single and double excitations (CCSD) was used with the cc-pVTZ (triple-ζ) basis set to calculate the two-dimensional potential energy surface (PES) governing the conformational dynamics along the ring-puckering and internal rotation coordinates. The two conformers with the hydrogen bonding lie about 300 cm⁻¹ (0.8 kcal/mole) lower in energy than the other four conformers. The lowest energy conformation has a calculated distance of 2.68 Å from the hydrogen atom on the OH group to the middle of the C=C double bond. For the other conformers, this distance is at least 0.3 Å longer. The infrared spectrum in the O-H stretching region agrees well with the predicted frequency differences between the conformers and shows the conformers with the hydrogen bonding to have the lowest values. The infrared spectra in other regions arise mostly from the two hydrogen-bonded species.

Quasi-Classical Trajectory Study of the CN + NH3 Reaction Based on a Global Potential Energy Surface

Based on a combination of valence-bond and molecular mechanics functions which were fitted to high-level ab initio calculations, we constructed an analytical full-dimensional potential energy surface, named PES-2020, for the hydrogen abstraction title reaction for the first time. This surface is symmetrical with respect to the permutation of the three hydrogens in ammonia, it presents numerical gradients and it improves the description presented by previous theoretical studies. In order to analyze its quality and accuracy, stringent tests were performed, exhaustive kinetics and dynamics studies were carried out using quasi-classical trajectory calculations, and the results were compared with the available experimental evidence. Firstly, the properties (geometry, vibrational frequency and energy) of all stationary points were found to reasonably reproduce the ab initio information used as input; due to the complicated topology with deep wells in the entrance and exit channels and a “submerged” transition state, the description of the intermediate complexes was poorer, although it was adequate to reasonably simulate the kinetics and dynamics of the title reaction. Secondly, in the kinetics study, the rate constants simulated the experimental data in the wide temperature range of 25–700 K, improving the description presented by previous theoretical studies. In addition, while previous studies failed in the description of the kinetic isotope effects, our results reproduced the experimental information. Finally, in the dynamics study, we analyzed the role of the vibrational and rotational excitation of the CN(v,j) reactant and product angular scattering distribution. We found that vibrational excitation by one quantum slightly increased reactivity, thus reproducing the only experimental measurement, while rotational excitation strongly decreased reactivity. The scattering distribution presented a forward-backward shape, associated with the presence of deep wells along the reaction path. These last two findings await experimental confirmation.

Instanton calculations of tunneling splittings in chiral molecules

We report the ground state tunneling splittings (ΔE_± ) of a number of axially chiral molecules using the ring-polymer instanton (RPI) method (J. Chem. Phys., 2011, 134, 054109). The list includes isotopomers of hydrogen dichalcogenides H₂ X₂ (X = O, S, Se, Te, and Po), hydrogen thioperoxide HSOH and dichlorodisulfane S₂ Cl₂ . Ab initio electronic-structure calculations have been performed on the level of second-order Møller-Plesset perturbation (MP2) theory either with split-valance basis sets or augmented correlation-consistent basis sets on H, O, S, and Cl atoms. Energy-consistent pseudopotential and corresponding triple zeta basis sets of the Stuttgart group are used on Se, Te, and Po atoms. The results are further improved using single point calculations performed at the coupled cluster level with iterative singles and doubles and perturbative triples amplitudes. When available for comparison, our computed values of ΔE_± are found to lie within the same order of magnitude as values reported in the literature, although RPI also provides predictions for H₂ Po₂ and S₂ Cl₂ , which have not previously been directly calculated. Since RPI is a single-shot method which does not require detailed prior knowledge of the optimal tunneling path, it offers an effective way for estimating the tunneling dynamics of more complex chiral molecules, and especially those with small tunneling splittings.

Spectroscopy and Scattering Studies Using Interpolated Ab Initio Potentials

The Born-Oppenheimer potential energy surface (PES) has come a long way since its introduction in the 1920s, both conceptually and in predictive power for practical applications. Nevertheless, nearly 100 years later-despite astonishing advances in computational power-the state-of-the-art first-principles prediction of observables related to spectroscopy and scattering dynamics is surprisingly limited. For example, the water dimer, (H₂O)₂, with only six nuclei and 20 electrons, still presents a formidable challenge for full-dimensional variational calculations of bound states and is considered out of reach for rigorous scattering calculations. The extremely poor scaling of the most rigorous quantum methods is fundamental; however, recent progress in development of approximate methodologies has opened the door to fairly routine high-quality predictions, unthinkable 20 years ago. In this review, in relation to the workflow of spectroscopy and/or scattering studies, we summarize progress and challenges in the component areas of electronic structure calculations, PES fitting, and quantum dynamical calculations.

Geometric analysis of anharmonic downward distortion following paths

A mathematical aspect of the anharmonic downward distortion following (ADDF) path is discussed. The ADDF method is utilized as an automated reaction path search method, which can explore transition state geometries on a potential energy surface from a potential minimum. We show that the maximum number of the ADD stationary paths intersecting the potential minimum is 2^f + 1  - 2, where f denotes the degree of freedom of the system. We also show that the bifurcation of the ADD stationary path is essential to detect all the transition states connected from a given minimum. The ADDF computation is demonstrated for a H₂ O molecule in which pitchfork bifurcation is observed.

Proton Transfer and Nitro Rotation Tuned Photoisomerization of Artificial Base Pair-ZP

Recently, the successful incorporation of artificial base pairs in genetics has made a significant progress in synthetic biology. The present work reports the proton transfer and photoisomerization of unnatural base pair ZP, which is synthesized from the pyrimidine analog 6-amino-5-nitro-3-(1-β-D-2'-deoxyribo-furanosyl)-2 (1H)-pyridone (Z) and paired with its Watson-Crick complement, the purine analog 2-amino-8-(1'-β-D-2'- deoxyribofuranosyl)-imidazo[1,2-a]-1,3,5-triazin-4(8H)-one (P). To explain the mechanism of proton transfer process, we constructed the relaxed potential energy surfaces (PESs) linking the different tautomers in both gas phase and solution. Our results show that the double proton transfer in the gas phase occurs in a concerted way both in S₀ and S₁ states, while the stepwise mechanism becomes more favorable in solution. The solvent effect can promote the single proton transfer, which undergoes a lower energy barrier in S₁ state due to the strengthened hydrogen bond. In contrast to the excited state ultrafast deactivation process of the natural bases, there is no conical intersection between S₀ and S₁ states along the proton transfer coordinate to activate the decay mechanism in ZP. Of particular relevance to the photophysical properties, charge-transfer character is obviously related to the nitro rotation in S₁ state. We characterized the molecular vibration effect on the electronic properties, which reveals the electronic excitation can be tuned by the rotation-induced structural distortion accompanied with the electron localization on nitro group.

Quantum Chemical Exploration of Intermolecular Reactions of Acetylene

Quantum chemical explorations of potential energy surfaces showed that acetylene produces various products starting from molecular arrays in short distances of 1.3-2.5 Å. Arrays of (C₂ H₂ )₂ gave cyclobutadiene, tetrahedrane, and acetylene dimers. Arrays of (C₂ H₂ )₃ gave benzene, prismane, benzvalene, Dewar benzene, and acetylene trimers. Arrays of (C₂ H₂ )₄ gave cubane, cyclooctatetranene, and acetylene tetramers. Different forms of initial arrays yielded different sets of products; a parallel array of two monomers gave cyclobutadiene, whereas a cross array gave tetrahedrane. Initial molecular arrays with unusually close contacts were estimated to require local forces of 1-9 × 10^-8  N. © 2019 Wiley Periodicals, Inc.

Chemical promenades: Exploring potential-energy surfaces with immersive virtual reality

The virtual-reality framework AVATAR (Advanced Virtual Approach to Topological Analysis of Reactivity) for the immersive exploration of potential-energy landscapes is presented. AVATAR is based on modern consumer-grade virtual-reality technology and builds on two key concepts: (a) the reduction of the dimensionality of the potential-energy surface to two process-tailored, physically meaningful generalized coordinates, and (b) the analogy between the evolution of a chemical process and a pathway through valleys (potential wells) and mountain passes (saddle points) of the associated potential energy landscape. Examples including the discovery of competitive reaction paths in simple A + BC collisional systems and the interconversion between conformers in ring-puckering motions of flexible rings highlight the innovation potential that augmented and virtual reality convey for teaching, training, and supporting research in chemistry.

The Generation of the Oxidant Agent of a Mononuclear Nonheme Fe(II) Biomimetic Complex by Oxidative Decarboxylation. A DFT Investigation

The oxidative decarboxylation of the iron(II) α-hydroxy acid (mandelic acid) complex model, biomimetic of Rieske dioxygenase, has been investigated at the density functional level. The explored mechanism sheds light on the role of the α-hydroxyl group on the dioxygen activation. The potential energy surfaces have been explored in different electronic spin states. The rate-determining step of the process is the proton transfer. The oxidative decarboxylation preferentially takes place on the quintet state.

Determinant of ESIPT Mechanism by the Structure Designed for Symmetrical and Unsymmetrical Molecules

In this work, we present that different structures lead to different excited state properties based on the investigations of the systems with symmetrical and unsymmetrical properties. Our work shows that the symmetrical and unsymmetrical compounds with a modified structure play key roles in regulating the excited state intramolecular proton transfer (PT) (ESIPT) process. For N,N'-bis(salicylidene)-p-phenyle-nediamine (BSP) and N,N'-disalicylidene-1,6-pyrenediamine (BSD), when they are excited, the torsion angels result in the rematch of the spectra. By analyzing the potential energy surfaces (PESs) of the torsion angle and PT process, we conclude that the size of the π stack largely affects the molecular properties in the excited states. For 2'-hydroxychalcone derivatives, which have important applications in biotransformation reactions, investigating the molecules of M1 and M2 could promote innovation in bioengineering. The results of molecular electrostatic potential (MEP) and the real space intramolecular interactions show the state and position of the hydrogen bond (HB) in both S₀ and S₁ states. The corresponding PT PESs show that the ESIPT reaction is much easier for the M1 system due to the lack of the side chain hydroxyl compared with the M2 system. This work is not only consistent with experimental results and explains its mechanism but also presents that symmetric and nonsymmetric structures modify their own potential regulation and controlling effects for ESIPT behaviors.

Computational Mechanistic Insights on the NO Oxidation Reaction Catalyzed by Non-Heme Biomimetic Cr-N-Tetramethylated Cyclam Complexes

The conversion reaction of NO to NO₃⁻ ion catalyzed by the end-on [Cr(III)(n-TMC)(O₂)(Cl)]⁺ superoxo and side-on [Cr(IV)(n-TMC)(O₂)(Cl)]⁺ peroxo non-heme complexes (n = 12, 13, 14 and 15), which are biomimetic systems of nitric oxide dioxygenases (NODs), has been explored using a computational protocol in the framework of density functional theory. Results show that the potential energy profiles for the studied reactions lie above the reagent energies, regardless of the used catalyst. Both the O-O bond breaking in the biomimetics and the NO₃⁻ ion formation require low energy barriers suggesting an efficient catalytic power of the studied systems. The rate-determining step depends on ligand size.

Augmenting Basin-Hopping With Techniques From Unsupervised Machine Learning: Applications in Spectroscopy and Ion Mobility

Evolutionary algorithms such as the basin-hopping (BH) algorithm have proven to be useful for difficult non-linear optimization problems with multiple modalities and variables. Applications of these algorithms range from characterization of molecular states in statistical physics and molecular biology to geometric packing problems. A key feature of BH is the fact that one can generate a coarse-grained mapping of a potential energy surface (PES) in terms of local minima. These results can then be utilized to gain insights into molecular dynamics and thermodynamic properties. Here we describe how one can employ concepts from unsupervised machine learning to augment BH PES searches to more efficiently identify local minima and the transition states connecting them. Specifically, we introduce the concepts of similarity indices, hierarchical clustering, and multidimensional scaling to the BH methodology. These same machine learning techniques can be used as tools for interpreting and rationalizing experimental results from spectroscopic and ion mobility investigations (e.g., spectral assignment, dynamic collision cross sections). We exemplify this in two case studies: (1) assigning the infrared multiple photon dissociation spectrum of the protonated serine dimer and (2) determining the temperature-dependent collision cross-section of protonated alanine tripeptide.

Mechanism and rate constants of the CH2 + CH2 CO reactions in triplet and singlet states: A theoretical study

Ab initio and density functional CCSD(T)-F12/cc-pVQZ-f12//B2PLYPD3/6-311G** calculations have been performed to unravel the reaction mechanism of triplet and singlet methylene CH₂ with ketene CH₂ CO. The computed potential energy diagrams and molecular properties have been then utilized in Rice-Ramsperger-Kassel-Marcus-Master Equation (RRKM-ME) calculations of the reaction rate constants and product branching ratios combined with the use of nonadiabatic transition state theory for spin-forbidden triplet-singlet isomerization. The results indicate that the most important channels of the reaction of ketene with triplet methylene lead to the formation of the HCCO + CH₃ and C₂ H₄ + CO products, where the former channel is preferable at higher temperatures from 1000 K and above. In the C₂ H₄ + CO product pair, the ethylene molecule can be formed either adiabatically in the triplet electronic state or via triplet-singlet intersystem crossing in the singlet electronic state occurring in the vicinity of the CH₂ COCH₂ intermediate or along the pathway of CO elimination from the initial CH₂ CH₂ CO complex. The predominant products of the reaction of ketene with singlet methylene have been shown to be C₂ H₄ + CO. The formation of these products mostly proceeds via a well-skipping mechanism but at high pressures may to some extent involve collisional stabilization of the CH₃ CHCO and cyclic CH₂ COCH₂ intermediates followed by their thermal unimolecular decomposition. The calculated rate constants at different pressures from 0.01 to 100 atm have been fitted by the modified Arrhenius expressions in the temperature range of 300-3000 K, which are proposed for kinetic modeling of ketene reactions in combustion. © 2018 Wiley Periodicals, Inc.

Singlet O2 Oxidation of a Deprotonated Guanine-Cytosine Base Pair and Its Entangling with Intra-Base-Pair Proton Transfer

We report an experimental and computational study on the ¹ O₂ oxidation of gas-phase deprotonated guanine-cytosine base pair [G ⋅ C-H]^- that is composed of 9HG ⋅ [C-H]^- and 7HG ⋅ [C-H]^- (pairing 9H- or 7H-guanine with N1-deprotonated cytosine), and 9HG ⋅ [C-H]^- _PT and 7HG ⋅ [C-H]^- _PT (formed by intra-base-pair proton transfer from the N1 of guanine to the N3 of [C-H]^- ). The conformer-averaged reaction product ions and cross section were measured over a center-of-mass collision energy range from 0.1 to 0.5 eV using a guided-ion-beam tandem mass spectrometer. To explore conformation-specific reactivity, collision dynamics of ¹ O₂ with each of the four [G ⋅ C-H]^- conformers was simulated at B3LYP/6-31G(d). Trajectories showed that the ¹ O₂ oxidation of the base pair entangles with intra-base-pair proton transfer, and prefers to occur in a collision when the base pair adopts a proton-transferred structure; trajectories also indicate that the 9HG-containing base pair favors stepwise formation of 4,8-endoperoxide of guanine, whereas the 7HG-containing base pair prefers concerted formation of guanine 5,8-endoperoxide. Using trajectory results as a guide, potential energy surfaces (PESs) along all possible reaction pathways were established using the approximately spin-projected ωB97XD/6-311++G(d,p)//B3LYP/6-311++G(d,p) method. PESs have not only rationalized trajectory findings but provided more accurate energetics and indicated that the proton-transferred base-pair conformers have lower activation barriers for oxidation than their non-proton-transferred counterparts.

Nucleophilic Substitution (SN 2): Dependence on Nucleophile, Leaving Group, Central Atom, Substituents, and Solvent

The reaction potential energy surface (PES), and thus the mechanism of bimolecular nucleophilic substitution (SN 2), depends profoundly on the nature of the nucleophile and leaving group, but also on the central, electrophilic atom, its substituents, as well as on the medium in which the reaction takes place. Here, we provide an overview of recent studies and demonstrate how changes in any one of the aforementioned factors affect the SN 2 mechanism. One of the most striking effects is the transition from a double-well to a single-well PES when the central atom is changed from a second-period (e. g. carbon) to a higher-period element (e.g, silicon, germanium). Variations in nucleophilicity, leaving group ability, and bulky substituents around a second-row element central atom can then be exploited to change the single-well PES back into a double-well. Reversely, these variations can also be used to produce a single-well PES for second-period elements, for example, a stable pentavalent carbon species.

Roads to pentazolate anion: a theoretical insight

The formation mechanism of pentazolate anion (PZA) is not yet clear. In order to present the possible formation pathways of PZA, the potential energy surfaces of phenylpentazole (PPZ), phenylpentazole radical (PPZ-R), phenylpentazole radical anion (PPZ-RA), PPZ and m-chloroperbenzoic acid (m-CPBA), p-pentazolylphenolate anion (p-PZPolA) and m-CPBA, and p-pentazolylphenol (p-PZPol) and m-CPBA were calculated by the computational electronic structure methods including the hybrid density functional, the double hybrid density functional and the coupled-cluster theories. At the thermodynamic point of view, the cleavages of C-N bonds of PPZ and PPZ-R need to absorb large amounts of heat. Thus, they are not feasible entrance for PZA formation at ambient condition. But excitation of PPZ and deprotonation of PPZ-RA probably happen before cleavage of C-N bond of PPZ at high-energy condition. As to the radical anion mechanism, the high accuracy calculations surveyed that the barrier of PZA formation is probably lower than that of dinitrogen evolution, but the small ionization potential of PPZ-RA gives rise to the unstable ionic pair of sodium PPZ at high temperature. In respect of oxidation mechanism, except for PPZ, the reactions of p-PZPolA and p-PZPol with m-CPBA can form PZA and quinone. The PZA formations have the barriers of about 20 kcal mol-1 which compete with the dinitrogen evolutions. The stabilities of PZA in both solid and gas phases were also studied herein. The proton prefers to transfer to pentazolyl group in the (N5)6(H3O)3(NH4)4Cl system which leads to the dissociation of pentazole ring. The ground states of M(N5)2(H2O)4 (M = Co, Fe and Mn) are high-spin states. The pentazolyl groups confined by the crystal waters in the coordinate compounds can improve the kinetic stability. As to the reactivity of PZA, it can be persistently oxidized by m-CPBA to oxo-PZA and 1,3-oxo-PZA with the barriers of about 20 kcal mol-1.

Probing the rate-determining region of the potential energy surface for a prototypical ion-molecule reaction

We report a joint experimental-theoretical study of the F^- + HCl → HF + Cl^- reaction kinetics. The experimental measurement of the rate coefficient at several temperatures was made using the selected ion flow tube method. Theoretical rate coefficients are calculated using the quasi-classical trajectory method on a newly developed global potential energy surface, obtained by fitting a large number of high-level ab initio points with augmentation of long-range electrostatic terms. In addition to good agreement between experiment and theory, analyses suggest that the ion-molecule reaction rate is significantly affected by shorter-range interactions, in addition to the traditionally recognized ion-dipole and ion-induced dipole terms. Furthermore, the statistical nature of the reaction is assessed by comparing the measured and calculated HF product vibrational state distributions to that predicted by the phase space theory.This article is part of the theme issue 'Modern theoretical chemistry'.

Automated reaction path searches for spin-forbidden reactions

Many catalytic and biomolecular reactions containing transition metals involve changes in the electronic spin state. These processes are referred to as "spin-forbidden" reactions within nonrelativistic quantum mechanics framework. To understand detailed reaction mechanisms of spin-forbidden reactions, one must characterize reaction pathways on potential energy surfaces with different spin states and then identify crossing points. Here we propose a practical computational scheme, where only the lowest mixed-spin eigenstate obtained from the diagonalization of the spin-coupled Hamiltonian matrix is used in reaction path search calculations. We applied this method to the ^6,4 FeO⁺  + H₂ → ^6,4 Fe⁺  + H₂ O, ^6,4 FeO⁺  + CH₄ → ^6,4 Fe⁺  + CH₃ OH, and ⁷ Mn⁺  + OCS → ⁵ MnS⁺  + CO reactions, for which crossings between the different spin states are known to play essential roles in the overall reaction kinetics. © 2018 Wiley Periodicals, Inc.

Development of a Flexible-Monomer Two-Body Carbon Dioxide Potential and Its Application to Clusters up to (CO2 )13

A flexible-monomer two-body potential energy function was developed that approaches the high level CCSD(T)/CBS potential energy surface (PES) of carbon dioxide (CO₂ ) systems. This function was generated by fitting the electronic energies of unique CO₂ monomers and dimers to permutationally invariant polynomials. More than 200,000 CO₂ configurations were used to train the potential function. Comparisons of the PESs of six orientations of flexible CO₂ dimers were evaluated to demonstrate the accuracy of the potential. Furthermore, the potential function was used to determine the minimum energy structures of CO₂ clusters containing as many as 13 molecules. For isomers of (CO₂ )₃ , the potential demonstrated energetic agreement with the M06-2X functional and structural agreement of the B2PLYP-D functional at substantially reduced computational costs. A separate function, fit to MP2/aug-cc-pVDZ reference energies, was developed to directly compare the two-body potential to the ab initio MP2 level of theory. © 2017 Wiley Periodicals, Inc.