Pith. sign in

REVIEW 4 major objections 4 minor 57 references

Mechanical Degradation of Unentangled Polymer Melts under Uniaxial Extensional Flow

T0 review · 4 major / 4 minor · reviewed 2026-08-07 · deepseek-v4-flash

Pith's one-line read This paper argues that unentangled polymer melts under uniaxial extension show five viscosity stages, ending in thinning from chain scission, and that bond-fracture kinetics turns non-first-order once bonds break before chains stretch.

desk verdict Careful NEMD study adding a scission-driven fifth thinning stage and a rate-dependent fracture-kinetics switch, but the scaling exponent is inconsistently reported and the atomic-SLLOD assumption needs a control before the kinetics are quantitative. read the letter →

arxiv 2506.03945 v1 pith:7KVWSO5B submitted 2025-06-04 cond-mat.soft cond-mat.other

classification cond-mat.softcond-mat.other
keywords unentangledpolymermeltuniaxialextensionalflownon-equilibriummoleculardynamicschainscissionfracturekineticsviscositythickeningandthinningquarticbond-breakingpotentialmechanochemistry
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

This paper asks what happens to unentangled polymer melts when strong uniaxial stretching is forceful enough to snap the chains themselves, not just deform them. Using non-equilibrium molecular dynamics with a bond potential that can irreversibly break, it shows that as the extension rate rises the melt's steady-state viscosity passes through five stages: Newtonian, thickening, thinning, thickening, and finally a second thinning caused by chain scission. The thickening-thinning-thickening sequence is tied to how chains orient, stretch, and then stretch their bonds, and the final thinning appears only when bonds are allowed to break. The paper also claims that fracture kinetics is first-order at moderately high rates, because chains stretch before they snap, but stops being first-order once bonds break before chains have time to stretch. A sympathetic reader would care because these are exactly the flows encountered in polymer processing, where chain scission or mechanophore activation degrades the manufactured material.

What carries the argument

The load-bearing object is the quartic bond-breaking potential that replaces the FENE spring in the coarse-grained bead-spring model, allowing each backbone bond to rupture irreversibly when stretched past a cutoff, with a parameter $B_2$ tuning the bond strength. Around it, the argument uses the atomic SLLOD equations of motion with a thermostat to impose uniaxial extensional flow, combined with generalized Kraynik-Reinelt boundary conditions (periodic remapping that keeps the simulation box from collapsing under steady extension), and monitors three structural observables, the chain orientation order parameter, the chain stretch ratio, and the average bond force, plus the intact-bond count $N_{\mathrm{bond}}$. These quantities connect each viscosity stage to a chain-conformation mechanism: alignment in the first thickening, friction suppression in the thinning, bond stretching in the second thickening, and chain shortening in the final thinning. The intact-bond count versus time is what reveals the change in fracture reaction order.

What would settle it

Run the identical uniaxial-extension simulation with the flow equations and thermostat applied at the molecular center-of-mass level instead of to individual beads, and check whether the five viscosity stages and the switch from first-order to non-first-order bond-fracture kinetics survive; if they shift or disappear, the reported stages and reaction orders are artifacts of the atomic-level forcing.

Watch

Extended reading notes

Core claim

On its own terms, the paper's central claim is that mechanical degradation of unentangled melts under uniaxial extensional flow is a staged, rate-dependent process with a distinctive kinetic signature. Below a critical extension rate the melt behaves as a Newtonian fluid; increasing the rate first thickens the flow as chains align and partially extend, then thins it as aligned chains slide past one another with reduced friction, then thickens again as bonds themselves are stretched toward their limit. When a quartic bond potential replaces the unbreakable FENE potential, a further rise in extension rate triggers irreversible bond fracture, shortening chains and producing a final viscosity decrease. The discovery emphasized in the abstract is about reaction order: fracture is first-order in the regime where chains are highly stretched before they snap, with the logarithm of the intact-bond count falling linearly in time, and it departs from first-order behavior when the extension rate is so high that bonds fracture before chains stretch. The authors read this as evidence that the kinetics of flow-induced degradation is controlled by the competition between chain stretching and bond rupture.

