FMO-MD simulations on the hydration of formaldehyde in water solution with constraint dynamics

Full-quantum mechanical fragment molecular orbital-based molecular dynamics (FMO-MD) simulations were applied to the hydration reaction of formaldehyde in water solution under neutral conditions. Two mechanisms, a concerted and a stepwise one, were considered with respect to the nucleophilic addition and the proton transfer. Preliminary molecular orbital calculations by means of polarized continuum model reaction field predicted that the hydration prefers a concerted mechanism. Because the calculated activation barriers were too high for free FMO-MD simulations to give reactive trajectories spontaneously, a More O'Ferrall-Jencks-type diagram was constructed from the statistical analysis of the FMO-MD simulations with constraint dynamics. The diagram showed that the hydration proceeds through a zwitterionic-like (ZW-like) structure. The free energy changes along the reaction coordinate calculated by means of the blue moon ensemble for the hydration and the amination of formaldehyde indicated that the hydration proceeds through a concerted process through the ZW-like structure, whereas the amination goes through a stepwise mechanism with a ZW intermediate. In inspection of the FMO-MD trajectories, water-mediated cyclic proton transfers were observed in both reactions on the way from the ZW-like structure to the product. These proton transfers also have an asynchronous character, in which deprotonation from the nucleophilic oxygen atom (or nitrogen atom for amination) precedes the protonation of the carbonyl oxygen atom. The results showed the strong advantage of the FMO-MD simulations to obtain detailed information at a molecular level for solution reactions.     

Vibrational overtone induced elimination reactions within hydrogen-bonded molecular clusters: the dynamics of water catalyzed reactions in CH2FOH.H2On

The dynamics of the light initiated OH-overtone induced elimination reactions CH(2)FOH.(H(2)O)(n) + hnu--> HF + CH(2)O + n(H(2)O), n = 1-3, are studied using classical trajectory simulations where the ab initio potential energy surface is computed "on-the-fly". Hydrogen bonding to the water is found to lower the barrier to reaction by over 20 kcal mol(-1) and modifies the mechanism to a concerted multiple H-atom transfer process. The reaction process is found to occur on a rapid timescale, <100 fs, and involves the hydronium ion as an intermediate. An essential aspect of dynamics is the successful competition of reaction with energy dissipation through water evaporation from the cluster.     

FMO‐MD Simulations on the Hydration of Formaldehyde in Water Solution with Constraint Dynamics

Full‐quantum mechanical fragment molecular orbital‐based molecular dynamics (FMO‐MD) simulations were applied to the hydration reaction of formaldehyde in water solution under neutral conditions. Two mechanisms, a concerted and a stepwise one, were considered with respect to the nucleophilic addition and the proton transfer. Preliminary molecular orbital calculations by means of polarized continuum model reaction field predicted that the hydration prefers a concerted mechanism. Because the calculated activation barriers were too high for free FMO‐MD simulations to give reactive trajectories spontaneously, a More O’Ferrall–Jencks‐type diagram was constructed from the statistical analysis of the FMO‐MD simulations with constraint dynamics. The diagram showed that the hydration proceeds through a zwitterionic‐like (ZW‐like) structure. The free energy changes along the reaction coordinate calculated by means of the blue moon ensemble for the hydration and the amination of formaldehyde indicated that the hydration proceeds through a concerted process through the ZW‐like structure, whereas the amination goes through a stepwise mechanism with a ZW intermediate. In inspection of the FMO‐MD trajectories, water‐mediated cyclic proton transfers were observed in both reactions on the way from the ZW‐like structure to the product. These proton transfers also have an asynchronous character, in which deprotonation from the nucleophilic oxygen atom (or nitrogen atom for amination) precedes the protonation of the carbonyl oxygen atom. The results showed the strong advantage of the FMO‐MD simulations to obtain detailed information at a molecular level for solution reactions.     

"Concerted" transition state, stepwise mechanism. Dynamics effects in C2-C6 enyne allene cyclizations

