Molecular dynamics simulation of thermomechanical properties of montmorillonite crystal. 3. montmorillonite crystals with PEO oligomer intercalates

We present the results of molecular dynamics (MD) simulation of the structure and thermomechanical behavior of Wyoming-type Na+-montmorillonite (MMT) with poly(ethylene oxide) (PEO) oligomer intercalates. Periodic boundary conditions in all three directions and simulation cells containing two MMT lamellae [Si248Al8][Al112Mg16]O640[OH]128 oriented parallel to the XY-plane were used. The interlamellar space, or gallery, between neighboring MMT lamellae was populated by 24 Na+ counterions and PEO macromolecules of different lengths, ranging from 2 up to 240 repeat units. We considered three different loadings of PEO within the gallery: 80, 160, and 240 repeat units, corresponding to 13, 23, and 31 wt % PEO based on total mass of the nanocomposite, respectively. In the cases of 13 and 23 wt %, the polymer chains formed one or two well-defined amorphous layers with interlayer distances of 1.35 and 1.8 nm, respectively. We have observed also formation of a wider monolayer gallery with interlayer distances of 1.6 nm. Three-layer PEO films formed in the case of 31 wt % loading. The thermal properties were analyzed over the range 300-400 K, and the isothermal linear compressibility, transversal moduli, and shear moduli were calculated at 300 K. These properties are compared with the results of our simulation of thermal and mechanical properties of MMT crystal with galleries filled by one or two water layers as well as with those of an isolated clay nanoplate.     

Molecular dynamics simulations of nanocomposites based on poly(epsilon-caprolactone) grafted on montmorillonite clay

Intercalated and exfoliated models of polymer nanocomposites based on poly(epsilon-caprolactone) and functionalized montmorillonite clay are studied by means of molecular dynamics simulations. Intercalated and exfoliated models are considered for probing the structural characteristics of the corresponding nanocomposites prepared by melt intercalation and in situ polymerization, respectively. In the exfoliated system, the organization of the polymer chains onto the clay surface is examined in terms of the density profiles and the order parameter function. A layered structure can clearly be seen to form near the surface with density maxima higher than in amorphous poly(epsilon-caprolactone). This can be viewed as an increase in effective particle thickness, which can contribute to the outstanding gas barrier properties of the exfoliated nanocomposites. The comparison of the structures and energetics of the intercalated model with those of a nanocomposite model based on a nonfunctionalized clay indicates nearly similar characteristics. Nevertheless, the slight differences observed for the interfacial polymer density and clay- and surfactant-polymer binding energies can account for the differences in rheological measurements. The results also suggest that the difference in morphology obtained for the nanocomposites prepared by the two synthetic approaches can be ascribed to both a difference in interfacial polymer density and the formation of bridging polymer chain structures that hinder the exfoliation process.     

Effects of montmorillonite modification on optical properties of heterogeneous nematic liquid crystal - clay mineral nanocomposites

Investigations of montmorillonite (MMT) clay mineral modification effects on electro-optical properties of nanocomposites, based on the nematic liquid crystal 4-pentyl-4'-cyanobiphenyl (5CB) and MMT have been carried out. Only the composite with MMT modified by an organic surface-active substance, dioctadecyldimethylammonium chloride, has been shown to manifest electro-optical memory effect and contrast. A polar dopant (acetone), added to the mixture, significantly increases composite homogeneity. Results of IR spectroscopy measurements lead to the conclusion, that there is a mutual influence of components on each other in organoclay systems, which appears as an alignment of near-surface layers of both the organic and inorganic components of the composite. Due to such interactions these systems show electro-optical contrast and memory effect. A composite with organophobic Na-MMT does not show these electro-optical properties, due to the absence of component interactions, as shown by IR spectroscopic data.     

Molecular dynamics simulation of LiTFSI-acetamide electrolytes: structural properties

