Atomistic molecular dynamics simulation of stress relaxation upon cessation of steady-state uniaxial elongational flow

A new approach is presented for predicting the linear viscoelastic properti es of a polymer melt through a series of molecular dynamics (MD) simulation s of the relaxation of well-equilibrated, preoriented, strained configurati ons. Such strained configurations have been accumulated (Mavrantzas and The odorou, 1998) by employing the end-bridging Monte Carlo (EBMC) algorithm in the presence of a small tensorial field alpha (xx) which orients the chain s. They are representative of a melt under conditions of steady-state, unia xial elongational flow. In the dynamic studies presented in this work, the tensorial field alpha (xx) is removed and the relaxation of the system back toward its field-free, equilibrium state is monitored with MD. All simulat ions are performed in the NTL(x)sigma (yy)sigma (zz) statistical ensemble ( Yang et al., 1997), where the following variables are kept constant: the nu mber of particles N, the temperature T, the length L-x of the simulation bo x in the direction of flow, and the average stress (sigma (yy) + sigma (zz) )/2 in the other two (lateral) directions. The physical experiment modeled through these NTL(x)sigma (yy)sigma (zz) MD simulations is one of stress re laxation upon cessation of a steady-state, uniaxial elongational flout. The relaxation of the melt is quantified by monitoring the temporal evolution of the normal stress sigma (xx) in the x direction, of the volume V, and of certain descriptors of the short- and long-length scale conformation of ch ains. These include the diagonal components (c) over tilde (xx), (c) over t ilde (yy), and (c) over tilde (zz) of the conformation tensor and the chain mean-square end-to-end distance [R-2], all functions of time t. Results fo r the aforementioned properties, accumulated as statistical averages over m any initial configurations subjected to NTL(x)sigma (yy)sigma (zz) MD simul ations, are presented for two PE melt systems, C-24 and C-78, both of which have been studied extensively in the past. By invoking the Rouse model, an alytical expressions are derived for the functions sigma (xx)(t) and (c) ov er tilde (xx)(t) corresponding to the experiment simulated. By mapping the simulation results on these expressions, the longest relaxation times tau ( R) of the melts are extracted in excellent agreement with previous equilibr ium computer experiments (Harmandaris et al.; 1998). The stress relaxation modulus G(t) is computed from the equilibrium shear stress autocorrelation function at short times and from the spectrum of relaxation times extracted by mapping the MD results on the Rouse model at long times, yielding consi stent and physically meaningful results.