Estimation of relaxation modulus of polymer melts by molecular dynamics simulations: Application of relaxation mode analysis

A framework for estimating the linear relaxation modulus of polymer melts by molecular dynamics (MD) simulations is presented on the basis of relaxation mode analysis (RMA). Conventional calculations of the relaxation modulus based on the Green-Kubo formula are computationally very expensive owing to long relaxation times and poor convergence of stress autocorrelation functions: In practice, reliable calculations usually require a time average over O(10τ₁-100τ₁), where τ₁ is the longest relaxation time of a chain in the melt. RMA is a method that systematically extracts relaxation modes and rates of a polymer chain from the time correlation functions of coordinates of polymer segments. In the present method, the relaxation modulus is evaluated by fitting the data of the stress autocorrelation functions to the generalized Maxwell model whose relaxation times are determined from the relaxation rate spectrum obtained by RMA.