Atomistic simulations of calcite nanoparticles and their interaction with water

Molecular dynamics (MD) simulations have been used to study the stability of calcite nanoparticles ranging in size from 18 to 324 f.u., both in vacuo and in the presence of explicit water molecules. In vacuo, the smallest particles become highly disordered during the MD simulation due to rotation and translation of the undercoordinated CO(3) (2-) anions at the edges of the particles. As the nanoparticle size increases, the influence of the fully coordinated bulk ions begins to dominate and long-range order is seen both in the Ca-C pair distribution functions and in the degree of rotational order of the CO(3) (2-) anions. However, when explicit water is added to the system, the molecules in the first hydration layer complete the coordination shell of the surface ions, preserving structural order even in the smallest of the nanoparticles. Close to particle surface, the structure of the water itself shows features similar to those seen close to planar periodic (1014) surfaces, although the molecules are far less tightly bound.     

A density-functional approach to polarizable models: a Kim-Gordon response density interaction potential for molecular simulations

A combined linear-response-frozen electron-density model has been implemented in a molecular-dynamics scheme derived from an extended Lagrangian formalism. This approach is based on a partition of the electronic charge distribution into a frozen region described by Kim-Gordon theory [J. Chem. Phys. 56, 3122 (1972); J. Chem. Phys. 60, 1842 (1974)] and a response contribution determined by the instantaneous ionic configuration of the system. The method is free from empirical pair potentials and the parametrization protocol involves only calculations on properly chosen subsystems. We apply this method to a series of alkali halides in different physical phases and are able to reproduce experimental structural and thermodynamic properties with an accuracy comparable to Kohn-Sham density-functional calculations.     

Liquid structures of water, methanol, and hydrogen fluoride at ambient conditions from first principles molecular dynamics simulations with a dispersion corrected density functional

Using first principles molecular dynamics simulations in the isobaric-isothermal ensemble (T = 300 K, p = 1 atm) with the Becke-Lee-Yang-Parr exchange/correlation functional and a dispersion correction due to Grimme, the hydrogen bonding networks of pure liquid water, methanol, and hydrogen fluoride are probed. Although an accurate density is found for water with this level of electronic structure theory, the average liquid densities for both hydrogen fluoride and methanol are overpredicted by 50 and 25%, respectively. The radial distribution functions indicate somewhat overstructured liquid phases for all three compounds. The number of hydrogen bonds per molecule in water is about twice as high as for methanol and hydrogen fluoride, though the ratio of cohesive energy over number of hydrogen bonds is lower for water. An analysis of the hydrogen-bonded aggregates revealed the presence of mostly linear chains in both hydrogen fluoride and methanol, with a few stable rings and chains spanning the simulation box in the case of hydrogen fluoride. Only an extremely small fraction of smaller clusters was found for water, indicating that its hydrogen bond network is significantly more extensive. A special form of water with on average about two hydrogen bonds per molecule yields a hydrogen-bonding environment significantly different from the other two compounds.     

Atomistic molecular dynamics simulations on the interaction of TEMPO-oxidized cellulose nanofibrils in water

Atomistic molecular dynamics simulations were carried out to obtain information on the rheological, aggregation and disintegration properties of carboxylated (TEMPO-oxidized) cellulose nanofibrils with different functionalization levels. The magnitude of the inter-fibril interaction was quantified for parallel nanofibrils using the umbrella sampling method. The obtained potential of mean force was found highly sensitive to the charge configuration for intermediate functionalization levels. This feature was further studied with an electrostatic model for similar charge configurations and system periodicity as in the case of the molecular dynamics simulations. The electrostatic contribution of the charged surfaces varied from repulsive to attractive depending on the distribution of the carboxylate groups and nearby counter-ions, as well as the distance between the fibrils. The simulated deviations from average behavior for single fibrils in both models suggest heterogeneity in their aggregation and disintegration behavior. This was seen in disintegration experiments, where the differences in disintegration energy and in the structural variation qualitatively agreed with the model predictions. As to aggregation behavior, the studied case with parallel fibrils reflects the upper boundary of the repulsive interaction.     

Theoretical study of intramolecular interaction energies during dynamics simulations of oligopeptides by the fragment molecular orbital-Hamiltonian algorithm method

