scholarly journals Polarizable force field for molecular dynamics simulations of silver nanoparticles

Contact of silver metal surfaces with water, ions and organic ligands experiences induced charges, leading to attractive polarization. These forces play an important role at inorganic/organic interfaces and complement other non-bonded surface interactions. Despite the importance of these interactions, it, however, remains difficult to implement polarization effects to classical molecular dynamics (MD) simulations. In this contribution, we first present an overview of two popular polarizable models, such as Drude oscillator and the rigid rod model, which are utilized to mimic the polarizability of bulk metals. Second, we implemented the rigid rod model to the polarizable force field (FF) for a silver atom, which was further adapted for atomistic MD simulations of silver nanoparticles (AgNPs) composed of 1397 atoms. In our model, induced charge polarization is represented by the displacement of a charge-carrying virtual site attached rigidly to an original Ag atom. To explore the role of polarization, we compared the performance of the classical nonpolarizable FF and the new polarizable model in the MD simulations of adsorption of water and ions onto quasi-spherical AgNP and the flat crystalline silver surface. The analysis of the radial distribution function of Ag-Ag atoms demonstrated that the introduction of the polarization effect had minor effects on face-centered cubic (fcc) packing of silver atoms of bare and water-solvated AgNPs. We found that the polarizable FF causes some increase in attractive interactions between the silver surface and water molecules and Na+ ions. As a crucial test of the developed polarizable model, the structure of adsorbed interfacial water molecules was analyzed. Our data suggest that the environment-induced polarization of the silver surface contributes significantly to the structure of adsorbed interfacial water layers and it also plays an important role in the adsorption of positive ions. However, it was also found out that the polarization effect has a rather short-range effect, so that a minor contribution of silver polarization was seen for adsorption of water molecules and ions from distant solvation shells.

1999 ◽  
Vol 54 (11) ◽  
pp. 896-902 ◽  
Author(s):  
Antonio Matas ◽  
Antonio Heredia

Abstract A theoretical molecular modelling study has been conducted for cutin, the biopolyester that forms the main structural component of the plant cuticle. Molecular dynamics (MD) simulations, extended over several ten picoseconds, suggests that cutin is a moderately flexible netting with motional constraints mainly located at the cross-link sites of functional ester groups. This study also gives structural information essentially in accordance with previously reported experimental data, obtained from X -ray diffraction and nuclear magnetic resonance experiments. MD calculations were also performed to simulate the diffusion of water mole­cules through the cutin biopolymer. The theoretical analysis gives evidence that water perme­ation proceedes by a “hopping mechanism”. Coefficients for the diffusion of the water molecules in cutin were obtained from their mean-square displacements yielding values in good agreement with experimental data.


2020 ◽  
Vol 73 (8) ◽  
pp. 734
Author(s):  
W. A. Monika Madhavi ◽  
Samantha Weerasinghe ◽  
Konstantin I. Momot

Rotational motion of molecules plays an important role in determining NMR spin relaxation properties of liquids. The textbook theory of NMR spin relaxation predominantly uses the assumption that the reorientational dynamics of molecules is described by a continuous time rotational diffusion random walk with a single rotational diffusion coefficient. Previously we and others have shown that reorientation of water molecules on the timescales of picoseconds is not consistent with the Debye rotational-diffusion model. In particular, multiple timescales of molecular reorientation were observed in liquid water. This was attributed to the hydrogen bonding network in water and the consequent presence of collective rearrangements of the molecular network. In order to better understand the origins of the complex reorientational behaviour of water molecules, we carried out molecular dynamics (MD) simulations of a liquid that has a similar molecular geometry to water but does not form hydrogen bonds: hydrogen sulfide. These simulations were carried out at T=208K and p=1 atm (~5K below the boiling point). Ensemble-averaged Legendre polynomial functions of hydrogen sulfide exhibited a Gaussian decay on the sub-picosecond timescale but, unlike water, did not exhibit oscillatory behaviour. We attribute these differences to hydrogen sulfide’s absence of hydrogen bonding.


