Mechanical stimuli in energetic materials initiate chemical reactions at shock fronts prior to detonation. Shock sensitivity measurements provide widely varying results, and quantum-mechanical calculations are unable to handle systems large enough to describe shock structure. Recent developments in reactive forcefield molecular dynamics (REAXFF-MD) combined with advances in parallel computing have paved the way to accurately simulate reaction pathways along with the structure of shock fronts. Our multimillionatom REAXFF-MD simulations of l,3,5-trinitro-l,3,5-triazine (RDX) reveal that detonation is preceded by a transition from a diffuse shock front with well-ordered molecular dipoles behind it to a disordered dipole distribution behind a sharp front. [12 -14]. A major unsolved problem in this area of mechanochemistry is how strains trigger chemical reactions that cause detonation. Recently, the role of a mechanical stimulus, such as bond bending, in initiating reaction at a shock front has attracted much attention. Manna et al. have shown on the basis of density functional calculation of a TATB crystal that bond-bending of nitro group in TATB molecule closes the HOMO-LUMO gap and have speculated that reaction may proceed at supersonic speed [15]. Atomistic-level understanding of the structure of shock fronts at the onset of detonation requires multimillion-atom simulations with reaction chemistry under high strain rate deformations.We have performed a series of molecular dynamics (MD) simulations of planar shock on a slab of an l,3,5-trinitro-l,3,5-triazine (RDX) [16,17] crystal. Each RDX molecule consists of 21 atoms [ Fig. 1(a)], and the unit cell of the RDX crystal contains 8 RDX molecules [ Fig. 1(b)]. The x, y, and z axes are aligned with the [100], [010], and [001] crystallographic orientations, respectively. We apply periodic boundary conditions in all directions after adding 2 nm of vacuum layers on both sides of the slab along the x axis. We have simulated a system of N 2 322 432 atoms with dimensions 318:48 284:08 271:76 A 3 and also a system of 145 152 atoms to confirm that major results do not change with the system size. Initially, the system is relaxed for 1 ps at 5 K and then quenched to 0 K.Simulations reported here are based on a scalable implementation of a reactive force-field (REAXFF) MD [18,19] on massively parallel computers [20]. The REAXFF potential energy function consists of bonding and nonbonding interactions. Bonding interactions comprise coordination energy, two-body stretching, three-body bending and fourbody torsion energy functions. Nonbonding interactions consists of van der Waals and Coulomb energies [21], where charge transfer is included by an electronegativity FIG. 1 (color). (a) An RDX molecule with carbon (yellow), hydrogen (white), oxygen (red), and nitrogen (blue) atoms. (b) The unit cell of an RDX crystal contains 8 RDX molecules, which are colored blue and red depending on whether the NO 2 groups faces away from (group 1) or faces towards (group 2) the shock plane. Th...