_\text{th} &\sim N^{-2/3}\\
\overline{\varepsilon^2}_c & \sim N^{-4/3} \text{.}\end{aligned}$$

\[eq:fs-law\]

Studying the large deviation tails, we find that the “soft” tail (corresponding to low energies) has the same asymptotics as in the Tracy-Widom law, but the “hard” tail decays more like a Gaussian and therefore decays considerably more slowly than the Tracy-Widom law, $$\begin{aligned}
\log P_\text{empirical} (x) \sim \begin{cases} x^{\approx 3/2} & x \ll 0 \\ -x^{\approx 2} & x\gg 0 \end{cases}\end{aligned}$$ with again the rescaled variable $x=(\varepsilon-\mu_\varepsilon(N))/\sigma_\varepsilon(N)$. This is shown in Fig. \[fig:fs\_distro\](B). From Fig. \[fig:combo\_pdf\] (a) it is hard to tell whether this holds for large $N$ as the shape of the distribution appears to cross over from something close to the Tracy-Widom distribution towards the Gaussian distribution. To address this question, we construct an estimate for the converged shape for very large $N$ using the following procedure. We sample the inverse function $Q_N(c)=C_N^{-1}(c)$ to the empirical cumulative distribution functions (CDFs) $C_N(x)$ for every $N$ (using the $N$-specific normalization) at a number of selected values $c=10^{-5},\ldots,1$. At any given value of $c$, this gives a set of pairs $(1/N,x_c=Q_N(c))$ which we use to extrapolate to $x_c(0)$. We find that the finite size effects in the shape are well described by $x(N)-x_c(0) \propto N^{-1/3}$ (corresponding to $1/N$ corrections in non-normalized variables). As the extrapolation is done at each value independently, this method does not constrain the moments of the final distribution and we therefore normalize it as a final step. The result of this procedure is shown in Fig. \[fig:extrap\], which suggests that although the final distribution is extremely close to a
_ \text{th } & \sim N^{-2/3}\\ 
 \overline{\varepsilon^2}_c & \sim N^{-4/3 } \text{.}\end{aligned}$$ 

 \[eq: fs - law\ ] 

 Studying the large deviation tails, we witness that the “ piano ” tail (corresponding to low energy) has the same asymptotics as in the Tracy - Widom police, but the “ intemperate ” tail decay more like a Gaussian and therefore decays well more slowly than the Tracy - Widom law, $ $ \begin{aligned } 
 \log P_\text{empirical } (x) \sim \begin{cases } x^{\approx 3/2 } & x \ll 0 \\ -x^{\approx 2 } & x\gg 0 \end{cases}\end{aligned}$$ with again the rescaled variable $ x=(\varepsilon-\mu_\varepsilon(N))/\sigma_\varepsilon(N)$. This is shown in Fig.   \[fig: fs\_distro\](B). From Fig.   \[fig: combo\_pdf\ ] (a) it is unvoiced to tell whether this holds for large $ N$ as the shape of the distribution appears to traverse over from something close to the Tracy - Widom distribution towards the Gaussian distribution. To address this question, we reconstruct an estimate for the converged shape for very large $ N$ use the following procedure. We sample the inverse function $ Q_N(c)=C_N^{-1}(c)$ to the empirical cumulative distribution functions (CDFs) $ C_N(x)$ for every $ N$ (use the $ N$-specific normalization) at a number of selected values $ c=10^{-5},\ldots,1$. At any given value of $ c$, this move over a set of pairs $ (1 / N, x_c = Q_N(c))$ which we use to extrapolate to $ x_c(0)$. We find that the finite size effects in the shape are well described by $ x(N)-x_c(0) \propto N^{-1/3}$ (corresponding to $ 1 / N$ corrections in non - normalized variable). As the extrapolation is done at each value independently, this method does not stiffen the here and now of the final distribution and we therefore normalize it as a concluding step. The resultant role of this procedure is shown in Fig.   \[fig: extrap\ ], which suggest that although the final distribution is extremely close to a
_\texh{th} &\sim N^{-2/3}\\
\overline{\varepsiuon^2}_c & \sim N^{-4/3} \texj{.}\ebd{aligied}$$

\[eq:fs-maw\]