The liquid structures of nonaqueous electrolytes composed of lithium bis(trifluoromethylsulfonyl)imide (LiTFSI) and acetamide, with LiTFSI/acetamide molar ratios of 1:2, 1:4, and 1:6, were studied by molecular dynamics simulations. The simulations indicate that the Li+ cations prefer to be six-coordinate by the sulfonyl oxygen atoms of the TFSI- anions and the carbonyl oxygen atoms of the acetamide molecules, rather than by the most electronegative nitrogen atom of the TFSI- anion. Therefore, close Li+-TFSI- contact pairs exist in the system. The TFSI- anion prefers to provide only one of four possible oxygen atoms to coordinate to the same Li+ cation. Three conformations (cis, trans, and gauche) of the TFSI- anions were found to coexist in the liquid electrolyte. At high salt concentrations, the TFSI- anions mainly adopt the gauche conformation in order to provide more oxygen atoms to coordinate to different Li+ cations, while simultaneously reducing the repulsion among the Li+ cations. On the other hand, the fraction of TFSI- anions adopting the cis conformation is largest for the system with the molar ratio of 1:6, in which many clusters, mainly composed of the Li+ cations and the TFSI- anions, are immersed in the acetamide molecules. The size and charge distribution of clusters were also investigated. In the system with the molar ratio of 1:2, nearly all of the ions in the PBC (periodic boundary conditions) box aggregate into a bulky cluster that gradually disassembles into small clusters with decreasing salt concentration. The addition of acetamide molecules was found to effectively relax the liquid electrolyte structure, and the system with the molar ratio of 1:4 was found to exhibit a more homogeneous liquid structure than the other two electrolyte systems with molar ratios of 1:2 and 1:6.     

Molecular dynamics simulation of uranyl(VI) adsorption equilibria onto an external montmorillonite surface

We used molecular dynamics simulations to study the adsorption of aqueous uranyl species (UO(2)(2+)) onto clay mineral surfaces in the presence of sodium counterions and carbonato ligands. The large system size (10,000 atoms) and long simulation times (10 ns) allowed us to investigate the thermodynamics of ion adsorption, and the atomistic detail provided clues for the observed adsorption behavior. The model system consisted of the basal surface of a low-charge Na-montmorillonite clay in contact with aqueous uranyl carbonate solutions with concentrations of 0.027 M, 0.081 M, and 0.162 M. Periodic boundary conditions were used in the simulations to better represent an aqueous solution interacting with an external clay surface. Uranyl adsorption tendency was found to decrease as the aqueous uranyl carbonate concentration was increased, while sodium adsorption remained constant. The observed behavior is explained by physical and chemical effects. As the ionic strength of the aqueous solution was increased, electrostatic factors prevented further uranyl adsorption once the surface charge had been neutralized. Additionally, the formation of aqueous uranyl carbonate complexes, including uranyl carbonato oligomers, contributed to the decreased uranyl adsorption tendency.     

Molecular dynamics simulation of polymer dispersed liquid crystal droplets under competing boundary conditions

This paper is devoted to the molecular dynamics simulation of structural organization inside a polydispersed liquid crystal (LC) droplet under competing boundary conditions. The droplet is assumed to be placed at the liquid crystal interface between two different regions of the solid polymer matrix, which accordingly separates the droplet into two hemispheres: the first of these is under radial boundary conditions; the second hemisphere is under bipolar boundary conditions. The droplet is considered as a jagged sphere filled with LC molecules, modelled as classical spins (unit vectors), whose centres of mass are associated with sites of a cubic lattice inside the cavity. The orienting action of the polymer matrix, and hence the resulting boundary conditions, are modelled by the interaction between the internal LC molecules (possessing only orientational degrees of freedom), and those of a delimiting surface layer (a jagged spherical shell), whose orientations are fixed, radial or bipolar, respectively. All interactions are modelled by the short range McMillan pair potential. The molecular orientation inside the LC droplet has been determined for various anchoring strengths of the interaction between internal spins and boundary layers. We have investigated the structure of the spherical defect resulting in the central region of the droplet, as well as of the boojum ‐ like defects existing near the poles of the droplet. It has been found that a change of relative radial and bipolar anchoring strengths can affect both central and boojum ‐ like defects. The effect of an external field on the molecular orientation inside the droplet has also been investigated. It has been found that a sufficiently strong external field increases the radius of the spherical defect placed in the central region of the droplet.     

Molecular dynamics simulation of polymer dispersed liquid crystal droplets under competing boundary conditions

