It is shown that the Hugoniot and the critical shear stress required to deform a metal plastically in shock compression can be obtained directly from molecular dynamics simulations without recourse to surface velocity profiles, or to details of the dislocation evolution. Specific calculations are shown for aluminum shocked along the [100] direction, and containing an initial distribution of microscopic defects. The presence of such defects has a minor effect on the Hugoniot and on the dynamic strength at high pressures. Computed results agree with experimental data.Introduction.-The carriers of plastic deformation in crystalline solids are dislocations. To set them in motion requires shear stresses above a critical value, called the yield or plastic flow strength, or strength for short. It depends on the type of crystal, its orientation with respect to the applied forces, on the temperature, on any previous plastic deformation, and above all on the rate of deformation enforced by the applied forces. In general, the strength increases as the rate of deformation increases. To explore strength at extreme plastic strain rates of 10 4 and beyond, dynamic forces must be applied by impact of solids, by chemical explosives, pulsed magnetic fields, or laser ablation. The rapid application of these extreme forces induce shock waves that propagate from the impacted surface through the material at speeds on the order of the sound velocity, typically several km/s. Compression of the material and the ensuing plastic flow occur behind the shock front [1]. When the shock reaches the exit surface, its velocity can be measured as a function of time, and from the analysis of this velocity profile, the degree of compression or the density, the pressure, and the dynamic strength can be inferred. However, direct measurements of all three parameters behind the shock front are not possible.It has therefore been a goal for many years to employ atomistic simulations and study strong shock waves in solids, visualize the dynamic processes during shock compression and deformation, and also "measure", as it were, the thermodynamic field variables such as temperature, density, pressure, flow velocities, and stresses. Holian, one of the prominent pioneers in this field, has reviewed the historical development and the remaining challenges to reach this goal [2]. While the molecular dynamic simulations have provided numerous visualizations of dynamic dislocation behavior during shock compression, the complete set of thermodynamic field variables has so far not been obtained.Here, we give a brief account of having accomplished this last step, namely a complete mapping of the spatial and temporal distributions of all the thermodynamic variables: temperature, stresses and their rates, strains and their rates, as well as derivative parameters such as pressure, density, von Mises stress, and plastic strain and strain rates. Counting all components of tensors, this is a total of 40 distribution functions. With this information in hand, it is now possi...