Each result cited below is universally quantified over the data in its own statement.
For μ ∈ P ( T d ) \mu\in\mathcal{P}(\mathbb{T}^{d}) μ ∈ P ( T d ) the class of a square-integrable Borel map G : R d → R d G:\mathbb{R}^{d}\to\mathbb{R}^{d} G : R d → R d in L 2 ( μ ; R d ) L^{2}(\mu;\mathbb{R}^{d}) L 2 ( μ ; R d ) is written [ G ] [G] [ G ] ; by The Space of Square-Integrable Random Vectors §classes (on the probability space ( R d , B ( R d ) , μ ) (\mathbb{R}^{d},\mathcal{B}(\mathbb{R}^{d}),\mu) ( R d , B ( R d ) , μ ) , see Optimal Transport on the Flat Torus: Standing Notation §fields ), [ G ] = [ H ] [G]=[H] [ G ] = [ H ] exactly when μ ( { y : G ( y ) = H ( y ) } ) = 1 \mu(\{y:G(y)=H(y)\})=1 μ ({ y : G ( y ) = H ( y )}) = 1 , and a [ G ] + b [ H ] = [ a G + b H ] a[G]+b[H]=[aG+bH] a [ G ] + b [ H ] = [ a G + b H ] for real a , b a,b a , b . A sequence ( μ n ) (\mu_{n}) ( μ n ) converges to μ \mu μ in ( P ( T d ) , W T ) (\mathcal{P}(\mathbb{T}^{d}),W_{\mathbb{T}}) ( P ( T d ) , W T ) exactly when the real sequence ( W T ( μ n , μ ) ) (W_{\mathbb{T}}(\mu_{n},\mu)) ( W T ( μ n , μ )) converges to 0 0 0 , by Convergent Sequence in a Metric Space .
Step 0 (a Lipschitz bound in the torus distance). Let u ∈ C p e r 1 u\in C^{1}_{\mathrm{per}} u ∈ C per 1 . We show that there is a real L u ≥ 0 L_{u}\ge0 L u ≥ 0 with ∣ u ( y ) − u ( x ) ∣ ≤ L u d T ( x , y ) |u(y)-u(x)|\le L_{u}\,d_{\mathbb{T}}(x,y) ∣ u ( y ) − u ( x ) ∣ ≤ L u d T ( x , y ) for all x , y ∈ R d x,y\in\mathbb{R}^{d} x , y ∈ R d . For each i ∈ [ d ] i\in[d] i ∈ [ d ] the function ∂ i u \partial_{i}u ∂ i u belongs to C p e r C_{\mathrm{per}} C per by Elementary Properties of Lattice-Periodic Functions §derivative , so there is a real M i ≥ 0 M_{i}\ge0 M i ≥ 0 with ∣ ∂ i u ∣ ≤ M i |\partial_{i}u|\le M_{i} ∣ ∂ i u ∣ ≤ M i on R d \mathbb{R}^{d} R d by Elementary Properties of Lattice-Periodic Functions §bounded ; put L u = M 1 + ⋯ + M d L_{u}=M_{1}+\dots+M_{d} L u = M 1 + ⋯ + M d . Let x , y ∈ R d x,y\in\mathbb{R}^{d} x , y ∈ R d , put v = ϖ ( y − x ) v=\varpi(y-x) v = ϖ ( y − x ) and z = x + v z=x+v z = x + v . By The Flat Torus Distance: Minimality of the Wrapped Displacement, Periodicity, the Metric on the Unit Cell and the Lipschitz Bound §range , y − z = ( y − x ) − ϖ ( y − x ) ∈ Z d y-z=(y-x)-\varpi(y-x)\in\mathbb{Z}^{d} y − z = ( y − x ) − ϖ ( y − x ) ∈ Z d , so u ( y ) = u ( z ) u(y)=u(z) u ( y ) = u ( z ) by the periodicity of u u u ; and ∥ v ∥ = d T ( x , y ) \lVert v\rVert=d_{\mathbb{T}}(x,y) ∥ v ∥ = d T ( x , y ) by The Wrapped Displacement and the Flat Torus Distance §distance . If v = 0 v=0 v = 0 , then u ( y ) = u ( x ) u(y)=u(x) u ( y ) = u ( x ) . Otherwise let F : [ 0 , 1 ] → R F:[0,1]\to\mathbb{R} F : [ 0 , 1 ] → R , F ( τ ) = u ( x + τ v ) F(\tau)=u(x+\tau v) F ( τ ) = u ( x + τv ) . Since u u u is of class C 1 C^{1} C 1 on the open set R d \mathbb{R}^{d} R d (claim 1 of Euclidean Space is Open in Itself, and C k C^k C k Maps are Continuous ), Derivatives Along a Segment for C^1 Functions on a Euclidean Open Set (claims 2 and 3, with J = [ 0 , 1 ] J=[0,1] J = [ 0 , 1 ] ) shows that F F F is continuous on [ 0 , 1 ] [0,1] [ 0 , 1 ] and differentiable at every τ ∈ ( 0 , 1 ) \tau\in(0,1) τ ∈ ( 0 , 1 ) with F ′ ( τ ) = ∑ i = 1 d ∂ i u ( x + τ v ) v i F'(\tau)=\sum_{i=1}^{d}\partial_{i}u(x+\tau v)\,v_{i} F ′ ( τ ) = ∑ i = 1 d ∂ i u ( x + τv ) v i . By the mean value theorem there is c ∈ ( 0 , 1 ) c\in(0,1) c ∈ ( 0 , 1 ) with u ( z ) − u ( x ) = F ′ ( c ) u(z)-u(x)=F'(c) u ( z ) − u ( x ) = F ′ ( c ) , whence ∣ u ( y ) − u ( x ) ∣ = ∣ F ′ ( c ) ∣ ≤ ∑ i M i ∣ v i ∣ ≤ L u ∥ v ∥ = L u d T ( x , y ) |u(y)-u(x)|=|F'(c)|\le\sum_{i}M_{i}|v_{i}|\le L_{u}\lVert v\rVert=L_{u}\,d_{\mathbb{T}}(x,y) ∣ u ( y ) − u ( x ) ∣ = ∣ F ′ ( c ) ∣ ≤ ∑ i M i ∣ v i ∣ ≤ L u ∥ v ∥ = L u d T ( x , y ) .
Consequently, if G : R d → R d G:\mathbb{R}^{d}\to\mathbb{R}^{d} G : R d → R d has components G 1 , … , G d ∈ C p e r 1 G_{1},\dots,G_{d}\in C^{1}_{\mathrm{per}} G 1 , … , G d ∈ C per 1 , then with L G = L G 1 + ⋯ + L G d L_{G}=L_{G_{1}}+\dots+L_{G_{d}} L G = L G 1 + ⋯ + L G d ,
∥ G ( y ) − G ( x ) ∥ ≤ ∑ i = 1 d ∣ G i ( y ) − G i ( x ) ∣ ≤ L G d T ( x , y ) ( x , y ∈ R d ) . \lVert G(y)-G(x)\rVert\le\sum_{i=1}^{d}|G_{i}(y)-G_{i}(x)|\le L_{G}\,d_{\mathbb{T}}(x,y)\qquad(x,y\in\mathbb{R}^{d}). ∥ G ( y ) − G ( x )∥ ≤ i = 1 ∑ d ∣ G i ( y ) − G i ( x ) ∣ ≤ L G d T ( x , y ) ( x , y ∈ R d ) .
Claim 1. Let g g g and h h h be gradient fields of φ \varphi φ .
Step 1 (agreement at heat-smoothed measures). Let μ ∈ P ( T d ) \mu\in\mathcal{P}(\mathbb{T}^{d}) μ ∈ P ( T d ) , let 0 < s ≤ 1 2 0<s\le\tfrac12 0 < s ≤ 2 1 , let S s S_{s} S s be the heat semigroup at time s s s , and put ν = S s μ \nu=S_{s}\mu ν = S s μ ; then ν ∈ P I ( T d ) ⊆ P ( T d ) \nu\in\mathcal{P}^{\mathcal{I}}(\mathbb{T}^{d})\subseteq\mathcal{P}(\mathbb{T}^{d}) ν ∈ P I ( T d ) ⊆ P ( T d ) by Heat Smoothing on the Torus: Distance to the Identity, Contraction, a Smooth Positive Periodic Density, Finite Entropy, and Finite Fisher Information §score . By Gradient Fields and Laplacian Test Functions on the Torus Wasserstein Space §represents , applied to g g g and to h h h , [ g ( ν , ⋅ ) ] = ∇ φ ( ν ) = [ h ( ν , ⋅ ) ] [g(\nu,\cdot)]=\nabla\varphi(\nu)=[h(\nu,\cdot)] [ g ( ν , ⋅ )] = ∇ φ ( ν ) = [ h ( ν , ⋅ )] , so the set D = { y ∈ R d : g ( ν , y ) ≠ h ( ν , y ) } D=\{y\in\mathbb{R}^{d}:g(\nu,y)\ne h(\nu,y)\} D = { y ∈ R d : g ( ν , y ) = h ( ν , y )} satisfies ν ( D ) = 0 \nu(D)=0 ν ( D ) = 0 . For each i i i the functions g i ( ν , ⋅ ) g_{i}(\nu,\cdot) g i ( ν , ⋅ ) and h i ( ν , ⋅ ) h_{i}(\nu,\cdot) h i ( ν , ⋅ ) lie in C p e r 1 C^{1}_{\mathrm{per}} C per 1 by Gradient Fields and Laplacian Test Functions on the Torus Wasserstein Space §regular , hence are continuous (C^k Maps on a Euclidean Open Set , clause 1) and periodic; so D D D , the union over i ∈ [ d ] i\in[d] i ∈ [ d ] of the sets where g i ( ν , ⋅ ) − h i ( ν , ⋅ ) ≠ 0 g_{i}(\nu,\cdot)-h_{i}(\nu,\cdot)\ne0 g i ( ν , ⋅ ) − h i ( ν , ⋅ ) = 0 , is open (hence Borel), and D + k = D D+k=D D + k = D for k ∈ Z d k\in\mathbb{Z}^{d} k ∈ Z d .
Suppose that D D D has a point y 0 y_{0} y 0 . By The Half-Open Unit Cell Tiles Euclidean Space §wrap , z 0 = π ( y 0 ) ∈ Q z_{0}=\pi(y_{0})\in Q z 0 = π ( y 0 ) ∈ Q and y 0 − z 0 ∈ Z d y_{0}-z_{0}\in\mathbb{Z}^{d} y 0 − z 0 ∈ Z d , so z 0 ∈ D z_{0}\in D z 0 ∈ D , and there is r > 0 r>0 r > 0 with the open ball B ( z 0 , r ) ⊆ D B(z_{0},r)\subseteq D B ( z 0 , r ) ⊆ D . Choose a real δ \delta δ with 0 < δ < 1 − ( z 0 ) i 0<\delta<1-(z_{0})_{i} 0 < δ < 1 − ( z 0 ) i for every i ∈ [ d ] i\in[d] i ∈ [ d ] and δ d < r \delta\sqrt{d}<r δ d < r , let c ∈ R d c\in\mathbb{R}^{d} c ∈ R d have coordinates c i = ( z 0 ) i + δ / 2 c_{i}=(z_{0})_{i}+\delta/2 c i = ( z 0 ) i + δ /2 , and let B ′ = B ( c , δ / 2 ) B'=B(c,\delta/2) B ′ = B ( c , δ /2 ) . For x ∈ B ′ x\in B' x ∈ B ′ and i ∈ [ d ] i\in[d] i ∈ [ d ] one has ∣ x i − c i ∣ ≤ ∥ x − c ∥ < δ / 2 |x_{i}-c_{i}|\le\lVert x-c\rVert<\delta/2 ∣ x i − c i ∣ ≤ ∥ x − c ∥ < δ /2 , so 0 ≤ ( z 0 ) i < x i < ( z 0 ) i + δ < 1 0\le(z_{0})_{i}<x_{i}<(z_{0})_{i}+\delta<1 0 ≤ ( z 0 ) i < x i < ( z 0 ) i + δ < 1 , whence x ∈ Q x\in Q x ∈ Q ; and ∥ x − z 0 ∥ ≤ ∥ x − c ∥ + ∥ c − z 0 ∥ < δ / 2 + δ d / 2 ≤ δ d < r \lVert x-z_{0}\rVert\le\lVert x-c\rVert+\lVert c-z_{0}\rVert<\delta/2+\delta\sqrt{d}/2\le\delta\sqrt{d}<r ∥ x − z 0 ∥ ≤ ∥ x − c ∥ + ∥ c − z 0 ∥ < δ /2 + δ d /2 ≤ δ d < r , whence x ∈ D x\in D x ∈ D . Thus B ′ ⊆ D ∩ Q B'\subseteq D\cap Q B ′ ⊆ D ∩ Q , and λ d ( B ′ ) > 0 \lambda_{d}(B')>0 λ d ( B ′ ) > 0 by Balls Have Positive Lebesgue Measure and Bounded Sets Have Finite Lebesgue Measure (claim 1). Let p p p be the positive periodic density of Heat Smoothing on the Torus: Distance to the Identity, Contraction, a Smooth Positive Periodic Density, Finite Entropy, and Finite Fisher Information §density , so that ν ( B ′ ) = ∫ 1 B ′ 1 Q p d λ d = ∫ 1 B ′ p d λ d \nu(B')=\int\mathbf{1}_{B'}\mathbf{1}_{Q}\,p\,d\lambda_{d}=\int\mathbf{1}_{B'}\,p\,d\lambda_{d} ν ( B ′ ) = ∫ 1 B ′ 1 Q p d λ d = ∫ 1 B ′ p d λ d . If this integral were 0 0 0 , then 1 B ′ p = 0 \mathbf{1}_{B'}p=0 1 B ′ p = 0 almost everywhere by The Lebesgue Integral and Null Sets: Almost-Everywhere Comparison, Markov's Inequality, and Dominated Convergence Almost Everywhere §vanishing , hence 1 B ′ = 0 \mathbf{1}_{B'}=0 1 B ′ = 0 almost everywhere since p > 0 p>0 p > 0 everywhere, hence λ d ( B ′ ) = ∫ 1 B ′ d λ d = 0 \lambda_{d}(B')=\int\mathbf{1}_{B'}\,d\lambda_{d}=0 λ d ( B ′ ) = ∫ 1 B ′ d λ d = 0 by the same clause, a contradiction. So 0 < ν ( B ′ ) ≤ ν ( D ) = 0 0<\nu(B')\le\nu(D)=0 0 < ν ( B ′ ) ≤ ν ( D ) = 0 , which is absurd. Therefore D D D is empty: g ( S s μ , ⋅ ) = h ( S s μ , ⋅ ) g(S_{s}\mu,\cdot)=h(S_{s}\mu,\cdot) g ( S s μ , ⋅ ) = h ( S s μ , ⋅ ) on R d \mathbb{R}^{d} R d .
Step 2 (general measures). Let μ ∈ P ( T d ) \mu\in\mathcal{P}(\mathbb{T}^{d}) μ ∈ P ( T d ) and x ∈ R d x\in\mathbb{R}^{d} x ∈ R d . For n ∈ N n\in\mathbb{N} n ∈ N put s n = 1 / ( 2 ( n + 1 ) ) ∈ ( 0 , 1 2 ] s_{n}=1/(2(n+1))\in(0,\tfrac12] s n = 1/ ( 2 ( n + 1 )) ∈ ( 0 , 2 1 ] and μ n = S s n μ \mu_{n}=S_{s_{n}}\mu μ n = S s n μ . By Heat Smoothing on the Torus: Distance to the Identity, Contraction, a Smooth Positive Periodic Density, Finite Entropy, and Finite Fisher Information §distance , W T ( μ n , μ ) 2 ≤ d s n W_{\mathbb{T}}(\mu_{n},\mu)^{2}\le d\,s_{n} W T ( μ n , μ ) 2 ≤ d s n , which tends to 0 0 0 ; so W T ( μ n , μ ) 2 → 0 W_{\mathbb{T}}(\mu_{n},\mu)^{2}\to0 W T ( μ n , μ ) 2 → 0 by claim 3 of Order Properties of Limits of Real Sequences (with L = 0 L=0 L = 0 ), and hence W T ( μ n , μ ) = ( W T ( μ n , μ ) 2 ) 1 / 2 → 0 1 / 2 = 0 W_{\mathbb{T}}(\mu_{n},\mu)=\bigl(W_{\mathbb{T}}(\mu_{n},\mu)^{2}\bigr)^{1/2}\to0^{1/2}=0 W T ( μ n , μ ) = ( W T ( μ n , μ ) 2 ) 1/2 → 0 1/2 = 0 by Properties of Real Powers of Nonnegative Real Numbers §continuity (the first equality by Properties of Real Powers of Nonnegative Real Numbers §agreement and Properties of Real Powers of Nonnegative Real Numbers §exponents , the last by Properties of Real Powers of Nonnegative Real Numbers §values ); so ( μ n ) (\mu_{n}) ( μ n ) converges to μ \mu μ . By Gradient Fields and Laplacian Test Functions on the Torus Wasserstein Space §joint-continuity , applied to g g g and to h h h with the constant sequence x n = x x_{n}=x x n = x , for each i ∈ [ d ] i\in[d] i ∈ [ d ] the sequences ( g i ( μ n , x ) ) n (g_{i}(\mu_{n},x))_{n} ( g i ( μ n , x ) ) n and ( h i ( μ n , x ) ) n (h_{i}(\mu_{n},x))_{n} ( h i ( μ n , x ) ) n converge to g i ( μ , x ) g_{i}(\mu,x) g i ( μ , x ) and h i ( μ , x ) h_{i}(\mu,x) h i ( μ , x ) respectively. By Step 1 the two sequences coincide term by term, and a real sequence has only one limit (claim 1 of Uniqueness of Limits and Boundedness of Convergent Real Sequences ), so g i ( μ , x ) = h i ( μ , x ) g_{i}(\mu,x)=h_{i}(\mu,x) g i ( μ , x ) = h i ( μ , x ) . Hence g = h g=h g = h , proving claim 1.
Claim 2. Let φ \varphi φ and g g g be as in the hypotheses of claim 2 of the statement; for μ ∈ P ( T d ) \mu\in\mathcal{P}(\mathbb{T}^{d}) μ ∈ P ( T d ) we have ∇ φ ( μ ) = [ g ( μ , ⋅ ) ] \nabla\varphi(\mu)=[g(\mu,\cdot)] ∇ φ ( μ ) = [ g ( μ , ⋅ )] by Gradient Fields and Laplacian Test Functions on the Torus Wasserstein Space §represents . We verify the three clauses of Intrinsic Test Functions on the Torus Wasserstein Space §test with K = P ( T d ) \mathcal{K}=\mathcal{P}(\mathbb{T}^{d}) K = P ( T d ) .
Step 3 (continuity). Let μ ∈ P ( T d ) \mu\in\mathcal{P}(\mathbb{T}^{d}) μ ∈ P ( T d ) and ε > 0 \varepsilon>0 ε > 0 , and put η = ∇ φ ( μ ) \eta=\nabla\varphi(\mu) η = ∇ φ ( μ ) . By Gradient Fields and Laplacian Test Functions on the Torus Wasserstein Space §represents and Differentiability Along Couplings of a Function on the Torus Wasserstein Space, and Its Gradient §differentiable , applied with the value 1 1 1 in place of ε \varepsilon ε , there is θ > 0 \theta>0 θ > 0 such that ∣ φ ( ν ) − φ ( μ ) − J T ( η , γ ) ∣ ≤ I T ( γ ) |\varphi(\nu)-\varphi(\mu)-\mathcal{J}_{\mathbb{T}}(\eta,\gamma)|\le\sqrt{I_{\mathbb{T}}(\gamma)} ∣ φ ( ν ) − φ ( μ ) − J T ( η , γ ) ∣ ≤ I T ( γ ) whenever γ ∈ Π ( μ , ν ) \gamma\in\Pi(\mu,\nu) γ ∈ Π ( μ , ν ) and I T ( γ ) < θ 2 I_{\mathbb{T}}(\gamma)<\theta^{2} I T ( γ ) < θ 2 . Let δ \delta δ be the smaller of θ \theta θ and ε / ( ∥ η ∥ μ + 2 ) \varepsilon/(\lVert\eta\rVert_{\mu}+2) ε / (∥ η ∥ μ + 2 ) , and let ν ∈ P ( T d ) \nu\in\mathcal{P}(\mathbb{T}^{d}) ν ∈ P ( T d ) with W T ( μ , ν ) < δ W_{\mathbb{T}}(\mu,\nu)<\delta W T ( μ , ν ) < δ . By The Torus Wasserstein Space is a Sequentially Compact Metric Space in which Optimal Couplings Exist §optimal there is an optimal γ ∈ Π ( μ , ν ) \gamma\in\Pi(\mu,\nu) γ ∈ Π ( μ , ν ) , so I T ( γ ) = W T ( μ , ν ) 2 < θ 2 I_{\mathbb{T}}(\gamma)=W_{\mathbb{T}}(\mu,\nu)^{2}<\theta^{2} I T ( γ ) = W T ( μ , ν ) 2 < θ 2 by Probability Measures on the Flat Torus, the Torus Cost of a Coupling, the Torus Wasserstein Distance and Optimal Couplings §optimal . By the triangle inequality and The Torus Displacement Pairing: Linearity, the Cost Bound, and Vanishing of a Field with First-Order Small Pairings §bound ,
∣ φ ( ν ) − φ ( μ ) ∣ ≤ ∣ J T ( η , γ ) ∣ + W T ( μ , ν ) ≤ ( ∥ η ∥ μ + 1 ) W T ( μ , ν ) < ( ∥ η ∥ μ + 1 ) δ < ε . |\varphi(\nu)-\varphi(\mu)|\le|\mathcal{J}_{\mathbb{T}}(\eta,\gamma)|+W_{\mathbb{T}}(\mu,\nu)\le(\lVert\eta\rVert_{\mu}+1)\,W_{\mathbb{T}}(\mu,\nu)<(\lVert\eta\rVert_{\mu}+1)\,\delta<\varepsilon . ∣ φ ( ν ) − φ ( μ ) ∣ ≤ ∣ J T ( η , γ ) ∣ + W T ( μ , ν ) ≤ (∥ η ∥ μ + 1 ) W T ( μ , ν ) < (∥ η ∥ μ + 1 ) δ < ε .
So φ \varphi φ is continuous at every μ \mu μ , which is Intrinsic Test Functions on the Torus Wasserstein Space §continuity .
Step 4 (differentiability). For every μ \mu μ , φ \varphi φ is differentiable along couplings at μ \mu μ by Gradient Fields and Laplacian Test Functions on the Torus Wasserstein Space §represents , and ∇ φ ( μ ) = [ g ( μ , ⋅ ) ] ∈ T μ \nabla\varphi(\mu)=[g(\mu,\cdot)]\in T_{\mu} ∇ φ ( μ ) = [ g ( μ , ⋅ )] ∈ T μ by hypothesis. This is Intrinsic Test Functions on the Torus Wasserstein Space §differentiability .
Step 5 (uniform convergence of the fields). Let ( μ n ) (\mu_{n}) ( μ n ) converge to μ \mu μ in ( P ( T d ) , W T ) (\mathcal{P}(\mathbb{T}^{d}),W_{\mathbb{T}}) ( P ( T d ) , W T ) . We show: for every ε > 0 \varepsilon>0 ε > 0 there is N ∈ N N\in\mathbb{N} N ∈ N with ∥ g ( μ n , x ) − g ( μ , x ) ∥ < ε \lVert g(\mu_{n},x)-g(\mu,x)\rVert<\varepsilon ∥ g ( μ n , x ) − g ( μ , x )∥ < ε for all n ≥ N n\ge N n ≥ N and all x ∈ R d x\in\mathbb{R}^{d} x ∈ R d . Suppose not. Then there is ε > 0 \varepsilon>0 ε > 0 such that for every N N N there are n ≥ N n\ge N n ≥ N and x ∈ R d x\in\mathbb{R}^{d} x ∈ R d with ∥ g ( μ n , x ) − g ( μ , x ) ∥ ≥ ε \lVert g(\mu_{n},x)-g(\mu,x)\rVert\ge\varepsilon ∥ g ( μ n , x ) − g ( μ , x )∥ ≥ ε . Choosing successively N = 0 N=0 N = 0 and then N = n k − 1 + 1 N=n_{k-1}+1 N = n k − 1 + 1 , we obtain a strictly increasing sequence ( n k ) k ∈ N (n_{k})_{k\in\mathbb{N}} ( n k ) k ∈ N in N \mathbb{N} N and points x k x_{k} x k with ∥ g ( μ n k , x k ) − g ( μ , x k ) ∥ ≥ ε \lVert g(\mu_{n_{k}},x_{k})-g(\mu,x_{k})\rVert\ge\varepsilon ∥ g ( μ n k , x k ) − g ( μ , x k )∥ ≥ ε . Put y k = π ( x k ) y_{k}=\pi(x_{k}) y k = π ( x k ) ; by The Half-Open Unit Cell Tiles Euclidean Space §wrap , y k ∈ Q y_{k}\in Q y k ∈ Q and x k − y k ∈ Z d x_{k}-y_{k}\in\mathbb{Z}^{d} x k − y k ∈ Z d , so by periodicity of the components (Gradient Fields and Laplacian Test Functions on the Torus Wasserstein Space §regular ) ∥ g ( μ n k , y k ) − g ( μ , y k ) ∥ ≥ ε \lVert g(\mu_{n_{k}},y_{k})-g(\mu,y_{k})\rVert\ge\varepsilon ∥ g ( μ n k , y k ) − g ( μ , y k )∥ ≥ ε . Every point of Q Q Q has coordinates in [ 0 , 1 ) [0,1) [ 0 , 1 ) , hence norm at most d \sqrt{d} d , so Q Q Q is bounded; by Bolzano-Weierstrass Theorem in Euclidean Space there are y ∗ ∈ R d y^{*}\in\mathbb{R}^{d} y ∗ ∈ R d and a strictly increasing sequence ( k j ) j ∈ N (k_{j})_{j\in\mathbb{N}} ( k j ) j ∈ N with y k j → y ∗ y_{k_{j}}\to y^{*} y k j → y ∗ . The sequence ( n k j ) j (n_{k_{j}})_{j} ( n k j ) j is strictly increasing, so ( μ n k j ) j (\mu_{n_{k_{j}}})_{j} ( μ n k j ) j converges to μ \mu μ by A Subsequence of a Convergent Sequence Has the Same Limit . By Gradient Fields and Laplacian Test Functions on the Torus Wasserstein Space §joint-continuity , applied to the sequences ( μ n k j ) j (\mu_{n_{k_{j}}})_{j} ( μ n k j ) j and ( y k j ) j (y_{k_{j}})_{j} ( y k j ) j , and again to the constant sequence ( μ ) j (\mu)_{j} ( μ ) j and ( y k j ) j (y_{k_{j}})_{j} ( y k j ) j , for each i ∈ [ d ] i\in[d] i ∈ [ d ] both g i ( μ n k j , y k j ) g_{i}(\mu_{n_{k_{j}}},y_{k_{j}}) g i ( μ n k j , y k j ) and g i ( μ , y k j ) g_{i}(\mu,y_{k_{j}}) g i ( μ , y k j ) converge to g i ( μ , y ∗ ) g_{i}(\mu,y^{*}) g i ( μ , y ∗ ) as j → ∞ j\to\infty j → ∞ . Hence ∥ g ( μ n k j , y k j ) − g ( μ , y k j ) ∥ ≤ ∑ i ∣ g i ( μ n k j , y k j ) − g i ( μ , y k j ) ∣ \lVert g(\mu_{n_{k_{j}}},y_{k_{j}})-g(\mu,y_{k_{j}})\rVert\le\sum_{i}|g_{i}(\mu_{n_{k_{j}}},y_{k_{j}})-g_{i}(\mu,y_{k_{j}})| ∥ g ( μ n k j , y k j ) − g ( μ , y k j )∥ ≤ ∑ i ∣ g i ( μ n k j , y k j ) − g i ( μ , y k j ) ∣ , which tends to 0 0 0 by claims 1 and 3 of Arithmetic of Limits of Real Sequences and claim 4 of Order Properties of Limits of Real Sequences , contradicting the lower bound ε \varepsilon ε .
Step 6 (gradient continuity). Let μ ∈ P ( T d ) \mu\in\mathcal{P}(\mathbb{T}^{d}) μ ∈ P ( T d ) , let ( μ n ) (\mu_{n}) ( μ n ) be a sequence in P ( T d ) \mathcal{P}(\mathbb{T}^{d}) P ( T d ) , and let γ n ∈ Π ( μ n , μ ) \gamma_{n}\in\Pi(\mu_{n},\mu) γ n ∈ Π ( μ n , μ ) with I T ( γ n ) → 0 I_{\mathbb{T}}(\gamma_{n})\to0 I T ( γ n ) → 0 . Since W T ( μ n , μ ) 2 ≤ I T ( γ n ) W_{\mathbb{T}}(\mu_{n},\mu)^{2}\le I_{\mathbb{T}}(\gamma_{n}) W T ( μ n , μ ) 2 ≤ I T ( γ n ) by Probability Measures on the Flat Torus, the Torus Cost of a Coupling, the Torus Wasserstein Distance and Optimal Couplings §distance , ( μ n ) (\mu_{n}) ( μ n ) converges to μ \mu μ . Use the representatives g ( μ n , ⋅ ) g(\mu_{n},\cdot) g ( μ n , ⋅ ) and g ( μ , ⋅ ) g(\mu,\cdot) g ( μ , ⋅ ) , and let A n A_{n} A n denote the n n n th term of the sequence in Intrinsic Test Functions on the Torus Wasserstein Space §gradient-continuity . By Step 0 there is L ≥ 0 L\ge0 L ≥ 0 with ∥ g ( μ , y ) − g ( μ , x ) ∥ ≤ L d T ( x , y ) \lVert g(\mu,y)-g(\mu,x)\rVert\le L\,d_{\mathbb{T}}(x,y) ∥ g ( μ , y ) − g ( μ , x )∥ ≤ L d T ( x , y ) for all x , y x,y x , y . Let ε > 0 \varepsilon>0 ε > 0 . First choose N 1 N_{1} N 1 as in Step 5 for this ε \varepsilon ε ; then, since I T ( γ n ) → 0 I_{\mathbb{T}}(\gamma_{n})\to0 I T ( γ n ) → 0 , choose N 2 N_{2} N 2 with L 2 I T ( γ n ) < ε 2 L^{2}I_{\mathbb{T}}(\gamma_{n})<\varepsilon^{2} L 2 I T ( γ n ) < ε 2 for n ≥ N 2 n\ge N_{2} n ≥ N 2 . For n n n at least the larger of N 1 N_{1} N 1 and N 2 N_{2} N 2 and w ∈ R d + d w\in\mathbb{R}^{d+d} w ∈ R d + d , with x = p r 1 ( w ) x=\mathrm{pr}_{1}(w) x = pr 1 ( w ) and y = p r 2 ( w ) y=\mathrm{pr}_{2}(w) y = pr 2 ( w ) , the second inequality of Pairs of Euclidean Points: Coordinate Projections, Pairings, the Product Measure on a Euclidean Space, Borel Norm Functions and Finite Sets §functions (applied to g ( μ n , x ) − g ( μ , y ) g(\mu_{n},x)-g(\mu,y) g ( μ n , x ) − g ( μ , y ) and g ( μ n , x ) − g ( μ , x ) g(\mu_{n},x)-g(\mu,x) g ( μ n , x ) − g ( μ , x ) ) gives
∥ g ( μ n , x ) − g ( μ , y ) ∥ 2 ≤ 2 ∥ g ( μ n , x ) − g ( μ , x ) ∥ 2 + 2 ∥ g ( μ , x ) − g ( μ , y ) ∥ 2 ≤ 2 ε 2 + 2 L 2 d T ( x , y ) 2 . \lVert g(\mu_{n},x)-g(\mu,y)\rVert^{2}\le2\lVert g(\mu_{n},x)-g(\mu,x)\rVert^{2}+2\lVert g(\mu,x)-g(\mu,y)\rVert^{2}\le2\varepsilon^{2}+2L^{2}d_{\mathbb{T}}(x,y)^{2}. ∥ g ( μ n , x ) − g ( μ , y ) ∥ 2 ≤ 2 ∥ g ( μ n , x ) − g ( μ , x ) ∥ 2 + 2 ∥ g ( μ , x ) − g ( μ , y ) ∥ 2 ≤ 2 ε 2 + 2 L 2 d T ( x , y ) 2 .
Integrating against the probability measure γ n \gamma_{n} γ n , by linearity and monotonicity of the integral (claim 1 of Linearity and Monotonicity of the Lebesgue Integral ) and Probability Measures on the Flat Torus, the Torus Cost of a Coupling, the Torus Wasserstein Distance and Optimal Couplings §cost , 0 ≤ A n ≤ 2 ε 2 + 2 L 2 I T ( γ n ) < 4 ε 2 0\le A_{n}\le2\varepsilon^{2}+2L^{2}I_{\mathbb{T}}(\gamma_{n})<4\varepsilon^{2} 0 ≤ A n ≤ 2 ε 2 + 2 L 2 I T ( γ n ) < 4 ε 2 . As ε > 0 \varepsilon>0 ε > 0 was arbitrary, A n → 0 A_{n}\to0 A n → 0 , which is Intrinsic Test Functions on the Torus Wasserstein Space §gradient-continuity .
Step 7 (conclusion of claim 2). By Steps 3, 4 and 6, φ \varphi φ is an intrinsic test function on P ( T d ) \mathcal{P}(\mathbb{T}^{d}) P ( T d ) , and g g g is a gradient field of φ \varphi φ ; so φ \varphi φ is a Laplacian test function by Gradient Fields and Laplacian Test Functions on the Torus Wasserstein Space §test .
Claim 3. Let χ = a φ + b ψ \chi=a\varphi+b\psi χ = a φ + b ψ , and put c = ∣ a ∣ + ∣ b ∣ + 1 c=|a|+|b|+1 c = ∣ a ∣ + ∣ b ∣ + 1 .
Step 8 (continuity). Let μ ∈ P ( T d ) \mu\in\mathcal{P}(\mathbb{T}^{d}) μ ∈ P ( T d ) and ε > 0 \varepsilon>0 ε > 0 . By Intrinsic Test Functions on the Torus Wasserstein Space §continuity for φ \varphi φ and ψ \psi ψ there are δ 1 , δ 2 > 0 \delta_{1},\delta_{2}>0 δ 1 , δ 2 > 0 with ∣ φ ( ν ) − φ ( μ ) ∣ < ε / c |\varphi(\nu)-\varphi(\mu)|<\varepsilon/c ∣ φ ( ν ) − φ ( μ ) ∣ < ε / c when W T ( μ , ν ) < δ 1 W_{\mathbb{T}}(\mu,\nu)<\delta_{1} W T ( μ , ν ) < δ 1 and ∣ ψ ( ν ) − ψ ( μ ) ∣ < ε / c |\psi(\nu)-\psi(\mu)|<\varepsilon/c ∣ ψ ( ν ) − ψ ( μ ) ∣ < ε / c when W T ( μ , ν ) < δ 2 W_{\mathbb{T}}(\mu,\nu)<\delta_{2} W T ( μ , ν ) < δ 2 . For W T ( μ , ν ) W_{\mathbb{T}}(\mu,\nu) W T ( μ , ν ) below the smaller of δ 1 , δ 2 \delta_{1},\delta_{2} δ 1 , δ 2 , ∣ χ ( ν ) − χ ( μ ) ∣ ≤ ∣ a ∣ ∣ φ ( ν ) − φ ( μ ) ∣ + ∣ b ∣ ∣ ψ ( ν ) − ψ ( μ ) ∣ ≤ ( ∣ a ∣ + ∣ b ∣ ) ε / c < ε |\chi(\nu)-\chi(\mu)|\le|a|\,|\varphi(\nu)-\varphi(\mu)|+|b|\,|\psi(\nu)-\psi(\mu)|\le(|a|+|b|)\varepsilon/c<\varepsilon ∣ χ ( ν ) − χ ( μ ) ∣ ≤ ∣ a ∣ ∣ φ ( ν ) − φ ( μ ) ∣ + ∣ b ∣ ∣ ψ ( ν ) − ψ ( μ ) ∣ ≤ ( ∣ a ∣ + ∣ b ∣ ) ε / c < ε . So χ \chi χ is continuous .
Step 9 (differentiability). Let μ ∈ K \mu\in\mathcal{K} μ ∈ K and put η = a ∇ φ ( μ ) + b ∇ ψ ( μ ) ∈ L 2 ( μ ; R d ) \eta=a\nabla\varphi(\mu)+b\nabla\psi(\mu)\in L^{2}(\mu;\mathbb{R}^{d}) η = a ∇ φ ( μ ) + b ∇ ψ ( μ ) ∈ L 2 ( μ ; R d ) . Let ε > 0 \varepsilon>0 ε > 0 ; by Intrinsic Test Functions on the Torus Wasserstein Space §differentiability and Differentiability Along Couplings of a Function on the Torus Wasserstein Space, and Its Gradient §differentiable , applied with ε / c \varepsilon/c ε / c , there are θ 1 , θ 2 > 0 \theta_{1},\theta_{2}>0 θ 1 , θ 2 > 0 for φ \varphi φ and ψ \psi ψ ; let θ \theta θ be the smaller. For ν ∈ P ( T d ) \nu\in\mathcal{P}(\mathbb{T}^{d}) ν ∈ P ( T d ) and γ ∈ Π ( μ , ν ) \gamma\in\Pi(\mu,\nu) γ ∈ Π ( μ , ν ) with I T ( γ ) < θ 2 I_{\mathbb{T}}(\gamma)<\theta^{2} I T ( γ ) < θ 2 , The Torus Displacement Pairing: Linearity, the Cost Bound, and Vanishing of a Field with First-Order Small Pairings §linear gives
χ ( ν ) − χ ( μ ) − J T ( η , γ ) = a ( φ ( ν ) − φ ( μ ) − J T ( ∇ φ ( μ ) , γ ) ) + b ( ψ ( ν ) − ψ ( μ ) − J T ( ∇ ψ ( μ ) , γ ) ) , \chi(\nu)-\chi(\mu)-\mathcal{J}_{\mathbb{T}}(\eta,\gamma)=a\bigl(\varphi(\nu)-\varphi(\mu)-\mathcal{J}_{\mathbb{T}}(\nabla\varphi(\mu),\gamma)\bigr)+b\bigl(\psi(\nu)-\psi(\mu)-\mathcal{J}_{\mathbb{T}}(\nabla\psi(\mu),\gamma)\bigr), χ ( ν ) − χ ( μ ) − J T ( η , γ ) = a ( φ ( ν ) − φ ( μ ) − J T ( ∇ φ ( μ ) , γ ) ) + b ( ψ ( ν ) − ψ ( μ ) − J T ( ∇ ψ ( μ ) , γ ) ) ,
whose absolute value is at most ( ∣ a ∣ + ∣ b ∣ ) ( ε / c ) I T ( γ ) ≤ ε I T ( γ ) (|a|+|b|)(\varepsilon/c)\sqrt{I_{\mathbb{T}}(\gamma)}\le\varepsilon\sqrt{I_{\mathbb{T}}(\gamma)} ( ∣ a ∣ + ∣ b ∣ ) ( ε / c ) I T ( γ ) ≤ ε I T ( γ ) . Hence χ \chi χ is differentiable along couplings at μ \mu μ with gradient η \eta η , so ∇ χ ( μ ) = a ∇ φ ( μ ) + b ∇ ψ ( μ ) \nabla\chi(\mu)=a\nabla\varphi(\mu)+b\nabla\psi(\mu) ∇ χ ( μ ) = a ∇ φ ( μ ) + b ∇ ψ ( μ ) by the uniqueness in Differentiability Along Couplings of a Function on the Torus Wasserstein Space, and Its Gradient §gradient ; and ∇ χ ( μ ) ∈ T μ \nabla\chi(\mu)\in T_{\mu} ∇ χ ( μ ) ∈ T μ because ∇ φ ( μ ) , ∇ ψ ( μ ) ∈ T μ \nabla\varphi(\mu),\nabla\psi(\mu)\in T_{\mu} ∇ φ ( μ ) , ∇ ψ ( μ ) ∈ T μ and T μ T_{\mu} T μ is a linear subspace by The Torus Tangent Space: Linearity of the Periodic Calculus, Closed Subspace, and Representation of Bounded Functionals on Gradients §subspace .
Step 10 (gradient continuity). Let μ \mu μ , ( μ n ) (\mu_{n}) ( μ n ) in K \mathcal{K} K and γ n ∈ Π ( μ n , μ ) \gamma_{n}\in\Pi(\mu_{n},\mu) γ n ∈ Π ( μ n , μ ) with I T ( γ n ) → 0 I_{\mathbb{T}}(\gamma_{n})\to0 I T ( γ n ) → 0 . Choose Borel representatives G n , G , H n , H G_{n},G,H_{n},H G n , G , H n , H of ∇ φ ( μ n ) , ∇ φ ( μ ) , ∇ ψ ( μ n ) , ∇ ψ ( μ ) \nabla\varphi(\mu_{n}),\nabla\varphi(\mu),\nabla\psi(\mu_{n}),\nabla\psi(\mu) ∇ φ ( μ n ) , ∇ φ ( μ ) , ∇ ψ ( μ n ) , ∇ ψ ( μ ) ; then a G n + b H n aG_{n}+bH_{n} a G n + b H n and a G + b H aG+bH a G + b H represent ∇ χ ( μ n ) \nabla\chi(\mu_{n}) ∇ χ ( μ n ) and ∇ χ ( μ ) \nabla\chi(\mu) ∇ χ ( μ ) by Step 9. For w ∈ R d + d w\in\mathbb{R}^{d+d} w ∈ R d + d , with x = p r 1 ( w ) x=\mathrm{pr}_{1}(w) x = pr 1 ( w ) , y = p r 2 ( w ) y=\mathrm{pr}_{2}(w) y = pr 2 ( w ) , the first inequality of Pairs of Euclidean Points: Coordinate Projections, Pairings, the Product Measure on a Euclidean Space, Borel Norm Functions and Finite Sets §functions (applied to a ( G n ( x ) − G ( y ) ) a(G_{n}(x)-G(y)) a ( G n ( x ) − G ( y )) and − b ( H n ( x ) − H ( y ) ) -b(H_{n}(x)-H(y)) − b ( H n ( x ) − H ( y )) ) gives
∥ ( a G n + b H n ) ( x ) − ( a G + b H ) ( y ) ∥ 2 ≤ 2 a 2 ∥ G n ( x ) − G ( y ) ∥ 2 + 2 b 2 ∥ H n ( x ) − H ( y ) ∥ 2 . \lVert(aG_{n}+bH_{n})(x)-(aG+bH)(y)\rVert^{2}\le2a^{2}\lVert G_{n}(x)-G(y)\rVert^{2}+2b^{2}\lVert H_{n}(x)-H(y)\rVert^{2}. ∥( a G n + b H n ) ( x ) − ( a G + b H ) ( y ) ∥ 2 ≤ 2 a 2 ∥ G n ( x ) − G ( y ) ∥ 2 + 2 b 2 ∥ H n ( x ) − H ( y ) ∥ 2 .
Integrating against γ n \gamma_{n} γ n (claim 1 of Linearity and Monotonicity of the Lebesgue Integral ), the n n n th term for χ \chi χ in Intrinsic Test Functions on the Torus Wasserstein Space §gradient-continuity lies between 0 0 0 and 2 a 2 2a^{2} 2 a 2 times the n n n th term for φ \varphi φ plus 2 b 2 2b^{2} 2 b 2 times the n n n th term for ψ \psi ψ ; the latter two sequences converge to 0 0 0 by that clause for φ \varphi φ and ψ \psi ψ , so, by claims 1 and 3 of Arithmetic of Limits of Real Sequences and the squeeze claim (claim 2) of Order Properties of Limits of Real Sequences , the former does too. With Steps 8 and 9, χ \chi χ is an intrinsic test function on K \mathcal{K} K with ∇ χ ( μ ) = a ∇ φ ( μ ) + b ∇ ψ ( μ ) \nabla\chi(\mu)=a\nabla\varphi(\mu)+b\nabla\psi(\mu) ∇ χ ( μ ) = a ∇ φ ( μ ) + b ∇ ψ ( μ ) .
Step 11 (the Laplacian case). Let φ , ψ \varphi,\psi φ , ψ be Laplacian test functions with gradient fields g , h g,h g , h . By Gradient Fields and Laplacian Test Functions on the Torus Wasserstein Space §test they are intrinsic test functions on P ( T d ) \mathcal{P}(\mathbb{T}^{d}) P ( T d ) , so by Steps 8 to 10 with K = P ( T d ) \mathcal{K}=\mathcal{P}(\mathbb{T}^{d}) K = P ( T d ) , χ \chi χ is one as well. We check that k = a g + b h k=ag+bh k = a g + bh is a gradient field of χ \chi χ . Clause (a): k i ( μ , ⋅ ) = a g i ( μ , ⋅ ) + b h i ( μ , ⋅ ) ∈ C p e r 1 k_{i}(\mu,\cdot)=a\,g_{i}(\mu,\cdot)+b\,h_{i}(\mu,\cdot)\in C^{1}_{\mathrm{per}} k i ( μ , ⋅ ) = a g i ( μ , ⋅ ) + b h i ( μ , ⋅ ) ∈ C per 1 by Elementary Properties of Lattice-Periodic Functions §algebra . Clause (b): χ \chi χ is differentiable along couplings at every μ \mu μ with ∇ χ ( μ ) = a [ g ( μ , ⋅ ) ] + b [ h ( μ , ⋅ ) ] = [ k ( μ , ⋅ ) ] \nabla\chi(\mu)=a[g(\mu,\cdot)]+b[h(\mu,\cdot)]=[k(\mu,\cdot)] ∇ χ ( μ ) = a [ g ( μ , ⋅ )] + b [ h ( μ , ⋅ )] = [ k ( μ , ⋅ )] , by Step 9 and Gradient Fields and Laplacian Test Functions on the Torus Wasserstein Space §represents for g g g and h h h . Clause (c): by claim 1 of Constants, Coordinate Functions, Sums and Products of C k C^k C k Functions on a Euclidean Open Set , ∂ j k i ( μ , x ) = a ∂ j g i ( μ , x ) + b ∂ j h i ( μ , x ) \partial_{j}k_{i}(\mu,x)=a\,\partial_{j}g_{i}(\mu,x)+b\,\partial_{j}h_{i}(\mu,x) ∂ j k i ( μ , x ) = a ∂ j g i ( μ , x ) + b ∂ j h i ( μ , x ) , so if μ n → μ \mu_{n}\to\mu μ n → μ and x n → x x_{n}\to x x n → x , then k i ( μ n , x n ) k_{i}(\mu_{n},x_{n}) k i ( μ n , x n ) and ∂ j k i ( μ n , x n ) \partial_{j}k_{i}(\mu_{n},x_{n}) ∂ j k i ( μ n , x n ) are the corresponding linear combinations of sequences that converge by Gradient Fields and Laplacian Test Functions on the Torus Wasserstein Space §joint-continuity for g g g and h h h , and so converge to k i ( μ , x ) k_{i}(\mu,x) k i ( μ , x ) and ∂ j k i ( μ , x ) \partial_{j}k_{i}(\mu,x) ∂ j k i ( μ , x ) . Hence χ \chi χ is a Laplacian test function with gradient field a g + b h ag+bh a g + bh by Gradient Fields and Laplacian Test Functions on the Torus Wasserstein Space §test . This proves claim 3.