The C2-C6 (Schmittel)/ene cyclization of enyne-allenes is studied by a combination of kinetic isotope effects, theoretical calculations, and dynamics trajectories. For the cyclization of allenol acetate 9, the isotope effect (k(CH3)/k(CD3) is approximately 1.43. The isotope effect is interpreted in terms of a highly asynchronous transition state near the concerted/stepwise boundary. This is supported by density functional theory calculations that locate a highly asynchronous transition structure for the concerted ene reaction. However, calculations of both the experimental system and a model reaction were unable to locate a transition structure for formation of the diradical intermediate of a stepwise mechanism. The stepwise mechanism and the asynchronous concerted mechanism start out geometrically similar, and the two pathways appear to have merged as far as the initial transition structure. For the model reaction, quasiclassical direct dynamics trajectories emanating from the initial transition structure afforded the diradical intermediate in 29 out of 101 trajectories. A large portion of the remaining trajectories completes hydrogen transfer before carbon-carbon bond formation, despite the advanced carbon-carbon bond formation in the asynchronous transition structure. Overall, the single minimum-energy path from starting material to product is inadequate to describe the reaction, and a consideration of dynamic effects is necessary to understand the mechanism. The implications of these observations toward questions of concert in other reactions are discussed.     

Three-Body Photodissociation of Thionyl Chloride

The ultrafast photodissociation dynamics of thionyl chloride after excitation with a 235 nm pump pulse has been studied using a femtosecond time-resolved mass spectroscopy. The observed parent transient suggests that the excited states initiated by 235 nm are very shortlived, that is, approximately 166 fs. Not only a stepwise three-body dissociation pathway but also a concerted three-body dissociation pathway is observed.     

Competition Between Concerted and Stepwise Dynamics in the Triplet Di‐π‐Methane Rearrangement

The molecular dynamics of the triplet‐state Zimmerman di‐π‐methane rearrangement of dibenzobarrelene were computed with B3LYP and M06‐2X density functionals. All productive quasiclassical trajectories involve sequential formation and cleavage of C? C bonds and an intermediate with lifetimes ranging from 13 to 1160 fs. Both dynamically concerted and stepwise trajectories are found. The average lifetime of this intermediate is significantly shorter than predicted by either transition‐state theory or the Rice–Ramsperger–Kassel–Marcus model, thus indicating the non‐statistical nature of the reaction mechanism.     

Electron and nuclear dynamics of molecular clusters in ultraintense laser fields. II. Electron dynamics and outer ionization of the nanoplasma

We explore electron dynamics in molecular (CD4)(1061) clusters and elemental Xen (n=249-2171) clusters, responding to ultraintense (intensity I=10(16)-10(19) W cm(-2)) laser fields. Molecular dynamics simulations (including magnetic field and relativistic effects) and analyses of high-energy electron dynamics and nuclear ion dynamics in a cluster interacting with a Gaussian shaped laser field (frequency 0.35 fs(-1), photon energy 1.44 eV, phase 0, temporal width 25 fs) elucidated the time dependence of inner ionization, the formation of a nanoplasma of unbound electrons within the cluster or its vicinity, and of outer ionization. We determined the cluster size and the laser intensity dependence of these three sequential-parallel electronic processes. The characteristic times for cluster inner ionization (tau(ii)) and for outer ionization (tau(oi)) fall in the femtosecond time domain, i.e., tau(ii)=2-9 fs and tau(oi)=4-15 fs for (CD4)(1061), tau(ii)=7-30 fs and tau(oi)=5-13 fs for Xe(n) (n=479,1061), with both tau(ii) and tau(oi) decreasing with increasing I, in accord with the barrier suppression ionization mechanism for inner ionization of the constituents and the cluster barrier suppression ionization mechanism for outer ionization. The positive delay times Deltatau(OI) between outer and inner ionization (e.g., Deltatau(OI)=6.5 fs for Xen at I=10(16) W cm(-2) and Deltatau(OI)=0.2 fs for (CD4)(1061) at I=10(19) W cm(-2)) demonstrate that the outer/inner ionization processes are sequential. For (CD4)(1061), tau(ii)tau(oi), reflecting on the energetic hierarchy in the ionization of the Xe atoms. Quasiresonance contributions to the outer ionization of the nanoplasma were established, as manifested in the temporal oscillations in the inner/outer ionization levels, and in the center of mass of the nanoplasma electrons. The formation characteristics, dynamics, and response of the nanoplasma in molecular or elemental clusters were addressed. The nanoplasma is positively charged, with a high-average electron density [rho(P)=(2-3)10(22) cm(-3)], being characterized by high-average electron energies epsilon(av) (e.g., in Xe(1061) clusters epsilon(av)=54 eV at I=10(16) W cm(-2) and epsilon(av)=0.56-0.37 keV at I=10(18) W cm(-2), with epsilon(av) proportional, variant I(1/2)). Beyond the cluster boundary the average electron energy markedly increases, reaching electron energies in the range of 1.2-40 keV for outer ionization of Xe(n) (n=249-2171) clusters. The nanoplasma exhibits spatial inhomogeneity and angular anisotropy induced by the laser field. Femtosecond time scales are predicted for the nanoplasma production (rise times 7-3 fs), for the decay (decay times approximately 5 fs), and for the persistence time (30-10 fs) of a transient nanoplasma at I=10(17)-10(18) W cm(-2). At lower intensities of I=10(16) W cm(-2) a persistent nanoplasma with a "long" lifetime of > 50 fs will prevail.     

Multiple time scales in solvation dynamics of DNA in aqueous solution: the role of water, counterions, and cross-correlations

Recent time domain experiments have explored solvation dynamics of a probe located inside a DNA duplex, in an effort to gain information, e.g., on the dynamics of water molecules in the DNA major and minor grooves and their environment. Multiple time constants in the range of a few picoseconds to several nanoseconds were obtained. We have carried out 15 ns long atomistic molecular dynamics simulations to study the solvation dynamics of bases of a 38 base-pair long DNA duplex in an aqueous solution containing counterions. We have computed the energy-energy time correlation function (TCF) of the four individual bases (A, T, G, and C) to characterize the solvation dynamics. All the TCFs display highly nonexponential decay with time. When the trajectories are analyzed with 100 fs time resolution, the TCF of each base shows initial ultrafast decay (with tau1 approximately equal 60-80 fs) followed by two intermediate components (tau2 approximately equal 1 ps, tau3 approximately equal 20-30 ps), in near complete agreement with a recent time domain experiment on DNA solvation. Interestingly, the solvation dynamics of each of the four different nucleotide bases exhibit rather similar time scales. To explore the existence of slow relaxation at longer times reported recently in a series of experiments, we also analyzed the solvation TCFs calculated with longer time trajectories and with a larger time resolution of 1 ps. In this case, an additional slow component with a time constant of the order of 250 ps is observed. Through an analysis of partial solvation TCFs, we find that the slow decay originates mainly from the interaction of the nucleotides with the dipolar water molecules and the counterions. An interesting negative cross-correlation between water and counterions is observed, which makes an important contribution to relaxation at intermediate to longer times.     

Molecular Dynamics Insights into Electron-Catalyzed Dissociation Repair of Cyclobutane Pyrimidine Dimer

Excess electrons are not only an important source of radiation damage, but also participate in the repair process of radiation damage such as cyclobutane pyrimidine dimer (CPD). Using ab initio molecular dynamics (AIMD) simulations, we reproduce the single excess electron stepwise catalytic CPD dissociation process in detail with an emphasis on the energy levels and molecular structure details associated with excess electrons. On the basis of the AIMD simulations on the CPD aqueous solution with two vertically added excess electrons, we exclude the early-proposed [2+2]-like concerted synchronous dissociation mechanism, and analyze the difference between the symmetry of the actual reaction and the symmetry of the frontier molecular orbitals which deeply impact the mechanism. Importantly, we propose a new model of the stepwise electron-catalyzed dissociation mechanism that conforms to the reality. This work not only provides dynamics insights into the excess electron catalyzed dissociation mechanism, but also reveals different roles of two excess electrons in two bond-cleavage steps (promoting versus inhibiting).     

Dynamic Effects Dictate the Mechanism and Selectivity of Dehydration–Rearrangement Reactions of Protonated Alcohols [Me2(R)CCH(OH2)Me]+ (R=Me,Et, iPr) in the Gas Phase

The gas‐phase dehydration–rearrangement (DR) reactions of protonated alcohols [Me₂(R)CCH(OH₂)Me]⁺ [R=Me ( ME ), Et ( ET ), and iPr ( I‐PR )] were studied by using static approaches (intrinsic reaction coordinate (IRC), Rice–Ramsperger–Kassel–Marcus theory) and dynamics (quasiclassical trajectory) simulations at the B3LYP/6‐31G(d) level of theory. The concerted mechanism involves simultaneous water dissociation and alkyl migration, whereas in the stepwise reaction pathway the dehydration step leads to a secondary carbocation intermediate followed by alkyl migration. Internal rotation (IR) can change the relative position of the migrating alkyl group and the leaving group (water), so distinct products may be obtained: [Me(R)CCH(Me)Me ??? OH₂]⁺ and [Me(Me)CCH(R)Me ??? OH₂]⁺. The static approach predicts that these reactions are concerted, with the selectivity towards these different products determined by the proportion of the conformers of the initial protonated alcohols. These selectivities are explained by the DR processes being much faster than IR. These results are in direct contradiction with the dynamics simulations, which indicate a predominantly stepwise mechanism and selectivities that depend on the alkyl groups and dynamics effects. Indeed, despite the lifetimes of the secondary carbocations being short (<0.5 ps), IR can take place and thus provide a rich selectivity. These different selectivities, particularly for ET and I‐PR , are amenable to experimental observation and provide evidence for the minor role played by potential‐energy surface and the relevance of the dynamics effects (non‐IRC pathways, IR) in determining the reaction mechanisms and product distribution (selectivity).     

Quantum tautomerization in porphycene and its isotopomers: Path-integral molecular dynamics simulations

Path-integral molecular dynamics simulations have been performed for porphycene and its isotopic variants in order to understand the effect of isotopic substitution of inner protons on the double proton transfer mechanism. We have used an on-the-fly direct dynamics technique at the semiempirical PM6 level combined with specific reaction parameterization. Our quantum simulations show that double proton transfer of the unsubstituted porphycene at T = 300 K mainly occurs via a so-called concerted mechanism through the D_2h second-order saddle point. In addition, we found that both isotopic substitution and temperature significantly affect the double proton transfer mechanism. For example, the contribution of the stepwise mechanism increases with a temperature increase. We have also carried out hypothetical simulations with the porphycene configurations being completely planar. It has been found that out-of-plane vibrational motions significantly decrease the contribution of the concerted proton transfer mechanism.     

Solvation of radical cations in water--reactive or unreactive solvation?

The solvation and reaction of ethylene radical cation in aqueous solution has been studied with Car-Parrinello molecular dynamics simulations. All ab initio simulations were performed using a system of 56 water and one ethylene molecule. Using a favorable symmetrically solvated radical cation as the starting point of the simulation a fast addition of water (within 90 fs) to the radical cation is observed. The primary addition product is rapidly deprotonated (within 100 fs) to yield the ethanol-2-yl radical. A second simulation was initiated through vertical ionization of neutral hydrated ethylene, representing a significantly less favorable situation for the addition process. No addition of water can be observed in this second simulation over a time span of 1.7 ps. Taken together the two simulations are indicative of a rearrangement of the solvent shell which represents the major part of the overall reaction barrier. Under these circumstances, the reaction rate of an otherwise spontaneous reaction is limited by the intrinsic solvent relaxation time. This interpretation of the reactivity of hydrated radical cations reconciles previously conflicting experimental condensed phase and theoretical gas phase studies.     

Vibrational spectral diffusion and hydrogen bond dynamics in heavy water from first principles

We present a first-principles theoretical study of vibrational spectral diffusion and hydrogen bond dynamics in heavy water without using any empirical model potentials. The calculations are based on ab initio molecular dynamics simulations for trajectory generation and a time series analysis using the wavelet method for frequency calculations. It is found that, in deuterated water, although a one-to-one relation does not exist between the instantaneous frequency of an OD bond and the distance of its associated hydrogen bond, such a relation does hold on average. The dynamics of spectral diffusion is investigated by means of frequency-time correlation and spectral hole dynamics calculations. Both of these functions are found to have a short-time decay with a time scale of approximately 100 fs corresponding to dynamics of intact hydrogen bonds and a slower long-time decay with a time constant of approximately 2 ps corresponding to lifetimes of hydrogen bonds. The connection of the slower time scale to the dynamics of local structural relaxation is also discussed. The dynamics of hydrogen bond making is shown to have a rather fast time scale of approximately 100 fs; hence, it can also contribute to the short-time dynamics of spectral diffusion. A damped oscillation is also found at around 150-200 fs, which is shown to have come from underdamped intermolecular vibrations of a hydrogen-bonded water pair. Such assignments are confirmed by independent calculations of power spectra of intermolecular motion and hydrogen bond kinetics using the population correlation function formalism. The details of the time constants of frequency correlations and spectral shifts are found to depend on the frequencies of chosen OD bonds and are analyzed in terms of the dynamics of hydrogen bonds of varying strengths and also of free non-hydrogen-bonded OD groups.     

Theoretical thermodynamic study of the adenine–thymine tautomeric equilibrium: Electronic structure calculations and steered molecular dynamic simulations

Free energy of the tautomeric equilibrium A‐T ? A*‐T* between the canonical and noncanonical DNA base equilibrium in aqueous solution was theoretically determined by applying electronic structure methods (at the M06‐2X‐PCM/6‐311++G(d,p) level) and steered molecular dynamic simulations. Concerted and stepwise mechanisms were considered for the double proton transfer in an effort to explain the anomalous behavior of this system where an unfavorable process without a transition state can be observed depending on the level of calculation used. Of the different mechanisms used in the simulations, the stepwise mechanism, in which the first step implies the transference of a proton from thymine to adenine, and a second step with the transference of a different proton from adenine to thymine, was the only one that showed two transition states and a reaction intermediate. However, a concerted and stepwise mechanism has similar kinetic and thermodynamic behavior, with similar reaction and activation energies. Simple proton transfer was more favorable for the transference of the hydrogen from the adenine to the thymine. The inclusion of an aqueous medium in this study only slightly modified these energies, but the barrier energy was higher when the solvent was described as a discrete medium. Transition states and intermediate structures were analyzed at molecular dynamic level.     

Increasing the time step with mass scaling in Born‐Oppenheimer ab initio QM/MM molecular dynamics simulations

Born‐Oppenheimer ab initio QM/MM molecular dynamics simulation with umbrella sampling is a state‐of‐the‐art approach to calculate free energy profiles of chemical reactions in complex systems. To further improve its computational efficiency, a mass‐scaling method with the increased time step in MD simulations has been explored and tested. It is found that by increasing the hydrogen mass to 10 amu, a time step of 3 fs can be employed in ab initio QM/MM MD simulations. In all our three test cases, including two solution reactions and one enzyme reaction, the resulted reaction free energy profiles with 3 fs time step and mass scaling are found to be in excellent agreement with the corresponding simulation results using 1 fs time step and the normal mass. These results indicate that for Born‐Oppenheimer ab initio QM/MM molecular dynamics simulations with umbrella sampling, the mass‐scaling method can significantly reduce its computational cost while has little effect on the calculated free energy profiles. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2009     

Dynamics of photoinduced processes in adenine and thymine base pairs

The excited-state dynamics of adenine and thymine dimers and the adenine-thymine base pair were investigated by femtosecond pump-probe ionization spectroscopy with excitation wavelengths of 250-272 nm. The base pairs showed a characteristic ultrafast decay of the initially excited pi pi* state to an n pi* state (lifetime tau(pi pi*) approximately 100 fs) followed by a slower decay of the latter with tau(n pi*) approximately 0.9 ps for (adenine)2, tau(n pi*) = 6-9 ps for (thymine)2, and tau(n pi*) approximately 2.4 ps for the adenine-thymine base pair. In the adenine dimer, a competing decay of the pi pi* state via the pi sigma* state greatly suppressed the n pi* state signals. Similarities of the excited-state decay parameters in the isolated bases and the base pairs suggest an intramonomer relaxation mechanism in the base pairs.     

A molecular spring for vision

Light absorption by the visual pigment rhodopsin leads to vision via a complex signal transduction pathway that is initiated by the ultrafast and highly efficient photoreaction of its chromophore, the retinal protonated Schiff base (RPSB). Here, we investigate this reaction in real time by means of unrestrained molecular dynamics simulations of the protein in a membrane mimetic environment, treating the chromophore at the density functional theory level. We demonstrate that a highly strained all-trans RPSB is formed starting from the 11-cis configuration (dark state) within approximately 100 fs by a minor rearrangement of the nuclei under preservation of the saltbridge with Glu113 and virtually no deformation of the binding pocket. Hence, the initial step of vision can be understood as the compression of a molecular spring by a minor change of the nuclear coordinates. This spring can then release its strain by altering the protein environment.     

Macrophomate synthase: QM/MM simulations address the Diels-Alder versus Michael-Aldol reaction mechanism

Macrophomate synthase (MPS) of the phytopathogenic fungus Macrophoma commelinae catalyzes the transformation of 2-pyrone derivatives to the corresponding benzoate analogues. The transformation proceeds through three separate chemical reactions, including decarboxylation of oxalacetate to produce pyruvate enolate, two C-C bond formations between 2-pyrone and pyruvate enolate that form a bicyclic intermediate, and final decarboxylation with concomitant dehydration. Although some evidence suggests that the second step of the reaction catalyzed by MPS is a Diels-Alder reaction, definite proof that the C-C bond formations occur via a concerted mechanism is still required. An alternative route for formation of the C-C bonds is a stepwise Michael-aldol reaction. In this work, mixed quantum and molecular mechanics (QM/MM) combined with Monte Carlo simulations and free-energy perturbation (FEP) calculations were performed to investigate the relative stabilities of the transition states (TS) for both reaction mechanisms. The key results are that the Diels-Alder TS model is 17.7 and 12.1 kcal/mol less stable than the Michael and aldol TSs in the stepwise route, respectively. Therefore, this work indicates that the Michael-aldol mechanism is the route used by MPS to catalyze the second step of the overall transformation, and that the enzyme is not a natural Diels-Alderase, as claimed by Ose and co-workers (Nature 2003, 422, 185-189; Acta Crystallogr. 2004, D60, 1187-1197). A modified link-atom treatment for the bonds at the QM/MM interface is also presented.     

Quasi-Hamiltonian equations of motion for internal coordinate molecular dynamics of polymers

Conventional molecular dynamics simulations of macromolecules require long computational times because the most interesting motions are very slow compared to the fast oscillations of bond lengths and bond angles that limit the integration time step. Simulation of dynamics in the space of internal coordinates, that is, with bond lengths, bond angles, and torsions as independent variables, gives a theoretical possibility of eliminating all uninteresting fast degrees of freedom from the system. This article presents a new method for internal coordinate molecular dynamics simulations of macromolecules. Equations of motion are derived that are applicable to branched chain molecules with any number of internal degrees of freedom. Equations use the canonical variables and they are much simpler than existing analogs. In the numerical tests the internal coordinate dynamics are compared with the traditional Cartesian coordinate molecular dynamics in simulations of a 56 residue globular protein. For the first time it was possible to compare the two alternative methods on identical molecular models in conventional quality tests. It is shown that the traditional and internal coordinate dynamics require the same time step size for the same accuracy and that in the standard geometry approximation of amino acids, that is, with fixed bond lengths, bond angles, and rigid aromatic groups, the characteristic step size is 4 fs, which is 2 times higher than with fixed bond lengths only. The step size can be increased up to 11 fs when rotation of hydrogen atoms is suppressed. © 1997 by John Wiley & Sons, Inc. J Comput Chem 18 : 1354–1364, 1997     

Quantum dynamics simulations of interfacial electron transfer in sensitized TiO2 semiconductors

Ab initio DFT molecular dynamics simulations are combined with quantum dynamics calculations of electronic relaxation to investigate the interfacial electron transfer in catechol/TiO(2)-anatase nanostructures under vacuum conditions. It is found that the primary process in the interfacial electron-transfer dynamics involves an ultrafast (tau(1) approximately 6 fs) electron-injection event that localizes the charge in the Ti(4+) surface ions next to the catechol adsorbate. The primary event is followed by charge delocalization (i.e., carrier diffusion) through the TiO(2)-anatase crystal, an anisotropic diffusional process that can be up to an order of magnitude slower along the [-101] direction than carrier relaxation along the [010] and [101] directions in the anatase crystal. It is shown that both the mechanism of electron injection and the time scales for interfacial electron transfer are quite sensitive to the symmetry of the electronic state initially populated in the adsorbate molecule. The results are particularly relevant to the understanding of surface charge separation in efficient mechanisms of molecular-based photovoltaic devices.