Polymer glasses have very appealing properties for a number of optical and other applications and their high resistance to plastic deformation, expressed as the ratio of low temperature yield stress to elastic modulus of order 0.05, is an additional bonus. Their normalized resistance to crack propagation, as measured by their ratio of fracture energy Gc to 2γ, twice the surface energy, is more than 5000 for polymer glasses whose molecular weight exceeds 10 Me, where Me is their entanglement molecular weight. These are potentially very tough materials, which nevertheless have an important flaw, their propensity to deform plastically by crazing, which competes with shear deformation. Because crazes are highly localized zones of plastic deformation that can break down to form cracks at very low overall strains, crazing is a source of brittleness. Suppressing crazing and encouraging less localized shear deformation is thus an important molecular design goal. It has been known for some 20 years that two glassy polymer properties strongly influence whether crazing or shear deformation dominates, the physical aging of the polymer glass and its entangled strand density νe = ρNAv/Me, where ρ is the polymer strand mass density and NAv is Avogadro's number. Physical aging of the polymer glass increases the shear yield stress, as measured for example in compression, above the flow stress, the stress that is required to continue shear deformation at moderately larger strains. The flow stress has been shown to be independent of physical aging, giving rise to so-called strain softening, a decrease in (true) stress as the deformation continues past yield until the flow stress is reached. Recent experiments demonstrate unequivocally that a mild plastic predeformation by rolling erases the effects of physical aging (“rejuvenation”) so that if the polymer glass is deformed immediately after, even in tension, the yield stress and flow stress are the same. Shear deformation is favored over crazing even in polymers such as polystyrene (PS), which would normally craze in preference to deform by shear.1 Aging at room temperature for 48 h is sufficient to restore the normal yield and crazing behavior of PS. Obviously an understanding of the changes in structure of the glass during the initial stages of plastic deformation by shear is crucial to achieving a fundamental understanding of the competition between brittle (crazing) and ductile (shear) modes of deformation of polymer glasses, but to date very little is known except from molecular dynamics simulations of model Lennard-Jones nonpolymeric glasses under conditions of constant density at 0 K.2, 3 These simulations, viewed from the standpoint of energy landscape models, show that the rejuvenation by deformation is due to irreversible creation of structural states with shallower energy minima, that resemble, but are distinct from, those of the glass freshly quenched from an equilibrium melt, and that physical aging corresponds to the relaxation of these structural states toward those with deeper energy minima. This relaxation at constant density restores the yield point and strain softening behavior. Although previously this relaxation has been taken to be due to a decrease in free volume, the simulations show it occurs even in the absence of such a decrease. In addition, there is no experimental evidence of a volume increase during and immediately after moderate shear deformation that produces rejuvenation. In fact, a volume decrease has been reported.1 Recent positron annihilation experiments support this conclusion.4, 5 Simulations that treat the aging and rejuvenation of polymeric glasses, at constant pressure rather than constant density and at finite temperature, will certainly help us to clarify these important issues. Strain hardening modulus GR for PPO/PS blends and cross-linked PS as well as the rubbery plateau modulus G of PPO/PS and modulus Gx of rubbery cross-linked PS versus the strand density νe of the entanglement and cross-link density. Although both sets of moduli scale linearly with νe, the hardening modulus is two orders of magnitude larger than the true modulus of the entangled or cross-linked elastomer. (Data taken from Ref.6.) Although clearly experiments that can monitor the changes in strand dimensions with λ can help us understand how to fundamentally explain strain hardening in polymer glasses and its dependence on νe, I think that simulations can perhaps play an even more important role. In this regard the work of Rottler, Robbins, and coworkers suggests that it should now be possible to arrive at a molecular level understanding of this phenomenon.7 I therefore issue this as a challenge to the simulation community and look forward to a much improved physical picture for this important property of polymer glasses.