Molecular dynamics simulations of D2O ice photodesorption

Molecular dynamics (MD) calculations have been performed to study the ultraviolet (UV) photodissociation of D(2)O in an amorphous D(2)O ice surface at 10, 20, 60, and 90 K, in order to investigate the influence of isotope effects on the photodesorption processes. As for H(2)O, the main processes after UV photodissociation are trapping and desorption of either fragments or D(2)O molecules. Trapping mainly takes place in the deeper monolayers of the ice, whereas desorption occurs in the uppermost layers. There are three desorption processes: D atom, OD radical, and D(2)O molecule photodesorption. D(2)O desorption takes places either by direct desorption of a recombined D(2)O molecule, or when an energetic D atom produced by photodissociation kicks a surrounding D(2)O molecule out of the surface by transferring part of its momentum. Desorption probabilities are calculated for photoexcitation of D(2)O in the top four monolayers and are compared quantitatively with those for H(2)O obtained from previous MD simulations of UV photodissociation of amorphous water ice at different ice temperatures [Arasa et al., J. Chem. Phys. 132, 184510 (2010)]. The main conclusions are the same, but the average D atom photodesorption probability is smaller than that of the H atom (by about a factor of 0.9) because D has lower kinetic energy than H, whereas the average OD radical photodesorption probability is larger than that of OH (by about a factor of 2.5-2.9 depending on ice temperature) because OD has higher translational energy than OH for every ice temperature studied. The average D(2)O photodesorption probability is larger than that of H(2)O (by about a factor of 1.4-2.3 depending on ice temperature), and this is entirely due to a larger contribution of the D(2)O kick-out mechanism. This is an isotope effect: the kick-out mechanism is more efficient for D(2)O ice, because the D atom formed after D(2)O photodissociation has a larger momentum than photogenerated H atoms from H(2)O, and D transfers momentum more easily to D(2)O than H to H(2)O. The total (OD + D(2)O) yield has been compared with experiments and the total (OH + H(2)O) yield from previous simulations. We find better agreement when we compare experimental yields with calculated yields for D(2)O ice than when we compare with calculated yields for H(2)O ice.     

Molecular dynamics simulations of cluster nucleation during inert gas condensation

Molecular dynamics simulations of vapor-phase nucleation of germanium in an argon atmosphere were performed and a unexpected channel of nucleation was observed. This channel, vapor-induced cluster splitting, is important for more refractory materials since the critical nucleus size can fall below the size of a dimer. As opposed to conventional direct vapor nucleation of the dimer, which occurs by three-body collisions, cluster-splitting nucleation is a second-order reaction. The most important cluster-splitting reaction is the collision of a vapor atom and a trimer that leads to the formation of two dimers. The importance of the cluster-splitting nucleation channel relative to the direct vapor nucleation channel is observed to increase with decreasing vapor density and increasing ratio of vapor to carrier gas atoms.     

Metadynamics simulations of ice nucleation and growth

The metadynamics method for accelerating rate events in molecular simulations is applied to the problem of ice freezing. We demonstrate homogeneous nucleation and growth of ice at 180 K in the isothermal-isobaric ensemble without the presence of external fields or surfaces. This result represents the first report of continuous and dynamic ice nucleation in a system of freely evolving density. Simulations are conducted using a variety of periodic simulation domains. In all cases the cubic polymorph ice I(c) is grown. The influence of boundary effects on estimates of the nucleation free energy barrier are discussed in relation to differences between this and earlier work.     

Molecular dynamics simulation of quasi-two-dimensional water clusters on ice nucleation protein

We performed molecular dynamics simulations of systems that consisted of the ice nucleation protein and the quasi-two-dimensional water cluster on it. The angle distributions, percolation probabilities, mean cluster sizes, cluster size distributions, and hydrogen bond relaxation times were analyzed. We concluded that the behavior of the water clusters on the ice nucleation protein was elaborately intertwined by the interaction between the ice nucleation protein and water, the interaction between the water molecules and the effect of temperature. The percolation probability and mean cluster size depended on the interactions and temperatures.     

