Brownian Motion on CUE
Brownian motion on the manifold of the unitary matrices can be constructed by infinitesimal generators from Gaussian unitary ensemble. The stationary distribution of this Brownian motion is then identical to the distribution of CUE.
In[1]:=
![Click for copyable input](assets.en/brownian-motion-on-cue/In_56.png)
mats = RandomVariate[GaussianUnitaryMatrixDistribution[0.1, 2],
100000];
mats = Table[MatrixExp[I mat], {mat, mats}];
Generate a Brownian path with initial point sampled from CUE.
In[2]:=
![Click for copyable input](assets.en/brownian-motion-on-cue/In_57.png)
initial = RandomVariate[CircularUnitaryMatrixDistribution[2]];
res = FoldList[#2.#1 &, initial, mats];
Compute the phases of the eigenvalues and compare them with the PDF of the eigenvalues of matrices from CUE.
In[3]:=
![Click for copyable input](assets.en/brownian-motion-on-cue/In_58.png)
phases = RandomSample /@ Arg[Eigenvalues /@ res];
show complete Wolfram Language input
Out[4]=
![](assets.en/brownian-motion-on-cue/O_27.png)