Hydrogen irradiation-driven computational surface chemistry of lithium oxide and hydroxide.

We have investigated, using molecular dynamics, the surface chemistry of hydrogen incident on the amorphous and crystalline lithium oxide and lithium hydroxide surfaces upon being slowed down by a collision cascade and retained in the amorphous surface of either Li2O or LiOH. We looked for the bonding of H to the resident structures in the surface to understand a possible chain of chemical reactions that can lead to surface transformation upon H atom impact. Our findings, using Density-Functional Theory (DFT) trained ReaxFF force field/electronegativity equalization method potentials, stress the importance of inclusion of polarization in the dynamics of a Li-O-H system, which is also illustrated by DFT energy minimization and quantum-classical molecular dynamics using tight binding DFT. The resulting polar-covalent chemistry of the studied systems is complex and very sensitive to the instantaneous positions of all atoms as well as the ratio of concentrations of various resident atoms in the surface.

ReaxFF reactive force field for molecular dynamics simulations of liquid Cu and Zr metals.

We develop a ReaxFF reactive force field used for the molecular dynamics simulations of thermophysical properties of liquid Cu and Zr metals. The ReaxFF parameters are optimized by fitting to the first-principles density-functional calculations on the equations of state for bulk crystal structures and surface energies. To validate the force field, we compare the ReaxFF results with those from experiments and embedded-atom-method (EAM) potentials. We demonstrate that the present ReaxFF force field well represents structural characteristics and diffusion behaviors of elemental Cu and Zr up to high-temperature liquid regions. It reasonably reproduces the thermodynamic processes associated with crystal-liquid interface. In particular, the equilibrium melting temperatures show better agreement with experimental measurements than the results from EAM potentials. The ReaxFF reactive force field method exhibits a good transferability to the nonreactive processes of liquid systems.

Study of thermal conductivity of ice clusters after impact deposition on the silica surfaces using the ReaxFF reactive force field.

During aircraft or spacecraft missions, ice accumulates on different parts of their surface elements. An important parameter affecting the ability to remove this ice from the surface is the heat transfer characteristics of the accumulated ice. The ice heat transfer is related to the process of ice formation and its density and internal structure. In this study we investigate the effects of the ice and silica structure and the ice cluster attachment mechanism to the silica surface on the thermal conductivity (TC) of the attached ice cluster using the ReaxFF reactive force field. The purpose of this study is to investigate the thermal transport in amorphous and crystalline ice after high-velocity deposition on the silica surfaces. A dual thermostat method has been applied for the calculation of TC values. The validity of this method has been verified by comparing the calculated values of TC for crystal and amorphous ice with available experimental values. Our calculations show that the TC values of both crystal and amorphous ice drop after deposition on the silica surfaces. This decrease in the TC is more significant for the ice deposition on suboxide silica surfaces. Furthermore, crystal ice shows higher TC values than amorphous ice after accumulation. However, when crystal ice impacts on the silica surface at 1 km s(-1) impact speed, the crystalline shape of the ice cluster is lost to a considerable level and the TC values obtained for the ice clusters in such cases are closer to amorphous ice TC values. We observed a decrease in the TC values when ionic species are added inside the ice clusters.

Reactive molecular dynamics study of the pH-dependent dynamic structure of α-helix.

We have studied the α-helix to random coil transition using ReaxFF reactive molecular dynamics as a function of pH. Urea binding to peptides and associated interference with backbone H-bonds and charged side chains interactions, which can both denature the helices, have been studied previously using nonreactive force fields (Topol, I. A. J. Am. Chem. Soc. 2001, 123, 6054-6060). This study reveals new proton-transfer mechanisms related to the denaturation of α-helical structures, which cannot be captured by nonreactive molecular dynamics. In addition, we show that proton transfer between the solution and the peptide can break the α-helix hydrogen bonds, and consequently, at extreme pHs, a significant amount of helix will unravel. We also compare the effects of temperature in the denaturation mechanism. The ReaxFF findings are in significantly better agreement with ab initio calculations than previous nonreactive force field results, indicating the relevance of the reactive component on helical loss.

