Self-similar chain conformations in polymer gels

We use molecular dynamics simulations to study the swelling of randomly end-cross-linked polymer networks in good solvent conditions. We find that the equilibrium degree of swelling saturates at Q_eq = N_e**(3/5) for mean strand lengths N_s exceeding the melt entanglement length N_e. The internal structure of the network strands in the swollen state is characterized by a new exponent nu=0.72. Our findings are in contradiction to de Gennes' c*-theorem, which predicts Q_eq proportional N_s**(4/5) and nu=0.588. We present a simple Flory argument for a self-similar structure of mutually interpenetrating network strands, which yields nu=7/10 and otherwise recovers the classical Flory-Rehner theory. In particular, Q_eq = N_e**(3/5), if N_e is used as effective strand length.Comment: 4 pages, RevTex, 3 Figure

Topological interactions in systems of mutually interlinked polymer rings

The topological interaction arising in interlinked polymeric rings such as DNA catenanes is considered. More specifically, the free energy for a pair of linked random walk rings is derived where the distance $R$ between two segments each of which is part of a different ring is kept constant. The topology conservation is imposed by the Gauss invariant. A previous approach (M.Otto, T.A. Vilgis, Phys.Rev.Lett. {\bf 80}, 881 (1998)) to the problem is refined in several ways. It is confirmed, that asymptotically, i.e. for large $R\gg R_G$ where $R_G$ is average size of single random walk ring, the effective topological interaction (free energy) scales $\propto R^4$.Comment: 16 pages, 3 figur

Variational bounds for the shear viscosity of gelling melts

We study shear stress relaxation for a gelling melt of randomly crosslinked, interacting monomers. We derive a lower bound for the static shear viscosity $\eta$, which implies that it diverges algebraically with a critical exponent $k\ge 2\nu-\beta$. Here, $\nu$ and $\beta$ are the critical exponents of percolation theory for the correlation length and the gel fraction. In particular, the divergence is stronger than in the Rouse model, proving the relevance of excluded-volume interactions for the dynamic critical behaviour at the gel transition. Precisely at the critical point, our exact results imply a Mark-Houwink relation for the shear viscosity of isolated clusters of fixed size.Comment: 5 pages; CHANGES: typos corrected, some references added; version as publishe

Tube Models for Rubber-Elastic Systems

In the first part of the paper we show that the constraining potentials introduced to mimic entanglement effects in Edwards' tube model and Flory's constrained junction model are diagonal in the generalized Rouse modes of the corresponding phantom network. As a consequence, both models can formally be solved exactly for arbitrary connectivity using the recently introduced constrained mode model. In the second part, we solve a double tube model for the confinement of long paths in polymer networks which is partially due to crosslinking and partially due to entanglements. Our model describes a non-trivial crossover between the Warner-Edwards and the Heinrich-Straube tube models. We present results for the macroscopic elastic properties as well as for the microscopic deformations including structure factors.Comment: 15 pages, 8 figures, Macromolecules in pres

Dynamic charge density correlation function in weakly charged polyampholyte globules

We study solutions of statistically neutral polyampholyte chains containing a large fraction of neutral monomers. It is known that, even if the quality of the solvent with respect to the neutral monomers is good, a long chain will collapse into a globule. For weakly charged chains, the interior of this globule is semi-dilute. This paper considers mainly theta-solvents, and we calculate the dynamic charge density correlation function g(k,t) in the interior of the globules, using the quadratic approximation to the Martin-Siggia-Rose generating functional. It is convenient to express the results in terms of dimensionless space and time variables. Let R be the blob size, and let T be the characteristic time scale at the blob level. Define the dimensionless wave vector q = R k, and the dimensionless time s = t/T. We find that for q<1, corresponding to length scales larger than the blob size, the charge density fluctuations relax according to g(q,s) = q^2(1-s^(1/2)) at short times s < 1, and according to g(q,s) = q^2 s^(-1/2) at intermediate times 1 < s 0.1, where entanglements are unimportant.Comment: 12 pages RevTex, 1 figure ps, PACS 61.25.Hq, reason replacement: Expression for dynamic corr. function g(k,t) in old version was incorrect (though expression for Fourier transform g(k,w) was correct, so the major part of the calculation remains.) Also major textual chang

