pith. sign in

arxiv: 2605.01261 · v1 · submitted 2026-05-02 · ❄️ cond-mat.soft · cond-mat.stat-mech· hep-th

Loop expansion in polymer field theory: application to phase separation

Kiyoharu Kawana , Kyosuke Adachi This is my paper

Pith reviewed 2026-05-09 18:37 UTC · model grok-4.3

classification ❄️ cond-mat.soft cond-mat.stat-mechhep-th
keywords polymer phase separationloop expansionrandom phase approximationfield theorybinodalcoexistence densitymolecular dynamics
0
0 comments X

The pith

Identifying inverse polymer density as the expansion parameter refines RPA predictions for polymer phase separation binodals.

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

The paper develops a loop expansion for homopolymer field theories by treating the inverse polymer density as the small parameter, analogous to Planck's constant in quantum field theory. This generates systematic corrections to the random phase approximation free energy, starting with the leading RPA+ term and the next RPA++ term. When the RPA+ binodal is compared to molecular dynamics simulations of bead-spring chains, it improves the predicted density of the dilute coexisting phase while the critical point error stays similar to plain RPA. Accurate binodals matter for modeling liquid-liquid phase separation in both synthetic polymers and biological systems such as protein condensates, where the basic RPA loses quantitative accuracy away from high densities.

Core claim

By identifying the inverse polymer density ρ^{-1} as ħ, we derive the leading and next-to-leading corrections to the RPA free energy, denoted RPA+ and RPA++. Testing the RPA+ binodal against molecular dynamics simulations of bead-spring chains with Gaussian pair interactions shows that it qualitatively improves the dilute-phase coexistence density over the RPA, while the critical point error remains comparable.

What carries the argument

The loop expansion in homopolymer field theory with inverse density ρ^{-1} identified as ħ, which systematically generates the RPA+ and RPA++ corrections to the random phase approximation free energy.

If this is right

  • The RPA+ correction improves the description of the dilute branch of the binodal curve compared with plain RPA.
  • The location of the critical point remains comparably accurate to the RPA level at this order.
  • The method supplies a perturbative series that can be extended order by order for further refinement of phase coexistence predictions.
  • The expansion applies to homopolymer systems with the Gaussian pair interactions used in the simulations.

Where Pith is reading between the lines

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

  • Computing the RPA++ term explicitly would test whether the series converges and whether higher orders can reduce the remaining critical-point discrepancy.
  • The same identification of inverse density as ħ could be applied to field theories that include specific interactions or polydispersity.
  • Direct comparison of the predicted free-energy corrections against simulation data at varying densities would clarify the radius of convergence of the expansion.

Load-bearing premise

The inverse polymer density provides a valid and controlled small parameter for a perturbative expansion around the RPA that captures the dominant corrections in the dilute regime without requiring resummation.

What would settle it

Comparing the RPA++ binodal prediction directly to the same molecular dynamics data to check whether the next order further reduces the dilute-phase density error beyond the improvement already seen with RPA+.

Figures

Figures reproduced from arXiv: 2605.01261 by Kiyoharu Kawana, Kyosuke Adachi.

Figure 1
Figure 1. Figure 1: FIG. 1. A graphical representation of the connected correlation view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. The LO diagrams contributing to the free energy. view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. Binodal lines for phase separation in a polymer solution view at source ↗
read the original abstract

Liquid-liquid phase separation underlies phenomena ranging from protein condensate formation to the phase coexistence of synthetic polymers. Although the random phase approximation (RPA) is widely used to predict such phase behavior, its quantitative accuracy for binodals of polymer solutions, particularly outside the high-density regime, remains incompletely characterized. Here, we develop a field theoretic loop expansion in homopolymer systems by identifying the inverse polymer density $\rho^{-1}$ as the Planck constant $\hbar$ in quantum field theory. We calculate the leading-order and next-to-leading-order corrections to the RPA free energy, denoted as RPA+ and RPA++, respectively. Testing the binodal predicted by the RPA+ against molecular dynamics simulations of bead-spring chains with Gaussian pair interactions, we find that the RPA+ qualitatively improves the dilute-phase coexistence density over the RPA, while the critical point error remains comparable to that of the RPA. Our results establish the loop expansion as a systematic route for refining the RPA-based binodal predictions for polymer phase separation.

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

2 major / 2 minor