This paper is devoted to the molecular dynamics simulation of structural organization inside a polydispersed liquid crystal (LC) droplet under competing boundary conditions. The droplet is assumed to be placed at the liquid crystal interface between two different regions of the solid polymer matrix, which accordingly separates the droplet into two hemispheres: the first of these is under radial boundary conditions; the second hemisphere is under bipolar boundary conditions. The droplet is considered as a jagged sphere filled with LC molecules, modelled as classical spins (unit vectors), whose centres of mass are associated with sites of a cubic lattice inside the cavity. The orienting action of the polymer matrix, and hence the resulting boundary conditions, are modelled by the interaction between the internal LC molecules (possessing only orientational degrees of freedom), and those of a delimiting surface layer (a jagged spherical shell), whose orientations are fixed, radial or bipolar, respectively. All interactions are modelled by the short range McMillan pair potential. The molecular orientation inside the LC droplet has been determined for various anchoring strengths of the interaction between internal spins and boundary layers. We have investigated the structure of the spherical defect resulting in the central region of the droplet, as well as of the boojum - like defects existing near the poles of the droplet. It has been found that a change of relative radial and bipolar anchoring strengths can affect both central and boojum - like defects. The effect of an external field on the molecular orientation inside the droplet has also been investigated. It has been found that a sufficiently strong external field increases the radius of the spherical defect placed in the central region of the droplet.     

Molecular simulation to discover rheological properties and soil-binding ability of PHPA polymer on montmorillonite surface

Synthetic polymer fluids are increasingly being applied to support excavations in deep foundations. As these fluids are molecularly engineered, their underlying microstructure interaction with in situ soils significantly affect excavation stability and soil dispersion. However, little molecular-scale research has been done on the rheological behavior of partially hydrolyzed polyacrylamides (PHPA) polymer fluids on the clay surface. Molecular models of the clay–polymer systems are constructed using PHPA on montmorillonite (MMT) clay surface. Initial rheological properties and soil-binding ability at different shear rates, temperatures, and polymer concentrations are first studied using molecular dynamics (MD) simulations. It is found that the functional groups of PHPA can interact with the MMT surface and form a viscous film under the atomic interaction of hydrogen bonds, water bridges, and electrostatic attraction. The shear stress, σ increases with the shear rate and follows the power-law model. And the viscosity, η decreases as the shear rate increases, which is consistent with the experimental trend. However, the σ and η decrease with the increase of temperature. And the action mode of PHPA concentration has been identified from the MD perspective. This work provides insight into the molecular mechanism for PHPA's rheology on the clay surface and their interaction.     

Molecular dynamics simulation of the interaction between methyl benzotriazole and Cu2O crystal

In circulating water system , methyl benzotriazole (TTA) is one of the common corrosion inhibitors for copper. But the inhibition mechanisms have not been clearly understood so far. In different number of water molecules, the interaction between TTA and Cu₂O (copper surface) was investigated with molecular dynamics (MD) method. The results showed that the MD simulation result with water was more consistent with the experiment results. In different number of water molecules, the sequence of the interaction energies between TTA and Cu₂O (001) was ? E ₁ (150H₂O) > ? E ₁(200H₂O) > ? E ₁(100H₂O) > ? E ₁(50H₂O) > ? E ₁(0H₂O). The number of water molecules had an important influence on the interaction between corrosion inhibitors and Cu₂O crystal. From non‐bond energy and pair correlation functions, the interaction energies of the model system were mainly contributed by the non‐bond interaction. Strong adsorption could be raised by the Coulomb interaction between the negatively charged functional groups in TTA and the positive copper ions in the Cu₂O (001) face, and further interaction between aggressive media and copper could be restricted. So, copper corrosion could be avoided. Chemical bonds and non‐bond interactions were formed between TTA and Cu₂O (001) in different number of water molecules. Water molecules could not be ignored during the MD simulation, too. The results obtained here may provide theoretical supports for developing new corrosion inhibitors.     

Molecular dynamics simulation of semirigid polymers

