Structure and ultrafast dynamics of liquid water: a quantum mechanics/molecular mechanics molecular dynamics simulations study

A quantum mechanics/molecular mechanics molecular dynamics simulation was performed for liquid water to investigate structural and dynamical properties of this peculiar liquid. The most important region containing a central reference molecule and all nearest surrounding molecules (first coordination shell) was treated by Hartree-Fock (HF), post-Hartree-Fock [second-order Moller-Plesset perturbation theory (MP2)], and hybrid density functional B3LYP [Becke's three parameter functional (B3) with the correlation functional of Lee, Yang, and Parr (LYP)] methods. In addition, another HF-level simulation (2HF) included the full second coordination shell. Site to site interactions between oxygen-oxygen, oxygen-hydrogen, and hydrogen-hydrogen atoms of all ab initio methods were compared to experimental data. The absence of a second peak and the appearance of a shoulder instead in the gO-O graph obtained from the 2HF simulation is notable, as this feature has been observed so far only for pressurized or heated water. Dynamical data show that the 2HF procedure compensates some of the deficiency of the HF one-shell simulation, reducing the difference between correlated (MP2) and HF results. B3LYP apparently leads to too rigid structures and thus to an artificial slow down of the dynamics.     

Optimal control of ultrafast laser driven many-electron dynamics in a polyatomic molecule: N-methyl-6-quinolone

We report time-dependent configuration interaction singles calculations for the ultrafast laser driven many-electron dynamics in a polyatomic molecule, N-methyl-6-quinolone. We employ optimal control theory to achieve a nearly state-selective excitation from the S(0) to the S(1) state, on a time scale of a few ( approximately 6) femtoseconds. The optimal control scheme is shown to correct for effects opposing a state-selective transition, such as multiphoton transitions and other, nonlinear phenomena, which are induced by the ultrashort and intense laser fields. In contrast, simple two-level pi pulses are not effective in state-selective excitations when very short pulses are used. Also, the dependence of multiphoton and nonlinear effects on the number of states included in the dynamical simulations is investigated.     

A theoretical study on quantum control of photodissociation dynamics by ultrashort chirped laser pulses

By a wavepacket propagation, we demonstrate the possibility of controlling the photodissociation branching ratio between two fragment channels by a single ultrashort linearly chirped laser pulse. It is found that a negatively chirped pulse of a moderate chirp rate completely prohibits the production of one of the photofragment channels. Two characteristics of chirped laser pulses contribute to this remarkable effect: the mechanism of adiabatic rapid passage (ARP) for the population transfer between the ground and excited states and the intrapulse pump‐dump process for determining the branching ratio. ©1999 John Wiley & Sons, Inc. Int J Quant Chem 72: 525–532, 1999     

Hybrid quantum mechanics/molecular mechanics-based molecular dynamics simulation of acid-catalyzed dehydration of polyols in liquid water

We use the conversion of protonated glycerol to acrolein for a case study of the mechanism of acid-catalyzed dehydration of polyols in aqueous environments. We employ hybrid Quamtum Mechanics/Molecular Mechanics Molecular Dynamics (QM/MM MD) simulations with biased sampling and perform free energy calculations for the elementary steps of the reaction. We investigate the effects of solvent dynamics and in particular the role of quantum mechanical water in the dehydration mechanism. We present results supporting a mechanism that proceeds via water-mediated proton transfers and thus through an enol intermediate. We find that the first dehydration may take place by two, low-energy pathways requiring, respectively, 20.9 and 18.8 kcal/mol of activation free energy. The second dehydration requires 19.9 kcal/mol of activation free energy while for the overall reaction we compute a free energy change of -8 kcal/mol.     

Computer simulation of quantum dynamics in a classical spin environment

In this paper, a formalism for studying the dynamics of quantum systems coupled to classical spin environments is reviewed. The theory is based on generalized antisymmetric brackets and naturally predicts open-path off-diagonal geometric phases in the evolution of the density matrix. It is shown that such geometric phases must also be considered in the quantum–classical Liouville equation for a classical bath with canonical phase space coordinates; this occurs whenever the adiabatics basis is complex (as in the case of a magnetic field coupled to the quantum subsystem). When the quantum subsystem is weakly coupled to the spin environment, non-adiabatic transitions can be neglected and one can construct an effective non-Markovian computer simulation scheme for open quantum system dynamics in classical spin environments. In order to tackle this case, integration algorithms based on the symmetric Trotter factorization of the classical-like spin propagator are derived. Such algorithms are applied to a model comprising a quantum two-level system coupled to a single classical spin in an external magnetic field. Starting from an excited state, the population difference and the coherences of this two-state model are simulated in time while the dynamics of the classical spin is monitored in detail. It is the author’s opinion that the numerical evidence provided in this paper is a first step toward developing the simulation of quantum dynamics in classical spin environments into an effective tool. In turn, the ability to simulate such a dynamics can have a positive impact on various fields, among which, for example, nanoscience.     

Non-adiabatic molecular dynamics simulation of ultrafast solar cell electron transfer

A non-adiabatic (NA) molecular dynamics (MD) simulation of the photoinduced electron transfer (ET) from a molecular electron donor to the TiO₂ acceptor is carried out. The system under study is typical of the dye sensitized semiconductor nanomaterials used in solar cells, photocatalysis and photoelectrolysis. The electronic structure of the dye-semiconductor system and the adiabatic dynamics are simulated by ab initio density functional theory MD, while the NA effects are incorporated by the quantum-classical mean-field approach. A novel procedure separating the NA and adiabatic ET pathways is developed. The simulation provides a detailed picture of the ET process. For the specific system under study, ET occurs on a 30 fs time scale, in agreement with the ultrafast experimental data. Both adiabatic and NA pathways for the ET are observed. The NA transfer entirely dominates at short times and can occur due to strong localized avoided crossing as well as extended regions of weaker NA coupling. Although the adiabatic ET contribution accumulates more slowly, it approaches that of the NA ET pathway asymptotically. The electron acceptor states are formed by the d-orbital of Ti atoms of the semiconductor and are localized within the first 3–4 layers of the surface. About 20% of the acceptor state density is localized on a single Ti atom of the first surface layer. The simulation predicts a complex non-single-exponential time dependence of the ET process.     

Finite element interpolation for combined classical/quantum mechanical molecular dynamics simulations

A method is presented to interpolate the potential energy function for a part of a system consisting of a few degrees of freedom, such as a molecule in solution. The method is based on a modified finite element (FE) interpolation scheme. The aim is to save computer time when expensive methods such as quantum-chemical calculations are used to determine the potential energy function. The expensive calculations are only carried out if the molecule explores new unknown regions of the conformation space. If the molecule resides in regions previously explored, a cheap interpolation is performed instead of an expensive calculation, using known neighboring points. We report the interpolation techniques for the energies and the forces of the molecule, the handling of the FE mesh, and an application to a simple test example in molecular dynamics (MD) simulations. Good performance of the method was obtained (especially for MD simulations with a preceding Monte Carlo mesh generation) without losing accuracy. © 1997 John Wiley & Sons, Inc. J Comput Chem 18 : 1484–1495, 1997     

Dynamic molecular interferometer: probe of inversion symmetry in I2(-) photodissociation

Time-resolved photoelectron imaging of negative ions is employed to examine 780-nm dissociation dynamics of I2(-), emphasizing the effects of interference in time-resolved photoelectron angular distributions obtained with 390-nm probe. No energetic changes are observed after about 700 fs, but the evolution of the photoelectron anisotropy persists for up to 2.5 ps, indicating that the electronic wave function of the dissociating anion continues to evolve long after the asymptotic energetic limit of the reaction has been effectively reached. The time scale of the anisotropy variation corresponds to a fragment separation of the same order of magnitude as the de Broglie wavelength of the emitted electrons (lambda=35 A). These findings are interpreted by considering the effect of I2(-) inversion symmetry and viewing the dissociating anion as a dynamic molecular-scale "interferometer," with the electron waves emitted from two separating centers. The predictions of the model are in agreement with the present experiment and shed new light on previously published results [A. V. Davis, R. Wester, A. E. Bragg, and D. M. Neumark, J. Chem. Phys. 118, 999 (2003)].     

Optimal control of classical molecular dynamics: A perturbation formulation and the existence of multiple solutions

This paper considers the prospect for there being multiple solutions to the control of classically modelled molecular dynamical systems. The research presented here follows up on a parallel study based on quantum mechanics. For polyatomic molecules it is generally expected that a classical mechanical model will be adequate and necessary as a means for designing optical fields for molecular control. The prospect for multiple control field solutions existing in this domain is important to establish in terms of ultimate laboratory realization of molecular control. A general formulation of the multiplicity problem is considered and the existence of a denumerably infinite number of solutions for the control field amplitude is shown to be the case under certain mild limitations on the physical variables.     

Hybrid quantum/classical studies of photodissociation and recombination of I2 (A) in rare gas matrices: A linear chain model

A hybrid quantum/classical model is developed for the photodissociation and recombination dynamics of an I₂ molecule in low-temperature rare-gas (Rg) matrices. The simplified model consists of an I₂ molecule embedded in a linear chain of Rg atoms. The aggregate is partitioned into a quantum system and a classical bath, which are self-consistently coupled. Two partitioning schemes are used. The first treats the I-I coordinate quantum mechanically and the Rg coordinates classically. The second and more reliable scheme includes in the quantum system both the I-I mode and the symmetric motion of the two nearest Rg atoms. Both models show substantial energy transfer from the dissociating I₂ to the solvent, followed by coherent vibrational motion of the recombined I₂. It is found that the one-dimensional quantum/classical scheme is consistent with its higher dimensional counterpart, although the latter shows much faster dephasing. © 1996 John Wiley & Sons, Inc.     

The structure of the hydrated electron. Part 2. A mixed quantum/classical molecular dynamics embedded cluster density functional theory: single-excitation configuration interaction study

Adiabatic mixed quantum/classical (MQC) molecular dynamics (MD) simulations were used to generate snapshots of the hydrated electron in liquid water at 300 K. Water cluster anions that include two complete solvation shells centered on the hydrated electron were extracted from the MQC MD simulations and embedded in a roughly 18 Ax18 Ax18 A matrix of fractional point charges designed to represent the rest of the solvent. Density functional theory (DFT) with the Becke-Lee-Yang-Parr functional and single-excitation configuration interaction (CIS) methods were then applied to these embedded clusters. The salient feature of these hybrid DFT(CIS)/MQC MD calculations is significant transfer (approximately 18%) of the excess electron's charge density into the 2p orbitals of oxygen atoms in OH groups forming the solvation cavity. We used the results of these calculations to examine the structure of the singly occupied and the lower unoccupied molecular orbitals, the density of states, the absorption spectra in the visible and ultraviolet, the hyperfine coupling (hfcc) tensors, and the infrared (IR) and Raman spectra of these embedded water cluster anions. The calculated hfcc tensors were used to compute electron paramagnetic resonance (EPR) and electron spin echo envelope modulation (ESEEM) spectra for the hydrated electron that compared favorably to the experimental spectra of trapped electrons in alkaline ice. The calculated vibrational spectra of the hydrated electron are consistent with the red-shifted bending and stretching frequencies observed in resonance Raman experiments. In addition to reproducing the visible/near IR absorption spectrum, the hybrid DFT model also accounts for the hydrated electron's 190-nm absorption band in the ultraviolet. Thus, our study suggests that to explain several important experimentally observed properties of the hydrated electron, many-electron effects must be accounted for: one-electron models that do not allow for mixing of the excess electron density with the frontier orbitals of the first-shell solvent molecules cannot explain the observed magnetic, vibrational, and electronic properties of this species. Despite the need for multielectron effects to explain these important properties, the ensemble-averaged radial wavefunctions and energetics of the highest occupied and three lowest unoccupied orbitals of the hydrated electrons in our hybrid model are close to the s- and p-like states obtained in one-electron models. Thus, one-electron models can provide a remarkably good approximation to the multielectron picture of the hydrated electron for many applications; indeed, the two approaches appear to be complementary.     

Exploring the role of decoherence in condensed-phase nonadiabatic dynamics: a comparison of different mixed quantum/classical simulation algorithms for the excited hydrated electron

Mixed quantum/classical (MQC) molecular dynamics simulation has become the method of choice for simulating the dynamics of quantum mechanical objects that interact with condensed-phase systems. There are many MQC algorithms available, however, and in cases where nonadiabatic coupling is important, different algorithms may lead to different results. Thus, it has been difficult to reach definitive conclusions about relaxation dynamics using nonadiabatic MQC methods because one is never certain whether any given algorithm includes enough of the necessary physics. In this paper, we explore the physics underlying different nonadiabatic MQC algorithms by comparing and contrasting the excited-state relaxation dynamics of the prototypical condensed-phase MQC system, the hydrated electron, calculated using different algorithms, including: fewest-switches surface hopping, stationary-phase surface hopping, and mean-field dynamics with surface hopping. We also describe in detail how a new nonadiabatic algorithm, mean-field dynamics with stochastic decoherence (MF-SD), is to be implemented for condensed-phase problems, and we apply MF-SD to the excited-state relaxation of the hydrated electron. Our discussion emphasizes the different ways quantum decoherence is treated in each algorithm and the resulting implications for hydrated-electron relaxation dynamics. We find that for three MQC methods that use Tully's fewest-switches criterion to determine surface hopping probabilities, the excited-state lifetime of the electron is the same. Moreover, the nonequilibrium solvent response function of the excited hydrated electron is the same with all of the nonadiabatic MQC algorithms discussed here, so that all of the algorithms would produce similar agreement with experiment. Despite the identical solvent response predicted by each MQC algorithm, we find that MF-SD allows much more mixing of multiple basis states into the quantum wave function than do other methods. This leads to an excited-state lifetime that is longer with MF-SD than with any method that incorporates nonadiabatic effects with the fewest-switches surface hopping criterion.     

Femtosecond pump-probe spectroscopy of I2 in a dense rare gas environment: a mixed quantum/classical study of vibrational decoherence

The process of decoherence of vibrational states of I2 in a dense helium environment is studied theoretically using the mixed quantum/classical method based on the Bohmian formulation of quantum mechanics [E. Gindensperger, C. Meier, and J. A. Beswick, J. Chem. Phys. 113, 9369 (2000)]. Specifically, the revival of vibrational wave packets is a quantum phenomena which depends sensitively on the coherence between the vibrational states excited by an ultrafast laser pulse. Its detection by a pump-probe setup as a function of rare gas pressure forms a very accurate way of detecting vibrational dephasing. Vibrational revivals of I2 in high pressure rare gas environments have been observed experimentally, and the very good agreement with the simulated spectra confirms that the method can accurately describe decoherence processes of quantum systems in interaction with an environment.     

Sulfur dioxide in water: structure and dynamics studied by an ab initio quantum mechanical charge field molecular dynamics simulation

An ab initio Quantum Mechanical Charge Field Molecular Dynamics Simulation (QMCF MD) was performed to investigate structure and dynamics behavior of hydrated sulfur dioxide (SO(2)) at the Hartree-Fock level of theory employing Dunning DZP basis sets for solute and solvent molecules. The intramolecular structural characteristics of SO(2), such as S═O bond lengths and O═S═O bond angle, are in good agreement with the data available from a number of different experiments. The structural features of the hydrated SO(2) were primarily evaluated in the form of S-O(wat) and O(SO(2))-H(wat) radial distribution functions (RDFs) which gave mean distances of 2.9 and 2.2 ?, respectively. The dynamical behavior characterizes the solute molecule to have structure making properties in aqueous solution or water aerosols, where the hydrated SO(2) can easily get oxidized to form a number of sulfur(VI) species, which are believed to play an important role in the atmospheric processes.     

Moving solvated electrons with light: nonadiabatic mixed quantum/classical molecular dynamics simulations of the relocalization of photoexcited solvated electrons in tetrahydrofuran (THF)

Motivated by recent ultrafast spectroscopic experiments [Martini et al., Science 293, 462 (2001)], which suggest that photoexcited solvated electrons in tetrahydrofuran (THF) can relocalize (that is, return to equilibrium in solvent cavities far from where they started), we performed a series of nonequilibrium, nonadiabatic, mixed quantum/classical molecular dynamics simulations that mimic one-photon excitation of the THF-solvated electron. We find that as photoexcited THF-solvated electrons relax to their ground states either by continuous mixing from the excited state or via nonadiabatic transitions, approximately 30% of them relocalize into cavities that can be over 1 nm away from where they originated, in close agreement with the experiments. A detailed investigation shows that the ability of excited THF-solvated electrons to undergo photoinduced relocalization stems from the existence of preexisting cavity traps that are an intrinsic part of the structure of liquid THF. This explains why solvated electrons can undergo photoinduced relocalization in solvents like THF but not in solvents like water, which lack the preexisting traps necessary to stabilize the excited electron in other places in the fluid. We also find that even when they do not ultimately relocalize, photoexcited solvated electrons in THF temporarily visit other sites in the fluid, explaining why the photoexcitation of THF-solvated electrons is so efficient at promoting recombination with nearby scavengers. Overall, our study shows that the defining characteristic of a liquid that permits the photoassisted relocalization of solvated electrons is the existence of nascent cavities that are attractive to an excess electron; we propose that other such liquids can be found from classical computer simulations or neutron diffraction experiments.     

Solvation of xyloglucan in water/alcohol systems by molecular dynamics simulation

Xyloglucan in water solution turns into a gel with addition of alcohol such as methanol and ethanol. In regard to this phenomenon, we investigated the adhesive property of alcohol to xyloglucan and proposed the mechanism of the gelation by molecular dynamics (MD) simulation of a xyloglucan in water, water/methanol, and water/ethanol solution for 10 ns. The alcohol molecules showed its adhesive property to the xyloglucan and made the swelling-shrinking motion of the xyloglucan slow. Alcohol molecules solvated to the xyloglucan mainly in hydrophobic way so as to fill the void of water hydration shell, resulting in reformation of the hydrogen-bond network of water molecules around the solute. We also found that alcohol molecules have strong tendency to hydrogen-bond on xylose O3 in xyloglucan. According to these results, we proposed the gelation mechanism of xyloglucan in water/alcohol solution.     

Ab initio nonadiabatic quantum dynamics of cyclohexadiene/hexatriene ultrafast photoisomerization

Reaction mechanisms of the ultrafast photoisomerization between cyclohexadiene and hexatriene have been elucidated by the quantum dynamics on the ab initio potential energy surfaces calculated by multireference configuration interaction method. In addition to the quantum wave-packet dynamics along the two-dimensional reaction coordinates, the semiclassical analyses have also been carried out to correctly estimate the nonadiabatic transition probabilities around conical intersections in the full-dimensional space. The reaction time durations of radiationless decays in the wave-packet dynamics are found to be generally consistent with the femtosecond time-resolution experimental observations. The nonadiabatic transition probabilities among the ground (S0), first (S1), and second (S2) excited states have been estimated by using the semiclassical Zhu-Nakamura formula considering the full-dimensional wave-packet density distributions in the vicinity of conical intersections under the harmonic normal mode approximation. The cyclohexadiene (CHD) ring-opening process proceeds descending on the S1(1 1B) potential after the photoexcitation. The major part of the wave-packet decays from S1(1 1B) to S1(2 1A) by the first seam line crossing along the C2-symmetry-breaking directions. The experimentally observed ultrafast S1-S0 decay can be explained by the dynamics through the S1-S0 conical intersection along the direction toward the five-membered ring. The CHD: hexatriene (HT) branching ratio is estimated to be approximately 5:5, which is in accordance with the experiment in solution. This branching ratio is found to be mainly governed by the location of the five-membered ring S1-S0 conical intersection along the ground state potential ridge between CHD and HT.     

IR spectra of N-methylacetamide in water predicted by combined quantum mechanical/molecular mechanical molecular dynamics simulations

We applied the combined quantum mechanical (QM)/molecular mechanical (MM) molecular dynamics (MD) simulation method in assessing IR spectra of N-methylacetamide and its deuterated form in aqueous solutions. The model peptide is treated at the Austin Model 1 (AM1) level and the induced dipole effects by the solvent are incorporated in fluctuating solute dipole moments, which are calculated using partial charges from Mulliken population analyses without resorting to any available high-level ab initio dipole moment data. Fourier transform of the solute dipole autocorrelation function produces in silico IR spectra, in which the relative peak intensities and bandwidths of major amide bands are quantitatively compatible with experimental results only when both geometric and electronic polarizations of the peptide by the solvent are dealt with at the same quantum-mechanical level. We cast light on the importance of addressing dynamic charge fluctuations of the solute in calculating IR spectra by comparing classical and QM/MM MD simulation results. We propose the adjustable scaling factors for each amide mode to be directly compared with experimental data.