Small, solubilized platinum nanocrystals consist of an ordered core surrounded by mobile surface atoms.

In situ structures of Platinum (Pt) nanoparticles (NPs) can be determined with graphene liquid cell transmission electron microscopy. Atomic-scale three-dimensional structural information about their physiochemical properties in solution is critical for understanding their chemical function. We here analyze eight atomic-resolution maps of small (<3 nm) colloidal Pt NPs. Their structures are composed of an ordered crystalline core surrounded by surface atoms with comparatively high mobility. 3D reconstructions calculated from cumulative doses of 8500 and 17,000 electrons/pixel, respectively, are characterized in terms of loss of atomic densities and atomic displacements. Less than 5% of the total number of atoms are lost due to dissolution or knock-on damage in five of the structures analyzed, whereas 10-16% are lost in the remaining three. Less than 5% of the atomic positions are displaced due to the increased electron irradiation in all structures. The surface dynamics will play a critical role in the diverse catalytic function of Pt NPs and must be considered in efforts to model Pt NP function computationally.

Obtaining QM/MM binding free energies in the SAMPL8 drugs of abuse challenge: indirect approaches.

Accurately predicting free energy differences is essential in realizing the full potential of rational drug design. Unfortunately, high levels of accuracy often require computationally expensive QM/MM Hamiltonians. Fortuitously, the cost of employing QM/MM approaches in rigorous free energy simulation can be reduced through the use of the so-called "indirect" approach to QM/MM free energies, in which the need for QM/MM simulations is avoided via a QM/MM "correction" at the classical endpoints of interest. Herein, we focus on the computation of QM/MM binding free energies in the context of the SAMPL8 Drugs of Abuse host-guest challenge. Of the 5 QM/MM correction coupled with force-matching submissions, PM6-D3H4/MM ranked submission proved the best overall QM/MM entry, with an RMSE from experimental results of 2.43 kcal/mol (best in ranked submissions), a Pearson's correlation of 0.78 (second-best in ranked submissions), and a Kendall [Formula: see text] correlation of 0.52 (best in ranked submissions).

The Extended Eighth-Shell method for periodic boundary conditions with rotational symmetry.

The Eighth-Shell method for parallelization of molecular dynamics simulations has previously been shown to be the most optimal for efficiency at large process counts. However, in its current formulation only the P1 space group is supported for periodic boundary conditions (PBC) and thus reflection and/or rotational crystal symmetries are not supported. In this work, we outline the development and implementation of the Extended Eighth-Shell (EES) method that allows rotational symmetry by using an extended import region compared to the ES method. It simulates only the asymmetric unit and communicates coordinates and forces with images that correspond to P21 PBC. The P21 PBC has application in lipid bilayer simulations as it can be used to allow lipids to switch leaftlets, thus rapidly balancing the chemical potential difference between the two layers. Our results show that the EES method scales efficiently over large number of processes and can be used for simulations with P21 symmetry in an orthorhombic crystal.

A replica exchange umbrella sampling (REUS) approach to predict host-guest binding free energies in SAMPL8 challenge.

In this study, we report binding free energy calculations of various drugs-of-abuse to Cucurbit-[8]-uril as part of the SAMPL8 blind challenge. Force-field parameters were obtained from force-matching with different quantum mechanical levels of theory. The Replica Exchange Umbrella Sampling (REUS) approach was used with a cylindrical restraint to enhance the sampling of host-guest binding. Binding free energy was calculated by pulling the guest molecule from one side of the symmetric and cylindrical host, then into and through the host, and out the other side (bidirectional) as compared to pulling only to the bound pose inside the cylindrical host (unidirectional). The initial results with force-matched MP2 parameter set led to RMSE of 4.68 [Formula: see text] from experimental values. However, the follow-up study with CHARMM generalized force field parameters and force-matched PM6-D3H4 parameters resulted in RMSEs from experiment of [Formula: see text] and [Formula: see text], respectively, which demonstrates the potential of REUS for accurate binding free energy calculation given a more suitable description of energetics. Moreover, we compared the free energies for the so called bidirectional and unidirectional free energy approach and found that the binding free energies were highly similar. However, one issue in the bidirectional approach is the asymmetry of profile on the two sides of the host. This is mainly due to the insufficient sampling for these larger systems and can be avoided by longer sampling simulations. Overall REUS shows great promise for binding free energy calculations.

