Influence of interatomic bonding potentials on detonation properties.

The dependences of the macroscopic detonation properties of a two-dimensional (2D) diatomic (AB) molecular system on the fundamental molecular properties were investigated. This includes examining the detonation velocity, reaction zone thickness, and critical width as functions of the exothermicity (Q) of the gas-phase reaction [AB --> (1/2)(A(2) + B(2))] and the gas-phase dissociation energy (D(e)(AB)) for AB --> A + B . Following previous work, molecular dynamics (MD) simulations with a reactive empirical bond-order potential were used to characterize the shock-induced response of a diatomic AB molecular solid, which exothermically reacts to produce A2 and B2 gaseous products.

Shock waves in polycrystalline iron.

The propagation of shock waves through polycrystalline iron is explored by large-scale atomistic simulations. For large enough shock strengths the passage of the wave causes the body-centered-cubic phase to transform into a close-packed phase with most structure being isotropic hexagonal-close-packed (hcp) and, depending on shock strength and grain orientation, some fraction of face-centered-cubic (fcc) structure. The simulated shock Hugoniot is compared to experiments.

Direct observation of the alpha-epsilon transition in shock-compressed iron via nanosecond x-ray diffraction.

In situ x-ray diffraction studies of iron under shock conditions confirm unambiguously a phase change from the bcc (alpha) to hcp (epsilon) structure. Previous identification of this transition in shock-loaded iron has been inferred from the correlation between shock-wave-profile analyses and static high-pressure x-ray measurements. This correlation is intrinsically limited because dynamic loading can markedly affect the structural modifications of solids.

Nanohydrodynamics simulations: an atomistic view of the Rayleigh-Taylor instability.

Nanohydrodynamics simulations, hydrodynamics on the nanometer and nanosecond scale by molecular dynamics simulations for up to 100 million particles, are performed on the latest generation of supercomputers. Such simulations exhibit Rayleigh-Taylor instability, the mixing of a heavy fluid on top of a light in the presence of a gravitational field, initiated by thermal fluctuations at the interface, leading to the chaotic regime in the long-time evolution of the mixing process. The early-time behavior is in general agreement with linear analysis of continuum theory (Navier-Stokes), and the late-time behavior agrees quantitatively with experimental observations.

Microscopic view of structural phase transitions induced by shock waves.

Multimillion-atom molecular-dynamics simulations are used to investigate the shock-induced phase transformation of solid iron. Above a critical shock strength, many small close-packed grains nucleate in the shock-compressed body-centered cubic crystal growing on a picosecond time scale to form larger, energetically favored grains. A split two-wave shock structure is observed immediately above this threshold, with an elastic precursor ahead of the lagging transformation wave.

Uniaxial Hugoniostat: a method for atomistic simulations of shocked materials.

An new equilibrium molecular-dynamics method (the uniaxial Hugoniostat) is proposed to study the energetics and deformation structures in shocked crystals. This method agrees well with nonequilibrium molecular-dynamics simulations used to study shock-wave propagation in solids and liquids.

Orientation dependence in molecular dynamics simulations of shocked single crystals.

We use multimillion-atom molecular dynamics simulations to study shock wave propagation in fcc crystals. As shown recently, shock waves along the <100> direction form intersecting stacking faults by slippage along 111 close-packed planes at sufficiently high shock strengths. We find even more interesting behavior of shocks propagating in other low-index directions: for the <111> case, an elastic precursor separates the shock front from the slipped (plastic) region.

Plasticity induced by shock waves in nonequilibrium molecular-dynamics simulations.

Nonequilibrium molecular-dynamics simulations of shock waves in three-dimensional 10-million atom face-centered cubic crystals with cross-sectional dimensions of 100 by 100 unit cells show that the system slips along all of the available 111 slip planes, in different places along the nonplanar shock front. Comparison of these simulations with earlier ones on a smaller scale not only eliminates the possibility that the observed slippage is an artifact of transverse periodic boundary conditions, but also reveals the richness of the nanostructure left behind. By introducing a piston face that is no longer perfectly flat, mimicking a line or surface inhomogeneity in the unshocked material, it is shown that for weaker shock waves (below the perfect-crystal yield strength), stacking faults can be nucleated by preexisting extended defects.

Large-scale molecular dynamics simulations of dislocation intersection in copper.

The results of massively parallel three-dimensional molecular dynamics simulations of the perpendicular intersection of extended dislocations in copper are reported. The intersection process, which involves three of the four possible 111 glide planes in the face-centered cubic lattice, begins with junction formation, followed by unzipping, partial dislocation bowing, cutting, and, finally, unit jog formation. The investigation provides insights into this complex atomistic process, which is currently not accessible to experimental investigation.