Global Reaction Route Mapping on Potential Energy Surfaces of Formaldehyde, Formic Acid, and Their Metal-Substituted Analogues

Observation of the Infrared-Induced Structural Change in the Microscopic Hydrogen Bond Network of Phenol-Methanol Cluster Cations in a Cold-Ion Trap

To gain insight into microscopic hydrogen bond networks, we measured ultraviolet photodissociation (UVPD) spectra of the phenol-methanol 1:3 cluster cation, [PhOH(MeOH)3]+ trapped in a variable temperature ion trap. At low temperatures, an isomer with a ring-type hydrogen bond structure dominates, whereas at higher temperatures the chain-type isomers become dominant due to the flexibility of their hydrogen bond structures. We also found a clear temperature dependence of the spectral features, such as band position and width. In addition to the above measurement, we observed the infrared (IR) induced isomerization of [PhOH(MeOH)3]+ to study the dynamical aspects of hydrogen bond networks. We succeeded in observing IR-induced isomerization from the ring to chain forms of [PhOH(MeOH)3]+ at low temperature. The isomerization was clearly identified as a change in the UVPD spectra. The time evolution of the UVPD spectra after IR excitation indicated that the IR-induced isomerization occurs within a nanosecond. The chain-type isomers produced by the IR-induced isomerization are then converted back to the ring-type form by collisions with cold He buffer gas in the trap. This backward isomerization proceeds with a time constant of 100 μs under our experimental conditions. In this study, we evaluated the temperatures of the chain isomers during the backward isomerization on the basis of the spectral features.

Theoretical Analysis of Coordination Geometries in Transition Metal-Histidine Complexes Using Quantum Chemical Calculations

A theoretical investigation utilizing density functional theory (DFT) calculations was conducted to explore the coordination complexes formed between histidine (His) ligands and various divalent transition metal ions (Mn2+, Fe2+, Co2+, Ni2+, Cu2+, and Zn2+). Conformational exploration of the His ligand was initially performed to assess its stability upon coordination. Both 1:1 and 1:2 of metal-to-ligand complexes were scrutinized to elucidate their structural features and the relative stability of the complexes. This study examined the ability of His to act as a bidentate or tridentate coordinating ligand, along with the differences in coordination geometry when solvent effects were incorporated. The reduced density gradient (RDG) analysis and local electron attachment energy (LEAE) analysis were employed to elucidate the interaction planes and the nucleophilic and electrophilic properties. The electronic properties were analyzed through electrostatic potential (ESP) maps and natural population analysis (NPA) of atomic charge distributions. This computational study provides valuable insights into the diverse coordination modes of His and its interactions with divalent transition metal ions, contributing to a better understanding of the role of this amino acid ligand in the formation of transition metal complexes. The findings can aid in the design and construction of self-assembled structures involving His-metal coordination.

Quantum Chemical Investigation into the Structural Analysis and Calculated Raman Spectra of Amylose Modeled with Linked Glucose Molecules

This study presents a quantum chemical investigation into the structural analysis and calculated Raman spectra of modeled amylose with varying units of linked glucose molecules. We systematically examined the rotation of hydroxymethyl groups and intramolecular hydrogen bonds within these amylose models. Our study found that as the number of linked glucose units increases, the linear structure becomes more complex, resulting in curled, cyclic, or helical structures facilitated by establishing various intramolecular interactions. The hydroxymethyl groups were confirmed to form interactions with oxygen atoms and with hydroxymethyl and hydroxyl groups from adjacent rings in the molecular structures. We identified distinct peaks and selected specific bands applicable in various analytical contexts by comparing their calculated Raman spectra. Representative vibrational modes within selected regions were identified across the different lengths of amylose models, serving as characteristic signatures for linear and more coiled structural conformations. Our findings contribute to a deeper understanding of amylose structures and spectroscopic signatures, with implications for theoretical studies and potential applications. This work provides valuable reference points for the detailed assignment of Raman peaks of amylose structure, facilitating their application in broader research on carbohydrate structures and their associated spectroscopic properties.

Proton Transfer Processes in 2-Butenenitrile Dimer Cation Studied by Mass-Selective Infrared Spectroscopy

2-Butenenitrile (2-Bu) is a recently discovered crucial interstellar molecule. Herein, an abnormal NH band was observed in the infrared spectrum of the 2-Bu dimer cation, suggestive of a proton transfer reaction within the cluster. Through a comprehensive theoretical analysis of the IR spectrum of (2-Bu)2+, we discovered not only the formation of a new C-N bond through the attachment of one 2-Bu to another but also the occurrence of a proton transfer reaction in the cluster. This proton was identified as originating from the methyl group of the attaching 2-Bu in the cluster based on the analysis of IR spectra of (2-Bu)+ and [2-Bu-acrylonitrile (AN)]+. Furthermore, the detailed reaction process of this ion-molecule reaction is examined with theoretical calculation. This finding contributes significantly to our deeper understanding of ion-molecule reactions in the gas phase and the formation of nitrogen-containing prebiotic molecules in the interstellar medium.

Infrared Spectroscopy of Neutral and Cationic Benzonitrile-Methanol Binary Clusters in Supersonic Jets

