i}E_{i}^{2}}dE_{1}dE_{2}\text{.}$$]{} By introducing $x_{1}=E_{1}-E_{2}$, $x_{2}=E_{1}+E_{2}$, we have$$\begin{aligned}
P\left( s\right) &\propto &2\int_{-\infty }^{\infty }\left\vert
x_{1}\right\vert ^{\nu }\delta \left( s-\left\vert x_{1}\right\vert \right)
e^{-\frac{A}{2}\sum_{i}x_{i}^{2}}dx_{1}dx_{2}  \notag \\
&=&Cs^{\nu }e^{-As^{2}/2}\text{.}\end{aligned}$$The constants $A,C$ can be determined by working out the integral about $%
x_{2}$, but it is more convenient to obtain by imposing the normalization condition$$\int_{0}^{\infty }P\left( s\right) ds=1\text{, }\int_{0}^{\infty }sP\left(
s\right) ds=1\text{.}  \label{equ:normalization}$$From which we can reach to the famous Wigner-Dyson distribution $$P(s)=\left\{
\begin{array}{ll}
\frac{\pi }{2}s\exp \big(-\frac{\pi }{4}s^{2}\big) & \nu =1\quad \text{GOE}
\\[1mm]
\frac{32}{\pi ^{2}}s^{2}\exp \big(-\frac{4}{\pi }s^{2}\big) & \nu =2\quad
\text{GUE} \\[1mm]
\frac{2^{18}}{3^{6}\pi ^{3}}s^{4}\exp \big(-\frac{64}{9\pi }s^{2}\big) & \nu
=4\quad \text{GSE}%
\end{array}%
\right.  \label{equ:nearest}$$

On the other hand, the levels are independent in Poisson ensemble, which means the occurrence of next level is independent of previous level, the nearest level spacings then follows a Poisson distribution $P\left( s\right)
=\exp \left( -s\right) $.

Although the Wigner surm
i}E_{i}^{2}}dE_{1}dE_{2}\text{.}$$ ] { } By introducing $ x_{1}=E_{1}-E_{2}$, $ x_{2}=E_{1}+E_{2}$, we have$$\begin{aligned } 
 P\left (s\right) & \propto & 2\int_{-\infty } ^{\infty } \left\vert 
 x_{1}\right\vert ^{\nu } \delta \left (s-\left\vert x_{1}\right\vert \right) 
 e^{-\frac{A}{2}\sum_{i}x_{i}^{2}}dx_{1}dx_{2 }   \notag \\ 
 & = & Cs^{\nu } e^{-As^{2}/2}\text{.}\end{aligned}$$The constants $ A, C$ can be determined by work out the built-in about $% 
 x_{2}$, but it is more convenient to receive by inflict the normalization condition$$\int_{0}^{\infty } P\left (s\right) ds=1\text {, } \int_{0}^{\infty } sP\left (
 s\right) ds=1\text {. }   \label{equ: normalization}$$From which we can reach to the celebrated Wigner - Dyson distribution $ $ P(s)=\left\ { 
 \begin{array}{ll } 
 \frac{\pi } { 2}s\exp \big(-\frac{\pi } { 4}s^{2}\big) & \nu = 1\quad \text{GOE } 
 \\[1 mm ] 
 \frac{32}{\pi ^{2}}s^{2}\exp \big(-\frac{4}{\pi } s^{2}\big) & \nu = 2\quad 
 \text{GUE } \\[1 mm ] 
 \frac{2^{18}}{3^{6}\pi ^{3}}s^{4}\exp \big(-\frac{64}{9\pi } s^{2}\big) & \nu 
 = 4\quad \text{GSE}% 
 \end{array}% 
 \right.   \label{equ: nearest}$$ 

 On the other hand, the levels are autonomous in Poisson corps de ballet, which means the occurrence of future level is independent of former level, the nearest degree spacings then follows a Poisson distribution $ P\left (s\right) 
 = \exp \left (-s\right) $. 

 Although the Wigner surm
i}E_{i}^{2}}fE_{1}dE_{2}\text{.}$$]{} By introducing $w_{1}=E_{1}-E_{2}$, $x_{2}=E_{1}+E_{2}$, we have$$\bggun{aligied}
P\leff( s\right) &\propto &2\int_{-\infty }^{\infty }\left\vxrt
x_{1}\eight\cert ^{\nu }\delta \left( s-\lewt\vert x_{1}\rpght\vert \eighu)
e^{-\frac{A}{2}\sum_{i}x_{i}^{2}}dx_{1}dx_{2}  \notag \\
&=&Gf^{\nu }s^{-Ws^{2}/2}\tert{.}\xnd{aligned}$$The cpnstants $A,W$ can be deterkiveb by working out the integral about $%
v_{2}$, but iy ls more convenyent eo ognaln by imposing the normalization conditpon$$\int_{0}^{\infty }P\left( s\right) ds=1\text{, }\int_{0}^{\infty }sP\peft(
d\right) ds=1\text{.}  \labfl{equ:normaluzatyin}$$From which we can reach to the fzmous Wigner-Dyson distribution $$O(s)=\lefc\{
\begin{array}{ol}
\draf{\ki }{2}s\exp \big(-\fcac{\pi }{4}f^{2}\big) & \nu =1\quad \text{GOA}
\\[1mm]
\frac{32}{\li ^{2}}s^{2}\exp \big(-\frag{4}{\pi }s^{2}\uig) & \nu =2\quad
\text{GUE} \\[1mm]
\frar{2^{18}}{3^{6}\pi ^{3}}s^{4}\exp \big(-\frac{64}{9\pi }s^{2}\fig) & \nu
=4\qudd \text{GSE}%
\end{array}%
\rugyt.  \ldbel{aqu:ndqrert}$$

Oh vhe other haid, the levems are independent in Poisson tnsqnble, which mezns thq jccurrence of next level is independent of previous level, the neaeest level spacings tjen follors a Poisson distribution $P\left( s\right)
=\exp \left( -s\rhght) $.

Elghobnm thd Wlgner surm
i}E_{i}^{2}}dE_{1}dE_{2}\text{.}$$]{} By introducing $x_{1}=E_{1}-E_{2}$, $x_{2}=E_{1}+E_{2}$, we have$$\begin{aligned} &\propto }^{\infty }\left\vert ^{\nu }\delta \left( \\ }e^{-As^{2}/2}\text{.}\end{aligned}$$The constants $A,C$ be determined by out the integral about $% x_{2}$, it is more convenient to obtain by imposing the normalization condition$$\int_{0}^{\infty }P\left( s\right) }\int_{0}^{\infty }sP\left( s\right) ds=1\text{.} \label{equ:normalization}$$From which we can reach to the famous Wigner-Dyson $$P(s)=\left\{ \frac{\pi \big(-\frac{\pi & \nu =1\quad \text{GOE} \\[1mm] \frac{32}{\pi ^{2}}s^{2}\exp \big(-\frac{4}{\pi }s^{2}\big) & \nu =2\quad \text{GUE} \\[1mm] \frac{2^{18}}{3^{6}\pi ^{3}}s^{4}\exp }s^{2}\big) & \nu =4\quad \text{GSE}% \end{array}% \right. \label{equ:nearest}$$ the other hand, the are independent in Poisson ensemble, means occurrence of level independent previous level, the level spacings then follows a Poisson distribution $P\left( s\right) =\exp \left( -s\right) $. Although the Wigner surm
i}E_{i}^{2}}dE_{1}dE_{2}\text{.}$$]{} By introducing $x_{1}=e_{1}-E_{2}$, $x_{2}=E_{1}+E_{2}$, we havE$$\begiN{alIgnEd}
p\lefT( s\riGht) &\propto &2\int_{-\inFTy }^{\inFty }\left\vert
x_{1}\right\vert ^{\nU }\deltA \lEFt( s-\lEFt\Vert x_{1}\Right\veRT \rIGHt)
e^{-\FrAc{a}{2}\suM_{i}X_{I}^{2}}dX_{1}dx_{2}  \noTag \\
&=&cs^{\nu }e^{-As^{2}/2}\Text{.}\end{aliGneD}$$THe constants $A,c$ CaN be determiNed By working out The IntegrAl AboUT $%
x_{2}$, but It iS more ConvenIEnt to oBtain by imPoSIng the NOrmalizATIoN conDition$$\int_{0}^{\infty }P\leFT( s\RIght) ds=1\text{, }\int_{0}^{\iNfty }sP\LeFT(
s\RIGht) Ds=1\tExt{.}  \label{eqU:nOrmalIZation}$$FROm WHICh wE Can reach to the Famous WigneR-dysOn distRiButIOn $$P(s)=\leFt\{
\begIn{ARraY}{ll}
\frac{\pi }{2}s\eXp \biG(-\frac{\pi }{4}s^{2}\bIg) & \nu =1\quAD \text{GOe}
\\[1Mm]
\frac{32}{\pI ^{2}}s^{2}\exp \bIg(-\fRac{4}{\Pi }s^{2}\bIG) & \nU =2\qUad
\TeXT{GUe} \\[1Mm]
\FraC{2^{18}}{3^{6}\Pi ^{3}}s^{4}\Exp \big(-\frAc{64}{9\Pi }S^{2}\big) & \nU
=4\quaD \TEXT{GSE}%
\End{ArraY}%
\righT.  \label{equ:nearEst}$$

on thE OthEr hanD, the lEvelS aRe indEpendeNt in POiSson ensemble, whiCh meAns the occUrrEnCe oF nExt leVEl is inDepEndEnt of prEvious lEVel, ThE NEArEst level spacings thEn FOLlOws a PoisSon disTRiBuTIon $P\left( S\rIghT)
=\exp \LEFt( -s\riGht) $.

ALThOugh the WIgner sURm
i}E_{i}^{2}}dE_{1}dE_{2}\t ext{.}$$]{ } Byint rod uc ing$x_{ 1}=E_{1}-E_{2} $ , $x _{2}=E_{1}+E_{2}$, wehave$ $\ b egin { al igned }
P\lef t (s \ rig ht )&\p ro p to  &2\i nt_ {-\inft y }^{\inft y } \l eft\vert
x_{ 1 }\ right\vert  ^{ \nu }\delta\le ft( s- \l eft \ vertx_{ 1}\ri ght\ve r t \rig ht)
e^{-\ fr a c{A}{2 } \sum_{i } x _{ i}^{ 2}}dx_{1}dx_{2} \ no t ag \\
&=&Cs^{\ nu }e^ {- A s^ { 2 }/2 }\t ext{.}\end {a ligne d }$$Thec on s t a nts $A,C$ can bedeterminedb y w orking  o utt he int egral  a b out  $%
x_{2}$,  but  it is mo re con v enientt o obtai n by i mpo sin g th e  n or mal iz a tio n  c ond i tio n$$\int_ {0 }^ {\inf ty } P \ l e ft(s\r ight ) ds= 1\text{, }\in t_{ 0}^{ \ inf ty }s P\lef t(
s \r ight)  ds=1\ text{ .}   \label{equ:no rmal ization}$ $Fr om  wh ic h wec an rea chtothe fam ous Wig n er- Dy s o n  d istribution $$P(s) =\ l e ft \{
\begi n{arra y }{ ll } 
\frac{\ pi  }{ 2}s\ e x p \bi g(-\ f ra c{\pi }{ 4}s^{2 } \b ig ) & \nu  = 1\quad  \ tex t{G OE}
\ \ [1mm ]
\fra c{32}{\p i ^{2 } }s^{2}\exp \bi g (-\frac{4}{\p i  } s ^ {2 } \big ) &  \nu =2\qua d
\t e xt{G UE}\ \[ 1mm ] 
\fra c{2^{ 18 } }{ 3 ^{6}\pi ^{3}}s^{4}\ ex p \big (-\fr ac{64}{9\pi } s^{2}\big) & \nu
=4\q uad\ te x t{GSE}%
\end{a rray} %
\right.\label{e qu:ne arest}$$ 

On theo t her hand , t helev els a re  independenti n  Poi ss on ense mbl e, whic h m ean s t heoc currenceof nextle ve lis  in depen d ent of p re vio us  le vel,t he nea restleve lsp a cin gs then fo l l owsaPo isso n d is tribu tion $P\ left( s \right)
= \ex p  \le ft (-s\righ t) $.

Althou gh  the Wigne rsur m
i}E_{i}^{2}}dE_{1}dE_{2}\text{.}$$]{} By_introducing $x_{1}=E_{1}-E_{2}$,_$x_{2}=E_{1}+E_{2}$, we have$$\begin{aligned}
P\left( s\right)_&\propto &2\int_{-\infty_}^{\infty_}\left\vert
x_{1}\right\vert ^{\nu_}\delta_\left( s-\left\vert x_{1}\right\vert_\right)
e^{-\frac{A}{2}\sum_{i}x_{i}^{2}}dx_{1}dx_{2}  \notag_\\
&=&Cs^{\nu }e^{-As^{2}/2}\text{.}\end{aligned}$$The constants $A,C$_can be determined_by_working out the integral about $%
x_{2}$, but it is more convenient to obtain by_imposing_the normalization_condition$$\int_{0}^{\infty_}P\left(_s\right) ds=1\text{, }\int_{0}^{\infty }sP\left(
s\right) ds=1\text{.}_ \label{equ:normalization}$$From which we can_reach to_the famous Wigner-Dyson distribution $$P(s)=\left\{
\begin{array}{ll}
\frac{\pi }{2}s\exp \big(-\frac{\pi }{4}s^{2}\big)_&_\nu =1\quad \text{GOE}
\\[1mm]
\frac{32}{\pi_^{2}}s^{2}\exp \big(-\frac{4}{\pi }s^{2}\big) & \nu =2\quad
\text{GUE} \\[1mm]
\frac{2^{18}}{3^{6}\pi ^{3}}s^{4}\exp \big(-\frac{64}{9\pi_}s^{2}\big) & \nu
=4\quad \text{GSE}%
\end{array}%
\right.  \label{equ:nearest}$$

On_the other hand,_the_levels_are independent in Poisson_ensemble, which means the occurrence of_next level is independent of previous_level, the nearest level spacings then follows_a Poisson distribution $P\left( s\right)
=\exp \left(_-s\right) $.

Although the Wigner surm