Temperature effects on spreading of water nano‐droplet on poly(methyl methacrylate): A molecular dynamics simulation study

In this article, the effect of temperature on the spreading behavior of a water nano‐droplet on poly(methyl methacrylate) substrate is investigated. The contact angle analysis illustrates that the spreading process occurs in a stage‐like manner and the increase in temperature causes a regime change from partial to total wetting. The interaction energy distributions show that there exist sites on the surface which could trap water molecules and provide a better path for other molecules to overcome the asperities. Estimations of the coefficients of self‐diffusivity suggest that temperature has a major effect in the reorientation stage, which results in the formation of the interfacial layer. In the second stage of spreading, temperature affects the process by providing sufficient energy for water molecules to overcome the interactions with the substrate. Therefore, this stage is controlled by the movement of water molecules on the surface and is highly influenced by their interaction with the surface asperities, strong interaction sites, and the carbonyl groups. © 2017 Wiley Periodicals, Inc. J. Polym. Sci., Part B: Polym. Phys. 2017 , 55, 1532–1541     

Tetracyclines adsorption onto alumina: A comparative experimental and molecular dynamics simulation study

Using alumina (Al₂O₃) as the adsorbent, a static adsorption experiment was carried out in this study. It comprehensively evaluated the factors including Al₂O₃ dosage, adsorption temperature, and pH that influence the adsorption capability of three tetracyclines (TCs), namely, tetracycline hydrochloride (TC), chlortetracycline hydrochloride (CTC) and oxytetracycline hydrochloride (OTC). The results demonstrate that the adsorption efficiency increases with Al₂O₃ dosage. In addition, low-acid or natural solution is benefit for the adsorption. The adsorption behavior is more reasonably described with the Freundlich isotherm, and fits well with the pseudo-second-order kinetic model (R²?>?0.999). The results of molecular dynamics (MD) simulation show that the structures of TCs deformed during the combining process. The values of binding energy of TCs follow the order as: CTC (88.45?kcal/mol)?>?OTC (73.54?kcal/mol)?>?TC (54.28?kcal/mol). The MD simulation results agree well with the adsorption experimental results, which indicates that the MD simulation is reliable and reasonable. The MD simulation will provide theoretical knowledge in understanding the adsorption mechanism and environmental behavior of TCs.     

Thermal analysis via molecular dynamics simulation

Thermal analysis by classical molecular dynamics simulations is discussed on hand of heat capacity of crystals of 9600 atoms. The differences between quantum mechanical and classical mechanical calculations are shown. Anharmonicity is proven to be an important factor. Finally, it is found that defects contribute to an increase in heat capacity before melting. The energy of conformational gauche defects within the crystal is only about 10% due to internal rotation. The other energy must be generated by cooperative strain. The conclusion is that the next generation of faster computers may permit wider use of molecular dynamics simulations in support of the interpretation of thermal analysis.Dedicated to Professor Bernhard Wunderlich on the occasion of his 65th birthday     

Nano‐adhesion and friction of multi‐asperity contact: a molecular dynamics simulation study

In our study, the contact and sliding processes between a flat plate and a substrate with multiple asperities are studied by molecular dynamics (MD) simulations, and how the number of asperities and asperity height influence the adhesion force and friction force are investigated thoroughly. The normal force versus the separation distance curve during contact processes is analyzed completely and from which the van der Waals (vdW) force (F_vdW) and the adhesion force (F_adh) are obtained and compared with the Katainen model. The adhesion force and the friction force increase linearly as the increase of the number of asperities (i.e. real contact area) with same asperity height. With the identical number of asperities, the adhesion force and the friction force decrease with the increase of the asperity height at first. However the reductions of the adhesion force and the friction force become less obvious, when the asperity height is larger than a critical value (20 Å for our simulation parameters). Copyright © 2015 John Wiley & Sons, Ltd.     

IBIsCO: a molecular dynamics simulation package for coarse-grained simulation

IBIsCO is a parallel molecular dynamics simulation package developed specially for coarse-grained simulations with numerical potentials derived by the iterative Boltzmann inversion (IBI) method (Reith et al., J Comput Chem 2003, 24, 1624). In addition to common features of molecular dynamics programs, the techniques of dissipative particle dynamics (Groot and Warren, J Chem Phys 1997, 107, 4423) and Lowe-Andersen dynamics (Lowe, Europhys Lett 1999, 47, 145) are implemented, which can be used both as thermostats and as sources of friction to compensate the loss of degrees of freedom by coarse-graining. The reverse nonequilibrium molecular dynamics simulation method (Müller-Plathe, Phys Rev E 1999, 59, 4894) for the calculation of viscosities is also implemented. Details of the algorithms, functionalities, implementation, user interfaces, and file formats are described. The code is parallelized using PE_MPI on PowerPC architecture. The execution time scales satisfactorily with the number of processors.     