The formation of nitrogen-containing organic interstellar molecules is of great importance to reveal chemical processes and the origin of life on Earth. Benzonitrile (BN) is one of the simplest nitrogen-containing aromatic molecules in the interstellar medium (ISM) that has been detected in recent years. Methanol (CH3OH) exists widely in interstellar space with high reactivity. Herein, we measured the infrared (IR) spectra of neutral and cationic BN-CH3OH clusters by vacuum ultraviolet (VUV) photoionization combined with time-of-flight mass spectrometry. Combining IR spectra with the density functional theory calculations, we reveal that the BN-CH3OH intends to form a cyclic H-bonded structure in neutral clusters. However, after the ionization of BN-CH3OH clusters, proton-shared N···H···O and N···H···C structures are confirmed to form between BN and CH3OH, with the minor coexistence of H-bond and O-π structures. The formation of the proton-shared structure expands our knowledge of the evolution of the life-related nitrogen-containing molecules in the universe and provides a possible pathway to the further study of biorelevant aromatic organic macromolecules.

Quantum Chemical Calculations of Flexible Tripeptide-Ni(II) Ion-Mediated Supramolecular Fragments and Comparative Analysis of Tripeptide Complexes with Various Metal(II) Ions

An automated conformational search method was employed to efficiently determine the stable conformers and weak hydrogen bonds of a flexible tripeptide coordinated with a solitary metal(II) ion in an aqueous environment. Quantum chemical calculations were performed to investigate the tendency of octahedral coordination formation between different metal(II) ions and various coordination models (ammonia molecule, chelate molecule, and flexible tripeptide). The octahedral coordination was analyzed by decomposing it into tridentate, bidentate, and monodentate coordination model complexes to assess their formation propensities and conformational properties. Additionally, population analysis, including electrostatic potential mapping and natural population analysis, was performed to identify the unique properties of the Ni(II) ion in forming octahedral coordination in crystals and to explore the potential of other metal(II) ions for self-assembling novel coordination configurations in peptide-metal compounds. Two common hydrogen bonding interactions were examined by using artificial forces to facilitate dissociation or reinforcement.

Single-Photon Ionization Induced New Covalent Bond Formation in Acrylonitrile(AN)-Pyrrole(Py) Clusters

The formation of nitrogen-containing organic compounds is crucial for understanding chemical evolution and the origin of life in the interstellar medium (ISM). In this study, we explore whether acrylonitrile (AN) and pyrrole (Py) can form new nitrogen-containing compounds after single-photon ionization in their gaseous clusters by vacuum ultraviolet (VUV)-infrared (IR) spectroscopy and theoretical calculations. The results show that a strong linear H-bond is formed in neutral AN-Py, while cyclic or bicyclic H-bonded networks are formed in the neutral AN-Py2 cluster. It is found that the structure containing a new C-C covalent bond between two moieties in (AN-Py)+ is formed besides the formation of H-bonded structures after AN-Py is ionized by VUV light. In (AN-Py2)+ cluster cations, new C-C or C-N covalent bonds tend to be formed between two Py, with (Py)2+ as the core in the cluster. The results reveal that new covalent bonds are more likely to be formed between two Py species when AN and Py are present in the cationic clusters. These results provide spectroscopic evidence of the formation of new nitrogen-containing organic compounds from AN and Py induced by VUV, which are helpful for our understanding of the formation of diverse prebiotic molecules in interstellar space.

Combined molecular dynamics and coordinate driving method for automatically searching complicated reaction pathways

The combined molecular dynamics and coordinate driving (MD/CD) method is updated and generalized in this work to broaden its applications in automatically searching reaction pathways for complicated reactions. In this updated version, MD simulations are performed with the GFN's family of methods to systematically sample conformers for almost any systems with atomic numbers Z ≤ 86. The improved CD procedure is greatly accelerated by applying a pre-screening stage at the semiempirical GFN2-xTB level. An automatic module based on the Marcus theory and its improved version (the Wolynes theory) is designed to include single electron transfer (SET) processes into reaction pathways. The capabilities of this method are demonstrated by exploring the most possible reaction pathways of three experimentally reported reactions: the organophosphine-catalyzed trans phosphinoboration, the Fe(II) complex-mediated C(sp2)-H borylation reaction, and the SET-triggered deaminative radical cross-coupling reaction. Comprehensive reaction networks are obtained for all three reactions with reasonable computational costs. Detailed mechanisms for these reactions can account for the reported experimental facts.

Designing main-group catalysts for low-temperature methane combustion by ozone

The catalytic combustion of methane at a low temperature is becoming increasingly key to controlling unburned CH₄ emissions from natural gas vehicles and power plants, although the low activity of benchmark platinum-group-metal catalysts hinders its broad application. Based on automated reaction route mapping, we explore main-group elements catalysts containing Si and Al for low-temperature CH₄ combustion with ozone. Computational screening of the active site predicts that strong Brønsted acid sites are promising for methane combustion. We experimentally demonstrate that catalysts containing strong Bronsted acid sites exhibit improved CH₄ conversion at 250 °C, correlating with the theoretical predictions. The main-group catalyst (proton-type beta zeolite) delivered a reaction rate that is 442 times higher than that of a benchmark catalyst (5 wt% Pd-loaded Al₂O₃) at 190 °C and exhibits higher tolerance to steam and SO₂. Our strategy demonstrates the rational design of earth-abundant catalysts based on automated reaction route mapping.

An automated reaction route mapping for the reaction of NO and active species on Ag4 clusters in zeolites

A computational investigation of the catalytic reaction on multinuclear sites is very challenging. Here, using an automated reaction route mapping method, the single-component artificial force induced reaction (SC-AFIR) algorithm, the catalytic reaction of NO and OH/OOH species over the Ag₄²⁺ cluster in a zeolite is investigated. The results of the reaction route mapping for H₂ + O₂ reveal that OH and OOH species are formed over the Ag₄²⁺ cluster via an activation barrier lower than that of OH formation from H₂O dissociation. Then, reaction route mapping is performed to examine the reactivity of the OH and OOH species with NO molecules over the Ag₄²⁺ cluster, resulting in the facile reaction path of HONO formation. With the aid of the automated reaction route mapping, the promotion effect of H₂ addition on the SCR reaction was computationally proposed (boosting the formation of OH and OOH species). In addition, the present study emphasizes that automated reaction route mapping is a powerful tool to elucidate the complicated reaction pathway on multi-nuclear clusters.

