Preliminaries. Write EK=∑k=0Kμk/k! for K∈N0 and E−1=E−2=0. By The Real Exponential Function, EK→exp(μ) as K→∞ (the partial sums of the defining series converge to exp(μ)), and hence also EK−1→exp(μ) and EK−2→exp(μ) as K→∞: for ε>0, if ∣EL−exp(μ)∣<ε for all L≥L0, then ∣EK−1−exp(μ)∣<ε and ∣EK−2−exp(μ)∣<ε for all K≥L0+2, which is the definition of the limit. By claim 2 of Basic Properties of the Exponential Function, exp(−μ)exp(μ)=1 and exp(−μ)>0. Since every term μk/k! is positive, EK is nondecreasing in K and μk/k!≤EK for k≤K; moreover EK≤exp(μ), by claim 1 of Order Properties of Limits of Real Sequences applied to the constant sequence with value EK (which converges to EK) and the sequence (EL)L≥K, which dominates it termwise and converges to exp(μ) (a tail of a convergent sequence converges to the same limit, by the argument just given for EK−1). Hence 0<pμ(k)≤exp(−μ)exp(μ)=1 for k∈N0, so pμ takes values in [0,1] and {pμ>0}=N0. Finally, each singleton {k}=⋂j∈N(k−1/j,k+1/j) is a Borel set (a countable intersection of open intervals), and by Poisson Distribution, Pμ({k})=exp(−μ)∑j∈N0∩{k}μj/j!=exp(−μ)μk/k!=pμ(k) for k∈N0, as asserted in the statement.
(P) Let f:R→[0,∞) vanish on R∖N0, suppose the partial sums sK=∑k=0Kf(k) converge to a real number s, and let A be a set with N0⊆A⊆R. Then the sum of f over A equals s (note s≥0: each sK≥0, so claim 1 of Order Properties of Limits of Real Sequences with the constant sequence 0 gives s≥0). Indeed, for a finite set F⊆A, either F∩N0=∅ and ∑x∈Ff(x)=0≤s (if F=∅ every term vanishes, so the sum is a sum of zeros over [∣F∣] along an enumeration by Sum over a Finite Index Set, hence 0 by Finite Sum Notation in a Field; if F=∅ the empty-sum convention applies), or K=max(F∩N0) exists and ∑x∈Ff(x)=∑x∈F∩N0f(x)≤sK≤s, using claims 4 and 3 of Peeling, Splitting, and Interchange for Sums over a Finite Index Set (terms vanishing off F∩N0; F∩N0⊆{0,…,K} and the omitted terms are nonnegative, the case F∩N0={0,…,K} being an equality) and then claim 1 of Order Properties of Limits of Real Sequences applied to the constant sequence sK and the nondecreasing sequence (sL)L≥K, which converges to s as a tail of a convergent sequence (by the argument given for EK−1 in the preliminaries). So s is an upper bound of the finite sums. Conversely each sK is itself a finite sum (over F={0,…,K}⊆A), so any upper bound b satisfies b≥sK for all K, whence b≥s by the same claim. Thus s is the least upper bound.
Claim 1. Apply (P) to f=pμ and A=R: sK=exp(−μ)EK→exp(−μ)exp(μ)=1 by claim 2 of Arithmetic of Limits of Real Sequences. So ∑x∈Rpμ(x)=1, and pμ is a discrete probability mass function by Discrete Probability Mass Function on Euclidean Space; the identities pμ(k)=Pμ({k}) and {pμ>0}=N0 were shown in the preliminaries.
Claim 2. Let k∈N0. If k≥1, then k!=k⋅(k−1)! by Factorial of a Natural Number (with 0!=1 when k=1) and μk=μ⋅μk−1, so
pμ(k)pμ(k−1)=μk/k!μk−1/(k−1)!=μk.
If k=0, then k−1=−1∈/N0, so pμ(k−1)=0=k/μ. In both cases ρpμ,a,w(k)=w1(1−pμ(k−a1)/pμ(k))=w1(1−k/μ) by Move Score of a Discrete Probability Mass Function. Finally k+a1=k+1∈N⊆N0 for every k∈N0 by Natural Numbers, so no k∈N0 has k+a1∈/N0.
Claim 3. By Move Information of a Discrete Probability Mass Function and claim 2, J(pμ;a,w) is the sum over N0 of f(k)=pμ(k)w12(1−k/μ)2=μ2w12pμ(k)(k−μ)2, and we extend f by 0 to R∖N0; by (P) with A=N0 it suffices to show that the partial sums ∑k=0Kpμ(k)(k−μ)2 converge to μ. Using (k−μ)2=k(k−1)+k−2μk+μ2 and, for k≥1, kμk/k!=μμk−1/(k−1)!, and for k≥2, k(k−1)μk/k!=μ2μk−2/(k−2)! (the terms with k=0, respectively k≤1, vanish),
k=0∑Kkk!μk=μEK−1,k=0∑Kk(k−1)k!μk=μ2EK−2,
so that
k=0∑Kpμ(k)(k−μ)2=exp(−μ)(μ2EK−2+μEK−1−2μ2EK−1+μ2EK).
By the preliminaries and claims 1 and 3 of Arithmetic of Limits of Real Sequences, the right side converges to exp(−μ)exp(μ)(μ2+μ−2μ2+μ2)=μ; multiplying by the constant w12/μ2 (claim 3 there), the partial sums of f converge to w12/μ. Hence J(pμ;a,w)=μ2w12⋅μ=μw12. ■