We analyzed the interaction energies between residues (fragments) in an oligopeptide occurring during dynamic simulations by using the fragment molecular orbital-Hamiltonian algorithm (FMO-HA) method, an ab initio MO-molecular dynamics technique. The FMO method enables not only calculation of large molecules based on ab initio MO but also easy evaluation of interfragment interaction energies. The glycine pentamer [(Gly)(5)] and decamer [(Gly)(10)] were divided into five and ten fragments, respectively. alpha-helix structures of (Gly)(5) and (Gly)(10) were stabilized by attractive interaction energies owing to intramolecular hydrogen bonds between fragments n and n+3 (and n+4), and beta-strand structures were characterized by repulsive interaction energies between fragments n and n+2. We analyzed interfragment interaction energies during dynamics simulations as the peptides' geometries changed from alpha helix to beta strand. Intramolecular hydrogen bonds between fragments 2-4 and 2-5 control the geometrical preference of (Gly)(5) for the beta-strand structure. The pitch of one turn of the alpha helix of (Gly)(10) elongated and thus weakened during dynamics due to a shifting of the intramolecular hydrogen bonds, and enabled the beta-strand structure to form. Changes in interaction energies due to the intramolecular hydrogen bonds controlled the tendency toward alpha-helix or beta-strand structure of (Gly)(5) and (Gly)(10). Evaluation of interfragment interaction energies during dynamics simulations thus enabled detailed analysis of the process of the geometrical changes occurring in oligopeptides.     

Pair interaction potentials with explicit polarization for molecular dynamics simulations of La(3+) in bulk water

Pair interaction potentials (IPs) were defined to describe the La(3+)-OH(2) interaction for simulating the La(3+) hydration in aqueous solution. La(3+)-OH(2) IPs are taken from the literature or parametrized essentially to reproduce ab initio calculations at the second-order Moller-Plesset level of theory on La(H(2)O)(8) (3+). The IPs are compared and used with molecular dynamics (MD) including explicit polarization, periodic boundary conditions of La(H(2)O)(216) (3+) boxes, and TIP3P water model modified to include explicit polarization. As expected, explicit polarization is crucial for obtaining both correct La-O distances (r(La-O)) and La(3+) coordination number (CN). Including polarization also modifies hydration structure up to the second hydration shell and decreases the number of water exchanges between the La(3+) first and second hydration shells. r(La-O) ((1))=2.52 A and CN((1))=9.02 are obtained here for our best potential. These values are in good agreement with experimental data. The tested La-O IPs appear to essentially account for the La-O short distance repulsion. As a consequence, we propose that most of the multibody effects are correctly described by the explicit polarization contributions even in the first La(3+) hydration shell. The MD simulation results are slightly improved by adding a-typically negative 1r(6)-slightly attractive contribution to the-typically exponential-repulsive term of the La-O IP. Mean residence times are obtained from MD simulations for a water molecule in the first (1082 ps) and second (7.6 ps) hydration shells of La(3+). The corresponding water exchange is a concerted mechanism: a water molecule leaving La(H(2)O)(9) (3+) in the opposite direction to the incoming water molecule. La(H(2)O)(9) (3+) has a slightly distorded "6+3" tricapped trigonal prism D(3h) structure, and the weakest bonding is in the medium triangle, where water exchanges take place.     

A configuration analysis for fragment interaction

We propose a modified version of configuration analysis (CA) for the fragment interaction in conjunction with Kitaura’s fragment molecular orbital (FMO) scheme. The proposal is abbreviated as CAFI. The MO sets of fragments are merged and then orthonormalized by the use of a weighted Löwdin orthonormalization. The energy calculation is performed with the concurrent electron relaxation functional (CERF). The relaxation energy is obtained in an orbital-wise fashion and is distinguished as the charge-transfer and the polarization. The utility of CAFI is demonstrated through test calculations on hydrogen-bonding systems.     

Computational simulations of the conformational behaviour of the adhesive proteins RGDS fragment