Direct diabatization and analytic representation of coupled potential energy surfaces and couplings for the reactive quenching of the excited 2Σ+ state of OH by molecular hydrogen.

We have employed extended multiconfiguration quasidegenerate perturbation theory, fourfold-way diabatic molecular orbitals, and configurational uniformity to develop a global three-state diabatic representation of the potential energy surfaces and their couplings for the electronically nonadiabatic reaction OH* + H2 → H2O + H, where * denotes electronic excitation to the A 2Σ+ state. To achieve sign consistency of the computed diabatic couplings, we developed a graphics processing unit-accelerated algorithm called the cluster-growing algorithm. Having obtained consistent signs of the diabatic couplings, we fit the diabatic matrix elements (which consist of the diabatic potentials and the diabatic couplings) to analytic representations. Adiabatic potential energy surfaces are generated by diagonalizing the 3 × 3 diabatic potential energy matrix. The comparisons between the fitted and computed diabatic matrix elements and between the originally computed adiabatic potential energy surfaces and those generated from the fits indicate that the current fit is accurate enough for dynamical studies, and it may be used for quantal or semiclassical dynamics calculations.

Anharmonicity of Coupled Torsions: The Extended Two-Dimensional Torsion Method and Its Use To Assess More Approximate Methods.

In this work we present the extended two-dimensional torsion (E2DT) method and use it to analyze the performance of several methods that incorporate torsional anharmonicity more approximately for calculating rotational-vibrational partition functions. Twenty molecules having two hindered rotors were studied for temperatures between 100 and 2500 K. These molecules present several kinds of situations; they include molecules with nearly separable rotors, molecules in which the reduced moments of inertia change substantially with the internal rotation, and molecules presenting compound rotation. Partition functions obtained by the rigid-rotor harmonic oscillator approximation, a method involving global separability of torsions and the multistructural methods without explicit potential coupling [MS-T(U)] and with explicit potential coupling [MS-T(C)] of torsions, are compared to those obtained with a quantized version - called the extended two-dimensional torsion (E2DT) method - of the extended hindered rotor approximation of Vansteenkiste et al. ( Vansteenkiste et al. J. Chem. Phys. 2006 , 124 , 044314 ). In the E2DT method, quantum effects due to the torsional modes were incorporated by the two-dimensional nonseparable method, which is a method that is based on the solution of the torsional Schrödinger equation and that includes full coupling in both the kinetic and potential energy. By comparing other methods to the E2DT method and to experimental thermochemical data, this study concludes that the harmonic approximation yields very poor results at high temperatures; the global separation of torsions from the rest of the degrees of freedom is not justified even when an accurate method to treat the torsions is employed; it is confirmed that methods based on less complete potential energy coupling of torsions, such as MS-T(U), are not accurate when dealing with rotors with different barrier heights, and more complete inclusion of torsional coupling to the method in MS-T(C) improves substantially the results in such a way that it could be used in cases where the E2DT method is unaffordable.

Computational Kinetics by Variational Transition-State Theory with Semiclassical Multidimensional Tunneling: Direct Dynamics Rate Constants for the Abstraction of H from CH3OH by Triplet Oxygen Atoms.

Rate constants and the product branching ratio for hydrogen abstraction from CH3OH by O(3P) were computed with multistructural variational transition-state theory including microcanonically optimized multidimensional tunneling. Benchmark calculations of the forward and reverse classical barrier heights and the reaction energetics have been carried out by using coupled cluster theory and multireference calculations to select the most reliable density functional method for direct dynamics computations of the rate constants. The dynamics calculations included the anharmonicity of the zero-point energies and partition functions, with specific-reaction-parameter scaling factors for reactants and transition states, and multistructural torsional anharmonicity was included for the torsion around the C-O bond in methanol and in the transition states. The resulting rate constants are presented over a wider range than they are available from experiment, but in the temperature range where experiments are available, they agree well with experimental values, which is encouraging for their reliability over the wider temperature range and for future computations of oxygen atom reaction rates. In contrast to a previous computational prediction, the branching ratio predicted by the present work shows that the formation of CH2OH + OH is the dominant channel over the whole range of temperature from 250 to 2000 K.

A quasiclassical trajectory study of the N2(X(1)Σ) + O((3)P) → NO(X(2)Π) + N((4)S) reaction.