Essential dynamics/factor analysis for the interpretation of molecular dynamics trajectories

Subject of this work is the analysis of molecular dynamics (MD) trajectories of neurophysins I (NPI) and II (NPII) and their complexes with the neurophyseal nonapeptide hormones oxytocin (OT) and vasopresssin (VP), respectively, simulated in water. NPs serve in the neurosecretory granules as carrier proteins for the hormones before their release to the blood. The starting data consisted of two pairs of different trajectories for each of the (NPII/VP)2 and (NPI/OT)2 heterotetramers and two more trajectories for the NPII2 and NPI2 homodimers (six trajectories in total). Using essential dynamics which, to our judgement, is equivalent to factor analysis, we found that only about 10 degrees of freedom per trajectory are necessary and sufficient to describe in full the motions relevant for the function of the protein. This is consistent with these motions to explain about 90% of the total variance of the system. These principal degrees of freedom represent slow anharmonic motional modes, clearly pointing at distinguished mobility of the atoms involved in the protein's functionality.     

Atomistic detailed free‐energy landscape of intrinsically disordered protein studied by multi‐scale divide‐and‐conquer molecular dynamics simulation

Calcineurin (CaN) is a eukaryotic serine/threonine protein phosphatase activated by both Ca²⁺ and calmodulin (CaM), including intrinsically disordered region (IDR). The region undergoes folding into an α‐helix form in the presence Ca²⁺‐loaded CaM. To sample the ordered structure of the IDR by conventional all atom model (AAM) molecular dynamics (MD) simulation, the IDR and Ca²⁺‐loaded CaM must be simultaneously treated. However, it is time‐consuming task because the coupled folding and binding should include repeated binding and dissociation. Then, in this study, we propose novel multi‐scale divide‐and‐conquer MD (MSDC‐MD), which combines AAM‐MD and coarse‐grained model MD (CGM‐MD). To speed up the conformation sampling, MSDC‐MD simulation first treats the IDR by CGM to sample conformations from wide conformation space; then, multiple AAM‐MD in a limited area is initiated using the resultant CGM conformation, which is reconstructed by homology modeling method. To investigate performance, we sampled the ordered conformation of the IDR using MSDC‐MD; the root‐mean‐square distance (RMSD) with respect to the experimental structure was 2.23 Å.     

Communication performance of d -meshes in molecular dynamics simulation

Communication algorithms, tailored for molecular dynamics simulation on d-meshes, are evaluated in terms of communication efficiency. It has been shown elsewhere that d-meshes are better than other regular topologies, e.g., hypercubes and standard toroidal 4-meshes, when compared in their diameter and average distance among nodes. Collective communication is needed in molecular dynamics simulation for the distribution of coordinates and calculation and distribution of new energies. We show that both collective communication patterns used in molecular dynamics can be efficiently solved with congestion-free algorithms for all-to-all communication based on store-and-forward routing and routing tables. Our results indicate that d-meshes compete with hypercubes in parallel computers. Therefore d-meshes can also be used as a communication upgrade of existing molecular dynamics simulation platforms and can be successfully applied to perform fast molecular dynamics simulation.     

A single‐molecule reaction cascade: First‐principles molecular dynamics simulation

The success of mechanochemistry is continued with the targeted organic synthesis of functional nano‐scale devices. In the present theoretical study, first principles molecular dynamics simulations are performed for a recently synthesized three‐ring system that functions as a molecular reaction cascade. Mechanochemical and photochemical reaction conditions are investigated. The system was designed in a way that three bonds would break consecutively in mechanochemical and sonochemical setups. We succeeded to simulate the reaction mechanisms with first‐principles molecular dynamics simulations and discuss the stereochemistry.     

Chemical reaction dynamics of PeCB and TCDD decomposition: A tight‐binding quantum chemical molecular dynamics study with first‐principles parameterization

The decomposition reaction dynamics of 2,3,4,4′,5‐penta‐chlorinated biphenyl (2,3,4,4′,5‐PeCB), 3,3′,4,4′,5‐penta‐chlorinated biphenyl (3,3′,4,4′,5‐PeCB), and 2,3,7,8‐tetra‐chlorinated dibenzo‐p‐dioxin (2,3,7,8‐TCDD) was clarified for the first time at atomic and electronic levels, using our novel tight‐binding quantum chemical molecular dynamics method with first‐principles parameterization. The calculation speed of our new method is over 5000 times faster than that of the conventional first‐principles molecular dynamics method. We confirmed that the structure, energy, and electronic states of the above molecules calculated by our new method are quantitatively consistent with those by first‐principles calculations. After the confirmation of our methodology, we investigated the decomposition reaction dynamics of the above molecules and the calculated dynamic behaviors indicate that the oxidation of the 2,3,4,4′,5‐PeCB, 3,3′,4,4′,5‐PeCB, and 2,3,7,8‐TCDD proceeds through an epoxide intermediate, which is in good agreement with the previous experimental reports and consistent with our static density functional theory calculations. These results proved that our new tight‐binding quantum chemical molecular dynamics method with first‐principles parameterization is an effective tool to clarify the chemical reaction dynamics at reaction temperatures. © 2004 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

