Fix Σ∈Δl and α∈Rm, and for each ordered pair (σ,γ)∈{1,…,l}2 with σ=γ abbreviate wσγ=Σσβ(σ,γ,Σ,α). By the definition of the probability simplex, Σσ≥0, and by clause 1 (bounds) of the definition of a transition-rate family, β(σ,γ,Σ,α)≥0; hence wσγ≥0.
Step 1 (Entries of the outer products). Fix an ordered pair (σ,γ) with σ=γ and set u=eγ−eσ, whose component up equals 1 for p=γ, equals −1 for p=σ, and equals 0 otherwise. By the definitions of the matrix product and the transpose, uu⊤ is the matrix with l rows and l columns whose entry in row p and column q is upuq: the row-p, column-q entry of the product of the one-column matrix u and the one-row matrix u⊤ is the product of the sole entry up of row p of u and the sole entry uq of column q of u⊤.
Step 2 (Jump representation). Write Θ′ for the right-hand side of part 1, so that, sums of matrices being entrywise,
(Θ′)pq=(σ,γ):σ=γ∑wσγupuq(p,q∈{1,…,l}),
with u=eγ−eσ depending on the pair as in Step 1.
First take q=p. The term of the pair (σ,γ) contributes wσγ(up)2, which equals wσγ when p∈{σ,γ} and 0 otherwise. The pairs with γ=p contribute ∑σ:σ=pwσp and the pairs with σ=p contribute ∑γ:γ=pwpγ, so, renaming the summation index to σ in both sums,
(Θ′)pp=σ:σ=p∑(Σσβ(σ,p,Σ,α)+Σpβ(p,σ,Σ,α)),
which is exactly the diagonal entry Θpp(Σ,α) prescribed by the definition of the aggregate fluctuation covariance.
Now fix p=q. The term of the pair (σ,γ) contributes wσγupuq, which vanishes unless both p and q lie in {σ,γ}; since p=q and σ=γ, this happens exactly for (σ,γ)=(p,q) and (σ,γ)=(q,p). In the first case upuq=(−1)⋅1=−1, and in the second upuq=1⋅(−1)=−1. Hence
(Θ′)pq=−wpq−wqp=−Σpβ(p,q,Σ,α)−Σqβ(q,p,Σ,α)=Θpq(Σ,α),
again by the definition of the aggregate fluctuation covariance. Thus Θ′=Θ(Σ,α), proving part 1.
Step 3 (Quadratic form). Let x∈Rl. By part 1, entrywise, and interchange of the finite sums,
p=1∑lq=1∑lΘpq(Σ,α)xpxq=(σ,γ):σ=γ∑wσγp=1∑lq=1∑lupuqxpxq=(σ,γ):σ=γ∑wσγ(p=1∑lupxp)2,
and ∑p=1lupxp=xγ−xσ by the description of u in Step 1. This proves part 2.
Step 4 (Positive semidefiniteness). Θ(Σ,α) is symmetric by its definition. Let x∈Rl. By the definitions of the matrix-vector product and the dot product,
x⋅(Θ(Σ,α)x)=p=1∑lxpq=1∑lΘpq(Σ,α)xq=p=1∑lq=1∑lΘpq(Σ,α)xpxq,
which by part 2 and wσγ≥0 is a finite sum of nonnegative terms, hence nonnegative. By the definition of a positive semidefinite matrix, Θ(Σ,α) is positive semidefinite, proving part 3.