-
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.
-
Ripple-like instability in the simulated gel phase of finite size phosphocholine bilayers
Authors:
Vivien Walter,
Céline Ruscher,
Adrien Gola,
Carlos M. Marques,
Olivier Benzerara,
Fabrice Thalmann
Abstract:
Atomistic molecular dynamics simulations have reached a degree of maturity that makes it possible to investigate the lipid polymorphism of model bilayers over a wide range of temperatures. However if both the fluid $L_α$ and tilted gel $L_{β'}$ states are routinely obtained, the $P_{β'}$ ripple phase of phosphatidylcholine lipid bilayers is still unsatifactorily described. Performing simulations o…
▽ More
Atomistic molecular dynamics simulations have reached a degree of maturity that makes it possible to investigate the lipid polymorphism of model bilayers over a wide range of temperatures. However if both the fluid $L_α$ and tilted gel $L_{β'}$ states are routinely obtained, the $P_{β'}$ ripple phase of phosphatidylcholine lipid bilayers is still unsatifactorily described. Performing simulations of lipid bilayers made of different numbers of DPPC (1,2-dipalmitoylphosphatidylcholine) molecules ranging from 32 to 512, we demonstrate that the tilted gel phase $L_{β'}$ expected below the pre-transition cannot be obtained for large systems ($>$ 94 DPPC molecules) through common simulations settings or temperature treatments. Large systems are instead found in a disordered gel phase which display configurations, topography and energies reminiscent from the ripple phase $P_{β'}$ observed between the pretransition and the main melting transition. We show how the state of the bilayers below the pretransition can be controlled and depends on thermal history and conditions of preparations. A mechanism for the observed topographic instability is suggested.
△ Less
Submitted 7 February, 2021;
originally announced February 2021.
-
A machine learning assessment of the two states model for lipid bilayer phase transitions
Authors:
Viven Walter,
Céline Ruscher,
Olivier Benzerara,
Carlos M. Marques,
Fabrice Thalmann
Abstract:
We have adapted a set of classification algorithms, also known as Machine Learning, to the identification of fluid and gel domains close to the main transition of dipalmitoyl-phosphatidylcholine (DPPC) bilayers. Using atomistic molecular dynamics conformations in the low and high temperature phases as learning sets, the algorithm was trained to categorize individual lipid configurations as fluid o…
▽ More
We have adapted a set of classification algorithms, also known as Machine Learning, to the identification of fluid and gel domains close to the main transition of dipalmitoyl-phosphatidylcholine (DPPC) bilayers. Using atomistic molecular dynamics conformations in the low and high temperature phases as learning sets, the algorithm was trained to categorize individual lipid configurations as fluid or gel, in relation with the usual two-states phenomenological description of the lipid melting transition. We demonstrate that our machine can learn and sort lipids according to their most likely state without prior assumption regarding the nature of the order parameter of the transition. Results from our machine learning approach provides strong support in favor of a two-states model approach of membrane fluidity.
△ Less
Submitted 21 September, 2019; v1 submitted 12 July, 2019;
originally announced July 2019.
-
The role of shape disorder in the collective behaviour of aligned fibrous matter
Authors:
Salvatore Salamone,
Nava Schulmann,
Olivier Benzerara,
Hendrik Meyer,
Thierry Charitat,
Carlos M. Marques
Abstract:
We study the compression of bundles of aligned macroscopic fibers with intrinsic shape disorder, as found in human hair and in many other natural and man-made systems. We show by a combination of experiments, numerical simulations and theory how the statistical properties of the shapes of the fibers control the collective mechanical behaviour of the bundles. This work paves the way for designing a…
▽ More
We study the compression of bundles of aligned macroscopic fibers with intrinsic shape disorder, as found in human hair and in many other natural and man-made systems. We show by a combination of experiments, numerical simulations and theory how the statistical properties of the shapes of the fibers control the collective mechanical behaviour of the bundles. This work paves the way for designing aligned fibrous matter with pre-required properties from large numbers of individual strands of selected geometry and rigidity.
△ Less
Submitted 4 March, 2019; v1 submitted 26 September, 2018;
originally announced September 2018.
-
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.
-
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.