We present quasiclassical trajectory simulations of the high-energy collision-induced O((3)P) + N2 reaction. The simulation results show that reactive scattering may proceed on both the (3)A' and (3)A″ potential energy surfaces, but the reaction on (3)A″ surface overwhelmingly dominates. At total energies of collision slightly above the threshold, the NO(X(2)Π) product is found to be scattered primarily in the forward direction, and at higher collision energies, the fraction of NO products that scatter in the backward direction increases. Two-dimensional joint distributions of the center-of-mass differential cross section and translational energy have been examined to unravel the details of the scattering as functions of collision energy.

An Ancient Fingerprint Indicates the Common Ancestry of Rossmann-Fold Enzymes Utilizing Different Ribose-Based Cofactors.

Nucleoside-based cofactors are presumed to have preceded proteins. The Rossmann fold is one of the most ancient and functionally diverse protein folds, and most Rossmann enzymes utilize nucleoside-based cofactors. We analyzed an omnipresent Rossmann ribose-binding interaction: a carboxylate side chain at the tip of the second ß-strand (ß2-Asp/Glu). We identified a canonical motif, defined by the ß2-topology and unique geometry. The latter relates to the interaction being bidentate (both ribose hydroxyls interacting with the carboxylate oxygens), to the angle between the carboxylate and the ribose, and to the ribose's ring configuration. We found that this canonical motif exhibits hallmarks of divergence rather than convergence. It is uniquely found in Rossmann enzymes that use different cofactors, primarily SAM (S-adenosyl methionine), NAD (nicotinamide adenine dinucleotide), and FAD (flavin adenine dinucleotide). Ribose-carboxylate bidentate interactions in other folds are not only rare but also have a different topology and geometry. We further show that the canonical geometry is not dictated by a physical constraint--geometries found in noncanonical interactions have similar calculated bond energies. Overall, these data indicate the divergence of several major Rossmann-fold enzyme classes, with different cofactors and catalytic chemistries, from a common pre-LUCA (last universal common ancestor) ancestor that possessed the ß2-Asp/Glu motif.

Potential energy surface of triplet N2O2.

We present a global ground-state triplet potential energy surface for the N2O2 system that is suitable for treating high-energy vibrational-rotational energy transfer and collision-induced dissociation. The surface is based on multi-state complete-active-space second-order perturbation theory/minimally augmented correlation-consistent polarized valence triple-zeta electronic structure calculations plus dynamically scaled external correlation. In the multireference calculations, the active space has 14 electrons in 12 orbitals. The calculations cover nine arrangements corresponding to dissociative diatom-diatom collisions of N2, O2, and nitric oxide (NO), the interaction of a triatomic molecule (N2O and NO2) with the fourth atom, and the interaction of a diatomic molecule with a single atom (i.e., the triatomic subsystems). The global ground-state potential energy surface was obtained by fitting the many-body interaction to 54 889 electronic structure data points with a fitting function that is a permutationally invariant polynomial in terms of bond-order functions of the six interatomic distances.

Multi-path variational transition state theory for chiral molecules: the site-dependent kinetics for abstraction of hydrogen from 2-butanol by hydroperoxyl radical, analysis of hydrogen bonding in the transition state, and dramatic temperature dependence of the activation energy.

The goal of the present work is modeling the kinetics of a key reaction involved in the combustion of the biofuel 2-butanol. To accomplish this we extended multi-path variational transition state theory (MP-VTST) with the small curvature tunneling (SCT) approximation to include multistructural anharmonicity factors for molecules with chiral carbons. We use the resulting theory to predict the site-dependent rate constants of the hydrogen abstraction from 2-butanol by hydroperoxyl radical. The generalized transmission coefficients were averaged over the four lowest-energy reaction paths. The computed forward reaction rate constants indicate that hydrogen abstraction from the C-2 site has the largest contribution to the overall reaction from 200 K to 2400 K, with a contribution ranging from 99.9988% at 200 K to 88.9% at 800 K to 21.2% at 3000 K, while hydrogen abstraction from the oxygen site makes the lowest contribution at all temperatures, ranging from 2.5 × 10-9% at 200 K to 0.65% at 800 K to 18% at 3000 K. This work highlights the importance of including the multiple-structure and torsional potential anharmonicity in the computation of the thermal rate constants. We also analyzed the role played by the hydrogen bond at the transition state, and we illustrated the risks of (a) considering only the lowest-energy conformations in the calculations of the rate constants or (b) ignoring the nonlinear temperature dependence of the activation energies. A hydrogen bond at the transition state can lower the enthalpy of activation, but raise the free energy of activation. We find an energy of activation that increases from 11 kcal mol-1 at 200 K to more than 36 kcal mol-1 at high temperature for this radical reaction with a biofuel molecule.

