Effect Of Water On Brittle Fracture Of Sio2 By Molecular Dynamics Study

The influence of water on the brittle behavior of beta-cristobalite is studied by means of molecular dynamics (MD) simulation With the TTAM potential. Crack extension of mode 1 type is observed as the crack opening is filled LIP With water. The critical stress intensity factor K-lc(MD) is used to characterize the crack extension of MD simulation. The surface energy of SiO2 covered with layers of water is calculated at temperature of 300 K. Based oil the Griffith fracture criterion, the critical stress intensity factor K-lc(Griffith) is calculated, and it is in good agreement with that of MD simulation. (C) 2008 Elsevier B.V. All rights reserved.

A molecular dynamics study of the microscopic properties of simple dense fluids

The microscopic properties of a two-dimensional model dense fluid of Lennard-Jones disks have been studied using the so-called "molecular dynamics" method. Analyses of the computer-generated simulation data in terms of "conventional" thermodynamic and distribution functions verify the physical validity of the model and the simulation technique.

The radial distribution functions g(r) computed from the simulation data exhibit several subsidiary features rather similar to those appearing in some of the g(r) functions obtained by X-ray and thermal neutron diffraction measurements on real simple liquids. In the case of the model fluid, these "anomalous" features are thought to reflect the existence of two or more alternative configurations for local ordering.

Graphical display techniques have been used extensively to provide some intuitive insight into the various microscopic phenomena occurring in the model. For example, "snapshots" of the instantaneous system configurations for different times show that the "excess" area allotted to the fluid is collected into relatively large, irregular, and surprisingly persistent "holes". Plots of the particle trajectories over intervals of 2.0 to 6.0 x 10^-12 sec indicate that the mechanism for diffusion in the dense model fluid is "cooperative" in nature, and that extensive diffusive migration is generally restricted to groups of particles in the vicinity of a hole.

A quantitative analysis of diffusion in the model fluid shows that the cooperative mechanism is not inconsistent with the statistical predictions of existing theories of singlet, or self-diffusion in liquids. The relative diffusion of proximate particles is, however, found to be retarded by short-range dynamic correlations associated with the cooperative mechanism--a result of some importance from the standpoint of bimolecular reaction kinetics in solution.

A new, semi-empirical treatment for relative diffusion in liquids is developed, and is shown to reproduce the relative diffusion phenomena observed in the model fluid quite accurately. When incorporated into the standard Smoluchowski theory of diffusion-controlled reaction kinetics, the more exact treatment of relative diffusion is found to lower the predicted rate of reaction appreciably.

Finally, an entirely new approach to an understanding of the liquid state is suggested. Our experience in dealing with the simulation data--and especially, graphical displays of the simulation data--has led us to conclude that many of the more frustrating scientific problems involving the liquid state would be simplified considerably, were it possible to describe the microscopic structures characteristic of liquids in a concise and precise manner. To this end, we propose that the development of a formal language of partially-ordered structures be investigated.

Controlling electronic structures by irradiation in single-walled SiC nanotubes: a first-principles molecular dynamics study

Using first-principles molecular dynamics simulations, the displacement threshold energy and defect configurations are determined in SiC nanotubes. The simulation results reveal that a rich variety of defect structures (vacancies, Stone-Wales defects and antisite defects) are formed with threshold energies from 11 to 64 eV. The threshold energy shows an anisotropic behavior and exhibits a dramatic decrease with decreasing tube diameter. The electronic structure can be altered by the defects formed by irradiation, which suggests that the electron irradiation may be a way to use defect engineering to tailor electronic properties of SiC nanotubes.

Molecular dynamics study on the phase diagrams of linear and branched chain molecules

The particle transfer molecular dynamics is used to study the phase equilibria of linear and branched chain molecules. The scaling of the critical temperature versus chain length is obtained and the critical densities are found to decrease with increasing chain length, which are in agreement with the results of experiment and theory. The phase diagrams of the linear and the branched chain molecules nearly overlap with each other. Moreover, the radial distribution functions of linear and branched chain molecules in gas phase are very similar, but in the liquid phase, they are different for different kinds of chains.

Docking and molecular dynamics study on the inhibitory activity of coumarins on aldose reductase

In order to explore the inhibitory mechanism of coumarins toward aldose reductase (ALR2), AutoDock and Gromacs software were used for docking and molecular dynamics studies on 14 coumarins (CM) and ALR2 protease. The docking results indicate that residues TYR48, HIS110, and TRP111 construct the active pocket of ALR2 and, besides van der Waals and hydrophobic interaction, CM mainly interact with ALR2 by forming hydrogen bonds to cause inhibitory behavior. Except for CM1, all the other coumarins take the lactone part as acceptor to build up the hydrogen bond network with active-pocket residues. Unlike CM3, which has two comparable binding modes with ALR2, most coumarins only have one dominant orientation in their binding sites. The molecular dynamics calculation, based on the docking results, implies that the orientations of CM in the active pocket show different stabilities. Orientation of CM1 and CM3a take an unstable binding mode with ALR2; their conformations and RMSDs relative to ALR2 change a lot with the dynamic process. While the remaining CM are always hydrogen-bonded with residues TYR48 and HIS110 through the carbonyl O atom of the lactone group during the whole process, they retain the original binding mode and gradually reach dynamic equilibrium.

Kinetics of phase transformations in a model with metastable fluid-fluid separation: A molecular dynamics study

