Molecular dynamics simulations of carbon nanotube/silicon interfacial thermal conductance

Using molecular dynamics simulations with Tersoff reactive many-body potential for Si-Si, Si-C, and C-C interactions, we have calculated the thermal conductance at the interfaces between carbon nanotube (CNT) and silicon at different applied pressures. The interfaces are formed by axially compressing and indenting capped or uncapped CNTs against 2 x 1 reconstructed Si surfaces. The results show an increase in the interfacial thermal conductance with applied pressure for interfaces with both capped and uncapped CNTs. At low applied pressure, the thermal conductance at interface with uncapped CNTs is found to be much higher than that at interface with capped CNTs. Our results demonstrate that the contact area or the number of bonds formed between the CNT and Si substrate is key to the interfacial thermal conductance, which can be increased by either applying pressure or by opening the CNT caps that usually form in the synthesis process. The temperature and size dependences of interfacial thermal conductance are also simulated. These findings have important technological implications for the application of vertically aligned CNTs as thermal interface materials.     

Molecular dynamics simulations of monoclinic calcium apatites: A universal equation of state

Molecular dynamics simulations of monoclinic (P2₁/b) hydroxy- and chlorapatite were undertaken in the range 498 K < T < 1298 K, and for pressures up to 7.5 GPa. The all-atom Born–Huggins–Mayer force field, that had been previously used to successfully describe the room temperature isotherms of both compounds, was also used in this work. The isothermal sets of p–V data generated by simulation were each fitted to the three-parameter form of the isothermal Parsafar and Mason equation of state (EoS) with an accuracy better than 0.07%. Taking the temperature dependence of the coefficients into account, it was found that the MD data are satisfactorily reproduced by the universal EoS. The isothermal compressibility coefficient dependence with pressure can be described by a linear relation.     

Molecular dynamics simulations for human CAR inverse agonists

Constitutive androstane receptor (CAR), along with pregnane x receptor (PXR), is an important metabolic sensor in the hepatocytes. Like all other nuclear receptors (NRs), CAR works in concert with coregulator proteins, coactivators, and corepressors which bind to the NRs. The main basis for the receptor to distinguish between coactivators and corepressors is the position of the C-terminal helix 12 (H12), which is determined by the bound NR ligand. CAR, having constitutive activity, can be repressed or further activated by its ligands. Crystal structure of human CAR bound to an agonist and a coactivator peptide is available, but no structural information on an inverse agonist-bound human CAR and a corepressor exists. In our previous molecular dynamics (MD) studies, no corepressor peptide was included. Therefore, probably due to the strong interactions which keep the relatively short H12 of CAR in the active position, the structural changes elicited by inverse agonists were very subtle, and H12 of CAR seemed to more or less retain its active conformation. Here, we have run a series of MD simulations to study the movement of H12 in the presence of both activating and repressing ligands as well as a corepressor peptide. The presence of the corepressor on the coregulator surface of CAR induced a clear shift of H12 of the inverse agonists-bound CAR. In general, H12 moved toward H10 and not away from the ligand binding domain, as seen in some other NRs. However, H12 of CAR is short enough that this movement seems to be adequate to accommodate the binding of the corepressor.     

Molecular dynamics simulations for water and ions in protein crystals

The spatial and temporal properties of water and ions in bionanoporous materials-protein crystals-have been investigated using molecular dynamics simulations. Three protein crystals are considered systematically with different morphologies and chemical topologies: tetragonal lysozyme, orthorhombic lysozyme, and tetragonal thermolysin. It is found that the thermal fluctuations of C(alpha) atoms in the secondary structures of protein molecules are relatively weak due to hydrogen bonding. The solvent-accessible surface area per residue is nearly identical in the three protein crystals; the hydrophobic and hydrophilic residues in each crystal possess approximately the same solvent-accessible surface area. Water distributes heterogeneously and has different local structures within the biological nanopores of the three protein crystals. The mobility of water and ions in the crystals is enhanced as the porosity increases and also by the fluctuations of protein atoms particularly in the two lysozyme crystals. Anisotropic diffusion is found preferentially along the pore axis, as experimentally observed. The anisotropy of the three crystals increases in the order: tetragonal thermolysin < tetragonal lysozyme < orthorhombic lysozyme.     

Molecular dynamics simulations of water droplets on polymer surfaces

