Equilibration strategies in molecular dynamics simulation of polyethers

A comparison of the NVT and NVE ensemble simulations of tetraglyme showed that, in terms of energy, temperature and most of the structural features the results were very similar. However, major differences were observed in dynamic properties, ie in the mean square displacement and in the O--O distances. A fast equilibration method suitable for amorphous polymer systems is also detailed. This was accomplished by the reassignment of the velocity distribution to the particles, after the interruption of an NVE simulation.

Study of growth mechanism of TiC cluster in ferrite via molecular dynamics simulation

Abstract Atomistic simulations were used to investigate the evolution process of titanium carbide clusters to mature precipitates in ferrite. The typical kinetic of carbide cluster growth was studied in detail through analyzing the atomic interactions of a carbide cluster with scattered carbon atoms. The driving force required for cluster growth was calculated along with the atomic diffusivity in the iron matrix, exploring the change in response as two main growth steps. The growth kinetic improved the understanding of precipitate evolution at the atomic level.

A molecular dynamics simulation of alloy carbide clusters formation

 The formation of alloy carbide cluster in ferrite was investigated via molecular dynamics simulation, which disclosed the cluster property and formation mechanism. These together provided a better fundamental understanding of the cluster formation and firm information for the evolution of cluster and precipitate in high-strength low-alloy steel.

Molecular dynamics study of response of liquid N,N-dimethylformamide to externally applied electric field using a polarizable force field

The behavior of Liquid N,N-dimethylformamide subjected to a wide range of externally applied electric fields (from 0.001 V/nm to 1 V/nm) has been investigated through molecular dynamics simulation. To approach the objective the AMOEBA polarizable force field was extended to include the interaction of the external electric field with atomic partial charges and the contribution to the atomic polarization. The simulation results were evaluated with quantum mechanical calculations. The results from the present force field for the liquid at normal conditions were compared with the experimental and molecular dynamics results with non-polarizable and other polarizable force fields. The uniform external electric fields of higher than 0.01 V/nm have a significant effect on the structure of the liquid, which exhibits a variation in numerous properties, including molecular polarization, local cluster structure, rotation, alignment, energetics, and bulk thermodynamic and structural properties.

Molecular dynamics modelling of diffusional formation of titanium carbide clusters in iron matrix

Molecular dynamics simulation was employed to study the atomic interactions in titanium carbides and iron matrix containing carbon and titanium, which are significant for understanding the formation of titanium carbide cluster during precipitate process. The atoms trajectory and diffusion coefficients of carbon in titanium carbide were analyzed to provide a vacancy-exchanging mechanism and clarify the carbon concentration dependence of carbon diffusion in titanium carbide. The dependence of the formation of titanium carbide cluster in iron matrix on carbon was determined from the study of atoms diffusivity, cluster formation and formation energy of titanium carbide cluster. The simulation results provided insight into the carbon diffusion process and improved the understanding of the formation of titanium carbide cluster.

Simulation of oleuropein structural conformation in vacuum, water and triolein-water systems using molecular dynamics

Oleuropein, the main phenolic compound of olive leaves, exhibits a unique blend of biological activities and has been shown to locate itself at the oil-water (O/W) interface. This behavior could influence the physico-chemical properties of dispersed systems such as emulsions. In this work, we study the effect of the microenvironment (vacuum, water, and triolein-water) on the conformational preferences of oleuropein using molecular dynamics (MD) simulations at 300K for at least 30ns. The seven torsions that describe the flexible skeleton of oleuropein were monitored together with the distance between the glucose (Glu) and hydroxytyrosol (Hyd) moieties (dglu-hyd) of the molecule. The obtained trajectories demonstrated that oleuropein adopts different conformations that depend on the environment. The preferential conformers in each system were analyzed for their molecular geometry and internal energy. In vacuum, the oleuropein preferential conformation is tight with the glucose moiety in close proximity with the hydroxytyrosol moiety. In water, oleuropein preferential conformers presented large differences in their structural properties, varying from a close like U form, and a semi-opened form, to an opened form characterized by high fluctuations in dglu-hyd values. In a triolein-water system, oleuropein tends to adopt a more open form where the glucose moiety could be approximately aligned with the hydroxytyrosol and elenolic acid moieties. Based on a calculation at the HF/6-31G* level, these flexibilities of oleuropein required energy of 19.14kcal/mol in order to adopt the conformation between water and triolein-water system. A radial distribution function (RDF) analysis showed that specific hydroxyl groups of Hyd and Glu interact with water molecules, enabling us to understand the amphiphilic character of oleuropein at the triolein-water interface. MD calculations together with interfacial tension measurements revealed that the oleuropein binding at O/W interface is an enthalpy driven mechanism.

Ion clustering in molecular dynamics simulations of sodium iodide solutions