Molecular dynamics simulations of methane hydrate using polarizable force fields

Molecular dynamics simulations of methane hydrate have been carried out using the polarizable AMOEBA and COS/G2 force fields. Properties calculated include the temperature dependence of the lattice constant, the OC and OO radial distribution functions, and the vibrational spectra. Both the AMOEBA and COS/G2 force fields are found to successfully account for the available experimental data, with overall somewhat better agreement with experiment being found for the AMOEBA model. Comparison is made with previous results obtained using TIP4P and SPC/E effective two-body force fields and the polarizable TIP4P-FQ force field, which allows for in-plane polarization only. Significant differences are found between the properties calculated using the TIP4P-FQ model and those obtained using the other models, indicating an inadequacy of restricting explicit polarization to in-plane only.     

Molecular dynamics simulations of nucleation from vapor to solid composed of Lennard-Jones molecules

We performed molecular dynamics (MD) simulations of nucleation from vapor at temperatures below the triple point for systems consisting of 10(4)-10(5) Lennard-Jones (L-J) type molecules in order to test nucleation theories at relatively low temperatures. Simulations are performed for a wide range of initial supersaturation ratio (S(0) ? 10-10(8)) and temperature (kT = 0.2-0.6ε), where ε and k are the depth of the L-J potential and the Boltzmann constant, respectively. Clusters are nucleated as supercooled liquid droplets because of their small size. Crystallization of the supercooled liquid nuclei is observed after their growth slows. The classical nucleation theory (CNT) significantly underestimates the nucleation rates (or the number density of critical clusters) in the low-T region. The semi-phenomenological (SP) model, which corrects the CNT prediction of the formation energy of clusters using the second virial coefficient of a vapor, reproduces the nucleation rate and the cluster size distributions with good accuracy in the low-T region, as well as in the higher-T cases considered in our previous study. The sticking probability of vapor molecules onto the clusters is also obtained in the present MD simulations. Using the obtained values of sticking probability in the SP model, we can further refine the accuracy of the SP model.     

Kinetics of NaCl nucleation in supercritical water investigated by molecular dynamics simulations

At system pressures between 17 MPa and 25 MPa the nucleation and growth of NaCl nanoparticles in water at supercritical conditions was investigated by molecular dynamics simulations at different system temperatures and system densities. Our results show that particle formation takes place within a few hundred picoseconds after the jump from ambient to supercritical conditions. After nucleation a phase of growth by adding monomers is followed by growth via cluster-cluster collisions. We present results on the time development of distributions of cluster sizes, cluster compositions, and cluster temperatures as well as radial distribution functions and nucleation rates.     

Molecular dynamics simulations of halide glasses

Summary Halide glasses have been extensively studied in recent years because of their potential application as infrared transmitting fibre optic materials. They are believed to be more ionic than glasses based on silica and should therefore be more amenable to molecular dynamics simulation using simple two-body potentials. Here the main features of structural models derived using such techniques are described and compared with available structural data. Possible future applications of this approach are outlined.     

Molecular dynamics simulations of peptide-surface interactions