Molecular dynamics simulations were used to study the wetting of polymer surfaces with water. Contact angles of water droplets on crystalline and two amorphous polyethylene (PE) and poly(vinyl chloride) (PVC) surfaces were extracted from atomistic simulations. Crystalline surfaces were produced by duplicating the unit cell of an experimental crystal structure, and amorphous surfaces by pressing the bulk polymer step by step at elevated temperature between two repulsive grid surfaces to a target density. Different-sized water droplets on the crystalline PE surface revealed a slightly positive line tension on the order of 10(-12)-10(-11) N, whereas droplets on crystalline PVC did not yield a definite line tension. Microscopic contact angles produced by the simple point charge (SPC) water model were mostly a few degrees smaller than those produced by the extended SPC model, which, as the model with lowest bulk energy, presents an upper boundary for contact angles. The macroscopic contact angle for the SPC model was 94 degrees on crystalline PVC and 113 degrees on crystalline PE. Amorphicity of the surface increased the water contact angle on PE but decreased it on PVC, for both water models. If the simulated contact angles on crystalline and amorphous surfaces are combined in proportion to the crystallinity of the polymer in question, simulated values in relatively good agreement with measured values are obtained.     

Molecular dynamics simulations of structure and dynamics of organic molecular crystals

A set of model compounds covering a range of polarity and flexibility have been simulated using GAFF, CHARMM22, OPLS and MM3 force fields to examine how well classical molecular dynamics simulations can reproduce structural and dynamic aspects of organic molecular crystals. Molecular structure, crystal structure and thermal motion, including molecular reorientations and internal rotations, found from the simulations have been compared between force fields and with experimental data. The MM3 force field does not perform well in condensed phase simulations, while GAFF, CHARMM and OPLS perform very similarly. Generally molecular and crystal structure are reproduced well, with a few exceptions. The atomic displacement parameters (ADPs) are mostly underestimated in the simulations with a relative error of up to 70%. Examples of molecular reorientation and internal rotation, observed in the simulations, include in-plane reorientations of benzene, methyl rotations in alanine, decane, isopropylcyclohexane, pyramidal inversion of nitrogen in amino group and rotation of the whole group around the C-N bond. Frequencies of such dynamic processes were calculated, as well as thermodynamic properties for reorientations in benzene and alanine. We conclude that MD simulations can be used for qualitative analysis, while quantitative results should be taken with caution. It is important to compare the outcomes from simulations with as many experimental quantities as available before using them to study or quantify crystal properties not available from experiment.     

Molecular dynamics simulations of charged dendrimers: low-to-intermediate half-generation PAMAMs

We have performed atomistic molecular dynamics simulations of PAMAM dendrimers of generations 0.5, 1.5, 2.5, 3.5, and 4.5. The simulated systems comprise the charged dendrimer and its counterions embedded in a dielectric continuum (i.e., without explicit solvent). Structural properties of these dendrimers, like the radius of gyration, the principal moments of inertia, and the segment density profiles, were evaluated from the simulations. The average radius of gyration obtained for the intermediate half-generations 2.5, 3.5, and 4.5 follows the same scaling law that was previously inferred from simulations of full-generation PAMAMs, Rg approximately M1/3, and is characteristic of space-filling objects. The low half-generations 0.5 and 1.5 deviate, however, to greater Rg values. The shape of the smaller dendrimers is approximately that of a prolate ellipsoid, which becomes more spherical for higher generations. The segment density profiles show features identical to those obtained in other simulations of flexible-chain dendrimers, like dendron-backfolding. Two slightly different configurations, in terms of size and shape, were identified for generation 2.5. The radial distributions of counterions extracted from the simulations compare well with the solutions of Poisson-Boltzmann cell model, and the dendrimer's effective charge was estimated using the Bjerrum criterion. The influence of electrostatic interactions in the dendrimer's conformation due to repulsion between the charged end-groups and its relation to counterion effects is discussed for the several generations simulated. The form factors calculated from the simulations are compared with the model of a homogeneous ellipsoid of revolution. The overall results are in agreement with the previously established morphological transition of PAMAM dendrimers toward a more spherical and compact conformation above generations 3 or 4.     

Molecular dynamics simulations of the enzyme catechol-O-methyltransferase: methodological issues

Results from extensive 70 ns all-atom molecular dynamics simulations of catechol-O-methyltransferase (COMT) enzyme are reported. The simulations were performed with explicit TIP3P water and Mg2+ ions. Four different crystal structures of COMT, with and without different ligands, were used. These simulations are among the most extensive of their kind and as such served as a stability test for such simulations. On the methodological side we found that the initial energy minimization procedure may be a crucial step: particular hydrogen bonds may break, and this can initiate an irreversible loss of protein structure that becomes observable in longer time scales of the order of tens of nanoseconds. This has important implications for both molecular dynamics and quantum mechanics-molecular mechanics simulations.     

