-
Strain correlation functions in isotropic elastic bodies: Large wavelength limit for two-dimensional systems
Authors:
J. P. Wittmer,
A. N. Semenov,
J. Baschnagel
Abstract:
Strain correlation functions in two-dimensional isotropic elastic bodies are shown both theoretically (using the general structure of isotropic tensor fields) and numerically (using a glass-forming model system) to depend on the coordinates of the field variable (position vector r in real space or wavevector q in reciprocal space) and thus on the direction of the field vector and the orientation o…
▽ More
Strain correlation functions in two-dimensional isotropic elastic bodies are shown both theoretically (using the general structure of isotropic tensor fields) and numerically (using a glass-forming model system) to depend on the coordinates of the field variable (position vector r in real space or wavevector q in reciprocal space) and thus on the direction of the field vector and the orientation of the coordinate system. Since the fluctuations of the longitudinal and transverse components of the strain field in reciprocal space are known in the long-wavelength limit from the equipartition theorem, all components of the correlation function tensor field are imposed and no additional physical assumptions are needed. An observed dependence on the field vector direction thus cannot be used as an indication for anisotropy or for a plastic rearrangement. This dependence is different for the associated strain response field containing also information on the localized stress perturbation
△ Less
Submitted 16 August, 2023; v1 submitted 29 March, 2023;
originally announced March 2023.
-
Correlations of tensor field components in isotropic systems with an application to stress correlations in elastic bodies
Authors:
J. P. Wittmer,
A. N. Semenov,
J. Baschnagel
Abstract:
Correlation functions of components of second-order tensor fields in isotropic systems can be reduced to an isotropic forth-order tensor field characterized by a few invariant correlation functions (ICFs). It is emphasized that components of this field depend in general on the coordinates of the field vector variable and thus on the orientation of the coordinate system. These angular dependencies…
▽ More
Correlation functions of components of second-order tensor fields in isotropic systems can be reduced to an isotropic forth-order tensor field characterized by a few invariant correlation functions (ICFs). It is emphasized that components of this field depend in general on the coordinates of the field vector variable and thus on the orientation of the coordinate system. These angular dependencies are distinct from those of ordinary anisotropic systems. As a simple example of the procedure to obtain the ICFs we discuss correlations of time-averaged stresses in isotropic glasses where only one ICF in reciprocal space becomes a finite constant e for large sampling times and small wavevectors. It is shown that e is set by the typical size of the frozen-in stress components normal to the wavevectors, i.e. it is caused by the symmetry breaking of the stress for each independent configuration. Using the presented general mathematical formalism for isotropic tensor fields this finding explains in turn the observed long-range stress correlations in real space. Under additional but rather general assumptions e is shown to be given by a thermodynamic quantity, the equilibrium Young modulus E. We thus relate for certain isotropic amorphous bodies the existence of finite Young or shear moduli to the symmetry breaking of a stress component in reciprocal space.
△ Less
Submitted 13 June, 2023; v1 submitted 29 March, 2023;
originally announced March 2023.
-
Spatial correlation functions for non-ergodic stochastic processes of macroscopic system
Authors:
J. P. Wittmer,
A. N. Semenov,
J. Baschnagel
Abstract:
Focusing on non-ergodic macroscopic systems we reconsider the variances of time averages time-series. The total variance (direct average over all time-series) is known to be the sum of an internal variance (fluctuations within the meta-basins) and an external variance (fluctuations between meta-basins). It is shown that whenever the time-averaged observable can be expressed as a volume average of…
▽ More
Focusing on non-ergodic macroscopic systems we reconsider the variances of time averages time-series. The total variance (direct average over all time-series) is known to be the sum of an internal variance (fluctuations within the meta-basins) and an external variance (fluctuations between meta-basins). It is shown that whenever the time-averaged observable can be expressed as a volume average of a local field the three variances can be written as volume averages of correlation functions with the total correlation function being the sum of an internal and an external correlation function. The dependences of the the different variancescan thus be traced back to the internal and the external correlation function. Various relations are illustrated using lattice spring models with spatially correlated spring constants.
△ Less
Submitted 18 August, 2022;
originally announced August 2022.
-
Numerical determination of shear stress relaxation modulus of polymer glasses
Authors:
I. Kriuchevskyi,
J. P. Wittmer,
O. Benzerara,
H. Meyer,
J. Baschnagel
Abstract:
Focusing on simulated polymer glasses well below the glass transition, we confirm the validity and the efficiency of the recently proposed simple-average expression $G(t) = μ_A - h(t)$ for the computational determination of the shear stress relaxation modulus $G(t)$. Here, $μ_A = G(0)$ characterizes the affine shear transformation of the system at $t=0$ and $h(t)$ the mean-square displacement of t…
▽ More
Focusing on simulated polymer glasses well below the glass transition, we confirm the validity and the efficiency of the recently proposed simple-average expression $G(t) = μ_A - h(t)$ for the computational determination of the shear stress relaxation modulus $G(t)$. Here, $μ_A = G(0)$ characterizes the affine shear transformation of the system at $t=0$ and $h(t)$ the mean-square displacement of the instantaneous shear stress as a function of time $t$. This relation is seen to be particulary useful for systems with quenched or sluggish transient shear stresses which necessarily arise below the glass transition. The commonly accepted relation $G(t)=c(t)$ using the shear stress auto-correlation function $c(t)$ becomes incorrect in this limit.
△ Less
Submitted 31 January, 2022;
originally announced January 2022.
-
Marginally compact hyperbranched polymer trees
Authors:
M. Dolgushev,
J. P. Wittmer,
A. Johner,
O. Benzerara,
H. Meyer,
J. Baschnagel
Abstract:
Assuming Gaussian chain statistics along the chain contour, we generate by means of a proper fractal generator hyperbranched polymer trees which are marginally compact. Static and dynamical properties, such as the radial intrachain pair density distribution or the shear-stress relaxation modulus, are investigated theoretically and by means of computer simulations. We emphasize that albeit the self…
▽ More
Assuming Gaussian chain statistics along the chain contour, we generate by means of a proper fractal generator hyperbranched polymer trees which are marginally compact. Static and dynamical properties, such as the radial intrachain pair density distribution or the shear-stress relaxation modulus, are investigated theoretically and by means of computer simulations. We emphasize that albeit the self-contact density diverges logarithmically with the total mass $N$, this effect becomes rapidly irrelevant with increasing spacer length $S$. In addition to this it is seen that the standard Rouse analysis must necessarily become inappropriate for compact objects for which the relaxation time $τ_p$ of mode $p$ must scale as $τ_p \sim (N/p)^{5/3}$ rather than the usual square power law for linear chains.
△ Less
Submitted 31 January, 2022;
originally announced January 2022.
-
Simple models for strictly non-ergodic stochastic processes of macroscopic systems
Authors:
G. George,
L. Klochko,
A. N. Semenov,
J. Baschnagel,
J. P. Wittmer
Abstract:
We investigate simple models for strictly non-ergodic stochastic processes $x_t$ ($t$ being the discrete time step) focusing on the expectation value $v$ and the standard deviation $δv$ of the empirical variance $v[x]$ of finite time series $x$. $x_t$ is averaged over a fluctuating field $σ_{r}$ ($r$ being the microcell position) characterized by a quenched spatially correlated Gaussian field. Due…
▽ More
We investigate simple models for strictly non-ergodic stochastic processes $x_t$ ($t$ being the discrete time step) focusing on the expectation value $v$ and the standard deviation $δv$ of the empirical variance $v[x]$ of finite time series $x$. $x_t$ is averaged over a fluctuating field $σ_{r}$ ($r$ being the microcell position) characterized by a quenched spatially correlated Gaussian field. Due to the quenched field $δv(Δt)$ becomes a finite constant, $Δ_{ne} > 0$, for large sampling times $Δt$. The volume dependence of the non-ergodicity parameter $Δ_{ne}$ is investigated for different spatial correlations. Models with marginally long-ranged $\fr$-correlations are successfully mapped on shear-stress data from simulated amorphous glasses of polydisperse beads.
△ Less
Submitted 22 November, 2021;
originally announced November 2021.
-
Fluctuations of non-ergodic stochastic processes
Authors:
G. George,
L. Klochko,
A. N. Semenov,
J. Baschnagel,
J. P. Wittmer
Abstract:
We investigate the standard deviation $δv(\tsamp)$ of the variance $v[\xbf]$ of time series $\xbf$ measured over a finite sampling time $\tsamp$ focusing on non-ergodic systems where independent "configurations" $c$ get trapped in meta-basins of a generalized phase space. It is thus relevant in which order averages over the configurations $c$ and over time series $k$ of a configuration $c$ are per…
▽ More
We investigate the standard deviation $δv(\tsamp)$ of the variance $v[\xbf]$ of time series $\xbf$ measured over a finite sampling time $\tsamp$ focusing on non-ergodic systems where independent "configurations" $c$ get trapped in meta-basins of a generalized phase space. It is thus relevant in which order averages over the configurations $c$ and over time series $k$ of a configuration $c$ are performed. Three variances of $v[\xbf_{ck}]$ must be distinguished: the total variance $\dvtot = \dvint + \dvext$ and its contributions $\dvint$, the typical internal variance within the meta-basins, and $\dvext$, characterizing the dispersion between the different basins. We discuss simplifications for physical systems where the stochastic variable $x(t)$ is due to a density field averaged over a large system volume $V$. The relations are illustrated for the shear-stress fluctuations in quenched elastic networks and low-temperature glasses formed by polydisperse particles and free-standing polymer films. The different statistics of $\svint$ and $\svext$ are manifested by their different system-size dependence
△ Less
Submitted 12 March, 2021;
originally announced March 2021.
-
Ensemble fluctuations matter for variances of macroscopic variables
Authors:
G. George,
L. Klochko,
A. N. Semenov,
J. Baschnagel,
J. P. Wittmer
Abstract:
Extending recent work on stress fluctuations in complex fluids and amorphous solids we describe in general terms the ensemble average $v(Δt)$ and the standard deviation $δv(Δt)$ of the variance $v[\mathbf{x}]$ of time series $\mathbf{x}$ of a stochastic process $x(t)$ measured over a finite sampling time $Δt$. Assuming a stationary, Gaussian and ergodic process, $δv$ is given by a functional…
▽ More
Extending recent work on stress fluctuations in complex fluids and amorphous solids we describe in general terms the ensemble average $v(Δt)$ and the standard deviation $δv(Δt)$ of the variance $v[\mathbf{x}]$ of time series $\mathbf{x}$ of a stochastic process $x(t)$ measured over a finite sampling time $Δt$. Assuming a stationary, Gaussian and ergodic process, $δv$ is given by a functional $δv_G[h]$ of the autocorrelation function $h(t)$. $δv(Δt)$ is shown to become large and similar to $v(Δt)$ if $Δt$ corresponds to a fast relaxation process. Albeit $δv = δv_G[h]$ does not hold in general for non-ergodic systems, the deviations for common systems with many microstates are merely finite-size corrections. Various issues are illustrated for shear-stress fluctuations in simple coarse-grained model systems.
△ Less
Submitted 17 November, 2020;
originally announced November 2020.
-
Glassy dynamics of a binary Voronoi fluid: A mode-coupling analysis
Authors:
Céline Ruscher,
Simone Ciarella,
Chengjie Luo,
Liesbeth M. C. Janssen,
Jean Farago,
Jörg Baschnagel
Abstract:
The binary Voronoi mixture is a fluid model whose interactions are local and many-body. Here we perform molecular-dynamics (MD) simulations of an equimolar mixture that is weakly polydisperse and additive. For the first time we study the structural relaxation of this mixture in the supercooled-liquid regime. From the simulations we determine the time- and temperature-dependent scattering functions…
▽ More
The binary Voronoi mixture is a fluid model whose interactions are local and many-body. Here we perform molecular-dynamics (MD) simulations of an equimolar mixture that is weakly polydisperse and additive. For the first time we study the structural relaxation of this mixture in the supercooled-liquid regime. From the simulations we determine the time- and temperature-dependent scattering functions for a large range of wave vectors, as well as the mean-square displacements of both particle species. We perform a detailed analysis of the dynamics by comparing the MD results with the first-principles-based idealized mode-coupling theory (MCT). To this end, we employ two approaches: fits to the asymptotic predictions of the theory, and fit-parameter-free binary MCT calculations based on static-structure-factor input from the simulations. We find that many-body interactions of the Voronoi mixture do not lead to strong qualitative differences relative to similar analyses carried out for simple liquids with pair-wise interactions. For instance, the fits give an exponent parameter $λ\approx 0.746$ comparable to typical values found for simple liquids, the wavevector dependence of the Kohlrausch relaxation time is in good agreement with literature results for polydisperse hard spheres, and the MCT calculations based on static input overestimate the critical temperature, albeit only by a factor of about 1.2. This overestimation appears to be weak relative to other well-studied supercooled-liquid models such as the binary Kob--Andersen Lennard-Jones mixture. Overall, the agreement between MCT and simulation suggests that it is possible to predict several microscopic dynamic properties with qualitative, and in some cases near-quantitative, accuracy based solely on static two-point structural correlations, even though the system itself is inherently governed by many-body interactions.
△ Less
Submitted 10 August, 2020;
originally announced August 2020.
-
Shear-stress fluctuations in free-standing polymer films
Authors:
G. George,
I. Kriuchevskyi,
H. Meyer,
J. Baschnagel,
J. P. Wittmer
Abstract:
Using molecular dynamics simulation of a polymer glass model we investigate free-standing polymer films focusing on the in-plane shear modulus $μ$ and the corresponding shear-stress relaxation modulus $G(t)$ as functions of temperature $T$, film thickness $H$ (tuned by means of the lateral box size $L$) and sampling time $Δt$. Various observables are seen to vary linearly with $1/H$ demonstrating…
▽ More
Using molecular dynamics simulation of a polymer glass model we investigate free-standing polymer films focusing on the in-plane shear modulus $μ$ and the corresponding shear-stress relaxation modulus $G(t)$ as functions of temperature $T$, film thickness $H$ (tuned by means of the lateral box size $L$) and sampling time $Δt$. Various observables are seen to vary linearly with $1/H$ demonstrating thus the (to leading order) linear superposition of bulk and surface properties. In agreement with recent studies on three-dimensional polymer glass-formers, $μ$ and $G(t)$ are found to decrease continuously with $T$. A jump-singularity is not observed. Confirming the time-translational invariance of our systems, the $Δt$-dependence of $μ$ is traced back to $G(t)$.
△ Less
Submitted 28 November, 2018; v1 submitted 7 September, 2018;
originally announced September 2018.
-
Voronoi Glass-Forming Liquids : A Structural Study
Authors:
Céline Ruscher,
Joerg Baschnagel,
Jean Farago
Abstract:
We introduce a new theoretical model of simple fluid, whose interactions, defined in terms of the Voronoi cells of the configurations, are local and many-body. The resulting system is studied both theoretically and numerically. We show that the fluid, though sharing the global features of other models of fluids with soft interactions, has several unusual characteristics, which are investigated and…
▽ More
We introduce a new theoretical model of simple fluid, whose interactions, defined in terms of the Voronoi cells of the configurations, are local and many-body. The resulting system is studied both theoretically and numerically. We show that the fluid, though sharing the global features of other models of fluids with soft interactions, has several unusual characteristics, which are investigated and discussed.
△ Less
Submitted 23 December, 2017;
originally announced December 2017.
-
Shear modulus and shear-stress fluctuations in polymer glasses
Authors:
I. Kriuchevskyi,
J. P. Wittmer,
H. Meyer,
J. Baschnagel
Abstract:
Using molecular dynamics simulation of a standard coarse-grained polymer glass model we investigate by means of the stress-fluctuation formalism the shear modulus $μ$ as a function of temperature $T$ and sampling time $Δt$. While the ensemble-averaged modulus $μ(T)$ is found to decrease continuously for all $Δt$ sampled, its standard deviation $δμ(T)$ is non-monotonous with a striking peak at the…
▽ More
Using molecular dynamics simulation of a standard coarse-grained polymer glass model we investigate by means of the stress-fluctuation formalism the shear modulus $μ$ as a function of temperature $T$ and sampling time $Δt$. While the ensemble-averaged modulus $μ(T)$ is found to decrease continuously for all $Δt$ sampled, its standard deviation $δμ(T)$ is non-monotonous with a striking peak at the glass transition. Confirming the effective time-translational invariance of our systems, $μ(Δt)$ can be understood using a weighted integral over the shear-stress relaxation modulus $G(t)$. While the crossover of $μ(T)$ gets sharper with increasing $Δt$, the peak of $δμ(T)$ becomes more singular. %
It is thus elusive to predict the modulus of a single configuration at the glass transition.
△ Less
Submitted 2 November, 2017;
originally announced November 2017.
-
Shear-stress fluctuations and relaxation in polymer glasses
Authors:
I. Kriuchevskyi,
J. P. Wittmer,
H. Meyer,
O. Benzerara,
J. Baschnagel
Abstract:
We investigate by means of molecular dynamics simulation a coarse-grained polymer glass model focusing on (quasi-static and dynamical) shear-stress fluctuations as a function of temperature T and sampling time $Δt$. The linear response is characterized using (ensemble-averaged) expectation values of the contributions (time-averaged for each shear plane) to the stress-fluctuation relation $μ_{sf}$…
▽ More
We investigate by means of molecular dynamics simulation a coarse-grained polymer glass model focusing on (quasi-static and dynamical) shear-stress fluctuations as a function of temperature T and sampling time $Δt$. The linear response is characterized using (ensemble-averaged) expectation values of the contributions (time-averaged for each shear plane) to the stress-fluctuation relation $μ_{sf}$ for the shear modulus and the shear-stress relaxation modulus $G(t)$. Using 100 independent configurations we pay attention to the respective standard deviations. While the ensemble-averaged modulus $μ_{sf}(T)$ decreases continuously with increasing T for all $Δt$ sampled, its standard deviation $δμ_{sf}(T)$ is non-monotonous with a striking peak at the glass transition. The question of whether the shear modulus is continuous or has a jump-singularity at the glass transition is thus ill-posed. Confirming the effective time-translational invariance of our systems, the $Δt$-dependence of $μ_{sf}$ and related quantities can be understood using a weighted integral over $G(t)$. This implies that the shear viscosity $η(T)$ may be readily obtained from the $1/Δt$-decay of $μ_{sf}$ above the glass transition.
△ Less
Submitted 2 November, 2017;
originally announced November 2017.
-
arXiv:1611.04860
[pdf]
cond-mat.soft
cond-mat.mes-hall
cond-mat.stat-mech
physics.comp-ph
physics.flu-dyn
Molecular dynamics simulation of the capillary leveling of viscoelastic polymer films
Authors:
Ioannis Tanis,
Hendrik Meyer,
Thomas Salez,
Elie Raphaël,
Anthony C. Maggs,
Jörg Baschnagel
Abstract:
Surface tension-driven flow techniques have recently emerged as an efficient means of shedding light into the rheology of thin polymer films. Motivated by experimental and theoretical approaches in films bearing a varying surface topography, we present results on the capillary relaxation of a square pattern at the free surface of a viscoelastic polymer film, using molecular dynamics simulations of…
▽ More
Surface tension-driven flow techniques have recently emerged as an efficient means of shedding light into the rheology of thin polymer films. Motivated by experimental and theoretical approaches in films bearing a varying surface topography, we present results on the capillary relaxation of a square pattern at the free surface of a viscoelastic polymer film, using molecular dynamics simulations of a coarse-grained polymer model. Height profiles are monitored as a function of time after heating the system above its glass-transition temperature and their time dependence is fitted to the theory of capillary leveling. Results show that the viscosity is not constant, but time dependent. In addition to providing a complementary insight about the local inner mechanisms, our simulations of the capillary-leveling process therefore probe the viscoelasticity of the polymer and not only its viscosity, in contrast to previous experimental approaches.
△ Less
Submitted 20 February, 2017; v1 submitted 15 November, 2016;
originally announced November 2016.
-
Shear-stress fluctuations in self-assembled transient elastic networks
Authors:
J. P. Wittmer,
I. Kriuchevskyi,
A. Cavallo,
H. Xu,
J. Baschnagel
Abstract:
Focusing on shear-stress fluctuations we investigate numerically a simple generic model for self-assembled transient networks formed by repulsive beads reversibly bridged by ideal springs. With $Δdt$ being the sampling time and $t_*(f) \sim 1/f$ the Maxwell relaxation time (set by the spring recombination frequency $f$) the dimensionless parameter $Δx = dt/t_*(f)$ is systematically scanned from th…
▽ More
Focusing on shear-stress fluctuations we investigate numerically a simple generic model for self-assembled transient networks formed by repulsive beads reversibly bridged by ideal springs. With $Δdt$ being the sampling time and $t_*(f) \sim 1/f$ the Maxwell relaxation time (set by the spring recombination frequency $f$) the dimensionless parameter $Δx = dt/t_*(f)$ is systematically scanned from the liquid limit ($Δdx \gg 1)$ to the solid limit ($Δx \ll 1$) where the network topology is quenched and an ensemble average over $m$ independent configurations is required. Generalizing previous work on permanent networks it is shown that the shear-stress relaxation modulus $G(t)$ may be efficiently determined for all $Δx$ using the simple-average expression $G(t) = μ_A - h(t)$ with $μ_A = G(0)$ characterizing the canonical-affine shear transformation of the system at $t=0$ and $h(t)$ the (rescaled) mean-square displacement of the instantaneous shear stress as a function of time $t$. This relation is compared to the standard expression $G(t) = C(t)$ using the (rescaled) shear-stress autocorrelation function $C(t)$. Lower bounds for the $m$ configurations required by both relations are given.
△ Less
Submitted 21 June, 2016;
originally announced June 2016.
-
Simple-average expressions for shear-stress relaxation modulus
Authors:
J. P. Wittmer,
H. Xu,
J. Baschnagel
Abstract:
Focusing on isotropic elastic networks we propose a novel simple-average expression $G(t) = μ_A - h(t)$ for the computational determination of the shear-stress relaxation modulus $G(t)$ of a classical elastic solid or fluid and its equilibrium modulus $\G_{eq} = \lim_{t \to \infty} G(t)$. Here, $μ_A = G(0)$ characterizes the shear transformation of the system at $t=0$ and $h(t)$ the (rescaled) mea…
▽ More
Focusing on isotropic elastic networks we propose a novel simple-average expression $G(t) = μ_A - h(t)$ for the computational determination of the shear-stress relaxation modulus $G(t)$ of a classical elastic solid or fluid and its equilibrium modulus $\G_{eq} = \lim_{t \to \infty} G(t)$. Here, $μ_A = G(0)$ characterizes the shear transformation of the system at $t=0$ and $h(t)$ the (rescaled) mean-square displacement of the instantaneous shear stress $\hatτ(t)$ as a function of time $t$. While investigating sampling time effects we also discuss the related expressions in terms of shear-stress autocorrelation functions. We argue finally that our key relation may be readily adapted for more general linear response functions.
△ Less
Submitted 6 October, 2015;
originally announced October 2015.
-
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.
-
Compressibility and pressure correlations in isotropic solids and fluids
Authors:
J. P. Wittmer,
H. Xu,
P. Polinsak,
C. Gillig,
J. Helfferich,
F. Weysser,
J. Baschnagel
Abstract:
Presenting simple coarse-grained models of isotropic solids and fluids in $d=1$, $2$ and $3$ dimensions we investigate the correlations of the instantaneous pressure and its ideal and excess contributions at either imposed pressure (NPT-ensemble, $λ=0$) or volume (NVT-ensemble, $λ=1$) and for more general values of the dimensionless parameter $λ$ characterizing the constant-volume constraint.
Presenting simple coarse-grained models of isotropic solids and fluids in $d=1$, $2$ and $3$ dimensions we investigate the correlations of the instantaneous pressure and its ideal and excess contributions at either imposed pressure (NPT-ensemble, $λ=0$) or volume (NVT-ensemble, $λ=1$) and for more general values of the dimensionless parameter $λ$ characterizing the constant-volume constraint.
△ Less
Submitted 15 August, 2015;
originally announced August 2015.
-
Fluctuation-dissipation relation between shear stress relaxation modulus and shear stress autocorrelation function revisited
Authors:
J. P. Wittmer,
H. Xu,
O. Benzerara,
J. Baschnagel
Abstract:
The shear stress relaxation modulus $G(t)$ may be determined from the shear stress $τ(t)$ after switching on a tiny step strain $γ$ or by inverse Fourier transformation of the storage modulus $G^{\prime}(ω)$ or the loss modulus $G^{\prime\prime}(ω)$ obtained in a standard oscillatory shear experiment at angular frequency $ω$. It is widely assumed that $G(t)$ is equivalent in general to the equilib…
▽ More
The shear stress relaxation modulus $G(t)$ may be determined from the shear stress $τ(t)$ after switching on a tiny step strain $γ$ or by inverse Fourier transformation of the storage modulus $G^{\prime}(ω)$ or the loss modulus $G^{\prime\prime}(ω)$ obtained in a standard oscillatory shear experiment at angular frequency $ω$. It is widely assumed that $G(t)$ is equivalent in general to the equilibrium stress autocorrelation function $C(t) = βV \langle δτ(t) δτ(0)\rangle$ which may be readily computed in computer simulations ($β$ being the inverse temperature and $V$ the volume). Focusing on isotropic solids formed by permanent spring networks we show theoretically by means of the fluctuation-dissipation theorem and computationally by molecular dynamics simulation that in general $G(t) = G_{eq} + C(t)$ for $t > 0$ with $G_{eq}$ being the static equilibrium shear modulus. A similar relation holds for $G^{\prime}(ω)$. $G(t)$ and $C(t)$ must thus become different for a solid body and it is impossible to obtain $G_{eq}$ directly from $C(t)$.
△ Less
Submitted 15 August, 2015;
originally announced August 2015.
-
Shear stress relaxation and ensemble transformation of shear stress autocorrelation functions revisited
Authors:
J. P. Wittmer,
H. Xu,
J. Baschnagel
Abstract:
We revisit the relation between the shear stress relaxation modulus $G(t)$, computed at finite shear strain $0 < γ\ll 1$, and the shear stress autocorrelation functions $C(t)|_γ$ and $C(t)|_τ$ computed, respectively, at imposed strain $γ$ and mean stress $τ$. Focusing on permanent isotropic spring networks it is shown theoretically and computationally that in general…
▽ More
We revisit the relation between the shear stress relaxation modulus $G(t)$, computed at finite shear strain $0 < γ\ll 1$, and the shear stress autocorrelation functions $C(t)|_γ$ and $C(t)|_τ$ computed, respectively, at imposed strain $γ$ and mean stress $τ$. Focusing on permanent isotropic spring networks it is shown theoretically and computationally that in general $G(t) = C(t)|_τ = C(t)|_γ + G_{eq}$ for $t > 0$ with $G_{eq}$ being the static equilibrium shear modulus. $G(t)$ and $C(t)|_γ$ thus must become different for solids and it is impossible to obtain $G_{eq}$ alone from $C(t)|_γ$ as often assumed. We comment briefly on self-assembled transient networks where $G_{eq}(f)$ must vanish for a finite scission-recombination frequency $f$. We argue that $G(t) = C(t)|_τ = C(t)|_γ$ should reveal an intermediate plateau set by the shear modulus $G_{eq}(f=0)$ of the quenched network.
△ Less
Submitted 15 August, 2015;
originally announced August 2015.
-
Shear-strain and shear-stress fluctuations in generalized Gaussian ensemble simulations of isotropic elastic networks
Authors:
J. P. Wittmer,
I. Kriuchevskyi,
J. Baschnagel,
H. Xu
Abstract:
Shear-strain and shear-stress correlations in isotropic elastic bodies are investigated both theoretically and numerically at either imposed mean shear-stress $τ$ ($λ=0$) or shear-strain $γ$ ($λ=1$) and for more general values of a dimensionless parameter $λ$ characterizing the generalized Gaussian ensemble. It allows to tune the strain fluctuations $μ_{γγ} \equiv βV \la δγ^2 \ra = (1-λ)/G_{eq}$ w…
▽ More
Shear-strain and shear-stress correlations in isotropic elastic bodies are investigated both theoretically and numerically at either imposed mean shear-stress $τ$ ($λ=0$) or shear-strain $γ$ ($λ=1$) and for more general values of a dimensionless parameter $λ$ characterizing the generalized Gaussian ensemble. It allows to tune the strain fluctuations $μ_{γγ} \equiv βV \la δγ^2 \ra = (1-λ)/G_{eq}$ with $β$ being the inverse temperature, $V$ the volume, $γ$ the instantaneous strain and $G_{eq}$ the equilibrium shear modulus. Focusing on spring networks in two dimensions we show, e.g., for the stress fluctuations $μ_{ττ} \equiv βV \la δτ^2 \ra$ ($τ$ being the instantaneous stress) that $μ_{ττ} = μ_{A} - λG_{eq}$ with $μ_{A} = μ_{ττ}|_{λ=0}$ being the affine shear-elasticity. For the stress autocorrelation function $c_{ττ}(t) \equiv βV \la δτ(t) δτ(0) \ra$ this result is then seen (assuming a sufficiently slow shear-stress barostat) to generalize to $c_{ττ}(t) = G(t) - λ\Geq$ with $G(t)$ being the shear-stress relaxation modulus.
△ Less
Submitted 15 August, 2015;
originally announced August 2015.
-
New conserved structural fields for supercooled liquids
Authors:
Jean Farago,
Alexander N. Semenov,
Stefan Frey,
Joerg Baschnagel
Abstract:
By considering Voronoi tessellations of the configurations of a fluid, we propose two new conserved fields, which provide structural information not fully accounted for by the usual 2-point density field fluctuations (structure factor). One of these fields is scalar and associated to the Voronoi cell volumes, whereas the other one, termed the "geometrical polarisation", is vectorial, related to th…
▽ More
By considering Voronoi tessellations of the configurations of a fluid, we propose two new conserved fields, which provide structural information not fully accounted for by the usual 2-point density field fluctuations (structure factor). One of these fields is scalar and associated to the Voronoi cell volumes, whereas the other one, termed the "geometrical polarisation", is vectorial, related to the very local anisotropy of the configurations. We study the static and dynamical properties of these fields in the supercooled regime of a model glass-forming liquid. We show in particular that the geometrical polarisation is both statically correlated to the force field and contrary to it develops a plateau regime when the temperature is lowered. We attribute this behaviour to the microsopic disorder of the underlying inherent structures (IS) which dictate the dynamics on time scales larger than the true microscopic time, in the strong supercooled regime. In this respect, this work raises the issue of to what extent the inter IS dynamics, intrinsically anisotropic and collective (cf. T.B. Schrøder et al. {\it J. of Chem. Phys.}, {\bf 112}, 9834 (2000)), could be related to their polarisation field.
△ Less
Submitted 12 December, 2013;
originally announced December 2013.
-
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.
-
Interchain monomer contact probability in two-dimensional polymer solutions
Authors:
N. Schulmann,
H. Meyer,
J. P. Wittmer,
A. Johner,
J. Baschnagel
Abstract:
Using molecular dynamics simulation of a standard bead-spring model we investigate the density crossover scaling of strictly two-dimensional self-avoiding polymer chains focusing on properties related to the contact exponent set by the intrachain subchain size distribution. Irrespective of the density sufficiently long chains are found to consist of compact packings of blobs of fractal perimeter d…
▽ More
Using molecular dynamics simulation of a standard bead-spring model we investigate the density crossover scaling of strictly two-dimensional self-avoiding polymer chains focusing on properties related to the contact exponent set by the intrachain subchain size distribution. Irrespective of the density sufficiently long chains are found to consist of compact packings of blobs of fractal perimeter dimension dp = 5/4.
△ Less
Submitted 20 December, 2011;
originally announced December 2011.
-
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.
-
Static Properties of Polymer Melts in Two Dimensions
Authors:
H. Meyer,
J. P. Wittmer,
T. Kreer,
A. Johner,
J. Baschnagel
Abstract:
Self-avoiding polymers in strictly two-dimensional ($d=2$) melts are investigated by means of molecular dynamics simulation of a standard bead-spring model with chain lengths ranging up to N=2048. % The chains adopt compact configurations of typical size $R(N) \sim N^ν$ with $ν=1/d$. % The precise measurement of various distributions of internal chain distances allows a direct test of the contact…
▽ More
Self-avoiding polymers in strictly two-dimensional ($d=2$) melts are investigated by means of molecular dynamics simulation of a standard bead-spring model with chain lengths ranging up to N=2048. % The chains adopt compact configurations of typical size $R(N) \sim N^ν$ with $ν=1/d$. % The precise measurement of various distributions of internal chain distances allows a direct test of the contact exponents $Θ_0=3/8$, $Θ_1=1/2$ and $Θ_2=3/4$ predicted by Duplantier. % Due to the segregation of the chains the ratio of end-to-end distance $\Rend(N)$ and gyration radius $\Rgyr(N)$ becomes $\Rend^2(N)/\Rgyr^2(N) \approx 5.3 < 6$ for $N \gg 100$ and the chains are more spherical than Gaussian phantom chains. % The second Legendre polynomial $P_2(s)$ of the bond vectors decays as $P_2(s) \sim 1/s^{1+νΘ_2}$ measuring thus the return probability of the chain after $s$ steps. % The irregular chain contours are shown to be characterized by a perimeter length $L(N) \sim R(N)^{\dc}$ of fractal line dimension $\dc = d-Θ_2 =5/4$. % % In agreement with the generalized Porod scattering of compact objects with fractal contour the Kratky representation of the intramolecular structure factor $F(q)$ reveals a strong non-monotonous behavior with $q^dF(q) \sim 1/(q R(N))^{Θ_2}$ in the intermediate regime of the wave vector $q$. This may allow to confirm the predicted contour fractality in a real experiment.
△ Less
Submitted 23 April, 2010;
originally announced April 2010.
-
Molecular dynamics simulations of glassy polymers
Authors:
Jean-Louis Barrat,
Joerg Baschnagel,
Alexei Lyulin
Abstract:
We review recent results from computer simulation studies of polymer glasses, from chain dynamics around the glass transition temperature Tg to the mechanical behaviour below Tg. These results clearly show that modern computer simulations are able to address and give clear answers to some important issues in the field, in spite of the obvious limitations in terms of length and time scales. In th…
▽ More
We review recent results from computer simulation studies of polymer glasses, from chain dynamics around the glass transition temperature Tg to the mechanical behaviour below Tg. These results clearly show that modern computer simulations are able to address and give clear answers to some important issues in the field, in spite of the obvious limitations in terms of length and time scales. In the present review we discuss the cooling rate effects, and dynamic slowing down of different relaxation processes when approaching Tg for both model and chemistry-specific polymer glasses. The impact of geometric confinement on the glass transition is discussed in detail. We also show that computer simulations are very useful tools to study structure and mechanical response of glassy polymers. The influence of large deformations on mechanical behaviour of polymer glasses in general, and strain hardening effect in particular are reviewed. Finally, we suggest some directions for future research, which we believe will be soon within the capabilities of state of the art computer simulations, and correspond to problems of fundamental interest.
△ Less
Submitted 10 February, 2010;
originally announced February 2010.
-
Non-extensivity of the chemical potential of polymer melts
Authors:
J. P. Wittmer,
A. Johner,
A. Cavallo,
P. Beckrich,
F. Crevel,
J. Baschnagel
Abstract:
Following Flory's ideality hypothesis the chemical potential of a test chain of length $n$ immersed into a dense solution of chemically identical polymers of length distribution P(N) is extensive in $n$. We argue that an additional contribution $δμ_c(n) \sim +1/ρ\sqrt{n}$ arises ($ρ$ being the monomer density) for all $¶(N)$ if $n \ll <N>$ which can be traced back to the overall incompressibilit…
▽ More
Following Flory's ideality hypothesis the chemical potential of a test chain of length $n$ immersed into a dense solution of chemically identical polymers of length distribution P(N) is extensive in $n$. We argue that an additional contribution $δμ_c(n) \sim +1/ρ\sqrt{n}$ arises ($ρ$ being the monomer density) for all $¶(N)$ if $n \ll <N>$ which can be traced back to the overall incompressibility of the solution leading to a long-range repulsion between monomers. Focusing on Flory distributed melts we obtain $δμ_c(n) \approx (1- 2 n/<N>) / ρ\sqrt{n}$ for $n \ll <N>^2$, hence, $δμ_c(n) \approx - 1/ρ\sqrt{n}$ if $n$ is similar to the typical length of the bath $<N>$. Similar results are obtained for monodisperse solutions. Our perturbation calculations are checked numerically by analyzing the annealed length distribution P(N) of linear equilibrium polymers generated by Monte Carlo simulation of the bond-fluctuation model. As predicted we find, e.g., the non-exponentiality parameter $K_p \equiv 1 - <N^>/p!<N>^p$ to decay as $K_p \approx 1 / \sqrt{<N>}$ for all moments $p$ of the distribution.
△ Less
Submitted 17 December, 2009;
originally announced December 2009.
-
Distance dependence of angular correlations in dense polymer solutions
Authors:
J. P. Wittmer,
A. Johner,
S. P. Obukhov,
H. Meyer,
A. Cavallo,
J. Baschnagel
Abstract:
Angular correlations in dense solutions and melts of flexible polymer chains are investigated with respect to the distance $r$ between the bonds by comparing quantitative predictions of perturbation calculations with numerical data obtained by Monte Carlo simulation of the bond-fluctuation model. We consider both monodisperse systems and grand-canonical (Flory-distributed) equilibrium polymers.…
▽ More
Angular correlations in dense solutions and melts of flexible polymer chains are investigated with respect to the distance $r$ between the bonds by comparing quantitative predictions of perturbation calculations with numerical data obtained by Monte Carlo simulation of the bond-fluctuation model. We consider both monodisperse systems and grand-canonical (Flory-distributed) equilibrium polymers. Density effects are discussed as well as finite chain length corrections. The intrachain bond-bond correlation function $P(r)$ is shown to decay as $P(r) \sim 1/r^3$ for $ξ\ll r \ll \r^*$ with $ξ$ being the screening length of the density fluctuations and $r^* \sim N^{1/3}$ a novel length scale increasing slowly with (mean) chain length $N$.
△ Less
Submitted 17 December, 2009;
originally announced December 2009.
-
A finite excluded volume bond-fluctuation model: Static properties of dense polymer melts revisited
Authors:
J. P. Wittmer,
A. Cavallo,
T. Kreer,
J. Baschnagel,
A. Johner
Abstract:
The classical bond-fluctuation model (BFM) is an efficient lattice Monte Carlo algorithm for coarse-grained polymer chains where each monomer occupies exclusively a certain number of lattice sites. In this paper we propose a generalization of the BFM where we relax this constraint and allow the overlap of monomers subject to a finite energy penalty $\overlap$. This is done to vary systematically…
▽ More
The classical bond-fluctuation model (BFM) is an efficient lattice Monte Carlo algorithm for coarse-grained polymer chains where each monomer occupies exclusively a certain number of lattice sites. In this paper we propose a generalization of the BFM where we relax this constraint and allow the overlap of monomers subject to a finite energy penalty $\overlap$. This is done to vary systematically the dimensionless compressibility $g$ of the solution in order to investigate the influence of density fluctuations in dense polymer melts on various s tatic properties at constant overall monomer density. The compressibility is obtained directly from the low-wavevector limit of the static structure fa ctor. We consider, e.g., the intrachain bond-bond correlation function, $P(s)$, of two bonds separated by $s$ monomers along the chain. It is shown that the excluded volume interactions are never fully screened for very long chains. If distances smaller than the thermal blob size are probed ($s \ll g$) the chains are swollen acc ording to the classical Fixman expansion where, e.g., $P(s) \sim g^{-1}s^{-1/2}$. More importantly, the polymers behave on larger distances ($s \gg g$) like swollen chains of incompressible blobs with $P(s) \si m g^0s^{-3/2}$.
△ Less
Submitted 11 August, 2009;
originally announced August 2009.
-
Perimeter Length and Form Factor of Two-Dimensional Polymer Melts
Authors:
H. Meyer,
T. Kreer,
M. Aichele,
A. Cavallo,
A. Johner,
J. Baschnagel,
J. P. Wittmer
Abstract:
Self-avoiding polymers in two-dimensional ($d=2$) melts are known to adopt compact configurations of typical size $R(N) \sim N^{1/d}$ with $N$ being the chain length. Using molecular dynamics simulations we show that the irregular shapes of these chains are characterized by a perimeter length $L(N) \sim R(N)^{\dpm}$ of fractal dimension $\dpm = d-Θ_2 =5/4$ with $Θ_2=3/4$ being a well-known conta…
▽ More
Self-avoiding polymers in two-dimensional ($d=2$) melts are known to adopt compact configurations of typical size $R(N) \sim N^{1/d}$ with $N$ being the chain length. Using molecular dynamics simulations we show that the irregular shapes of these chains are characterized by a perimeter length $L(N) \sim R(N)^{\dpm}$ of fractal dimension $\dpm = d-Θ_2 =5/4$ with $Θ_2=3/4$ being a well-known contact exponent. Due to the self-similar structure of the chains, compactness and perimeter fractality repeat for subchains of all arc-lengths $s$ down to a few monomers. The Kratky representation of the intramolecular form factor $F(q)$ reveals a strong non-monotonous behavior with $q^2F(q) \sim 1/(qN^{1/d})^{Θ_2}$ in the intermediate regime of the wavevector $q$. Measuring the scattering of labeled subchains %($s F(q) \sim L(s)$) the form factor may allow to test our predictions in real experiments.
△ Less
Submitted 7 May, 2009;
originally announced May 2009.
-
Structural and conformational dynamics of supercooled polymer melts: Insights from first-principles theory and simulations
Authors:
Song-Ho Chong,
Martin Aichele,
Hendrik Meyer,
Matthias Fuchs,
Jörg Baschnagel
Abstract:
We report on quantitative comparisons between simulation results of a bead-spring model and mode-coupling theory calculations for the structural and conformational dynamics of a supercooled, unentangled polymer melt. We find semiquantitative agreement between simulation and theory, except for processes that occur on intermediate length scales between the compressibility plateau and the amorphous…
▽ More
We report on quantitative comparisons between simulation results of a bead-spring model and mode-coupling theory calculations for the structural and conformational dynamics of a supercooled, unentangled polymer melt. We find semiquantitative agreement between simulation and theory, except for processes that occur on intermediate length scales between the compressibility plateau and the amorphous halo of the static structure factor. Our results suggest that the onset of slow relaxation in a glass-forming melt can be described in terms of monomer-caging supplemented by chain connectivity. Furthermore, a unified atomistic description of glassy arrest and of conformational fluctuations that (asymptotically) follow the Rouse model, emerges from our theory.
△ Less
Submitted 5 December, 2007;
originally announced December 2007.
-
Static Rouse Modes and Related Quantities: Corrections to Chain Ideality in Polymer Melts
Authors:
H. Meyer,
J. P. Wittmer,
T. Kreer,
P. Beckrich,
A. Johner,
J. Farago,
J. Baschnagel
Abstract:
Following the Flory ideality hypothesis intrachain and interchain excluded volume interactions are supposed to compensate each other in dense polymer systems. Multi-chain effects should thus be neglected and polymer conformations may be understood from simple phantom chain models. Here we provide evidence against this phantom chain, mean-field picture. We analyze numerically and theoretically th…
▽ More
Following the Flory ideality hypothesis intrachain and interchain excluded volume interactions are supposed to compensate each other in dense polymer systems. Multi-chain effects should thus be neglected and polymer conformations may be understood from simple phantom chain models. Here we provide evidence against this phantom chain, mean-field picture. We analyze numerically and theoretically the static correlation function of the Rouse modes. Our numerical results are obtained from computer simulations of two coarse-grained polymer models for which the strength of the monomer repulsion can be varied, from full excluded volume (`hard monomers') to no excluded volume (`phantom chains'). For nonvanishing excluded volume we find the simulated correlation function of the Rouse modes to deviate markedly from the predictions of phantom chain models. This demonstrates that there are nonnegligible correlations along the chains in a melt. These correlations can be taken into account by perturbation theory. Our simulation results are in good agreement with these new theoretical predictions.
△ Less
Submitted 4 December, 2007;
originally announced December 2007.
-
Soft particle model for block copolymers
Authors:
F. Eurich,
A. Karatchentsev,
J. Baschnagel,
W. Dieterich,
P. Maass
Abstract:
A soft particle model for diblock (AB) copolymer melts is proposed. Each molecule is mapped onto two soft spheres built by Gaussian A- and B-monomer distributions. An approximate analytical expression for the joint distribution function for the distance between both spheres and their radii of gyration is derived, which determines the entropic contribution to the intramolecular free energy. Addin…
▽ More
A soft particle model for diblock (AB) copolymer melts is proposed. Each molecule is mapped onto two soft spheres built by Gaussian A- and B-monomer distributions. An approximate analytical expression for the joint distribution function for the distance between both spheres and their radii of gyration is derived, which determines the entropic contribution to the intramolecular free energy. Adding a mean-field expression for the intermolecular interactions, we obtain the total free energy of the system. Based on this free energy, Monte Carlo simulations are carried out to study the kinetics of microphase ordering in the bulk and its effect on molecular diffusion. This is followed by an analysis of thin films with emphasis on pattern transfer from walls with a periodic structure. It is shown that the level of coarse graining in the soft particle model is suitable to describe structural and kinetic properties of copolymers on mesoscopic scales.
△ Less
Submitted 8 June, 2007;
originally announced June 2007.
-
Bond formation and slow heterogeneous dynamics in adhesive spheres with long--ranged repulsion: Quantitative test of Mode Coupling Theory
Authors:
O. Henrich,
A. M. Puertas,
M. Sperl,
J. Baschnagel,
M. Fuchs
Abstract:
A colloidal system of spheres interacting with both a deep and narrow attractive potential and a shallow long-ranged barrier exhibits a prepeak in the static structure factor. This peak can be related to an additional mesoscopic length scale of clusters and/or voids in the system. Simulation studies of this system have revealed that it vitrifies upon increasing the attraction into a gel-like sol…
▽ More
A colloidal system of spheres interacting with both a deep and narrow attractive potential and a shallow long-ranged barrier exhibits a prepeak in the static structure factor. This peak can be related to an additional mesoscopic length scale of clusters and/or voids in the system. Simulation studies of this system have revealed that it vitrifies upon increasing the attraction into a gel-like solid at intermediate densities. The dynamics at the mesoscopic length scale corresponding to the prepeak represents the slowest mode in the system. Using mode coupling theory with all input directly taken from simulations, we reveal the mechanism for glassy arrest in the system at 40% packing fraction. The effects of the low-q peak and of polydispersity are considered in detail. We demonstrate that the local formation of physical bonds is the process whose slowing down causes arrest.
It remains largely unaffected by the large-scale heterogeneities, and sets the clock for the slow cluster mode. Results from mode-coupling theory without adjustable parameters agree semi-quantitatively with the local density correlators but overestimate the lifetime of the mesoscopic structure (voids).
△ Less
Submitted 4 May, 2007;
originally announced May 2007.
-
Intramolecular long-range correlations in polymer melts: The segmental size distribution and its moments
Authors:
J. P. Wittmer,
P. Beckrich,
H. Meyer,
A. Cavallo,
A. Johner,
J. Baschnagel
Abstract:
Presenting theoretical arguments and numerical results we demonstrate long-range intrachain correlations in concentrated solutions and melts of long flexible polymers which cause a systematic swelling of short chain segments. They can be traced back to the incompressibility of the melt leading to an effective repulsion $u(s) \approx s/ρR^3(s) \approx ce/\sqrt{s}$ when connecting two segments tog…
▽ More
Presenting theoretical arguments and numerical results we demonstrate long-range intrachain correlations in concentrated solutions and melts of long flexible polymers which cause a systematic swelling of short chain segments. They can be traced back to the incompressibility of the melt leading to an effective repulsion $u(s) \approx s/ρR^3(s) \approx ce/\sqrt{s}$ when connecting two segments together where $s$ denotes the curvilinear length of a segment, $R(s)$ its typical size, $ce \approx 1/ρbe^3$ the ``swelling coefficient", $be$ the effective bond length and $ρ$ the monomer density. The relative deviation of the segmental size distribution from the ideal Gaussian chain behavior is found to be proportional to $u(s)$. The analysis of different moments of this distribution allows for a precise determination of the effective bond length $be$ and the swelling coefficient $ce$ of asymptotically long chains. At striking variance to the short-range decay suggested by Flory's ideality hypothesis the bond-bond correlation function of two bonds separated by $s$ monomers along the chain is found to decay algebraically as $1/s^{3/2}$. Effects of finite chain length are considered briefly.
△ Less
Submitted 24 May, 2007; v1 submitted 12 April, 2007;
originally announced April 2007.
-
Why polymer chains in a melt are not random walks
Authors:
J. P. Wittmer,
P. Beckrich,
J. Johner,
A. N. Semenov,
S. P. Obukhov,
H. Meyer,
J. Baschnagel
Abstract:
A cornerstone of modern polymer physics is the `Flory ideality hypothesis' which states that a chain in a polymer melt adopts `ideal' random-walk-like conformations. Here we revisit theoretically and numerically this pivotal assumption and demonstrate that there are noticeable deviations from ideality. The deviations come from the interplay of chain connectivity and the incompressibility of the…
▽ More
A cornerstone of modern polymer physics is the `Flory ideality hypothesis' which states that a chain in a polymer melt adopts `ideal' random-walk-like conformations. Here we revisit theoretically and numerically this pivotal assumption and demonstrate that there are noticeable deviations from ideality. The deviations come from the interplay of chain connectivity and the incompressibility of the melt, leading to an effective repulsion between chain segments of all sizes $s$. The amplitude of this repulsion increases with decreasing $s$ where chain segments become more and more swollen. We illustrate this swelling by an analysis of the form factor $F(q)$, i.e. the scattered intensity at wavevector $q$ resulting from intramolecular interferences of a chain. A `Kratky plot' of $q^2F(q)$ {\em vs.} $q$ does not exhibit the plateau for intermediate wavevectors characteristic of ideal chains. One rather finds a conspicuous depression of the plateau, $δ(F^{-1}(q)) = |q|^3/32ρ$, which increases with $q$ and only depends on the monomer density $ρ$.
△ Less
Submitted 11 January, 2007; v1 submitted 13 November, 2006;
originally announced November 2006.
-
On the Dynamics and Disentanglement in Thin and Two-Dimensional Polymer Films
Authors:
H. Meyer,
T. Kreer,
A. Cavallo,
J. P. Wittmer,
J. Baschnagel
Abstract:
We present results from molecular dynamics simulations of strictly two-dimensional (2D) polymer melts and thin polymer films in a slit geometry of thickness of the order of the radius of gyration. We find that the dynamics of the 2D melt is qualitatively different from that of the films. The 2D monomer mean-square displacement shows a $t^{8/15}$ power law at intermediate times instead of the…
▽ More
We present results from molecular dynamics simulations of strictly two-dimensional (2D) polymer melts and thin polymer films in a slit geometry of thickness of the order of the radius of gyration. We find that the dynamics of the 2D melt is qualitatively different from that of the films. The 2D monomer mean-square displacement shows a $t^{8/15}$ power law at intermediate times instead of the $t^{1/2}$ law expected from Rouse theory for nonentangled chains. In films of finite thickness, chain entanglements may occur. The impact of confinement on the entanglement length $N_\mathrm{e}$ has been analyzed by a primitive path analysis. The analysis reveals that $N_\mathrm{e}$ increases strongly with decreasing film thickness.
△ Less
Submitted 6 September, 2006;
originally announced September 2006.
-
Monte Carlo SImulation of Polymers: Coarse-Grained Models
Authors:
J. Baschnagel,
J. P. Wittmer,
H. Meyer
Abstract:
A coarse-grained simulation model eliminates microscopic degrees of freedom and represents a polymer by a simplified structure. A priori, two classes of coarse-grained models may be distinguished: those which are designed for a specific polymer and reflect the underlying atomistic details to some extent, and those which retain only the most basic features of a polymer chain (chain connectivity,…
▽ More
A coarse-grained simulation model eliminates microscopic degrees of freedom and represents a polymer by a simplified structure. A priori, two classes of coarse-grained models may be distinguished: those which are designed for a specific polymer and reflect the underlying atomistic details to some extent, and those which retain only the most basic features of a polymer chain (chain connectivity, short-range excluded-volume interactions, etc.).
In this review we mainly focus on the second class of generic polymer models, while the first class of specific coarse-grained models is only touched upon briefly.
△ Less
Submitted 28 July, 2004;
originally announced July 2004.
-
Long Range Bond-Bond Correlations in Dense Polymer Solutions
Authors:
J. P. Wittmer,
H. Meyer,
J. Baschnagel,
A. Johner,
S. Obukhov,
L. Mattioni,
M. Mueller,
A. N. Semenov
Abstract:
The scaling of the bond-bond correlation function $C(s)$ along linear polymer chains is investigated with respect to the curvilinear distance, $s$, along the flexible chain and the monomer density, $ρ$, via Monte Carlo and molecular dynamics simulations. % Surprisingly, the correlations in dense three dimensional solutions are found to decay with a power law $C(s) \sim s^{-ω}$ with $ω=3/2$ and t…
▽ More
The scaling of the bond-bond correlation function $C(s)$ along linear polymer chains is investigated with respect to the curvilinear distance, $s$, along the flexible chain and the monomer density, $ρ$, via Monte Carlo and molecular dynamics simulations. % Surprisingly, the correlations in dense three dimensional solutions are found to decay with a power law $C(s) \sim s^{-ω}$ with $ω=3/2$ and the exponential behavior commonly assumed is clearly ruled out for long chains. % In semidilute solutions, the density dependent scaling of $C(s) \approx g^{-ω_0} (s/g)^{-ω}$ with $ω_0=2-2ν=0.824$ ($ν=0.588$ being Flory's exponent) is set by the number of monomers $g(ρ)$ contained in an excluded volume blob of size $ξ$. % Our computational findings compare well with simple scaling arguments and perturbation calculation. The power-law behavior is due to self-interactions of chains on distances $s \gg g$ caused by the connectivity of chains and the incompressibility of the melt. %
△ Less
Submitted 21 June, 2004; v1 submitted 20 April, 2004;
originally announced April 2004.
-
Static Properties of a Simulated Supercooled Polymer Melt: Structure Factors, Monomer Distributions Relative to the Center of Mass, and Triple Correlation Functions
Authors:
Martin Aichele,
Song-Ho Chong,
Jörg Baschnagel,
Matthias Fuchs
Abstract:
We analyze structural and conformational properties in a simulated bead-spring model of a non-entangled, supercooled polymer melt. We explore the statics of the model via various structure factors, involving not only the monomers, but also the center of mass (CM). We find that the conformation of the chains and the CM-CM structure factor, which is well described by a recently proposed approximat…
▽ More
We analyze structural and conformational properties in a simulated bead-spring model of a non-entangled, supercooled polymer melt. We explore the statics of the model via various structure factors, involving not only the monomers, but also the center of mass (CM). We find that the conformation of the chains and the CM-CM structure factor, which is well described by a recently proposed approximation [Krakoviack et al., Europhys. Lett. 58, 53 (2002)], remain essentially unchanged on cooling toward the critical glass transition temperature of mode-coupling theory. Spatial correlations between monomers on different chains, however, depend on temperature, albeit smoothly. This implies that the glassy behavior of our model cannot result from static intra-chain or CM-CM correlations. It must be related to inter-chain correlations at the monomer level. Additionally, we study the dependence of inter-chain correlation functions on the position of the monomer along the chain backbone. We find that this site-dependence can be well accounted for by a theory based on the polymer reference interaction site model (PRISM). We also analyze triple correlations by means of the three-monomer structure factors for the melt and for the chains. These structure factors are compared with the convolution approximation that factorizes them into a product of two-monomer structure factors. For the chains this factorization works very well, indicating that chain connectivity does not introduce special triple correlations in our model. For the melt deviations are more pronounced, particularly at wave vectors close to the maximum of the static structure factor.
△ Less
Submitted 25 February, 2004;
originally announced February 2004.
-
Dynamical Properties of the Slithering Snake Algorithm: A numerical test of the activated reptation hypothesis
Authors:
L. Mattioni,
J. P. Wittmer,
J. Baschnagel,
J. -L. Barrat,
E. Luijten
Abstract:
The correlations in the motion of reptating polymers in their melt are investigated by means of kinetic Monte Carlo simulations of the three dimensional slithering snake version of the bond-fluctuation model.
Surprisingly, the slithering snake dynamics becomes inconsistent with classical reptation predictions at high chain overlap (either chain length $N$ or volume fraction $φ$) where the rela…
▽ More
The correlations in the motion of reptating polymers in their melt are investigated by means of kinetic Monte Carlo simulations of the three dimensional slithering snake version of the bond-fluctuation model.
Surprisingly, the slithering snake dynamics becomes inconsistent with classical reptation predictions at high chain overlap (either chain length $N$ or volume fraction $φ$) where the relaxation times increase much faster than expected.
This is due to the anomalous curvilinear diffusion in a finite time window whose upper bound $τ_+$ is set by the chain end density $φ/N$. Density fluctuations created by passing chain ends allow a reference polymer to break out of the local cage of immobile obstacles created by neighboring chains.
The dynamics of dense solutions of snakes at $t \ll τ_+$ is identical to that of a benchmark system where all but one chain are frozen. We demonstrate that it is the slow creeping of a chain out of its correlation hole which causes the subdiffusive dynamical regime.
Our results are in good qualitative agreement with the activated reptation scheme proposed recently by Semenov and Rubinstein [Eur. Phys. J. B, {\bf 1} (1998) 87].
Additionally, we briefly comment on the relevance of local relaxation pathways within a slithering snake scheme. Our preliminary results suggest that a judicious choice of the ratio of local to slithering snake moves is crucial to equilibrate a melt of long chains efficiently.
△ Less
Submitted 22 March, 2003; v1 submitted 17 December, 2002;
originally announced December 2002.
-
Gaussian ellipsoid model for confined polymer systems
Authors:
F. Eurich,
P. Maass,
J. Baschnagel
Abstract:
Polymer systems in slab geometries are studied on the basis of the recently presented Gaussian Ellipsoid Model [J. Chem. Phys. 114, 7655 (2001)].The potential of the confining walls has an exponential shape. For homogeneous systems in thermodynamic equilibrium we discuss density, orientation and deformation profiles of the polymers close to the walls. For strongly segregated mixtures of polymer…
▽ More
Polymer systems in slab geometries are studied on the basis of the recently presented Gaussian Ellipsoid Model [J. Chem. Phys. 114, 7655 (2001)].The potential of the confining walls has an exponential shape. For homogeneous systems in thermodynamic equilibrium we discuss density, orientation and deformation profiles of the polymers close to the walls. For strongly segregated mixtures of polymer components A and B equilibrium profiles are studied near a planar interface separating A and B rich regions. Spinodal decomposition processes of the mixtures in the presence of neutral walls show upon strong confinement an increase of the lateral size of A and B rich domains and a slowing down of the demixing kinetics. These findings are in agreement with predictions from time dependent Ginzburg--Landau theory. In the case, where one wall periodically favors one of the two mixture components over the other, different equilibrium structures emerge and lead to different kinetic pathways of spinodal decomposition processes in such systems.
△ Less
Submitted 25 February, 2002;
originally announced February 2002.
-
Simulation of models for the glass transition: Is there progress?
Authors:
Kurt Binder,
Jörg Baschnagel,
Walter Kob,
Wolfgang Paul
Abstract:
The glass transition of supercooled fluids is a particular challenge for computer simulation, because the (longest) relaxation times increase by about 15 decades upon approaching the transition temperature T_g. Brute-force molecular dynamics simulations, as presented here for molten SiO_2 and coarse-grained bead-spring models of polymer chains, can yield very useful insight about the first few d…
▽ More
The glass transition of supercooled fluids is a particular challenge for computer simulation, because the (longest) relaxation times increase by about 15 decades upon approaching the transition temperature T_g. Brute-force molecular dynamics simulations, as presented here for molten SiO_2 and coarse-grained bead-spring models of polymer chains, can yield very useful insight about the first few decades of this slowing down. Hence this allows to access the temperature range around T_c of the so-called mode coupling theory, whereas the dynamics around the experimental glass transition is completely out of reach. While methods such as ``parallel tempering'' improve the situation somewhat, a method that allows to span a significant part of the region T_g\leq T\leq T_c is still lacking. Only for abstract models such as the infinite range 10-state Potts glass with a few hundred spins this region can be explored. However this model suffers from very strong finite size effects thus making it difficult to extrapolate the results obtained for the finite system sizes to the thermodynamic limit. For the case of polymer melts, two different strategies to use lattice models instead of continuum models are discussed: In the first approach, a mapping of an atomistically realistic model of polyethylene to the bond fluctuation model with suitable effective potentials and a temperature-dependent time rescaling factor is attempted. In the second approach, devoted to a test of the entropy theory, moves that are artificial but which lead to a faster relaxation (``slithering snake'' algorithm) are used, to get at least static properties at somewhat lower temperatures than possible with a ``realistic'' dynamics. The merits and shortcomings of all these approaches are discussed.
△ Less
Submitted 20 February, 2002;
originally announced February 2002.
-
Recent MD Results on Supercooled Thin Polymer Films
Authors:
Fathollah Varnik,
Joerg Baschnagel,
Kurt Binder
Abstract:
The dynamic and static properties of a supercooled (non-entangled) polymer melt are investigated via molecular dynamics (MD) simulations. The system is confined between two completely smooth and purely repulsive walls. The wall-to-wall separation (film thickness), $D$, is varied from about 3 to about 14 times the bulk radius of gyration. Despite the geometric confinement, the supercooled films e…
▽ More
The dynamic and static properties of a supercooled (non-entangled) polymer melt are investigated via molecular dynamics (MD) simulations. The system is confined between two completely smooth and purely repulsive walls. The wall-to-wall separation (film thickness), $D$, is varied from about 3 to about 14 times the bulk radius of gyration. Despite the geometric confinement, the supercooled films exhibit many qualitative features which were also observed in the bulk and could be analyzed in terms of mode-coupling theory (MCT). Examples are the two-step relaxation of the incoherent intermediate scattering function, the time-temperature superposition property of the late time $α$-process and the space-time factorization of the scattering function on the intermediate time scale of the MCT $β$-process. An analysis of the temperature dependence of the $α$-relaxation time suggests that the critical temperature, $\Tc$, of MCT decreases with $D$. If the confinement is not too strong ($D \ge 10 {monomer diameter}$) the static structure factor of the film coincides with that of the bulk when compared for the same distance, $T-\Tc(D)$, to the critical temperature. This suggests that $T-\Tc(D)$ is an important temperature scale of our model both in the bulk and in the films.
△ Less
Submitted 11 September, 2001;
originally announced September 2001.