Dynamical Bifurcation in Gas‐Phase XH− + CH3Y SN2 Reactions: The Role of Energy Flow and Redistribution in Avoiding the Minimum Energy Path

The gas‐phase reactions of XH^? (X=O, S) + CH₃Y (Y=F, Cl, Br) span nearly the whole range of S_N2 pathways, and show an intrinsic reaction coordinate (IRC) (minimum energy path) with a deep well owing to the CH₃XH???Y^? (or CH₃S^????HF) hydrogen‐bonded postreaction complex. MP2 quasiclassical‐type direct dynamics starting at the [HX???CH₃???Y]^? transition‐state (TS) structure reveal distinct mechanistic behaviors. Trajectories that yield the separated CH₃XH+Y^? (or CH₃S^?+HF) products directly are non‐IRC, whereas those that sample the CH₃XH???Y^? (or CH₃S^????HF) complex are IRC. The IRCIRC/non‐IRC ratios of 90:10, 40:60, 25:75, 2:98, 0:100, and 0:100 are obtained for (X, Y)=(S, F), (O, F), (S, Cl), (S, Br), (O, Cl), and (O, Br), respectively. The properties of the energy profiles after the TS cannot provide a rationalization of these results. Analysis of the energy flow in dynamics shows that the trajectories cross a dynamical bifurcation, and that the inability to follow the minimum energy path arises from long vibration periods of the X?C???Y bending mode. The partition of the available energy to the products into vibrational, rotational, and translational energies reveals that if the vibrational contribution is more than 80 %, non‐IRC behavior dominates, unless the relative fraction of the rotational and translational components is similar, in which case a richer dynamical mechanism is shown, with an IRC/non‐IRC ratio that correlates to this relative fraction.     

RDX分子内振动能量重分配的动力学研究

基于从头算分子动力学(Born-oppenheimer molecular dynamics, BOMD)模拟, 构建了环硝胺六氢-1,3,5-三硝基-1,3,5-三嗪(RDX)单分子不同振动模式之间的耦合矩阵, 并计算了在不同加载能量下从低频振动模式到高频振动模式的最优能量传输路径. 结果表明, RDX单分子中—NNO₂基团更有利于能量局域化, 振动模式v₃和v₄在从低频振动模式到高频振动模式的能量传输过程中扮演着重要角色. 通过对v₃和v₄两个振动模式的进一步分析发现, 加载能量的不同会导致RDX单分子能量传输路径的不同. 当加载能量较低时, RDX单分子倾向于从低频振动模式到中频振动模式再到高频振动模式的能量传输路径; 当加载能量较高时, 能量更倾向于从低频振动模式直接传输到高频振动模式上. 揭示了RDX分子内振动耦合能量转移的微观机制, 为进一步探索RDX将“机械能”转化为“化学能”的微观过程提供了理论基础.     

Proton dissociation and transfer in hydrated phosphoric acid clusters

The dynamics and mechanisms of proton dissociation and transfer in hydrated phosphoric acid (H₃PO₄) clusters under excess proton conditions were studied based on the concept of presolvation using the H₃PO₄–H₃O⁺–nH₂O complexes (n = 1–3) as the model systems and ab initio calculations and Born–Oppenheimer molecular dynamics (BOMD) simulations at the RIMP2/TZVP level as model calculations. The static results showed that the smallest, most stable intermediate complex for proton dissociation (n = 1) is formed in a low local‐dielectric constant environment (e.g., ε = 1), whereas proton transfer from the first to the second hydration shell is driven by fluctuations in the number of water molecules in a high local‐dielectric constant environment (e.g., ε = 78) through the Zundel complex in a linear H‐bond chain (n = 3). The two‐dimensional potential energy surfaces (2D‐PES) of the intermediate complex (n = 1) suggested three characteristic vibrational and ¹H NMR frequencies associated with a proton moving on the oscillatory shuttling and structural diffusion paths, which can be used to monitor the dynamics of proton dissociation in the H‐bond clusters. The BOMD simulations over the temperature range of 298–430 K validated the proposed proton dissociation and transfer mechanisms by showing that good agreement between the theoretical and experimental data can be achieved with the proposed rate‐determining processes. The theoretical results suggest the roles played by the polar solvent and iterate that insights into the dynamics and mechanisms of proton transfer in the protonated H‐bond clusters can be obtained from intermediate complexes provided that an appropriate presolvation model is selected and that all of the important rate‐determining processes are included in the model calculations. © 2015 Wiley Periodicals, Inc.     