Load-bearing premise

The load-bearing premise is that forcing the flow and the thermostat directly on individual monomer beads, rather than on whole molecules, barely distorts the dense melt; if that premise fails, the computed viscosity stages and scission kinetics would be unreliable.

Editorial extensions

If this is right

  • Below a rate-dependent critical extension rate the melt is safe from degradation; above it, chain scission shortens the chains and permanently lowers the melt's viscosity.
  • Longer chains fracture at lower extension rates, so higher-molecular-weight melts are the most vulnerable to flow-induced damage during processing.
  • Because fracture is first-order only when stretching precedes breaking, the molecular-weight distribution left after a strong flow will look different in the moderate-rate and very-high-rate regimes.
  • The same stretch-before-break logic should govern mechanophore activation, meaning flow conditions could in principle be chosen to either avoid or deliberately trigger mechanochemical reactions.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • A natural dimensionless criterion, the ratio of the chain-stretching time to the bond-fracture time, might collapse the first-order/non-first-order crossover onto a single curve for all chain lengths and bond strengths; the paper does not propose such a master parameter.
  • Because scission shortens chains and shorter chains resist further scission, the model implies self-limiting degradation: at a fixed extension rate the melt evolves toward a polydisperse mixture whose longest chains sit just below the fracture threshold, so the final molecular-weight distribution could be predicted from the initial one.
  • For real processing flows that mix shear and extension, such as injection-molding gates and electrospinning jets, the strain-rate history should matter more than the steady-state rate, since the simulations show fracture initiates during the transient stretching stage.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

4 major / 4 minor

Summary. This paper reports nonequilibrium molecular dynamics (NEMD) simulations of unentangled bead-spring polymer melts under uniaxial extensional flow (UEF), comparing a non-breakable FENE bond potential with a breakable QUARTIC potential. The steady-state extensional viscosity is mapped over roughly four decades of extension rate for chain lengths N = 20–80, revealing five stages: Newtonian, thickening, thinning, thickening, and a final thinning attributed to bond scission. The authors connect these stages to the steady-state orientation order parameter, chain stretch ratio, and average bond force. They further report the critical fracture rate as a function of chain length, an Arrhenius-like rate constant for bond fracture in the first-order regime, and a transition to non-first-order kinetics at very high extension rates.

Significance. The main value of the paper is its systematic computational survey of flow-induced chain scission in unentangled melts and the connection between rheological stages and molecular conformations. If the claims hold, the identification of a fifth thinning stage caused by bond breaking, and the distinction between first-order and non-first-order scission kinetics, are useful for interpreting mechanophore activation and polymer degradation in processing flows. The study is strengthened by the use of three independent initial configurations, time-step convergence checks, multiple chain lengths and bond strengths, and comparison with an earlier molecular-SLLOD simulation for intact chains. The central quantitative relations, however, need correction and additional validation before the conclusions can be fully accepted.