Statistical mechanics of triangulated ribbons

We use computer simulations and scaling arguments to investigate statistical and structural properties of a semiflexible ribbon composed of isosceles triangles. We study two different models, one where the bending energy is calculated from the angles between the normal vectors of adjacent triangles, the second where the edges are viewed as semiflexible polymers so that the bending energy is related to the angles between the tangent vectors of next-nearest neighbor triangles. The first model can be solved exactly whereas the second is more involved. It was recently introduced by Liverpool and Golestanian Phys.Rev.Lett. 80, 405 (1998), Phys.Rev.E 62, 5488 (2000) as a model for double-stranded biopolymers such as DNA. Comparing observables such as the autocorrelation functions of the tangent vectors and the bond-director field, the probability distribution functions of the end-to-end distance, and the mean squared twist we confirm the existence of local twist correlation, but find no indications for other predicted features such as twist-stretch coupling, kinks, or oscillations in the autocorrelation function of the bond-director field.Comment: 10 pages, 13 figures. submitted to PRE, revised versio

Screening of Hydrodynamic Interactions in Semidilute Polymer Solutions: A Computer Simulation Study

We study single-chain motion in semidilute solutions of polymers of length N = 1000 with excluded-volume and hydrodynamic interactions by a novel algorithm. The crossover length of the transition from Zimm (short lengths and times) to Rouse dynamics (larger scales) is proportional to the static screening length. The crossover time is the corresponding Zimm time. Our data indicate Zimm behavior at large lengths but short times. There is no hydrodynamic screening until the chains feel constraints, after which they resist the flow: "Incomplete screening" occurs in the time domain.Comment: 3 figure

Variational theory for a single polyelectrolyte chain revisited

We reconsider the electrostatic contribution to the persistence length, $\ell_e$, of a single, infinitely long charged polymer in the presence of screening. A Gaussian variational method is employed, taking $\ell_e$ as the only variational parameter. For weakly charged and flexible chains, crumpling occurs at small length scales because conformational fluctuations overcome electrostatic repulsion. The electrostatic persistence length depends on the square of the screening length, $\ell_e\sim\kappa^{-2}$, as first argued by Khokhlov and Khachaturian by applying the Odijk-Skolnick-Fixman (OSF) theory to a string of crumpled blobs. We compare our approach to previous theoretical works (including variational formulations) and show that the result $\ell_e\sim\kappa^{-1}$ found by several authors comes from the improper use of a cutoff at small length scales. For highly charged and stiff chains, crumpling does not occur; here we recover the OSF result and validate the perturbative calculation for slightly bent rods.Comment: 11 pages, 6 figure

DNA uptake into nuclei: Numerical and analytical results

The dynamics of polymer translocation through a pore has been the subject of recent theoretical and experimental works. We have considered theoretical estimates and performed computer simulations to understand the mechanism of DNA uptake into the cell nucleus, a phenomenon experimentally investigated by attaching a small bead to the free end of the double helix and pulling this bead with the help of an optical trap. The experiments show that the uptake is monotonous and slows down when the remaining DNA segment becomes very short. Numerical and analytical studies of the entropic repulsion between the DNA filament and the membrane wall suggest a new interpretation of the experimental observations. Our results indicate that the repulsion monotonically decreases as the uptake progresses. Thus, the DNA is pulled in (i) either by a small force of unknown origin, and then the slowing down can be interpreted only statistically; (ii) or by a strong but slow ratchet mechanism, which would naturally explain the observed monotonicity, but then the slowing down requires additional explanations. Only further experiments can unambiguously distinguish between these two mechanisms.Comment: 12 pages, 6 figures, submitted to J. Phys. Cond. Ma

Force-Extension Relation and Plateau Modulus for Wormlike Chains

We derive the linear force-extension relation for a wormlike chain of arbitrary stiffness including entropy elasticity, bending and thermodynamic buckling. From this we infer the plateau modulus $G^0$ of an isotropic entangled solution of wormlike chains. The entanglement length $L_e$ is expressed in terms of the characteristic network parameters for three different scaling regimes in the entangled phase. The entanglement transition and the concentration dependence of $G^0$ are analyzed. Finally we compare our findings with experimental data.Comment: 5 pages, 1 eps-figure, to appear in PR