Simulation of silicon clusters from “quantum” Langevin molecular dynamics

We present a simulation method to determine from first principles the structure of low symmetry atomic systems. Our method is based on Langevin molecular dynamics and quantum mechanical interactions derived fromab initio pseudopotential calculations. The molecular dynamics time step with this approach can be one to two orders of magnitude larger than in the Car-Parrinello method, compensating for the time required for self-consistency at each step. Moreover, because the simulation is constrained to reside on the Born-Oppenheimer surface, this method can be used for insulating as well as metallic and charged systems. Application will be made to small silicon clusters.     

Molecular dynamics simulation study on stabilities and reactivities of NADH cytochrome B5 reductase

Binding free energies between coenzyme (FAD and NADH) and the apoenzyme of NADH-cytochrome b5 reductase (b5R) were estimated by applying the continuum Poisson-Boltzmann (PB) model to structures sampled from molecular dynamics simulations in explicit water molecules. Important residues for the enzymatic catalysis were clarified using a computational alanine scanning method. The binding free energies calculated by applying an alanine scanning method can successfully reproduce the trends of the measured steady-state enzymatic activities kcatNADH/KmNADH. Significant decreases in the binding free energy are expected when one of the four residues Arg91, Lys110, Ser127, and Thr181 is mutated into Ala. According to the results of the molecular dynamics simulation, Thr181 is considered to be one of the key residues that helps NADH to approach the isoalloxazine in FAD. Finally, we have constructed very simplified model systems and carried out density functional theory calculations using B3LYP/LANL2DZ//ROHF(or RHF)/LANL2DZ level of theory in order to elucidate a realistic and feasible mechanism of the hydride-ion transfer from NADH to FAD affected by HEME(Fe3+) as an electron acceptor. Our calculated results suggest that the electron and/or hydride-ion transfer reaction from NADH to FAD can be accelerated in the presence of HEME(Fe3+).     

Clustering of Lennard-Jones particles in water: temperature and pressure effects

While the hydrophobic effect is, for many systems, one of the most relevant interactions, it may be said that in the case of biological systems this effect becomes of determinant importance. Although the matter has been analyzed extensively, certain aspects are yet to be elucidated. Hence, the study on the behavior of the hydrophobic effect with temperature, and particularly with pressure deserves further investigation; model systems may help us in the task. We have analyzed the behavior of Lennard-Jones particles in water by means of molecular dynamics simulation under different conditions of size, concentration, temperature, and pressure. Following the formation of particle aggregates we can observe an increase of the hydrophobic effect with temperature and a strong weakening of the effect at high pressures. The results agree with the experimental evidence and show the ability of molecular dynamics simulation to account for the behavior of nonpolar substances under different conditions, provided that the intermolecular interactions used are adequate.     

Discontinuous molecular dynamics for rigid bodies: applications

Event-driven molecular dynamics simulations are carried out on two rigid-body systems which differ in the symmetry of their molecular mass distributions. First, simulations of methane in which the molecules interact via discontinuous potentials are compared with simulations in which the molecules interact through standard continuous Lennard-Jones potentials. It is shown that under similar conditions of temperature and pressure, the rigid discontinuous molecular dynamics method reproduces the essential dynamical and structural features found in continuous-potential simulations at both gas and liquid densities. Moreover, the discontinuous molecular dynamics approach is demonstrated to be between 3 and 100 times more efficient than the standard molecular dynamics method depending on the specific conditions of the simulation. The rigid discontinuous molecular dynamics method is also applied to a discontinuous-potential model of a liquid composed of rigid benzene molecules, and equilibrium and dynamical properties are shown to be in qualitative agreement with more detailed continuous-potential models of benzene. The few qualitative differences in the angular dynamics of the two models are related to the relatively crude treatment of variations in the discontinuous repulsive interactions as one benzene molecule rotates by another.     

Implementation of a molecular dynamics approach with rigid fragments to simulation of chemical reactions in biomolecular systems

