Using bonding to guide transition state optimization

Locating Transition States by Variational Reaction Path Optimization with an Energy-Derivative-Free Objective Function

Locating transition states is essential for understanding molecular reactions. We propose a double-ended transition state search method by revisiting a variational reaction path optimization method known as the MaxFlux method. Although its original purpose is to add temperature effects to reaction paths, we conversely let the temperature approach zero to obtain an asymptotically exact minimum energy path and its corresponding transition state in variational formalism with an energy-derivative-free objective function. Using several numerical techniques to directly optimize the objective function, the present method reliably finds transition states with low computational cost. In particular, only three force evaluations per iteration are sufficient. This is confirmed on a variety of molecular reactions where the nudged elastic band method often fails. The present method is implemented in Python using the Atomic Simulation Environment and is available on GitHub.

Elastic Image Pair Method for Finding Transition States on Potential Energy Surfaces Using Only First Derivatives

Herein, an elastic image pair (EIP) method is proposed to search transition states between two potential-energy minima using only first derivatives. In this method, two images are generated, and the spring forces are added to the images to control the distance between the two images. Transition states are reached when the force and the distance of the image pair are both converged. A set of test molecules is optimized using the EIP method, which shows its efficiency in transition state searching compared to other methods. This new method is more stable and reliable in finding transition states with much less computations.

Evaluation of inhibitive corrosion potential of symmetrical hydrazine derivatives containing nitrophenyl moiety in 1M HCl for C38 steel: experimental and theoretical studies

The exploitation of cost-effective, sustainable, green and efficient compounds is a renewed science and a demanding mission for today's chemists and technologists. In this view, the inhibitive corrosion properties of some hydrazine derivatives named (1E,2E)-bis(1-(2-nitrophenyl)ethylidene)hydrazine (SSBO), (1E,2E)-bis(1-(3-nitrophenyl)ethylidene)hydrazine (SSBM) and (1E,2E)-bis(1-(4-nitrophenyl)ethylidene)hydrazine (SSBP) on the C38 steel corrosion in 1M HCl media has been evaluated by different techniques like electrochemical impedance spectroscopy (EIS), potentiodynamic polarization (PDP), scanning electron microscopy (SEM) and energy-dispersive X-ray spectroscopy. The EIS results showed that SSBM is the greatest inhibitor (η>93%) among the three tested compounds. The SSBM gives considerable inhibition efficiency against corrosion of steel compared to the previous studies. The PDP curves indicated that the studied inhibitors were a mixed-type inhibitor with a predominantly cathodic control. Quantum calculations of some descriptors derived from the density functional theory (DFT), the transition state theory (TST), the quantum theory of atoms in molecules (QTAIM) and molecular dynamics simulation have delivered helpful information regarding electron transfer and mechanism during adsorption of inhibitors on C38 steel surface.

Nudged Elastic Band Method for Molecular Reactions Using Energy-Weighted Springs Combined with Eigenvector Following

The climbing image nudged elastic band method (CI-NEB) is used to identify reaction coordinates and to find saddle points representing transition states of reactions. It can make efficient use of parallel computing as the calculations of the discretization points, the so-called images, can be carried out simultaneously. In typical implementations, the images are distributed evenly along the path by connecting adjacent images with equally stiff springs. However, for systems with a high degree of flexibility, this can lead to poor resolution near the saddle point. By making the spring constants increase with energy, the resolution near the saddle point is improved. To assess the performance of this energy-weighted CI-NEB method, calculations are carried out for a benchmark set of 121 molecular reactions. The performance of the method is analyzed with respect to the input parameters. Energy-weighted springs are found to greatly improve performance and result in successful location of the saddle points in less than a thousand energy and force evaluations on average (about a hundred per image) using the same set of parameter values for all of the reactions. Even better performance is obtained by stopping the calculation before full convergence and complete the saddle point search using an eigenvector following method starting from the location of the climbing image. This combination of methods, referred to as NEB-TS, turns out to be robust and highly efficient as it reduces the average number of energy and force evaluations down to a third, to 305. An efficient and flexible implementation of these methods has been made available in the ORCA software.

Simple position and orientation preconditioning scheme for minimum energy path calculations

Minimum-energy path (MEP) calculations, such as those typified by the nudged elastic band method, require input of reactant and product molecular configurations at initialization. In the case of reactions involving more than one molecule, generating initial reactant and product configurations requires careful consideration of the relative position and orientations of the reactive molecules in order to ensure that the resulting MEP calculation proceeds without converging on an alternative reaction-path, and without requiring excessive numbers of optimization iterations; as such, this initial system set-up is most commonly performed "by hand," with an expert user arranging reactive molecules in space to ensure that the following MEP calculation runs smoothly. In this Article, we introduce a simple preconditioning scheme which replaces this labor-intensive, human-knowledge-based step with an automated deterministic computational scheme. In our approach, initial reactant and product configurations are generated such that steric hindrance between reactive molecules is minimized in the reactant and product configurations, while also simultaneously requiring minimal structural differences between the reactants and products. The method is demonstrated using a benchmark test-set of >3400 organic molecular reactions, where comparison of the reactant/product configurations generated using our approach compare very well to initial configurations which were generated on an ad hoc basis.

