Molecular dynamics simulations are performed for a polymer melt composed of short chains in quiescent and sheared conditions. The stress relaxation function G(t) exhibits a stretched exponential form in a relatively early stage and ultimately follows the Rouse function in quiescent supercooled state. Transient stress evolution after application of shear obeys the linear growth t 0 dt ′ G(t ′ ) for strain less than 0.1 and then saturates into a non-Newtonian viscosity. In steady states, strong shear-thinning and elongation of chains into ellipsoidal shapes are found at extremely small shear. A glassy component of the stress is much enhanced in these examples.PACS numbers: 83.10. Nn, 83.20.Jp, 83.50.By, 64.70.Pf Stress and dielectric relaxations of glassy polymer melts occur from microscopic to macroscopic time scales in very complicated manners [1,2]. Experiments have shown that the stress relaxation function G(t) exhibits a glassy stretched exponential decay, a glass-rubber transition, a rubbery plateau, and a terminal decay, in this order over many decades of time. Such hierarchical relaxation behavior arises from rearrangements of jammed atomic configurations and subsequent evolution of chain conformations described by the Rouse or reptation dynamics [3,4]. The stress-optical relation between birefringence and stress has also been reported to be violated as the temperature T is approached the glass transition temperature T g [5][6][7], obviously owing to enhancement of a glassy part of the stress.Recent simulations on glassy polymer melts have mainly treated self-motions of particles in quiescent states [8][9][10]. However, not enough theoretical efforts have been made on the rheological properties of glassy polymers. Hence, we will first study linear rheology of a model short chain system via very long molecular dynamics simulations. Then, we will demonstrate that chains are very easily elongated at extremely small shear rateγ on the order of the inverse Rouse time. In our model all the bead particles interact with a Lennard-Jones potential of the form [4,8,10], U LJ (r) = 4ǫ[(σ/r) 12 − (σ/r) 6 ] + ǫ. It is cut off at the minimum distance 2 1/6 σ, so we use its repulsive part only to prevent spatial overlap of particles. Consecutive beads on each chain are connected by an anharmonic spring of the form, U F (r) = − 1 2 k c R 2 0 ln[1 − (r/R 0 ) 2 ] with k c = 30ǫ/σ 2 and R 0 = 1.5σ, so the bond length cannot exceed R 0 . In a cubic box with length L = 10σ under the periodic boundary condition, we put M = 100 chains composed of N = 10 beads. The number density is fixed at a very high value of n = N M/V = 1/σ 3 , which results in severely jammed configurations at low T . We will measure space and time in units of σ and τ 0 = (mσ 2 /ǫ) 1/2 with m being the mass of a bead. The temperature T will be measured in units of ǫ/k B . Simulations were performed in normal (T = 1.0) and supercooled (T = 0.4 and 0.2) states with and without shear. The bond lengths b j = |R j − R j+1 | between consecutive beads on each cha...