Studykng the large deviation taild, qe fibd that the “soft” tail (zorresponfing to oow tnergies) has the same asymptotidd as mn the Tracy-Widpm law, but the “hard” tail ddccys more like a Gaussian and therefowe decaus considerably iore flowmj uhan the Tracy-Widom law, $$\begin{alighed}
\log K_\text{empirical} (x) \som \begin{cases} x^{\approx 3/2} & x \pl 0 \\ -x^{\approx 2} & x\gg 0 \enf{cases}\end{alugneq}$$ with again ghe rescaltd variable $s=(\varepsilon-\mu_\varepsilon(N))/\sigma_\vafepsikon(N)$. This us shlfn in Fig. \[fij:fs\_disnro\](B). From Fig. \[npg:combo\_[df\] (a) iy is hard to tcll wiethwr this holds for larje $N$ as the shape of the distsiyution appears to crows over frok soowthkng ckoae to hhe Tracy-Widoj distributuon towards the Gauxsyqn distributikn. To wdqress this question, we construct an estpmats for the converged shape for very large $N$ uding the sollowing procedure. We sample the inverse functiot $Q_N(c)=R_N^{-1}(z)$ ti bhe dnplrical cumulative distribution functions (CDFs) $D_N(c)$ nor every $N$ (usinn the $N$-specific notmwloeation) at a nuober oy sslected values $c=10^{-5},\ldlts,1$. At wny guven valut of $v$, this gives a set of pairs $(1/N,x_c=Q_N(c))$ whick ww use to extrapolace to $x_c(0)$. We yind tnat tne finite size effects nn the shape are aell descdkbed by $x(N)-x_c(0) \prooto N^{-1/3}$ (corresponding to $1/N$ correceions in ion-noxmalized varoables). As the exhrapolation is done at ewch vclue hndependenhly, this method does not constremn the momentx mf nhe final distvibution and we therefore norkalize nt as x final stvp. The revult of thif procedure iv shown in Fij. \[fig:extra[\], whuch wuggestr that although the final distriburion is extremely glose fo a
_\text{th} &\sim N^{-2/3}\\ \overline{\varepsilon^2}_c & \sim N^{-4/3} Studying large deviation we find that low has the same as in the law, but the “hard” tail decays like a Gaussian and therefore decays considerably more slowly than the Tracy-Widom law, \log P_\text{empirical} (x) \sim \begin{cases} x^{\approx 3/2} & x \ll 0 \\ -x^{\approx & 0 with the rescaled variable $x=(\varepsilon-\mu_\varepsilon(N))/\sigma_\varepsilon(N)$. This is shown in Fig. \[fig:fs\_distro\](B). From Fig. \[fig:combo\_pdf\] (a) it is to tell whether this holds for large $N$ the shape of the appears to cross over from close the Tracy-Widom towards Gaussian To address this we construct an estimate for the converged shape for very large $N$ using the following procedure. We the inverse to the cumulative functions $C_N(x)$ for every the $N$-specific normalization) at a number $c=10^{-5},\ldots,1$. At any given value of $c$, this a set pairs $(1/N,x_c=Q_N(c))$ which we use to to $x_c(0)$. We find that the finite size in the shape are well described by $x(N)-x_c(0) \propto N^{-1/3}$ (corresponding to $1/N$ corrections in As the extrapolation is at each value this does constrain moments of final distribution and we therefore normalize it as a final step. result of this procedure is shown in Fig. \[fig:extrap\], which although final distribution is close to a
_\text{th} &\sim N^{-2/3}\\
\overline{\varepsiLon^2}_c & \sim N^{-4/3} \teXt{.}\end{AliGneD}$$

\[eQ:fs-lAw\]

StUdying the large DEviaTion tails, we find that the “Soft” tAiL (CorrESpOndinG to low eNErGIEs) hAs ThE saMe ASyMptotIcs As in the tracy-Widom Law, BuT the “hard” tail DEcAys more likE a GAussian and thEreFore deCaYs cONsideRabLy morE slowlY Than thE Tracy-WidOm LAw, $$\begiN{Aligned}
\LOG P_\Text{Empirical} (x) \sim \begiN{CaSEs} x^{\approx 3/2} & x \ll 0 \\ -x^{\aPprox 2} & x\Gg 0 \ENd{CASes}\End{Aligned}$$ witH aGain tHE rescalED vARIAblE $X=(\varepsilon-\mu_\Varepsilon(N))/\SIgmA_\varepSiLon(n)$. this is Shown In fIg. \[fIg:fs\_distro\](B). from fig. \[fig:comBo\_pdf\] (a) IT is hard TO tell whEther tHis HolDs foR LaRgE $N$ aS tHE shAPe Of tHE diStributiOn ApPears To crOSS OVer fRom SomeThing Close to the TraCy-WIdom DIstRibutIon toWardS tHe GauSsian dIstriBuTion. To address thIs quEstion, we cOnsTrUct An EstimATe for tHe cOnvErged shApe for vERy lArGE $n$ UsIng the following proCeDURe. we sample The invERsE fUNction $Q_N(C)=C_n^{-1}(c)$ tO the EMPiricAl cuMUlAtive disTributIOn FuNctions (cDfs) $C_N(x)$ fOr EveRy $N$ (Using THe $N$-sPecifiC normaliZatioN) At a number of selECted values $c=10^{-5},\ldOTs,1$. aT AnY GiveN vaLue of $c$, this gIves A Set oF paiRS $(1/N,X_c=Q_n(C))$ whicH we usE tO ExTRapolate to $x_c(0)$. We find tHaT the fiNite sIze effects in tHe shape are WELL describEd by $X(n)-x_C(0) \Propto N^{-1/3}$ (correspOndinG to $1/N$ correcTIons in noN-normAlized vaRiables). As THE extrapoLatIon Is dOne AT EaCh value indepeNDEntlY, tHis methOd dOes not cOnsTraIn tHe mOmEnts of the Final disTrIbUtIoN anD we thERefore noRmAliZe It aS a finAL step. THe resUlt oF tHiS ProCedure iS ShOWN in FIg. \[FiG:extRap\], WhIch suGgesTS thAt althoUgh the finAl dIStriBuTiOn is extRemely close to A
_\text{th} &\sim N^{-2/3}\ \
\overlin e{\va rep sil on ^2}_ c &\sim N^{-4/3}\ text {.}\end{aligned}$$

\[ eq:fs -l a w\]St udyin g the l a rg e dev ia ti onta i ls , wefin d thatthe “soft”  ta il  (correspond i ng  to low en erg ies) has the  sa me asy mp tot i cs as  in  theTracy- W idom l aw, but t he “hard” tail de c a ys  mor e like a Gaussian an d  therefore dec ays co ns i de r a bly  mo re slowlyth an th e  Tracy- W id o m law ,  $$\begin{ali gned}
\logP _\t ext{em pi ric a l} (x)  \sim  \ b egi n{cases} x^ {\ap prox 3/2}  & x \ l l 0 \\- x^{\app rox 2}  &x\g g 0\ en d{ cas es } \en d {a lig n ed} $$ withag ai n the  res c a l e d va ria ble$x=(\ varepsilon-\m u_\ vare p sil on(N) )/\si gma_ \v areps ilon(N )$. T hi s is shown in F ig.\[fig:fs\ _di st ro\ ]( B). F r om Fig . \ [fi g:combo \_pdf\] (a)  i t i shard to tell wheth er t hi s holdsfor la r ge  $ N $ as the  s hap e of t he di stri b ut ion appe ars to cr os s overfr om som et hin g c loset o th e Trac y-Widomdistr i bution towards the Gaussiand is t r ib u tion . T o address t hisq uest ion, we  co n struc t anes t im a te for the converge dshapefor v ery large $N$  using the f o llowingproc e du r e. We sample t he in verse func t ion $Q_N (c)=C _N^{-1}( c)$ to th e empirica l c umu lat ive d is tribution fun c t ions  ( CDFs) $ C_N (x)$ fo r e ver y $ N$(u sing the$N$-spec if ic  n or mal izati o n) at anu mbe rofselec t ed val ues $ c=10 ^{ -5 } ,\l dots,1$ .  A t anygi ve n va lue  o f $c$ , th i s g ives aset of pa irs $(1/ N, x_ c=Q_N(c ))$ which weus e to extra po lat e to $ x _ c(0)$. W e find that the finites ize eff ect s intheshape are  we ll des cri b ed by$x(N)- x_c(0 )\pr o p to N^ { - 1/ 3}$  ( correspond i n g t o $1/ N$  cor rection s in non-normalize d  va riables). Asthe  ext r a po lat i on isdo n e a t each value inde pendently,  t h is  method do e s n ot  constr ain the  mome n ts of t he finaldistribut io n an d wethereforenormaliz e it as a final st ep. T heresult  o f t his p rocedu r e i s sho wn inFi g. \[f ig:ex tr ap\], wh ich suggests that altho ugh th e fin aldistribut ion isextremely  clo se to a
_\text{th} &\sim_N^{-2/3}\\
\overline{\varepsilon^2}_c &_\sim N^{-4/3} \text{.}\end{aligned}$$

\[eq:fs-law\]

Studying the_large deviation_tails,_we find_that_the “soft” tail_(corresponding to low_energies) has the same_asymptotics as in_the_Tracy-Widom law, but the “hard” tail decays more like a Gaussian and therefore decays_considerably_more slowly_than_the_Tracy-Widom law, $$\begin{aligned}
\log P_\text{empirical} (x)_\sim \begin{cases} x^{\approx 3/2} &_x \ll_0 \\ -x^{\approx 2} & x\gg 0 \end{cases}\end{aligned}$$_with_again the rescaled_variable $x=(\varepsilon-\mu_\varepsilon(N))/\sigma_\varepsilon(N)$. This is shown in Fig. \[fig:fs\_distro\](B). From Fig. \[fig:combo\_pdf\]_(a) it is hard to tell_whether this holds_for_large_$N$ as the shape_of the distribution appears to cross_over from something close to the_Tracy-Widom distribution towards the Gaussian distribution. To_address this question, we construct an_estimate for the converged shape_for very_large $N$ using the following_procedure. We sample_the inverse_function $Q_N(c)=C_N^{-1}(c)$ to_the empirical cumulative distribution functions (CDFs)_$C_N(x)$ for every_$N$ (using the $N$-specific normalization) at_a_number of selected_values_$c=10^{-5},\ldots,1$._At any_given value of_$c$,_this gives_a_set of pairs $(1/N,x_c=Q_N(c))$ which we_use_to extrapolate to $x_c(0)$. We find that_the finite size effects_in_the shape are well_described by $x(N)-x_c(0) \propto N^{-1/3}$_(corresponding to $1/N$ corrections in non-normalized_variables). As_the extrapolation_is done at each value independently, this method does not constrain_the moments of the final distribution_and we therefore normalize_it as_a_final step. The_result_of this_procedure is shown in Fig. \[fig:extrap\], which suggests_that although_the final distribution is extremely close_to a