Structure and dynamics of surfactant and hydrocarbon aggregates on graphite: a molecular dynamics simulation study

We have examined the structure and dynamics of sodium dodecyl sulfate (SDS) and dodecane (C12) molecular aggregates at varying surface coverages on the basal plane of graphite via classical molecular dynamics simulations. Our results suggest that graphite-hydrocarbon chain interactions favor specific molecular orientations at the single-molecule level via alignment of the tail along the crystallographic directions. This orientational bias is reduced greatly upon increasing the surface coverage for both molecules due to intermolecular interactions, leading to very weak bias at intermediate surface coverages. Interestingly, for complete monolayers, we find a re-emergent orientational bias. Furthermore, by comparing the SDS behavior with C12, we demonstrate that the charged head group plays a key role in the aggregate structures: SDS molecules display a tendency to form linear file-like aggregates while C12 forms tightly bound planar ones. The observed orientational bias for SDS molecules is in agreement with experimental observations of hemimicelle orientation and provides support for the belief that an initial oriented layer governs the orientation of hemimicellar aggregates.     

Adsorption of a single polymer chain on a surface: A molecular dynamics simulation study

We performed simulations of the physical adsorption of a single globular chain on a surface of hemispherical shape by means of molecular dynamics simulations. For the chain, we took advantage of a united atom model. Interactions within the chain were limited to stretching, bending, and torsional as well as nonbonded interactions between the nonadjacent atoms. The interaction between each chain element and the surface formation are reigned by a Lennard–Jones potential. In this article, we focused on differences in the behavior of the adsorbed globule to the free unadsorbed one particularly in two different zones of the immediate vicinity of the surface. There were strong indications for a localized acceleration of the dynamics as compared with the bulk that appears in an increase of trans–gauche switches. For explanation we came up with an adsorption scenario. Special attention was given to the shift of the percentage of trans and gauche conformations within the globule in dependence on the strength of the adsorption potential that might be related to crystallization or glass transition. © 2001 John Wiley & Sons, Inc. J Polym Sci Part B: Polym Phys 39: 2333–2339, 2001     

Statics and dynamics of ethane molecules in AlPO(4)-5: a molecular dynamics simulation study.

From an experimental perspective, there has been disagreement among researchers on whether ethane would display single-file or normal diffusive behavior in the channels of AlPO(4)-5. Pulsed field gradient nuclear magnetic resonance measurements implied single-file diffusion, while quasielastic neutron scattering showed normal diffusion. In this paper we present the results of extensive classical molecular dynamics simulations of the diffusion of ethane molecules adsorbed in AlPO(4)-5. Our aim is to provide microscopic details of the static and dynamic properties of the adsorbed molecules in order to verify whether the conditions for the single-file regime can be achieved in a nondefective AlPO(4)-5 crystal structure.     

Structure,thermodynamics and diffusion in asymmetric binary mixtures: a molecular dynamics simulation study

Structural and thermodynamic properties as well as diffusion coefficients of binary fluid mixtures with asymmetry in mass, size, charge and their combinations have been studied using classical molecular dynamics simulations. The fluid mixture is modelled as spherical particles interacting via the Weeks–Chandler–Andersen and Coulomb potential. The diameter, charge and mass of the fluid particles are in the range 6–60 Å, 1–10e and 1—500 amu, respectively. Systematic variations in pair-correlation functions, thermodynamic properties as well as the self-diffusion coefficient are found with the size, charge and mass ratio of the particles. The self-diffusion coefficient for systems having more than one type of asymmetry is calculated and expressed in terms of diffusion coefficients of systems with only one type of asymmetry.     

Thermophoresis in liquids: a molecular dynamics simulation study

Thermophoresis in liquids is studied by molecular dynamics simulation (MD). A theory is developed that divides the problem in the way consistent with the characteristic scales. MD is then conducted to obtain the solution of each problem, which is to be all combined for macroscopic predictions. It is shown that when the temperature gradient is applied to the nonconducting liquid bath that contains neutral particles, there occurs a pressure gradient tangential to the particle surface at the particle-liquid interface. This may induce the flow in the interfacial region and eventually the particle to move. This applies to the material system that interacts through van der Waals forces and may be a general source of the thermophoresis phenomenon in liquids. The particle velocity is linearly proportional to the temperature gradient. And, in a large part of the given temperature range, the particle motion is in the direction toward the cold end and decreases with respect to the temperature. It is also shown that the particle velocity decreases or even reverses its sign in the lowest limit of the temperature range or with a particle of relatively weak molecular interactions with the liquid. The characteristics of the phenomenon are analyzed in molecular details.     

