The theory of rubber elasticity

Initial Crystallization Effects in Coarse-Grained Polyethylene Systems After Uni- and Biaxial Stretching in Blow-Molding Cooling Scenarios

This study investigates the initial stage of the thermo-mechanical crystallization behavior for uni- and biaxially stretched polyethylene. The models are based on a mesoscale molecular dynamics approach. We take constraints that occur in real-life polymer processing into account, especially with respect to the blowing stage of the extrusion blow-molding process. For this purpose, we deform our systems using a wide range of stretching levels before they are quenched. We discuss the effects of the stretching procedures on the micro-mechanical state of the systems, characterized by entanglement behavior and nematic ordering of chain segments. For the cooling stage, we use two different approaches which allow for free or hindered shrinkage, respectively. During cooling, crystallization kinetics are monitored: We precisely evaluate how the interplay of chain length, temperature, local entanglements and orientation of chain segments influence crystallization behavior. Our models reveal that the main stretching direction dominates microscopic states of the different systems. We are able to show that crystallization mainly depends on the (dis-)entanglement behavior. Nematic ordering plays a secondary role.

A Ribbon Model for Nematic Polymer Networks

We present a theory of deformation of ribbons made of nematic polymer networks (NPNs). These materials exhibit properties of rubber and nematic liquid crystals, and can be activated by external stimuli of heat and light. A two-dimensional energy for a sheet of such a material has already been derived from the celebrated neo-classical energy of nematic elastomers in three space dimensions. Here, we use a dimension reduction method to obtain the appropriate energy for a ribbon from the aforementioned sheet energy. We also present an illustrative example of a rectangular NPN ribbon that undergoes in-plane serpentine deformations upon activation under an appropriate set of boundary conditions.

Effect of Nonisoprene Degradation and Naturally Occurring Network during Maturation on the Properties of Natural Rubber

It well-known that the superior performance of natural rubber (NR) compared to its synthetic counterpart mainly derives from nonisoprene components and naturally occurring network, which varies during the progress of the maturation and thereby results in technically graded rubber with different properties. However, identifying the roles of these two factors in the forming of excellent performance of NR is still a challenge as they change simultaneously during the maturation process. Here, influences of naturally occurring networking and nonisoprene degradation on the components, structures and properties of NR were systematically investigated by tailored treatments of maturation. It was found that the maturation-induced formation of natural network structure contributes to the increase in initial plastic value, Mooney viscosity and gel content for un-crosslinked NR, while the decomposition of nonisoprene components plays a dominant role in improving the mechanical properties of vulcanized NR. Stress-strain curve and Mooney-Rivlin analysis demonstrate that the biodegradation of the nonisoprene components significantly boost the vulcanization process, which significantly increases the number of chemical cross-link networks and effective cross-link density of the material, greatly improving the mechanical properties of NR vulcanizates. This resulted in the tensile strength of TSR 10CV being able to reach 22.6 MPa, which is significantly improved compared to 15.8 MPa of TSR 3CV. Evidenced by tubular model fitting, the increase in chemical cross-linking points effectively reduces the movable radius of the molecular chain under dynamic loading, making the molecular chain more difficult to move, which suppresses the entropy change under dynamic loading and consequently endows NR excellent dynamic mechanical properties. This resulted in a significant decrease in the temperature rising of TSR 10CV to 3.3 °C, while the temperature rising of TSR 3CV was still as high as 14.5 °C. As a minor factor, the naturally occurring network improves the mechanical properties of vulcanizates in the form of sacrificial bonds.

Investigation of Crystallization and Relaxation Effects in Coarse-Grained Polyethylene Systems after Uniaxial Stretching

In this study, we investigate the thermo-mechanical relaxation and crystallization behavior of polyethylene using mesoscale molecular dynamics simulations. Our models specifically mimic constraints that occur in real-life polymer processing: After strong uniaxial stretching of the melt, we quench and release the polymer chains at different loading conditions. These conditions allow for free or hindered shrinkage, respectively. We present the shrinkage and swelling behavior as well as the crystallization kinetics over up to 600 ns simulation time. We are able to precisely evaluate how the interplay of chain length, temperature, local entanglements and orientation of chain segments influences crystallization and relaxation behavior. From our models, we determine the temperature dependent crystallization rate of polyethylene, including crystallization onset temperature.

Tensile behavior of non-crosslinked networks of athermal fibers in the presence of entanglements and friction