A conserved asparagine in a ubiquitin-conjugating enzyme positions the substrate for nucleophilic attack

The mechanism used by the ubiquitin-conjugating enzyme, Ubc13, to catalyze ubiquitination is probed with three computational techniques: Born–Oppenheimer molecular dynamics, single point quantum mechanics/molecular mechanics energies, and classical molecular dynamics. These simulations support a long-held hypothesis and show that Ubc13-catalyzed ubiquitination uses a stepwise, nucleophilic attack mechanism. Furthermore, they show that the first step—the formation of a tetrahedral, zwitterionic intermediate—is rate limiting. However, these simulations contradict another popular hypothesis that supposes that the negative charge on the intermediate is stabilized by a highly conserved asparagine (Asn79 in Ubc13). Instead, calculated reaction profiles of the N79A mutant illustrate how charge stabilization actually increases the barrier to product formation. Finally, an alternate role for Asn79 is suggested by simulations of wild-type, N79A, N79D, and H77A Ubc13: it stabilizes the motion of the electrophile prior to the reaction, positioning it for nucleophilic attack. © 2019 Wiley Periodicals, Inc.     

Multiscale electrostatic embedding simulations for modeling structure and dynamics of molecules in solution: A tutorial review

The main concepts and important technical details of electrostatic embedding quantum mechanics/molecular mechanics (QM/MM) simulations are explained and illustrated with the intent of assisting newcomers in performing and gauging the accuracy of such simulations, focused on smaller molecules in solution. Beginners are advised on how to increase the reliability and accuracy of the simulations through benchmarking. Central considerations on methodologies for QM/MM Molecular Dynamics (MD) simulations are presented, alongside technical fundamentals regarding the construction and manipulation of simulation systems using the python-based Atomic Simulation Environment (ASE). A worked example of QM/MM Born–Oppenheimer MD is included, and a flowchart summarizing the most salient decisions and tasks within the methodology is presented.     

Ab initio molecular dynamics study of the reactions of allene cation induced by intense 7 micron laser pulses

ABSTRACT

The isomerisation and fragmentation of allene cation (H₂C=C=CH₂⁺) by short, intense laser pulses were simulated by Born-Oppenheimer molecular dynamics (BOMD) on the ground state potential energy surface using the B3LYP/6-31?+?G(d,p) level of theory and a 10 cycle 7?µm cosine squared pulse with a maximum field strength of 0.07?au. Laser fields polarised along the C=C=C axis deposits an average of 150?kcal/mol in the molecule, compared to only 25 and 51?kcal/mol for perpendicular polarisations. Approximately 90% of the trajectories with the field aligned with the C=C=C axis underwent one or more structural rearrangement steps to form H₂C=CH–CH⁺ (15%), H₃CCCH⁺ (4%), cyclopropene cation (6%), and allene cation with rearranged hydrogens and carbons (47%). In addition, a variety of fragments including H₂CCCH⁺?+?H (10%), c-C₃H₃⁺?+?H (7%), and HCCCH⁺?+?H₂ (2%) trajectories were produced after isomerisation. With the same amount of thermal energy, field-free BOMD shows good agreements with the BOMD with the field. However, RRKM calculations favour isomerisation to propyne cation and dissociation to HCCCH⁺?+?H₂. This suggests that for molecules in intense laser fields the energy in the intermediate isomers is not distributed statistically.     

Dynamics of proton exchange in a model phosphonic acid‐functionalized polymer

The dynamics and mechanism of proton exchange in phosphonic acid‐functionalized polymers were studied using poly(vinyl‐phosphonic acid) (PVPA) as a model system along with quantum chemical calculations and Born–Oppenheimer molecular dynamics (BOMD) simulations at the B3LYP/TZVP level as model calculations. This theoretical study began with searching for the smallest, most active polymer segments and their intermediate conformations which could be involved in the local proton‐exchange process. The B3LYP/TZVP results confirmed that a low local dielectric environment and excess proton conditions are required to generate the intermediate conformations, and the shapes of the potential energy curves of the proton exchange between the two phosphonic acid functional groups are sensitive to the local conformational changes. In contrast, a high local dielectric environment increases the energy barriers, thereby preventing the proton from returning to the original functional group. Based on the static results, a mechanism for the proton exchange between the two functional groups involving fluctuations in the local dielectric environment and a local conformational change was proposed. The BOMD results confirmed the proposed mechanism by showing that the activation energies for the proton exchange in the hydrogen bond between two immobilized phosphonic acid moieties, obtained from the exponential relaxation behaviors of the envelopes of the velocity autocorrelation functions and the ¹H Nuclear Magnetic Resonance (NMR) line‐shape analyses, are too low to be the rate‐determining process. Instead, coupled librational motion in the backbone which leads to the interconversion between the two intermediate conformations possesses higher activation energy, and therefore represents one of the most important rate‐determining processes. These findings suggested that the rate of the proton exchange in the model phosphonic acid‐functionalized polymer is determined by the polymer mobility which, in this case, is the large‐amplitude librational motion of the vinyl backbone. © 2015 Wiley Periodicals, Inc.     