Summary. The paper develops a loop expansion in homopolymer field theory by identifying the inverse polymer density ρ^{-1} as the Planck constant ħ. It computes the first correction (RPA+) and next-to-leading correction (RPA++) to the RPA free energy, then compares the RPA+ binodal predictions to molecular dynamics simulations of bead-spring chains with Gaussian interactions. The abstract reports that RPA+ qualitatively improves the dilute-phase coexistence density relative to RPA while the critical-point error remains comparable, and concludes that the loop expansion provides a systematic route to refine RPA-based binodal predictions.

Significance. If the expansion parameter is demonstrably small and the series controlled, the work would supply a field-theoretic framework for systematically improving mean-field predictions of polymer phase separation, with relevance to both synthetic polymers and protein condensates. Direct comparison to independent MD simulations is a strength that grounds the claims in concrete data. However, the current presentation leaves open whether the reported improvements arise from a controlled perturbative correction or from an ad-hoc adjustment, which limits the immediate significance.

major comments (2)
  1. [Abstract] Abstract (and the definition of the loop expansion): identifying ħ ≡ ρ^{-1} makes the expansion parameter O(1) or larger precisely in the dilute regime where RPA+ is claimed to improve the binodal. The manuscript does not show that the RPA++ shift is parametrically smaller than the RPA+ correction at the reported coexistence densities, so the observed change cannot yet be attributed to a controlled perturbative series rather than an uncontrolled correction.
  2. [Abstract] Abstract (comparison to MD): the claim of qualitative improvement in dilute-phase density is stated without error bars on the simulation data, without quantitative metrics (e.g., relative error or χ²), and without explicit values of the expansion parameter ρ^{-1} at the binodal points. This leaves the support for the central claim of systematic improvement incomplete.
minor comments (2)
  1. Define the precise meaning of RPA+ and RPA++ (free-energy corrections or effective potentials) at first use and ensure consistent notation throughout.
  2. Add a short paragraph or table listing the numerical values of ρ^{-1} at the critical point and at representative dilute binodal points to allow readers to assess the size of the expansion parameter directly.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and constructive feedback on our manuscript. We address each major comment below and indicate the revisions planned for the next version.

read point-by-point responses

  1. Referee: [Abstract] Abstract (and the definition of the loop expansion): identifying ħ ≡ ρ^{-1} makes the expansion parameter O(1) or larger precisely in the dilute regime where RPA+ is claimed to improve the binodal. The manuscript does not show that the RPA++ shift is parametrically smaller than the RPA+ correction at the reported coexistence densities, so the observed change cannot yet be attributed to a controlled perturbative series rather than an uncontrolled correction.

    Authors: We agree that ρ^{-1} is O(1) or larger in the dilute regime, which raises a legitimate question about the control of the perturbative series there. The loop expansion remains formally systematic, but we acknowledge that the improvement cannot be attributed to a controlled expansion without evidence of convergence. In the revised manuscript we will add a direct comparison (table or plot) of the RPA+ and RPA++ corrections evaluated at the coexistence densities, together with the numerical values of ρ^{-1} at those points, so that readers can assess the relative size of successive terms. revision: partial

  2. Referee: [Abstract] Abstract (comparison to MD): the claim of qualitative improvement in dilute-phase density is stated without error bars on the simulation data, without quantitative metrics (e.g., relative error or χ²), and without explicit values of the expansion parameter ρ^{-1} at the binodal points. This leaves the support for the central claim of systematic improvement incomplete.

    Authors: We accept this criticism. The revised manuscript will include error bars on the MD coexistence densities (computed from multiple independent runs), quantitative measures of improvement such as relative errors in the dilute-phase density for both RPA and RPA+, and the explicit values of ρ^{-1} at the reported binodal points. These additions will make the comparison more rigorous and allow a clearer evaluation of the claimed improvement. revision: yes

Circularity Check

0 steps flagged

No significant circularity; derivation is self-contained with external validation.

full rationale

The paper introduces the identification ρ^{-1} ≡ ħ to enable a standard loop expansion in the homopolymer field theory, then applies established diagrammatic rules to compute the RPA+ and RPA++ corrections to the free energy. Binodal predictions are tested directly against independent molecular dynamics simulations of bead-spring chains, providing an external benchmark. No quoted step reduces a claimed prediction to a fitted input by construction, a self-citation chain, or a definitional tautology; the central claim of systematic refinement rests on the computed corrections and their comparison to simulation data rather than internal equivalence.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on the validity of the field-theoretic identification and perturbative expansion in the polymer context, plus the representativeness of the chosen MD model for validation.

axioms (1)
  • domain assumption The polymer system admits a field-theoretic description in which the inverse density ρ^{-1} functions as the small expansion parameter analogous to ħ in quantum field theory.
    This identification is the foundational step that converts the loop expansion into a practical calculation for homopolymers.