Mechanical properties and fracture dynamics of silicene membranes.

As graphene has become one of the most important materials, there is renewed interest in other similar structures. One example is silicene, the silicon analogue of graphene. It shares some of the remarkable graphene properties, such as the Dirac cone, but presents some distinct ones, such as a pronounced structural buckling. We have investigated, through density functional based tight-binding (DFTB), as well as reactive molecular dynamics (using ReaxFF), the mechanical properties of suspended single-layer silicene. We calculated the elastic constants, analyzed the fracture patterns and edge reconstructions. We also addressed the stress distributions, unbuckling mechanisms and the fracture dependence on the temperature. We analysed the differences due to distinct edge morphologies, namely zigzag and armchair.

Reactive molecular dynamics simulation on the disintegration of Kapton, POSS polyimide, amorphous silica, and teflon during atomic oxygen impact using the ReaxFF reactive force-field method.

Atomic oxygen (AO) is the most abundant element in the low Earth orbit (LEO). It is the result of the dissociation of molecular oxygen by ultraviolet radiation from the sun. In the LEO, it collides with the materials used on spacecraft surfaces and causes degradation of these materials. The degradation of the materials on the surface of spacecrafts at LEO has been a significant problem for a long time. Kapton polyimide, polyhedral oligomeric silsesquioxane (POSS), silica, and Teflon are the materials extensively used in spacecraft industry, and like many other materials used in spacecraft industry, AO collision degradation is an important issue in their applications on spacecrafts. To investigate the surface chemistry of these materials in exposure to space AO, a computational chemical evaluation of the Kapton polyimide, POSS, amorphous silica, and Teflon was performed in separate simulations under similar conditions. For performing these simulations, the ReaxFF reactive force-field program was used, which provides the computational speed required to perform molecular dynamics (MD) simulations on system sizes sufficiently large to describe the full chemistry of the reactions. Using these simulations, the effects of AO impact on different materials and the role of impact energies, the content of material, and temperature of material on the behavior of the materials are studied. The ReaxFF results indicate that Kapton is less resistant than Teflon toward AO damage. These results are in good agreement with experiment. These simulations indicate that the amorphous silica shows the highest stability among these materials before the start of the highly exothermic silicon oxidation. We have verified that adding silicon to the bulk of the Kapton structure enhances the stability of the Kapton against AO impact. Our canonical MD simulations demonstrate that an increase in the heat transfer in materials during AO impact can provide a considerable decrease in the disintegration of the material. This effect is especially relevant in silica AO collision. Considerable experimental efforts have been undertaken to minimize such AO-based degradations. As our simulations demonstrate, ReaxFF can provide a cost-effective screening tool for future material optimization.

Defect healing and enhanced nucleation of carbon nanotubes by low-energy ion bombardment.

Structural defects inevitably appear during the nucleation event that determines the structure and properties of single-walled carbon nanotubes. By combining ion bombardment experiments with atomistic simulations we reveal that ion bombardment in a suitable energy range allows these defects to be healed resulting in an enhanced nucleation of the carbon nanotube cap. The enhanced growth of the nanotube cap is explained by a nonthermal ion-induced graphene network restructuring mechanism.

Graphene to fluorographene and fluorographane: a theoretical study.

We report here a fully reactive molecular dynamics study on the structural and dynamical aspects of the fluorination of graphene membranes (fluorographene). Our results show that fluorination tends to produce defective areas on the graphene membranes with significant distortions of carbon-carbon bonds. Depending on the amount of incorporated fluorine atoms, large membrane holes were observed due to carbon atom losses. These results may explain the broad distribution of the structural lattice parameter values experimentally observed. We have also investigated the effects of mixing hydrogen and fluorine atoms on the graphene functionalization. Our results show that, when in small amounts, the presence of hydrogen atoms produces a significant decrease in the rate of fluorine incorporation onto the membrane. On the other hand, when fluorine is the minority element, it produces a significant catalytic effect on the rate of hydrogen incorporation. We have also observed the spontaneous formation of new hybrid structures with different stable configurations (chair-like, zigzag-like and boat-like) which we named fluorographane.

