Each result cited is universally quantified over the data in its own statement.
Claim 1. The function u is continuous by condition (d) of Test Functions on the Wasserstein Space: the Intrinsic Gradient and the Translation Hessian (i), and E is lower semicontinuous on D relative to D by (E1). By Basic Properties of the Delta-Envelopes on the Wasserstein Space (iii), applied for each positive δ∈R, the function u is bounded above near each point and bounded below near each point of P2(Rd), and uδ−=u−δE and uδ+=u+δE on D.
Claim 2. Let δ, φ and μ^ be as stated, and put χ=u−φ. By Sums, Real Multiples and Differences of Test Functions on the Wasserstein Space (ii) the function χ is a test function on P2(Rd) with
∇χ(ν)=∇u(ν)−∇φ(ν),Hχ(ν)=Hu(ν)−Hφ(ν)(ν∈P2(Rd)).
By claim 1, for μ∈D one has uδ−(μ)−φ(μ)=u(μ)−δE(μ)−φ(μ)=χ(μ)−δE(μ), so μ^ is a point of D at which χ−δE has a local maximum relative to D. By (E2), μ^∈DΣ, and by First- and Second-Order Conditions at a Penalised Extremum of a Test Function on the Wasserstein Space (i),
∇u(μ^)−∇φ(μ^)=δΣ(μ^),Hu(μ^)−Hφ(μ^)⪯0d.
The first identity gives ∇φ(μ^)+δΣ(μ^)=∇u(μ^) by Elementary Identities in a Vector Space, and the second gives Hu(μ^)⪯Hφ(μ^) by claim 1 of Comparison with the Zero Matrix in the Positive Semidefinite Ordering. By claim 1 again, uδ−(μ^)+δE(μ^)=u(μ^). Hence, by The Bundle of Vector Fields over a Set of Measures, Second-Order Equation Operators on the Wasserstein Space, and Their Delta-Shifts (iii),
Fδ−(μ^,uδ−(μ^),∇φ(μ^),Hφ(μ^))=F(μ^,u(μ^),∇u(μ^),Hφ(μ^))≤F(μ^,u(μ^),∇u(μ^),Hu(μ^))≤0,
the first inequality by Degenerate Elliptic Second-Order Equation Operators on the Wasserstein Space (i) applied to Hu(μ^)⪯Hφ(μ^), and the second because u is a classical subsolution of F on DΣ and μ^∈DΣ (Classical Sub- and Supersolutions of a Second-Order Equation on the Wasserstein Space (i)).
Claim 3. By claim 1 the function u is bounded above near each point of P2(Rd), so the envelopes uδ− are defined. Let δ∈R be positive, let φ be a test function on P2(Rd), let μ^∈D be a point at which the function with value uδ−(μ)−φ(μ) at μ has a local maximum relative to D, and let ε∈R be positive. By claim 2, μ^∈DΣ. Take
ν=μ^,π=(id,id)#μ^,s=uδ−(μ^),q=∇φ(μ^),Y=Hφ(μ^),
where id is the identity map of Rd. By Couplings on Euclidean Space: Product Coupling, Swap, Finiteness of the Cost, Push-Forward Couplings, Modifying One Marginal, Quantisation, Gluing over a Finitely Supported Measure, and the Lipschitz Bound (iv), applied with S=T=id, one has π∈Π(μ^,μ^) and
I(π)=∫Rd∥id−id∥2dμ^=0.
The number ε2 is positive by claim 5 of Elementary Order Arithmetic in an Ordered Field, so I(π)<ε2. Since ν=μ^ and s=uδ−(μ^), the differences uδ−(ν)−uδ−(μ^) and s−uδ−(μ^) are 0, of absolute value 0<ε by claim 1 of Properties of the Absolute Value in an Ordered Field. By The Discrepancy of Two Square-Integrable Vector Fields Along a Coupling of Their Base Measures (ii) the discrepancy of q and ∇φ(μ^) along π is ∥q−∇φ(μ^)∥μ^2=0<ε2, the norm of the zero element vanishing by Elementary Identities in a Real Inner Product Space (iii). Also Y−Hφ(μ^)=0d by Difference of Real Matrices, so ∥Y−Hφ(μ^)∥=0<ε by claim 4 of Properties of the Norm of a Symmetric Real Matrix. Finally, claim 2 gives Fδ−(ν,s,q,Y)≤0≤ε. Every requirement of Viscosity Subsolution, Supersolution and Solution of a Second-Order Equation on the Wasserstein Space (i) is therefore met, and u is a viscosity subsolution of F relative to the penalty pair.
Claim 4. Let δ, φ and μ^ be as stated and put χ=u−φ as in claim 2. By claim 1, uδ+(μ)−φ(μ)=χ(μ)+δE(μ) for μ∈D, so μ^ is a point of D at which χ+δE has a local minimum relative to D. By (E2), μ^∈DΣ, and by First- and Second-Order Conditions at a Penalised Extremum of a Test Function on the Wasserstein Space (ii),
∇u(μ^)−∇φ(μ^)=−δΣ(μ^),0d⪯Hu(μ^)−Hφ(μ^).
The first identity gives ∇φ(μ^)−δΣ(μ^)=∇u(μ^) by Elementary Identities in a Vector Space, and the second gives Hφ(μ^)⪯Hu(μ^) by claim 2 of Comparison with the Zero Matrix in the Positive Semidefinite Ordering. By claim 1, uδ+(μ^)−δE(μ^)=u(μ^). Hence, by The Bundle of Vector Fields over a Set of Measures, Second-Order Equation Operators on the Wasserstein Space, and Their Delta-Shifts (iii),
Fδ+(μ^,uδ+(μ^),∇φ(μ^),Hφ(μ^))=F(μ^,u(μ^),∇u(μ^),Hφ(μ^))≥F(μ^,u(μ^),∇u(μ^),Hu(μ^))≥0,
the first inequality by Degenerate Elliptic Second-Order Equation Operators on the Wasserstein Space (i) applied to Hφ(μ^)⪯Hu(μ^), and the second because u is a classical supersolution of F on DΣ and μ^∈DΣ (Classical Sub- and Supersolutions of a Second-Order Equation on the Wasserstein Space (ii)).
Claim 5. With the data of claim 3, taken now at a local minimum of the function with value uδ+(μ)−φ(μ) at μ and with s=uδ+(μ^), the same computations give I(π)<ε2, ∣uδ+(ν)−uδ+(μ^)∣=0<ε, ∣s−uδ+(μ^)∣=0<ε, discrepancy 0<ε2 and ∥Y−Hφ(μ^)∥=0<ε, while claim 4 gives −ε≤0≤Fδ+(ν,s,q,Y), the first inequality by claim 4 of Elementary Order Arithmetic in an Ordered Field. Every requirement of Viscosity Subsolution, Supersolution and Solution of a Second-Order Equation on the Wasserstein Space (ii) is met, so u is a viscosity supersolution of F relative to the penalty pair.
Claim 6. A classical solution of F on DΣ is both a classical subsolution and a classical supersolution there, by Classical Sub- and Supersolutions of a Second-Order Equation on the Wasserstein Space (iii); by claims 3 and 5 it is then both a viscosity subsolution and a viscosity supersolution of F relative to the penalty pair, hence a viscosity solution by Viscosity Subsolution, Supersolution and Solution of a Second-Order Equation on the Wasserstein Space (iii), the required local bounds holding by claim 1.