major comments (4)
  1. [Section 3.3, Figure 6, Conclusions] In Section 3.3 and Figure 6 the authors report a fitted scaling relation ε̇_c ~ N^(-0.6) and then state that 'Such a scaling relation matches that of transient extensional flow ε̇_c ~ N^(-1).' These two exponents are not consistent, so the claim as written is internally contradictory. The same -0.6 exponent is repeated in the Conclusions. Please correct either the fitted exponent or the stated comparison, and report the fit uncertainty. If -0.6 is indeed the fitted value, discuss explicitly why it differs from the -1 theoretical prediction for transient extensional flow.
  2. [Section 2, Eq. (4)-(5)] The manuscript acknowledges that applying atomic SLLOD equations and atomic thermostats to monomer beads 'can induce intramolecular stresses' and can misinterpret angular streaming as thermal energy, but asserts that in a dense polymer this 'should be minimal' without a quantitative test. Because the reported Stage II thickening and the scission kinetics depend on intramolecular bond tension and chain stretch, please provide a control simulation using molecular SLLOD and a molecular thermostat for intact chains (where molecular identity is fixed) to quantify the artifact, or cite a direct validation for this model. This is needed to establish that the viscosity stages and the first-order/non-first-order transition are not distorted by the flow-driving algorithm.
  3. [Section 3.3, Figure 8] The Arrhenius relation is written as k_f ~ exp[-(E_f - f ε̇)/kT], but E_f, f, and T are not defined in terms of the model, and the fitted dashed lines in Figure 8 are not accompanied by the fitted parameters or their uncertainties. Please state the exact fitting function, the fitted coefficients, and the chain-length dependence so the claim is falsifiable. Also clarify what f represents physically and whether it can be related to the maximum force of the QUARTIC bond.
  4. [Section 3.3, Figure 7] The classification of bond-fracture kinetics as first-order versus non-first-order is made by visual inspection of linearity in ln(N_bond) versus t. Please provide a quantitative criterion, such as an R² threshold, a comparison with a stretched-exponential fit, or a statistical test, to support the claimed transition and its location in extension rate.
minor comments (4)
  1. [Abstract] The phrase 'shear thickening-thinning-thickening stages' should be 'extensional thickening-thinning-thickening stages' because the flow studied is uniaxial extension, not shear; using 'shear' may mislead readers.
  2. [Figure 8 caption] The fitted scaling relation in the caption is written as k_f~e^(-(a-b ε̇)); please define the constants a and b or use the same symbols as in the main text.
  3. [Section 3.3, paragraph 2] The sentence 'fewer polymer chains are strongly stretched with lower λc before bond fracture' is ambiguous; please rephrase to clarify that the average chain stretch at fracture decreases with increasing extension rate.
  4. [Figure 7 caption] The caption contains a duplicated word: 'show show results at higher extension rate'; please correct this typo.

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity: all reported stages, scaling relations, and kinetics emerge from simulation outputs rather than being imposed by the model inputs or fits.

full rationale

This is an empirical NEMD simulation study. The flow stages (Newtonian, thickening, thinning, thickening, thinning) are read off steady-state viscosities computed from Eq. (6), and the chain-conformation descriptors (stretch ratio, order parameter, bond force) are independently measured from the simulated trajectories. The scaling relation eps_dot_c ~ N^{-1.6} in Figure 6 and the Arrhenius form k_f ~ exp(-(E_f - f eps_dot)/k_B T) in Figure 8 are fits to the simulated fracture data, not inputs that define those outputs. The first-order to non-first-order kinetic transition is diagnosed from the shape of ln(N_bond) versus time curves, which is an emergent property of the breakable-bond model. The only caveat the authors flag, the use of atomic SLLOD and atomic thermostats, is a modeling assumption about possible artifacts and is explicitly acknowledged in the Computational Methods section; it is a correctness/robustness concern, not a circularity. Citations to prior work are used for model validation or comparison (e.g., Ref. 54 for Stage IV thickening, Refs. 23, 50, 53 for Stages II/III behavior), and the self-citation of the SLLOD equations (Ref. 42) is a standard equations-of-motion reference, not a load-bearing derivation of the results. Therefore no circular step is present.

Assumptions & free parameters 3 free parameters · 6 assumptions · 0 invented entities

The model potentials (WCA, FENE, QUARTIC) and their parameters are taken from prior literature (refs 26,29,38) and are standard coarse-grained models. The key fitted quantities are scaling exponents and kinetic parameters extracted from the simulation output. No new physical entities are postulated.

free parameters (3)
  • B2 (quartic bond strength) = 1.25, 1.55, 1.75
    Controls maximum bond force in the quartic bond potential (Eq. 3). Chosen by hand as a control variable; not fitted to external data.
  • Critical fracture rate scaling exponent = -0.6
    Exponent fitted to the epsilon_crit vs N_mol data in Figure 6. This is an output characterization of the simulation results, not an input parameter.
  • Arrhenius prefactor and force coupling (E_f, f)
    Rate constants k_f in Figure 8 are fit to the exponential form k_f ~ exp(-(E_f - f epsilon_dot)/k_B T); the parameters are not tabulated.