pith-pipeline@v0.9.0 · 5474 in / 1263 out tokens · 39833 ms · 2026-05-09T18:37:41.820358+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

34 extracted references · 34 canonical work pages

  1. [1]

    C. P. Brangwynne, P. Tompa, and R. V . Pappu, Polymer physics of intracellular phase transitions, Nat. Phys.11, 899 (2015)

    work page 2015

  2. [2]

    S. F. Banani, H. O. Lee, A. A. Hyman, and M. K. Rosen, Biomolecular condensates: organizers of cellular biochemistry, Nat. Rev. Mol. Cell Biol.18, 285 (2017)

    work page 2017

  3. [3]

    Berry, C

    J. Berry, C. P. Brangwynne, and M. Haataja, Physical princi- ples of intracellular organization via active and passive phase transitions, Rep. Prog. Phys.81, 046601 (2018)

    work page 2018

  4. [4]

    J.-M. Choi, A. S. Holehouse, and R. V . Pappu, Physical princi- ples underlying the complex biology of intracellular phase tran- sitions, Annu. Rev. Biophys.49, 107 (2020)

    work page 2020

  5. [5]

    E. W. Martin, A. S. Holehouse, I. Peran, M. Farag, J. J. Incicco, A. Bremer, C. R. Grace, A. Soranno, R. V . Pappu, and T. Mit- tag, Valence and patterning of aromatic residues determine the phase behavior of prion-like domains, Science367, 694 (2020)

    work page 2020

  6. [6]

    Bremer, M

    A. Bremer, M. Farag, W. M. Borcherds, I. Peran, E. W. Martin, R. V . Pappu, and T. Mittag, Deciphering how naturally occur- ring sequence features impact the phase behaviours of disor- dered prion-like domains, Nat. Chem.14, 196 (2022)

    work page 2022

  7. [8]

    Rizuan, S

    A. Rizuan, S. Rekhi, Y . C. Kim, S. Najafi, J.-E. Shea, and J. Mit- tal, Computational modeling of biomolecular phase separation: Current progress and open challenges, Annu. Rev. Phys. Chem. 77, 557 (2026)

    work page 2026

  8. [9]

    G. H. Fredrickson, V . Ganesan, and F. Drolet, Field-theoretic computer simulation methods for polymers and complex fluids, Macromolecules35, 16 (2002)

    work page 2002

  9. [10]

    G. H. Fredrickson and K. T. Delaney,Field-Theoretic Simula- tions in Soft Matter and Quantum Fluids, International Series of Monographs on Physics (Oxford University Press, Oxford, 2023)

    work page 2023

  10. [11]

    McCarty, K

    J. McCarty, K. T. Delaney, S. P. O. Danielsen, G. H. Fredrick- son, and J.-E. Shea, Complete phase diagram for liquid–liquid phase separation of intrinsically disordered proteins, J. Phys. Chem. Lett.10, 1644 (2019)

    work page 2019

  11. [12]

    G. H. Fredrickson,The Equilibrium Theory of Inhomogeneous Polymers, International Series of Monographs on Physics, V ol. 134 (Oxford University Press, Oxford, 2006)

    work page 2006

  12. [13]

    Y .-H. Lin, J. D. Forman-Kay, and H. S. Chan, Sequence-specific polyampholyte phase separation in membraneless organelles, Phys. Rev. Lett.117, 178101 (2016)

    work page 2016

  13. [14]

    S. F. Edwards, The theory of polymer solutions at intermediate concentration, Proc. Phys. Soc.88, 265 (1966)

    work page 1966

  14. [15]

    V . Y . Borue and I. Y . Erukhimovich, A statistical theory of weakly charged polyelectrolytes: fluctuations, equation of state and microphase separation, Macromolecules21, 3240 (1988)

    work page 1988

  15. [16]

    K. A. Mahdi and M. Olvera de la Cruz, Phase diagrams of salt- free polyelectrolyte semidilute solutions, Macromolecules33, 7649 (2000)

    work page 2000

  16. [17]

    Y . O. Popov, J. Lee, and G. H. Fredrickson, Field-theoretic sim- ulations of polyelectrolyte complexation, J. Poly. Sci. B45, 3223 (2007)

    work page 2007

  17. [18]

    K. T. Delaney and G. H. Fredrickson, Theory of polyelectrolyte complexation—complex coacervates are self-coacervates, J. Chem. Phys.146, 224902 (2017)

    work page 2017

  18. [19]

    Y .-H. Lin, J. P. Brady, J. D. Forman-Kay, and H. S. Chan, Charge pattern matching as a “fuzzy” mode of molecular recog- nition for the functional phase separations of intrinsically dis- ordered proteins, New J. Phys.19, 115003 (2017)

    work page 2017

  19. [20]

    Y .-H. Lin, J. P. Brady, H. S. Chan, and K. Ghosh, A unified ana- lytical theory of heteropolymers for sequence-specific phase be- haviors of polyelectrolytes and polyampholytes, J. Chem. Phys. 10 152, 045102 (2020)

    work page 2020

  20. [21]

    Wess ´en, S

    J. Wess ´en, S. Das, T. Pal, and H. S. Chan, Analytical For- mulation and Field-Theoretic Simulation of Sequence-Specific Phase Separation of Protein-Like Heteropolymers with Short- and Long-Spatial-Range Interactions, J. Phys. Chem. B126, 9222 (2022)

    work page 2022

  21. [22]

    Doi and S

    M. Doi and S. F. Edwards,The Theory of Polymer Dynamics (Oxford University Press, Oxford, 1986)

    work page 1986

  22. [23]

    Muthukumar and S

    M. Muthukumar and S. F. Edwards, Extrapolation formulas for polymer solution properties, J. Chem. Phys.76, 2720 (1982)

    work page 1982

  23. [24]

    Ohta and A

    T. Ohta and A. Nakanishi, Theory of semi-dilute polymer solu- tions. I. Static property in a good solvent, J. Phys. A16, 4155 (1983)

    work page 1983

  24. [25]

    Yamakawa, Modern theory of polymer solutions (1971)

    H. Yamakawa, Modern theory of polymer solutions (1971)

    work page 1971

  25. [26]

    M. E. Peskin and D. V . Schroeder,An Introduction to quantum field theory(Addison-Wesley, Reading, USA, 1995)

    work page 1995

  26. [27]

    Y .-H. Lin, J. Wess´en, T. Pal, S. Das, and H. S. Chan, Numerical techniques for applications of analytical theories to sequence- dependent phase separations of intrinsically disordered pro- teins, Methods Mol. Biol.2563, 51 (2023)

    work page 2023

  27. [28]

    Stukowski, Visualization and analysis of atomistic simula- tion data with ovito–the open visualization tool, Model

    A. Stukowski, Visualization and analysis of atomistic simula- tion data with ovito–the open visualization tool, Model. Simul. Mater. Sci. Eng.18, 015012 (2009)

    work page 2009

  28. [29]

    Ruelle,Statistical Mechanics: Rigorous Results(World Sci- entific, 1999)

    D. Ruelle,Statistical Mechanics: Rigorous Results(World Sci- entific, 1999)

    work page 1999

  29. [30]

    J. A. Anderson, J. Glaser, and S. C. Glotzer, HOOMD-blue: A Python package for high-performance molecular dynamics and hard particle Monte Carlo simulations, Comput. Mater. Sci. 173, 109363 (2020)

    work page 2020

  30. [31]

    F. J. Blas, L. G. MacDowell, E. de Miguel, and G. Jackson, Vapor-liquid interfacial properties of fully flexible Lennard- Jones chains, J. Chem. Phys.129, 144703 (2008)

    work page 2008

  31. [32]

    G. L. Dignon, W. Zheng, Y . C. Kim, R. B. Best, and J. Mittal, Sequence determinants of protein phase behavior from a coarse- grained model, PLOS Comput. Biol.14, e1005941 (2018)

    work page 2018

  32. [33]

    Tesei, T

    G. Tesei, T. K. Schulze, R. Crehuet, and K. Lindorff-Larsen, Accurate model of liquid-liquid phase behavior of intrinsically disordered proteins from optimization of single-chain proper- ties, Proc. Natl. Acad. Sci. U.S.A.118, e2111696118 (2021)

    work page 2021

  33. [34]

    J. A. Joseph, A. Reinhardt, A. Aguirre, P. Y . Chew, K. O. Rus- sell, J. R. Espinosa, A. Garaizar, and R. Collepardo-Guevara, Physics-driven coarse-grained model for biomolecular phase separation with near-quantitative accuracy, Nat. Comput. Sci. 1, 732 (2021)

    work page 2021

  34. [35]

    Feito, I

    A. Feito, I. Sanchez-Burgos, I. Tejero, E. Sanz, A. Rey, R. Collepardo-Guevara, A. R. Tejedor, and J. R. Espinosa, Benchmarking residue-resolution protein coarse-grained mod- els for simulations of biomolecular condensates, PLOS Com- put. Biol.21, e1012737 (2025)

    work page 2025