Measuring molecular force fields: Terahertz,inelastic neutron scattering,Raman, FTIR,DFT, and BOMD molecular dynamics of solid l-serine

A vibrational analysis of polycrystalline l-serine is provided using experimental terahertz, FTIR, Raman and inelastic neutron scattering (INS) spectra, calculated INS spectra – and Born–Oppenheimer molecular dynamics (BOMD) simulations from which the power spectra for the electronegative elements are compared to the THz spectra. Corrections are made to density functional theory (DFT) calculations for van der Waals interactions. Assignments and potential energy distributions are included for all 3N = 336 normal modes of an eight molecule supercell, including those for 48 non-bonded whole molecule translating and rotating vibrations, of which three are acoustic modes, usually not considered. Calculated and observed frequencies differ by an average 3 cm⁻¹ (s = 4). The INS spectrum of these modes below 100 cm⁻¹, calculated from energy second derivatives, show a remarkable similarity to the experimental 10 K spectra. The calculated low frequency modes are insensitive to small changes in cell parameters and geometry. THz intensities are represented by power spectra and not calculated explicitly. Nevertheless, power spectra of 13 ps BOMD trajectories at classical temperatures of 20 K, 400 K, and 500 K are markedly similar to the experimental terahertz spectra at 77 K and 298 K. Calculations on a serine crystal supercell 2 × 2 × 2 molecules deep appear to include, in a crude but fortuitously accurate way, enough of the principle out of phase dispersion to yield a match with experimental frequencies and intensities.     

Direct MP2 molecular dynamics studies of H atom reaction with CD4 and CH4

The mechanism and dynamics of the H + CD₄ → CD₃ + HD (I) and H + CH₄ → CH₃ + H₂ (II) reactions have been investigated by electronic structure methods. The minimum‐energy path and vibrational frequencies along the intrinsic reaction coordinate are calculated at MP2/cc‐pVDZ level. Energy distributions of the products are also obtained by the direct classical trajectory calculations at the MP2/ cc‐pVDZ level. It is found that most of the available energy appears as product translational energy, and very little of the available energy is partitioned into internal excitation of the HD (H₂) product for reaction I (II), which is in agreement with the experimental evidence. The results indicate that the experimental results could be reproduced by the direct MP2 molecular dynamics calculations. The rotational state distributions of the products show the HD (H₂) products are formed with lower rotational quantum numbers than the CD₃ (CH₃) products. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Inelastic neutron scattering,Raman, vibrational analysis with anharmonic corrections,and scaled quantum mechanical force field for polycrystalline l-alanine

A scaled quantum mechanical harmonic force field (SQMFF) corrected for anharmonicity is obtained for the 23 K l-alanine crystal structure using van der Waals corrected periodic boundary condition density functional theory (DFT) calculations with the PBE functional. Scale factors are obtained with comparisons to inelastic neutron scattering (INS), Raman, and FT-IR spectra of polycrystalline l-alanine at 15–23 K. Calculated frequencies for all 153 normal modes differ from observed frequencies with a standard deviation of 6 wavenumbers. Non-bonded external k = 0 lattice modes are included, but assignments to these modes are presently ambiguous. The extension of SQMFF methodology to lattice modes is new, as are the procedures used here for providing corrections for anharmonicity and van der Waals interactions in DFT calculations on crystals. First principles Born–Oppenheimer molecular dynamics (BOMD) calculations are performed on the l-alanine crystal structure at a series of classical temperatures ranging from 23 K to 600 K. Corrections for zero-point energy (ZPE) are estimated by finding the classical temperature that reproduces the mean square displacements (MSDs) measured from the diffraction data at 23 K. External k = 0 lattice motions are weakly coupled to bonded internal modes.