Accounting for conformational flexibility and torsional anharmonicity in the H + CH3CH2OH hydrogen abstraction reactions: a multi-path variational transition state theory study.

This work reports a detailed theoretical study of the hydrogen abstraction reactions from ethanol by atomic hydrogen. The calculated thermal rate constants take into account torsional anharmonicity and conformational flexibility, in addition to the variational and tunneling effects. Specifically, the kinetics calculations were performed by using multi-path canonical variational transition state theory with least-action path tunneling corrections, to which we have added the two-dimensional non-separable method to take into account torsional anharmonicity. The multi-path thermal rate constant is expressed as a sum over conformational reaction channels. Each of these channels includes all the transition states that can be reached by internal rotations. The results show that, in the interval of temperatures between 250 and 2500 K, the account for multiple paths leads to higher thermal rate constants with respect to the single path approach, mainly at low and at high temperatures. In addition, torsional anharmonicity enhances the slope of the Arrhenius plot in this range of temperatures. Finally, we show that the incorporation of tunneling into the hydrogen abstraction reactions substantially changes the contribution of each of the transition states to the conformational reaction channel.

A convergent approach to the dioxaadamantane core of (±)-tetrodotoxin.

A fully stereocontrolled 1,3-diol orthoesterification and a water-promoted intramolecular Henry addition, combined with the previously reported formal (3 + 3) annulation of α-nitro-α,ß-enals and 2,2-dimethyl-1,3-dioxan-5-one, provided for a short convergent pathway to the dioxaadamantane core of (±)-tetrodotoxin.

Prediction of experimentally unavailable product branching ratios for biofuel combustion: the role of anharmonicity in the reaction of isobutanol with OH.

Isobutanol is a prototype biofuel, and sorting out the mechanism of its combustion is an important objective where theoretical modeling can provide information that is unavailable and not easily obtained by experiment. In the present work the rate constants and branching ratios for the hydrogen abstraction reactions from isobutanol by hydroxyl radical have been calculated using multi-path variational transition-state theory with small-curvature tunneling. We use hybrid degeneracy-corrected vibrational perturbation theory to show that it is critical to consider the anharmonicity difference of high-frequency modes between reactants and transition states. To obtain accurate rate constants, we must apply different scaling factors to the calculated harmonic vibrational frequencies at the reactants and at the transition states. The factors determining the reaction rate constants have been analyzed in detail, including variational effects, tunneling contributions, the effect of multiple reaction paths on transmission coefficients, and anharmonicities of low- and high-frequency vibrational modes. The analysis quantifies the uncertainties in the rate calculations. A key result of the paper is a prediction for the site dependence of hydrogen abstraction from isobutanol by hydroxyl radical. This is very hard to measure experimentally, although it is critical for combustion mechanism modeling. The present prediction differs considerably from previous theoretical work.

Including Tunneling in Non-Born-Oppenheimer Simulations.

For electronically nonadiabatic processes in all but the simplest systems, the most practical multidimensional simulation method is a semiclassical approximation in which a trajectory or the center of a wave packet follows a classical path governed by an effective potential energy function. Here, we show how such simulations can be made more realistic by including tunneling by the army ants tunneling method. We illustrate the theory by calculations with model potential energy surfaces; one model study is in the adiabatic limit, and the other one has nonadiabatic transitions between two electronic states during the tunneling event. The army ants tunneling algorithm is used to efficiently sample tunneling events in the trajectories in both cases. This work makes it possible to simulate complex nonadiabatic chemical processes by efficiently including the important quantum effect of tunneling.

New pathways for formation of acids and carbonyl products in low-temperature oxidation: the Korcek decomposition of γ-ketohydroperoxides.