We describe a new implementation of the molecular dynamics method aimed at simulation of the properties of biomolecular systems in which chemical reactions are possible. The quantum mechanical/molecular mechanical method based on the effective fragment potential theory is used for calculating the energies and forces along trajectories. Due to specific features of the effective fragment theory, the behavior of the molecular mechanical subsystem is described by rigid body dynamics. The method has been applied to simulation of proton transfer along the chain of water molecules inside the gramicidin channel.     

Toward the correction of effective electrostatic forces in explicit-solvent molecular dynamics simulations: restraints on solvent-generated electrostatic potential and solvent polarization

Despite considerable advances in computing power, atomistic simulations under nonperiodic boundary conditions, with Coulombic electrostatic interactions and in systems large enough to reduce finite-size associated errors in thermodynamic quantities to within the thermal energy, are still not affordable. As a result, periodic boundary conditions, systems of microscopic size and effective electrostatic interaction functions are frequently resorted to. Ensuing artifacts in thermodynamic quantities are nowadays routinely corrected a posteriori, but the underlying configurational sampling still descends from spurious forces. The present study addresses this problem through the introduction of on-the-fly corrections to the physical forces during an atomistic molecular dynamics simulation. Two different approaches are suggested, where the force corrections are derived from special potential energy terms. In the first approach, the solvent-generated electrostatic potential sampled at a given atom site is restrained to a target value involving corrections for electrostatic artifacts. In the second approach, the long-range regime of the solvent polarization around a given atom site is restrained to the Born polarization, i.e., the solvent polarization corresponding to the ideal situation of a macroscopic system under nonperiodic boundary conditions and governed by Coulombic electrostatic interactions. The restraints are applied to the explicit-water simulation of a hydrated sodium ion, and the effect of the restraints on the structural and energetic properties of the solvent is illustrated. Furthermore, by means of the calculation of the charging free energy of a hydrated sodium ion, it is shown how the electrostatic potential restraint translates into the on-the-fly consideration of the corresponding free-energy correction terms. It is discussed how the restraints can be generalized to situations involving several solute particles. Although the present study considers a very simple system only, it is an important step toward the on-the-fly elimination of finite-size and approximate-electrostatic artifacts during atomistic molecular dynamics simulations.     

Ab initio based polarizable force field parametrization

Experimental and simulation studies of anion-water systems have pointed out the importance of molecular polarization for many phenomena ranging from hydrogen-bond dynamics to water interfaces structure. The study of such systems at molecular level is usually made with classical molecular dynamics simulations. Structural and dynamical features are deeply influenced by molecular and ionic polarizability, which parametrization in classical force field has been an object of long-standing efforts. Although when classical models are compared to ab initio calculations at condensed phase, it is found that the water dipole moments are underestimated by approximately 30%, while the anion shows an overpolarization at short distances. A model for chloride-water polarizable interaction is parametrized here, making use of Car-Parrinello simulations at condensed phase. The results hint to an innovative approach in polarizable force fields development, based on ab initio simulations, which do not suffer for the mentioned drawbacks. The method is general and can be applied to the modeling of different systems ranging from biomolecular to solid state simulations.     

Isotropic periodic sum: a method for the calculation of long-range interactions

This work presents an accurate and efficient approach to the calculation of long-range interactions for molecular modeling and simulation. This method defines a local region for each particle and describes the remaining region as images of the local region statistically distributed in an isotropic and periodic way, which we call isotropic periodic images. Different from lattice sum methods that sum over discrete lattice images generated by periodic boundary conditions, this method sums over the isotropic periodic images to calculate long-range interactions, and is referred to as the isotropic periodic sum (IPS) method. The IPS method is not a lattice sum method and eliminates the need for a reciprocal space sum. Several analytic solutions of IPS for commonly used potentials are presented. It is demonstrated that the IPS method produces results very similar to that of Ewald summation, but with three major advantages, (1) it eliminates unwanted symmetry artifacts raised from periodic boundary conditions, (2) it can be applied to potentials of any functional form and to fully and partially homogenous systems as well as finite systems, and (3) it is more computationally efficient and can be easily parallelized for multiprocessor computers. Therefore, this method provides a general approach to an efficient calculation of long-range interactions for various kinds of molecular systems.     