Reactive molecular dynamics of the initial oxidation stages of Ni111 in pure water: effect of an applied electric field.

Corrosion processes occurring in aqueous solutions are critically dependent upon the interaction between the metal electrode and the solvent. In this work, the interaction of a nickel substrate with water molecules has been investigated using reactive force field (ReaxFF) molecular dynamics simulations. This approach was originally developed by van Duin and co-workers to study hydrocarbon chemistry and the catalytic properties of organic compounds. To our knowledge, this method has not previously been used to study the corrosion of nickel. In this work, we studied the interaction of 480 molecules of water (ρ = 0.99 g·cm(-3)) with Ni(111) surfaces at 300 K. The results showed that a water "bilayer" was adsorbed on the nickel surface. In the absence of an applied electric field, no dissociation of water was observed. However, the nickel atoms at the surface were charged positively, whereas the first water layer was charged negatively, indicating the formation of an electric double layer. To study the corrosion of nickel in pure water, we introduced an external electric field between the metal and the solution. The electric field intensity varied between 10 and 20 MeV/cm. The presence of this electric field led to oxidation of the metal surface. The structural and morphological differences associated with the growth of this oxide film in the presence of the electric field were evaluated. The simulated atomic trajectories were used to analyze the atomic displacement during the reactive process. The growth of the oxide scale on the nickel surface was primarily due to the movement of anions toward the interior of the metal substrate and the migration of nickel toward the free surface. We found that increasing the electric field intensity sped up the corrosion of nickel. The results also showed that the oxide film thickness increased linearly with increasing electric field intensity.

Can amorphization take place in nanoscale interconnects?

The trend of miniaturization has highlighted the problems of heat dissipation and electromigration in nanoelectronic device interconnects, but not amorphization. While amorphization is known to be a high pressure and/or temperature phenomenon, we argue that defect density is the key factor, while temperature and pressure are only the means. For nanoscale interconnects carrying modest current density, large vacancy concentrations may be generated without the necessity of high temperature or pressure due to the large fraction of grain boundaries and triple points. To investigate this hypothesis, we performed in situ transmission electron microscope (TEM) experiments on 200 nm thick (80 nm average grain size) aluminum specimens. Electron diffraction patterns indicate partial amorphization at modest current density of about 10(5) A cm(-2), which is too low to trigger electromigration. Since amorphization results in drastic decrease in mechanical ductility as well as electrical and thermal conductivity, further increase in current density to about 7 × 10(5) A cm(-2) resulted in brittle fracture failure. Our molecular dynamics (MD) simulations predict the formation of amorphous regions in response to large mechanical stresses (due to nanoscale grain size) and excess vacancies at the cathode side of the thin films. The findings of this study suggest that amorphization can precede electromigration and thereby play a vital role in the reliability of micro/nanoelectronic devices.

Reactive molecular dynamics study on the first steps of DNA damage by free hydroxyl radicals.

We employ a large scale molecular simulation based on bond-order ReaxFF to simulate the chemical reaction and study the damage to a large fragment of DNA molecule in the solution by ionizing radiation. We illustrate that the randomly distributed clusters of diatomic OH radicals that are primary products of megavoltage ionizing radiation in water-based systems are the main source of hydrogen abstraction as well as formation of carbonyl and hydroxyl groups in the sugar moiety that create holes in the sugar rings. These holes grow up slowly between DNA bases and DNA backbone, and the damage collectively propagates to a DNA single and double strand break.

Structure and dynamics of shock-induced nanobubble collapse in water.

Shock-induced collapse of nanobubbles in water is investigated with molecular dynamics simulations based on a reactive force field. We observe a focused jet at the onset of bubble shrinkage and a secondary shock wave upon bubble collapse. The jet length scales linearly with the nanobubble radius, as observed in experiments on micron-to-millimeter size bubbles. Shock induces dramatic structural changes, including an ice-VII-like structural motif at a particle velocity of 1 km/s. The incipient ice VII formation and the calculated Hugoniot curve are in good agreement with experimental results.