This application keeps both the electron relative position and a quantum source coordinate. Its first calculation concerns a fixed Coulomb reference, whose exact spectral measure can be evaluated. The second restores the actual source interaction and feedback. These are different estimates: a continuum-only coefficient of the fixed reference cannot remove the poles or replace the coupled resolvent of the actual parent.
The Coulomb operator is self-adjoint on H2(R3) with form domain H1, ground energy −1/2 and normalized ground ϕ. Its excited bound energies are En=−1/(2n2), n≥2, with continuum [0,∞). We use the usual Coulomb spectrum, not a newly inferred set of levels. Only the isolated forcing f has the exact angular selection l=1,m=0; the actual coupled parent below retains every angular sector.
Theorem 9.1 (Full spectral measure of the dipole forcing)
Its total mass and first (h+1) form moment are both one:
∫dμf=1,∫Edμf=0,∫(E+1)dμf=1.(9.5)
On the whole true continuum,
ρC(E)≤ρ∗,(1+E)ρC(E)≤ρ∗,ρ∗=3e4256.(9.6)
Proof
With ψ=u(r)Y10(r)/r, the radial forcing and Friedrichs operator are
uf=32r2e−r,h1=−21∂r2+r−2−r−1.
Integration of r4e−2r gives ∥f∥2=1. Direct differentiation gives h1uf=(−1/2+1/r)uf, whence ⟨f,hf⟩=0.
We specify the continuum normalization because an energy cross-section convention could otherwise introduce a wrong Jacobian. The regular Coulomb function is normalized by its unit-amplitude large-radius oscillation. In the standard convention [10],
Its derivative jump is −2, as required by the coefficient −1/2 in h1. Its imaginary part is (2/k)F1(r)F1(r′). Stone's formula therefore gives the measure ukukdk, or ukuk/k per unit energy E=k2/2. The Coulomb asymptotics used to choose H1+ are asymptotic statements at large radius, not exact finite-radius equalities.
Here is the overlap calculation. Put a=2+i/k, s=1+ik, z=2ik/(1+ik). Laplace integration of the confluent hypergeometric series, first in its absolute-convergence region and then by analytic continuation, gives
∫0∞r4e−sr1F1(a;4;2ikr)dr=24s−52F1(a,5;4;z).
The coefficient identity (5)j/(4)j=1+j/4 implies
2F1(a,5;4;z)=(1−z)−a(1+4(1−z)az).
The last bracket is 1/[2(1−ik)]. Consequently the integral is
(1+ik)5(1−ik)12(1+ik1−ik)−2−i/k.
The continuous logarithm of the ratio is −2iarctank. Its contribution to the squared modulus is e−4arctan(k)/k. Multiplication by the normalized radial prefactors gives ∣⟨uk,uf⟩∣2=kρC(k2/2). The factor 1/k from Stone's energy measure proves (9.4).
The normalized negative-energy eigenfunction in this channel is
un1(r)=n34r2(n+1)!(n−2)!e−r/nLn−23(2r/n).
Insert Lj3=(4)j1F1(−j;4;⋅)/j! into the same Laplace identity. The bracket is now n/[2(n−1)] and gives (9.3); in particular d22=32768/59049.
For completeness, the radial regular/decaying solutions have only the simple negative poles −1/(2n2), n≥2. Their residues give precisely these normalized eigenprojections. On every compact positive-energy interval, (9.7) has continuous boundary values against the exponentially decaying forcing, so its forcing measure there is absolutely continuous. A remaining singular measure could only be supported at zero. The regular zero-energy solution is proportional to rJ3(8r); its large-radius oscillatory magnitude is of order r1/4 and it is not L2. Thus zero is not an eigenvalue and supplies no atom. A singular continuous measure cannot be supported on a single point. This proves completeness of (9.2). The already evaluated norm and first energy expectation now give (9.5); numerical quadrature is not used to establish them.
To prove the whole-continuum ceiling, let A(k)=arctank−k/(1+k2). Since
dkd{3(1+k2)2k3−A(k)}=3(1+k2)22k4≥0,
we have A(k)≤2k3/[3(1+k2)]. The logarithmic derivative of (1+k2/2)ρC(k2/2) is consequently at most
−3(1+k2)16k+k2(e2π/k−1)2π.
For k≤1, use ex−1≥x3/6 to bound the positive term by 3k/(2π2); this leaves a strictly negative margin. For k≥1, use ex−1>x to bound it by 1/k, whereas the negative magnitude is at least 8/(3k). The threshold limit is ρ∗. This proves the weighted ceiling and hence the unweighted one.
□
The bound lines are substantial. The formula gives n3dn2⟶ρ∗ and En+1−En∼n−3, matching the threshold density. For n≥100 one also has
dn2≤4n37,n>N∑dn2≤8N27.
For example the logarithm of the factor multiplying 256/(3n3) is at most −4+10/n≤−3.9, and 256e−3.9/3<7/4 by a positive exponential Taylor sum. These tails can check a truncated calculation; they do not turn the discrete spectrum into a continuum.
Extend the time input by zero outside [0,t]. Hilbert-valued Plancherel, integrated first on the whole frequency line, gives the factor 2π for this unnormalized Fourier integral. Restricting to E>0 and using (9.6) gives both estimates in (9.9). The source form commutes its own propagation, so the same argument applied to Hs1/2a gives (9.10).
□
For a=gYp, the last input is gHs1/2Yp, in that order. It is not gYHs1/2p. The source interaction picture here is used only for the norm proof; it has not changed the physical source coordinate or its current. The bound poles are excluded from this proposition and must be retained separately: a Poisson-smoothed atom contributes dn2/(πη) at its center. There is no uniform η-independent full-spectrum density ceiling.