Throughout, B2=B(R)⊗B(R) is the product σ-algebra on R2, and a function on R2 is called jointly Borel if it is measurable from (R2,B2) to (R,B(R)), as in Joint Distribution, Expectations, and Block Independence for Independent Random Variables. Let g(x)=exp(−x2/2) and c=∫Rgdλ∈(0,∞) be as in Standard Normal Distribution, with λ Lebesgue measure, which is σ-finite, so the product measure λ⊗λ on B2 is defined. Set G(x,y)=g(x)g(y).
Step 0: Measurability inventory. By Claim 4 of Joint Distribution, Expectations, and Block Independence for Independent Random Variables, the coordinate projections π1,π2 are jointly Borel, and t∘φ is jointly Borel whenever φ is jointly Borel and t is Borel; since g is Borel (recorded in Standard Normal Distribution), g∘π1 and g∘π2 are jointly Borel, and their product G is jointly Borel by the closure of measurable functions under products recorded in the preliminaries of Square-Integrable Random Variables and the Mean-Square Inner Product, applied on the measurable space (R2,B2). Next, a map T:R2→R2 whose two coordinate functions are jointly Borel is measurable from (R2,B2) to itself: the preimage of a measurable rectangle B×B′ is the intersection of the coordinate preimages, which lie in B2, and rectangles generate B2 by Product Sigma-Algebra, so the generator criterion applies. Finally, for a real α the function hα(x,y)=x+αy is jointly Borel: the map φα(x,y)=(x,αy) has jointly Borel coordinates (π1, and t∘π2 with the Borel map t(u)=αu), hence is measurable by the rectangle argument, and hα is the composition of the jointly Borel addition map of Claim 4 of Joint Distribution, Expectations, and Block Independence for Independent Random Variables with φα, and compositions of measurable maps are measurable (preimages compose).
Step 1: The joint law of (Z1,Z2) has density G. By Claim 1 of Joint Distribution, Expectations, and Block Independence for Independent Random Variables, the map V=(Z1,Z2):Ω→R2 is measurable and its distribution is the product measure N⊗N, where N is the standard normal distribution. Define ρ:B2→[0,∞] by
ρ(C)=c−2∫R21CGd(λ⊗λ).
Then ρ is a measure: ρ(∅)=0, and for pairwise disjoint (Ck)k∈N the partial sums of ∑k1CkG increase pointwise to 1∪kCkG, so countable additivity follows from Monotone Convergence Theorem together with the additivity of the integral in Linearity and Monotonicity of the Lebesgue Integral. On a measurable rectangle B×B′, the Tonelli theorem and the fact that constants factor out of the integral (by Linearity and Monotonicity of the Lebesgue Integral) give
ρ(B×B′)=c−2∫R1B′(y)g(y)(∫R1B(x)g(x)dλ(x))dλ(y)=N(B)N(B′).
Since ρ and N⊗N are two measures on B2 assigning N(B)N(B′) to every measurable rectangle, the uniqueness assertion of Existence and Uniqueness of the Product Measure yields ρ=N⊗N. Hence
P(V∈C)=c−2∫R21CGd(λ⊗λ)(C∈B2).(1)
Step 2: Shears preserve the integral. For a real α define the shears Sα(x,y)=(x+αy,y) and Sα′(x,y)=(x,y+αx); by Step 0 both are measurable from (R2,B2) to itself. We claim that for every B2-measurable f:R2→R with f≥0,
∫R2f∘Sαd(λ⊗λ)=∫R2fd(λ⊗λ),
and likewise for Sα′. Indeed, f∘Sα is B2-measurable as a composition, and by the Tonelli theorem, integrating first in x,
∫R2f∘Sαd(λ⊗λ)=∫R(∫Rf(x+αy,y)dλ(x))dλ(y).
For each fixed y, the inner integrand is the translate by t=αy of the section x↦f(x,y), which is measurable by the section assertion of Tonelli and Fubini Theorems; Claim 2 of Translation Invariance of Lebesgue Measure and the Lebesgue Integral shows the inner integral equals ∫Rf(x,y)dλ(x). Applying the Tonelli theorem once more gives the claim. The claim for Sα′ is proved symmetrically, integrating first in y and using the other iterated-integral identity of Tonelli and Fubini Theorems.
Step 3: Factorization of the rotation into shears. Let R(x,y)=(ax+by,−bx+ay); its coordinates are jointly Borel by Step 0 (sums of the jointly Borel functions t∘πi are jointly Borel, again by the preliminaries of Square-Integrable Random Variables and the Mean-Square Inner Product on (R2,B2)), so R is measurable. Suppose first b=0, and set
α=γ=b1−a,β=−b.
Then R=Sα∘Sβ′∘Sγ. Indeed, applying Sγ, then Sβ′, then Sα to (x,y) yields the pair with second coordinate βx+(1+βγ)y and first coordinate (1+αβ)x+(γ+α(1+βγ))y; here 1+αβ=1+βγ=1−(1−a)=a, and
γ+αa=b(1−a)(1+a)=b1−a2=bb2=b
using a2+b2=1, so the pair is exactly (ax+by,−bx+ay). Consequently, for B2-measurable f≥0, writing f∘R=((f∘Sα)∘Sβ′)∘Sγ and applying Step 2 three times,
∫R2f∘Rd(λ⊗λ)=∫R2fd(λ⊗λ).(2)
Step 4: The density is rotation invariant. By the multiplicative property of the exponential function (Basic Properties of the Exponential Function), G(x,y)=exp(−(x2+y2)/2), and
(ax+by)2+(−bx+ay)2=(a2+b2)(x2+y2)=x2+y2,
since the cross terms 2abxy and −2abxy cancel. Hence G∘R=G pointwise on R2.
Step 5: Conclusion for b=0. The functions W1,W2 are random variables by the closure preliminaries of Square-Integrable Random Variables and the Mean-Square Inner Product, and (W1,W2)=R∘V pointwise on Ω. Fix Borel sets B,B′ and put C=R−1(B×B′)∈B2. Then {W1∈B}∩{W2∈B′}=V−1(C), and 1C=1B×B′∘R, so by Step 4, 1CG=(1B×B′G)∘R. Combining the identity (1) of Step 1, the invariance (2) of Step 3 applied to f=1B×B′G, and Step 1 again,
P(W1∈B,W2∈B′)=c−2∫R2(1B×B′G)∘Rd(λ⊗λ)=c−2∫R21B×B′Gd(λ⊗λ)=N(B)N(B′).
Taking B′=R gives P(W1∈B)=N(B) for every Borel B, so W1 is standard normal; taking B=R gives the same for W2. The displayed identity then reads P(W1∈B,W2∈B′)=P(W1∈B)P(W2∈B′) for all Borel B,B′, which is independence of W1 and W2.
Step 6: The case b=0. Then a2=1. If a=1, then (W1,W2)=(Z1,Z2) and there is nothing to prove. If a=−1, then (W1,W2)=(−Z1,−Z2): each is standard normal by Claim 3 of Reflection Invariance of Lebesgue Measure and Symmetry of the Standard Normal Distribution, and for Borel B,B′ the events {−Z1∈B}={Z1∈−B} and {−Z2∈B′}={Z2∈−B′} involve the Borel sets −B,−B′ by Claim 1 of Reflection Invariance of Lebesgue Measure and Symmetry of the Standard Normal Distribution, so the product identity follows from the independence of Z1 and Z2. ■