Proteins, which are bioactive molecules, adsorb on implants placed in the body through complex and poorly understood mechanisms and directly influence biocompatibility. Molecular dynamics modeling using empirical force fields provides one of the most direct methods of theoretically analyzing the behavior of complex molecular systems and is well-suited for the simulation of protein adsorption behavior. To accurately simulate protein adsorption behavior, a force field must correctly represent the thermodynamic driving forces that govern peptide residue-surface interactions. However, since existing force fields were developed without specific consideration of protein-surface interactions, they may not accurately represent this type of molecular behavior. To address this concern, we developed a host-guest peptide adsorption model in the form of a G(4)-X-G(4) peptide (G is glycine, X is a variable residue) to enable determination of the contributions to adsorption free energy of different X residues when adsorbed to functionalized Au-alkanethiol self-assembled monolayers (SAMs). We have previously reported experimental results using surface plasmon resonance (SPR) spectroscopy to measure the free energy of peptide adsorption for this peptide model with X = G and K (lysine) on OH and COOH functionalized SAMs. The objectives of the present research were the development and assessment of methods to calculate adsorption free energy using molecular dynamics simulations with the GROMACS force field for these same peptide adsorption systems, with an oligoethylene oxide (OEG) functionalized SAM surface also being considered. By comparing simulation results to the experimental results, the accuracy of the selected force field to represent the behavior of these molecular systems can be evaluated. From our simulations, the G(4)-G-G(4) and G(4)-K-G(4) peptides showed minimal to no adsorption to the OH SAM surfaces and the G(4)-K-G(4) showed strong adsorption to the COOH SAM surface, which is in agreement with our SPR experiments. Contrary to our experimental results, however, the simulations predicted a relatively strong adsorption of G(4)-G-G(4) peptide to the COOH SAM surface. In addition, both peptides were unexpectedly predicted to adsorb to the OEG surface. These findings demonstrate the need for GROMACS force field parameters to be rebalanced for the simulation of peptide adsorption behavior on SAM surfaces. The developed methods provide a direct means of assessing, modifying, and validating force field performance for the simulation of peptide and protein adsorption to surfaces, without which little confidence can be placed in the simulation results that are generated with these types of systems.     

Molecular dynamics simulations of lipid bilayers

Computer simulation methods are becoming increasingly widespread as tools for studying the structure and dynamics of lipid bilayer membranes. The length scale and time scale accessible to atomic-level molecular dynamics simulations are rapidly increasing, providing insight into the relatively slow motions of molecular reorientation and translation and demonstrating that effects due to the finite size of the simulation cell can influence simulation results. Additionally, significant advances have been made in the complexity of membrane systems studied, including bilayers with cholesterol, small solute molecules, and lipid-protein and lipid-DNA complexes. Especially promising is the progress that continues to be made in the comparison of simulation results with experiment, both to validate the simulation algorithms and to aid in the interpretation of existing experimental data.     

Molecular dynamics simulations of polyelectrolyte adsorption

We have performed molecular dynamics simulations of polyelectrolyte adsorption at oppositely charged surfaces from dilute polyelectrolyte solutions. In our simulations, polyelectrolytes were modeled by chains of charged Lennard-Jones particles with explicit counterions. We have studied the effects of the surface charge density, surface charge distribution, solvent quality for the polymer backbone, strength of the short-range interactions between polymers and substrates on the polymer surface coverage, and the thickness of the adsorbed layer. The polymer surface coverage monotonically increases with increasing surface charge density for almost all studied systems except for the system of hydrophilic polyelectrolytes adsorbing at hydrophilic surfaces. In this case the polymer surface coverage saturates at high surface charge densities. This is due to additional monomer-monomer repulsion between adsorbed polymer chains, which becomes important in dense polymeric layers. These interactions also preclude surface overcharging by hydrophilic polyelectrolytes at high surface charge densities. The thickness of the adsorbed layer shows monotonic dependence on the surface charge density for the systems of hydrophobic polyelectrolytes for both hydrophobic and hydrophilic surfaces. Thickness is a decreasing function of the surface charge density in the case of hydrophilic surfaces while it increases with the surface charge density for hydrophobic substrates. Qualitatively different behavior is observed for the thickness of the adsorbed layer of hydrophilic polyelectrolytes at hydrophilic surfaces. In this case, thickness first decreases with increasing surface charge density, then it begins to increase.     

Molecular dynamics simulations of biotin carboxylase

