pith. sign in

arxiv: 2509.13146 · v3 · submitted 2025-09-16 · ❄️ cond-mat.soft · cond-mat.mtrl-sci· cond-mat.stat-mech

Multi-scale Modeling of the Electro-viscoelasticity of Charged Polymers in Combined Flow and Electric Fields

Zachary Wolfgram , Jeffrey G. Ethier , Matthew Grasinger This is my paper

Pith reviewed 2026-05-18 16:08 UTC · model grok-4.3

classification ❄️ cond-mat.soft cond-mat.mtrl-scicond-mat.stat-mech
keywords charged polymerselectro-viscoelasticityupper-convected Maxwell modelRouse modelmolecular dynamicsviscosity scalingelectric field dyadicKremer-Grest chains
0
0 comments X

The pith

The upper-convected time derivative of the electric field dyadic must be included in stress evolution to reproduce observed viscosity scaling for charged polymers in combined flow and electric fields.

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

The paper extends the Rouse bead-spring model by distributing charge density along polymer chains and derives the resulting viscoelastic stress under homogeneous shear and electric fields. Viscosity rises quadratically with select electric field components and depends on the relative orientation of flow and field, modulated by charge sequence relaxation time and dielectric constant. These results motivate the upper-convected electro-Maxwell (UCEM) continuum model, which adds upper-convected derivatives to the electric field dyadic so that polarization stresses properly track the stretching and rotation of charge pairs. Coarse-grained molecular dynamics simulations of Kremer-Grest chains with prescribed charge sequences confirm separate relaxation timescales for chain motion and charge redistribution, and show that standard models lacking the upper-convected terms miss the observed viscosity scaling.

Core claim

Standard continuum formulations of electro-viscoelasticity fail to capture the viscosity scaling observed in both the extended Rouse model and molecular dynamics simulations unless the upper-convected time derivative of the electric field dyadic is retained in the stress evolution equation; this term accounts for the convective stretching and rotation of charge pairs in flow, and the resulting UCEM model satisfies the second law while reproducing the quadratic field dependence and orientation effects.

What carries the argument

The upper-convected electro-Maxwell (UCEM) model, in which polarization stresses are expressed via an electric field dyadic that evolves under upper-convected time derivatives, analogous to the treatment of the conformation tensor in the upper-convected Maxwell fluid model.

If this is right

  • Viscosity increase scales quadratically with select electric field components and depends on field-flow orientation.
  • The UCEM model satisfies the second law of thermodynamics for the constitutive responses examined.
  • Distinct relaxation timescales separate overall chain dynamics from charge redistribution along the chain.
  • The model yields analytic constitutive relations for several canonical flows and field configurations.

Where Pith is reading between the lines

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

  • The same upper-convected treatment of the electric dyadic could be tested in extensional flows or inhomogeneous fields to check consistency with manufacturing processes.
  • The requirement for convective derivatives on the electric field suggests analogous corrections may be needed in other continuum models of charged macromolecules under flow.
  • Direct simulation of charge-pair trajectories in flow could provide an independent check on whether the upper-convected derivative is the minimal term needed.

Load-bearing premise

The derivation assumes homogeneous shear and electric fields together with a defined charge sequence on the polymer chains that produces distinct relaxation timescales for overall chain dynamics versus charge redistribution.

What would settle it

Viscosity measurements on charged polymer solutions in simple shear flow at fixed electric field strength, comparing data against predictions from the UCEM model versus a standard Maxwell model without upper-convected electric dyadic terms.

Figures

Figures reproduced from arXiv: 2509.13146 by Jeffrey G. Ethier, Matthew Grasinger, Zachary Wolfgram.