Summary Many adhesive proteins present in extracellular matrices and in blood contain the tetrapeptide sequence -Arg-Gly-Asp-Ser- (or RGDS) at their cell recognition site. Since this sequence, or similar ones, was found in many proteins involved in major biological mechanisms, conformational investigations were performed on the RGDS fragment. A preliminary review of available crystal structures indicates that the RxDy sequences exhibit 3 well-defined structural patterns: one corresponding to a strong interaction between the Arg and Asp ionic side chains which are only about 4 Å apart, one with the ions separated by about 8 Å, and another in which the side chains are further apart (about 11 Å).The conformational behaviour of the isolated RGDS fragment was next tackled using sequential building, Monte Carlo and molecular dynamics computational techniques. Analysis of the RGDS sequence conformational possibilities, as simulated in vacuum and in water solution, indicates that they can be classified into several conformational classes, which correspond roughly to the behaviour of the RGDS fragment as observed in protein matrices. This suggests the possibility of understanding the biological role of the RGDS or parent sequences in recognition processes.     

Polarizable interaction potential for molecular dynamics simulations of actinoids(III) in liquid water

In this work, we have developed a polarizable classical interaction potential to study actinoids(III) in liquid water. This potential has the same analytical form as was recently used for lanthanoid(III) hydration [M. Duvail, P. Vitorge, and R. Spezia, J. Chem. Phys. 130, 104501 (2009)]. The hydration structure obtained with this potential is in good agreement with the experimentally measured ion-water distances and coordination numbers for the first half of the actinoid series. In particular, the almost linearly decreasing water-ion distance found experimentally is replicated within the calculations, in agreement with the actinoid contraction behavior. We also studied the hydration of the last part of the series, for which no structural experimental data are available, which allows us to provide some predictive insights on these ions. In particular we found that the ion-water distance decreases almost linearly across the series with a smooth decrease of coordination number from nine to eight at the end.     

Towards an assessment of the accuracy of density functional theory for first principles simulations of water. II

A series of 20 ps ab initio molecular dynamics simulations of water at ambient density and temperatures ranging from 300 to 450 K are presented. Car-Parrinello (CP) and Born-Oppenheimer (BO) molecular dynamics techniques are compared for systems containing 54 and 64 water molecules. At 300 K, an excellent agreement is found between radial distribution functions (RDFs) obtained with BO and CP dynamics, provided an appropriately small value of the fictitious mass parameter is used in the CP simulation. However, we find that the diffusion coefficients computed from CP dynamics are approximately two times larger than those obtained with BO simulations for T>400 K, where statistically meaningful comparisons can be made. Overall, both BO and CP dynamics at 300 K yield overstructured RDFs and slow diffusion as compared to experiment. In order to understand these discrepancies, the effect of proton quantum motion is investigated with the use of empirical interaction potentials. We find that proton quantum effects may have a larger impact than previously thought on structure and diffusion of the liquid.     

Analytic gradient and molecular dynamics simulations using the fragment molecular orbital method combined with effective potentials

The completely analytic energy gradients are derived and implemented for the two-body fragment molecular orbital (FMO2) method combined with the model core potentials (MCP) and effective fragment potentials (EFP). The many-body terms in EFP require solving coupled-perturbed Hartree-Fock equations, which are derived and implemented. The molecular dynamics (MD) simulations are performed using the FMO2/MCP method for the capped alanine decamer and with the FMO2/EFP method for the zwitterionic conformer of glycine tetramer immersed in the water layer of 6.0 Å (135 water molecules). The results of the MD simulations using the FMO2/EFP and FMO2/MCP gradients show that the total energy is conserved at the time steps less than 1 fs.     

Charge transfer interaction in the effective fragment potential method

An approximate formula is derived and implemented in the general effective fragment potential (EFP2) method to model the intermolecular charge transfer interaction. This formula is based on second order intermolecular perturbation theory and utilizes canonical molecular orbitals and Fock matrices obtained with preparative self-consistent field calculations. It predicts charge transfer energies that are in reasonable agreement with the reduced variational space energy decomposition analysis. The formulas for the charge transfer gradients with respect to EFP translational and rotational displacements are also derived and implemented.     

Chemiluminescence in the interaction of organoboron peroxides with water

