Each result cited is universally quantified over the data in its own statement.
Fix n∈N until Step 5. For maps T:S1→S2 and R:S2→S3 between measurable spaces, each measurable, and a measure λ on S1, one has (R∘T)−1(C)=T−1(R−1(C)) for every measurable C⊆S3, so that the image measures of claim 1 of Image Measures, Measures with Densities, and Change of Variables satisfy
(R∘T)#λ=R#(T#λ).(†)
Let ρ be the product metric on Rn×X, so that ρ((u,w),(u′,w′))=max{∥u−u′∥,∣w−w′∣}, using dE(u,u′)=∥u−u′∥ from Euclidean Space and Lebesgue Measure: Standing Notation §space; its Borel σ-algebra is B(Rn)⊗B(X), as recalled in the statement. On X×X we use the norm, distance and coordinate maps π1,π2 of Borel Probability Measures on a Real Hilbert Space with an Orthonormal Basis: Standing Notation §pairs; by Properties of the Product of Two Real Inner Product Spaces §coordinates, π1 and π2 are linear with ∣πiz∣≤∣z∣, and by Borel Sets of a Hilbert Space with an Orthonormal Basis: Coordinates, Determination by Finite-Dimensional Projections, and Pairs §product-sigma, B(X×X)=B(X)⊗B(X) and π1,π2 are Borel. The function q:X→R, q(y)=∣y∣2, is Borel, as recorded in The Second Moment of a Borel Probability Measure on a Hilbert Space and the Probability Measures with Finite Second Moment.
Step 1 (a coupling). The measures μ,γc∈P(X) have finite mass, so by Borel Sets of a Hilbert Space with an Orthonormal Basis: Coordinates, Determination by Finite-Dimensional Projections, and Pairs §product-measure the product measure λ=μ⊗γc of Existence and Uniqueness of the Product Measure is a Borel measure on X×X with (π1)#λ=γc(X)μ=μ and (π2)#λ=μ(X)γc=γc; it is a probability measure, since λ(X×X)=μ(X)γc(X)=1.
Define jn:X×X→Rn×X by jn(z)=(pn(π1z),Qnπ2z). For z,z′∈X×X, Borel Sets of a Hilbert Space with an Orthonormal Basis: Coordinates, Determination by Finite-Dimensional Projections, and Pairs §continuity (the maps pn and Qn are Lipschitz with constant 1) and the linearity of πi give
ρ(jn(z),jn(z′))≤max{∣π1(z−z′)∣,∣π2(z−z′)∣}≤∣z−z′∣,
so jn is continuous (take δ=ε), hence Borel with respect to B(X×X) and B(Rn)⊗B(X) by claim 3 of Borel Measurability and Bounded Integration on a Metric Space. Let Tn=Φn∘jn:X×X→X, so that
Tn(x,g)=pn∗(pn(x))+Qng=Pnx+Qng(x,g∈X),
using Pn=pn∗∘pn from Borel Probability Measures on a Real Hilbert Space with an Orthonormal Basis: Standing Notation §coordinates. Since Φn is Borel by Head and Tail of a Diagonal Gaussian Measure on a Hilbert Space are Independent §synthesis, Tn is Borel by claim 4 of Borel Measurability and Bounded Integration on a Metric Space. Let Ψn=(π1,Tn):X×X→X×X, that is, Ψn(x,g)=(x,Pnx+Qng); it is Borel by Borel Sets of a Hilbert Space with an Orthonormal Basis: Coordinates, Determination by Finite-Dimensional Projections, and Pairs §pairing. Let π(n)=(Ψn)#λ, a Borel probability measure on X×X by claim 1 of Image Measures, Measures with Densities, and Change of Variables, so π(n)∈P(X×X).
Its first marginal: π1∘Ψn=π1, so by (†), (π1)#π(n)=(π1)#λ=μ.
Its second marginal: π2∘Ψn=Tn=Φn∘jn, so by (†), (π2)#π(n)=(Φn)#((jn)#λ). For A∈B(Rn) and B∈B(X) we have jn−1(A×B)=pn−1(A)×Qn−1(B), a measurable rectangle (pn, Qn being Borel), so
(jn)#λ(A×B)=μ(pn−1(A))γc(Qn−1(B))=((pn)#μ)(A)((Qn)#γc)(B).
The measure (jn)#λ on B(Rn)⊗B(X) therefore has the defining rectangle values of the product of the probability (hence σ-finite) measures (pn)#μ and (Qn)#γc, and by the uniqueness in Existence and Uniqueness of the Product Measure, (jn)#λ=(pn)#μ⊗(Qn)#γc. Hence (π2)#π(n)=μ(n). As μ(n) is a Borel probability measure on X, μ(n)∈P(X), and π(n) is a coupling of μ and μ(n): π(n)∈Π(μ,μ(n)).
Step 2 (the cost of the coupling). Let an=q∘Qn:X→R, an(y)=∣Qny∣2, Borel by claim 4 of Borel Measurability and Bounded Integration on a Metric Space, with 0≤an(y)≤∣y∣2 because ∣Qny∣≤∣y∣ by Borel Sets of a Hilbert Space with an Orthonormal Basis: Coordinates, Determination by Finite-Dimensional Projections, and Pairs §continuity. Let κ(z)=∣π1z−π2z∣2 be the Borel integrand of the quadratic cost. For (x,g)∈X×X,
κ(Ψn(x,g))=∣x−Pnx−Qng∣2=∣Qnx−Qng∣2≤(∣Qnx∣+∣Qng∣)2≤2an(x)+2an(g),
by the triangle inequality and (s+t)2≤2s2+2t2 for real s,t. By claim 2 of Image Measures, Measures with Densities, and Change of Variables (change of variables for nonnegative functions), the monotonicity and additivity of the integral of nonnegative measurable functions (Linearity and Monotonicity of the Lebesgue Integral, in force by Measure Spaces and the Lebesgue Integral: Standing Notation §background), and claim 2 of Image Measures, Measures with Densities, and Change of Variables once more together with the marginals of Step 1,
I(π(n))=∫X×Xκ∘Ψndλ≤2∫X×Xan∘π1dλ+2∫X×Xan∘π2dλ=2∫Xandμ+2∫Xandγc.
Write εn for the right-hand side. Since an≤q, monotonicity gives εn≤2M2(μ)+2M2(γc), with the second moments of The Second Moment of a Borel Probability Measure on a Hilbert Space and the Probability Measures with Finite Second Moment §moment; here M2(μ)<∞ because μ∈P2(X), and M2(γc)<∞ because γc∈P2(X) by Moments of a Diagonal Gaussian Measure on a Hilbert Space: Coordinate Covariances, Finite Second Moment, and Exponential Moments of the Squared Norm §moment. Thus I(π(n))≤εn<∞.
Step 3 (claim 1). Since μ∈P2(X), π(n)∈Π(μ,μ(n)) and I(π(n))<∞, the converse part of Couplings on a Hilbert Space: Product Coupling, Swap, Finiteness of the Cost, Push-Forward Couplings, Modifying One Marginal, Quantisation, Gluing over a Finitely Supported Measure, the Lipschitz Bound and the Moment Bound §cost-finite gives μ(n)∈P2(X).
For the head, note that pn is linear (each coordinate x↦⟨x,ek⟩ is linear) and that pn(pn∗(y))=y for y∈Rn by Borel Sets of a Hilbert Space with an Orthonormal Basis: Coordinates, Determination by Finite-Dimensional Projections, and Pairs §continuity; hence pn(Pny)=pn(y) and pn(Qny)=pn(y)−pn(Pny)=0 for y∈X. Consequently pn(Tn(x,g))=pn(x)+0=pn(x), that is, pn∘π2∘Ψn=pn∘π1 on X×X. By (†), applied twice to each three-fold composite, and Step 1,
(pn)#μ(n)=(pn∘π2∘Ψn)#λ=(pn∘π1)#λ=(pn)#((π1)#λ)=(pn)#μ.
This proves claim 1 for the arbitrary n fixed above.
Step 4 (a bound on W2). Both μ(n) and μ belong to P2(X), so W2(μ(n),μ) is defined. With σ(x,y)=(y,x) the swap of X×X, Couplings on a Hilbert Space: Product Coupling, Swap, Finiteness of the Cost, Push-Forward Couplings, Modifying One Marginal, Quantisation, Gluing over a Finitely Supported Measure, the Lipschitz Bound and the Moment Bound §swap gives σ#π(n)∈Π(μ(n),μ) and I(σ#π(n))=I(π(n)). By The Quadratic Wasserstein Distance on a Hilbert Space §distance and Step 2,
W2(μ(n),μ)2≤I(σ#π(n))=I(π(n))≤εn.
Step 5 (claim 2). Let η∈{μ,γc}. The function q is integrable with respect to η, being nonnegative with ∫Xqdη=M2(η)<∞ (Measure Spaces and the Lebesgue Integral: Standing Notation §integral). For every y∈X, by Borel Probability Measures on a Real Hilbert Space with an Orthonormal Basis: Standing Notation §coordinates the sequence (Xm)m∈N is exhausting, Pm is the orthogonal projection onto Xm and Qmy=y−Pmy, so Exhausting Sequences of Finite-Dimensional Subspaces in a Separable Real Hilbert Space, and Their Projections §tail gives that (Qmy)m∈N converges to 0X, that is, ∣Qmy∣→0, and hence am(y)=∣Qmy∣2→0. Since ∣am(y)∣≤q(y) for all m and y, claim 3 of Dominated Convergence Theorem gives ∫Xamdη→0 as m→∞. Hence εm→0.
Now let ε>0 be real and choose m0∈N with εm<ε2 for every m≥m0. By Step 4, applied with each m≥m0 in place of n, W2(μ(m),μ)2<ε2, and as W2(μ(m),μ)≥0 this gives W2(μ(m),μ)<ε. Thus limmW2(μ(m),μ)=0. Since (μ(m))m∈N is a sequence in P2(X) by claim 1 and μ∈P2(X), Wasserstein Convergence on a Hilbert Space: Weak Convergence, Integrals of Continuous Functions of Quadratic Growth, Convergence from Weak Convergence with Uniformly Integrable Second Moments, and Compactness §weak gives μ(m)⇒μ. This proves claim 2.