Many biological and soft artificial materials contain a random network of non-crosslinked fibers as their main structural component. The excluded volume interactions (contact forces) at fiber contacts control the mechanical behavior of these systems. This physics has been studied extensively in compression, but little is known about the relation between network structure and its mechanical response in tension. In particular, although occasionally used conjecturally, the notion of fiber entanglements in athermal networks is not well defined, nor is it clear what role entanglements play in athermal network mechanics. The primary contribution of this work is the introduction of a measure of the degree of entanglement of a system of random athermal fibers, and the definition of its relationship with the mechanical behavior of the network. Entanglements confine the fibers during tensile loading, reduce the auxetic effect in mat-like networks, and maintain the inter-fiber contact density. In the absence of this contribution, reduction of the contact density during tensile loading due to auxeticity results in stress reduction. Entanglements stabilize the network via a tensegrity mechanism similar to that operating in woven materials and lead to network stiffening. The relation between the proposed measure of entanglements and the fiber volume fraction is defined. The effect of inter-fiber friction on the mechanics of entangled mat-like non-crosslinked fiber networks is also evaluated.

Ridge energy for thin nematic polymer networks

Minimizing the elastic free energy of a thin sheet of nematic polymer network among smooth isometric immersions is the strategy purported by the mainstream theory. In this paper, we broaden the class of admissible spontaneous deformations: we consider ridged isometric immersions, which can cause a sharp ridge in the immersed surfaces. We propose a model to compute the extra energy distributed along such ridges. This energy comes from bending; it is shown under what circumstances it scales quadratically with the sheet's thickness, falling just in between stretching and bending energies. We put our theory to the test by studying the spontaneous deformation of a disk on which a radial hedgehog was imprinted at the time of crosslinking. We predict the number of folds that develop in terms of the degree of order induced in the material by external agents (such as heat and illumination).

A blend of stretching and bending in nematic polymer networks

Nematic polymer networks are (heat and light) activable materials, which combine the features of rubber and nematic liquid crystals. When only the stretching energy of a thin sheet of nematic polymer network is minimized, the intrinsic (Gaussian) curvature of the shape it takes upon (thermal or optical) actuation is determined. This, unfortunately, produces a multitude of possible shapes, for which we need a selection criterion, which may only be provided by a correcting bending energy depending on the extrinsic curvatures of the deformed shape. The literature has so far offered approximate corrections depending on the mean curvature. In this paper, we derive the appropriate bending energy for a sheet of nematic polymer network from the celebrated neo-classical energy of nematic elastomers in three space dimensions. This task is performed via a dimension reduction based on a modified Kirchhoff-Love hypothesis, which withstands the criticism of more sophisticated analytical tools. The result is a surface elastic free-energy density where stretching and bending are blended together; they may or may not be length-separated, and should be minimized together. The extrinsic curvatures of the deformed shape not only feature in the bending energy through the mean curvature, but also through the relative orientation of the nematic director in the frame of the directions of principal curvatures.

Topological Models for Open-Knotted Protein Chains Using the Concepts of Knotoids and Bonded Knotoids

In this paper we introduce a method that offers a detailed overview of the entanglement of an open protein chain. Further, we present a purely topological model for classifying open protein chains by also taking into account any bridge involving the backbone. To this end, we implemented the concepts of planar knotoids and bonded knotoids. We show that the planar knotoids technique provides more refined information regarding the knottedness of a protein when compared to established methods in the literature. Moreover, we demonstrate that our topological model for bonded proteins is robust enough to distinguish all types of lassos in proteins.

Microscopic Dynamics and Topology of Polymer Rings Immersed in a Host Matrix of Longer Linear Polymers: Results from a Detailed Molecular Dynamics Simulation Study and Comparison with Experimental Data

We have performed molecular dynamics (MD) simulations of melt systems consisting of a small number of long ring poly(ethylene oxide) (PEO) probes immersed in a host matrix of linear PEO chains and have studied their microscopic dynamics and topology as a function of the molecular length of the host linear chains. Consistent with a recent neutron spin echo spectroscopy study (Goossen et al., Phys. Rev. Lett. 2015, 115, 148302), we have observed that the segmental dynamics of the probe ring molecules is controlled by the length of the host linear chains. In matrices of short, unentangled linear chains, the ring probes exhibit a Rouse-like dynamics, and the spectra of their dynamic structure factor resemble those in their own melt. In striking contrast, in matrices of long, entangled linear chains, their dynamics is drastically altered. The corresponding dynamic structure factor spectra exhibit a steep initial decay up to times on the order of the entanglement time τe of linear PEO at the same temperature but then they become practically time-independent approaching plateau values. The plateau values are different for different wavevectors; they also depend on the length of the host linear chains. Our results are supported by a geometric analysis of topological interactions, which reveals significant threading of all ring molecules by the linear chains. In most cases, each ring is simultaneously threaded by several linear chains. As a result, its dynamics at times longer than a few τe should be completely dictated by the release of the topological restrictions imposed by these threadings (interpenetrations). Our topological analysis did not indicate any effect of the few ring probes on the statistical properties of the network of primitive paths of the host linear chains.