Concertedness and Activation Energy Control by Distal Methyl Group during Ring Contraction/Expansion in Scalarane-Type Sesterterpenoid Biosynthesis

Salmahyritisol A, similan A, and hippospongide A, which are scalarane-type sesterterpenoids, feature 6/6/5/7/5 pentacyclic skeletons. Although their biosyntheses have been previously proposed to involve a unique skeletal rearrangement reaction, the detailed reaction mechanism remains unclear as none of the corresponding biosynthetic enzymes for this reaction have been reported. Herein, this skeletal rearrangement reaction was investigated using computational techniques, which revealed the following four key features: (i) the distal 24-Me substituent controls both the concertedness and activation energy of this transformation, (ii) enzymes are not responsible for the observed regioselectivity of C12-C20 bond formation, (iii) stereoselectivity is enzyme-regulated, and (iv) protonation is a key step in this skeletal rearrangement process. These new findings provide insight into the C-ring-contraction and D-ring-expansion mechanisms in scalarane-type sesterterpenoid biosyntheses.

Global Search for Stable C4H5NO Compounds-Guinness Molecules and Stability Islands

Global reaction route mapping (GRRM) analysis for compounds with the formula C₄H₅NO allowed for the detection of the corresponding "Guinness molecules" 000 and 001, as well as around 150 other stable minima of the same composition. The results suggest that compounds of similar functionality form a kind of "Stability Island" with their free energies of formation falling within s relatively limited range.

Calculation of the permeability coefficients of small molecules through lipid bilayers by free-energy reaction network analysis following the explicit treatment of the internal conformation of the solute

Biomembrane permeation represents a major barrier to pharmacokinetics. During preclinical drug discovery, the coefficients of the permeation of molecules through lipid bilayers account for a valuable property of such molecules. Therefore, the control of the permeation of molecules through lipid bilayers is an essential factor in drug design, and the estimation of the permeation phenomena is a crucial study in pharmacy. Thus, there are many published studies on the theoretical estimations of permeation coefficients. Here, we propose a molecular dynamics (MD) simulation method for estimating the permeation of small molecules through lipid bilayers based on the free-energy reaction network (FERN) analysis. In this method, the collective variables (CVs) of the free energy calculations explicitly include the conformational changes in the rotational bonds of the solute molecules. The advantages of the present method over the other method are that it is possible to estimate reaction pathways and their reaction rates, i.e., permeation coefficients or passage times, in multidimensional space spanned by CVs though conventional methods such as the umbrella sampling method and target MDs often dealt with a few degrees of freedom. To demonstrate the efficacy of our method, we calculate the coefficients of the permeation of three small aromatic peptides, namely N-acetylphenylalanineamide (Ac-Phe-NH2 or NAFA), N-acetyltyrosineamide (Ac-Tyr-NH2 or NAYA), and N-acetyltryptophanamide (Ac-Trp-NH2 or NATA), through a 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) lipid bilayer. In these cases we adopted one CV for the permeation direction and four CVs for the internal rotational coordinates. The results reveal that the permeation coefficients of NAFA, NAYA, and NATA are 1.7 × 10^-2, 0.51 × 10^-4, and 5.7 × 10^-4 cm s^-1, respectively. Compared with the experimental data, our simulation results followed the same trend, i.e., NAFA > NATA > NAYA. By analyzing the structures of metastable points of the solute molecules, our simulation result reveals that the aforementioned trend is caused by the differences in stability among their rotamers. Furthermore, we evaluate the statistical fluctuation of the rotamers, and the time scale of flipping the side chain reveals that the structures rigidify as the ligand moves deeper into the membrane.

Molecular Conformer Search with Low-Energy Latent Space

Identifying low-energy conformers with quantum mechanical accuracy for molecules with many degrees of freedom is challenging. In this work, we use the molecular dihedral angles as features and explore the possibility of performing molecular conformer search in a latent space with a generative model named variational auto-encoder (VAE). We bias the VAE towards low-energy molecular configurations to generate more informative data. In this way, we can effectively build a reliable energy model for the low-energy potential energy surface. After the energy model has been built, we extract local-minimum conformations and refine them with structure optimization. We have tested and benchmarked our low-energy latent-space (LOLS) structure search method on organic molecules with 5–9 searching dimensions. Our results agree with previous studies.

Spectroscopic evidence of S∴N and S∴O hemibonds in heterodimer cations

Computational and condensed phase experimental evidence for the existence of S∴N and S∴O hemibonded structures has been reported previously, but no gas phase experimental evidence has been reported. To experimentally explore the existence of the S∴N and S∴O hemibonds in the gas phase, we recorded the infrared photodissociation action spectra of four cationic clusters: [CH₃SH-NH₃]⁺, [CH₃SCH₃-NH₃]⁺, [CH₃SCH₃-H₂O]⁺, and [CH₃OCH₃-H₂O]⁺. Combined with the calculation results, it is found that the S∴N hemibonded structure is competitive with the S⋯HN H-bonded structure, though only the latter structure is actually observed in [CH₃SH-NH₃]⁺. The spectral and theoretical results show that hemibonds can form between the second- (oxygen or nitrogen) and the third-period elements (sulfur) in the heterodimer clusters of [CH₃SCH₃-NH₃]⁺ and [CH₃SCH₃-H₂O]⁺. However, the S∴N and S∴O hemibonded structures are found competitive with the C⋯HN and CH⋯O H-bonded structures, respectively, and both the structures coexist. On the other hand, the O∴O hemibonded structure is much less stable than other hydrogen bonded (H-bonded) structures in [CH₃OCH₃-H₂O]⁺, and it shows no clear contribution to the observed spectrum. This study provides direct spectroscopic evidence for the existence of S∴N and S∴O hemibonds in the gas phase and their competition with the H-bonds, which may be also fundamentally important in biological processes.

