Atomistic force field for alumina fit to density functional theory.

We present a force field for bulk alumina (Al2O3), which has been parametrized by fitting the energies, forces, and stresses of a large database of reference configurations to those calculated with density functional theory (DFT). We use a functional form that is simpler and computationally more efficient than some existing models of alumina parametrized by a similar technique. Nevertheless, we demonstrate an accuracy of our potential that is comparable to those existing models and to DFT.

Melting slope of MgO from molecular dynamics and density functional theory.

We combine density functional theory (DFT) with molecular dynamics simulations based on an accurate atomistic force field to calculate the pressure derivative of the melting temperature of magnesium oxide at ambient pressure--a quantity for which a serious disagreement between theory and experiment has existed for almost 15 years. We find reasonable agreement with previous DFT results and with a very recent experimental determination of the slope. We pay particular attention to areas of possible weakness in theoretical calculations and conclude that the long-standing discrepancy with experiment could only be explained by a dramatic failure of existing density functionals or by flaws in the original experiment.

Giant wave-drag enhancement of friction in sliding carbon nanotubes.

Molecular dynamics simulations of coaxial carbon nanotubes in relative sliding motion reveal a striking enhancement of friction when phonons whose group velocity is close to the sliding velocity of the nanotubes are strongly excited. The effect is analogous to the dramatic increase in air drag experienced by aircraft flying close to the speed of sound but differs in that it can occur in multiple velocity ranges with varying magnitude, depending on the atomic level structures of the nanotubes. The phenomenon is a general one that may occur in other nanoscale mechanical systems.

On the theory underlying the Car-Parrinello method and the role of the fictitious mass parameter.

The theory underlying the Car-Parrinello extended-Lagrangian approach to ab initio molecular dynamics (CPMD) is reviewed and reexamined using "heavy" ice as a test system. It is emphasized that the adiabatic decoupling in CPMD is not a decoupling of electronic orbitals from the ions but only a decoupling of a subset of the orbital vibrational modes from the rest of the necessarily coupled system of orbitals and ions. Recent work [J.

Structural transformations of carbon nanotubes under hydrostatic pressure.

We used simulations with a classical force field to study the transformation under hydrostatic pressure of isolated single-walled nanotubes (SWNT) from a circular to a collapsed cross section. Small-diameter SWNTs deform continuously under pressure, whereas larger-diameter SWNTs display hysteresis and undergo a first-order-like transformation. The different behavior is due to the changing proportions in the total energy of the wall-curvature energy and the van der Waals attraction between opposite walls of the tube.

Temperature dependence of the band gap of semiconducting carbon nanotubes.

The temperature dependence of the band gap of semiconducting single-wall carbon nanotubes (SWNTs) is calculated by direct evaluation of electron-phonon couplings within a "frozen-phonon" scheme. An interesting diameter and chirality dependence of E(g)(T) is obtained, including nonmonotonic behavior for certain tubes and distinct "family" behavior. These results are traced to a strong and complex coupling between band-edge states and the lowest-energy optical phonon modes in SWNTs.

Dynamic sliding friction between concentric carbon nanotubes.

Molecular dynamics simulations are used to study mechanical energy dissipation in carbon nanotube oscillators of lengths of tens of nanometers. The principal source of friction is found to be the ends of the tubes and hence dynamical friction is virtually independent of the overlap area between tubes. As a result of this, tube commensuration does not lead to significantly increased frictional forces.

Pressure-induced structural changes in liquid SiO2 from Ab initio simulations.

First-principles molecular dynamics simulations at constant pressure have been used to investigate the mechanisms of compression of liquid SiO2. Liquid SiO2 is found to become denser than quartz at a pressure of about 6 GPa, in agreement with extrapolations of lower pressure experimental data. The high compressibility of the liquid is traced to medium-range changes in the topology of the atomic network.