Chemiluminescence (CL) was detected in the reactions of organoboron peroxides, intermediate products of the oxidation of 9-borabicyclo[3.3.1]nonane (1) by water. Triplet cyclooctanone emitting at λ=420±15 nm is the emitter of CL. A mechanism for the formation of the excited ketone was suggested. The interaction of water with other components of the reaction mixture (the initial boron hydride I, end products of its oxidation, and a number of intermediate radicals) are not accompanied by CL. Translated fromIzvestiya Akademii Nauk. Seriya Khimicheskaya, No. 9, pp. 1692–1695, September, 1997.     

Implementation of the analytic energy gradient for the combined time-dependent density functional theory/effective fragment potential method: application to excited-state molecular dynamics simulations

Excited-state quantum mechanics/molecular mechanics molecular dynamics simulations are performed, to examine the solvent effects on the fluorescence spectra of aqueous formaldehyde. For that purpose, the analytical energy gradient has been derived and implemented for the linear-response time-dependent density functional theory (TDDFT) combined with the effective fragment potential (EFP) method. The EFP method is an efficient ab initio based polarizable model that describes the explicit solvent effects on electronic excitations, in the present work within a hybrid TDDFT/EFP scheme. The new method is applied to the excited-state MD of aqueous formaldehyde in the n-π* state. The calculated π*→n transition energy and solvatochromic shift are in good agreement with other theoretical results.     

Fitting properties from density functional theory based molecular dynamics simulations to parameterize a rigid water force field

In the quest towards coarse-grained potentials and new water models, we present an extension of the force matching technique to parameterize an all-atom force field for rigid water. The methodology presented here allows to improve the matching procedure by first optimizing the weighting exponents present in the objective function. A new gauge for unambiguously evaluating the quality of the fit has been introduced; it is based on the root mean square difference of the distributions of target properties between reference data and fitted potentials. Four rigid water models have been parameterized; the matching procedure has been used to assess the role of the ghost atom in TIP4P-like models and of electrostatic damping. In the former case, burying the negative charge inside the molecule allows to fit better the torques. In the latter, since short-range interactions are damped, a better fit of the forces is obtained. Overall, the best performing model is the one with a ghost atom and with electrostatic damping. The approach shown in this paper is of general validity and could be applied to any matching algorithm and to any level of coarse graining, also for non-rigid molecules.     

Many-body interaction in glycine-(water)3 complex using density functional theory method

Various configurations were investigated to find the most stable structures of glycine-(water)3 complex. Five different optimized conformers of glycine-(water)3 complex are obtained from density functional theory calculations using 6-311++G* basis set. Relaxation energy and many body interaction energies (two, three, and four body) are also calculated for these conformers. Out of the five conformers, the most stable conformer has the BSSE corrected total energy -513.917 967 7 Hartree and binding energy -27.28 Kcal/mol. It has been found that the relaxation energies, two body energies and three body energies have significant contribution to the total binding energy whereas four body energies are very small. The chemical hardness and chemical potential also confirmed the stability of the conformer having lowest total energy.     

Auger decay calculations with core-hole excited-state molecular-dynamics simulations of water

We report a new theoretical procedure for calculating Auger decay transition rates including effects of core-hole excited-state dynamics. Our procedure was applied to the normal and first resonant Auger processes of gas-phase water and compared to high-resolution experiments. In the normal Auger decay, calculated Auger spectra were found to be insensitive to the dynamics, while the repulsive character of the first resonant core-excited state makes the first resonantly excited Auger decay spectra depend strongly on the dynamics. The ultrafast dissociation of water upon O(1s)-->4a(1) excitation was analyzed and found to be very sensitive to initial vibrational distortions in the ground state which furthermore affect the excitation energy. Our calculated spectra reproduce the experimental Auger spectra except for the Franck-Condon vibrational structure which is not included in the procedure. We found that the Auger decay of OH and O fragments contributes to the total intensity, and that the contribution from these fragments increases with increasing excitation energy.     

Effective fragment potential study of the interaction of DNA bases

Hydrogen-bonded and stacked structures of adenine-thymine and guanine-cytosine nucleotide base pairs, along with their methylated analogues, are examined with the ab inito based general effective fragment potential (EFP2) method. A comparison of coupled cluster with single, double, and perturbative triple (CCSD(T)) energies is presented, along with an EFP2 energy decomposition to illustrate the components of the interaction energy.