Biotin carboxylase catalyzes the ATP-dependent carboxylation of biotin and is one component of the multienzyme complex acetyl-CoA carboxylase that catalyzes the first committed step in fatty acid synthesis in all organisms. Biotin carboxylase from Escherichia coli, whose crystal structures with and without ATP bound have been determined, has served as a model system for this component of the acetyl-CoA carboxylase complex. The two crystal structures revealed a large conformational change of one domain relative to the other domains when ATP is bound. Unfortunately, the crystal structure with ATP bound was obtained with an inactive site-directed mutant of the enzyme. As a consequence the structure with ATP bound lacked key structural information such as for the Mg2+ ions and contained altered conformations of key active-site residues. Therefore, nanosecond molecular dynamics studies of the wild-type biotin carboxylase were undertaken to supplant and amend the results of the crystal structures. Specifically, the protein-metal interactions of the two catalytically critical Mg2+ ions bound in the active site are presented along with a reevaluation of the conformations of active-site residues bound to ATP. In addition, the regions of the polypeptide chain that serve as hinges for the large conformational change were identified. The results of the hinge analysis complemented a covariance analysis that identified the individual structural elements of biotin carboxylase that change their conformation in response to ATP binding.     

Molecular dynamics simulations of cardiolipin bilayers

Cardiolipin is a key lipid component in the inner mitochondrial membrane, where the lipid is involved in energy production, cristae structure, and mechanisms in the apoptotic pathway. In this article we used molecular dynamics computer simulations to investigate cardiolipin and its effect on the structure of lipid bilayers. Three cardiolipin/POPC bilayers with different lipid compositions were simulated: 100, 9.2, and 0% cardiolipin. We found strong association of sodium counterions to the carbonyl groups of both lipid types, leaving in the case of 9.2% cardiolipin virtually no ions in the aqueous compartment. Although binding occurred primarily at the carbonyl position, there was a preference to bind to the carbonyl groups of cardiolipin. Ion binding and the small headgroup of cardiolipin gave a strong ordering of the hydrocarbon chains. We found significant effects in the water dipole orientation and water dipole potential which can compensate for the electrostatic repulsion that otherwise should force charged lipids apart. Several parameters relevant for the molecular structure of cardiolipin were calculated and compared with results from analyses of coarse-grained simulations and available X-ray structural data.     

Molecular dynamics simulations of lithium alkyl carbonates

A quantum chemistry-based many-body polarizable force field has been developed for two model solid-electrolyte interphase (SEI) components: dilithium ethylene dicarbonate (Li(2)EDC) and lithium methyl carbonate (LiMC). Molecular dynamics (MD) simulations of amorphous Li(2)EDC and LiMC were performed at temperatures from 393 to 600 K. Simulations reveal that Li(+) is coordinated by approximately 4.6 oxygen atoms from -COO(-) groups coming from different alkyl carbonate molecules. Charge transport in Li(2)EDC was found to be almost entirely due to Li(+). The temperature dependence of the ionic conductivity of the SEI model compounds Li(2)EDC and LiMC was found to be significantly stronger than that of liquid electrolytes (e.g., ethylene carbonate + LiTFSI), yielding extrapolated conductivities of the Li(2)EDC on order of 10(-10) S/cm at -30 degrees C.     

Molecular dynamics simulations of peptide carboxylate hydration

Aqueous solvation of carboxylate groups, as present in the glycine zwitterion and the dipeptide aspartylalanine, is studied employing a force-field that includes distributed multipole electrostatics and induction contributions (Amoebapro: P. Ren and J. W. Ponder, J. Comput. Chem., 2002, 23, 1497; P. Ren and J. W. Ponder, J. Phys. Chem. B, 2003, 107, 5933; J. W. Ponder and D. A. Case, Adv. Protein Chem., 2003, 66, 27). Radial and orientation distribution functions, as well as hydration numbers, are calculated and compared with existing simulation data derived from Car-Parrinello molecular dynamics (CPMD), and also distributed-charge force-fields. Connections are also made with experimental data for solvation of carboxylates in water. Our findings show that Amoebapro yields carboxylate solvation properties in very good agreement with CPMD results, significantly closer agreement than can be obtained from traditional force-fields. We also demonstrate that the influence of solvation on the conformation of the dipeptide is markedly different using Amoebapro compared with the other force-fields.     