2008 ◽  
Vol 128 (17) ◽  
pp. 174705 ◽  
Author(s):  
Meng-Hsiung Weng ◽  
Wen-Jay Lee ◽  
Shin-Pon Ju ◽  
Chien-Hsiang Chao ◽  
Nan-Kai Hsieh ◽  
...  

Author(s):  
Fabian Keller ◽  
Andreas Heuer ◽  
Hans-Joachim Galla ◽  
Jens Smiatek

The interactions between DPPC lipid bilayers in different phases with ectoine, amino ectoine and water molecules are studied by means of atomistic molecular dynamics (MD) simulations and conceptual density functional theory (DFT) calculations.


2017 ◽  
Vol 19 (11) ◽  
pp. 7476-7480 ◽  
Author(s):  
Tian Sheng ◽  
Jin-Yu Ye ◽  
Wen-Feng Lin ◽  
Shi-Gang Sun

In this work, we have studied methanol oxidation mechanisms on RuO2(100) by using density functional theory (DFT) calculations and ab initio molecular dynamics (MD) simulations with some explicit interfacial water molecules.


2020 ◽  
Vol 4 (1) ◽  
pp. 20
Author(s):  
Thi Tuong Vy Phan

The voltage-gated proton channel (Hv1) plays the important role in proton extrusion, pH homeostasis, sperm motility, and cancer progression. The closed-state structure of Hv1 was recently revealed by X-ray crystallography. However, the opened-state structure has not been captured yet. To investigate the mechanism of proton transfer in Hv1, molecular dynamics (MD) simulations were performed with the closed-state structure of Hv1 under electric field and pH conditions. The residues arrangement on the closed-state structure revealed that the selectivity filter (Asp108) which is located in the hydrophobic layer (consists of two Phe residues 146 and 179) might prevent water penetration. In molecular dynamics simulations, we observed that the channel opened by moving 3 Arg up on the S4 helix and a continuous hydrogen-bonded chain of water molecules (a “water wire”) went through the channel when it opened. During simulations, the open channel allowed water molecules to pass through the channel but excluded other ions. This indicates the Hv1 channel is highly selective for protons. Our results clearly showed the Hv1 channel is voltage-and pH-gradient sensing.


2011 ◽  
Vol 233-235 ◽  
pp. 1872-1877 ◽  
Author(s):  
Run Liang Zhu ◽  
Thomas V. Shapley ◽  
Marco Molinari ◽  
Ge Fei ◽  
Stephen C. Parker

Molecular dynamics (MD) simulations have been used to investigate the interlayer structure of water saturated organoclays. The basal spacing values of cetyltrimethylammonium (CTMA) intercalated montmorillonite (CTMA-Mont) in dry and water saturated states were detected using XRD. Then the results were compared with simulation results of dry CTMA-Mont. The MD simulations show that the CTMA cations form layer structures on siloxane surface and aggregate in the interlayer space. Water molecules can access part of the siloxane surface and form H-bonds with surface oxygen atoms by donating one or two of the hydrogen atoms. Thus, the water molecules close to the surface have a preferred orientation with the dipole pointing towards the surface, while in the interlayer space, the water molecules aggregate to form large clusters. The H-bonds between surface oxygen and water molecules are shown to be slightly weaker than those between water molecules. Although water molecules within interlayer space can form strong H-bonds as in bulk water the number of H-bond for each water molecule is reduced. Our results indicate that MD simulations represent a useful tool for exploring the microstructure of water saturated organoclays.


2018 ◽  
Vol 209 ◽  
pp. 249-257 ◽  
Author(s):  
Rob D. Coalson

Atomistic Molecular Dynamics (MD) simulations provide numerous insights into the process whereby water is driven through a narrow nanopore (diameter on the order of a few water molecules) by application of hydrostatic pressure.


Sign in / Sign up

Export Citation Format

Share Document