pith. sign in

arxiv: 2604.05338 · v1 · submitted 2026-04-07 · ❄️ cond-mat.soft · cond-mat.stat-mech

Linear Viscoelasticity of Semidilute Unentangled Flexible Polymer Solutions

Amit Varakhedkar , P. Sunthar , J. Ravi Prakash This is my paper

Pith reviewed 2026-05-10 19:21 UTC · model grok-4.3

classification ❄️ cond-mat.soft cond-mat.stat-mech
keywords linear viscoelasticitysemidilute unentangled polymersBrownian dynamicsRouse dynamicsZimm dynamicssuccessive fine-grainingstorage modulusloss modulus
0
0 comments X

The pith

Brownian dynamics simulations show a crossover to Rouse-like dynamics in semidilute unentangled polymer solutions from screened interactions.

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

The paper investigates the linear viscoelastic response of flexible polymer solutions in dilute and semidilute unentangled regimes using Brownian dynamics simulations of bead-spring chains. In the dilute limit the simulations recover Zimm-like behavior with solvent-quality-dependent power-law scaling. In the semidilute unentangled regime a systematic crossover to Rouse-like dynamics occurs with rising concentration because excluded-volume and hydrodynamic interactions become screened. The computed storage modulus agrees well with experiments in both regimes while the loss modulus matches at low and intermediate frequencies; high-frequency deviations are removed by successive fine-graining extrapolation to yield quantitative predictions in the infinite-chain-length limit.

Core claim

In the semidilute unentangled regime the simulations reveal a systematic crossover to Rouse-like dynamics with increasing concentration due to the screening of excluded volume and hydrodynamic interactions. This is confirmed by excellent agreement with experimental storage modulus across both concentration regimes and loss modulus at low and intermediate frequencies, with high-frequency deviations arising from finite-chain discretization effects that are removed by the successive fine-graining technique to enable quantitative prediction of the loss modulus in the infinite-chain-length limit.

What carries the argument

The Brownian dynamics bead-spring model with hydrodynamic interactions and excluded-volume potentials, combined with the successive fine-graining extrapolation to remove finite-chain-length artifacts.

Load-bearing premise

The chosen bead-spring model and hydrodynamic screening rules faithfully reproduce the real dynamics of flexible polymer chains, and successive fine-graining correctly removes discretization artifacts without introducing new biases.

What would settle it

High-frequency loss-modulus data measured on extremely long chains in the semidilute unentangled regime that deviate from the infinite-chain extrapolation obtained by successive fine-graining.

Figures

Figures reproduced from arXiv: 2604.05338 by Amit Varakhedkar, J. Ravi Prakash, P. Sunthar.