The chain rigidity of poly(p-hydroxybenzoate) was estimated through the theoretical evaluation of its persistence length (L_p). A non-Brownian molecular dynamics (MD) simulation of an isolated chain with 20 monomeric units was performed. The sampled conformational population was analyzed and the orientational correlation function between monomeric units along the chain was calculated. An algorithm based on the worm-like chain model was applied to evaluate the persistence length. The results were compared with those obtained from equilibrium models like the freely-rotating-chain and the rotational-matrix method with fluctuations. Equilibrium models give different results depending on the degree of accuracy used in describing the monomeric unit. The inclusion of thermal fluctuations is crucial to obtain realistic results. These coincide with those given by MD simulation when only nearest-neighbour orientational correlations are taken into account: inclusion of higher-order correlation terms leads to lower values of the persistence length. The origin of this discrepancy was investigated. The MD simulation results are characterized by an overrepresentation of conformations with a short end-to-end distance resulting from an anomalous energy concentration in the first bending mode of the chain. In analogy with previous simulation results from systems characterized by a week coupling amoung their degrees of freedom, failure in the energy equipartition is proposed as a likely explanation of the anomalous dynamical behaviour.     

液态水的分子动力学模拟

用分子动力学(MD)模拟方法在150~376K的温度范围内对液态水的微正则系统进行了研究。考察了液态水的结构及其性质。模拟采用了由从头算得出的柔性水-水相互作用势MCYL。对时间和空间的平均得出了液态中水分子几何构型及温度改变所引起的液态水结构变化。对径向分布函数gOH, gOO, gHH及配位数的分析表明, 在所考察的温度范围内, 每个水分子与相邻分子形成的氢键数为2~3, 水分子在参与的2个氢键中同时作为授受体。结合对振动谱的研究表明在低温时液态水形成的网络结构可能随温度的升高而形成小的簇结构。     

Molecular dynamics simulation of liquid trimethylphosphine

Structural and dynamical properties of liquid trimethylphosphine (TMP), (CH(3))(3)P, as a function of temperature is investigated by molecular dynamics (MD) simulations. The force field used in the MD simulations, which has been proposed from molecular mechanics and quantum chemistry calculations, is able to reproduce the experimental density of liquid TMP at room temperature. Equilibrium structure is investigated by the usual radial distribution function, g(r), and also in the reciprocal space by the static structure factor, S(k). On the basis of center of mass distances, liquid TMP behaves like a simple liquid of almost spherical particles, but orientational correlation due to dipole-dipole interactions is revealed at short-range distances. Single particle and collective dynamics are investigated by several time correlation functions. At high temperatures, diffusion and reorientation occur at the same time range as relaxation of the liquid structure. Decoupling of these dynamic properties starts below ca. 220 K, when rattling dynamics of a given TMP molecules due to the cage effect of neighbouring molecules becomes important.     

Molecular dynamics simulation of surface morphology and thermodynamic properties of polyaniline nanostructured film

To develop predictive models in nanostructured films, there is an ongoing research to validate molecular dynamics (MD) simulation results with experimental data. The morphology and surface topography of polyaniline (PANI) nanostructured film coated on a TiO₂ nanocrystalline surface were investigated by scanning electron microscopy and atomic force microscopy, respectively. The atomistic model of the simulated PANI was generated using energy minimization with a condensed‐phase optimized molecular potential for atomistic studies force field function to reach a thermodynamic equilibrium state. Various parameters of PANI such as density, energy, cavity size, and free volume distributions are calculated. MD simulation has also been used to obtain specific volume (V) as a function of temperature (T). It is demonstrated that this V–T curve can be used to determinate glass transition temperature T_g, reliably. Although experimental data available for the PANI film are very limited, simulation results such as density and T_g are in good agreement with the experimental values reported in the literature. Comparison of the surface topography of PANI demonstrates a reasonable trend between atomic force microscopy image analysis and the MD simulation results at various temperatures. Copyright © 2014 John Wiley & Sons, Ltd.     

Formability and properties of self-standing clay film by montmorillonite with different interlayer cations