Structural fluctuation and concerted motions in F1‐ATPase: A molecular dynamics study

F₁‐ATPase is an adenosine tri‐phosphate (ATP)‐driven rotary motor enzyme. We investigated the structural fluctuations and concerted motions of subunits in F₁‐ATPase using molecular dynamics (MD) simulations. An MD simulation for the α₃β₃γ complex was carried out for 30 ns. Although large fluctuations of the N‐terminal domain observed in simulations of the isolated β_E subunit were suppressed in the complex simulation, the magnitude of fluctuations in the C‐terminal domain was clearly different among the three β subunits (β_E, β_TP, and β_DP). Despite fairly similar conformations of the β_TP and β_DP subunits, the β_DP subunit exhibits smaller fluctuations in the C‐terminal domain than the β_TP subunit due to their dissimilar interface configurations. Compared with the β_TP subunit, the β_DP subunit stably interacts with both the adjacent α_DP and α_E subunits. This sandwiched configuration in the β_DP subunit leads to strongly correlated motions between the β_DP and adjacent α subunits. The β_DP subunit exhibits an extensive network of highly correlated motions with bound ATP and the γ subunit, as well as with the adjacent α subunits, suggesting that the structural changes occurring in the catalytically active β_DP subunit can effectively induce movements of the γ subunit. © 2010 Wiley Periodicals, Inc. J Comput Chem 2010     

Coarse‐grained molecular dynamics simulation study on spherical and tube‐like vesicles formed by amphiphilic copolymers

Molecular‐level understanding of the vesicular structure and formation process is beneficial for potential vesicle applications, especially in drug delivery. In this article, coarse‐grained molecular dynamics simulation was used to study the self‐assembly behavior of amphiphilic poly(acrylic acid)‐b‐polystyrene copolymers in water at different concentrations and PS/PAA block ratios. It was found that various spherical and tube‐like vesicles formed at PS/PAA 3:3 and 4:2. For spherical vesicles, analysis of vesicular structure indicated that the cavity size was influenced by copolymer concentration and wall thickness by the block ratio. Tube‐like vesicle was formed via the fusion of two spherical vesicles, and a key factor for this morphology is polymer movements between inner and outer layer. This simulation study identifies the key factors governing vesicle formation and structure, and provides a guidance to design and prepare various vesicles for wide applications in drug delivery. © 2017 Wiley Periodicals, Inc. J. Polym. Sci., Part B: Polym. Phys. 2017 , 55, 1220–1226     

Short-time Fourier transform analysis of ab initio molecular dynamics simulation: collision reaction between CN and C4H6

Collision reactions between cyano radical (CN) and dimethylacetylene (C4H6) are thought to occur in the atmosphere of Saturn's moon Titan. However, it is difficult to reproduce reactions occurring in unique environments to study their dynamical processes. In this study, collision reactions between CN and C4H6 were investigated using ab initio molecular dynamics (AIMD) simulations. The simulation results were categorized into three kinds: nonreactive collision, incorporation, and substitution. Short-time Fourier transform analysis of velocity autocorrelation functions obtained by the AIMD simulations, which has been recently developed by our research group, was performed to examine the nonequilibrium condition of the vibrational states. Spectrograms, which correspond to the time evolution of power spectra, clarify the relationship between the three reaction channels and the dynamical changes of the vibrational states.     

GridMAT‐MD: A grid‐based membrane analysis tool for use with molecular dynamics

GridMAT‐MD is a new program developed to aid in the analysis of lipid bilayers from molecular dynamics simulations. It reads a GROMACS coordinate file and generates two types of data: a two‐dimensional contour plot depicting membrane thickness, and a polygon‐based tessellation of the individual lipid headgroups. GridMAT‐MD can also account for proteins or small molecules within the headgroups of the lipids, closely approximating their occupied lateral area. The program requires no installation, is fast, and is freely available. © 2008 Wiley Periodicals, Inc. J Comput Chem, 2009     

Electrode polarization effects on interfacial kinetics of ionic liquid at graphite surface: An extended lagrangian-based constant potential molecular dynamics simulation study