Influence of molecular architecture on the entanglement network: topological analysis of linear, long- and short-chain branched polyethylene melts via Monte Carlo simulations

We present detailed results on the effect of chain branching on the topological properties of entangled polymer melts via an advanced connectivity-altering Monte Carlo (MC) algorithm. Eleven representative model linear, short-chain branched (SCB), and long-chain branched (LCB) polyethylene (PE) melts were employed, based on the total chain length and/or the longest linear chain dimension. Directly analyzing the entanglement [or the primitive path (PP)] network of the system via the Z-code, we quantified several important topological measures: (a) the PP contour length Lpp, (b) the number of entanglements Zes per chain, (c) the end-to-end length of an entanglement strand des, (d) the number of carbon atoms per entanglement strand Nes, and (e) the probability distribution for each of these quantities. The results show that the SCB polymer melts have significantly more compact overall chain conformations compared to the linear polymers, exhibiting, relative to the corresponding linear analogues, (a) ∼20% smaller values of 〈Lpp〉 (the statistical average of Lpp), (b) ∼30% smaller values of 〈Zes〉, (c) ∼20% larger values of 〈des〉, and (d) ∼50% larger values of 〈Nes〉. In contrast, despite the intrinsically smaller overall chain dimensions than those of the linear analogues, the LCB (H-shaped and A3AA3 multiarm) PE melts exhibit relatively (a) 7-8% larger values of 〈Lpp〉, (b) 6-11% larger values of 〈Zes〉 for the H-shaped melt and ∼2% smaller values of 〈Zes〉 for the A3AA3 multiarm, (c) 2-5% smaller values of 〈des〉, and (d) 7-11% smaller values of 〈Nes〉. Several interesting features were also found in the results of the probability distribution functions P for each topological measure.

Topological Constraints in Directed Polymer Melts

Polymers in a melt may be subject to topological constraints, as in the example of unlinked polymer rings. How to do statistical mechanics in the presence of such constraints remains a fundamental open problem. We study the effect of topological constraints on a melt of directed polymers, using simulations of a simple quasi-2D model. We find that fixing the global topology of the melt to be trivial changes the polymer conformations drastically. Polymers of length L wander in the transverse direction only by a distance of order (lnL)^{ζ} with ζ≃1.5. This is strongly suppressed in comparison with the Brownian L^{1/2} scaling which holds in the absence of the topological constraint. It is also much smaller than the predictions of standard heuristic approaches-in particular the L^{1/4} of a mean-field-like "array of obstacles" model-so our results present a sharp challenge to theory. Dynamics are also strongly affected by the constraints, and a tagged monomer in an infinite system performs logarithmically slow subdiffusion in the transverse direction. To cast light on the suppression of the strands' wandering, we analyze the topological complexity of subregions of the melt: the complexity is also logarithmically small, and is related to the wandering by a power law. We comment on insights the results give for 3D melts, directed and nondirected.

Long range correlations and folding angle with applications to α-helical proteins

The conformational complexity of chain-like macromolecules such as proteins and other linear polymers is much larger than that of point-like atoms and molecules. Unlike particles, chains can bend, twist, and even become knotted. Thus chains might also display a much richer phase structure. Unfortunately, it is not very easy to characterize the phase of a long chain. Essentially, the only known attribute is the radius of gyration. The way how it changes when the degree of polymerization becomes different, and how it evolves when the ambient temperature and solvent properties change, is commonly used to disclose the phase. But in any finite length chain there are corrections to scaling that complicate the detailed analysis of the phase structure. Here we introduce a quantity that we call the folding angle to identify and scrutinize the phase structure, as a complement to the radius of gyration. We argue for a mean-field level relationship between the folding angle and the scaling exponent in the radius of gyration. We then estimate the value of the folding angle in the case of crystallographic α-helical protein structures in the Protein Data Bank. We also show how the experimental value of the folding angle can be obtained computationally, using a semiclassical Born-Oppenheimer description of α-helical chiral chains.

From a melt of rings to chromosome territories: the role of topological constraints in genome folding