Fragment-Based Ab Initio Molecular Dynamics Simulation for Combustion

We develop a fragment-based ab initio molecular dynamics (FB-AIMD) method for efficient dynamics simulation of the combustion process. In this method, the intermolecular interactions are treated by a fragment-based many-body expansion in which three- or higher body interactions are neglected, while two-body interactions are computed if the distance between the two fragments is smaller than a cutoff value. The accuracy of the method was verified by comparing FB-AIMD calculated energies and atomic forces of several different systems with those obtained by standard full system quantum calculations. The computational cost of the FB-AIMD method scales linearly with the size of the system, and the calculation is easily parallelizable. The method is applied to methane combustion as a benchmark. Detailed reaction network of methane reaction is analyzed, and important reaction species are tracked in real time. The current result of methane simulation is in excellent agreement with known experimental findings and with prior theoretical studies.     

A Dirac–Fock self-consistent field method for closed-shell molecules including Breit interaction

We present a method for including the Breit interaction in relativistic self-consistent field calculations for closed-shell molecular systems using atomic basis spinors of kinetically balanced Gaussian-type functions. The method extends the formalism described in a previous paper [A. Mohanty and E. Clementi, Int. J. Quantum Chem. 39 , 487–517 (1991)] that dealt with the two-electron effect due to Coulomb interaction only. It is shown that both frequency-dependent and frequency-independent Breit interactions can be treated on equal footing, and the corresponding matrix elements are evaluated following the well-known Fourier transform technique applied to electron repulsion integral evaluation in nonrelativistic molecular calculations.     

Hydrated Electron Transfer to Nucleobases in Aqueous Solutions Revealed by Ab Initio Molecular Dynamics Simulations

We present an ab initio molecular dynamics (AIMD) simulation study into the transfer dynamics of an excess electron from its cavity‐shaped hydrated electron state to a hydrated nucleobase (NB)‐bound state. In contrast to the traditional view that electron localization at NBs (G/A/C/T), which is the first step for electron‐induced DNA damage, is related only to dry or prehydrated electrons, and a fully hydrated electron no longer transfers to NBs, our AIMD simulations indicate that a fully hydrated electron can still transfer to NBs. We monitored the transfer dynamics of fully hydrated electrons towards hydrated NBs in aqueous solutions by using AIMD simulations and found that due to solution‐structure fluctuation and attraction of NBs, a fully hydrated electron can transfer to a NB gradually over time. Concurrently, the hydrated electron cavity gradually reorganizes, distorts, and even breaks. The transfer could be completed in about 120–200 fs in four aqueous NB solutions, depending on the electron‐binding ability of hydrated NBs and the structural fluctuation of the solution. The transferring electron resides in the π*‐type lowest unoccupied molecular orbital of the NB, which leads to a hydrated NB anion. Clearly, the observed transfer of hydrated electrons can be attributed to the strong electron‐binding ability of hydrated NBs over the hydrated electron cavity, which is the driving force, and the transfer dynamics is structure‐fluctuation controlled. This work provides new insights into the evolution dynamics of hydrated electrons and provides some helpful information for understanding the DNA‐damage mechanism in solution.     

Long-lived charge-transfer states in compact donor-acceptor dyads.

Photoinduced electron transfer is a widely applied method to convert photon energy into a useful (electro)chemical potential, both in nature and in artificial devices. There is a continuing effort to develop molecular systems in which the charge-transfer state, populated by photoinduced electron transfer, survives sufficiently long to tap the energy stored in it. In general this has been found to require the construction of rather complex molecular systems, but more recently a few approaches have been reported that allow the use of much more simple and relatively small electron donor-acceptor dyads for this purpose. The most successful examples of such systems seem to be those that apply "electron spin control" to slow down the spontaneous decay of the charge-transfer state, and these are reviewed in this minireview, with a discussion of the underlying principles and a critical evaluation of some of the claims made with regard to using a pronounced "inverted-region effect" as an alternative method to prolong the lifetime of charge-transfer states.     