Nascent crystallization of a growing chain on a catalyst surface: a nonequilibrium molecular dynamics simulation study

The growing chain molecular dynamics (GCMD) simulation method, a new nonequilibrium molecular dynamics code, is proposed to simulate the polymer chain aggregation behavior during polymerization on a catalyst surface. We found that the growing chain crystallizes on the surface in two stages: the nucleation stage and the crystal growth stage. In the first part of the nucleation period, the short polymerizing chain first absorbs on the surface and can be in either an ordered or disordered structure. Still in the nucleation period, when the chain reaches a degree of polymerization, about 100 bonds, the chain folds into a stable nucleus on the substrate with 3-5 stems. In the crystal growth stage where the polymerization also proceeds, we observed a stem elongation process in combination with a chain folding process. In the stem elongation step, the number of stems in the nucleus remains constant, and all the stems expand together to a length of ca. 5-25 ns. In the subsequent chain folding step, the stem length decreases about 20 bonds within a period of ca. 0.1-0.5 ns. During chain growth, the elongation process and the folding process occur in an alternating and repeated fashion. The crystallization mechanism of the polymerizing chain was discussed.     

Structure and diffusion in amorphous aluminum silicate: a molecular dynamics computer simulation

The amorphous aluminum silicate (Al2O3)2(SiO2) [AS2] is investigated by means of large scale molecular dynamics computer simulations. We consider fully equilibrated melts in the temperature range 6100 K> or =T> or =2300 K as well as glass configurations that were obtained from cooling runs from T=2300 to 300 K with a cooling rate of about 10(12) K/s. Already at temperatures as high as 4000 K, most of the Al and Si atoms are fourfold coordinated by oxygen atoms. Thus, the structure of AS2 is that of a disordered tetrahedral network. The packing of AlO4 tetrahedra is very different from that of SiO4 tetrahedra in that Al is involved with a relatively high probability in small-membered rings and in triclusters in which an O atom is surrounded by four cations. We find as typical configurations two-membered rings with two Al atoms in which the shared O atoms form a tricluster. On larger length scales, the system shows a microphase separation in which the Al-rich network structure percolates through the SiO2 network. The latter structure gives rise to a prepeak in the static structure factor at a wave number q=0.5 A(-1). A comparison of experimental x-ray data with the results from the simulation shows good agreement for the structure function. The diffusion dynamics in AS2 is found to be much faster than in SiO2. We show that the self-diffusion constants for O and Al are very similar and that they are by a factor of 2-3 larger than the one for Si.     

Ion condensation behavior and dynamics of water molecules surrounding the sodium poly(methacrylic acid) chain in water: a molecular dynamics study

All-atom molecular dynamics simulations are used to study the condensation behavior of monovalent (Na(+)) and multivalent (Ca(2+)) salt counterions associated with the co-ions (Cl(-)) surrounding the charged poly(methacrylic acid) (PMAA) chain in water. The study is extended to the influences on chain conformation, local arrangement, and dynamics of water in the highly diluted aqueous solutions. We find that even when the salt ions are monovalent, they attract more than one charged monomer and act as a bridging agent within the chain, as the multivalent salt ions. In principle, the salt ions bridge between not only the "non-adjacent" but also the "adjacent" charged monomers, leading to a more coil-like and a locally stretched conformation, respectively. With an increase in the salt concentration, the amount of coiled-type condensed ions increase and reach a maximum when the chain conformation becomes the most collapsed; whereas, the stretched-type shows an opposite trend. Our results show that the attractive interactions through the condensed salt ions between the non-adjacent monomers are responsible for the conformational collapse. When the salt concentration increases high enough, a significant increase for the stretched-type condensed ions makes an expansion effect on the chain. These stretched-type salt ions, followed by the adsorption of the co-ions and water molecules, tend to form a multilayer organization outside surrounding the PMAA chain. Thus, the expansion degree of the chain conformation is greatly limited. When only the monovalent Na(+) ions are present in the solutions, water molecules are primarily adsorbed into either the condensed Na(+) ions or the COO(-) groups. These adsorbed water molecules form hydrogen bonds with each other and enhance the local bridging behavior associated with the Na(+) condensation on the resultant chain conformation. With an increase in the amount of multivalent Ca(2+) salt ions, more water molecules are bonded directly with the condensed Ca(2+) ions. In this case, only the condensed Ca(2+) ions provide a strong bridging effect within the polymer chain. We observe a significant shift towards a higher frequency of the oxygen vibration spectrum and only a slight shift towards a higher frequency of the hydrogen spectrum for the water molecules associated with the ion condensation.     

