-
Dynamics and Rheology of Polymer Melts via Hierarchical Atomistic, Coarse-grained, and Slip-spring Simulations
Authors:
Alireza F. Behbahani,
Ludwig Schneider,
Anastassia Rissanou,
Anthony Chazirakis,
Petra Bačová,
Pritam Kumar Jana,
Wei Li,
Manolis Doxastakis,
Patrycja Polińska,
Craig Burkhart,
Marcus Müller,
Vagelis A. Harmandaris
Abstract:
A hierarchical (triple scale) simulation methodology is presented for the prediction of the dynamical and rheological properties of high molecular weight entangled polymer melts. The methodology consists of atomistic, moderately coarse-grained (mCG), and highly coarse-grained slip-spring (SLSP) simulations. At the mCG level, a few chemically bonded atoms are lumped into one coarse-grained bead. At…
▽ More
A hierarchical (triple scale) simulation methodology is presented for the prediction of the dynamical and rheological properties of high molecular weight entangled polymer melts. The methodology consists of atomistic, moderately coarse-grained (mCG), and highly coarse-grained slip-spring (SLSP) simulations. At the mCG level, a few chemically bonded atoms are lumped into one coarse-grained bead. At this level, the chemical identity of the atomistic system, and the interchain topological constraints (entanglements) are preserved. The mCG potentials are derived by matching local structural distributions of the mCG model to those of the atomistic model through iterative Boltzmann inversion. For matching mCG and atomistic dynamics, the mCG time is scaled by a time scaling factor, which compensates for the lower monomeric friction coefficient of the mCG model than that of the atomistic one. At the SLSP level, multiple Kuhn segments of a polymer chain are represented by one coarse-grained bead. The very soft nonbonded interactions between beads do not prevent chain crossing and, hence, can not capture entanglements. The topological constraints are represented by slip-springs. A compensating pair potential is used in the SLSP model, to keep the static macromolecular properties unaltered upon the introduction of slip-springs. The static and kinetic parameters of the SLSP model are determined based on the lower level simulation models. Particularly, matching the orientational autocorrelation of the end-to-end vector, we determine the number of slip-springs and calibrate the timescale of the SLSP model. As the test case, the hierarchical methodology is applied to $cis$-1,4-polybutadiene (cPB) at 413 K. Dynamical properties of cPB melts are calculated for a broad range of molecular weights. The calculations are compared, and found in good agreement, with experimental data from the literature.
△ Less
Submitted 22 January, 2021;
originally announced January 2021.
-
Hyperbranched polymer stars with Gaussian chain statistics revisited
Authors:
P. Polinska,
C. Gillig,
J. P. Wittmer,
J. Baschnagel
Abstract:
Conformational properties of regular dendrimers and more general hyperbranched polymer stars with Gaussian statistics for the spacer chains between branching points are revisited numerically. We investigate the scaling for asymptotically long chains especially for fractal dimensions $d_f = 3$ (marginally compact) and $d_f = 2.5$ (diffusion limited aggregation). Power-law stars obtained by imposing…
▽ More
Conformational properties of regular dendrimers and more general hyperbranched polymer stars with Gaussian statistics for the spacer chains between branching points are revisited numerically. We investigate the scaling for asymptotically long chains especially for fractal dimensions $d_f = 3$ (marginally compact) and $d_f = 2.5$ (diffusion limited aggregation). Power-law stars obtained by imposing the number of additional arms per generation are compared to truly self-similar stars. We discuss effects of weak excluded volume interactions and sketch the regime where the Gaussian approximation should hold in dense solutions and melts for sufficiently large spacer chains.
△ Less
Submitted 15 August, 2015;
originally announced August 2015.
-
Pressure fluctuations in isotropic solids and fluids
Authors:
J. P. Wittmer,
H. Xu,
P. Polińska,
F. Weysser,
J. Baschnagel
Abstract:
Comparing isotropic solids and fluids at either imposed volume or pressure we investigate various correlations of the instantaneous pressure and its ideal and excess contributions. Focusing on the compression modulus K it is emphasized that the stress fluctuation representation of the elastic moduli may be obtained directly (without a microscopic displacement field) by comparing the stress fluctua…
▽ More
Comparing isotropic solids and fluids at either imposed volume or pressure we investigate various correlations of the instantaneous pressure and its ideal and excess contributions. Focusing on the compression modulus K it is emphasized that the stress fluctuation representation of the elastic moduli may be obtained directly (without a microscopic displacement field) by comparing the stress fluctuations in conjugated ensembles. This is made manifest by computing the Rowlinson stress fluctuation expression K_row of the compression modulus for NPT-ensembles. It is shown theoretically and numerically that K_row|P = P_id (2 - P_id/K) with P_id being the ideal pressure contribution.
△ Less
Submitted 7 May, 2013;
originally announced May 2013.
-
Shear modulus of simulated glass-forming model systems: Effects of boundary condition, temperature and sampling time
Authors:
J. P. Wittmer,
H. Xu,
P. Polińska,
F. Weysser,
J. Baschnagel
Abstract:
The shear modulus G of two glass-forming colloidal model systems in d=3 and d=2 dimensions is investigated by means of, respectively, molecular dynamics and Monte Carlo simulations. Comparing ensembles where either the shear strain gamma or the conjugated (mean) shear stress tau are imposed, we compute G from the respective stress and strain fluctuations as a function of temperature T while keepin…
▽ More
The shear modulus G of two glass-forming colloidal model systems in d=3 and d=2 dimensions is investigated by means of, respectively, molecular dynamics and Monte Carlo simulations. Comparing ensembles where either the shear strain gamma or the conjugated (mean) shear stress tau are imposed, we compute G from the respective stress and strain fluctuations as a function of temperature T while keeping a constant normal pressure P. The choice of the ensemble is seen to be highly relevant for the shear stress fluctuations mu_F(T) which at constant tau decay monotonously with T following the affine shear elasticity mu_A(T), i.e. a simple two-point correlation function. At variance, non-monotonous behavior with a maximum at the glass transition temperature T_g is demonstrated for mu_F(T) at constant gamma. The increase of G below T_g is reasonably fitted for both models by a continuous cusp singularity, G(T) is proportional to (1-T/T_g)^(1/2), in qualitative agreement with some recent replica calculations. It is argued, however, that longer sampling times may lead to a sharper transition. The additive jump discontinuity predicted by mode-coupling theory and other replica calculations thus cannot ultimately be ruled out.
△ Less
Submitted 7 December, 2012;
originally announced December 2012.
-
Impulsive correction to the elastic moduli obtained using the stress-fluctuation formalism in systems with truncated pair potential
Authors:
H. Xu,
J. P. Wittmer,
P. Polińska,
J. Baschnagel
Abstract:
The truncation of a pair potential at a distance r_cut is well-known to imply in general an impulsive correction to the pressure and other moments of the first derivatives of the potential. That depending on r_cut the truncation may also be of relevance to higher derivatives is shown theoretically for the Born contributions to the elastic moduli obtained using the stress-fluctuation formalism in d…
▽ More
The truncation of a pair potential at a distance r_cut is well-known to imply in general an impulsive correction to the pressure and other moments of the first derivatives of the potential. That depending on r_cut the truncation may also be of relevance to higher derivatives is shown theoretically for the Born contributions to the elastic moduli obtained using the stress-fluctuation formalism in d dimensions. Focusing on isotropic liquids for which the shear modulus G must vanish by construction, the predicted corrections are tested numerically for binary mixtures and polydisperse Lennard-Jones beads in, respectively, d=3 and d=2 dimensions.
△ Less
Submitted 2 August, 2012;
originally announced August 2012.
-
Strictly two-dimensional self-avoiding walks: Thermodynamic properties revisited
Authors:
N. Schulmann,
H. Meyer,
P. Polińska,
J. Baschnagel,
J. P. Wittmer
Abstract:
The density crossover scaling of various thermodynamic properties of solutions and melts of self-avoiding and highly flexible polymer chains without chain intersections confined to strictly two dimensions is investigated by means of molecular dynamics and Monte Carlo simulations of a standard coarse-grained bead-spring model. In the semidilute regime we confirm over an order of magnitude of the mo…
▽ More
The density crossover scaling of various thermodynamic properties of solutions and melts of self-avoiding and highly flexible polymer chains without chain intersections confined to strictly two dimensions is investigated by means of molecular dynamics and Monte Carlo simulations of a standard coarse-grained bead-spring model. In the semidilute regime we confirm over an order of magnitude of the monomer density rho the expected power-law scaling for the interaction energy between different chains e_inter\simρ^(21/8), the total pressure P\simρ^3 and the dimensionless compressibility gT=lim(q->0)(S(q)\sim1/ρ^2). Various elastic contributions associated to the affine and non-affine response to an infinitesimal strain are analyzed as functions of density and sampling time. We show how the size xi(rho) of the semidilute blob may be determined experimentally from the total monomer structure factor S(q) characterizing the compressibility of the solution at a given wavevector q. We comment briefly on finite persistence length effects.
△ Less
Submitted 31 August, 2012; v1 submitted 28 June, 2012;
originally announced June 2012.
-
Note: Scale-free center-of-mass displacement correlations in polymer films without topological constraints and momentum conservation
Authors:
J. P. Wittmer,
N. Schulmann,
P. Polińska,
J. Baschnagel
Abstract:
We present here computational work on the center-of-mass displacements in thin polymer films of finite width without topological constraints and without momentum conservation obtained using a well-known lattice Monte Carlo algorithm with chain lengths ranging up to N=8192. Computing directly the center-of-mass displacement correlation function C_N(t) allows to make manifest the existence of scale-…
▽ More
We present here computational work on the center-of-mass displacements in thin polymer films of finite width without topological constraints and without momentum conservation obtained using a well-known lattice Monte Carlo algorithm with chain lengths ranging up to N=8192. Computing directly the center-of-mass displacement correlation function C_N(t) allows to make manifest the existence of scale-free colored forces acting on a reference chain. As suggested by the scaling arguments put forward in a recent work on three-dimensional melts, we obtain a negative algebraic decay C_N(t) \sim -1/(Nt) for times t << T_N with T_N being the chain relaxation time. This implies a logarithmic correction to the related center-of-mass mean square-displacement h_N(t) as has been checked directly.
△ Less
Submitted 31 October, 2011;
originally announced October 2011.
-
Scale-free static and dynamical correlations in melts of monodisperse and Flory-distributed homopolymers: A review of recent bond-fluctuation model studies
Authors:
J. P. Wittmer,
A. Cavallo,
H. Xu,
J. E. Zabel,
P. Polińska,
N. Schulmann,
H. Meyer,
J. Farago,
A. Johner,
S. P. Obukhov,
J. Baschnagel
Abstract:
It has been assumed until very recently that all long-range correlations are screened in three-dimensional melts of linear homopolymers on distances beyond the correlation length $ξ$ characterizing the decay of the density fluctuations. Summarizing simulation results obtained by means of a variant of the bond-fluctuation model with finite monomer excluded volume interactions and topology violating…
▽ More
It has been assumed until very recently that all long-range correlations are screened in three-dimensional melts of linear homopolymers on distances beyond the correlation length $ξ$ characterizing the decay of the density fluctuations. Summarizing simulation results obtained by means of a variant of the bond-fluctuation model with finite monomer excluded volume interactions and topology violating local and global Monte Carlo moves, we show that due to an interplay of the chain connectivity and the incompressibility constraint, both static and dynamical correlations arise on distances $r \gg ξ$. These correlations are scale-free and, surprisingly, do not depend explicitly on the compressibility of the solution. Both monodisperse and (essentially) Flory-distributed equilibrium polymers are considered.
△ Less
Submitted 22 July, 2011;
originally announced July 2011.
-
Scale-free center-of-mass displacement correlations in dense polymer solutions and melts without topological constraints and momentum conservation: A bond-fluctuation model study
Authors:
Joachim P. Wittmer,
Patrycia Polińska,
Anna Cavallo,
Hendrik Meyer,
Jean Farago,
Albert Johner,
Jörg Baschnagel
Abstract:
By Monte Carlo simulations of a variant of the bond-fluctuation model without topological constraints we examine the center-of-mass (COM) dynamics of polymer melts in $d=3$ dimensions. Our analysis focuses on the COM displacement correlation function $\CN(t) \approx \partial_t^2 \MSDcmN(t)/2$, measuring the curvature of the COM mean-square displacement $\MSDcmN(t)$. We demonstrate that…
▽ More
By Monte Carlo simulations of a variant of the bond-fluctuation model without topological constraints we examine the center-of-mass (COM) dynamics of polymer melts in $d=3$ dimensions. Our analysis focuses on the COM displacement correlation function $\CN(t) \approx \partial_t^2 \MSDcmN(t)/2$, measuring the curvature of the COM mean-square displacement $\MSDcmN(t)$. We demonstrate that $\CN(t) \approx -(\RN/\TN)^2 (\rhostar/ρ) \ f(x=t/\TN)$ with $N$ being the chain length ($16 \le N \le 8192$), $\RN\sim N^{1/2}$ the typical chain size, $\TN\sim N^2$ the longest chain relaxation time, $ρ$ the monomer density, $\rhostar \approx N/\RN^d$ the self-density and $f(x)$ a universal function decaying asymptotically as $f(x) \sim x^{-ω}$ with $ω= (d+2) \times α$ where $α= 1/4$ for $x \ll 1$ and $α= 1/2$ for $x \gg 1$. We argue that the algebraic decay $N \CN(t) \sim - t^{-5/4}$ for $t \ll \TN$ results from an interplay of chain connectivity and melt incompressibility giving rise to the correlated motion of chains and subchains.
△ Less
Submitted 27 May, 2011; v1 submitted 24 February, 2011;
originally announced February 2011.