Model systems of sodium iodide dissolved in dimethyl ether or 1,2-dimethoxyethane (glyme) were studied in order to investigate the structural and dynamic properties of ionic solutions in small and polymeric ethers. Full molecular dynamics simulations were performed at a range of different salt concentrations. An algorithm was designed which assigns ions to clusters and then calculates all the terms which contribute to ionic conductivity. In dilute solutions, free ions are the most common ionic species, followed by ion pairs. As the concentration increases, pairs become the most common species, with significant concentrations of clusters with 3 through 6 ions. Changing the solvent from dimethyl ether to glyme significantly decreases the ion clustering due to the chelate effect in which the two oxygens on a solvent stabilize an associated cation. The conductivity in stable systems is shown to be primarily the result of the movement of free ions and the relative movement of ions within neutral pairs. The Nernst-Einstein relation, commonly used in the discussion of polymer electrolytes, is shown to be inadequate to quantitatively describe conductivity in the model systems.

Characterization of structure and energy of TixCy cluster at early formation stage in iron matrix by molecular dynamics

The structure, energy and bonding property of TixCy clusters formed in iron matrix were studied through molecular dynamics (MD) simulation method. The initial clusters with 1D-linear, 2D-ring, and 3D-tetrahedral structures were determined and their stability was calculated. The effect of temperature on the stability of the clusters was also discussed. In addition, the dissociation path of TiC clusters in iron matrix and the corresponding energy variation were analyzed. © 2014 Elsevier B.V. All rights reserved.

Molecular dynamics approach to the structural characterization and transport properties of poly(acrylonitrile)/N,N-dimethylformamide solutions

Poly(acrylonitrile) (PAN) in N,N-dimethylformamide (DMF) is a popular solution for producing large variety of polymer products. To precisely describe the behaviours of PAN and DMF in the synthesis processes, it is significant to call for more details about the structure, some thermodynamic and dynamical properties of PAN-DMF solutions. A PAN-DMF solution was simulated via molecular dynamics with an all-atom OPLS type potential in both the NPT and NVT ensembles. The simulation results were evaluated with quantum mechanical calculations (MP2/6-311 ++G(d,p) and counterpoise procedure) and were compared with available experimental results. The liquid structure was illustrated with pair correlation functions and transport and dynamics properties were calculated with the mean-square displacements MSD and the velocity autocorrelation functions. The strong H-bonds of C≡N « H-C=O, CH » O=C-H and CH2 O=C-H, with distances of 2.55 Å, 2.55 Å and 2.65 Å, respectively, were found. The largest interaction energy of - 7.157 kcal/mol between DMF molecules and PAN molecules was found at 4.9 Å center-of-mass distance. A potential profile of intermolecular interaction of DMF with PAN along the interaction distance was presented, clearly showing an increase of DMF vaporisation heat when it getting close to PAN molecules. This provided very useful information to analyse the vaporisation behaviours of DMF at the microscopic level, which is essential to comprehensively understand molecular rearrangements towards the design of synthetic processes. The impact of the presence of the PAN on the DMF solution properties were also benchmarked with pure DMF solution.

Parallelization of dissipative particle dynamics simulation

The dissipative particle dynamics simulation is usually used to study polymer in mesoscopic space. The traditional methods are resource intensive, especially when the scale of research is large. Therefore, improving computing efficiency is a key point in this research area. Two major issues are addressed in this paper. First, the DPD methods are analysed and the most time-consuming parts are identified: conservative force, dissipative force and random force. Second, we describe how to parallelize the existing serial application in the Windows Compute Cluster Server (WCCS) platform. The results show that the parallel algorithm not only effectively shortens the computing time, but also improves the resource utilization rate.

Molecular dynamics simulations of highly concentrated salt solutions : structural and transport effects in polymer electrolytes

Structural, thermodynamic and transport properties have been calculated in concentrated non-aqueous NaI solutions using molecular dynamics simulations. Although the solvent has been represented by a simplistic Stockmayer fluid (spherical particles with point dipoles), the general trends observed are still a useful indication of the behavior of real non-aqueous electrolyte systems. Results indicate that in low dielectric media, significant ion pairing and clustering occurs. Contact ion pairs become more prominent at higher temperatures, independent of the dielectric strength of the solvent. Thermodynamic analysis shows that this temperature behavior is predominantly entropically driven. Calculation of ionic diffusivities and conductivities in the NaI/ether system confirms the clustered nature of the salt, with the conductivities significantly lower than those predicted from the Nernst-Einstein relation. In systems where the solvent-ion interactions increase relative to ion-ion interactions (lower charge or higher solvent dipole moment), less clustering is observed and the transport properties indicate independent motion of the ions, with higher calculated conductivities. The solvent in this system is the most mobile species, in comparison with the polymer electrolytes where the solvent is practically immobile.

Highly concentrated salt solutions : molecular dynamics simulations of structure and transport

Molecular dynamics (MD) simulations in NaI solutions, where the solvent has been represented by the Stockmayer fluid, were performed as a function of temperature, salt concentration, and solvent dipole strength. At higher temperatures contact ion pairs become more prevalent, regardless of solvent strength. An examination of the temperature dependence of the potential of mean force demonstrates the entropic nature of this effect. The transport properties calculated in the simulations are dependent on the balance between solvent dielectric constant and ion charge. In systems with a large solvent dipole moment, the ions appear to be independently mobile, and deviations from Nernst–Einstein behavior are small. In systems of smaller solvent dipole moment or greater ion charge, the ions form clusters, and large deviations from Nernst–Einstein behavior are observed.