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

Recognition: unknown

Loop expansion in polymer field theory: application to phase separation

Kiyoharu Kawana , Kyosuke Adachi

Authors on Pith no claims yet

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 comments

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: FIG. 1. A graphical representation of the connected correlation view at source ↗

**Figure 2.** Figure 2: FIG. 2. The LO diagrams contributing to the free energy. view at source ↗

**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.

Desk Editor's Note private letter to a colleague

The paper maps 1/ρ to ħ to organize loop corrections around RPA for homopolymer binodals, shows the first correction moves the dilute coexistence density closer to MD, but the expansion parameter is not small in the regime where that shift occurs.

read the letter

The core of this paper is a field-theoretic loop expansion for homopolymer phase separation. They identify the inverse density ρ^{-1} as ħ, treat the RPA free energy as the tree-level term, and compute the one-loop (RPA+) and two-loop (RPA++) corrections. When they extract binodals from the RPA+ free energy and compare to bead-spring MD simulations with Gaussian interactions, the dilute-phase density improves qualitatively while the critical-point location stays roughly as inaccurate as plain RPA.

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)

[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.
[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)

Define the precise meaning of RPA+ and RPA++ (free-energy corrections or effective potentials) at first use and ensure consistent notation throughout.
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

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
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)

Reference graph

Works this paper leans on

34 extracted references

[1]

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

2015
[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)

2017
[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)

2018
[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)

2020
[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)

2020
[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)

2022
[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)

2026
[9]

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

2002
[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)

2023
[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)

2019
[12]

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

2006
[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)

2016
[14]

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

1966
[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)

1988
[16]

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

2000
[17]

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

2007
[18]

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

2017
[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)

2017
[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)

2020
[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)

2022
[22]

Doi and S

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

1986
[23]

Muthukumar and S

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

1982
[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)

1983
[25]

Yamakawa, Modern theory of polymer solutions (1971)

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

1971
[26]

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

1995
[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)

2023
[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)

2009
[29]

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

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

1999
[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)

2020
[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)

2008
[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)

2018
[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)

2021
[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)

2021
[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)

2025