{"id":"75c7b996-d74b-4b07-9c9e-ae31e577f381","arxiv_id":"2412.15753","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Oppositely charged polyelectrolytes starting from a homogeneous solution form percolated networks that coarsen as t^{1/2} with hydrodynamic interactions, not the usual t^{1/3} droplet growth.","lead":"Simulations of oppositely charged polymers show they can form a connected, web-like network instead of round droplets when the mixture starts well mixed, and this network thickens over time with a distinctive square-root growth law. This suggests that the round droplets seen in many coacervate experiments may be a result of imperfect initial mixing, and that gentle mixing should produce network-shaped condensates.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The late-time segment of the claimed ⟨q⟩ ∼ t^{-1/2} window approaches the box-size floor: in L=69.2σ, the implied wavelength reaches 50–100σ, so the asymptotic exponent is not established beyond the finite-size regime.","rationale":"The reader's CONDITIONAL verdict is appropriate, and my concern does not change it (UNCHANGED). I differ from the reader on which assumption is most load-bearing. The constancy of DP/ϕd is the stated pivot of the scaling derivation, but the derivation is a dimensional rationalization: a diffusion equation with a single length yields t^{1/2} by construction, so the decisive question is whether the observed window actually exhibits that law. That is where the data are least secure: the late-time values of ⟨q⟩ approach the box-scale fundamental mode, so the asymptotic segment of the 'unique growth law' may be a finite-size artifact. This is a data-quality concern, directly checkable by a larger-box run, rather than an internal inconsistency.\n\nI credit the paper for what is solid: the FPD/BD contrast (−1/2 vs −1/3) is a falsifiable, mechanism-relevant comparison; the morphology crossover with volume fraction (network at ϕ = 1.2%, droplets at ϕ = 0.6%; shifted threshold for N = 92) is a clean physical result; the robustness to Bjerrum length (lB = 1.1–3σ) and to added LJ attraction (ε = 1–2 kBT) supports generality; and the authors explicitly flag the transient-network limitation and the coarse-grained-solvent caveats. The remaining load is therefore not the mechanism's plausibility but the empirical reach of the exponent.\n\nSecondary observations, not load-bearing here: S(q) is computed from the total density field including counterions (Methods), so ⟨q⟩ formally mixes counterion redistribution with polymer-network coarsening; and exponents are quoted without confidence intervals from four runs. Neither settles the verdict, but both should be addressed in a larger-box replicate study or in revision.","tokens_in":23726,"tokens_out":16236,"duration_ms":148085,"concrete_test":"Run the charge-symmetric FPD case (Nc = Na = 40, lB = 2σ, ϕ = 2.3%) in a doubled box L ≈ 138σ with all other parameters fixed, recording ⟨q⟩(t) and the dense-phase volume fraction ϕd(t) (from the Voronoi local-density distribution) at intervals of ~10^2 τBD. If ⟨q⟩ continues along the t^{-1/2} guide line below ⟨q⟩σ/2π ≈ 0.02 (wavelength > 100σ) without flattening onto the 2π/L floor, and if ϕd stays within a few percent of 0.4, both the finite-size and DP-constancy concerns are resolved. If ⟨q⟩ flattens or ϕd drifts substantially, the asymptotic t^{1/2} law and its poroelastic interpretation are not established.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is the unique coarsening law ℓ ∼ t^{1/2} (⟨q⟩ ∼ t^{-1/2}) for network-forming PE coacervation, attributed to poroelastic control of viscoelastic phase separation. Two conditions must hold: (i) the measured window actually displays the asymptotic law, and (ii) the poroelastic diffusivity DP in Eq. (1) is effectively time-independent, i.e., the dense-phase volume fraction ϕd stays fixed.\n\nCondition (i) is not settled by the presented data. Fig. 2 extends to t ≈ 10^4–10^5 τBD, where ⟨q⟩σ/2π falls to ~0.01–0.02. With σ = 1 and box size L = 69.2σ, the implied length ℓ = 2π/⟨q⟩ is 50–100σ, comparable to or larger than L; the lowest nonzero wavevector of the box corresponds to qσ/2π = 1/69.2 ≈ 0.0145. The last decade of the drawn slope −1/2 therefore sits close to the fundamental lattice mode, where ⟨q⟩ is pinned by periodic wraparound of a single spanning network. The FPD/BD exponent difference is visible at intermediate times, but the asymptotic t^{1/2} law is a finite-window observation; the paper itself concedes the network is transient and its long-time fate unresolved.\n\nCondition (ii) is asserted but not demonstrated: the text says ϕd 'remains nearly constant' but shows no time-resolved ϕd(t); the scaled-structure-factor collapse in Fig. 3b covers only t = 4000–12000 τBD (a factor of 3), not the full range over which the exponent is claimed. If ϕd drifts, DP changes and the t^{1/2} scaling from Eq. (1) no longer follows. The absence of confidence intervals on the exponents and the use of only four independent runs compound the uncertainty.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper uses fluid particle dynamics (FPD) simulations with explicit electrostatics and hydrodynamic interactions, plus free-draining Brownian dynamics (BD) controls, to study phase separation of oppositely charged polyelectrolytes starting from a homogeneous mixed state. The central observation is that charge-symmetric semi-dilute solutions form a space-spanning network that coarsens as ⟨q⟩ ∼ t^{-1/2} in the presence of hydrodynamic interactions, whereas the same model in BD gives ⟨q⟩ ∼ t^{-1/3}. The authors interpret the t^{-1/2} law as viscoelastic phase separation controlled by poroelastic solvent permeation, with τα ≫ τd. They also report that charge asymmetry decelerates coarsening, that network formation requires a volume fraction above roughly 1% (lower for longer chains), and that droplet-like coacervates formed at lower volume fraction are initially irregular and only slowly become spherical. The paper explicitly acknowledges that the network is transient and that its long-time fate is unresolved.","tokens_in":24140,"tokens_out":5522,"duration_ms":53032,"significance":"If the central claim is correct, the paper significantly expands the standard droplet picture of polyelectrolyte coacervation: it identifies a protocol-dependent network morphology, a distinct growth exponent ν=1/2, and a crucial role for hydrodynamic interactions. The strength of the work lies in the simulation design: FPD with Ewald electrostatics, internal BD comparison, multiple Bjerrum lengths, error bars from four independent runs, and a structure-factor collapse supporting dynamic self-similarity. The paper also states its own limitations clearly, including the transient nature of the network and the coarse-grained treatment of water. However, the asymptotic status of the exponent and the quantitative support for the poroelastic mechanism are not yet fully established, which is why the manuscript needs revision before the claims can be accepted at face value.","major_comments":[{"comment":"The claimed asymptotic exponent ν=1/2 is not established beyond the finite-size regime. With the box size L=69.2σ given in Methods, the smallest nonzero wavevector is q_minσ/2π = 1/69.2 ≈ 0.0145. In Fig. 2a, 2c, and 2e, ⟨q⟩σ/2π reaches values of order 0.01–0.02 near t ≈ 10^4–10^5 τBD, so the last decade of the drawn −1/2 slope sits close to the fundamental lattice mode, where periodic wraparound of a single spanning network can pin the first moment. The FPD/BD exponent difference is visible at intermediate times, but the asymptotic t^{1/2} law is a finite-window observation; the manuscript itself states that the network is transient and its long-time fate is unresolved. Please provide larger-box simulations or restrict the exponent claim to a window clearly above the box-size floor, and report uncertainties on the fitted exponents rather than only error bars on ⟨q⟩.","section":""},{"comment":"The derivation of ℓ ∼ t^{1/2} from Eq. (1) rests on the assumption that D_P is time-independent because the dense-phase volume fraction ϕd remains nearly constant during coarsening. No time-resolved ϕd(t) during network coarsening is shown; the values ϕ ≈ 0.38 and 0.42 quoted in Methods come from separate bulk equilibrium simulations, not from the coarsening network. The structural self-similarity evidence in Fig. 3b covers only t = 4000–12000 τBD, a factor of 3, whereas the exponent is claimed over a much longer interval. The poroelastic mechanism is therefore plausible but not demonstrated: if ϕd drifts, or if the dense phase cannot be described by a single poroelastic diffusivity, the t^{1/2} scaling from Eq. (1) does not follow. I would like to see ϕd(t) during coarsening and, if possible, a direct test of the diffusive relaxation implied by Eq. (1), or a clearly softened claim that the mechanism is inferred by analogy with previous work.","section":""},{"comment":"The central morphological claim that the system forms a percolated, space-spanning network, and the reported morphology transition at ϕ* between 0.6% and 1.2%, are based on representative snapshots and visual inspection. No percolation probability, largest-cluster size, or connectivity order parameter is reported, despite the availability of four independent runs. Because the abstract and conclusion emphasize network formation and the crossover volume fraction, a quantitative percolation analysis is needed to support the 'space-spanning network' assertion and the claimed dependence of ϕ* on chain length.","section":""}],"minor_comments":[{"comment":"The abstract calls the growth law 'unique', but the same ν=1/2 exponent was previously reported for viscoelastic phase separation of neutral low-molecular-weight polymer solutions and colloidal suspensions (Refs. [52,55]); the novelty is the polyelectrolyte system and the persistence of self-similarity, so the wording should be adjusted to avoid overstatement.","section":""},{"comment":"The panel labels inside Fig. 2 are difficult to parse in the present version (for example the text fragments near panels a and c); please redraw the figure with clear panel labels and consistent axis annotations.","section":""},{"comment":"The Methods section states that the strain rate is obtained from the linear increase of |ε_bulk| and |ε_shear| with δt, but it does not specify the range of δt used for the linear fit; please give this range explicitly, since the extracted τd depends on it.","section":""},{"comment":"The caption for panels g–i mentions data at lB=2σ and lB=3σ, but panel g is not annotated to show which curve corresponds to which Bjerrum length; please label the curves directly in the figure.","section":""}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Hi — quick take on Yuan and Tanaka, arXiv:2412.15753.\n\nThe paper is worth a serious look. What's genuinely new is that starting from a homogeneous mixture, oppositely charged polyelectrolytes form a space-spanning network rather than droplets, even at a few percent volume fraction, and that the network coarsens with an exponent near 1/2 when hydrodynamic interactions are included, versus 1/3 in free-draining BD. The contrast between FPD and BD is clean and important: it directly shows that hydrodynamics change the coarsening mechanism. The charge-asymmetry result — slowing due to net surface charge buildup, not chain stretching — is also a meaningful step beyond Chen and Wang. The authors are honest about the limits: the network is transient and the long-time fate is unresolved.\n\nWhere I'd push back: the exponent. The late-time segment of the claimed t^{-1/2} window sits right at the box-size floor. In a box of L=69.2σ, the fundamental mode is qσ/2π = 0.0145, and the data go down to about 0.01–0.02. That means the last decade or so of the drawn slope is close to where periodic wraparound pins the first moment. The asymptotic law is not established beyond the finite-size regime. The paper's own words — 'transient network' — concede that. Second, the poroelastic scaling derivation assumes D_P is time-independent because φ_d stays roughly constant, but there is no time-resolved φ_d(t) shown. The structure-factor collapse covers only a factor of three in time. That's a gap, not a fatal one. The 1/2 exponent is plausible, and the FPD/BD difference at intermediate times supports the mechanism, but the supporting evidence is weaker than the abstract's tone.\n\nMinor: only four independent runs, no confidence intervals on exponents, no code/data deposit. The extrapolation to experimental droplets as 'imperfect mixing' artifacts is speculative but clearly labeled as a suggestion, not a conclusion.\n\nWho should read this: anyone working on coacervate dynamics, viscoelastic phase separation, or condensate morphology. It deserves referee time — the question of protocol-dependent morphology is important and the simulation evidence is substantial even if the asymptotic claim is not airtight. I'd recommend a proper review, with the authors asked to address the finite-size issue and provide φ_d(t) if they can.","headline":"A solid simulation study with a genuinely new initial-condition story, but the headline t^{1/2} coarsening law is not as firmly established as the abstract implies.","tokens_in":24695,"tokens_out":2322,"would_cite":true,"duration_ms":21237,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["64.75.Gh","82.35.Rs"],"model":"deepseek-v4-flash","headline":"Oppositely charged polyelectrolytes coacervate into a percolating network, not droplets, when phase separation starts from a homogeneous solution.","keywords":["coacervation","polyelectrolyte complex","viscoelastic phase separation","network formation","domain coarsening","hydrodynamic interactions","charge asymmetry","interfacial tension"],"falsifier":"Measure the dense-phase volume fraction and the first moment $\\langle q \\rangle$ of the structure factor during salt-jump coarsening of a charge-symmetric polyelectrolyte solution: if $\\langle q \\rangle$ decays as $t^{-1/3}$ rather than $t^{-1/2}$, or if $\\phi_d$ changes by more than a few percent over the coarsening window, the poroelastic time-independent-$D_P$ derivation fails.","tokens_in":23505,"feed_emoji":"🕸️","tokens_out":6755,"duration_ms":51226,"temperature":0.7,"pith_summary":"This paper argues that the standard picture of coacervates as spherical droplets obeying classical liquid–liquid phase separation is incomplete. Starting from a thoroughly mixed semi-dilute solution of symmetric oppositely charged polyelectrolytes, the simulations show spontaneous formation of a space-spanning network whose characteristic length grows as $t^{1/2}$, an exponent distinct from the $t^{1/3}$ droplet law. The mechanism is viscoelastic phase separation: the polymer-rich phase relaxes more slowly than the domain deforms, and coarsening is limited by solvent permeation through the dense network. The paper also shows that charge asymmetry slows coarsening through electrostatic repulsion at the network surface, and that droplets formed at low concentration are irregularly shaped because interfacial tension is weak in good solvents. The result matters because it suggests observed droplet morphologies may often be artifacts of imperfect initial mixing, and because it connects coacervate formation to the physics of gels and porous materials.","feed_headline":"Polyelectrolyte coacervates coarsen as networks with a square-root law","feed_subtitle":"Simulations show charge-matched polyelectrolyte pairs form transient gels, not droplets; charge asymmetry can stabilize them.","key_machinery":"The central object is viscoelastic phase separation (VPS), the dynamically asymmetric regime in which the dense phase's structural relaxation time $\\tau_\\alpha$ is much longer than the domain deformation time $\\tau_d$, so the slow phase cannot follow the fast deformation and a transient network forms instead of droplets. The rate-controlling step is solvent permeation through the dense network, captured by the poroelastic equation $\\partial \\varepsilon/\\partial t = D_P \\nabla^2 \\varepsilon$, whose scaling analysis yields the characteristic length $\\ell \\sim t^{1/2}$ provided the poroelastic diffusivity $D_P$ can be treated as constant because the dense-phase volume fraction $\\phi_d$ barely changes during coarsening. The simulations compare models with hydrodynamic interactions (fluid particle dynamics) and free-draining Brownian dynamics, and the presence of the $t^{1/2}$ versus $t^{1/3}$ difference is used to show that hydrodynamic interactions are essential to this mechanism.","core_discovery":"On its own terms, the central discovery is that oppositely charged polyelectrolytes, when mixed from a homogeneous semi-dilute solution and allowed to phase separate, first form a transient percolated network that coarsens self-similarly with growth exponent $\\nu = 1/2$ in the presence of hydrodynamic interactions, rather than the classical $\\nu = 1/3$ of droplet coarsening. The exponent is the same as the mechanical-relaxation-limited coarsening seen in viscoelastic phase separation of neutral polymer solutions and colloidal suspensions, and it is traced to a poroelastic diffusion equation with a time-independent diffusivity because the dense-phase volume fraction stays nearly constant. A distinctive electrostatic ingredient is that the attractions between polycations and polyanions in good solvents are weak and long-ranged, arising from spatial charge inhomogeneity under global charge neutrality, which keeps the dense phase loosely packed and interfacial tension very low. Under charge asymmetry, net surface charge accumulates and decelerates coarsening, eventually causing dynamic slowing down of electrostatic origin.","pith_inferences":["By the paper's logic, bidisperse or partially charged polyelectrolyte mixtures in cells should be expected to form transient networks whenever chain dynamics is slow relative to local deformation; this could explain mesh-like condensates observed in vivo without invoking specific cross-linkers.","A direct test of the $1/2$ law would be a stopped-flow salt-jump experiment on a charge-matched polyelectrolyte pair, tracking $\\langle q \\rangle$ from small-angle scattering: observing $t^{-1/2}$ over two or more decades would confirm the poroelastic mechanism, while $t^{-1/3}$ would indicate classical droplet coarsening.","If hydrodynamic interactions are indeed required for the $t^{1/2}$ law, then experiments in highly viscous or confined environments (where hydrodynamic coupling is screened) should show a crossover back toward the $t^{-1/3}$ exponent, a prediction not stated in the paper.","The claim that droplets are kinetically trapped in irregular shapes suggests that coacervate 'roundness' could be used as a proxy for interfacial tension, and that measured interfacial tensions from droplet shape relaxation may systematically overestimate equilibration."],"forward_implications":["Spherical coacervate droplets seen in many experiments may partly result from imperfect initial mixing; initiation from a homogeneous state produces transient networks instead.","The coarsening exponent for network-forming coacervates should be $t^{1/2}$ with hydrodynamic interactions and $t^{1/3}$ without them, a testable signature.","Charge asymmetry provides a control knob: a small excess of one chain length slows coarsening and can stabilize the network against breakup, relevant for porous hydrogels and mesh-like biological condensates.","Because the dense phase is loose and interfacial tension is low, droplets are irregular and only slowly round; in good solvents they may remain nonspherical for very long times.","Longer polyelectrolyte chains lower the volume fraction needed for network formation, so percolated coacervates should be common at low concentration for long polymers."],"supporting_citations":[{"why":"Provides the viscoelastic phase separation concept used to interpret network formation.","marker":"[27]"},{"why":"Establishes the mechanical-slowing-down mechanism and the $t^{1/2}$ coarsening law for neutral polymer networks that this paper extends to polyelectrolytes.","marker":"[52]"},{"why":"Supplies the poroelastic scaling derivation for power-law coarsening governed by mechanical relaxation.","marker":"[55]"},{"why":"Gives the droplet coacervate baseline ($t^{-1/3}$ and charge-asymmetry slowing) that this paper contrasts with $t^{-1/2}$ network growth.","marker":"[17]"},{"why":"The fluid particle dynamics method used to incorporate hydrodynamic interactions.","marker":"[63]"},{"why":"Poroelasticity theory behind Eq. (1), the diffusion equation for local volume deformation.","marker":"[57]"}],"fun_headline_variants":["Coacervates form transient networks, not droplets, with t^1/2 coarsening","Charge-matched polyelectrolytes coarsen as gels with square-root growth","Network coarsening in coacervates: exponent 1/2, not 1/3","Percolating polyelectrolyte networks challenge droplet view of coacervates","Weak interfacial tension makes coacervate droplets irregular at first"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The scaling derivation treats the poroelastic diffusivity of the dense phase as time-independent, which requires the dense-phase volume fraction to remain nearly constant while the network coarsens; if that density drifts, the $t^{1/2}$ law and the mechanical-relaxation mechanism do not follow.","fun_headline_variants_meta":{"raw":{"variants":["Coacervates form transient networks, not droplets, with t^1/2 coarsening","Charge-matched polyelectrolytes coarsen as gels with square-root growth","Network coarsening in coacervates: exponent 1/2, not 1/3","Percolating polyelectrolyte networks challenge droplet view of coacervates","Weak interfacial tension makes coacervate droplets irregular at first"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000534,"raw_usage":{"total_tokens":2579,"prompt_tokens":966,"completion_tokens":1613,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":582,"completion_tokens_details":{"reasoning_tokens":1502}},"tokens_in":582,"tokens_out":1613,"duration_ms":9568,"temperature":1.0,"reasoning_tokens":1502,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T11:08:01.299187+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the dense-phase volume fraction and the first moment $\\langle q \\rangle$ of the structure factor during salt-jump coarsening of a charge-symmetric polyelectrolyte solution: if $\\langle q \\rangle$ decays as $t^{-1/3}$ rather than $t^{-1/2}$, or if $\\phi_d$ changes by more than a few percent over the coarsening window, the poroelastic time-independent-$D_P$ derivation fails.","supporting_citations":[{"cited_title":"Viscoelastic phase separation,","cited_arxiv_id":null,"evidence_quote":"Provides the viscoelastic phase separation concept used to interpret network formation."},{"cited_title":"Me- chanical slowing down of network-forming phase separa- tion of polymer solutions,","cited_arxiv_id":null,"evidence_quote":"Establishes the mechanical-slowing-down mechanism and the $t^{1/2}$ coarsening law for neutral polymer networks that this paper extends to polyelectrolytes."},{"cited_title":"Power-law coars- ening in network-forming phase separation governed by mechanical relaxation,","cited_arxiv_id":null,"evidence_quote":"Supplies the poroelastic scaling derivation for power-law coarsening governed by mechanical relaxation."},{"cited_title":"Charge asym- metry suppresses coarsening dynamics in polyelectrolyte complex coacervation,","cited_arxiv_id":null,"evidence_quote":"Gives the droplet coacervate baseline ($t^{-1/3}$ and charge-asymmetry slowing) that this paper contrasts with $t^{-1/2}$ network growth."},{"cited_title":"Simulation method of colloidal suspensions with hydrodynamic interactions: Fluid particle dynamics,","cited_arxiv_id":null,"evidence_quote":"The fluid particle dynamics method used to incorporate hydrodynamic interactions."},{"cited_title":"General theory of three-dimensional consolidation,","cited_arxiv_id":null,"evidence_quote":"Poroelasticity theory behind Eq. (1), the diffusion equation for local volume deformation."}],"review_version":1}