We present new reaction pathways relevant to low-temperature oxidation in gaseous and condensed phases. The new pathways originate from γ-ketohydroperoxides (KHP), which are well-known products in low-temperature oxidation and are assumed to react only via homolytic O-O dissociation in existing kinetic models. Our ab initio calculations identify new exothermic reactions of KHP forming a cyclic peroxide isomer, which decomposes via novel concerted reactions into carbonyl and carboxylic acid products. Geometries and frequencies of all stationary points are obtained using the M06-2X/MG3S DFT model chemistry, and energies are refined using RCCSD(T)-F12a/cc-pVTZ-F12 single-point calculations. Thermal rate coefficients are computed using variational transition-state theory (VTST) calculations with multidimensional tunneling contributions based on small-curvature tunneling (SCT). These are combined with multistructural partition functions (Q(MS-T)) to obtain direct dynamics multipath (MP-VTST/SCT) gas-phase rate coefficients. For comparison with liquid-phase measurements, solvent effects are included using continuum dielectric solvation models. The predicted rate coefficients are found to be in excellent agreement with experiment when due consideration is made for acid-catalyzed isomerization. This work provides theoretical confirmation of the 30-year-old hypothesis of Korcek and co-workers that KHPs are precursors to carboxylic acid formation, resolving an open problem in the kinetics of liquid-phase autoxidation. The significance of the new pathways in atmospheric chemistry, low-temperature combustion, and oxidation of biological lipids are discussed.

Chloroform as a hydrogen atom donor in Barton reductive decarboxylation reactions.

The utility of chloroform as both a solvent and a hydrogen atom donor in Barton reductive decarboxylation of a range of carboxylic acids was recently demonstrated (Ko, E. J. et al. Org. Lett. 2011, 13, 1944). In the present work, a combination of electronic structure calculations, direct dynamics calculations, and experimental studies was carried out to investigate how chloroform acts as a hydrogen atom donor in Barton reductive decarboxylations and to determine the scope of this process. The results from this study show that hydrogen atom transfer from chloroform occurs directly under kinetic control and is aided by a combination of polar effects and quantum mechanical tunneling. Chloroform acts as an effective hydrogen atom donor for primary, secondary, and tertiary alkyl radicals, although significant chlorination was also observed with unstrained tertiary carboxylic acids.

Tunneling and conformational flexibility play critical roles in the isomerization mechanism of vitamin D.

The thermal isomerization reaction converting previtamin D to vitamin D is an intramolecular [1,7]-sigmatropic hydrogen shift with antarafacial stereochemistry. We have studied the dynamics of this reaction by means of the variational transition-state theory with multidimensional corrections for tunneling in both gas-phase and n-hexane environments. Two issues that may have important effects on the dynamics were analyzed in depth, i.e., the conformations of previtamin D and the quantum effects associated with the hydrogen-transfer reaction. Of the large number of conformers of previtamin D that were located, there are 16 that have the right disposition to react. The transition-state structures associated with these reaction paths are very close in energy, so all of them should be taken into account for an accurate calculation of both the thermal rate constants and the kinetic isotope effects. This issue is particularly important because the contribution of each of the reaction paths to the total thermal rate constant is quite sensitive to the environment. The dynamics results confirm that tunneling plays an important role and that model systems that were considered previously to study the hydrogen shift reaction cannot mimic the complexity introduced by the flexibility of the rings of previtamin D. Finally, the characterization of the conformers of both previtamin D and vitamin D allowed the calculation of the thermal equilibrium constants of the isomerization process.

High-level direct-dynamics variational transition state theory calculations including multidimensional tunneling of the thermal rate constants, branching ratios, and kinetic isotope effects of the hydrogen abstraction reactions from methanol by atomic hydrogen.

We report a detailed theoretical study of the hydrogen abstraction reaction from methanol by atomic hydrogen. The study includes the analysis of thermal rate constants, branching ratios, and kinetic isotope effects. Specifically, we have performed high-level computations at the MC3BB level together with direct dynamics calculations by canonical variational transition state theory (CVT) with the microcanonically optimized multidimensional tunneling (µOMT) transmission coefficient (CVT/µOMT) to study both the CH(3)OH+H→CH(2)OH+H(2) (R1) reaction and the CH(3)OH+H→CH(3)O+H(2) (R2) reaction. The CVT/µOMT calculations show that reaction R1 dominates in the whole range 298≤T (K)≤2500 and that anharmonic effects on the torsional mode about the C-O bond are important, mainly at high temperatures. The activation energy for the total reaction sum of R1 and R2 reactions changes substantially with temperature and, therefore, the use of straight-line Arrhenius plots is not valid. We recommend the use of new expressions for the total R1 + R2 reaction and for the R1 and R2 individual reactions.