We review pro and contra of the hypothesis that generic polymer properties of topological constraints are behind many aspects of chromatin folding in eukaryotic cells. For that purpose, we review, first, recent theoretical and computational findings in polymer physics related to concentrated, topologically simple (unknotted and unlinked) chains or a system of chains. Second, we review recent experimental discoveries related to genome folding. Understanding in these fields is far from complete, but we show how looking at them in parallel sheds new light on both.

Identifying knots in proteins

Polypeptide chains form open knots in many proteins. How these knotted proteins fold and finding the evolutionary advantage provided by these knots are among some of the key questions currently being studied in the protein folding field. The detection and identification of protein knots are substantial challenges. Different methods and many variations of them have been employed, but they can give different results for the same protein. In the present article, we review the various knot identification algorithms and compare their relative strengths when applied to the study of knots in proteins. We show that the statistical approach based on the uniform closure method is advantageous in comparison with other methods used to characterize protein knots.

Topological versus rheological entanglement length in primitive-path analysis protocols, tube models, and slip-link models

We show that the front factor appearing in the shear modulus of a phantom network, G(ph) = (1-2/f)(ρk(B)T)/N(s), also controls the ratio of the strand length, N(s), and the number of monomers per Kuhn length of the primitive paths, N(ph)(PPKuhn), characterizing the average network conformation. In particular, N(ph)(PPKuhn) = N(s)/(1-2/f) and G(ph) = (ρk(B)T)/N(ph)(PPKuhn). Neglecting the difference between cross-links and slip-links, these results can be transferred to entangled systems and the interpretation of primitive path analysis data. In agreement with the tube model, the analogy to phantom networks suggest that the rheological entanglement length, N(e)(rheo) = (ρk(B)T)/G(e), should equal N(e)(PPKuhn). Assuming binary entanglements with f = 4 functional junctions, we expect that N(e)(rheo) should be twice as large as the topological entanglement length, N(e)(topo). These results are in good agreement with reported primitive path analysis results for model systems and a wide range of polymeric materials. Implications for tube and slip-link models are discussed.

A gaussian model for substrates of entangled cross-linked poly(ethylene glycol) in biomedical applications

Cells usually spread on a synthetic substrate through bonds between receptors and chemical groups on the substrate (ligands). Therefore, it is valuable to study the effects of the average number density of these chemical groups and the average distance between them to model and predict the cell behavior. Poly(ethylene glycol) [PEG] modified with peptide groups has been used widely in biomedical applications as a substrate material. In this study, a coarse-grained model is proposed for PEG to predict the average number density of ligands and the average distance between them. Molecular information such as initial molecular weight distribution, average molecular weight between cross-links, and average molecular weight between entanglements is used as input parameters. Based on simulation results, it is concluded that both entanglement and cross-link densities are required to create a network structure. The results suggest that an average initial molecular weight 2-3 times the average molecular weight between entanglements and a moderate cross-link density are sufficient to create a closed network structure with a high ligand density and a small average distance between them.

Topological analysis of polymeric melts: chain-length effects and fast-converging estimators for entanglement length

Primitive path analyses of entanglements are performed over a wide range of chain lengths for both bead spring and atomistic polyethylene polymer melts. Estimators for the entanglement length N_{e} which operate on results for a single chain length N are shown to produce systematic O(1/N) errors. The mathematical roots of these errors are identified as (a) treating chain ends as entanglements and (b) neglecting non-Gaussian corrections to chain and primitive path dimensions. The prefactors for the O(1/N) errors may be large; in general their magnitude depends both on the polymer model and the method used to obtain primitive paths. We propose, derive, and test new estimators which eliminate these systematic errors using information obtainable from the variation in entanglement characteristics with chain length. The new estimators produce accurate results for N_{e} from marginally entangled systems. Formulas based on direct enumeration of entanglements appear to converge faster and are simpler to apply.

Analytic expressions for the statistics of the primitive-path length in entangled polymers

An analytic expression is proposed for the primitive-path length of entangled polymer chains. The expression is derived from statistical mechanics of a chain that is a random walk with randomly scattered entanglements. The only parameters are the number of Kuhn steps in the chain and a dimensionless parameter beta that contains information about the entanglement density and Kuhn step size. The expression is found to compare very favorably with numerical results recently found from examining topological constraints in microscopic simulations. The comparison also predicts well the plateau modulus of polyethylene, suggesting that the slip-link model is a viable intermediate in the search for true ab initio rheology predictions. Since the expression is analytic, it can be used to make predictions where the simulations cannot reach, and hence is applicable for coarse graining.