DFT Study on the Biosynthesis of Asperterpenol and Preasperterpenoid Sesterterpenoids: Exclusion of Secondary Carbocation Intermediates and Origin of Structural Diversification

The biosynthetic pathway to asperterpenol, a sesterterpenoid featuring a 6/6/8/5 tetracyclic ring system, was proposed to involve three secondary (2°) carbocation intermediates (B, D, and I), but it remains controversial whether or not these are viable. Further, the proposed 11/6/5 tricyclic intermediate C has the same "ChemDraw" structure as an intermediate in the biosynthesis of preasperterpenoid, which has a very different 5/7/(3)6/5 pentacyclic skeleton. Here, we present a detailed scrutiny of the asperterpenol/preasperterpenoid biosynthetic pathways based on comprehensive DFT calculations.

Relative Populations and IR Spectra of Cu38 Cluster at Finite Temperature Based on DFT and Statistical Thermodynamics Calculations

The relative populations of Cu₃₈ isomers depend to a great extent on the temperature. Density functional theory and nanothermodynamics can be combined to compute the geometrical optimization of isomers and their spectroscopic properties in an approximate manner. In this article, we investigate entropy-driven isomer distributions of Cu₃₈ clusters and the effect of temperature on their IR spectra. An extensive, systematic global search is performed on the potential and free energy surfaces of Cu₃₈ using a two-stage strategy to identify the lowest-energy structure and its low-energy neighbors. The effects of temperature on the populations and IR spectra are considered via Boltzmann factors. The computed IR spectrum of each isomer is multiplied by its corresponding Boltzmann weight at finite temperature. Then, they are summed together to produce a final temperature-dependent, Boltzmann-weighted spectrum. Our results show that the disordered structure dominates at high temperatures and the overall Boltzmann-weighted spectrum is composed of a mixture of spectra from several individual isomers.

New Tools for Taming Complex Reaction Networks: The Unimolecular Decomposition of Indole Revisited

The level of detail attained in the computational description of reaction mechanisms can be vastly improved through tools for automated chemical space exploration, particularly for systems of small to medium size. Under this approach, the unimolecular decomposition landscape for indole was explored through the automated reaction mechanism discovery program AutoMeKin. Nevertheless, the sheer complexity of the obtained mechanisms might be a hindrance regarding their chemical interpretation. In this spirit, the new Python library amk-tools has been designed to read and manipulate complex reaction networks, greatly simplifying their overall analysis. The package provides interactive dashboards featuring visualizations of the network, the three-dimensional (3D) molecular structures and vibrational normal modes of all chemical species, and the corresponding energy profiles for selected pathways. The combination of the joined mechanism generation and postprocessing workflow with the rich chemistry of indole decomposition enabled us to find new details of the reaction (obtained at the CCSD(T)/aug-cc-pVTZ//M06-2X/MG3S level of theory) that were not reported before: (i) 16 pathways leading to the formation of HCN and NH₃ (via amino radical); (ii) a barrierless reaction between methylene radical and phenyl isocyanide, which might be an operative mechanism under the conditions of the interstellar medium; and (iii) reaction channels leading to both hydrogen cyanide and hydrogen isocyanide, of potential astrochemical interest as the computed HNC/HCN ratios greatly exceed the calculated equilibrium value at very low temperatures. The reported reaction networks can be very valuable to supplement databases of kinetic data, which is of remarkable interest for pyrolysis and astrochemical studies.

In vitro and in silico studies of potential coronavirus-specific 3C-like protease inhibitors

We investigated similar compounds to ebselen and tideglusib, which exhibit strong activity against severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), using MACCS keys. Four candidate compounds were identified. One of them, phenyl-benzothiazol-3-one, showed coronavirus-specific 3C-like (3CL) protease inhibitory activity. The results indicated that a similarity score above 0.81 is a good indicator of activity for ebselen-and-tideglusib-like compounds. Subsequently, we simulated the ring-cleavage Michael reaction of ebselen at the Se center, which is responsible for its 3CL protease inhibitory activity, and determined the activation free energy of the reaction. The results showed that reaction simulation is a useful tool for estimating the activity of inhibitory compounds that undergo Michael addition reactions with the relevant cysteine S atom of 3CL proteases.

Theoretical Study on the Mechanism of Spirocyclization in Spiroviolene Biosynthesis

Spiroviolene is a spirocyclic triquinane diterpene produced by Streptomyces violens. Recently, a biosynthetic pathway that includes secondary carbocation intermediates and a complicated concerted skeletal rearrangement was proposed for spiroviolene, based upon careful labeling experiments. On the basis of density functional theory (DFT) calculations, we propose a revised pathway for spiroviolene biosynthesis, involving a multistep carbocation cascade that bypasses the formation of unstable secondary carbocations by breaking the adjacent C-C bond to form a more stable tertiary carbocation (IM3) and by Wagner-Meerwein 1,2-methyl rearrangement (IM7).

Predicting potentially hazardous chemical reactions using an explainable neural network