Structure and dynamics of water at a clay surface from molecular dynamics simulation

We report a molecular dynamics study of the structure and dynamics of water at a clay surface. The negative charge of the surface and the presence of surface oxygen atoms perturbs water over two to three molecular layers, while the nature of the counterions (Na(+)or Cs(+)) has only a small effect. In the first molecular layer, approximately half of the water molecules are H-bonded to the surface. We also analyze the H-bond network between surface water molecules. The diffusion of water molecules along the surface is slowed down compared to the bulk case. As far as the orientational order and dynamics of the water dipole are concerned, only the component normal to the clay surface is perturbed. We investigate the surface H-bond formation and dissociation dynamics and their coupling to the release of molecules from the first molecular layer. We introduce a simple kinetic model in the spirit of Luzar and Chandler [Nature, 1996, 379, 55] to allow for a comparison with bulk water dynamics. This model semi-quantitatively reproduces the molecular simulation results and suggests that H-bond formation is faster with the surface than in the bulk, while H-bond dissociation is slower.     

Local chain dynamics of a model polycarbonate near glass transition temperature: A molecular dynamics simulation

Constant pressure constant temperature molecular dynamics method is employed to investigate the atomistic scale dynamics of a model Bisphenol A polycarbonate in the vicinity of its glass transition temperature. First, the glass transition temperature and the thermal expansion coefficients of the polymer are predicted by performing simulations at different temperatures. To explore the significance of different modes of motion, various types of time correlation functions are utilized in analyzing the trajectories. In these nanosecond scale simulations, the motion of the chain segments is found to be highly localized with little reorientation of the vectors representing these segments. Detailed analysis of trajectories and the correlation functions of the backbone dihedrals and side methyl groups indicates that they exhibit numerous conformational transitions. The activation energies of the conformational transitions obtained from the simulation are generally larger than the potential barriers for the rotations of these dihedrals, however, both show the same trend. We also have estimated the phenylene ring flip activation energy as 12.6 kcal/mol and the flip frequency as 0.77 MHz at 300 K. These values fall either fall within the range determined by various NMR spectroscopy experiments or slightly out of the range. The study shows that the conformational transitions between the adjacent dihedrals are strongly correlated. Three basic cooperative modes are identified from the simulation. They are: a positive synchronous rotation of two phenylene rings, a negative synchronous rotation of two phenylene rings, and a carbonate group rotation. Above the glass transition temperature, the large scale cooperative motions become much more significant.     

Adsorption of water molecules inside a Au nanotube: a molecular dynamics study

A molecular dynamics simulation of water molecules through a Au nanotube with a diameter of 20 A at bulk densities 0.8, 1, and 1.2 gcm(3) has been carried out. The water molecules inside a nanoscale tube, unlike those inside a bulk tube, have a confined effect. The interaction energy of the Au nanotube wall has a direct influence on the distribution of water molecules inside the Au tube in that the adsorption of the water molecules creates shell-like formations of water. Moreover, the high number of adsorbed molecules has already achieved saturation at the wall of the Au nanotube at three bulk densities. This work compares the distribution percentage profiles of hydrogen bonds for different regions inside the tube. The structural characteristics of water molecules inside the tube have also been studied. The results reveal that the numbers of hydrogen bonds per water molecule influence the orientational order parameter q. In addition, the phenomenon of a group of molecules bonded inside the tube can be observed as the number of hydrogen bonds increase.     

Glass transition in single poly(ethylene oxide) chain: A molecular dynamics simulation study

Local dynamics of single poly(ethylene oxide) chain in various environments (bulk, film, and isolated systems) has been characterized by the reorientation functions of various backbone bond vectors. Within any observation time, the variations of these reorientation functions with the temperature can be well described by the Kohlrausch?Williams?Watts (KWW) like equation, in which the fitted temperature parameter is identified as the glass transition temperature (T _g). The so‐obtained T _g for that polymer faithfully reveals the effects of the observation time, chain flexibility and vector range on the local dynamics. Furthermore, it is found that the KWW like relation is also applicable to the temperature‐dependence of the fraction of frozen atoms or torsions defined by the trajectory radii of gyration or the conformational transitions. Consequently, different motions lead to different values of T _g for the same system. Despite all, the consistent trend can be yielded, namely, T _g (bulk) > T _g (film) > T _g (isolated), which captures the effects of free surfaces on enhanced dynamics. In addition, dynamics heterogeneity in the systems can be quantitatively revealed. The newly proposed method holds a bright promise to predict the T _g values of complex polymers especially for comparisons. © 2016 Wiley Periodicals, Inc. J. Polym. Sci., Part B: Polym. Phys. 2017 , 55 , 178–188     