Non-Condon效应对电子转移速率影响的理论模拟及其在有机半导体中的应用

随着Condon近似下各种电子转移理论的不断发展与完善和人们对non-Condon效应在电子转移过程中重要作用认识的逐步深入,已建立了几个理论模型来研究这种效应对电子转移速率的影响.本文主要总结了近两年来我们在non-Condon效应电子转移理论方面的工作,首先阐述了指数型、高斯型以及直线型non-Condon电子转移速率的全量子表达式,然后运用该理论模型以及分子动力学模拟计算了二噻吩四硫富瓦烯（DT-TTF）有机半导体的迁移率.此外,还进一步利用数值模拟详细研究了这三种线型的non-Condon效应在量子尺度上对电子转移速率的影响.     

The midpoint method for parallelization of particle simulations

The evaluation of interactions between nearby particles constitutes the majority of the computational workload involved in classical molecular dynamics (MD) simulations. In this paper, we introduce a new method for the parallelization of range-limited particle interactions that proves particularly suitable to MD applications. Because it applies not only to pairwise interactions but also to interactions involving three or more particles, the method can be used for evaluation of both nonbonded and bonded forces in a MD simulation. It requires less interprocessor data transfer than traditional spatial decomposition methods at all but the lowest levels of parallelism. It gains an additional practical advantage in certain commonly used interprocessor communication networks by distributing the communication burden more evenly across network links and by decreasing the associated latency. When used to parallelize MD, it further reduces communication requirements by allowing the computations associated with short-range nonbonded interactions, long-range electrostatics, bonded interactions, and particle migration to use much of the same communicated data. We also introduce certain variants of this method that can significantly improve the balance of computational load across processors.     

Allosteric Control of Naphthalene Diimide Encapsulation and Electron Transfer in Porphyrin Containers: Photophysical Studies and Molecular Dynamics Simulation

The complexation processes of N,N’-dibutyl-1,4,5,8-naphthalene diimide ( NDI ) into two types of π-electron-rich molecular containers consisting of two Zn(II)-porphyrins connected by four flexible linkers of two different lengths, were characterized by means of absorption and emission spectroscopies and molecular dynamics simulation. Notably, the addition of NDI leads to a strong quenching of the fluorescence of both cages only when they are in an open conformation suitable for guest encapsulation, a situation triggered by silver(I) ions binding to the lateral triazoles. Molecular dynamics simulations confirm the fast binding of NDI , likely assisted by NDI -silver(I) interactions. Upon NDI complexation, the two porphyrin macrocycles get closer, with an optimized face to face orientation, suggesting an induced-fit mechanism through π–π interactions with the NDI aromatic cycle. Ultrafast transient absorption experiments allowed to identify the process of quenching of the Zn-porphyrin fluorescence as an efficient photoinduced electron transfer reaction between the cage porphyrin and the included NDI guest. The process occurs on fast and ultrafast time scales in the two complexes (1.5 ps and ≤300 fs) leading to a short-lived charge separated state (charge recombination lifetimes in the order of 30–40 ps). The combined computational and experimental approach used here is able to furnish a reliable model of the NDI -cage complexation mechanism and of the corresponding electron transfer reaction, attesting the allosteric control of both processes by the silver(I) ions.     

Enveloping distribution sampling: a method to calculate free energy differences from a single simulation