Predicting potentially dangerous chemical reactions is a critical task for laboratory safety. However, a traditional experimental investigation of reaction conditions for possible hazardous or explosive byproducts entails substantial time and cost, for which machine learning prediction could accelerate the process and help detailed experimental investigations. Several machine learning models have been developed which allow the prediction of major chemical reaction products with reasonable accuracy. However, these methods may not present sufficiently high accuracy for the prediction of hazardous products which particularly requires a low false negative result for laboratory safety in order not to miss any dangerous reactions. In this work, we propose an explainable artificial intelligence model that can predict the formation of hazardous reaction products in a binary classification fashion. The reactant molecules are transformed into substructure-encoded fingerprints and then fed into a convolutional neural network to make the binary decision of the chemical reaction. The proposed model shows a false negative rate of 0.09, which can be compared with 0.47-0.66 using the existing main product prediction models. To provide explanations for what substructures of the given reactant molecules are important to make a decision for target hazardous product formation, we apply an input attribution method, layer-wise relevance propagation, which computes the contributions of individual inputs per input data. The computed attributions indeed match some of the existing chemical intuitions and mechanisms, and also offer a way to analyze possible data-imbalance issues of the current predictions based on relatively small positive datasets. We expect that the proposed hazardous product prediction model will be complementary to existing main product prediction models and experimental investigations.

DFT Study on the Biosynthesis of Verrucosane Diterpenoids and Mangicol Sesterterpenoids: Involvement of Secondary-Carbocation-Free Reaction Cascades

Some experimental observations indicate that a sequential formation of secondary (2°) carbocations might be involved in some biosynthetic pathways, including those of verrucosane-type diterpenoids and mangicol-type sesterterpenoids, but it remains controversial whether or not such 2° cations are viable intermediates. Here, we performed comprehensive density functional theory calculations of these biosynthetic pathways. The results do not support previously proposed pathways/mechanisms: in particular, we find that none of the putative 2° carbocation intermediates is involved in either of the biosynthetic pathways. In verrucosane biosynthesis, the proposed 2° carbocations (II and IV) in the early stage are bypassed by the formation of the adjacent 3° carbocations and by unusual skeletal rearrangement reactions, and in the later stage, the putative 2° carbocation intermediates (VI, VII, and VIII) are not present as the proposed forms but as nonclassical structures between homoallyl and cyclopropylcarbinyl cations. In the mangicol biosynthesis, one of the two proposed 2° carbocations (X) is bypassed by a C-C bond-breaking reaction to generate a 3° carbocation with a C=C bond, while the other (XI) is bypassed by a strong hyperconjugative interaction leading to a nonclassical carbocation. We propose new biosynthetic pathways/mechanisms for the verrucosane-type diterpenoids and mangicol-type sesterterpenoids. These pathways are in good agreement with the findings of previous biosynthetic studies, including isotope-labeling experiments and byproducts analysis, and moreover can account for the biosynthesis of related terpenes.

AutoMeKin2021: An open-source program for automated reaction discovery

AutoMeKin2021 is an updated version of tsscds2018, a program for the automated discovery of reaction mechanisms (J. Comput. Chem. 2018, 39, 1922). This release features a number of new capabilities: rare-event molecular dynamics simulations to enhance reaction discovery, extension of the original search algorithm to study van der Waals complexes, use of chemical knowledge, a new search algorithm based on bond-order time series analysis, statistics of the chemical reaction networks, a web application to submit jobs, and other features. The source code, manual, installation instructions and the website link are available at: https://rxnkin.usc.es/index.php/AutoMeKin.

Application of Reaction Path Search Calculations to Potential Energy Surface Fits

There has been significant progress in recent years in the use of machine learning techniques to model high-dimensional reactive potential energy surfaces using large-scale data obtained from ab initio electronic structure calculations. In these methods, the strategy used to gather data becomes a key issue as the molecular size increases. In this work, we examine the applicability of the reaction path search algorithm implemented in the Global Reaction Route Mapping (GRRM) code as a data-gathering approach. The electronic energies and gradients sampled by using the GRRM calculation are directly used in potential energy surface fitting to a permutationally invariant polynomial function. This simple approach was applied to the HNS and HCNO reaction systems, and we found that the fitted potential energy surfaces reasonably reproduce the features of the electronic structure calculations used in the GRRM calculations. This suggests that the GRRM sampling scheme can be used to construct an initial potential energy surface.

Proton Transfer in Nitromethane-Ammonia Clusters under VUV Single-Photon Ionization Explored by Infrared Spectroscopy and Theoretical Calculations

It is known that the acidity and reactivity of the CH bond can be enhanced after ionization. Also, this property plays a pivotal role in proton transfer reaction and in the formation of new molecules. Herein, infrared spectroscopy and high-precision quantum chemical calculations are used to study the neutral and cationic clusters of nitromethane-ammonia (CH₃NO₂-NH₃). It is found that in the neutral cluster, CH₃NO₂ and NH₃ are mainly bonded by three intermolecular hydrogen bonds, in which electrostatic contribution plays a major role. After vacuum ultraviolet (VUV) single-photon ionization of CH₃NO₂-NH₃, the positive charge redistributes from the ionized nitrogen atom of NH₃ to the CH₃NO₂ molecule immediately. Then, the proton of CH₃NO₂ transfers to NH₃ to form a proton-transferred type structure CH₂NO₂-NH₄⁺, without any effective energy barrier, due to the positive hyperconjugation of cationic nitromethane. A closed loop of positive charge transfer takes place in the CH₃NO₂-NH₃ cluster after VUV ionization. The present work demonstrates that both the proton transfer reaction and charge transfer process have occurred in the ionized CH₃NO₂-NH₃ cluster. Moreover, it is found that the proton transfer reaction is a result of the highly acidic CH bond caused by hyperconjugation between the σ (CH) bond and π orbital.