Figure 1
Figure 1. Figure 1: FIG. 1. Non-dimensional relaxation modulus [PITH_FULL_IMAGE:figures/full_fig_p005_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Non-dimensional relaxation modulus [PITH_FULL_IMAGE:figures/full_fig_p006_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Non-dimensional dynamic moduli, [PITH_FULL_IMAGE:figures/full_fig_p006_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Comparison of dynamic moduli with the experimental data for polystyrene in an infinitely dilute solution in (a) [PITH_FULL_IMAGE:figures/full_fig_p007_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. Finite-chain effects in the loss modulus [PITH_FULL_IMAGE:figures/full_fig_p007_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. Successive fine-graining extrapolation of finite chain data for the loss modulus [PITH_FULL_IMAGE:figures/full_fig_p008_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7. Loss modulus obtained by applying the successive fine [PITH_FULL_IMAGE:figures/full_fig_p009_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8. Successive fine-graining extrapolation of finite size data for the loss modulus [PITH_FULL_IMAGE:figures/full_fig_p010_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: FIG. 9. Loss modulus [PITH_FULL_IMAGE:figures/full_fig_p011_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: FIG. 10. Comparison of dynamic moduli with the experimental data for poly(AN-co-IA) at concentrations: (a) [PITH_FULL_IMAGE:figures/full_fig_p012_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: FIG. 11. Successive fine-graining extrapolation of finite chain data for the loss modulus [PITH_FULL_IMAGE:figures/full_fig_p013_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: FIG. 12. Comparison of the loss modulus [PITH_FULL_IMAGE:figures/full_fig_p014_12.png] view at source ↗
read the original abstract

The linear viscoelastic response of flexible polymer solutions in the dilute and semidilute unentangled regimes is investigated using Brownian dynamics simulations. The relaxation modulus and dynamic moduli are computed over a wide range of concentrations and chain discretizations for both $\theta$ and good solvents to establish the connection between microscopic chain dynamics and macroscopic viscoelastic response. In the dilute limit, the simulations recover the expected Zimm-like behavior with solvent-quality-dependent power-law scaling in the intermediate time and frequency regimes, while in the semidilute unentangled regime a systematic crossover to Rouse-like dynamics is observed with increasing concentration due to the screening of excluded volume and hydrodynamic interactions. Comparison with experimental measurements shows excellent agreement for the storage modulus across both concentration regimes and for the loss modulus at low and intermediate frequencies, with deviations at high frequencies as a result of finite-chain discretization effects. These finite-chain length effects are systematically accounted for using the successive fine-graining technique, enabling quantitative prediction of the loss modulus in the infinite-chain length limit.

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 / 1 minor

Summary. The manuscript uses Brownian dynamics simulations of bead-spring chains to compute the relaxation modulus and dynamic moduli (G' and G'') of flexible polymer solutions in the dilute and semidilute unentangled regimes. It reports recovery of Zimm-like scaling in dilute solutions (solvent-quality dependent), a systematic crossover to Rouse-like dynamics in the semidilute regime due to screening of excluded-volume and hydrodynamic interactions, and quantitative agreement with experimental data for G' across regimes and for G'' at low-to-intermediate frequencies. High-frequency deviations in G'' are attributed to finite-chain discretization and are corrected via successive fine-graining to obtain the infinite-chain limit.

Significance. If the central claims hold, the work supplies a direct microscopic-to-macroscopic connection for linear viscoelasticity in unentangled polymer solutions, with the fine-graining correction offering a route to infinite-chain predictions that can be compared to experiment. The explicit benchmarking against independent experimental moduli and the coverage of both θ and good solvents are positive features.

major comments (1)
  1. [Section describing successive fine-graining / high-frequency extrapolation] The headline quantitative result—the extrapolated infinite-chain G''(ω) at high frequencies—rests on the successive fine-graining procedure. The manuscript must demonstrate that the assumed extrapolation form isolates only discretization artifacts and converges to a result independent of the specific bead-spring potential, hydrodynamic tensor, and chosen N values; without such a test the corrected high-ω loss modulus may retain systematic bias (see abstract and the section describing the fine-graining technique).
minor comments (1)
  1. The abstract states that finite-chain effects are 'systematically accounted for' but does not quote the explicit functional form or fitting procedure used for the extrapolation; this should be stated clearly in the methods or results to permit independent verification.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their careful reading and constructive feedback on our manuscript. We address the single major comment below and have revised the manuscript to incorporate additional validation of the fine-graining procedure.

read point-by-point responses

  1. Referee: [Section describing successive fine-graining / high-frequency extrapolation] The headline quantitative result—the extrapolated infinite-chain G''(ω) at high frequencies—rests on the successive fine-graining procedure. The manuscript must demonstrate that the assumed extrapolation form isolates only discretization artifacts and converges to a result independent of the specific bead-spring potential, hydrodynamic tensor, and chosen N values; without such a test the corrected high-ω loss modulus may retain systematic bias (see abstract and the section describing the fine-graining technique).

    Authors: We agree that explicit validation of the extrapolation's robustness is required. The successive fine-graining procedure simulates progressively finer discretizations (larger N) and extrapolates the high-frequency G''(ω) using a form derived from the known high-frequency asymptotics of bead-spring chains. In the revised manuscript we have added a new subsection presenting results obtained with both Hookean and FENE spring potentials, with Oseen and Rotne-Prager hydrodynamic tensors, and across multiple overlapping ranges of N. In all cases the extrapolated high-frequency loss modulus converges to the same limiting curve, confirming that the correction isolates discretization effects without residual dependence on these modeling choices. The abstract has also been updated to note this validation. revision: yes

Circularity Check

0 steps flagged

No significant circularity; derivation chain is self-contained via simulation and external benchmarking.

full rationale

The paper computes relaxation and dynamic moduli from Brownian dynamics bead-spring simulations across concentrations and discretizations, observes the expected Zimm-to-Rouse crossover due to screening, and directly compares storage and loss moduli to independent experimental measurements. Finite-N effects are corrected via successive fine-graining extrapolation to the infinite-chain limit, but this is framed as a post-processing technique rather than a fit to the target observables. No load-bearing steps reduce by construction to fitted inputs, self-citations, or ansatzes imported from the authors' prior work; the central claims rest on the model's physics and external data agreement rather than internal redefinition.

Axiom & Free-Parameter Ledger

2 free parameters · 2 axioms · 0 invented entities

The central claim rests on the standard assumptions of Brownian dynamics for bead-spring chains and on the validity of the successive fine-graining extrapolation procedure; no new physical entities are introduced.

free parameters (2)
  • Number of beads per chain
    Discretization parameter that is varied and then extrapolated to the infinite-chain limit.
  • Concentration range
    Set of simulated concentrations used to observe the crossover.
axioms (2)
  • domain assumption Brownian dynamics with bead-spring chains and appropriate hydrodynamic and excluded-volume interactions accurately models real flexible polymer solutions.
    Core modeling assumption invoked throughout the simulation study.
  • ad hoc to paper The successive fine-graining technique removes finite-chain discretization artifacts without introducing systematic bias in the high-frequency regime.
    Method used to enable the infinite-chain prediction.

pith-pipeline@v0.9.0 · 5481 in / 1561 out tokens · 73864 ms · 2026-05-10T19:21:32.444332+00:00 · methodology

discussion (0)

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

Reference graph

Works this paper leans on

54 extracted references · 54 canonical work pages

  1. [1]

    Stress Tensor The dynamic properties such as relaxation modulus, zero- shear rate viscosity and dynamic moduli investigated in this work can be defined in terms of the components of the stress- tensor [38] for the polymer solution. The non-dimensional contribution to the stress tensor is given by the Kramers- Kirkwood expression, τp = 1 Nc * Nc ∑ ξ=1 Nb ∑...

    work page

  2. [2]

    Relaxation Modulus The relaxation modulusG(t)is obtained from equilibrium simulations using the Green-Kubo relation that relatesG(t) to the autocorrelation function of the stress tensor [39]. At equilibrium, the stress tensor is isotropic, hence the relaxation modulus is given by the expression, G(t) = 1 3 (Gxy(t) +G xz(t) +G yz(t))(14) where the componen...

    work page

  3. [3]

    Hereη p,0 is calculated from equilibrium simulations by inte- grating the relaxation modulusG(t)[40, 41], ηp,0 = Z ∞ 0 G(t)dt(17)

    Zero Shear Rate Viscosity While the study of shear viscosity at moderately high shear rates ˙γis important in non-linear rheology, the study of linear viscoelasticity primarily focuses on the polymeric component of the zero-shear rate viscosity, defined as:η p,0 =lim ˙γ→0 ηp. Hereη p,0 is calculated from equilibrium simulations by inte- grating the relaxa...

    work page

  4. [4]

    These proper- ties are typically obtained from oscillatory shear flow exper- iments

    Dynamic Moduli The elastic and viscous response of a viscoelastic fluid is generally characterised by the storageG′ and the lossG ′′ mod- uli, together referred to as the dynamic moduli. These proper- ties are typically obtained from oscillatory shear flow exper- iments. In the current simulations, in the limit of very small strain amplitude,G ′ andG ′′ a...

    work page

  5. [5]

    M. Doi, S. F. Edwards, and S. F. Edwards,The Theory of Poly- mer Dynamics, V ol. 73 (Oxford University Press, 1988)

    work page 1988

  6. [6]

    Rubinstein and R

    M. Rubinstein and R. H. Colby,Polymer Physics(Oxford uni- versity press, 2003)

    work page 2003

  7. [7]

    Rouse, Prince E., J

    J. Rouse, Prince E., J. Chem. Phys.21, 1272 (1953)

    work page 1953

  8. [8]

    B. H. Zimm, J. Chem. Phys.24, 269 (1956)

    work page 1956

  9. [9]

    Tschoegl, J

    N. Tschoegl, J. Chem. Phys.40, 473 (1964)

    work page 1964

  10. [10]

    Sahouani and T

    H. Sahouani and T. P. Lodge, Macromolecules25, 5632 (1992)

    work page 1992

  11. [11]

    De Gennes,Scaling concepts in polymer physics(Cornell university press, 1979)

    P.-G. De Gennes,Scaling concepts in polymer physics(Cornell university press, 1979)

    work page 1979

  12. [12]

    J. R. Prakash, Curr. Opin. Colloid Interface Sci.43, 63 (2019)

    work page 2019

  13. [13]

    Ahlrichs, R

    P. Ahlrichs, R. Everaers, and B. Dünweg, Phys. Rev. E64, 040501 (2001)

    work page 2001

  14. [14]

    X. Zhu, X. Chen, H. Saba, Y . Zhang, and H. Wang, Eur. Polym. J.48, 597 (2012)

    work page 2012

  15. [15]

    W. Paul, K. Binder, D. W. Heermann, and K. Kremer, J. Phys. II1, 37 (1991)

    work page 1991

  16. [16]

    Adam and M

    M. Adam and M. Delsanti, J. Phys. France45, 1513 (1984)

    work page 1984

  17. [17]

    Inoue, Y

    T. Inoue, Y . Yamashita, and K. Osaki, Macromolecules35, 9169 (2002)

    work page 2002

  18. [18]

    S. Pan, D. At Nguyen, T. Sridhar, P. Sunthar, and J. Ravi Prakash, J. Rheol.58, 339 (2014)

    work page 2014

  19. [19]

    Heo and R

    Y . Heo and R. G. Larson, J. Rheol.49, 1117 (2005)

    work page 2005

  20. [20]

    Clasen, J

    C. Clasen, J. P. Plog, W.-M. Kulicke, M. Owens, C. Macosko, L. E. Scriven, M. Verani, and G. H. McKinley, J. Rheol.50, 849 (2006)

    work page 2006

  21. [21]

    X. Chen, S. Liang, S.-W. Wang, and R. H. Colby, J. Rheol.62, 81 (2018)

    work page 2018

  22. [22]

    R. M. Johnson, J. L. Schrag, and J. D. Ferry, Polym. J.1, 742 (1970)

    work page 1970

  23. [23]

    Huang, R

    C.-C. Huang, R. G. Winkler, G. Sutmann, and G. Gompper, Macromolecules43, 10107 (2010)

    work page 2010

  24. [24]

    Varakhedkar, P

    A. Varakhedkar, P. Sunthar, and J. R. Prakash, Soft Matter22, 369 (2026)

    work page 2026

  25. [25]

    T. T. Pham, P. Sunthar, and J. R. Prakash, J. Non-Newtonian Fluid Mech.149, 9 (2008)

    work page 2008

  26. [26]

    Sunthar and J

    P. Sunthar and J. R. Prakash, Macromolecules38, 617 (2005)

    work page 2005

  27. [27]

    Prabhakar, J

    R. Prabhakar, J. R. Prakash, and T. Sridhar, J. Rheol.48, 1251 (2004)

    work page 2004

  28. [28]

    Sasmal, K.-W

    C. Sasmal, K.-W. Hsiao, C. M. Schroeder, and J. Ravi Prakash, J. Rheol.61, 169 (2016)

    work page 2016

  29. [29]

    A. Jain, B. Dünweg, and J. R. Prakash, Phys. Rev. Lett.109, 088302 (2012)

    work page 2012

  30. [30]

    Saadat and B

    A. Saadat and B. Khomami, J. Rheol.59, 1507 (2015)

    work page 2015

  31. [31]

    H. C. Öttinger, J. Chem. Phys.86, 3731 (1987)

    work page 1987

  32. [32]

    H. C. Öttinger, J. Chem. Phys.90, 463 (1989)

    work page 1989

  33. [33]

    J. R. Prakash and H. C. Öttinger, J. Non-Newtonian Fluid Mech. 71, 245 (1997)

    work page 1997

  34. [34]

    J. R. Prakash, J. Rheol.46, 1353 (2002)

    work page 2002

  35. [35]

    Kröger, A

    M. Kröger, A. Alba-Pérez, M. Laso, and H. C. Öttinger, J. Chem. Phys.113, 4767 (2000)

    work page 2000

  36. [36]

    H. C. Öttinger,Stochastic processes in polymeric fluids: tools and examples for developing simulation algorithms(Springer Science & Business Media, 2012)

    work page 2012

  37. [37]

    Soddemann, B

    T. Soddemann, B. Dünweg, and K. Kremer, Eur. Phys. J. E6, 409 (2001)

    work page 2001

  38. [38]

    Santra, K

    A. Santra, K. Kumari, R. Padinhateeri, B. Dünweg, and J. R. Prakash, Soft Matter15, 7876 (2019)

    work page 2019

  39. [39]

    J. A. Anderson, J. Glaser, and S. C. Glotzer, Comput. Mater. Sci173, 109363 (2020)

    work page 2020

  40. [40]

    A. M. Fiore, F. Balboa Usabiaga, A. Donev, and J. W. Swan, J. Chem. Phys.146, 124116 (2017)

    work page 2017

  41. [41]

    D. Robe, A. Santra, G. H. McKinley, and J. R. Prakash, Macro- molecules57, 4220 (2024)

    work page 2024

  42. [42]

    R. B. Bird, C. F. Curtiss, R. C. Armstrong, and O. Hassager, Dynamics of Polymeric Liquids, Volume 2: Kinetic Theory(Wi- ley, 1987)

    work page 1987

  43. [43]

    Wittmer, H

    J. Wittmer, H. Xu, O. Benzerara, and J. Baschnagel, Mol. Phys. 113, 2881 (2015)

    work page 2015

  44. [44]

    Fixman, Macromolecules14, 1710 (1981)

    M. Fixman, Macromolecules14, 1710 (1981)

    work page 1981

  45. [45]

    W. B. Lee and K. Kremer, Macromolecules42, 6270 (2009)

    work page 2009

  46. [46]

    H. C. Öttinger and Y . Rabin, J. Non-Newtonian Fluid Mech.33, 53 (1989)

    work page 1989

  47. [47]

    R. G. Larson, inConstitutive Equations for Polymer Melts and Solutions, Butterworths Series in Chemical Engineering, edited by R. G. Larson (Butterworth-Heinemann, 1988) pp. 219–270

    work page 1988

  48. [48]

    C. S. Landry, Ph.D. thesis, University of Wisconsin (1985)

    work page 1985

  49. [49]

    Sunthar and J

    P. Sunthar and J. R. Prakash, Europhys. Lett.75, 77 (2006)

    work page 2006

  50. [50]

    J. R. Prakash and H. C. Öttinger, Macromolecules32, 2028 (1999)

    work page 2028

  51. [51]

    Garcia De la Torre, M

    J. Garcia De la Torre, M. C. Lopez Martinez, and M. M. Tirado, Macromolecules17, 2715 (1984)

    work page 1984

  52. [52]

    J. J. Freire, A. Rey, and J. Garcia de la Torre, Macromolecules 19, 457 (1986)

    work page 1986

  53. [53]

    Osaki, Macromolecules5, 141 (1972)

    K. Osaki, Macromolecules5, 141 (1972)

    work page 1972

  54. [54]

    W. H. Press, S. A. Teukolsky, W. T. Vetterling, and B. P. Flan- nery, The art of scientific computing (1992)

    work page 1992