assumptions (6)
  • domain assumption The Kremer-Grest bead-spring model with WCA and FENE potentials is a valid representation of unentangled polymer melts.
    Used throughout; equilibrium properties match previous simulations (ref 46).
  • domain assumption The quartic bond potential (Eq. 3) from ref 29 captures irreversible bond scission under flow.
    Replaces FENE after equilibration; parameters from Yin et al. 2020.
  • standard math The generalized Kraynik-Reinelt method maintains periodic boundary conditions without artifacts under steady-state UEF.
    Implements the method of Dobson (ref 44) as in the LAMMPS UEF package.
  • domain assumption Atomic SLLOD and atomic thermostat artifacts are negligible in dense melts.
    Acknowledged as a limitation by the authors; if not negligible, results could be biased.
  • domain assumption Simulations reach steady state within the stated durations (3x10^3 to 3x10^4 time units).
    Based on plateaus in viscosity time series (Figures S4-S5).
  • domain assumption Chains with N=20-80 are unentangled since N_e ~ 85.
    Uses entanglement length estimate from ref 37.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Mechanical Degradation of Unentangled Polymer Melts under Uniaxial Extensional Flow." pith.science (2026). https://pith.science/paper/7KVWSO5B

@misc{pith2026250603945,
  author       = {Pith},
  title        = {Pith review of: Mechanical Degradation of Unentangled Polymer Melts under Uniaxial Extensional Flow},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/7KVWSO5B}},
  note         = {Machine review of arXiv:2506.03945}
}
read the original abstract

Complex flow fields govern the deformation of polymers in various manufacturing processes. However, high flow rates may trigger reaction events (i.e., bond breaking or undesirable reaction of mechanophores) in raw polymeric materials, leading to the mechanical or functional debasement of manufactured structures. Additionally, it is difficult to fully characterize such molecular-level flow in the laboratory due to time- and length-scale limits. In this study, we perform non-equilibrium molecular dynamics (NEMD) simulations to explore the rheological and mechanical degradation of unentangled polymer melts under uniaxial extensional flow (UEF), allowing for chain breaking. Our simulations demonstrate shear thickening-thinning-thickening stages with the increase of UEF extension rates, resulting from flow-induced changes of chain conformation. With further increasing UEF extension rates, a bond-breaking potential leads to another flow thinning stage. Interestingly, fracture kinetics is originally first-order owing to the need for highly stretched polymer chains before bond fracture. It is no longer first-order when bond fracture is instigated before chains are stretched. Our computational work provides insight into the optimal design of the manufacturing process for polymeric materials.

Figures

Figures reproduced from arXiv: 2506.03945 by the authors.

Figure 1
Figure 1. (a) Schematic diagram of the unentangled polymer melt under uniaxial extensional flow. Extension is applied along z axis, and contraction is applied along the x and y axes. (b-c) The predicted average magnitudes of (b) end-to-end vector 𝐑!"! and (c) mean-square end-to-end vector 𝐑!"! # of equilibrium polymer melts with chain lengths of N = 20 – 80. The scattered symbols are results simulated using the LJ+QUARTIC pot… view at source ↗
Figure 2
Figure 2. The average steady-state values of UEF viscosity 𝜂$%& '"!()* for unentangled polymer melts with different chain lengths of (a) N = 20, (b) N = 40, (c) N = 60 and (d) N = 80, as a function of UEF extension rates, 𝜀̇. Solid circles and squares show the calculated 𝜂$%& '"!()* using LJ+FENE and LJ+QUARTIC (B2 = 1.25) potentials, respectively. Stages I – V represent the typical UEF-dependent rheological behavior: I – New… view at source ↗
Figure 3
Figure 3. The average steady-state values of the order parameter 𝑆+ '"!()*, chain stretch ratio 𝜆' '"!()* and bond force 𝐹,-.) '"!()* for unentangled polymer melts simulated by the LJ+FENE potential, as a function of UEF extension rates 𝜀̇. Different chain lengths: (a) N = 20, (b) N = 40, (c) N = 60 and (d) N = 80 were simulated. The ranges of 𝜀̇filled with different background colors represent the flow stages illustrated in … view at source ↗