Spectroscopic evidence of the C-N covalent bond formed between two interstellar molecules (ISM): acrylonitrile and ammonia

Acrylonitrile (AN) and ammonia (NH3) are two important nitrogen-containing interstellar molecules in outer space, especially on Titan. Herein, we measured infrared (IR) spectra of neutral and cationic AN-NH3 complexes by VUV single-photon ionization combined with time-of-flight mass spectrometry. On combining IR spectra with the theoretical calculations, we found that the molecules prefer to form a single-ring cyclic H-bonded structure in the neutral AN-NH3 and (AN)2-NH3 clusters. However, after ionization of AN-NH3 and (AN)2-NH3 clusters, a new C-N-covalent bond is confirmed to form directly between AN and NH3, without any energy barrier in the cationic complexes. Moreover, in the ionized (AN)2-NH3 cluster, the covalent C-N bond prefers to form between AN and NH3 rather than the two AN groups. These results provide spectroscopic evidence of AN forming a new molecule with NH3, induced by VUV radiation. The formation of the new C-N bond broadens our knowledge on the evolution of the prebiotic nitrogen-containing molecules in space.

Application of Optimization Algorithms in Clusters

The structural characterization of clusters or nanoparticles is essential to rationalize their size and composition-dependent properties. As experiments alone could not provide complete picture of cluster structures, so independent theoretical investigations are needed to find out a detail description of the geometric arrangement and corresponding properties of the clusters. The potential energy surfaces (PES) are explored to find several minima with an ultimate goal of locating the global minima (GM) for the clusters. Optimization algorithms, such as genetic algorithm (GA), basin hopping method and its variants, self-consistent basin-to-deformed-basin mapping, heuristic algorithm combined with the surface and interior operators (HA-SIO), fast annealing evolutionary algorithm (FAEA), random tunneling algorithm (RTA), and dynamic lattice searching (DLS) have been developed to solve the geometrical isomers in pure elemental clusters. Various model or empirical potentials (EPs) as Lennard-Jones (LJ), Born-Mayer, Gupta, Sutton-Chen, and Murrell-Mottram potentials are used to describe the bonding in different type of clusters. Due to existence of a large number of homotops in nanoalloys, genetic algorithm, basin-hopping algorithm, modified adaptive immune optimization algorithm (AIOA), evolutionary algorithm (EA), kick method and Knowledge Led Master Code (KLMC) are also used. In this review the optimization algorithms, computational techniques and accuracy of results obtained by using these mechanisms for different types of clusters will be discussed.

High performance global exploration of isomers and isomerization channels on quantum chemical potential energy surface of H5 C2 NO2

High performance global exploration of isomers and isomerization channels on the quantum chemical potential energy surface (PES) is performed for H₅ C₂ NO₂ by using the scaled hypersphere search-anharmonic downward distortion following (SHS-ADDF) method. A multi-node operation, NeoGRRM, has achieved high performance exploration calculations for the large system by submitting SHS-ADDF sub-jobs into many cores in parallel and unifying the results of sub-jobs into the total lists of the main-job. Global exploration of equilibrium (EQ) and transition-state structures at the level of B3LYP/6-31G(d) gave 3210 EQs and 23278 TSs. Nine compounds were found in the low energy regions of 0-100 kJ/mol; the lowest energy compound is N-methylcarbamic acid, the second is methyl carbamate, and the third is glycine (the most fundamental amino acid). Interconversion pathways between the conformers of each of the low energy compounds were surveyed. Isomerization channels around glycine were explored in detail. The lowest energy barriers around some of the EQs turned to be negative after zero-point energy corrections. This indicates that those structures cannot exist as independent structures because they spontaneously collapse into more stable structures. The global PES search showed various interesting dissociating channels which indicate synthon reaction pathways in the reverse directions.

Exploration of Free Energy Surface and Thermal Effects on Relative Population and Infrared Spectrum of the Be6B11- Flux-Ional Cluster

The starting point to understanding cluster properties is the putative global minimum and all the nearby local energy minima; however, locating them is computationally expensive and difficult. The relative populations and spectroscopic properties that are a function of temperature can be approximately computed by employing statistical thermodynamics. Here, we investigate entropy-driven isomers distribution on Be₆B₁₁⁻ clusters and the effect of temperature on their infrared spectroscopy and relative populations. We identify the vibration modes possessed by the cluster that significantly contribute to the zero-point energy. A couple of steps are considered for computing the temperature-dependent relative population: First, using a genetic algorithm coupled to density functional theory, we performed an extensive and systematic exploration of the potential/free energy surface of Be₆B₁₁⁻ clusters to locate the putative global minimum and elucidate the low-energy structures. Second, the relative populations’ temperature effects are determined by considering the thermodynamic properties and Boltzmann factors. The temperature-dependent relative populations show that the entropies and temperature are essential for determining the global minimum. We compute the temperature-dependent total infrared spectra employing the Boltzmann factor weighted sums of each isomer’s infrared spectrum and find that at finite temperature, the total infrared spectrum is composed of an admixture of infrared spectra that corresponds to the spectra of the lowest-energy structure and its isomers located at higher energies. The methodology and results describe the thermal effects in the relative population and the infrared spectra.

Cu(I)/sucrose-catalyzed hydroxylation of arenes in water: the dual role of sucrose

