[Paper Review] Convergence of the Mayer Series via Cauchy Majorant Method with Application to the Yukawa Gas in the Region of Collapse

Leonardo F. Guidi, Domingos H. U. Marchetti|ArXiv.org|Oct 14, 2003

Fractional Differential Equations Solutions17 references3 citations

TL;DR

This paper establishes the convergence of the Mayer series for the Yukawa gas in the region of collapse using the Cauchy majorant method applied to a nonlinear first-order PDE derived from Ursell function flow equations. It proves Benfatto, Gallavotti, and Nicolò's conjecture up to β < 8π by introducing Lagrange multipliers and improved stability bounds for non-neutral n-particle configurations, showing convergence when leading terms are omitted and an enhanced stability condition holds.

ABSTRACT

We construct majorant functions $Φ(t,z)$ for the Mayer series of pressure satisfying a nonlinear differential equation of first order which can be solved by the method of characteristics. The domain $| z|

Motivation & Objective

• To establish convergence of the Mayer series for the Yukawa gas in the region of collapse, where standard methods fail due to instability.
• To resolve BGN's conjecture regarding finite radius of convergence for β ∈ [4π, 8π) after omitting leading terms.
• To develop a PDE-based majorant method using the method of characteristics to analyze the radius of convergence.
• To improve stability bounds for non-neutral n-particle configurations to extend convergence results beyond previous limits.
• To demonstrate numerically that the improved stability condition is feasible for small n, supporting the convergence proof.

Proposed method

• Formulate a majorant pressure function Φ(t,z) satisfying a nonlinear first-order PDE derived from the Ursell integral equation.
• Apply the method of characteristics to solve the PDE, with t parameterizing scale decomposition of the Yukawa potential.
• Introduce Lagrange multipliers to modify the PDE, improving the linear coefficient in the ODE system for majorant coefficients.
• Use the Lambert W-function to express the solution for non-negative potentials, linking it to exponential generating functions of rooted trees.
• Construct majorant series for Mayer series coefficients by comparing with solutions of modified PDEs with improved stability bounds.
• Verify numerically that the improved stability condition holds for small n (n ≤ 5), particularly for non-neutral configurations with minimal energy.

Experimental results

Research questions

• RQ1 Does the Mayer series for the Yukawa gas converge in the region of collapse for β < 8π when leading terms are omitted?
• RQ2 Can the Cauchy majorant method combined with PDE techniques establish a finite radius of convergence for the Mayer series in this regime?
• RQ3 What improved stability condition on non-neutral n-particle configurations ensures convergence of the Mayer series?
• RQ4 How does the inclusion of Lagrange multipliers enhance the convergence radius compared to standard majorant methods?
• RQ5 Can numerical evidence confirm that the required stability bounds are satisfied for small n in the Yukawa gas model?

Key findings

• The Mayer series for the Yukawa gas converges for β < 8π when the first 2r terms are omitted, with the radius of convergence given by |z| < (β_{2r+1} - β)/(β_{2r+1} β e), where β_{2r+1} = 8π(1 + (2r+1)^{-1})^{-1}.
• The solution to the majorant PDE is expressed in terms of the Lambert W-function, indicating a combinatorial origin of the non-physical singularity.
• For n=3, the minimal energy of non-neutral configurations yields δ₃ ≈ 0.8837, significantly exceeding 1/3, which improves the convergence bound.
• For n=4 and n=5, the minimal energy configurations yield δ₄ ≈ 1.8837 and δ₅ ≈ 0.8121, respectively, both exceeding their respective 1/n thresholds.
• Numerical calculations confirm that the improved stability condition is satisfied for small n, supporting the convergence result in the infinite volume limit.
• The method successfully extends the convergence radius beyond previous results, proving BGN's conjecture up to β < 8π under the improved stability condition.