Discussion (0). Sign in to comment.

Reference graph

Works this paper leans on

57 extracted references · 55 canonical work pages

  1. [1]

    A.; Daly, R

    Willis-Fox, N.; Rognin, E.; Aljohani, T. A.; Daly, R. Polymer mechanochemistry: Manufacturing is now a force to be reckoned with. Chem 2018, 4, 2499

  2. [2]

    A.; Basu, A.; Behrou, R.; Boechler, N.; Boydston, A

    Ghanem, M. A.; Basu, A.; Behrou, R.; Boechler, N.; Boydston, A. J.; Craig, S. L.; Lin, Y.; Lynde, B. E.; Nelson, A.; Shen, H.; et al. The role of polymer mechanochemistry in responsive materials and additive manufacturing. Nat. Rev. Mater. 2021, 6, 84

  3. [3]

    A.; Hamilton, A.; Yang, J.; Cremar, L

    Davis, D. A.; Hamilton, A.; Yang, J.; Cremar, L. D.; Van Gough, D.; Potisek, S. L.; Ong, M. T.; Braun, P. V.; Martínez, T. J.; White, S. R.; et al. Force-induced activation of covalent bonds in mechanoresponsive polymeric materials. Nature 2009, 459, 68

  4. [4]

    A.; Lamuta, C.; Kim, H.; Leal, C.; Sottos, N

    Kim, T. A.; Lamuta, C.; Kim, H.; Leal, C.; Sottos, N. R. Interfacial force-focusing effect in mechanophore-linked nanocomposites. Adv. Sci. 2020, 7, 1903464

  5. [5]

    Ramirez, A. L. B.; Kean, Z. S.; Orlicki, J. A.; Champhekar, M.; Elsakr, S. M.; Krause, W. E.; Craig, S. L. Mechanochemical strengthening of a synthetic polymer in response to typically destructive shear forces. Nat. Chem. 2013, 5, 757

  6. [6]

    Self-strengthening of cross-linked elastomers via the use of dynamic covalent macrocyclic mechanophores

    Kida, J.; Aoki, D.; Otsuka, H. Self-strengthening of cross-linked elastomers via the use of dynamic covalent macrocyclic mechanophores. ACS Macro Lett. 2021, 10, 558

  7. [7]

    A novel design of multi-mechanoresponsive and mechanically strong hydrogels

    Chen, H.; Yang, F.; Chen, Q.; Zheng, J. A novel design of multi-mechanoresponsive and mechanically strong hydrogels. Adv. Mater. 2017, 29, 1606900

  8. [8]

    P.; Creton, C

    Ducrot, E.; Chen, Y.; Bulters, M.; Sijbesma, R. P.; Creton, C. Toughening elastomers with sacrificial bonds and watching them break. Science 2014, 344, 186. 23