A protocol for the hydroxylation of aryl halides catalyzed by copper(i) and sucrose in neat water has been developed. The dual role of sucrose, the reaction pathway, and the high selectivity for hydroxylation were investigated using a combination of experimental and theoretical techniques.

A comprehensive infrared spectroscopic and theoretical study on phenol-ethyldimethylsilane dihydrogen-bonded clusters in the S0 and S1 states

To investigate microscopic characters of Si-H⋯H-O type dihydrogen bonds, we observed OH and SiH stretch bands in both the S₀ and S₁ states of phenol-ethyldimethylsilane (PhOH-EDMS) clusters by infrared (IR)-ultraviolet (UV) and UV-IR double resonance spectroscopies. Density functional theory (DFT) calculations and energy decomposition analysis were also performed. Structures of two isomers identified were unambiguously determined through the analysis of IR spectra and DFT calculations. To discuss the strength of dihydrogen bond in various systems, we performed theoretical calculations on clusters of EDMS with several acidic molecules in addition to PhOH. It was revealed that charge-transfer interaction energies from a bonding σ orbital of SiH bond to an anti-bonding σ^* orbital of OH bond well reflect strengths of dihydrogen bonds. Additionally, it was found that the red shift of SiH stretch frequencies can be used as a crude measure of the strength of dihydrogen bonds. Relationship between the red shifts of OH/SiH stretch frequencies and various electrostatic components of the interaction energy was examined. In the S₁ state, large increases in red shifts were observed for both the OH and SiH stretch frequencies. Since the EDMS moiety is not associated with the electronic excitation in a cluster, the strength of dihydrogen bonds in the S₁ and S₀ states was able to be directly compared based on the red shifts of the SiH stretch bands. A significant increase in the red shift of SiH stretch frequency indicates a strengthening of the dihydrogen bonds during the electronic excitation of the PhOH moiety.

Water Radical Cations in the Gas Phase: Methods and Mechanisms of Formation, Structure and Chemical Properties

Water radical cations, (H₂O)_n^+•, are of great research interest in both fundamental and applied sciences. Fundamental studies of water radical reactions are important to better understand the mechanisms of natural processes, such as proton transfer in aqueous solutions, the formation of hydrogen bonds and DNA damage, as well as for the discovery of new gas-phase reactions and products. In applied science, the interest in water radicals is prompted by their potential in radiobiology and as a source of primary ions for selective and sensitive chemical ionization. However, in contrast to protonated water clusters, (H₂O)_nH⁺, which are relatively easy to generate and isolate in experiments, the generation and isolation of radical water clusters, (H₂O)_n^+•, is tremendously difficult due to their ultra-high reactivity. This review focuses on the current knowledge and unknowns regarding (H₂O)_n^+• species, including the methods and mechanisms of their formation, structure and chemical properties.

Analytical Method Using a Scaled Hypersphere Search for High-Dimensional Metadynamics Simulations

Metadynamics (MTD) is one of the most effective methods for calculating the free energy surface and finding rare events. Nevertheless, numerous studies using MTD have been carried out using 3D or lower dimensional collective variables (CVs), as higher dimensional CVs require costly computational resources and the obtained results are too complex to understand the important events. The latter issue can be conveniently solved by utilizing the free energy reaction network (FERN), which is a graph structure consisting of edges of minimum free energy paths (MFEPs), nodes of equation (EQ) points, and transition state (TS) points. In the present article, a new method for exploring FERNs on high-dimensional CVs using MTD and the scaled hypersphere search (SHS) method is described. A test calculation based on the MTD-SHS simulation of met-enkephalin in explicit water with 7 CVs was conducted. As a result, 889 EQ points and 1805 TS points were found. The MTD-SHS approach can find MFEPs exhaustively; therefore, the FERNs can be estimated without any a priori knowledge of the EQ and TS points.

An effective Hamiltonian analysis of a Franck-Condon-like pattern in the IR spectra of phenol-alkylsilane dihydrogen-bonded clusters in the S1 state

Infrared (IR) spectra in a region of the OH stretch band of phenol (PhOH)-ethyldimethylsilane (EDMS), phenol (PhOH)-triethylsilane (TES), and phenol (PhOH)-t-butyldimethylsilane (BDMS) dihydrogen-bonded clusters in the S₁ state were observed. All of the species exhibited unconventional band patterns in which many combination bands appeared with comparable intensities to those of allowed bands. Such a behavior is sometimes called a Franck-Condon-like pattern. In the case of the PhOH-BDMS, one intermolecular vibrational mode is involved in this behavior. The observed IR spectra were well reproduced based on the concept of the Franck-Condon-like behavior. As an alternative treatment, we analyzed the band patterns on the concept of intensity borrowing due to the vibrational anharmonic interaction. The analysis was based on an effective Hamiltonian involving an anharmonic interaction between the OH stretch and intermolecular vibrational modes. Two treatments provided the same results. Thus, it was confirmed that the Franck-Condon-like behavior originates from vibrational anharmonic interactions. In the cases of the PhOH-EDMS and PhOH-TES, we carried out a two-dimensional Franck-Condon and an effective Hamiltonian analysis to interpret the Franck-Condon-like patterns. We examined vibrational wave functions obtained by the latter analysis. Shapes of the wave functions suggest that a recombination of the intermolecular vibrational modes occurs during the excitation of OH stretch mode in these clusters, which is a similar behavior to the Duschinsky effect in the electronic transition.

DFT Study on the Biosynthesis of Preasperterpenoid A: Role of Secondary Carbocations in the Carbocation Cascade

