[Paper Review] Topological Analysis of Linear Polymer Melts

Christos Tzoumanekas, Doros N. Theodorou|arXiv (Cornell University)|Feb 23, 2006

Advanced Theoretical and Applied Studies in Material Sciences and Geometry3 citations

TL;DR

This paper introduces the CReTA algorithm to extract primitive paths from atomistic polymer melts, enabling topological analysis of entanglement networks in polyethylene and cis-1,4-polybutadiene. The study reveals a universal scaling behavior in entanglement distributions, demonstrating that both melts exhibit statistically identical topological structure when scaled by mesh length, challenging the tube model’s assumption of uniform mesh size.

ABSTRACT

We introduce an algorithm for the reduction of computer generated atomistic polymer samples to networks of primitive paths. By examining network ensembles of Polyethylene and cis-1,4 Polybutadiene melts, we quantify the underlying topologies through the radial distribution function of entanglements and the distribution of the number of monomers between entanglements. A suitable scaling of acquired data leads to a unifying microscopic topological description of both melts.

Motivation & Objective

• To develop a microscopic, algorithmic method for extracting primitive paths from atomistic polymer configurations to reveal topological constraints.
• To quantify the statistical distribution of entanglements in flexible polymer melts, focusing on the number of monomers between entanglements.
• To test whether topological properties of chemically distinct polymers (PE and PB) can be unified under a single scaling framework.
• To validate the topological network representation against experimental data, establishing its physical relevance.
• To provide a foundation for constructing mesoscopic simulation models of polymer melts and glasses based on microscopically derived entanglement statistics.

Proposed method

• Application of the CReTA (Contour Reduction Topological Analysis) algorithm to reduce atomistic chain configurations to networks of primitive paths with fixed chain ends.
• Use of random aligning string moves and hardcore repulsion to prevent chain crossings during contour shortening, preserving topological constraints.
• Treatment of CHx monomers as united-atom hard spheres of diameter σ ≈ 3.5 Å, with chains modeled as fused spheres to avoid unphysical overlaps.
• Implementation of periodic boundary conditions and iterative reduction of unentangled loops to straight segments connected at nodal entanglements.
• Calculation of radial distribution functions (RDFs) of entanglements, distinguishing between spatial neighbors (RDF-0), topological neighbors (RDF-1, RDF-2), and total RDFs.
• Scaling of data using the average mesh length 〈d_ES〉 to test for universality across different polymer chemistries.

Experimental results

Research questions

• RQ1 Can a consistent, microscopic topological description be derived for linear polymer melts using primitive path networks?
• RQ2 Do the distributions of monomers between entanglements in polyethylene and cis-1,4-polybutadiene exhibit universal scaling when properly normalized?
• RQ3 What is the spatial correlation structure of entanglements in polymer melts, and how do topological versus spatial neighbors differ in their radial distribution?
• RQ4 To what extent do the observed entanglement distributions deviate from the assumptions of the classical tube model, particularly regarding uniform mesh size?
• RQ5 Can the stochastic process underlying entanglement placement along chains be modeled as a renewal process with effective repulsion between successive entanglements?

Key findings

• The distribution P(n) of monomers between entanglements exhibits a broad exponential tail, indicating significant fluctuations in mesh size and challenging the tube model’s assumption of uniformity.
• The downturn in P(n) at small n (n→0) suggests an effective repulsion between entanglements in monomer sequence space, implying a minimum separation requirement for entanglement formation.
• The radial distribution function (RDF) of entanglements for both PE and PB melts collapses onto a single curve when scaled by 〈d_ES〉, indicating universal spatial network correlations.
• RDF-0 (spatial neighbors) is below RDF-1 (topological neighbors) for r′ < 0.7, indicating that topological neighbors persist in close proximity, while spatial neighbors are excluded from the immediate vicinity of a node.
• The total RDF shows short-range correlations at r′ ≲ 0.8 and approaches ideal gas behavior at r′ ≳ 0.8, suggesting that long-range entanglement structure is dilute and uncorrelated.
• The exponential tail of P(n) is consistent with a stochastic model involving two uncorrelated Poisson processes, where one creates unentangled sequences and the other places entanglements, supporting a renewal process interpretation.