Systematic Study of the Glutathione Reactivity of N-Phenylacrylamides: 2. Effects of Acrylamide Substitution

A comprehensive understanding of structure-reactivity relationships is critical to the design and optimization of cysteine-targeted covalent inhibitors. Herein, we report glutathione (GSH) reaction rates for N-phenyl acrylamides with varied substitutions at the α- and β-positions of the acrylamide moiety. We find that the GSH reaction rates can generally be understood in terms of the electron donating or withdrawing ability of the substituent. When installed at the β-position, aminomethyl substituents with amine pK_a's > 7 accelerate, while those with pK_a's < 7 slow the rate of GSH addition at pH 7.4, relative to a hydrogen substituent. Although a computational model was able to only approximately capture experimental reactivity trends, our calculations do not support a frequently invoked mechanism of concerted amine/thiol proton transfer and C-S bond formation and instead suggest that protonated aminomethyl functions as an electron-withdrawing group to reduce the barrier for thiolate addition to the acrylamide.

Modelling of an aza-Michael reaction from crystalline naphthalene derivatives containing peri–peri interactions: very long N–C bonds?

A correlation between N–C bond formation and CC bond breaking is constructed from the structures of a family of peri-naphthalenes with a second set of peri substituents.

Minimum Energy Paths and Transition States by Curve Optimization

Transition states and minimum energy paths are essential to understand and predict chemical reactivity. Double-ended methods represent a standard approach for their determination. We introduce a new double-ended method that optimizes reaction paths described by curves. Unlike other methods, our approach optimizes the curve parameters rather than distinct structures along the path. With molecular paths represented as continuous curves, the optimization can benefit from the advantages of an integral-based formulation. We call this approach ReaDuct and demonstrate its applicability for molecular paths parametrized by B-spline curves.

Reliable and efficient reaction path and transition state finding for surface reactions with the growing string method

The computational challenge of fast and reliable transition state and reaction path optimization requires new methodological strategies to maintain low cost, high accuracy, and systematic searching capabilities. The growing string method using internal coordinates has proven to be highly effective for the study of molecular, gas phase reactions, but difficulties in choosing a suitable coordinate system for periodic systems has prevented its use for surface chemistry. New developments are therefore needed, and presented herein, to handle surface reactions which include atoms with large coordination numbers that cannot be treated using standard internal coordinates. The double-ended and single-ended growing string methods are implemented using a hybrid coordinate system, then benchmarked for a test set of 43 elementary reactions occurring on surfaces. These results show that the growing string method is at least 45% faster than the widely used climbing image-nudged elastic band method, which also fails to converge in several of the test cases. Additionally, the surface growing string method has a unique single-ended search method which can move outward from an initial structure to find the intermediates, transition states, and reaction paths simultaneously. This powerful explorative feature of single ended-growing string method is demonstrated to uncover, for the first time, the mechanism for atomic layer deposition of TiN on Cu(111) surface. This reaction is found to proceed through multiple hydrogen-transfer and ligand-exchange events, while formation of H-bonds stabilizes intermediates of the reaction. Purging gaseous products out of the reaction environment is the driving force for these reactions. © 2017 Wiley Periodicals, Inc.

Path optimization by a variational reaction coordinate method. II. Improved computational efficiency through internal coordinates and surface interpolation

Reaction path optimization is being used more frequently as an alternative to the standard practice of locating a transition state and following the path downhill. The Variational Reaction Coordinate (VRC) method was proposed as an alternative to chain-of-states methods like nudged elastic band and string method. The VRC method represents the path using a linear expansion of continuous basis functions, allowing the path to be optimized variationally by updating the expansion coefficients to minimize the line integral of the potential energy gradient norm, referred to as the Variational Reaction Energy (VRE) of the path. When constraints are used to control the spacing of basis functions and to couple the minimization of the VRE with the optimization of one or more individual points along the path (representing transition states and intermediates), an approximate path as well as the converged geometries of transition states and intermediates along the path are determined in only a few iterations. This algorithmic efficiency comes at a high per-iteration cost due to numerical integration of the VRE derivatives. In the present work, methods for incorporating redundant internal coordinates and potential energy surface interpolation into the VRC method are described. With these methods, the per-iteration cost, in terms of the number of potential energy surface evaluations, of the VRC method is reduced while the high algorithmic efficiency is maintained.

The computational prediction of Raman and ROA spectra of charged histidine tautomers in aqueous solution

Histidine is a key component of a number of enzymatic mechanisms, and undertakes many functionalities in biochemical systems.