By molecular dynamics (MD) simulations we study the crystallization process in a model system whose particles interact by a spherical pair potential with a narrow and deep attractive well adjacent to a hard repulsive core. The phase diagram of the model displays a solid-fluid equilibrium, with a metastable fluid-fluid separation. Our computations are restricted to fairly small systems (from 2592 to 10368 particles) and cover long simulation times, with constant energy trajectories extending up to 76x10(6) MD steps. By progressively reducing the system temperature below the solid-fluid line, we first observe the metastable fluid-fluid separation, occurring readily and almost reversibly upon crossing the corresponding line in the phase diagram. The nucleation of the crystal phase takes place when the system is in the two-fluid metastable region. Analysis of the temperature dependence of the nucleation time allows us to estimate directly the nucleation free energy barrier. The results are compared with the predictions of classical nucleation theory. The critical nucleus is identified, and its structure is found to be predominantly fcc. Following nucleation, the solid phase grows steadily across the system, incorporating a large number of localized and extended defects. We discuss the relaxation processes taking place both during and after the crystallization stage. The relevance of our simulation for the kinetics of protein crystallization under normal experimental conditions is discussed. (C) 2002 American Institute of Physics.

Neutron diffraction, NMR and molecular dynamics study of glucose dissolved in the ionic liquid 1-ethyl-3-methylimidazolium acetate

beta-D-glucose dissolved in the ionic liquid 1-ethyl-3-methylimidazolium acetate in a 6 : 1 molar ratio (ionic liquid : glucose) has been studied by neutron scattering, NMR and molecular dynamics simulations. Good agreement was found between simulated neutron scattering profiles generated for isotopically substituted liquid systems and those experimentally determined as well as between simulated and experimental diffusion coefficients obtained by Pulsed Field Gradient NMR spectroscopy. The overriding glucose-ionic liquid interactions in the liquid are hydrogen-bonding between acetate oxygens and sugar hydroxyl groups. The ionic liquid cation was found to play only a minor role in the solvation of the sugar and does not participate in hydrogen-bonding with the sugar to any significant degree. NOESY experiments lend further evidence that there is no direct interaction between sugar hydroxyl groups and acidic hydrogens on the ionic liquid cation.

A Molecular Dynamics Study of the Early Stages of Calcium Carbonate Growth

The precipitation of calcium carbonate in water has been examined using a combination of molecular dynamics and umbrella sampling. During 20 ns molecular dynamics trajectories at elevated calcium carbonate concentrations, amorphous particles are observed to form and appear to be composed of misaligned domains of vaterite and aragonite. The addition of further calcium ions to these clusters is found to be energetically favorable and virtually barrierless. By contrast, there is a large barrier to the addition of calcium to small calcite crystals. Thus, even though calcite nanocrystals are stable in solution, at high supersaturations, particles of amorphous material form because this material grows much faster than ordered calcite nanocrystals.

Protection of DNA against low-energy electrons by amino acids: A first-principles molecular dynamics study

Using first-principles molecular dynamics simulations, we have investigated the notion that amino acids can play a protective role when DNA is exposed to excess electrons produced by ionizing radiation. In this study we focus on the interaction of glycine with the DNA nucleobase thymine. We studied thymine-glycine dimers and a condensed phase model consisting of one thymine molecule solvated in amorphous glycine. Our results show that the amino acid acts as a protective agent for the nucleobase in two ways. If the excess electron is initially captured by the thymine, then a proton is transferred in a barrier-less way from a neighboring hydrogen-bonded glycine. This stabilizes the excess electron by reducing the net partial charge on the thymine. In the second mechanism the excess electron is captured by a glycine, which acts as a electron scavenger that prevents electron localization in DNA. Both these mechanisms introduce obstacles to further reactions of the excess electron within a DNA strand, e.g. by raising the free energy barrier associated with strand breaks.

Elementary steps in aqueous proton transfer reactions : a first principles molecular dynamics study

La nature des acides dans un environnement aqueux est primordiale dans de nombreux aspects de la chimie et de la biologie. La caractéristique principale d'un acide est sa capacité à transférer un proton vers une molécule d'eau ou vers n'importe quelle base, mais ce procédé n'est pas aussi simple qu'il y paraît. Il peut au contraire être extrêmement complexe et dépendre de manière cruciale de la solvatation des différents intermédiaires de réaction impliqués. Cette thèse décrit les études computationnelles basées sur des simulations de dynamique moléculaire ab initio qui ont pour but d'obtenir une description à l'échelle moléculaire des divers procédés de transferts de proton entre acide et bases dans un milieu aqueux. Pour cela, nous avons étudié une serie de système, dont l'acide hydrofluorique aqueux, l'acide trifluoroacétique aqueux, et un système modèle constitué d'un phénol et d'une entité carboxylate reliés entre eux par une molécule d'eau en solution aqueuse. Deux états intermédiaires ont été identifiés pour le transfert d'un proton depuis un acide. Ces intermédiaires apparaissent stabilisés par un motif local de solvatation via des ponts H. Leurs signatures spectroscopiques ont été caractérisées au moyen de la spectroscopie infrarouge, en utilisant le formalisme de la dynamique moléculaire ab initio, qui inclut l'effet quantique nucléaire de manière explicite. Cette étude a aussi identifié trois chemins de réaction élémentaire, qui sont responsable pour le transfert d'un proton d'un acide à une base, ainsi que leurs échelles de temps caractéristiques. Les conclusions tirées de ces études sont discutées dans les détails, au niveau moléculaire, avec une emphase sur les comparaisons entre les résultats théoriques et les mesures expérimentales obtenues dans a littérature ou via des collaborateurs.