Figure 1
Figure 1. Figure 1: FIG. 1. A visual representation of the linear polymer chain made us [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Polarization stress as shown in equations 35-37 for a constant [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. A normalization of the pressure-driven velocity profile for [PITH_FULL_IMAGE:figures/full_fig_p006_4.png] view at source ↗
read the original abstract

The behavior of polymers in combined flow and electric fields underlies many manufacturing processes but remains poorly understood. To address this, we model charged polymers across scales. We extend the original Rouse model for a bead-spring chain to include a charge density distributed along the polymer chain, and derive the viscoelastic stress under homogeneous shear and electric fields. The viscosity increase depends on field-flow orientation and scales quadratically with select components of the electric field strength, modulated by the effective charge sequence relaxation time and dielectric constant. Inspired by this result, a new continuum model--the upper-convected electro-Maxwell (UCEM) model--is proposed, resembling an upper-convected Maxwell model with polarization stresses expressed through an electric field dyadic subject to upper-convected time derivatives. We analyze the constitutive response for several flows and electric field strengths, discussing limitations and demonstrating compliance with the second law of thermodynamics. Lastly, coarse-grained molecular dynamics (MD) simulations of Kremer-Grest chains with a defined charge sequence confirm the existence of distinct relaxation timescales for overall chain dynamics versus charge redistribution, consistent with the UCEM model predictions. Critically, we demonstrate that the upper-convected time derivative of the electric field dyadic is required in the evolution of stress to account for stretching and rotation of the charge pairs in flow, reproducing the viscosity scaling observed in both the Rouse and MD results; whereas standard continuum formulations without these terms fail to capture this observed scaling.

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

1 major / 3 minor

Summary. The manuscript develops a multi-scale framework for understanding the electro-viscoelastic response of charged polymers in combined flow and electric fields. It begins by extending the classical Rouse model to incorporate a charge density distributed along the bead-spring chain and derives an expression for the viscoelastic stress tensor under conditions of homogeneous shear and electric fields. This leads to a prediction that the viscosity increases with field-flow orientation and scales quadratically with select components of the electric field, modulated by an effective charge sequence relaxation time and the dielectric constant. Motivated by this microscopic result, the authors introduce the upper-convected electro-Maxwell (UCEM) continuum model, which modifies the standard upper-convected Maxwell model by expressing polarization stresses through an electric field dyadic that is evolved using upper-convected time derivatives. The constitutive behavior is examined for several flow types and field strengths, with a demonstration of thermodynamic consistency via the second law. Coarse-grained molecular dynamics simulations using Kremer-Grest chains with a prescribed charge sequence are employed to confirm the presence of distinct relaxation timescales for overall chain motion versus charge redistribution. The pivotal result is that inclusion of the upper-convected derivative is essential for capturing the advection, stretching, and rotation of charge pairs in flow, therebyre

Significance. If validated, this contribution is significant in establishing a continuum-level description of electro-viscoelasticity that is directly informed by and consistent with a microscopic polymer model. The analytical derivation from the extended Rouse model provides a clear mechanistic origin for the quadratic scaling, while the MD simulations offer independent support through the observation of separate relaxation timescales. A particular strength is the explicit contrast showing that the upper-convected terms are necessary to match the scaling, offering a concrete test of the model's validity. The UCEM model also satisfies thermodynamic requirements, enhancing its potential utility in engineering simulations of charged polymer systems under electric fields.

major comments (1)
  1. [Model assumptions and MD validation] The central claim depends on the assumption of homogeneous fields and a defined charge sequence producing distinct relaxation timescales for chain dynamics and charge redistribution (as noted in the model assumptions). The manuscript would benefit from a sensitivity analysis or explicit discussion of how deviations from these assumptions might affect the necessity of the upper-convected derivative in the UCEM model, since this underpins the comparison to standard formulations.
minor comments (3)
  1. [Abstract] The abstract mentions 'a defined charge sequence' but provides limited detail on its exact form; a reference to the specific sequence or the section where it is defined would improve clarity for readers.
  2. [Figures] Figure captions comparing Rouse, UCEM, and MD results should explicitly distinguish the data sources and note any parameters (such as the effective relaxation time) used in the comparisons.
  3. [Thermodynamic analysis] The thermodynamic compliance section would be strengthened by including a brief reference to the specific dissipation inequality or equation that demonstrates non-negative entropy production in the UCEM model.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the positive assessment of our manuscript and for the constructive comment. We address the major point below and will revise the manuscript to incorporate additional discussion of the model assumptions as suggested.

read point-by-point responses

  1. Referee: [Model assumptions and MD validation] The central claim depends on the assumption of homogeneous fields and a defined charge sequence producing distinct relaxation timescales for chain dynamics and charge redistribution (as noted in the model assumptions). The manuscript would benefit from a sensitivity analysis or explicit discussion of how deviations from these assumptions might affect the necessity of the upper-convected derivative in the UCEM model, since this underpins the comparison to standard formulations.

    Authors: We agree that further elaboration on the role of these assumptions would strengthen the presentation. The homogeneous-field assumption is required for the closed-form analytical derivation of the stress tensor from the extended Rouse model, which directly yields the quadratic scaling with electric-field components and isolates the necessity of the upper-convected derivative. Likewise, the prescribed charge sequence is essential for the separation of relaxation timescales observed in both the Rouse analysis and the Kremer-Grest MD simulations. In the revised manuscript we will add a new paragraph in the Conclusions section that explicitly discusses deviations from these assumptions. We will note that, under inhomogeneous fields, the UCEM constitutive relation remains locally applicable but must be solved numerically together with the electrostatic equations; the upper-convected terms continue to be required to capture advection and rotation of charge pairs, although the precise quadratic scaling may be modulated. A comprehensive sensitivity analysis via additional MD runs lies outside the scope of the present work due to computational cost; however, we will provide qualitative arguments, supported by the existing data, showing that the necessity of the upper-convected derivative persists whenever flow induces stretching or rotation of charge pairs. This revision directly responds to the referee's suggestion while preserving the core claims of the paper. revision: partial

Circularity Check

0 steps flagged

Derivation chain is self-contained with independent MD validation

full rationale

The paper first extends the Rouse bead-spring model by distributing charge density along the chain and derives the viscoelastic stress tensor under homogeneous shear and electric fields, yielding a quadratic viscosity scaling with electric field components modulated by the charge-sequence relaxation time. This analytic result then motivates the proposal of the UCEM continuum constitutive model, which incorporates upper-convected derivatives on the electric dyadic to reproduce the same scaling. Coarse-grained MD simulations with an identical defined charge sequence independently confirm the existence of distinct relaxation timescales for chain dynamics versus charge redistribution and the necessity of the upper-convected terms for matching the observed scaling; standard models without those terms fail. Because the Rouse derivation is first-principles, the MD validation is external to the continuum ansatz, and no load-bearing step reduces to a fitted parameter renamed as a prediction or to a self-citation chain, the overall argument remains non-circular.

Axiom & Free-Parameter Ledger

2 free parameters · 1 axioms · 1 invented entities

The work extends the Rouse model with distributed charge, introduces the UCEM constitutive relation as a new entity, and employs effective parameters for charge-sequence relaxation time and dielectric constant to recover the observed quadratic viscosity scaling.

free parameters (2)
  • effective charge sequence relaxation time
    Modulates the magnitude of viscosity increase with electric field strength in the derived scaling relation.
  • dielectric constant
    Scales the quadratic dependence of viscosity on selected electric field components.
axioms (1)
  • domain assumption The Rouse bead-spring model assumptions remain valid when a charge density is distributed along the chain under homogeneous shear and electric fields.
    This extension forms the starting point for the viscoelastic stress derivation.
invented entities (1)
  • upper-convected electro-Maxwell (UCEM) model no independent evidence
    purpose: Continuum constitutive model expressing polarization stresses through an electric field dyadic subject to upper-convected time derivatives.
    Proposed to bridge molecular-scale Rouse results to macroscopic flows and shown to match MD observations only when the upper-convected terms are retained.

pith-pipeline@v0.9.0 · 5808 in / 1628 out tokens · 64939 ms · 2026-05-18T16:08:20.789638+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

28 extracted references · 28 canonical work pages

  1. [1]

    merlin.mbs aapmrev4-1.bst 2010-07-25 4.21a (PWD, AO, DPC) hacked

    FUNCTION id.bst "merlin.mbs aapmrev4-1.bst 2010-07-25 4.21a (PWD, AO, DPC) hacked" ENTRY address archive archivePrefix author bookaddress booktitle chapter collaboration doi edition editor eid eprint howpublished institution isbn issn journal key language month note number organization pages primaryClass publisher school SLACcitation series title translat...

    work page 2010

  2. [2]

    merlin.mbs aipauth4-1.bst 2010-07-25 4.21a (PWD, AO, DPC) hacked

    FUNCTION id.bst "merlin.mbs aipauth4-1.bst 2010-07-25 4.21a (PWD, AO, DPC) hacked" ENTRY address archive archivePrefix author bookaddress booktitle chapter collaboration doi edition editor eid eprint howpublished institution isbn issn journal key language month note number organization pages primaryClass publisher school SLACcitation series title translat...

    work page 2010

  3. [3]

    Modeling of electro–viscoelastic dielectric elastomer: A continuum mechanics approach

    Subrat Kumar Behera, Deepak Kumar, and Somnath Sarangi. Modeling of electro–viscoelastic dielectric elastomer: A continuum mechanics approach. European Journal of Mechanics - A/Solids , 90:104369, 2021

    work page 2021

  4. [4]

    Journal of the Mechanics and Physics of Solids , 157:104625, 2021

    A complete thermo-electro-viscoelastic characterization of dielectric elastomers, part ii: Continuum modeling approach. Journal of the Mechanics and Physics of Solids , 157:104625, 2021

    work page 2021

  5. [5]

    Modelling of Electro - Viscoelastic Materials through Rate Equations

    Claudio Giorgi and Angelo Morro. Modelling of Electro - Viscoelastic Materials through Rate Equations . Materials , 16(10):3661, May 2023

    work page 2023

  6. [6]

    Larrondo and R

    L. Larrondo and R. St. John Manley. Electrostatic fiber spinning from polymer melts. I . Experimental observations on fiber formation and properties. Journal of Polymer Science: Polymer Physics Edition , 19(6):909--920, 1981

    work page 1981

  7. [7]

    W. M. Winslow. Induced Fibration of Suspensions . Journal of Applied Physics , 20(12):1137--1140, December 1949

    work page 1949

  8. [8]

    Electrorheological fluids with shear dependent viscosities: Steady flows

    Michael Růžička. Electrorheological fluids with shear dependent viscosities: Steady flows. In Electrorheological Fluids : Modeling and Mathematical Theory , volume 1748, pages 61--103. Springer Berlin Heidelberg, Berlin, Heidelberg, 2000

    work page 2000

  9. [9]

    Vít Průša and K. R. Rajagopal. Flow of an electrorheological fluid between eccentric rotating cylinders. Theoretical and Computational Fluid Dynamics , 26(1-4):1--21, January 2012

    work page 2012

  10. [10]

    Modeling, Mathematical and Numerical Analysis of Electrorheological Fluids

    Michael Růžička. Modeling, Mathematical and Numerical Analysis of Electrorheological Fluids . Applications of Mathematics , 49(6):565--609, December 2004

    work page 2004

  11. [11]

    Enhancing mechanical strength of carbon fiber-epoxy interface through electrowetting of fiber surface

    Xiaoming Chen, Kaiqiang Wen, Chunjiang Wang, Siyi Cheng, Shuo Wang, Hechuan Ma, Hongmiao Tian, Jie Zhang, Xiangming Li, and Jinyou Shao. Enhancing mechanical strength of carbon fiber-epoxy interface through electrowetting of fiber surface. Composites Part B: Engineering , 234:109751, April 2022

    work page 2022

  12. [12]

    Electrorheology of polymers melts with ions along the backbone

    Zak Wolfgram, Jeffrey Ethier, and Matthew Grasinger. Electrorheology of polymers melts with ions along the backbone. (in preparation)

    work page

  13. [13]

    Electric Field Enhances Shear Resistance of Polymer Melts via Orientational Polarization in Microstructures

    Miao Huo and Yunlong Guo. Electric Field Enhances Shear Resistance of Polymer Melts via Orientational Polarization in Microstructures . Polymers , 12(2):335, February 2020

    work page 2020

  14. [14]

    Basim Abu-Jdayil and Peter O. Brunn. Effects of electrode morphology on the slit flow of an electrorheological fluid. Journal of Non-Newtonian Fluid Mechanics , 63(1):45--61, March 1996

    work page 1996

  15. [15]

    Thompson, H

    Aidan P. Thompson, H. Metin Aktulga, Richard Berger, Dan S. Bolintineanu, W. Michael Brown, Paul S. Crozier, Pieter J. in 't Veld, Axel Kohlmeyer, Stan G. Moore, Trung Dac Nguyen, Ray Shan, Mark J. Stevens, Julien Tranchida, Christian Trott, and Steven J. Plimpton. LAMMPS - a flexible simulation tool for particle-based materials modeling at the atomic, me...

    work page 2022

  16. [16]

    Kurt Kremer and Gary S. Grest. Dynamics of entangled linear polymer melts: A molecular-dynamics simulation. The Journal of Chemical Physics , 92(8):5057--5086, April 1990

    work page 1990

  17. [17]

    Smook and Sissi de Beer

    Leon A. Smook and Sissi de Beer. Electrical Chain Rearrangement : What Happens When Polymers in Brushes Have a Charge Gradient ? Langmuir , 40(8):4142--4151, February 2024

    work page 2024

  18. [18]

    Aoyagi and M

    T. Aoyagi and M. Doi. Molecular dynamics simulation of entangled polymers in shear flow. Computational and Theoretical Polymer Science , 10(3-4):317--321, June 2000

    work page 2000

  19. [19]

    O’Connor, and Lars Pastewka

    Jan Mees, Thomas C. O’Connor, and Lars Pastewka. Entropic stress of grafted polymer chains in shear flow. The Journal of Chemical Physics , 159(9):094902, September 2023

    work page 2023

  20. [20]

    Likhtman

    Jing Cao and Alexei E. Likhtman. Simulating startup shear of entangled polymer melts. ACS Macro Letters , 4(12):1376--1381, 2015

    work page 2015

  21. [21]

    Introduction to polymer physics

    Masao Doi. Introduction to polymer physics . Oxford university press, 1996

    work page 1996

  22. [22]

    C.W. Macosko. Rheology: Principles, Measurements, and Applications . Wiley, 1994

    work page 1994

  23. [23]

    A theory of the linear viscoelastic properties of dilute solutions of coiling polymers

    Prince E Rouse, Jr. A theory of the linear viscoelastic properties of dilute solutions of coiling polymers. The Journal of Chemical Physics , 21(7):1272, 1953

    work page 1953

  24. [24]

    Elements of nonequilibrium statistical mechanics , volume 3

    Venkataraman Balakrishnan. Elements of nonequilibrium statistical mechanics , volume 3. Springer, 2008

    work page 2008

  25. [25]

    Stochastic processes in polymeric fluids: tools and examples for developing simulation algorithms

    Hans C \"O ttinger. Stochastic processes in polymeric fluids: tools and examples for developing simulation algorithms . Springer Science & Business Media, 2012

    work page 2012

  26. [26]

    merlin.mbs aipnum4-1.bst 2010-07-25 4.21a (PWD, AO, DPC) hacked

    FUNCTION id.bst "merlin.mbs aipnum4-1.bst 2010-07-25 4.21a (PWD, AO, DPC) hacked" ENTRY address archive archivePrefix author bookaddress booktitle chapter collaboration doi edition editor eid eprint howpublished institution isbn issn journal key language month note number organization pages primaryClass publisher school SLACcitation series title translati...

    work page 2010

  27. [27]

    merlin.mbs apsrev4-1.bst 2010-07-25 4.21a (PWD, AO, DPC) hacked

    FUNCTION id.bst "merlin.mbs apsrev4-1.bst 2010-07-25 4.21a (PWD, AO, DPC) hacked" ENTRY address archive archivePrefix author bookaddress booktitle chapter collaboration doi edition editor eid eprint howpublished institution isbn issn journal key language month note number organization pages primaryClass publisher school SLACcitation series title translati...

    work page 2010

  28. [28]

    merlin.mbs apsrmp4-1.bst 2010-07-25 4.21a (PWD, AO, DPC) hacked

    FUNCTION id.bst "merlin.mbs apsrmp4-1.bst 2010-07-25 4.21a (PWD, AO, DPC) hacked" ENTRY address archive archivePrefix author bookaddress booktitle chapter collaboration doi edition editor eid eprint howpublished institution isbn issn journal key language month note number organization pages primaryClass publisher school SLACcitation series title translati...

    work page 2010