Preasperterpenoid A, featuring a 5/7/(3)6/5 pentacyclic structure, is a C₂₅ sesterterpenoid produced by Penicillium verruculosum. The results of density functional calculations on putative biosynthetic carbocation cyclization/rearrangements leading to preasperterpenoid A revealed a highly concerted four-step cyclization mechanism. Interestingly, two secondary carbocation structures were obtained as minima, but appeared almost as shoulders in the energy profile, and may represent essentially transient structures during the highly concerted reaction.

Spin-inversion mechanisms in O2 binding to a model heme complex revisited by density function theory calculations

Spin-inversion mechanisms in O₂ binding to a model heme complex, consisting of Fe(II)-porphyrin and imidazole, were investigated using density-functional theory calculations. First, we applied the recently proposed mixed-spin Hamiltonian method to locate spin-inversion structures between different total spin multiplicities. Nine spin-inversion structures were successfully optimized for the singlet-triplet, singlet-quintet, triplet-quintet, and quintet-septet spin-inversion processes. We found that the singlet-triplet spin-inversion points are located around the potential energy surface region at short Fe-O distances, whereas the singlet-quintet and quintet-septet spin-inversion points are located at longer Fe-O distances. This suggests that both narrow and broad crossing models play roles in O₂ binding to the Fe-porphyrin complex. To further understand spin-inversion mechanisms, we performed on-the-fly Born-Oppenheimer molecular dynamics calculations. The reaction coordinates, which are correlated to the spin-inversion dynamics between different spin multiplicities, are also discussed.

Migrations and Catalytic Action of Water Molecules in the Ionized Formamide-(H2O)2 Cluster

Isomerization dynamics involving the migrations, proton transfer reaction, and catalytic actions of water molecules upon vertical ionization of the formamide (FA)-(H₂O)₂ cluster is investigated by the infrared spectroscopy and theoretical reaction path search calculation. The infrared spectroscopic result indicates the [FA-(H₂O)₂]⁺ cation has the hydrogen-bonded structure of the enol isomer cation of formamide and the water dimer. This structure is formed by proton transfer from the CH bond to the carbonyl group through the catalytic action of the water molecules. The isomerization paths involving this enolization in ionized FA-(H₂O)₂ are explored by using the anharmonic downward distortion following method. We found multiple enolization paths which accompany proton exchanges among the formamide moiety and water molecules through the catalytic actions of the water molecules.

Toward strong self-healing polyisoprene elastomers with dynamic ionic crosslinks

To compromise high mechanical strength and efficient self-healing capability in an elastomer with dynamic crosslinks, optimization of the molecular structure is crucial in addition to the tuning of the dynamic properties of the crosslinks. Herein, we studied the effects of molecular weight, content of carboxy groups, and neutralization level of ionically crosslinked polyisoprene (PI) elastomers on their morphology, network rearrangement behavior, and self-healing and mechanical properties. In this PI elastomer, nanosized sphere-shaped ionic aggregates are formed by both neutralized and non-neutralized carboxy groups that act as stickers. The number density of the ionic aggregates that act as physical crosslinks increased with increase in the stickers' concentration, although the size of the ionic aggregates was independent of the molecular weight and the stickers' concentration. The ionic network was dynamically rearranged by the stickers' hopping between the ionic aggregates, and the rearrangement was accelerated by decreasing the neutralization level. We found that the 2Rg of the PI must be significantly larger than the average distance between the ionic aggregates to obtain a mechanically strong PI elastomer. We also found that further increase in the molecular weight is effective to enhance the dimensional stability of the elastomer. However, this approach reduced the elastomer's self-healing rate at the same time because the diffusion and randomization of the polymer chains between the damaged faces were reduced. In this work, we clearly demonstrated the principle in the optimization of the molecular structure for the ionically crosslinked PI elastomers to tune the mechanical and autonomous self-healing properties.

Conformational Landscapes and Infrared Spectra of Gas-phase Interstellar Molecular Clusters [(C3H3N)(CH3OH)n, n = 1-4]

Acrylonitrile (A) is one of the important interstellar molecules, which is considered closely related to the origin of life. And methanol (M) is one of the commonly used solvents, which is also found in outer space. Herein, we obtained the infrared (IR) spectra of size-selected AM_n (n = 1-4) clusters in supersonic jet by monitoring their fragments of H⁺AM_n-1 (n = 1-4) with vacuum ultraviolet single-photon soft ionization/IR-depletion technique. IR spectra of AM_n (n = 1-4) clusters were recorded in the CH and OH vibration bands in the range of 2700-3800 cm^-1. Spectra of AM_n (n = 1-4) clusters are similar in the CH stretching regions, while those show significant variations in the OH stretching regions with the increase of methanol molecules. Calculated IR spectra, which were predicted with the B3LYP-D3(BJ)/aug-cc-pVDZ method, were employed to compare with the experimental results. For AM, AM₂, and AM₃, the structures with the methanol cyclic hydrogen bonded with [N1-C4(H6)] of acrylonitrile are more stable than the other H-bonded structures. For the most stable structures of AM₄, however, the results show that the acrylonitrile is binding to a H-bonded ring formed by OH groups of four methanol molecules. The AM, AM₂, and AM₃ conformers with the single ring on the C₁ side of acrylonitrile, such as C₁-AM-a, C₁-AM₂-a, and C₁-AM₃-a, are dominant in the gas phase, while the C₂-AM₄-a conformer with the H-bonded ring formed by the OH groups on the C₂ side of acrylonitrile is more stable than that of CM₄-A-a in our experimental conditions (>130 K). These findings may provide valuable insight into the microsolvation process of the interstellar molecules and other biomolecules in gas phase.