pith. sign in

arxiv: 2604.03424 · v1 · submitted 2026-04-03 · ❄️ cond-mat.soft

Shear Banding in Simulations of Polymer Melts

Lucas L. Nelson , Gary S. Grest , Peter D. Olmsted This is my paper

Pith reviewed 2026-05-13 18:11 UTC · model grok-4.3

classification ❄️ cond-mat.soft
keywords shear bandingpolymer meltsentanglement numberconvected constraint releaseRolie-Poly modeltube dynamicsmolecular dynamicsrheology
0
0 comments X

The pith

A model coupling Rolie-Poly tube dynamics to CCR predicts shear banding only when polymer entanglement number Z exceeds a β-dependent critical value.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper compares molecular dynamics simulations of polymer melts in shear flow against a model that combines Rolie-Poly-like tube dynamics with entanglement evolution controlled by convected constraint release. The model shows shear banding emerges only when the equilibrium entanglement number Z surpasses a critical threshold Z_c that depends on the CCR strength parameter β, which is extracted from the simulations via chain bending stiffness. Simulations match the model predictions at a semi-quantitative level. The work offers a route to identify which real polymers are likely to exhibit shear banding during processing flows.

Core claim

The central claim is that a model coupling Rolie-Poly-like tube dynamics to CCR-mediated entanglement dynamics predicts shear banding for equilibrium entanglement numbers Z that exceed a critical value Z_c depending on the CCR parameter β, and that the simulations agree with this prediction in a semi-quantitative manner, with observed deviations attributed to approximations in the model and the limitations of tube descriptions for partially disentangled liquids under strong flow.

What carries the argument

The Rolie-Poly-like tube model coupled to convected constraint release (CCR) entanglement dynamics, with the parameter β controlling CCR strength and obtained from the bending stiffness k_θ measured in the simulations.

If this is right

  • Shear banding occurs only for Z greater than the β-dependent Z_c.
  • Simulations and model predictions agree semi-quantitatively.
  • Deviations are due to model approximations and tube-model limits in strong flow.
  • The results indicate which physical polymers are candidates for shear banding.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Processing flows for polymers could be tuned by chain length and stiffness to cross or avoid the Z_c threshold and thereby control banding.
  • Similar thresholds might appear in other inhomogeneous flows such as extension or mixed shear-extension.
  • Direct comparison with rheological experiments on model polymers with independently measured Z and β would provide a sharper test.

Load-bearing premise

The Rolie-Poly-like tube dynamics coupled to CCR remain adequate for describing entanglement evolution in strong shear, even though the paper notes their inadequacy for partially-disentangled liquids.

What would settle it

A direct simulation or experiment on a polymer melt with known Z and β showing the absence of shear banding well above the predicted Z_c or the presence of banding well below it would falsify the central prediction.

Figures

Figures reproduced from arXiv: 2604.03424 by Gary S. Grest, Lucas L. Nelson, Peter D. Olmsted.