Photoionization dynamics of glycine adsorbed on a silicon cluster: "on-the-fly" simulations

Dynamics of glycine chemisorbed on the surface of a silicon cluster is studied for a process that involves single-photon ionization, followed by recombination with the electron after a selected time delay. The process is studied by "on-the-fly" molecular dynamics simulations, using the semiempirical parametric method number 3 (PM3) potential energy surface. The system is taken to be in the ground state prior to photoionization, and time delays from 5 to 50 fs before the recombination are considered. The time evolution is computed over 10 ps. The main findings are (1) the positive charge after ionization is initially mostly distributed on the silicon cluster. (2) After ionization the major structural changes are on the silicon cluster. These include Si-Si bond breaking and formation and hydrogen transfer between different silicon atoms. (3) The transient ionization event gives rise to dynamical behavior that depends sensitively on the ion state lifetime. Subsequent to 45 fs evolution in the charged state, the glycine molecule starts to rotate on the silicon cluster. Implications of the results to various processes that are induced by transient transition to a charged state are discussed. These include inelastic tunneling in molecular devices, photochemistry on conducting surfaces, and electron-molecule scattering.     

Molecular dynamics simulations of grafted polyelectrolytes on two apposing walls

Molecular dynamics simulations of polyelectrolytes grafted to two apposing surfaces were performed. Bead-spring polymer models are used to treat flexible chains [e.g., sodium poly(styrene sulfonate)] and stiff chains (double-stranded DNA). The counterions are explicitly treated. The effect of the surface density of the grafted polymer, the chain length, and the gap width on the structure and the pressure were studied. Results are compared to experimental measurements and to simulations of polyelectrolyte brushes on a single surface. The density profiles exhibit a maximum not found in single surface data. The maximum is due to the brushes shrinking to avoid interpenetration.     

Molecular dynamics simulations of polyelectrolyte multilayering on a charged particle

Molecular dynamics simulations of polyelectrolyte multilayering on a charged spherical particle revealed that the sequential adsorption of oppositely charged flexible polyelectrolytes proceeds with surface charge reversal and highlighted electrostatic interactions as the major driving force of layer deposition. Far from being completely immobilized, multilayers feature a constant surge of chain intermixing during the deposition process, consistent with experimental observations of extensive interlayer mixing in these films. The formation of multilayers as well as the extent of layer intermixing depends on the degree of polymerization of the polyelectrolyte chains and the fraction of charge on its backbone. The presence of ionic pairs between oppositely charged macromolecules forming layers seems to play an important role in stabilizing the multilayer film.     

Molecular dynamics simulations of pressure effects on hydrophobic interactions.

We report results on the pressure effects on hydrophobic interactions obtained from molecular dynamics simulations of aqueous solutions of methanes in water. A wide range of pressures that is relevant to pressure denaturation of proteins is investigated. The characteristic features of water-mediated interactions between hydrophobic solutes are found to be pressure-dependent. In particular, with increasing pressure we find that (1) the solvent-separated configurations in the solute-solute potential of mean force (PMF) are stabilized with respect to the contact configurations; (2) the desolvation barrier increases monotonically with respect to both contact and solvent-separated configurations; (3) the locations of the minima and the barrier move toward shorter separations; and (4) pressure effects are considerably amplified for larger hydrophobic solutes. Together, these observations lend strong support to the picture of the pressure denaturation process proposed previously by Hummer et al. (Proc. Natl. Acad. Sci. U.S.A. 1998, 95, 1552): with increasing pressure, the transfer of water into protein interior becomes key to the pressure denaturation process, leading to the dissociation of close hydrophobic contacts and subsequent swelling of the hydrophobic protein interior through insertions of water molecules. The pressure dependence of the PMF between larger hydrophobic solutes shows that pressure effects on the interaction between hydrophobic amino acids may be considerably amplified compared to those on the methane-methane PMF.     

Molecular dynamics simulations of liquid water: Voronoi polyhedra and network topology

We have carried out various runs at different temperatures simulating liquid water via the TIP4P model, a well tested and widely used 4-site intermolecular potential. Voronoi polyhedra have been constructed and the distributions of various quantities analyzed. The higher the temperature, the more the distributions are similar to those of a normal liquid: as a reference for a corresponding liquid without hydrogen bonds, we have examined hydrogen sulphide. The typical shape of the Voronoi polyhedra is not too far from the polyhedra of a perfect ice I_h crystal, although the numbers of vertices and faces are almost doubled. Indeed the actual number of vertices increases rapidly toward the liquid value as the crystal is perturbed by thermal motions. It is remarkable that the total surface and volume are always found closely correlated. It denotes a quasi-constant shape, i.e. the various shapes are not dependent on the volume. Indeed a direct analysis shows that there is no correlation between volume and shape. A quantitative characterization of the Voronoi polyhedra shape is proposed, through an “asphericity” parameter, whose distribution is found to be greatly affected by temperature: a net shift of the centre of the curve toward the value for the ice is found as the temperature is lowered. Finally it is found that there is a positive correlation between the potential energy of a molecule and the local volume, thus indicating that the hypothesis that the four-bonded molecules have a local density lower than average is incorrect.     

