1. Anomalous moments from 1/(2pi/alpha - 1 +/- 4/sqrt3)
   Williams-Olsen 1979 (Aspden's input)   a_e = 0.00115964650448 ( -4.89 ppm vs 2023)   a_mu = 0.001165891225 (-25.3 ppm, -201.0 sigma vs Fermilab 2025)
   Rb atom recoil 2020                    a_e = 0.00115964619856 ( -5.16 ppm vs 2023)   a_mu = 0.001165890916 (-25.5 ppm, -203.0 sigma vs Fermilab 2025)
   Cs atom recoil 2018                    a_e = 0.00115964619991 ( -5.16 ppm vs 2023)   a_mu = 0.001165890917 (-25.5 ppm, -203.0 sigma vs Fermilab 2025)
   CODATA 2022                            a_e = 0.0011596461988 ( -5.16 ppm vs 2023)   a_mu = 0.001165890916 (-25.5 ppm, -203.0 sigma vs Fermilab 2025)
   Aspden's own lattice value             a_e = 0.00115964691175 ( -4.54 ppm vs 2023)   a_mu = 0.001165891637 (-24.9 ppm, -198.0 sigma vs Fermilab 2025)
   1979-83 data: a_e 0.0011596522, a_mu 0.001165924;  2023-25 data: a_e 0.00115965218059(1.3e-13), a_mu 0.001165920705(1.47e-10)
   muon in 1983: Aspden 0.00116589123 vs CERN 1979 0.001165924(8.5e-9): -3.86 sigma
   muon minus electron: Aspden 6.2447 x 10^-6, measured 6.2685 x 10^-6
   series coefficients of (alpha/pi)^n, n = 1..4:  Aspden ['0.5', '-0.32735', '0.214316', '-0.140313']  QED ['0.5', '-0.328479', '1.181241', '-1.912246']
   QED with alpha(Rb): a_e = 0.001159652180251, minus measured = -0.339 x 10^-12
   Aspden with alpha(Rb): minus measured = -5982.0 x 10^-12

2. 'Absolute gravitational potential' phi from eq. (3): phi = 2c^2 [(1+a_meas)/(1+a_Aspden) - 1]
   1979 g-factor, alpha Williams-Olsen 1979 (Aspden's input)   phi =   10.23 x 10^12 cm^2/s^2
   1979 g-factor, alpha Rb atom recoil 2020                    phi =   10.78 x 10^12 cm^2/s^2
   1979 g-factor, alpha Cs atom recoil 2018                    phi =   10.77 x 10^12 cm^2/s^2
   1979 g-factor, alpha Aspden's own lattice value             phi =   9.495 x 10^12 cm^2/s^2
   2023 g-factor, alpha Williams-Olsen 1979 (Aspden's input)   phi =   10.19 x 10^12 cm^2/s^2
   2023 g-factor, alpha Rb atom recoil 2020                    phi =   10.74 x 10^12 cm^2/s^2
   2023 g-factor, alpha Cs atom recoil 2018                    phi =   10.74 x 10^12 cm^2/s^2
   2023 g-factor, alpha Aspden's own lattice value             phi =    9.46 x 10^12 cm^2/s^2
   solar-system potential at the Earth's surface: 9.497 x 10^12 (Aspden: 9.49)
   3 K route, Aspden 1983: T = 2.7 K, m* = 0.04 m_e                phi = 10.23 x 10^12
   3 K route, T = 2.72548 K (COBE/FIRAS), m* = 0.04 m_e            phi = 10.33 x 10^12
   3 K route, T = 2.72548 K, m* = m_e/24.52 (his lattice particle) phi = 10.13 x 10^12
   Milky Way: Sun orbits at ~230 km/s, so |phi_galaxy| > v_c^2 = 529.0 x 10^12 = 55.7 x the solar-system value
   expected pull toward the galactic centre v^2/R0 = 2.09 x 10^-10 m/s^2; Gaia EDR3 measured 2.32 +/- 0.16 x 10^-10 m/s^2 toward the galactic centre
   annual change of phi from Earth's eccentricity: +/- 0.148 x 10^12 -> g/2 changes +/- 8.24e-11, versus the 2023 uncertainty 1.3e-13 (634.0 x)
   QED leaves no room: |a_QED - a_meas| ~ 1e-12 -> phi/(2c^2) < ~1e-12 -> phi < 0.0018 x 10^12

3. Lamb shift, hydrogen 2S1/2-2P1/2, eq. (18)
   Aspden's inputs (1983)   eq.(17) 1070.999 MHz;  eq.(18) 1057.8479 MHz
   CODATA 2022 inputs       eq.(17) 1071.0 MHz;  eq.(18) 1057.8491 MHz
   vs Lundeen & Pipkin 1981  1057.845 +/- 0.009 MHz: difference 4.138 kHz = 0.46 sigma
   vs Bezginov et al. 2019   1057.8298 +/- 0.0032 MHz: difference 19.34 kHz = 6.04 sigma
   the two correction factors move eq.(17) by -1.23 %
   step (13) uses dE = 2 hcR dn0 / n0^3 with n0 = 1; Bohr-Sommerfeld gives 2 hcR dn / n^3 with n = 2 -> 8x smaller: 132.23 MHz
   He+ (Z = 2), same formula scaled Z^4: 16932.5 MHz (with 1+2a) or 16692.4 MHz (with 1+2Za); measured 14041.9 +/- 1.5 MHz