Structure of short polymers at interfaces: a combined simulation and theoretical study

The structure of polymers confined between surfaces is studied using computer simulation and a density functional approach. The simple model system considers the polymer molecule as a pearl necklace of freely jointed hard spheres, having attractions among the beads, confined between attractive surfaces. This approach uses the universality of the free-energy functional to obtain the self-consistent field required in the single chain simulation. The second-order direct correlation function for the uniform bulk fluid required as input has been calculated from the reference interaction site model integral equation theory using mean spherical approximation. The theoretical results are shown to compare well with the Monte Carlo simulation results for varying densities, chain lengths, and with different attractive interaction parameters. The simulation results on the conformational properties give important indications regarding the behavior of chains as they approach the surfaces.     

Characterizing hydrophobicity at the nanoscale: a molecular dynamics simulation study

We use molecular dynamics (MD) simulations of water near nanoscopic surfaces to characterize hydrophobic solute-water interfaces. By using nanoscopic paraffin like plates as model solutes, MD simulations in isothermal-isobaric ensemble have been employed to identify characteristic features of such an interface. Enhanced water correlation, density fluctuations, and position dependent compressibility apart from surface specific hydrogen bond distribution and molecular orientations have been identified as characteristic features of such interfaces. Tetrahedral order parameter that quantifies the degree of tetrahedrality in the water structure and an orientational order parameter, which quantifies the orientational preferences of the second solvation shell water around a central water molecule, have also been calculated as a function of distance from the plate surface. In the vicinity of the surface these two order parameters too show considerable sensitivity to the surface hydrophobicity. The potential of mean force (PMF) between water and the surface as a function of the distance from the surface has also been analyzed in terms of direct interaction and induced contribution, which shows unusual effect of plate hydrophobicity on the solvent induced PMF. In order to investigate hydrophobic nature of these plates, we have also investigated interplate dewetting when two such plates are immersed in water.     

Effects of solute-solvent proton exchange on polypeptide chain dynamics: a constant-pH molecular dynamics study

A method for performing implicit-solvent molecular dynamics simulations at constant pH was applied to a pentapeptide acetyl-Ala-Asp-Ala-Lys-Ala-amide at pH 4. As a reference, molecular dynamics simulations were done for the same peptide with two variants of its fixed protonation patterns expected to dominate at pH 4, i.e., with a protonated and a deprotonated side chain of the Asp residue and the protonated Lys residue in both cases. The dynamic trajectories of the peptide were used to discuss the problem of the significance of the solute-solvent proton exchange phenomena for the dynamics and structural distributions of the polypeptide chain. The Asp-Lys distance was used as a probe of the overall molecular structure of the investigated pentapeptide. To characterize the dynamics, distributions of the "waiting" times for a transition from a "short" distance conformation to a "long" distance conformation were constructed, based on the generated molecular dynamics trajectories. We show that the relaxation time for the transitions, derived from the constant-pH simulations, is very close to the relaxation time characterizing a permanently protonated molecule, although the average protonation probability of the short-distance conformation is close to zero. However, the distribution of the Asp-Lys distances obtained from constant-pH simulations cannot be reproduced as a linear combination of the distributions resulting from the simulations with fixed protonation states.     

Structure and dynamics of methanol in water: a quantum mechanical charge field molecular dynamics study

An ab initio quantum mechanical charge field molecular dynamics simulation was carried out for one methanol molecule in water to analyze the structure and dynamics of hydrophobic and hydrophilic groups. It is found that water molecules around the methyl group form a cage-like structure whereas the hydroxyl group acts as both hydrogen bond donor and acceptor, thus forming several hydrogen bonds with water molecules. The dynamic analyses correlate well with the structural data, evaluated by means of radial distribution functions, angular distribution functions, and coordination number distributions. The overall ligand mean residence time, τ identifies the methanol molecule as structure maker. The relative dynamics data of hydrogen bonds between hydroxyl of methanol and water molecules prove the existence of both strong and weak hydrogen bonds. The results obtained from the simulation are in excellent agreement with the experimental results for dilute solution of CH(3)OH in water. The overall hydration shell of methanol consists in average of 18 water molecules out of which three are hydrogen bonded.     