Molecular dynamics simulations for pure epsilon-CL-20 and epsilon-CL-20-based PBXs

Molecular dynamics has been employed to simulate the well-known high energy density compound epsilon-CL-20 (hexanitrohexaazaisowurtzitane) crystal and 12 epsilon-CL-20-based PBXs (polymer bonded explosives) with four kinds of typical fluorine polymers, i.e., polyvinylidenedifluoride, polychlorotrifluoroethylene, fluorine rubber (F(2311)), and fluorine resin (F(2314)) individually. The elastic coefficients, isotropic mechanical properties (tensile moduli, bulk moduli, shear moduli, and poission's ratios), and bonding energy are first reported for epsilon-CL-20 crystal and epsilon-CL-20-based polymer bonded explosives (PBXs). The mechanical properties of epsilon-CL-20 can be effectively improved by blending with a small amount of fluorine polymers, and the whole effect of the adding fluorine polymers to improve mechanical properties of PBXs along the three crystalline surfaces of epsilon-CL-20 is found to be (100) approximately (001) > (010). The interaction between each of the crystalline surfaces and each of the fluorine polymers is different, and the ordering of binding energy for the three surfaces is (001) > (100) > (010); F(2314) always has the strongest binding ability with the three different surfaces. F(2314) can best improve the ductibility and tenacity of PBX when it is positioned on epsilon-CL-20 (001) crystal surface. The calculations on detonation performances for pure epsilon-CL-20 crystal and the four epsilon-CL-20-based PBXs show that adding a small amount of fluorine polymer into pure epsilon-CL-20 will lower detonation performance, but each detonation parameter of the obtained PBXs is still excellent.     

Molecular dynamics implementation in MSINDO: study of silicon clusters

Born-Oppenheimer molecular dynamics is implemented in the semiempirical self-consistent field molecular orbital method MSINDO. The method is employed for the investigation of the structure and dynamics of silicon clusters of various sizes. The reliability of the present parameterization for silicon compounds is demonstrated by a comparison of the results of simulated annealing and of density functional calculations of Si(n) clusters (n = 5-7). The melting behavior of the Si(7) cluster is investigated and the MSINDO results are compared to previous high-level calculations. The efficiency of the present approach for the treatment of large systems is demonstrated by an extensive simulated annealing study of the Si(45) and Si(60) clusters. New Si(45) and Si(60) structures are found and evaluated. The relative stability of various energy minimum structures is compared with density functional calculations and available literature data.     

Molecular dynamics simulations of the enzyme Catechol-O-Methyltransferase: methodological issues1

Results from extensive 70 ns all-atom molecular dynamics simulations of catechol-O-methyltransferase (COMT) enzyme are reported. The simulations were performed with explicit TIP3P water and Mg^2?+ ions. Four different crystal structures of COMT, with and without different ligands, were used. These simulations are among the most extensive of their kind and as such served as a stability test for such simulations. On the methodological side we found that the initial energy minimization procedure may be a crucial step: particular hydrogen bonds may break, and this can initiate an irreversible loss of protein structure that becomes observable in longer time scales of the order of tens of nanoseconds. This has important implications for both molecular dynamics and quantum mechanics–molecular mechanics simulations.     

分子动力学模拟研究聚乳酸/聚酰胺11共混物的相容性

采用分子动力学(MD)模拟方法在COMPASS力场下,研究了不同质量比(10/90,30/70,50/50,70/30和90/10)聚乳酸(PLA)/聚酰胺11(PA11)共混物的相容性.研究结果表明:不同比例下PLA/PA11共混物的Gibbs自由能变化均大于零,其共混物很难形成均相体系;共混体系结合能的计算以及不同组分分子间C—C原子对径向分布函数的分析揭示了PLA和PA11的相互作用主要源自其分子间的范德华力;此外,模拟得到的所有比例下共混物的Flory-Huggins相互作用参数(χ)均大于临界Flory-Huggins相互作用参数(χcritical),进一步证明PLA与PA11不能形成相容体系。