skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.

Explore Research Products in the PAR
It may take a few hours for recently added research products to appear in PAR search results.


Title: Chemical Potential of a Flexible Polymer Liquid in a Coarse-Grained Representation
While the excess chemical potential is the key quantity in determining phase diagrams, its direct computation for high-density liquids of long polymer chains has posed a significant challenge. Computationally, the excess chemical potential is calculated using the Widom insertion method, which involves monitoring the change in internal energy as one incrementally introduces individual molecules into the liquid. However, when dealing with dense polymer liquids, inserting long chains requires generating trial configurations with a bias that favors those at low energy on a unit-by-unit basis: a procedure that becomes more challenging as the number of units increases. Thus, calculating the excess chemical potential of dense polymer liquids using this method becomes computationally intractable as the chain length exceeds N ≥ 30. Here, we adopt a coarsegrained model derived from the integral equation theory for which inserting long polymer chains becomes feasible. The integral equation theory of coarse graining (IECG) represents a polymer as a sphere or a collection of blobs interacting through a soft potential. We employ the IECG approach to compute the excess chemical potential using Widom’s method for polymer chains of increasing lengths, extending up to N = 720 monomers, and at densities reaching up to ρ = 0.767 g/cm3. From a fundamental perspective, we demonstrate that the excess chemical potentials remain nearly constant across various levels of coarse graining, offering valuable insights into the consistency of this type of procedure. Ultimately, we argue that current Monte Carlo algorithms, originally designed for atomistic simulations, such as configurational bias Monte Carlo (CBMC) methods, can significantly benefit from the integration of the IECG approach, thereby enhancing their performance in the study of phase diagrams of polymer liquids.  more » « less
Award ID(s):
2154999
PAR ID:
10520057
Author(s) / Creator(s):
; ; ;
Publisher / Repository:
ACS
Date Published:
Journal Name:
The Journal of Physical Chemistry B
Volume:
128
ISSN:
1520-6106
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; (, Soft Matter)
    The integral equation coarse-graining (IECG) approach is a promising high-level coarse-graining (CG) method for polymer melts, with variable resolution from soft spheres to multi CG sites, which preserves the structural and thermodynamical consistencies with the related atomistic simulations. When compared to the atomistic description, the procedure of coarse-graining results in smoother free energy surfaces, longer-ranged potentials, a decrease in the number of interaction sites for a given polymer, and more. Because these changes have competing effects on the computational efficiency of the CG model, care needs to be taken when studying the effect of coarse-graining on the computational speed-up in CG molecular dynamics simulations. For instance, treatment of long-range CG interactions requires the selection of cutoff distances that include the attractive part of the effective CG potential and force. In particular, we show how the complex nature of the range and curvature of the effective CG potential, the selection of a suitable CG timestep, the choice of the cutoff distance, the molecular dynamics algorithms, and the smoothness of the CG free energy surface affect the efficiency of IECG simulations. By direct comparison with the atomistic simulations of relatively short chain polymer melts, we find that the overall computational efficiency is highest for the highest level of CG (soft spheres), with an overall improvement of the computational efficiency being about 10 6 –10 8 for various CG levels/resolutions. Therefore, the IECG method can have important applications in molecular dynamics simulations of polymeric systems. Finally, making use of the standard spatial decomposition algorithm, the parallel scalability of the IECG simulations for various levels of CG is presented. Optimal parallel scaling is observed for a reasonably large number of processors. Although this study is performed using the IECG approach, its results on the relation between the level of CG and the computational efficiency are general and apply to any properly-constructed CG model. 
    more » « less
  2. ; (, Minerals)
    We explore the structural properties of Mg, MgO, and MgSiO3 liquids from ab initio computer simulations at conditions that are relevant for the interiors of giant planets, stars, shock compression measurements, and inertial confinement fusion experiments. Using path-integral Monte Carlo and density functional theory molecular dynamics, we derive the equation of state of magnesium-rich liquids in the regime of condensed and warm dense matter, with densities ranging from 0.32 to 86.11 g cm−3 and temperatures from 20,000 K to 5 × 108 K. We study the electronic structure of magnesium as a function of density and temperature and the correlations of the atomic motion, finding an unexpected local maximum in the pair correlation functions that emerges at high densities which decreases the coordination number of elemental magnesium and reveals a higher packing. This phenomenon is not observed in other magnesium liquids, which maintain a rather constant coordination number. 
    more » « less
  3. ; ; ; (, Molecular Systems Design & Engineering)
    In semicrystalline conjugated polymer thin films, the mobility of charges depends on the arrangement of the individual polymer chains. In particular, the ordering of the polymer backbones affects the charge transport within the film, as electron transfer generally occurs along the backbones with alternating single and double bonds. In this paper, we demonstrate that polymer ordering should be discussed not only in terms of structural but also energetic ordering of polymer chains. We couple data from molecular dynamics simulations and quantum chemical calculations to quantify both structural and energetic ordering of polymer chains. We leverage a graph-based representation of the polymer chains to quantify the transport pathways in a computationally efficient way. Next, we formulate the morphological descriptors that correlate well with hole mobility determined using kinetic Monte Carlo simulations. We show that the shortest and fastest path calculations are predictive of mobility in equilibrated morphologies. In this sense, we leverage graph-based descriptors to provide a basis for the quantitative structure-property relationships. 
    more » « less
  4. ; (, The Journal of Chemical Physics)
    The Boyle temperature, TB, for an n-segment polymer in solution is the temperature where the second osmotic virial coefficient, A2, is zero. This characteristic is of interest for its connection to the polymer condensation critical temperature, particularly for n → ∞. TB can be measured experimentally or computed for a given model macromolecule. For the latter, we present and examine two approaches, both based on the Mayer-sampling Monte Carlo (MSMC) method, to calculate Boyle temperatures as a function of model parameters. In one approach, we use MSMC calculations to search for TB, as guided by the evaluation of temperature derivatives of A2. The second approach involves numerical integration of an ordinary differential equation describing how TB varies with a model parameter, starting from a known TB. Unlike general MSMC calculations, these adaptations are appealing because they neither invoke a reference for the calculation nor use special averages needed to avoid bias when computing A2 directly. We demonstrate these methods by computing TB lines for off-lattice linear Lennard-Jones polymers as a function of chain stiffness, considering chains of length n ranging from 2 to 512 monomers. We additionally perform calculations of single-molecule radius of gyration Rg and determine the temperatures Tθ, where linear scaling of Rg2 with n is observed, as if the polymers were long random-walk chains. We find that Tθ and TB seem to differ by 6% in the n → ∞ limit, which is beyond the statistical uncertainties of our computational methodology. However, we cannot rule out systematic error relating to our extrapolation procedure as being the source of this discrepancy. 
    more » « less
  5. ; ; (, Soft Matter)
    We propose an interpolation formula (the EV-WLC relation) for the force–extension behavior of wormlike chains in the presence of hard-core excluded volume interactions, analogous to the classic interpolation formula from Marko and Siggia for ideal wormlike chains. Using pruned-enriched Rosenbluth method (PERM) simulations of asymptotically long, discrete wormlike chains in an external force, we show that the error in the EV-WLC interpolation formula to describe discrete wormlike chains is systematically smaller than the error in the Marko–Siggia interpolation formula, except for the saturation region in which both formulas have the same limiting behavior. We anticipate that the EV-WLC interpolation formula will prove useful in the coarse-graining of wormlike chain models for dynamic simulations. Related results for the excess free energy due to excluded volume provide strong support for the physical basis of the Pincus regime. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email