Influences of exchangeable interlayer cations were investigated on self-standing film formability, film morphology, and properties of the clay films such as flexibility and gas barrier property. Ion-exchanged montmorillonite samples were prepared by a cation exchange from naturally bearing cation, mostly Na⁺, to Li⁺, Mg²⁺, Ca²⁺, Al³⁺, and Fe^{2+, 3+}. Self-standing films were prepared from aqueous colloidal dispersions of these montmorillonite samples with no additives. The montmorillonite samples with monovalent or divalent cation formed flat self-standing films while the Al-montmorillonite sample produced a distorted film. The Fe-montmorillonite sample formed many separated reddish-brown rod-shaped pieces. Clay film microstructures were different with interlayer cations. The films with monovalent interlayer cations were constructed by the stacking of units with delicately waved thin clay sheets in the whole film, but other films show different morphologies between the upper side and lower side; the upper side is laminated with thin sheets; the lower side is laminated with large thick sheets.The self-standing films’ flexibility and gas barrier property differed according to the interlayer cations. These properties were good in cases of samples with monovalent cations. The innumerable short wave and sheet thinness are considered to foster good flexibility and gas barrier properties. The differences in film formability and properties of the films are attributable to different swellability among samples with different interlayer cations. The montmorillonite samples with monovalent cations swell sufficiently by water, but those with polyvalent cations swell poorly. In the latter case, clay crystals aggregate in water, then the aggregate grows into large particles, creating a film with large particles.     

Molecular dynamics simulation of a lipid diamond cubic phase.

This paper presents the first atomistic simulation of a cubic membrane phase. Using the molecular dynamics simulation technique both the global and the local organization of glycerolmonoolein molecules inside the diamond cubic phase are studied. Multinanosecond simulations reveal that the center of the cubic bilayer remains close to the infinite periodic minimal surface that describes the diamond geometry. We further show that the equilibrium structure of the surfactant molecules inside the cubic phase is very similar to their structure inside a simulated lamellar bilayer. The small differences arise from the packing constraints of the surfactants within the cubic phase which has an area per surfactant that increases toward the bilayer center.     

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.     

Liquid crystal transitions induced by spherical particles in anisotropic polymer brushes: Molecular dynamics simulation

A change in the structure of a polymer brush formed by chains composed of anisotropic monomers brought about by filling the volume of the system with free spherical particles was studied by means of molecular dynamics simulation. Microphase separation takes place at a certain concentration of spherical particles inside the brush. A microphase that contains a small amount of isotropic particles and is characterized by a high density and high ordering of the monomers is formed near the wall. A disordered microphase with a low monomer density and a high concentration of spherical particles is located at the periphery of the brush. The microphases are separated from one another by a distinct boundary. As the concentration of isotropic particles further increases, the ordered microphase expands until the entire volume of the brush becomes occupied. It was also found that the brush in the ordered state forms different LC structures, namely, smectic or nematic entities at relatively low concentrations or at higher concentrations of spheres, respectively.     

Structural and electrical properties of doped polypyrrole and its composite with montmorillonite clay

Polypyrrole is synthesized and doped with Dodecylbenzenesulphonic acid. The doped Polypyrrole with Dodecylbenzenesulphonic acid is intercalated into the layers of Montmorillonite clay successfully by in situ polymerization. The structural properties of synthesized doped Polypyrrole and intercalated doped Polypyrrole were studied by XRD analysis. The crystallinity of intercalated doped Polypyrrole into the layers of Montmorillonite clay is confirmed by means of X-rays diffraction studies, which is more than the doped Polypyrrole. Enhanced d-spacing of Montmorillonite confirmed that doped Polypyrrole is interclated into the layers of Montmorillonite clay at nanoscale. The scanning electron micrographs also confirm the formation of dual phase of platelet as well as of flaky structure of intercalated doped Polypyrrole. Temperature dependant conductivity showed three dimensional variable ranges hopping model. Activation energy, density of states and hopping length are calculated and found to be influenced by intercalating Doped Polypyrole into the layers of Montmorillonite clay.     

Molecular dynamics simulation of the crystal nucleation behavior of a single chain touching a substrate surface

Molecular dynamics simulation has been used in exploring the crystal nucleation behaviour of a single chain touching a substrate surface. It shows that a polyethylene chain (980 CH₂) changed its overall shape from an isotropic coil to an oriented one in the case of touching a substrate surface of amorphous carbons at 300 K. Most repeats of the chain were aligned and ordered in a zigzag package. Surprisingly, the direction of the package is not parallel to the plane of the substrate, but almost perpendicular to it. This is in accordance with experimental observations.