First principles and classical molecular dynamics simulations of solvated benzene

We have performed extensive ab initio and classical molecular dynamics (MD) simulations of benzene in water in order to examine the unique solvation structures that are formed. Qualitative differences between classical and ab initio MD simulations are found and the importance of various technical simulation parameters is examined. Our comparison indicates that nonpolarizable classical models are not capable of describing the solute-water interface correctly if local interactions become energetically comparable to water hydrogen bonds. In addition, a comparison is made between a rigid water model and fully flexible water within ab initio MD simulations which shows that both models agree qualitatively for this challenging system.     

Matching-pursuit/split-operator Fourier-transform simulations of nonadiabatic quantum dynamics

A rigorous and practical approach for simulations of nonadiabatic quantum dynamics is introduced. The algorithm involves a natural extension of the matching-pursuitsplit-operator Fourier-transform (MPSOFT) method [Y. Wu and V. S. Batista, J. Chem. Phys. 121, 1676 (2004)] recently developed for simulations of adiabatic quantum dynamics in multidimensional systems. The MPSOFT propagation scheme, extended to nonadiabatic dynamics, recursively applies the time-evolution operator as defined by the standard perturbation expansion to first-, or second-order, accuracy. The expansion is implemented in dynamically adaptive coherent-state representations, generated by an approach that combines the matching-pursuit algorithm with a gradient-based optimization method. The accuracy and efficiency of the resulting propagation method are demonstrated as applied to the canonical model systems introduced by Tully for testing simulations of dual curve-crossing nonadiabatic dynamics.     

Trajectory-based solution of the nonadiabatic quantum dynamics equations: an on-the-fly approach for molecular dynamics simulations

The non-relativistic quantum dynamics of nuclei and electrons is solved within the framework of quantum hydrodynamics using the adiabatic representation of the electronic states. An on-the-fly trajectory-based nonadiabatic molecular dynamics algorithm is derived, which is also able to capture nuclear quantum effects that are missing in the traditional trajectory surface hopping approach based on the independent trajectory approximation. The use of correlated trajectories produces quantum dynamics, which is in principle exact and computationally very efficient. The method is first tested on a series of model potentials and then applied to study the molecular collision of H with H(2) using on-the-fly TDDFT potential energy surfaces and nonadiabatic coupling vectors.     

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     

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.     

The structures of ozone and HOx radicals in aqueous solution from combined quantum/classical molecular dynamics simulations

Ozone in aqueous solution decomposes through a complex mechanism that involves initial reaction with a hydroxide ion followed by formation of a variety of oxidizing species such as HO, HO(2), and HO(3) radicals. Though a number of hydrogen-bonded complexes have been described in the gas phase, both theoretically and experimentally, the structures of ozone and HO(x) in liquid water remain uncertain. In this work, combined quantum/classical computer simulations of aqueous solutions of these species have been reported. The results show that ozone undergoes noticeable electron polarization but it does not participate in hydrogen bonds with liquid water. The main contribution of the solvation energy comes from dispersion forces. In contrast, HO(x) radicals form strong hydrogen bonds. They are better proton donors but weaker proton acceptors than water. Their electronic and geometrical structures are significantly modified by the solvent, especially in the case of HO(3). In all cases, fluctuations in amplitudes of electronic properties are considerable, suggesting that solvent effects might play a crucial role on oxidation mechanisms initiated by ozone in liquid water. These mechanisms are important in a broad range of domains, such as atmospheric processes, plant response to ambient ozone, and medical and industrial applications.     

Influence of 3A molecular sieve on tetrahydrofuran (THF) hydrate formation

Visual observation of the THF hydrate formation process in the presence of a 3A molecular sieve has been made at normal atmosphere and below a temperature of zero by microscopy. The results indicate that a 3A molecular sieve can induce the nucleation of the THF hydrate and promote the THF hydrate growth. With the existence of a 3A molecular sieve, the growth rate of THF hydrate is between 0.01 and 0.05 μm/s. In comparison with the system without any 3A molecular sieve, the growth rate increases about 4 nm/s. After the THF hydrate grows into megacryst, the crystals will recombine and partially change under the same condition.     

Energy conservation in molecular dynamics simulations of classical systems

Classical Newtonian dynamics is analytic and the energy of an isolated system is conserved. The energy of such a system, obtained by the discrete "Verlet" algorithm commonly used in molecular dynamics simulations, fluctuates but is conserved in the mean. This is explained by the existence of a "shadow Hamiltonian" H [S. Toxvaerd, Phys. Rev. E 50, 2271 (1994)], i.e., a Hamiltonian close to the original H with the property that the discrete positions of the Verlet algorithm for H lie on the analytic trajectories of H. The shadow Hamiltonian can be obtained from H by an asymptotic expansion in the time step length. Here we use the first non-trivial term in this expansion to obtain an improved estimate of the discrete values of the energy. The investigation is performed for a representative system with Lennard-Jones pair interactions. The simulations show that inclusion of this term reduces the standard deviation of the energy fluctuations by a factor of 100 for typical values of the time step length. Simulations further show that the energy is conserved for at least one hundred million time steps provided the potential and its first four derivatives are continuous at the cutoff. Finally, we show analytically as well as numerically that energy conservation is not sensitive to round-off errors.     

Optimal laser control of ultrafast photodissociation of I2- in water: mixed quantum/classical molecular dynamics simulation