Computational models including electrode polarization can be essential to study electrode/electrolyte interfacial phenomena more realistically. We present here a constant-potential classical molecular dynamics simulation method based on the extended Lagrangian formulation where the fluctuating electrode atomic charges are treated as independent dynamical variables. The method is applied to a graphite/ionic liquid system for the validation and the interfacial kinetics study. While the correct adiabatic dynamics is achieved with a sufficiently small fictitious mass of charge, static properties have been shown to be almost insensitive to the fictitious mass. As for the kinetics study, electrical double layer (EDL) relaxation and ion desorption from the electrode surface are considered. We found that the polarization slows EDL relaxation greatly whereas it has little impact on the ion desorption kinetics. The findings suggest that the polarization is essential to estimate the kinetics in nonequilibrium processes, not in equilibrium. © 2019 Wiley Periodicals, Inc.     

Efficient 3D kinetic monte carlo method for modeling of molecular structure and dynamics

Self‐assembly of molecular systems is an important and general problem that intertwines physics, chemistry, biology, and material sciences. Through understanding of the physical principles of self‐organization, it often becomes feasible to control the process and to obtain complex structures with tailored properties, for example, bacteria colonies of cells or nanodevices with desired properties. Theoretical studies and simulations provide an important tool for unraveling the principles of self‐organization and, therefore, have recently gained an increasing interest. The present article features an extension of a popular code MBN Explorer (MesoBioNano Explorer) aiming to provide a universal approach to study self‐assembly phenomena in biology and nanoscience. In particular, this extension involves a highly parallelized module of MBN Explorer that allows simulating stochastic processes using the kinetic Monte Carlo approach in a three‐dimensional space. We describe the computational side of the developed code, discuss its efficiency, and apply it for studying an exemplary system. © 2014 Wiley Periodicals, Inc.     

兔子朊蛋白结构稳定性的分子动力学研究(英文)

朊蛋白病是一种能够对人类和动物带来致命影响,并具有高度传染性的神经退行性疾病.兔子是目前已经报道的哺乳类动物中对朊蛋白病免疫的少数几个物种之一.我们将分子动力学和操控式分子动力学模拟相结合,研究了兔子正常朊蛋白的结构稳定性;同时讨论了蛋白结构的收敛性及刚性分布,并揭示了兔子朊蛋白中关键二级结构的动力学以及受力各向异性特征,证实了兔子朊蛋白结构的稳定性特征.     

基于分子模拟技术煤焦分子模型构建

煤、焦是过程工业的重要原料。因此,有必要深入了解煤、焦分子结构以揭示其反应性。采用Materials Studio 7.0软件,从分子层次研究煤、焦的分子结构。根据已报道的文献,构建煤、焦的初始结构;基于分子力学原理对这些结构进行优化,使得模型物性与煤、焦物性相符;基于退火模拟算法对模型进行优化,从而使得模型密度、元素分析数据与真实值吻合;基于能量最小化原理,对煤、焦模型再次优化,从而获得其稳定、真实的分子构型。由计算结果发现,模型的估算密度、元素组成与已报道一致,说明构建的模型是有效、合理的;在模型优化过程中,相对于其他能量而言,库伦能和范德华能起着重要的作用。因此可以推断在煤、焦热加工过程中,弱键占据主要地位。另外,本文采用分子模拟技术构建煤、焦模型的方法对于构建其他复杂大分子结构有着重要的借鉴作用。     

Wavelet transform analysis of ab initio molecular dynamics simulation: application to core-excitation dynamics of BF3

We propose a novel analysis method of ab initio molecular dynamics (AIMD) simulation using a continuous wavelet transform (c-WT) technique. The c-WT technique, one of the time-frequency signal analysis methods, provides a clear view of the dynamical information in time developments. Combined with the auto-correlation function of velocity by AIMD simulation, c-WT analysis enables us to well understand dynamical distribution, such as the vibrational properties following a change of electronic structure in a molecular system. As a practical application, AIMD simulation of core-excited BF(3) (B1s --> 2a(2) (')) is illustrated. AIMD simulation leads to the change of vibrational motion as well as structural deformation by core-excitation. The c-WT analysis clarifies the relationship between structural deformation and the related significant vibrational modes in core-excitation within 50 fs.     

Hybrid particle‐field molecular dynamics simulations: Parallelization and benchmarks

The parallel implementation of a recently developed hybrid scheme for molecular dynamics (MD) simulations (Milano and Kawakatsu, J Chem Phys 2009, 130, 214106) where self‐consistent field theory (SCF) and particle models are combined is described. Because of the peculiar formulation of the hybrid method, considering single particles interacting with density fields, the most computationally expensive part of the hybrid particle‐field MD simulation can be efficiently parallelized using a straightforward particle decomposition algorithm. Benchmarks of simulations, including comparisons of serial MD and MD‐SCF program profiles, serial MD‐SCF and parallel MD‐SCF program profiles, and parallel benchmarks compared with efficient MD program GROMACS 4.5.4 are tested and reported. The results of benchmarks indicate that the proposed parallelization scheme is very efficient and opens the way to molecular simulations of large scale systems with reasonable computational costs. © 2012 Wiley Periodicals, Inc.