The authors present a method to calculate free energy differences between two states A and B "on the fly" from a single molecular dynamics simulation of a reference state R. No computer time has to be spent on the simulation of intermediate states. Only one state is sampled, i.e., the reference state R which is designed such that the subset of phase space important to it is the union of the parts of phase space important to A and B. Therefore, an accurate estimate of the relative free energy can be obtained by construction. The authors applied the method to four test systems (dipole inversion, van der Waals interaction perturbation, charge inversion, and water to methanol conversion) and compared the results to thermodynamic integration estimates. In two cases, the enveloping distribution sampling calculation was straightforward. However, in the charge inversion and the water to methanol conversion, Hamiltonian replica-exchange molecular dynamics of the reference state was necessary to observe transitions in the reference state simulation between the parts of phase space important to A and B, respectively. This can be explained by the total absence of phase space overlap of A and B in these two cases.     

Charge carrier dynamics in phonon-induced fluctuation systems from time-dependent wavepacket diffusion approach

A time-dependent wavepacket diffusion method is proposed to deal with charge transport in organic crystals. The electron-phonon interactions in both site energies and electronic couplings are incorporated by the time-dependent fluctuations which are generated from the corresponding spectral density functions. The numerical demonstrations reveal that the present approach predicts the consistent charge carrier dynamics with the rigorous quantum approaches. In addition, the diffusion coefficients obtained from the Marcus formula are well reproduced at the weak electronic coupling and high temperature limits. It is also found that the charge mobility feature of the crossover from the band-like to the hopping-type cannot be predicted from the fluctuations induced by the linear electron-phonon interactions with an Ohmic spectral density; however, it indeed appears as the electronic coupling fluctuation exponentially depends on the nuclear coordinates. Finally, it should be noted that although the present approach neglects the imaginary fluctuation, it essentially incorporates the coherent motion of the charge carrier and quantum effect of the phonon motion with a broad regime of the fluctuations for symmetric systems. Besides, the approach can easily be applied to systems having thousands of sites, which allows one to investigate charge transport in nanoscale organic crystals.     

Photochemical initiation of electron transfer reactions.

It is quite apparent that the use of photoinitiated electron transfer has become a powerful, if not dominating, technique in the study of biological electron transfer. It provides a means to measure directly very fast processes and, through the choice of approach (flavin semiquinones or related, metal substitution in hemes or modification with ruthenium) and experimental conditions, provides the ability to probe different features of the electron transfer mechanism. Nevertheless, much remains to be done to fully understand biological electron transfer. The use of photoinitiated electron transfer has clearly established a role for a number of factors involved in controlling the kinetics of electron transfer, including driving force, distance, intervening media, dynamics (conformational gating) and orientation of redox centers. However, we have only scratched the surface in regard to understanding in molecular terms the details of electron transfer in physiologically relevant systems. Thus, even relatively simple and well characterized systems like cytochrome c-cytochrome c peroxidase remain obscure in terms of the through-protein electron paths (intervening media) and the role of protein dynamics in controlling electron transfer kinetics. Indeed, it is the through-protein paths and conformational gating that are unique to biological systems and provide nature with the capability of modulating electron transfer kinetics to optimize biological function. Of the techniques described here, the use of flavin semiquinones is clearly the least invasive in that there is no evidence that flavin semiquinones bind to or perturb physiologically relevant systems. However, this approach is constrained in that precise distances and orientations are not always known, and the range of driving forces available is limited. In contrast, metal substitution and ruthenation allow the positions of interacting redox centers to be reasonably well defined and can provide a very large range of driving force. This latter point is particularly important since it provides a means to discriminate between rate limiting electron transfer and conformational gating. Nevertheless, chemically modifying redox proteins runs the risk of structural and electrostatic alterations which can be difficult to detect but have profound effects on the redox kinetics. Moreover, the intrinsic protein dynamics can be affected, resulting, in the worst case, in changes in conformational gating which cannot be resolved from rate limiting electron transfer. Given the early stage of development of photo-initiated electron transfer, substantial progress can be expected in the next few years. No doubt new approaches will be developed and existing approaches further refined. Especially important, the theoretical basis for interpreting and understanding electron transfer will continue to evolve.(ABSTRACT TRUNCATED AT 400 WORDS)