Figure 2
Figure 2. Figure 2: FIG. 2. Boundaries [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 1
Figure 1. Figure 1: FIG. 1. Stress-strain rate curves (solid lines) and relative en [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Disentanglement as a function of Rouse Weissenberg [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: shows how the inferred value of β depends on chain stiffness, quantified by the ratio of the Kuhn length bK to the packing length p. Dolata et al. 56 suggested that β should be a simple function of bK/p, which they ap￾proximated by a power law β ∼ (bK/p) 1.9 . We find that a slightly less steep power law β ∝ (bK/p) 1.6 agrees with our KG simulations and the results from Mohagheghi and Khomami. However, Dol… view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. Steady-state velocity profiles of shear-banded [PITH_FULL_IMAGE:figures/full_fig_p005_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. Material properties of simulations that do (solid [PITH_FULL_IMAGE:figures/full_fig_p006_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7. Range of reptation Weissenberg numbers for which [PITH_FULL_IMAGE:figures/full_fig_p007_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8. Normal stress ratios vs reptation Weissenberg num [PITH_FULL_IMAGE:figures/full_fig_p008_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: FIG. 9. Range of reptation Weissenberg numbers for which [PITH_FULL_IMAGE:figures/full_fig_p010_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: FIG. 10. Effect of varying [PITH_FULL_IMAGE:figures/full_fig_p011_10.png] view at source ↗
Figure 12
Figure 12. Figure 12: FIG. 12. Monotonicity boundaries in [PITH_FULL_IMAGE:figures/full_fig_p012_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: FIG. 13. Steady-state velocity profiles with different shear [PITH_FULL_IMAGE:figures/full_fig_p012_13.png] view at source ↗
read the original abstract

Results from numerical simulations of polymers under shear flow are compared with predictions for shear banding based on a model coupling Rolie-Poly-like tube dynamics to the entanglement dynamics mediated by Convected Constraint Release (CCR). CCR is controlled by a parameter $\beta$, whose dependence on the bending stiffness $k_{\theta}$ is calculated from the simulations. The model predicts shear banding for polymers whose equilibrium entanglement number $Z$ exceeds a critical value $Z_{c}$ that depends on $\beta$. The simulations are in semi-quantitative agreement with the model, with deviations that are attributed to approximations inherent in the model, and the inadequacy of tube models to describe partially-disentangled liquids in strong flow. These results may help determine which physical polymers could undergo shear banding.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

3 major / 2 minor

Summary. The paper compares molecular dynamics simulations of entangled polymer melts under shear to predictions from a constitutive model that couples Rolie-Poly-like tube dynamics to Convected Constraint Release (CCR) with a single parameter β. β is determined from the simulated dependence on bending stiffness k_θ. The model predicts the onset of shear banding when the equilibrium entanglement number Z exceeds a critical value Z_c(β). The simulations are reported to be in semi-quantitative agreement with these predictions, with deviations ascribed to model approximations and the known limitations of tube models for partially disentangled liquids in strong flow. The results are suggested to help identify which physical polymers may exhibit shear banding.

Significance. If the semi-quantitative agreement can be placed on a firmer footing, the work would usefully connect simulation-derived parameters to a continuum criterion for shear banding in entangled polymers. The explicit mapping of β(k_θ) and the identification of a Z_c threshold constitute concrete, testable outputs that could guide both further modeling and experimental design. However, the acknowledged breakdown of the underlying tube dynamics precisely in the strong-flow regime of interest reduces the confirmatory power of the comparison.

major comments (3)
  1. [Abstract] Abstract: β is extracted from the same set of simulations whose banding behavior is then compared against the model predictions. This creates a circularity in which the Z_c(β) threshold is calibrated, at least in part, against the very data used to test it, weakening the independence of the validation.
  2. [Abstract] Abstract: the paper states that deviations are due in part to 'the inadequacy of tube models to describe partially-disentangled liquids in strong flow.' Because shear banding onset and steady-state banding occur precisely in this regime, the self-acknowledged breakdown of the Rolie-Poly + CCR equations undermines the claim that the observed semi-quantitative agreement constitutes supportive evidence for the model.
  3. [Results] Results section (comparison of Z_c and simulated banding): no quantitative error bars, full data tables, or statistical measures of agreement are provided, making it impossible to assess whether the reported semi-quantitative match is robust or merely consistent within large uncertainties.
minor comments (2)
  1. [Model and Methods] The manuscript should clarify whether the reported Z_c(β) curve is obtained by solving the constitutive equations analytically or numerically, and should state the precise criterion used to identify the onset of banding in both the model and the simulations.
  2. [Figures] Figure captions and text should explicitly note the range of shear rates and entanglement numbers explored, so that readers can judge how far into the strong-flow regime the comparisons extend.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for their constructive comments, which have helped us clarify the presentation and strengthen the manuscript. We address each major point below and have revised the text accordingly where changes were warranted.

read point-by-point responses

  1. Referee: [Abstract] Abstract: β is extracted from the same set of simulations whose banding behavior is then compared against the model predictions. This creates a circularity in which the Z_c(β) threshold is calibrated, at least in part, against the very data used to test it, weakening the independence of the validation.

    Authors: We agree that the shared use of simulation data requires clearer separation. β is obtained from the k_θ dependence of the constraint-release rate measured in the linear-response regime of the same trajectories, while the banding comparison uses the nonlinear shear response at high Weissenberg number. The Z_c(β) curve is generated entirely from the constitutive model once β(k_θ) is fixed; it is not refitted to the banding data. We have added an explicit paragraph in the revised Methods and Results sections that isolates the linear-regime extraction step from the nonlinear banding test, thereby reducing the appearance of circularity. revision: partial

  2. Referee: [Abstract] Abstract: the paper states that deviations are due in part to 'the inadequacy of tube models to describe partially-disentangled liquids in strong flow.' Because shear banding onset and steady-state banding occur precisely in this regime, the self-acknowledged breakdown of the Rolie-Poly + CCR equations undermines the claim that the observed semi-quantitative agreement constitutes supportive evidence for the model.

    Authors: The referee correctly highlights a genuine limitation of the tube-model framework. We have revised the abstract and the concluding discussion to state that the semi-quantitative match is offered only as evidence that the CCR-augmented Rolie-Poly model captures the correct trend with Z and β, while acknowledging that quantitative discrepancies are expected once the tube picture itself fails at high strain rates. The revised text now frames the comparison as a useful diagnostic rather than a definitive validation. revision: yes

  3. Referee: [Results] Results section (comparison of Z_c and simulated banding): no quantitative error bars, full data tables, or statistical measures of agreement are provided, making it impossible to assess whether the reported semi-quantitative match is robust or merely consistent within large uncertainties.

    Authors: We accept this criticism. The revised Results section now includes error bars on all simulated banding-onset points (obtained from three independent runs), a supplementary table listing Z_c(model) versus Z_c(sim) for each k_θ together with the relative deviation, and a brief statement of the maximum observed deviation (approximately 15 %). These additions allow readers to judge the robustness of the reported agreement directly. revision: yes

Circularity Check

1 steps flagged

β fitted from simulations then used to predict Z_c(β) for comparison against same data

specific steps
  1. fitted input called prediction [Abstract]
    "CCR is controlled by a parameter β, whose dependence on the bending stiffness k_θ is calculated from the simulations. The model predicts shear banding for polymers whose equilibrium entanglement number Z exceeds a critical value Z_c that depends on β. The simulations are in semi-quantitative agreement with the model"

    β is extracted from the identical set of simulations whose banding behavior is later compared to the model's Z_c(β) prediction; the reported agreement therefore incorporates the fitted parameter by construction.

full rationale

The paper calculates β directly from the simulations and inserts it into the Rolie-Poly+CCR model to obtain a Z_c(β) threshold; the same simulations are then reported as showing semi-quantitative agreement with that threshold. This constitutes a fitted-input-called-prediction pattern. The underlying constitutive equations are independent of the present data, so the circularity is partial rather than total. The paper's own acknowledgment of tube-model breakdown in strong flow is noted but does not itself create an additional circularity step.

Axiom & Free-Parameter Ledger

1 free parameters · 2 axioms · 0 invented entities

The central claim rests on standard tube-model assumptions for entangled polymers plus the specific CCR mechanism; β is not a free parameter but is computed from the simulations themselves.

free parameters (1)
  • β
    CCR strength parameter whose value and dependence on bending stiffness k_θ are extracted from the simulations rather than chosen independently.
axioms (2)
  • domain assumption Rolie-Poly-like tube dynamics govern chain orientation and stretch under flow
    Invoked to couple tube relaxation to the CCR term in the model.
  • domain assumption Convected constraint release is the dominant entanglement-loss mechanism in shear
    Central modeling choice that controls the β dependence.

pith-pipeline@v0.9.0 · 5424 in / 1388 out tokens · 44413 ms · 2026-05-13T18:11:24.700635+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

5 extracted references · 5 canonical work pages · 1 internal anchor

  1. [1]

    The DO arXiv:2604.03424v1 [cond-mat.soft] 3 Apr 2026 2 model predicts banding for sufficiently-entangled poly- mersZ > Z c(β, α), where the critical entanglement num- berZ c is larger for larger CCR parameterβand smaller relaxation anisotropyα(defined below). The region of banding found in the simulations agrees surprisingly well with the model, given tha...

    work page internal anchor Pith review Pith/arXiv arXiv 2026

  2. [2]

    1a demonstrates the mechanism: stronger CCR increases disentanglement (smallerν) when ˙γ∼1/τ R, which decreases polymer alignment and increasesσ xy

    Fig. 1a demonstrates the mechanism: stronger CCR increases disentanglement (smallerν) when ˙γ∼1/τ R, which decreases polymer alignment and increasesσ xy. This value of ˙γis typically near the local stress minimum, so a sufficient increase of σxy there can eliminate the nonmonotonicity. Significantly, a larger relaxation anisotropyαpromotes banding. As evi...

    work page 2000

  3. [3]

    Lines are least- squares fits of the DO model to each of these systems, with the fit value ofβgiven

    Disentanglement as a function of Rouse Weissenberg number for four KG systems (triangles) and simulations by Mohagheghi and Khomami 47 (pentagons). Lines are least- squares fits of the DO model to each of these systems, with the fit value ofβgiven. Hollow triangles are data from Dolata et al. 56 computed using the Z1 code on a smaller system sus- ceptible...

    work page 2000

  4. [4]

    would be more likely to exhibit banding. ACKNOWLEDGMENTS We thank Joe Peterson, Ben Dolata, Thomas O’Connor, Marco Galvani Cunha, and Shi-Qing Wang for discussions; and Martin Kr¨ oger for generous help with and advice about the Z1+ code. PDO and LLN thank the National Science Foundation (DMREF-2118769), and PDO thanks Georgetown University and the Ives F...

    work page 2000

  5. [5]

    Reptation of a polymer chain in the presence of fixed obstacles,

    Steady-state velocity profiles with different shear histories for a system of 3200 chains of 400 beads withk θ = 2.0. In each part, two steady-state profiles are shown for the same shear rate (γ= 10 −5 and 10 −6τ −1 starting either from rest (red, hollow symbols) or from a steady-state configuration at a different shear rates (10 −3 and 10 −5τ −1) (black,...

    work page doi:10.1063/1.2192775 1988