Search

We investigate generic inequalities of various contributions to the shear modulus $\mu$ in ensembles of amorphous elastic bodies. We focus first on a simple elastic network model with connectivity matrices (CMs) which are either annealed or quenched, at or out of equilibrium. The stress-fluctuation formalism relation for $\mu$ is rewritten as $\mu = \mu_1 + \mu_a$ with $\mu_1 \ge 0$ characterizing the variance of the quenched shear stresses and $\mu_a$ being a simple average over all states and CMs. For equilibrium CM-distributions $\mu_a$ becomes equivalent to the shear modulus of annealed systems, i.e. $\mu_a \ge 0$, while more generally $\mu_a$ may become strongly negative as shown by considering a temperature quench and a scalar active two-temperature model. Consistent relations are also found for glass-forming colloids where $\mu-\mu_1=\mu_a=0$ for equilibrium ensembles, i.e. $\mu$ is set by the quenched shear stresses, while $\mu_a$ becomes again negative otherwise., Comment: 7 pages, 6 figures

Generalizing recent work on isotropic tensor fields in isotropic and achiral condensed matter systems from two to arbitrary dimensions we address both mathematical aspects assuming perfectly isotropic systems and applications focusing on correlation functions of displacement and strain field components in amorphous solids. Various general points are exemplified using simulated polydisperse Lennard-Jones particles in two dimensions. It is shown that the strain components in reciprocal space have essentially a complex circularly-symmetric Gaussian distribution albeit weak non-Gaussianity effects become visible for large wavenumbers q where also anisotropy effects become relevant. The dynamical strain correlation functions are strongly non-monotonic with respect to q with a minimum roughly at the breakdown of the continuum limit., Comment: 22 pages, 11 figures, 1 tables

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, Comment: 17 pages, 7 figures. arXiv admin note: text overlap with arXiv:2303.16571

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., Comment: 21 pages, 9 figures

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., Comment: 15 pages, 14 figures

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) = \mu_A - h(t)$ for the computational determination of the shear stress relaxation modulus $G(t)$. Here, $\mu_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., Comment: 5 pages, 2 figures. arXiv admin note: text overlap with arXiv:1711.00725

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 $\tau_p$ of mode $p$ must scale as $\tau_p \sim (N/p)^{5/3}$ rather than the usual square power law for linear chains., Comment: 14 pages, 11 figures

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 $\delta v$ of the empirical variance $v[x]$ of finite time series $x$. $x_t$ is averaged over a fluctuating field $\sigma_{r}$ ($r$ being the microcell position) characterized by a quenched spatially correlated Gaussian field. Due to the quenched field $\delta v(\Delta t)$ becomes a finite constant, $\Delta_{ne} > 0$, for large sampling times $\Delta t$. The volume dependence of the non-ergodicity parameter $\Delta_{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., Comment: 11 pages, 8 figures

We investigate the standard deviation $\delta 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, Comment: 15 pages, 10 figures

Extending recent work on stress fluctuations in complex fluids and amorphous solids we describe in general terms the ensemble average $v(\Delta t)$ and the standard deviation $\delta v(\Delta t)$ of the variance $v[\mathbf{x}]$ of time series $\mathbf{x}$ of a stochastic process $x(t)$ measured over a finite sampling time $\Delta t$. Assuming a stationary, Gaussian and ergodic process, $\delta v$ is given by a functional $\delta v_G[h]$ of the autocorrelation function $h(t)$. $\delta v(\Delta t)$ is shown to become large and similar to $v(\Delta t)$ if $\Delta t$ corresponds to a fast relaxation process. Albeit $\delta v = \delta 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., Comment: 19 pages, 18 figures, submitted to EPJE

Using molecular dynamics simulation of a polymer glass model we investigate free-standing polymer films focusing on the in-plane shear modulus $\mu$ 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 $\Delta 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, $\mu$ 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 $\Delta t$-dependence of $\mu$ is traced back to $G(t)$., Comment: 15 pages, 14 figures, accepted to Phys. Rev. E (Nov 2018) https://journals.aps.org/pre/accepted/b9077R77J041fb19b5843389c816aa554b5d02468

Using molecular dynamics simulation of a standard coarse-grained polymer glass model we investigate by means of the stress-fluctuation formalism the shear modulus $\mu$ as a function of temperature $T$ and sampling time $\Delta t$. While the ensemble-averaged modulus $\mu(T)$ is found to decrease continuously for all $\Delta t$ sampled, its standard deviation $\delta \mu(T)$ is non-monotonous with a striking peak at the glass transition. Confirming the effective time-translational invariance of our systems, $\mu(\Delta t)$ can be understood using a weighted integral over the shear-stress relaxation modulus $G(t)$. While the crossover of $\mu(T)$ gets sharper with increasing $\Delta t$, the peak of $\delta \mu(T)$ becomes more singular. % It is thus elusive to predict the modulus of a single configuration at the glass transition., Comment: 5 pages, 6 figure, published at PRL

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 $\Delta 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 $\mu_{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 $\mu_{sf}(T)$ decreases continuously with increasing T for all $\Delta t$ sampled, its standard deviation $\delta \mu_{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 $\Delta t$-dependence of $\mu_{sf}$ and related quantities can be understood using a weighted integral over $G(t)$. This implies that the shear viscosity $\eta(T)$ may be readily obtained from the $1/\Delta t$-decay of $\mu_{sf}$ above the glass transition., Comment: 21 pages, 21 figures, submitted to PRE

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 $\Delta 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 $\Delta x = dt/t_*(f)$ is systematically scanned from the liquid limit ($\Delta dx \gg 1)$ to the solid limit ($\Delta 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 $\Delta x$ using the simple-average expression $G(t) = \mu_A - h(t)$ with $\mu_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., Comment: 12 pages, 11 figures

Focusing on isotropic elastic networks we propose a novel simple-average expression $G(t) = \mu_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, $\mu_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{\tau}(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., Comment: 6 pages, 4 figures

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., Comment: 13 pages, 14 figures

The shear stress relaxation modulus $G(t)$ may be determined from the shear stress $\tau(t)$ after switching on a tiny step strain $\gamma$ or by inverse Fourier transformation of the storage modulus $G^{\prime}(\omega)$ or the loss modulus $G^{\prime\prime}(\omega)$ obtained in a standard oscillatory shear experiment at angular frequency $\omega$. It is widely assumed that $G(t)$ is equivalent in general to the equilibrium stress autocorrelation function $C(t) = \beta V \langle \delta \tau(t) \delta \tau(0)\rangle$ which may be readily computed in computer simulations ($\beta$ 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}(\omega)$. $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)$., Comment: 13 pages, 15 figures

We revisit the relation between the shear stress relaxation modulus $G(t)$, computed at finite shear strain $0 < \gamma \ll 1$, and the shear stress autocorrelation functions $C(t)|_{\gamma}$ and $C(t)|_{\tau}$ computed, respectively, at imposed strain $\gamma$ and mean stress $\tau$. Focusing on permanent isotropic spring networks it is shown theoretically and computationally that in general $G(t) = C(t)|_{\tau} = C(t)|_{\gamma} + G_{eq}$ for $t > 0$ with $G_{eq}$ being the static equilibrium shear modulus. $G(t)$ and $C(t)|_{\gamma}$ thus must become different for solids and it is impossible to obtain $G_{eq}$ alone from $C(t)|_{\gamma}$ 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)|_{\tau} = C(t)|_{\gamma}$ should reveal an intermediate plateau set by the shear modulus $G_{eq}(f=0)$ of the quenched network., Comment: 8 pages, 4 figures

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, $\lambda=0$) or volume (NVT-ensemble, $\lambda=1$) and for more general values of the dimensionless parameter $\lambda$ characterizing the constant-volume constraint., Comment: 17 pages, 9 figures

Shear-strain and shear-stress correlations in isotropic elastic bodies are investigated both theoretically and numerically at either imposed mean shear-stress $\tau$ ($\lambda=0$) or shear-strain $\gamma$ ($\lambda=1$) and for more general values of a dimensionless parameter $\lambda$ characterizing the generalized Gaussian ensemble. It allows to tune the strain fluctuations $\mu_{\gamma\gamma} \equiv \beta V \la \delta \gamma^2 \ra = (1-\lambda)/G_{eq}$ with $\beta$ being the inverse temperature, $V$ the volume, $\gamma$ 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 $\mu_{\tau\tau} \equiv \beta V \la \delta\tau^2 \ra$ ($\tau$ being the instantaneous stress) that $\mu_{\tau\tau} = \mu_{A} - \lambda G_{eq}$ with $\mu_{A} = \mu_{\tau\tau}|_{\lambda=0}$ being the affine shear-elasticity. For the stress autocorrelation function $c_{\tau\tau}(t) \equiv \beta V \la \delta \tau(t) \delta \tau(0) \ra$ this result is then seen (assuming a sufficiently slow shear-stress barostat) to generalize to $c_{\tau\tau}(t) = G(t) - \lambda \Geq$ with $G(t)$ being the shear-stress relaxation modulus., Comment: 17 pages, 15 figures

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.

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.

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.

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\rho^(21/8), the total pressure P\sim\rho^3 and the dimensionless compressibility gT=lim(q->0)(S(q)\sim1/\rho^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.

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.

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.

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 $\xi$ 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 \xi$. 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., Comment: 60 pages, 49 figures

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^{\nu}$ with $\nu=1/d$. % The precise measurement of various distributions of internal chain distances allows a direct test of the contact exponents $\Theta_0=3/8$, $\Theta_1=1/2$ and $\Theta_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+\nu\Theta_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-\Theta_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))^{\Theta_2}$ in the intermediate regime of the wave vector $q$. This may allow to confirm the predicted contour fractality in a real experiment., Comment: 14 figures, accepted at J.Chem.Phys.

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 $\delta \mu_c(n) \sim +1/\rho\sqrt{n}$ arises ($\rho$ being the monomer density) for all $\P(N)$ if $n \ll $ 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 $\delta \mu_c(n) \approx (1- 2 n/) / \rho \sqrt{n}$ for $n \ll ^2$, hence, $\delta \mu_c(n) \approx - 1/\rho \sqrt{n}$ if $n$ is similar to the typical length of the bath $$. 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 - /p!^p$ to decay as $K_p \approx 1 / \sqrt{}$ for all moments $p$ of the distribution., Comment: 14 pages, 6 figures, submitted to EPJE

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 $\xi \ll r \ll \r^*$ with $\xi$ 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$., Comment: 17 pages, 5 figures, accepted for publication at Macromolecules

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}$., Comment: 46 pages, 12 figures

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-\Theta_2 =5/4$ with $\Theta_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})^{\Theta_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., Comment: 4 pages, 4 figures, accepted for PRE Rapid Communications

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., Comment: 9 pages, 7 figures, accepted for publication in EPJE

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/\rho 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/\rho be^3$ the ``swelling coefficient", $be$ the effective bond length and $\rho$ 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., Comment: 50 pages, 13 figures, submitted to PRE

The question about the existence of a structural glass transition in two dimensions is studied using mode coupling theory (MCT). We determine the explicit d-dependence of the memory functional of mode coupling for one-component systems. Applied to two dimensions we solve the MCT equations numerically for monodisperse hard discs. A dynamic glass transition is found at a critical packing fraction phi_c^{d=2} = 0.697 which is above phi_c^{d=3} = 0.516 by about 35%. phi^d_c scales approximately with phi^d_{\rm rcp} the value for random close packing, at least for d=2, 3. Quantities characterizing the local, cooperative 'cage motion' do not differ much for d=2 and d=3, and we e.g. find the Lindemann criterion for the localization length at the glass transition. The final relaxation obeys the superposition principle, collapsing remarkably well onto a Kohlrausch law. The d=2 MCT results are in qualitative agreement with existing results from MC and MD simulations. The mean squared displacements measured experimentally for a quasi-two-dimensional binary system of dipolar hard spheres can be described satisfactorily by MCT for monodisperse hard discs over four decades in time provided the experimental control parameter Gamma (which measures the strength of dipolar interactions) and the packing fraction phi are properly related to each other., Comment: 14 pages, 15 figures

We discuss theoretically and numerically the intramolecular form factor $F(q)$ in dense polymer systems. Following Flory's ideality hypothesis, chains in the melt adopt Gaussian configurations and their form factor is supposed to be given by Debye's formula. At striking variance to this, we obtain noticeable (up to 20%) non-monotonic deviations which can be traced back to the incompressibility of dense polymer solutions beyond a local scale. The Kratky plot ($q^2F(q)$ {\it vs.} wavevector $q$) does not exhibit the plateau expected for Gaussian chains in the intermediate $q$-range. One rather finds a significant decrease according to the correction $\delta(F^{-1}(q)) = q^3/32\rho$ that only depends on the concentration $\rho$ of the solution, but neither on the persistence length or the interaction strength., Comment: 29 pages, 11 figures, submitted to Macromolecules

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, $\delta(F^{-1}(q)) = |q|^3/32\rho$, which increases with $q$ and only depends on the monomer density $\rho$., Comment: 4 pages, 4 figures, EPL, accepted January 2007

It is commonly accepted that in concentrated solutions or melts high-molecular weight polymers display random-walk conformational properties without long-range correlations between subsequent bonds. This absence of memory means, for instance, that the bond-bond correlation function, $P(s)$, of two bonds separated by $s$ monomers along the chain should exponentially decay with $s$. Presenting numerical results and theoretical arguments for both monodisperse chains and self-assembled (essentially Flory size-distributed) equilibrium polymers we demonstrate that some long-range correlations remain due to self-interactions of the chains caused by the chain connectivity and the incompressibility of the melt. Suggesting a profound analogy with the well-known long-range velocity correlations in liquids we find, for instance, $P(s)$ to decay algebraically as $s^{-3/2}$. Our study suggests a precise method for obtaining the statistical segment length \bstar in a computer experiment., Comment: 4 pages, 3 figures

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., Comment: 6 pages, 3 figures, proceedings 3rd International Workshop on Dynamics in Confinement (CONFIT 2006)

Self-assembled linear structures like giant cylindrical micelles or discotic molecules in solution stacked in flexible columns are systems reminiscent of polydisperse polymer solutions.These supramolecular polymers have an equilibrium length distribution, the result of a competition between the random breakage of chains and the fusion of chains to generate longer ones. In the present work, we review the basic theoretical concepts of these ``equilibrium polymers" and some of the important results obtained by simulation approaches. We propose a new version of a mesoscopic model in continuous space based on the bead and FENE spring polymer model which is treated by Brownian Dynamics and Monte-Carlo binding/unbinding reversible changes for adjacent monomers in space, characterized by an attempt frequency parameter. For a dilute and a moderately semi-dilute state-points which both correspond to dynamically unentangled regimes, the dynamic properties are found to depend upon $\omega$ through the effective life time $\tau_b$ of the average size chain which, in turn, yields the kinetic reaction coefficients of the mean-field kinetic model proposed by Cates. Simple kinetic theories seem to work for times $ t \geq \tau_b$ while at shorter time, strong dynamical correlation effects are observed., Comment: 17 figures; Lecture notes in Physics, Springer, to be published

Extending recent numerical studies on two dimensional amorphous bodies, we characterize the approach of elastic continuum limit in three dimensional (weakly polydisperse) Lennard-Jones systems. While performing a systematic finite-size analysis (for two different quench protocols) we investigate the non-affine displacement field under external strain, the linear response to an external delta force and the low-frequency harmonic eigenmodes and their density distribution. Qualitatively similar behavior is found as in two dimensions. We demonstrate that the classical elasticity description breaks down below an intermediate length scale $\xi$, which in our system is approximately 23 molecular sizes. This length characterizes the correlations of the non-affine displacement field, the self-averaging of external noise with distance from the source and gives the lower wave length bound for the applicability of the classical eigenfrequency calculations. We trace back the "Boson-peak" of the density of eigenfrequencies (obtained from the velocity auto-correlation function) to the inhomogeneities on wave lengths smaller than $\xi$., Comment: 27 pages, 11 figures, submitted to Phys. Rev. B

Conformational properties of polymer melts confined between two hard structureless walls are investigated by Monte Carlo simulation of the bond-fluctuation model. Parallel and perpendicular components of chain extension, bond-bond correlation function and structure factor are computed and compared with recent theoretical approaches attempting to go beyond Flory's and Silberberg's hypotheses. We demonstrate that for ultrathin films where the thickness, $H$, is smaller than the excluded volume screening length (blob size), $\xi$, the chain size parallel to the walls diverges logarithmically, $R^2/2N \approx b^2 + c \log(N)$ with $c \sim 1/H$. The corresponding bond-bond correlation function decreases like a power law, $C(s) = d/s^{\omega}$ with $s$ being the curvilinear distance between bonds and $\omega=1$. % Upon increasing the film thickness, $H$, we find -- in contrast to Flory's hypothesis -- the bulk exponent $\omega=3/2$ and, more importantly, an {\em decreasing} $d(H)$ that gives direct evidence for an {\em enhanced} self-interaction of chain segments reflected at the walls. Systematic deviations from the Kratky plateau as a function of $H$ are found for the single chain form factor parallel to the walls in agreement with the {\em non-monotonous} behaviour predicted by theory. This structure in the Kratky plateau might give rise to an erroneous estimation of the chain extension from scattering experiments. For large $H$ the deviations are linear with the wave vector, $q$, but are very weak. In contrast, for ultrathin films, $H<\xi$, very strong corrections are found (albeit logarithmic in $q$) suggesting a possible experimental verification of our results., Comment: 16 pages, 7 figures. Dedicated to L. Sch\"afer on the occasion of his 60th birthday

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., Comment: 57 pages, 19 figures, conference proceedings: Computational Soft Matter: From Synthetic Polymers to Proteins, edited by N. Attig et al. (NIC Series, Volume 23, Juelich, 2004), pp. 83-140

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, $\rho$, 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^{-\omega}$ with $\omega=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^{-\omega_0} (s/g)^{-\omega}$ with $\omega_0=2-2\nu=0.824$ ($\nu=0.588$ being Flory's exponent) is set by the number of monomers $g(\rho)$ contained in an excluded volume blob of size $\xi$. % 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. %, Comment: 4 pages, 4 figures

The linear response of two-dimensional amorphous elastic bodies to an external delta force is determined in analogy with recent experiments on granular aggregates. For the generated forces, stress and displacement fields, we find strong relative fluctuations of order one close to the source, which, however, average out readily to the classical predictions of isotropic continuum elasticity. The stress fluctuations decay (essentially) exponentially with distance from the source. Only beyond a surprisingly large distance, $b \approx 30$ interatomic distances, self-averaging dominates, and the quenched disorder becomes irrelevant for the response of an individual configuration. We argue that this self-averaging length $b$ sets also the lower wavelength bound for the applicability of classical eigenfrequency calculations.Particular attention is paid to the displacements of the source, allowing a direct measurement of the local rigidity. The algebraic correlations of these displacements demonstrate the existence of domains of slightly different rigidity without, however, revealing a characteristic length scale, at least not for the system sizes we are able to probe., Comment: 27 pages, 10 figures

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 $\phi$) 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 $\tau_+$ is set by the chain end density $\phi/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 \tau_+$ 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., Comment: 24 pages, 18 figures, submitted to EPJ E

The approach of the elastic continuum limit in small amorphous bodies formed by weakly polydisperse Lennard-Jones beads is investigated in a systematic finite-size study. We show that classical continuum elasticity breaks down when the wavelength of the sollicitation is smaller than a characteristic length of approximately 30 molecular sizes. Due to this surprisingly large effect ensembles containing up to N=40,000 particles have been required in two dimensions to yield a convincing match with the classical continuum predictions for the eigenfrequency spectrum of disk-shaped aggregates and periodic bulk systems. The existence of an effective length scale \xi is confirmed by the analysis of the (non-gaussian) noisy part of the low frequency vibrational eigenmodes. Moreover, we relate it to the {\em non-affine} part of the displacement fields under imposed elongation and shear. Similar correlations (vortices) are indeed observed on distances up to \xi~30 particle sizes., Comment: 28 pages, 13 figures, 3 tables

Structures involving solid particles of nanometric dimensions play an increasingly important role in material sciences. These structures are often characterized through the vibrational properties of their constituent particles, which can be probed by spectroscopic methods. Interpretation of such experimental data requires an extension of continuum elasticity theory down to increasingly small scales. Using numerical simulation and exact diagonalization for simple models, we show that continuum elasticity, applied to disordered system, actually breaks down below a length scale of typically 30 to 50 molecular sizes. This length scale is likely related to the one which is generally invoked to explain the peculiar vibrational properties of glassy systems., Comment: 4 pages, 5 figures, LATEX, Europhysics Letters accepted

The osmotic pressure $P$ in equilibrium polymers (EP) in good solvent is investigated by means of a three dimensional off-lattice Monte Carlo simulation. Our results compare well with real space renormalisation group theory and the osmotic compressibility $K \propto \phi \upd \phi/\upd P$ from recent light scattering study of systems of long worm-like micelles. We confirm the scaling predictions for EP based on traditional physics of quenched monodisperse polymers in the dilute and semidilute limit. Specifically, we find $P\propto \phi^{2.3}$ and, hence, $K \propto \phi^{-0.3}$ in the semidilute regime --- in agreement with both theory and experiment. At higher concentrations where the semidilute blobs become too small and hard-core interactions and packing effects become dominant, a much stronger increase % $\log(P/\phi)\approx \log(\Nav^2/\phi) \propto \phi$ is evidenced and, consequently, the compressibility decreases much more rapidly with $\phi$ than predicted from semidilute polymer theory, but again in agreement with experiment., Comment: 7 pages, 4 figures, LATEX

We investigate by means of a number of different dynamical Monte Carlo simulation methods the self-assembly of equilibrium polymers in dilute, semidilute and concentrated solutions under good-solvent conditions. In our simulations, both linear chains and closed loops compete for the monomers, expanding on earlier work in which loop formation was disallowed. Our findings show that the conformational properties of the linear chains, as well as the shape of their size distribution function, are not altered by the formation of rings. Rings only seem to deplete material from the solution available to the linear chains. In agreement with scaling theory, the rings obey an algebraic size distribution, whereas the linear chains conform to a Schultz--Zimm type of distribution in dilute solution, and to an exponentional distribution in semidilute and concentrated solution. A diagram presenting different states of aggregation, including monomer-, ring- and chain-dominated regimes, is given.