Show all 57 references
  1. [9]

    Matsuda, T.; Kawakami, R.; Namba, R.; Nakajima, T.; Gong, J. P. Mechanoresponsive self - growing hydrogels inspired by muscle training. Science 2019, 363, 504

  2. [10]

    Wang, J.; Piskun, I.; Craig, S. L. Mechanochemical strengthening of a multi -mechanophore benzocyclobutene polymer. ACS Macro Lett. 2015, 4, 834

  3. [11]

    Let spiropyran help polymers feel force! Prog

    Li, M.; Zhang, Q.; Zhou, Y.-N.; Zhu, S. Let spiropyran help polymers feel force! Prog. Polym. Sci. 2018, 79, 26

  4. [12]

    H.; Chang, C.-C.; Craig, S

    Lin, Y.; Barbee, M. H.; Chang, C.-C.; Craig, S. L. Regiochemical effects on mechanophore activation in bulk materials. J. Am. Chem. Soc. 2018, 140, 15969

  5. [13]

    A.; Moore, J

    May, P. A.; Moore, J. S. Polymer mechanochemistry: Techniques to generate molecular force via elongational flows. Chem. Soc. Rev. 2013, 42, 7497

  6. [14]

    A.; Göstl, R.; Daly, R

    Willis-Fox, N.; Rognin, E.; Baumann, C.; Aljohani, T. A.; Göstl, R.; Daly, R. Going with the flow: Tunable flow-induced polymer mechanochemistry. Adv. Funct. Mater. 2020, 30, 2002372

  7. [15]

    G.; Desai, P

    Larson, R. G.; Desai, P. S. Modeling the rheology of polymer melts and solutions. Annu. Rev. Fluid Mech. 2015, 47, 47

  8. [16]

    D.; Daivis, P

    Hajizadeh, E.; Todd, B. D.; Daivis, P. J. Nonequilibrium molecular dynamics simulation of dendrimers and hyperbranched polymer melts undergoing planar elongational flow. J. Rheol. 2014, 58, 281

  9. [17]

    Elongational viscosity of weakly entangled polymer melt via coarse-grained molecular dynamics simulation

    Murashima, T.; Hagita, K.; Kawakatsu, T. Elongational viscosity of weakly entangled polymer melt via coarse-grained molecular dynamics simulation. Nihon Reoroji Gakkaishi 2018, 46, 207

  10. [18]

    C.; Ge, T.; Rubinstein, M.; Grest, G

    O’Connor, T. C.; Ge, T.; Rubinstein, M.; Grest, G. S. Topological linking drives anomalous thickening of ring polymers in weak extensional flows. Phys. Rev. Lett. 2020, 124, 027801

  11. [19]

    Interchain pressure effect in extensional flows of entangled polymer melts

    Marrucci, G.; Ianniruberto, G. Interchain pressure effect in extensional flows of entangled polymer melts. Macromolecules 2004, 37, 3934

  12. [20]

    Structure and dynamics of dilute polymer solutions under shear flow via nonequilibrium molecular dynamics

    Aust, C.; Kröger, M.; Hess, S. Structure and dynamics of dilute polymer solutions under shear flow via nonequilibrium molecular dynamics. Macromolecules 1999, 32, 5660

  13. [21]

    J.; Evans, D

    Daivis, P. J.; Evans, D. J.; Morriss, G. P. Computer simulation study of the comparative rheology of branched and linear alkanes. J. Chem. Phys. 1992, 97, 616

  14. [22]

    Castrejón-González, O.; Castillo-Tejas, J.; Manero, O.; Alvarado, J. F. J. Structure factor and rheology of chain molecules from molecular dynamics. J. Chem. Phys. 2013, 138, 184901. 24

  15. [23]

    Structural mechanism for viscosity of semiflexible polymer melts in shear flow

    Xu, X.; Chen, J. Structural mechanism for viscosity of semiflexible polymer melts in shear flow. ACS Macro Lett. 2017, 6, 331

  16. [24]

    H.; Edwards, B

    Nafar Sefiddashti, M. H.; Edwards, B. J.; Khomami, B. Steady shearing flow of a moderately entangled polyethylene liquid. J. Rheol. 2016, 60, 1227

  17. [25]

    M.; Baig, C

    Cho, S.; Kim, J. M.; Baig, C. Scaling characteristics of rotational dynamics and rheology of linear polymer melts in shear flow. Macromolecules 2020, 53, 3030

  18. [26]

    Kremer, K.; Grest, G. S. Dynamics of entangled linear polymer melts: a molecular ‐ dynamics simulation. J. Chem. Phys. 1990, 92, 5057

  19. [27]

    Stevens, M. J. Interfacial fracture between highly cross-linked polymer networks and a solid surface: effect of interfacial bond density. Macromolecules 2001, 34, 2710

  20. [28]

    Stevens, M. J. Manipulating connectivity to control fracture in network polymer adhesives. Macromolecules 2001, 34, 1411

  21. [29]

    Topological origin of strain induced damage of multi-network elastomers by bond breaking

    Yin, Y.; Bertin, N.; Wang, Y.; Bao, Z.; Cai, W. Topological origin of strain induced damage of multi-network elastomers by bond breaking. Extreme Mech. Lett. 2020, 40, 100883

  22. [30]

    S.; Silberstein, M

    Manivannan, M. S.; Silberstein, M. N. Computational investigation of shear driven mechanophore activation at interfaces. Extreme Mech. Lett. 2016, 8, 6

  23. [31]

    Atomistic modeling framework for a cyclobutane-based mechanophore-embedded nanocomposite for damage precursor detection

    Koo, B.; Chattopadhyay, A.; Dai, L. Atomistic modeling framework for a cyclobutane-based mechanophore-embedded nanocomposite for damage precursor detection. Comput. Mater. Sci. 2016, 120, 135

  24. [32]

    Dimeric anthracene-based mechanophore for damage precursor detection in epoxy -based thermoset polymer matrix: Characterization and atomistic modeling

    Koo, B.; Nofen, E.; Chattopadhyay, A.; Dai, L. Dimeric anthracene-based mechanophore for damage precursor detection in epoxy -based thermoset polymer matrix: Characterization and atomistic modeling. Comput. Mater. Sci. 2017, 133, 167

  25. [33]

    G.; Khomami, B.; Sureshkumar, R

    Sim, H. G.; Khomami, B.; Sureshkumar, R. Flow –induced chain scission in dilute polymer solutions: Algorithm development and results for scission dynamics in elongational flow. J. Rheol. 2007, 51, 1223

  26. [34]

    A.; Daly, R

    Rognin, E.; Willis -Fox, N.; Aljohani, T. A.; Daly, R. A multiscale model for the rupture of linear polymers in strong flows. J. Fluid Mech. 2018, 848, 722

  27. [35]

    Fast parallel algorithms for short-range molecular dynamics

    Plimpton, S. Fast parallel algorithms for short-range molecular dynamics. J. Comput. Phys. 1995, 117, 1

  28. [36]

    Rubin, R. J. Random‐walk model of chain‐polymer adsorption at a surface. J. Chem. Phys. 2004, 43, 2392. 25

  29. [37]

    S.; Foteinopoulou, K.; Kröger, M

    Hoy, R. S.; Foteinopoulou, K.; Kröger, M. Topological analysis of polymeric melts: Chain - length effects and fast-converging estimators for entanglement length. Phys. Rev. E 2009, 80, 031803

  30. [38]

    D.; Chandler, D.; Andersen, H

    Weeks, J. D.; Chandler, D.; Andersen, H. C. Perturbation theory of the thermodynamic properties of simple liquids. J. Chem. Phys. 2003, 55, 5422

  31. [39]

    P.; Tildesley, D

    Allen, M. P.; Tildesley, D. J. Computer simulation of liquids; Oxford University Press (Oxford) 2017

  32. [40]

    A unified formulation of the constant temperature molecular dynamics methods

    Nosé, S. A unified formulation of the constant temperature molecular dynamics methods. J. Chem. Phys. 1984, 81, 511

  33. [41]

    Hoover, W. G. Canonical dynamics: Equilibrium phase-space distributions. Phys. Rev. A 1985, 31, 1695

  34. [42]

    J.; Searles, D

    Evans, D. J.; Searles, D. J.; Williams, S. R. Fundamentals of classical statistical thermodynamics: Dissipation, relaxation, and fluctuation theorems; John Wiley & Sons 2016

  35. [43]

    M.; Reinelt, D

    Kraynik, A. M.; Reinelt, D. A. Extensional motions of spatially periodic lattices. Int. J. Multiphase Flow 1992, 18, 1045

  36. [44]

    Periodic boundary conditions for long-time nonequilibrium molecular dynamics simulations of incompressible flows

    Dobson, M. Periodic boundary conditions for long-time nonequilibrium molecular dynamics simulations of incompressible flows. J. Chem. Phys. 2014, 141, 184103

  37. [45]

    A.; Rutledge, G

    Nicholson, D. A.; Rutledge, G. C. Molecular simulation of flow -enhanced nucleation in n- eicosane melts under steady shear and uniaxial extension. J. Chem. Phys. 2016, 145, 244903

  38. [46]

    K.; Keblinski, P

    Sen, S.; Kumar, S. K.; Keblinski, P. Viscoelastic properties of polymer melts from equilibrium molecular dynamics simulations. Macromolecules 2005, 38, 650

  39. [47]

    R.; Long, R.; Tang, T

    Lavoie, S. R.; Long, R.; Tang, T. Modeling the mechanics of polymer chains with deformable and active bonds. J. Phys. Chem. B 2020, 124, 253

  40. [48]

    R.; Silberstein, M

    Buche, M. R.; Silberstein, M. N. Statistical mechanical constitutive theory of polymer networks: The inextricable links between distribution, behavior, and ensemble. Phys. Rev. E 2020, 102, 012501

  41. [49]

    Treloar, L. R. G. The physics of rubber elasticity; Oxford University Press (Oxford) 1975

  42. [50]

    Nonlinear elongational rheology of unentangled polystyrene and poly(p -tert-butylstyrene) melts

    Matsumiya, Y.; Watanabe, H.; Masubuchi, Y.; Huang, Q.; Hassager, O. Nonlinear elongational rheology of unentangled polystyrene and poly(p -tert-butylstyrene) melts. Macromolecules 2018, 51, 9710. 26

  43. [51]

    Primitive chain network simulation of elongational flows of entangled linear chains: Stretch/orientation - induced reduction of monomeric friction

    Yaoita, T.; Isaki, T.; Masubuchi, Y.; Watanabe, H.; Ianniruberto, G.; Marrucci, G. Primitive chain network simulation of elongational flows of entangled linear chains: Stretch/orientation - induced reduction of monomeric friction. Macromolecules 2012, 45, 2773

  44. [52]

    Modeling unentangled polystyrene melts in fast elongational flows

    Ianniruberto, G.; Brasiello, A.; Marrucci, G. Modeling unentangled polystyrene melts in fast elongational flows. Macromolecules 2019, 52, 4610

  45. [53]

    Narimissa, E.; Wagner, M. H. Modeling nonlinear rheology of unentangled polymer melts based on a single integral constitutive equation. J. Rheol. 2019, 64, 129

  46. [54]

    J.; Matin, M

    Daivis, P. J.; Matin, M. L.; Todd, B. D. Nonlinear shear and elongational rheology of model polymer melts by non-equilibrium molecular dynamics. J. Non-Newton. Fluid Mech. 2003, 111, 1

  47. [55]

    Q.; Kausch, H.-H

    Nguyen, T. Q.; Kausch, H.-H. Chain scission in transient extensional flow kinetics and molecular weight dependence. J. Non-Newton. Fluid Mech. 1988, 30, 125

  48. [56]

    A.; Keller, A.; Rabin, Y

    Odell, J. A.; Keller, A.; Rabin, Y. Flow‐induced scission of isolated macromolecules. J. Chem. Phys. 1988, 88, 4022

  49. [57]

    # represents the transient state of the polymer melts, which is followed by a plateau in 𝜂!

    Reese, H. R.; Zimm, B. H. Fracture of polymer chains in extensional flow: Experiments with DNA, and a molecular-dynamics simulation. J. Chem. Phys. 1990, 92, 2650. S-1 Supporting Information for Mechanical Degradation of Unentangled Polymer Melts under Uniaxial Extensional Flo...

Pith tools

Reviewed August 7, 2026 · model on record in the stance chip above.