A linearized optimal control method in combination with mixed quantum/classical molecular dynamics simulation is used for numerically investigating the possibility of controlling photodissociation wave packets of I(2)(-) in water. Optimal pulses are designed using an ensemble of photodissociation samples, aiming at the creation of localized dissociation wave packets. Numerical results clearly show the effectiveness of the control although the control achievement is reduced with an increase in the internuclear distance associated with a target region. We introduce effective optimal pulses that are designed using a statistically averaged effective dissociation potential, and show that they semiquantitatively reproduce the control achievements calculated by using optimal pulses. The control mechanisms are interpreted from the time- and frequency-resolved spectra of the effective optimal pulses.     

Elastic bag model for molecular dynamics simulations of solvated systems: application to liquid water and solvated peptides

The fluctuating elastic boundary (FEB) model for molecular dynamics has recently been developed and validated through simulations of liquid argon. In the FEB model, a flexible boundary which consists of particles connected by springs is used to confine the solvated system, thereby eliminating the need for periodic boundary conditions. In this study, we extend this model to the simulation of bulk water and solvated alanine dipeptide. Both the confining potential and boundary particle interaction functions are modified to preserve the structural integrity of the boundary and prevent the leakage of the solute-solvent system through the boundary. A broad spectrum of structural and dynamic properties of liquid water are computed and compared with those obtained from conventional periodic boundary condition simulations. The applicability of the model to biomolecular simulations is investigated through the analysis of conformational population distribution of solvated alanine dipeptide. In most cases we find remarkable agreement between the two simulation approaches.     

Ejection of solvated ions from electrosprayed methanol/water nanodroplets studied by molecular dynamics simulations

The ejection of solvated small ions from nanometer-sized droplets plays a central role during electrospray ionization (ESI). Molecular dynamics (MD) simulations can provide insights into the nanodroplet behavior. Earlier MD studies have largely focused on aqueous systems, whereas most practical ESI applications involve the use of organic cosolvents. We conduct simulations on mixed water/methanol droplets that carry excess NH(4)(+) ions. Methanol is found to compromise the H-bonding network, resulting in greatly increased rates of ion ejection and solvent evaporation. Considerable differences in the water and methanol escape rates cause time-dependent changes in droplet composition. Segregation occurs at low methanol concentration, such that layered droplets with a methanol-enriched periphery are formed. This phenomenon will enhance the partitioning of analyte molecules, with possible implications for their ESI efficiencies. Solvated ions are ejected from the tip of surface protrusions. Solvent bridging prior to ion secession is more extensive for methanol/water droplets than for purely aqueous systems. The ejection of solvated NH(4)(+) is visualized as diffusion-mediated escape from a metastable basin. The process involves thermally activated crossing of a ~30 kJ mol(-1) free energy barrier, in close agreement with the predictions of the classical ion evaporation model.     

Solvation of excess electrons in LiF ionic pair matrix: evidence for a solvated dielectron from ab initio molecular dynamics simulations and calculations

Ab initio molecular dynamics simulations and first-principles calculations reveal the existence of a solvated dielectron species, (2e)s, in an LiF ionic matrix. The nature of the solvation mechanism and the stability of the species was explored. In addition to electrostatic interactions, a hole-orbital coupling among solvent molecules may significantly enhance the stability of the solvated electrons and govern the extent of electron solvation. This hole-orbital coupling is different from either an electrostatic coupling or conventional chemical bonding, and it may be described as a transition between them.     

Electrogenerated chemiluminescence with solvated electrons in hexamethylphosphoroamide (HMPA)

The electrogenerated chemiluminescence of N-toluenesulfonyl carbazole, 9-chlorofluorene derivatives and N,N,N′,N′-tetramethyl-p-phenylenediamine (TMPD) with solvated electrons in HMPA are discussed. The fairly bright emissions occurred simultaneously when solvated electrons were generated electrochemically. In the case of N-toluenesulfonyl carbazole and 9-chlorofluorene derivatives, the respective carbanions formed by “dissociative electron transfer reactions” are the fluorescents. The singlet state of TMPD seems to be directly formed by the electron transfer reactions between radical cations of TMPD and solvated electrons.     

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.     

Elastic bag model for molecular dynamics simulations of solvated systems: application to liquid argon

A new approach is developed to study the dynamics of the localized process in solutions and other condensed phase systems. The approach employs a fluctuating elastic boundary (FEB) model which encloses the simulated system in an elastic bag that mimics the effects of the bulk solvent. This alleviates the need for periodic boundary conditions and allows for a reduction in the number of solvent molecules that need to be included in the simulation. The boundary bag is modeled as a mesh of quasi-particles connected by elastic bonds. The FEB model allows for volume and density fluctuations characteristic of the bulk system, and the shape of the boundary fluctuates during the course of the simulation to adapt to the configuration fluctuations of the explicit solute-solvent system inside. The method is applied to the simulation of a Lennard-Jones model of liquid argon. Various structural and dynamical quantities are computed and compared with those obtained from conventional periodic boundary simulations. The agreement between the two is excellent in most cases, thus validating the viability of the FEB method.     

Anionic polymerization of styrene in tetrahydrofuran (THF) with cumylcaesium as initiator

The kinetics of the anionic polymerization of styrene initiated by cumylcaesium have been reinvestigated in tetrahydrofuran. Quite unlike earlier investigations, the ion pair rate constant $\text{k}{}_{\mathrm{p}}\text{(±)}$ was found to be 130 l/mole/sec at 25° and the dissociation constant of polystyrylcaesium to be about 5 × 10^?10 mole/l. In the presence of caesiumtriphenylcyanoborate used as a common ion salt addition, $\text{k}{}_{\mathrm{p}}\text{(±)}$ was measured over a wide range of temperature. The results are compared to the temperature behaviour of the ion pair rate constants obtained with polystyrylsodium in polar solvents; they indicate that the existence of solvent separated ion pairs cannot be excluded for polystyrylcaesium.