Molecular dynamics simulations of supramolecular polymer rheology

Using equilibrium and nonequilibrium molecular dynamics simulations, we studied the equilibrium and rheological properties of dilute and semidilute solutions of head-to-tail associating polymers. In our simulation model, a spontaneous complementary reversible association between the donor and the acceptor groups at the ends of oligomers was achieved by introducing a combination of truncated pseudo-Coulombic attractive potential and Lennard Jones repulsive potential between donor, acceptor, and neighboring groups. We have calculated the equilibrium properties of supramolecular polymers, such as the ring/chain equilibrium, average molecular weight, and molecular weight distribution of self-assembled chains and rings, which all agree well with previous analytical and computer modeling results. We have investigated shear thinning of solutions of 8- and 20-bead associating oligomers with different association energies at different temperatures and oligomer volume fractions. All reduced viscosity data for a given oligomer length can be collapsed into one master curve, exhibiting two power-law regions of shear-thinning behavior with an exponent of -0.55 at intermediate ranges of the reduced shear rate β and -0.8 (or -0.9) at larger shear rates. The equilibrium viscosity of supramolecular solutions with different oligomer lengths and associating energies is found to obey a power-law scaling dependence on oligomer volume fraction with an exponent of 1.5, in agreement with the experimental observations for several dilute or semidilute solutions of supramolecular polymers. This implies that dilute and semidilute supramolecular polymer solutions exhibit high polydispersity but may not be sufficiently entangled to follow the reptation mechanism of relaxation.     

Molecular dynamics simulations applied to electric field induced second harmonic generation in dipolar chromophore solutions

Electric field induced second harmonic generation (EFISH) is an important experimental technique in extracting the first hyperpolarizability of an organic chromophore molecule. Such experiments are carried out in solutions with chromophore molecules dissolved in some common solvents. A known fact is that the first hyperpolarizabilities extracted from EFISH experiments are subject to the use of local field factors. In this work, we apply simulations to study the EFISH properties of chromophore solutions. By combining quantum chemistry calculations with the results derived from molecular dynamics simulations, we show how macroscopic EFISH properties can be modeled, using 4-(dimethylamino)-4'-nitroazobenzene dissolved in chloroform as a demonstration case. The focus of the study is on deriving accurate local field factors. We find that the local field approach applies very well to dipolar solutions, such as the one studied here, but that the local field factors derived are much smaller than the commonly used Onsager or Lorentz local field factors. Our study indicates that many of the reported first hyperpolarizabilities for dipolar molecules from EFISH experiments are most probably underestimated because the Onsager/Lorentz approach, commonly used in extracting the molecular first hyperpolarizability, neglects the effects of the shapes of dipolar chromophore molecules on the local field factors.     

Molecular dynamics simulations of local field factors

In the present work, the authors evaluate a scheme based on molecular dynamics to derive local field factors. These are given without any assumption of a cavity by fitting the Langevin functions to the order parameters obtained from the molecular dynamics simulations. The local field factors so obtained, with the detailed chromophore-solvent interactions and solvent structures taken into account, are much smaller than those calculated from the conventional Onsager and Lorentz models. A numerical demonstration is given for two typical organic chromophore molecules, p-nitroaniline and p-nitro-N,N-dimethylaniline dissolved in chloroform.     

Molecular dynamics study of electric polarisation

Non-equilibrium molecular dynamics is used to calculate the frequency-dependent dielectric function for a polar fluid. Absorption due to rotational resonance, and strong power absorption due to vibrational resonance, are observed. The former is partially accounted for in terms of a memory function model. Non-linear dielectric saturation at zero frequency is well described by a Langevin function.