files / problems / simon-ionization-excess / compute / q1 / BOUNDS.md

3.6 KB · on GitHub

q1 bounds — HPS / Nam replay

Compute notes for compute/q1/. dent / residue as in AGENTS.md.

How to run

python3 replay_hps.py
python3 tighten_hps.py
python3 replay_nam_beta.py

or ./run_all.sh from this folder (invokes the three scripts if present). Needs numpy, scipy, mpmath. Writes certs/hps_replay.json, certs/hps_tight.json, certs/nam_beta.json.

Record replayed

Opened 2026-08-27: HPS 2504.18487v1 HTML, Nam 1009.2367v3 HTML, FHJN 1808.09017 HTML.

HPS §7 chain with 1.456 from FHJN Theorem 1 (not replaced). A later LT claim 1.44655 (arXiv:2403.04347) was not replayed and is not used.

printed independent enclosure valid?
b(2)∈(1.2071,1.2072)b(2)\in(1.2071,1.2072) 1.20710678118…1.20710678118\ldots yes
b(3)∈(1.1184,1.1185)b(3)\in(1.1184,1.1185) 1.11843379920…1.11843379920\ldots yes
1.11851.1185 upper on b(3)b(3) yes yes
a≤2.953a\le 2.953 then 2.962.96 a(5/2)=2.95203835…a(5/2)=2.95203835\ldots both yes
c<1.5855c<1.5855 1.58543698…1.58543698\ldots yes
a1<3.893a_1<3.893 on [β3−1,5/2][\beta_3^{-1},5/2] sup⁡=a1(5/2)=3.89949515…\sup=a_1(5/2)=3.89949515\ldots no
a1<3.90a_1<3.90 on that interval yes yes
extras 0.01340.0134, 0.1840.184, 0.01960.0196 0.013314680.01331468, 0.183615030.18361503, 0.019591740.01959174 yes
Nam β≥0.8218\beta\ge 0.8218 g(0.843)=0.82180392…g(0.843)=0.82180392\ldots yes
Nam β<0.8705\beta<0.8705 115/81−(1/2)ln⁡3=0.87044694…115/81-(1/2)\ln 3=0.87044694\ldots yes

Closed forms for b(2)b(2), b(3)b(3) match a numerical max of (1+ts−1)/(1+ts)(1+t^{s-1})/(1+t^s) on [0,1][0,1]. C1−1=2.34131154…C_1^{-1}=2.34131154\ldots and C2−1/2=2.21537683…C_2^{-1/2}=2.21537683\ldots match Lemma A.6 / (A.18), the displayed closed forms, and the three Lieb integrals for fpf_p.

HPS said the supremum of a1(x)a_1(x) on [β3−1,5/2][\beta_3^{-1},5/2] is at the left endpoint. Independently it is a minimum in the middle and the max is at x=5/2x=5/2. So 3.8933.893 is only an enclosure of the left-endpoint formula, not of their stated supremum. 3.903.90 still covers x=5/2x=5/2.

Nam g(λ)g(\lambda) at the maximiser λ≈0.8434764\lambda\approx 0.8434764 equals 0.82180662586…0.82180662586\ldots. The AM-GM equality geometry sits in the triangle, so inf⁡Wλ/(∣x∣+∣y∣)=g(λ)\inf W_\lambda/(|x|+|y|)=g(\lambda) there. A (r,s,(r,s, angle)) grid agrees to a few 10−710^{-7} and is not a lower bound. 1/0.8218=1.21684108…1/0.8218=1.21684108\ldots, which Nam rounds to 1.221.22.

What tightened (dent)

Same HPS theorem, tighter arithmetic, interval enclosures.

  1. Exact b(3)b(3) in a1a_1, not 1.11851.1185.
  2. Lieb N/Z<2+1/ZN/Z<2+1/Z for the Prop. 2.5 range Z≥4Z\ge 4, so N/Z<9/4N/Z<9/4 not 5/25/2. Then a1(x)a_1(x) on [b(3),9/4][b(3),9/4] is maximised at the left, and a1<3.892a_1<3.892.
  3. Evaluating the actual a1(x)a_1(x) is what shows (2). At x=b(3)x=b(3) the crude Z<Nβ3Z<N\beta_3 is equality.
  4. Scanning λ\lambda does not move a1a_1: HPS already minimises the left-endpoint coefficient.
  5. Intervals on every claimed decimal.

Certified inequalities (see certs/hps_tight.json):

r=λ(Nβ3)−1/3<0.5r=\lambda(N\beta_3)^{-1/3}<0.5 still holds on the contradiction set N≥β3−1ZN\ge\beta_3^{-1}Z, Z≥4Z\ge 4.

What did not

No new argument for β3\beta_3, so no dent of the leading 1.11851.1185. A second Nam trial (64 atoms on [1,10][1,10]) gives β≤0.87021\beta\le 0.87021 and sits next to Nam's own 0.87020.8702 remark; that tightens only the upper end of β\beta, which does not improve an ionization upper bound. Optimising a single power r−α1[1,L]r^{-\alpha}1_{[1,L]} recovers Nam's (α,L)=(3/2,9)(\alpha,L)=(3/2,9) and does not beat 115/81−(1/2)ln⁡3115/81-(1/2)\ln 3.