A mimetic porous carbon model by quench molecular dynamics simulation

A mimetic porous carbon model is generated using quench molecular dynamics simulations that reproduces experimental radial distribution functions of activated carbon. The resulting structure is composed of curved and defected graphene sheets. The curvature is induced by nonhexagonal rings. The quench conditions are systematically varied and the final porous structure is scrutinized in terms of its pore size distribution, pore connectivity, and fractal dimension. It is found that the initial carbon density affects the fractal dimension but only causes a minor shift in the pore size distribution. On the other hand, the quench rate affects the pore size distribution but only causes a minor shift in the fractal dimension.     

Structure and ultrafast dynamics of liquid water: a quantum mechanics/molecular mechanics molecular dynamics simulations study

A quantum mechanics/molecular mechanics molecular dynamics simulation was performed for liquid water to investigate structural and dynamical properties of this peculiar liquid. The most important region containing a central reference molecule and all nearest surrounding molecules (first coordination shell) was treated by Hartree-Fock (HF), post-Hartree-Fock [second-order Moller-Plesset perturbation theory (MP2)], and hybrid density functional B3LYP [Becke's three parameter functional (B3) with the correlation functional of Lee, Yang, and Parr (LYP)] methods. In addition, another HF-level simulation (2HF) included the full second coordination shell. Site to site interactions between oxygen-oxygen, oxygen-hydrogen, and hydrogen-hydrogen atoms of all ab initio methods were compared to experimental data. The absence of a second peak and the appearance of a shoulder instead in the gO-O graph obtained from the 2HF simulation is notable, as this feature has been observed so far only for pressurized or heated water. Dynamical data show that the 2HF procedure compensates some of the deficiency of the HF one-shell simulation, reducing the difference between correlated (MP2) and HF results. B3LYP apparently leads to too rigid structures and thus to an artificial slow down of the dynamics.     

Anionic clays containing anti-inflammatory drug molecules: comparison of molecular dynamics simulation and measurements

Three representative nonsteroidal anti-inflammatory drug molecules, Ibuprofen, Diclofenac, and Indomethacin, have been intercalated within the galleries of an anionic clay, Mg-Al layered double hydroxide (LDH). X-ray diffraction, IR and Raman vibrational spectroscopy and (13)C cross-polarization magic-angle spinning NMR have been used to characterize the confined drug molecules, while molecular dynamics (MD) simulations were used to probe the interlayer structure, arrangement, orientation, and geometry of the intercalated species. All three drug molecules are arranged as bilayers in the interlamellar space of the anionic clay. But while the structure of the intercalated Ibuprofen is identical to that of the molecule outside the layers, spectroscopy as well as MD simulation shows that there is a change in the geometry of Diclofenac and Indomethacin upon confinement within the galleries of the LDH. The change in geometry of Diclofenac and Indomethacin upon intercalation is shown to originate from the electrostatic interaction between the electronegative chlorine atoms on the drug molecule and the positively charged metal hydroxide sheets of the anionic clay. It is shown that these changes in the geometry of the intercalated drug molecules allow for the observed interlayer spacing to be realized without the bilayers having to interdigitate, which would otherwise have been necessary if the structure of the drug molecules had remained identical to that outside the layers. Comparisons of experimental measurements with simulation have provided a more detailed understanding of the geometry and organization of flexible drug molecules confined in the anionic clay.     

The properties of a single polymer chain in solvent confined in a slit: A molecular dynamics simulation

A single polymer chain in solvent confined in a slit formed by two parallel plates is studied by using molecular dynamics simulation method. The square radii of gyration and diffusion behaviors of polymers are greatly affected by the distance between the two plates, but they do not follow the same way. The chain size decays drastically with increasing h (h is the distance between two plates), until a basin occurs, and a universal h/〈R _g〉₀ dependence for polymer chains with different degrees of polymerization can be obtained. While, for the chain’s diffusion coefficient, it decays monotonously and there is no such basin-like behavior. Furthermore, we studied the radial distribution function of confined polymer chains to explain the reason why there is a difference for the decay behaviors between dynamic properties and static properties. Besides, we also give the degree of confinement dependence of the static scaling exponent for a single polymer chain. Our work provides an efficient way to estimate the dynamics and static properties of confined polymer chains, and also helps us to understand the behavior of polymer chains under confinement.