{"id":"40aa6f34-c226-449a-95c4-6b3d3b510eb0","arxiv_id":"2501.05061","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Exact equilibrium droplet and electrostatic energy are derived for the pre-critical (overlapping caps) phase of the two-charge spherical Coulomb gas.","lead":"This paper computes the exact shape of the charged-particle cloud on a sphere when two heavy fixed charges are close enough on the surface that their influence zones overlap. It provides the energy of that phase and links the result through duality to a one-dimensional random matrix ensemble.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The pre-critical droplet boundary is assumed, not derived, to be the rational map (2.46); if that classification does not cover the asymmetric two-charge field, the central claim fails.","rationale":"The reader's weakest assumption is exactly the same point: the rational conformal map (2.46) and the Stieltjes-transform representation (2.48) are imported from [23,25] rather than derived here. This is indeed the single most load-bearing step, because the parameter equations and the energy formula are all downstream of that representation. The paper's internal checks (the symmetric Q0=Q1 ellipse limit, the Q0→∞ planar limit, and the critical-boundary matching of K_pre with K_post in Remark 2.6.4) provide substantial evidence but do not cover the general asymmetric case. The computer-algebra reduction to the quartic is a secondary but real verifiability gap: Appendix B's discriminant argument is not written coherently, so the claim that the quartic has exactly two positive roots and that the smallest one is the correct branch rests on an unshown computation. A related minor slip is Remark 2.6.1, where the preimage of Z=w is stated as u_Z=w instead of u_Z=v0; the surrounding calculations use v0, so this appears typographical. None of this establishes a contradiction, and the numerical figure and critical matching suggest the formulas are likely correct. The appropriate response is therefore to keep the reader's CONDITIONAL verdict and require the numerical or independent-symbolic verification described above before upgrading to ACCEPT.","tokens_in":33117,"tokens_out":20117,"duration_ms":195244,"concrete_test":"Run a numerical equilibrium solver for the asymmetric pre-critical case, e.g. Q0=4, Q1=2, w=0.3 and one value closer to w_cri: discretize the energy functional (1.10), minimize over the droplet support on a fine grid, and compare the computed support boundary with ζ(|u|=1) using α from the smallest positive root of (2.63). Then evaluate the energy of the numerical equilibrium measure and compare it with (2.89), using u_Z=v0 in (2.101) for Z=w. Agreement to a few percent would confirm the imported parametrization; disagreement would break the central claim.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim depends on Eq. (2.46): the droplet boundary in the pre-critical phase is taken to be the image of the unit circle under ζ(u) = R/u (1-bu)/(1-au), with [25, Lemma 2] then used to write the Cauchy transform (2.48). This parametrization is imported from [23,25] and is not re-derived for the two-charge spherical potential. The paper verifies two limits (Q0=Q1 giving the known ellipse, and Q0→∞ giving the planar Ginibre result), but for general Q0≠Q1 there is no independent check that the true equilibrium support has this rational form. If the classification hypotheses fail, equations (2.55)–(2.63) describe a different measure and the energy (2.89) is not the pre-critical free energy. In addition, the reduction to the quartic (2.63) and the root count in Appendix B are delegated to computer algebra, and Appendix B's discriminant text is internally garbled ('roots of h are all positive'), so the choice of the smallest positive root is not human-verifiable as written. This is an addressable gap rather than a demonstrated contradiction, but it is the most load-bearing point in the paper.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the one-component Coulomb gas on the sphere with N mobile unit charges and two macroscopic external charges of strengths Q0N and Q1N. In the post-critical phase, when the two associated spherical caps do not overlap, the equilibrium measure and free energy are already known, and the paper relies on those results. The main new contribution is the pre-critical (overlapping-caps) phase: the droplet boundary is represented by the conformal map ζ(u) = R/u (1-bu)/(1-au), with the parameters R, a, b expressed through the smallest positive root α of the quartic (2.63), and the leading-order electrostatic energy is evaluated in closed form as K_pre^N in (2.89) via Propositions 2.8 and 2.9. The paper also uses a random-matrix duality to prove the post-critical expansion of Proposition 1.1, to compute the critical scaling regime of Corollary 3.1, and to derive an identity, (4.10), between the pre-critical energy difference and a constrained Jacobi unitary ensemble rate function. The derivation contains explicit cross-checks against known limits: the equal-charge ellipse of §2.2, the large-w disk limit (2.75), the phase-boundary equality (2.108), and the Ginibre limit of Remark 2.4.","tokens_in":33333,"tokens_out":8072,"duration_ms":80640,"significance":"If the derivation is correct, the paper supplies the exact equilibrium measure and the exact leading-order electrostatic energy for a nontrivial two-charge phase transition on the sphere, including a parameter-free conformal-map description and closed-form energy integrals. The manuscript has several concrete strengths: explicit formulas throughout, consistency with the known Q0=Q1 ellipse (Proposition 2.6 and Corollary 2.1), agreement in the large-w limit (2.75) and in the Q0→∞ planar Ginibre limit, the phase-boundary matching (2.108), and numerical support in Figures 1, 3, and 5. However, the central ansatz (2.46) is imported from the classification results in [23,25] rather than derived for the two-charge spherical potential, and the root selection for the quartic (2.63) is only justified through computer algebra, with Appendix B containing an internally inconsistent sentence. These issues make the significance conditional: the paper is likely correct in its main claims, but the load-bearing classification and root-selection steps need to be made verifiable before the results can be accepted as established.","major_comments":[{"comment":"The rational-form ansatz for the droplet boundary is the load-bearing assumption of the paper. The manuscript states that (2.46) follows from [23, Eq. (27) with N=1] and [25, Th. 5], and it uses [25, Lemma 2] to write the Stieltjes transform (2.48); however, it does not verify that the hypotheses of these classification theorems hold for the general Q0 ≠ Q1 two-charge spherical field. The checks for Q0=Q1 (§2.2, Proposition 2.6, Corollary 2.1) and for the Q0→∞ Ginibre limit (Remark 2.4) are strong consistency tests, but they do not exclude the possibility that the true support for generic Q0, Q1, w lies outside this rational family. Because equations (2.55)–(2.63) and the energy (2.89) are all derived from (2.46), this gap is directly load-bearing for the central claim. I recommend adding either a verification that the [23,25] classification applies in this setting, an independent derivation of the two-pole/one-zero form from the equilibrium equation (2.9), or a numerical check that the predicted boundary satisfies (2.9) for a representative non-symmetric parameter set.","section":"§2.3, Eq. (2.46)"},{"comment":"The reduction of the residue-matching equation (2.61) to the quartic (2.63) and the subsequent root counting are delegated to computer algebra, and Appendix B contains a contradictory sentence: after arguing that h(x) has one negative root and that the remaining two roots form a complex conjugate pair, it concludes that “roots of h are all positive.” The needed conclusion is the opposite—that h has no positive roots, so Disc_α(p) ≠ 0 and all roots of p are real and distinct. As written, Appendix B does not support the “smallest positive root” selection in Proposition 2.7, and the large-w expansion (2.72) alone does not identify the root unless the non-crossing property and the positivity of the other root are established. I recommend replacing the CAS-only statements with an explicit discriminant computation or a corrected, self-consistent argument, and stating precisely which criterion selects α.","section":"§2.3 and Appendix B"}],"minor_comments":[{"comment":"The sentence “which justifies (2.4) by continuity of a, b with respect to w” appears to refer to the inequality 0 < a < b stated in (2.47), not to equation (2.4); the cross-reference should be corrected.","section":"Remark 2.4.2"},{"comment":"In addition to the substantive issue raised above, the sentence “Thus we have the required result that roots of h are all positive” is a typographical inversion of the required statement, and it should be rewritten to say that h has no positive roots.","section":"Appendix B"},{"comment":"The phrase “determinisation of the electrostatic energy” should read “determination of the electrostatic energy.”","section":"§1.3"},{"comment":"The notation α(Q1, Q2) in the text before equation (3.18) is inconsistent with the definition α(Q0, Q1) in (3.20) and should be corrected.","section":"§3.4"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the scope of the journal and I see no novelty-disclosure concern. The main risk is the imported classification of the droplet boundary; if the authors can supply the missing verification or a corrected, human-verifiable root-count argument, the paper would be a solid contribution. I do not think the issues rise to rejection, but they are load-bearing enough to require a revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Bottom line: this is a serious paper that delivers the exact equilibrium measure and electrostatic energy for the pre-critical (overlapping caps) phase of the two-charge spherical Coulomb gas for general Q0≠Q1. The result is new and the internal consistency checks are convincing. I would send it to review, with a request to address the parametrization assumption and clean up the computer-algebra endpoints. What's actually new: the quartic (2.63) specifying the conformal map parameter α, and the exact energy K_pre (2.89)–(2.101). The equal-charge ellipse and the post-critical phase were known; the asymptotic checks against those limits, plus the Q0→∞ Ginibre limit and the large-w disk limit, are the right tests and they pass. The numerics in Figures 3 and 5 support the phase-boundary and the energy identity (4.10). The duality-based partition function expansion in §3 is a nice application, and the relation (4.10) between the sphere energy and a constrained Jacobi log-gas is a genuine byproduct. Soft spots: the rational form of the droplet boundary, (2.46), is assumed from [23,25] rather than derived for this two-charge problem. That is the load-bearing point. The paper does test the assumption in important limits, and the limits are consistent, so I don't think it's wrong, but a referee should ask the authors to make the domain of validity of the classification explicit, or to provide a self-contained derivation. The reduction to the quartic and the discriminant analysis in Appendix B sit on computer algebra; that's common in this area and acceptable, but the text is not human-verifiable as written. Also, Appendix B contains a typo that should be fixed: it says 'roots of h are all positive' where it clearly means no positive roots. The duality identity used for the large-N results is cited to an unpublished preprint [29]; not a fatal issue, but it should be stated which parts depend on that and whether a published reference exists. Overall: the central claim looks right, the gaps are addressable, and the paper deserves a serious referee. It will be useful to people working on log-gases, potential theory on the sphere, and Coulomb gas phase transitions.","headline":"Solid exact solution of the pre-critical droplet for the two-charge spherical Coulomb gas; the main caveat is the imported rational-form assumption for the droplet boundary and the heavy use of computer algebra, both addressable.","tokens_in":33896,"tokens_out":2421,"would_cite":true,"duration_ms":23457,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["31A15","60B20","82B05","30C20"],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper establishes the exact equilibrium measure and leading electrostatic energy of the one-component Coulomb gas on a sphere in the pre-critical phase, where the two external-charge caps overlap.","keywords":["one-component Coulomb gas","spherical ensemble","equilibrium measure","conformal map","pre-critical phase","electrostatic energy","Jacobi unitary ensemble","duality identity"],"falsifier":"Numerically minimise the energy functional (1.10) for parameter triples $(Q_0,Q_1,w)$ in the overlapping regime, compare the boundary of the minimiser with the curve $\\zeta(e^{i\\theta})$ obtained by solving (2.63) and using (2.54)-(2.60), and check the energy formula (2.89) by independent quadrature; any systematic discrepancy would falsify the claimed equilibrium measure and energy. Alternatively, compute both sides of the identity (4.10) numerically for several triples: the left side uses the droplet energies, the right side uses the explicit rate function (4.8).","tokens_in":32905,"feed_emoji":"⚛️","tokens_out":7585,"duration_ms":64002,"temperature":0.7,"pith_summary":"The paper studies the one-component Coulomb gas on a sphere: $N$ mobile unit charges repel logarithmically in the presence of two fixed charges of strengths $Q_0 N$ and $Q_1 N$. When the spherical caps naturally associated with the fixed charges overlap, earlier post-critical formulas cease to apply; this paper determines the equilibrium measure in that pre-critical phase. It asserts that the droplet boundary is the image of the unit circle under the rational conformal map $\\zeta(u)=\\frac{R}{u}\\frac{1-bu}{1-au}$, with all parameters fixed by the smallest positive root of a quartic, and derives the exact leading-order electrostatic energy. The payoff is that the exact free energy of the Coulomb gas is now known on both sides of the phase transition, and for $\\beta=2$ the partition function connects to a Jacobi unitary ensemble gap probability via a duality identity.","feed_headline":"Exact droplet boundary found for overlapping charges on a sphere","feed_subtitle":"A rational conformal map fixed by one quartic root gives the equilibrium measure and leading energy in the pre-critical phase.","key_machinery":"The central object is the rational conformal map $\\zeta(u)$ from the unit disk to the exterior of the droplet, whose form (2.46) is imported from the classification results [23,25]. The derivation runs through the Stieltjes transform representation (2.48): requiring that the interior pole $v_0$ of the first term satisfies $\\zeta(v_0)=w$ and that its residue cancels the pole of the charge term yields (2.55)-(2.57), reducing everything to the quartic (2.63). The energy computation is carried by two explicit contour-integral evaluations: Proposition 2.8 deforms a logarithmic contour around the branch point $v_0$ to compute $\\int_{\\tilde\\Omega_d}\\log(1+|z|^2)\\mu(z)\\,d^2z$, and Proposition 2.9 integrates the Stieltjes transform (2.100) to compute the logarithmic potential $W(Z)$. For the partition function at $\\beta=2$, the load-bearing identity is the duality (3.7)/(3.11), which equates the spherical Coulomb-gas configuration integral with a gap probability of the Jacobi unitary ensemble; the Wachter density (3.13) and its constrained analogue (4.6) then deliver the post-critical, critical, and pre-critical asymptotics.","core_discovery":"In the pre-critical phase, specified by $w>w_{\\rm cri}$ in (2.74), the equilibrium measure is supported on the droplet whose boundary is the image of $|u|=1$ under $\\zeta(u)=\\frac{R}{u}\\frac{1-bu}{1-au}$, with $0<a<b$, $a=R\\alpha$, $b=\\beta/R$, $\\beta=(1+Q_1)\\alpha/Q_0$, and $R^2$ given by (2.56). The parameter $\\alpha$ is the smallest positive root of the quartic (2.63), selected by its large-$w$ behaviour. The authors fix these parameters by imposing pole and residue cancellation in the Stieltjes transform representation (2.48), and they evaluate the background, one-body, and two-body energies to obtain the closed-form pre-critical Boltzmann constant $K_{\\rm pre}^N$ via (2.89), (2.90), and (2.101). At the phase boundary the pre-critical energy matches the post-critical value (1.5). For $\\beta=2$, the duality identity (3.7) yields the large-$N$ expansion (1.8) in the post-critical phase, the soft-edge Painlev\\'e II correction (3.21) in the critical scaling window, and, in the pre-critical phase, the electrostatic-energy identity (4.10).","pith_inferences":["Beyond the paper: the quartic-root parameterization should be verifiable by direct numerical minimisation of the energy functional (1.10); agreement of the support boundary would confirm the imported classification input.","Beyond the paper: the explicit formula for $K_{\\rm pre}^N$ is a natural anchor for computing subleading $1/N$ corrections in the pre-critical phase, which the paper defers.","Beyond the paper: the Appendix C comparison suggests the same energy formulas may transfer to truncated-unitary and Poincar\\'e-disk Coulomb systems, where a pre-critical energy identity of similar shape should hold.","Beyond the paper: the contour-integral technique for $W(Z)$ may extend to three or more external charges once the corresponding rational conformal-map classification is available."],"forward_implications":["The pre-critical droplet is exactly determined for general $Q_0\\neq Q_1$, extending the equal-charge ellipse result, and it reproduces the centred disk of radius $1/\\sqrt{Q_0+Q_1}$ as $w\\to\\infty$.","The closed-form $K_{\\rm pre}^N$ gives the leading $N^2$ free energy in the overlapping-caps phase, matching the post-critical constant at the phase boundary.","For $\\beta=2$, the post-critical partition function is independent of the charge position $w$ to all inverse powers of $N$, while the critical window adds $E_2^{\\rm soft}(0;(Q_0^{-2/3}s,\\infty))$.","In the pre-critical phase the duality identity forces (4.10), equating the difference of sphere-system energies with the constrained Jacobi large-deviation rate function.","The identity (4.10) lets the leading pre-critical gap probability of the Jacobi unitary ensemble be read off from the sphere droplet energy."],"supporting_citations":[{"why":"Supplies the rational form (2.46) of the droplet boundary conformal map used throughout.","marker":"[23]"},{"why":"Supplies the classification and the Stieltjes transform representation (2.48) (Lemma 2) that fix the parameter equations.","marker":"[25]"},{"why":"Treats the equal-charge symmetric case whose ellipse result is reclaimed by Proposition 2.2 and used as a consistency check.","marker":"[21]"},{"why":"Provides the planar analogue with a macroscopic charge, including the conformal-map parameterization and free-energy expansions that guide the pre-critical derivation.","marker":"[17]"},{"why":"Gives the post-critical Boltzmann factor (1.5) and the spherical Coulomb-gas calculation that the new pre-critical constant must match.","marker":"[26]"},{"why":"Supplies the duality identity (3.4)/(3.7) that connects the spherical ensemble to the Jacobi ensemble gap probability.","marker":"[29]"},{"why":"Provides the log-gas/Coulomb-gas formalism, stereographic projection, and normalisation identities underlying the energy computation.","marker":"[27]"},{"why":"Gives the constrained Jacobi spectral density (4.6) and the large-deviation rate function used in the pre-critical energy identity (4.10).","marker":"[44]"}],"fun_headline_variants":["Overlapping charges: exact droplet boundary on sphere","Quartic root determines droplet for overlapping charges","Exact energy for pre-critical Coulomb gas on sphere","Sphere Coulomb gas: overlapping-phase droplet solved","Pre-critical droplet from quartic root and conformal map"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the droplet boundary has the rational conformal-map form (2.46) imported from [23,25]; if the true equilibrium measure in the overlapping-caps phase is not of that form, the parameter equations (2.55)-(2.63) would not describe it, and the reduction to the quartic relies on computer algebra rather than a fully human-verifiable derivation.","fun_headline_variants_meta":{"raw":{"variants":["Overlapping charges: exact droplet boundary on sphere","Quartic root determines droplet for overlapping charges","Exact energy for pre-critical Coulomb gas on sphere","Sphere Coulomb gas: overlapping-phase droplet solved","Pre-critical droplet from quartic root and conformal map"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000619,"raw_usage":{"total_tokens":2917,"prompt_tokens":1038,"completion_tokens":1879,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":654,"completion_tokens_details":{"reasoning_tokens":1807}},"tokens_in":654,"tokens_out":1879,"duration_ms":12708,"temperature":1.0,"reasoning_tokens":1807,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T21:20:08.123342+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Numerically minimise the energy functional (1.10) for parameter triples $(Q_0,Q_1,w)$ in the overlapping regime, compare the boundary of the minimiser with the curve $\\zeta(e^{i\\theta})$ obtained by solving (2.63) and using (2.54)-(2.60), and check the energy formula (2.89) by independent quadrature; any systematic discrepancy would falsify the claimed equilibrium measure and energy. Alternatively, compute both sides of the identity (4.10) numerically for several triples: the left side uses the droplet energies, the right side uses the explicit rate function (4.8).","supporting_citations":[{"cited_title":"Crowdy and M","cited_arxiv_id":null,"evidence_quote":"Supplies the rational form (2.46) of the droplet boundary conformal map used throughout."},{"cited_title":"Etingof and X","cited_arxiv_id":null,"evidence_quote":"Supplies the classification and the Stieltjes transform representation (2.48) (Lemma 2) that fix the parameter equations."},{"cited_title":"Criado del Rey and A.B.J","cited_arxiv_id":null,"evidence_quote":"Treats the equal-charge symmetric case whose ellipse result is reclaimed by Proposition 2.2 and used as a consistency check."},{"cited_title":"Fischmann and P.J","cited_arxiv_id":null,"evidence_quote":"Gives the post-critical Boltzmann factor (1.5) and the spherical Coulomb-gas calculation that the new pre-critical constant must match."},{"cited_title":"Forrester,Dualities in random matrix theory, preprint","cited_arxiv_id":null,"evidence_quote":"Supplies the duality identity (3.4)/(3.7) that connects the spherical ensemble to the Jacobi ensemble gap probability."},{"cited_title":"Ramli, E","cited_arxiv_id":null,"evidence_quote":"Gives the constrained Jacobi spectral density (4.6) and the large-deviation rate function used in the pre-critical energy identity (4.10)."}],"review_version":1}