Related Researcher

Author's Photo

Baig, Chunggi
Theoretical and Computational study of Polymers & Nanomaterials Lab
Research Interests
  • Multiscale simulation, Polymer rheology, Nonequilibrium molecular dynamics/monte carlo


Analysis of the configurational temperature of polymeric liquids under shear and elongational flows using nonequilibrium molecular dynamics and Monte Carlo simulations

Cited 2 times inthomson ciCited 0 times inthomson ci
Analysis of the configurational temperature of polymeric liquids under shear and elongational flows using nonequilibrium molecular dynamics and Monte Carlo simulations
Baig, ChunggiEdwards, Brian J
Amorphous polyethylenes; Bending modes; Bond-stretching; Configurational entropy; Configurational temperature; Critical value; Dynamical effects; Elongational flow; Flow strength; Heat capacities; Interaction modes; Lennard jones; Linear polyethylene; Molecular collisions; MONTE CARLO; Monte Carlo Simulation; Non equilibrium; Nonequilibrium molecular dynamics; Nonequilibrium state; Nonmonotonic behaviors; Polymeric liquids; Polymeric material; Set-point temperatures; SHORT structures; Simulation data; Structural change; Structural characteristics
Issue Date
JOURNAL OF CHEMICAL PHYSICS, v.132, no.18, pp.184906-1 - 184906-19
We present a detailed analysis of the configurational temperature (T(conf)) for its application to polymeric materials using nonequilibrium molecular dynamics (NEMD) and nonequilibrium Monte Carlo (NEMC) methods. Simulations were performed of linear polyethylene liquid C(78)H(158) undergoing shear and elongational flows. At equilibrium, T(conf) is equal to the set point temperature of the simulation. An aphysically large decrease in T(conf) is observed in the NEMD simulations for both flows, especially at strong flow fields. By analyzing separately the individual contributions of the different potential interaction modes to the configurational temperature, it is found that the bonded modes (which constitutes almost 99.5% of the total) dominate the total T(conf) over the nonbonded ones; i.e., bond-stretching (approximate to 86.5%), bond-bending (approximate to 11.8%), bond-torsional (approximate to 1.2%), nonbonded intermolecular (approximate to 0.4%), and intramolecular (approximate to 0.1%) Lennard-Jones. The configurational temperature of the individual modes generally exhibits a nonmonotonic behavior with the flow strength and a dramatic change beyond a critical value of flow strength; this is mainly attributed to the dynamical effect of strong molecular collisions occurring at strong flow fields. In contrast, no such behavior is observed in the NEMC simulations where such dynamical effects are absent. Based on the principal physical concept of the configurational temperature, which represents the large-scale structural characteristics of the system, we propose to exclude the dynamical effects exhibited by the individual interaction modes, in obtaining a physically meaningful T(conf) as the configurational entropy of the system should not be affected by such factors. Since (a) the main difference between equilibrium and nonequilibrium states lies in the change in the overall (global) structure (represented by the bond torsional and nonbonded modes), and (b) the local, very short structure (represented by the bond-stretching and bond-bending modes) is barely changing between equilibrium and nonequilibrium states and its contribution to the total system configurational entropy is negligible compared to the large-scale structural changes, in order to accurately describe the structural changes occurring at nonequilibrium states by use of the configurational temperature, we further propose that only the contributions from the bond-torsional and nonbonded modes to Delta T(conf) between equilibrium and nonequilibrium states should be taken into account to generate a physically meaningful Delta T(conf). Applying the above hypothesis to the analysis of the simulation data, good agreement between the NEMD and NEMC simulations (and between NEMD simulations for different flows) is observed. Furthermore, the configurational temperature obtained in such way is found to match remarkably well with the heat capacity of amorphous polyethylene liquids and the flow-enhanced melting-point elevation reported in experiment.
Go to Link
Appears in Collections:
ECHE_Journal Papers
Files in This Item:
2-s2.0-77952730510.pdf Download

find_unist can give you direct access to the published full text of this article. (UNISTARs only)

Show full item record


  • mendeley


Items in DSpace are protected by copyright, with all rights reserved, unless otherwise indicated.