Split supersymmetry under GUT and current dark matter constraints

There are four types of two-Higgs doublet models under a discrete \(Z_2\) symmetry imposed to avoid tree-level flavor-changing neutral current, i.e. type-I, type-II, type-X, and type-Y models. We investigate the possibility to discriminate the four models in the light of the flavor physics data, including \(B_s\) – \(\bar{B}_s\) mixing, \(B_{s,d} \rightarrow \mu ^+ \mu ^-\) , \(B\rightarrow \tau \nu \) and \(\bar{B} \rightarrow X_s \gamma \) decays, the recent LHC Higgs data, the direct search for charged Higgs at LEP, and the constraints from perturbative unitarity and vacuum stability. After deriving the combined constraints on the Yukawa interaction parameters, we have shown that the correlation between the mass eigenstate rate asymmetry \(A_{\Delta \Gamma }\) of \(B_{s} \rightarrow \mu ^+ \mu ^-\) and the ratio \(R=\mathcal{B}(B_{s} \rightarrow \mu ^+ \mu ^-)_\mathrm{exp}/ \mathcal{B}(B_{s} \rightarrow \mu ^+ \mu ^-)_\mathrm{SM}\) could be a sensitive probe to discriminate the four models with future precise measurements of the observables in the \(B_{s} \rightarrow \mu ^+ \mu ^-\) decay at LHCb.     

Analyticity of \eta \pi isospin-violating form factors and the \tau \rightarrow \eta \pi \nu second-class decay

We consider the evaluation of the \(\eta \pi \) isospin-violating vector and scalar form factors relying on a systematic application of analyticity and unitarity, combined with chiral expansion results. It is argued that the usual analyticity properties do hold (i.e. no anomalous thresholds are present) in spite of the instability of the \(\eta \) meson in QCD. Unitarity relates the vector form factor to the \(\eta \pi \rightarrow \pi \pi \) amplitude: we exploit progress in formulating and solving the Khuri–Treiman equations for \(\eta \rightarrow 3\pi \) and in experimental measurements of the Dalitz plot parameters to evaluate the shape of the \(\rho \) -meson peak. Observing this peak in the energy distribution of the \(\tau \rightarrow \eta \pi \nu \) decay would be a background-free signature of a second-class amplitude. The scalar form factor is also estimated from a phase dispersive representation using a plausible model for the \(\eta \pi \) elastic scattering \(S\) -wave phase shift and a sum rule constraint in the inelastic region. We indicate how a possibly exotic nature of the \(a_0(980)\) scalar meson manifests itself in a dispersive approach. A remark is finally made on a second-class amplitude in the \(\tau \rightarrow \pi \pi \nu \) decay.     

Pesin’s Entropy Formula for Systems Between C^1 and C^{1+\alpha }

In this article we give a new observation of Pesin’s entropy formula, motivated from Mañé’s proof of (Ergod Theory Dyn Sys 1:95–102, 1981). Let \(M\) be a compact Riemann manifold and \(f:\,M\rightarrow M\) be a \(C^1\) diffeomorphism on \(M\) . If \(\mu \) is an \(f\) -invariant probability measure which is absolutely continuous relative to Lebesgue measure and nonuniformly-H \(\ddot{\text {o}}\) lder-continuous(see Definition 1.1), then we have Pesin’s entropy formula, i.e., the metric entropy \(h_\mu (f)\) satisfies $$\begin{aligned} h_{\mu }(f)=\int \sum _{\lambda _i(x)> 0}\lambda _i(x)d\mu , \end{aligned}$$ where \(\lambda _1(x)\ge \lambda _2(x)\ge \cdots \ge \lambda _{dim\,M}(x)\) are the Lyapunov exponents at \(x\) with respect to \(\mu .\) Nonuniformly-H \(\ddot{\text {o}}\) lder-continuous is a new notion from probabilistic perspective weaker than \(C^{1+\alpha }.\)      

Tracy–Widom at High Temperature

We investigate the marginal distribution of the bottom eigenvalues of the stochastic Airy operator when the inverse temperature \(\beta \) tends to \(0\) . We prove that the minimal eigenvalue, whose fluctuations are governed by the Tracy–Widom \(\beta \) law, converges weakly, when properly centered and scaled, to the Gumbel distribution. More generally we obtain the convergence in law of the marginal distribution of any eigenvalue with given index \(k\) . Those convergences are obtained after a careful analysis of the explosion times process of the Riccati diffusion associated to the stochastic Airy operator. We show that the empirical measure of the explosion times converges weakly to a Poisson point process using estimates proved in Dumaz and Virág (Ann Inst H Poincaré Probab Statist 49(4):915–933, 2013). We further compute the empirical eigenvalue density of the stochastic Airy ensemble on the macroscopic scale when \(\beta \rightarrow 0\) . As an application, we investigate the maximal eigenvalues statistics of \(\beta _N\) -ensembles when the repulsion parameter \(\beta _N\rightarrow 0\) when \(N\rightarrow +\infty \) . We study the double scaling limit \(N\rightarrow +\infty , \beta _N \rightarrow 0\) and argue with heuristic and numerical arguments that the statistics of the marginal distributions can be deduced following the ideas of Edelman and Sutton (J Stat Phys 127(6):1121–1165, 2007) and Ramírez et al. (J Am Math Soc 24:919–944, 2011) from our later study of the stochastic Airy operator.     

Optimal N-Point Configurations on the Sphere: “Magic” Numbers and Smale’s 7th Problem

This paper inquires into the concavity of the map \(N\mapsto v_s(N)\) from the integers \(N\ge 2\) into the minimal average standardized Riesz pair-energies \(v_s(N)\) of \(N\) -point configurations on the sphere \(\mathbb {S}^2\) for various \(s\in \mathbb {R}\) . The standardized Riesz pair-energy of a pair of points on \(\mathbb {S}^2\) a chordal distance \(r\) apart is \(V_s(r)= s^{-1}\left( r^{-s}-1 \right) \) , \(s \ne 0\) , which becomes \(V_0(r) = \ln \frac{1}{r}\) in the limit \(s\rightarrow 0\) . Averaging it over the \(\left( \begin{array}{c} N\\ 2\end{array}\right) \) distinct pairs in a configuration and minimizing over all possible \(N\) -point configurations defines \(v_s(N)\) . It is known that \(N\mapsto v_s(N)\) is strictly increasing for each \(s\in \mathbb {R}\) , and for \(s<2\) also bounded above, thus “overall concave.” It is (easily) proved that \(N\mapsto v_{-2}^{}(N)\) is even locally strictly concave, and that so is the map \(2n\mapsto v_s(2n)\) for \(s<-2\) . By analyzing computer-experimental data of putatively minimal average Riesz pair-energies \(v_s^x(N)\) for \(s\in \{-1,0,1,2,3\}\) and \(N\in \{2,\ldots ,200\}\) , it is found that the map \(N\mapsto {v}_{-1}^x(N)\) is locally strictly concave, while \(N\mapsto {v}_s^x(N)\) is not always locally strictly concave for \(s\in \{0,1,2,3\}\) : concavity defects occur whenever \(N\in {\mathcal {C}}^{x}_+(s)\) (an \(s\) -specific empirical set of integers). It is found that the empirical map \(s\mapsto {\mathcal {C}}^{x}_+(s),\ s\in \{-2,-1,0,1,2,3\}\) , is set-theoretically increasing; moreover, the percentage of odd numbers in \({\mathcal {C}}^{x}_+(s),\ s\in \{0,1,2,3\}\) is found to increase with \(s\) . The integers in \({\mathcal {C}}^{x}_+(0)\) are few and far between, forming a curious sequence of numbers, reminiscent of the “magic numbers” in nuclear physics. It is conjectured that these new “magic numbers” are associated with optimally symmetric optimal-log-energy \(N\) -point configurations on \(\mathbb {S}^2\) . A list of interesting open problems is extracted from the empirical findings, and some rigorous first steps toward their solutions are presented. It is emphasized how concavity can assist in the solution to Smale’s \(7\) th Problem, which asks for an efficient algorithm to find near-optimal \(N\) -point configurations on \(\mathbb {S}^2\) and higher-dimensional spheres.     

Comprehensive Bayesian analysis of rare (semi)leptonic and radiative \varvec{B} decays

The available data on \(|\Delta B| = |\Delta S| = 1\) decays are in good agreement with the Standard Model when permitting subleading power corrections of about \(15\,\%\) at large hadronic recoil. Constraining new-physics effects in \(\mathcal {C}_{7}^{\mathrm {}}\) , \(\mathcal {C}_{9}^{\mathrm {}}\) , \(\mathcal {C}_{10}^{\mathrm {}}\) , the data still demand the same size of power corrections as in the Standard Model. In the presence of chirality-flipped operators, all but one of the power corrections reduce substantially. The Bayes factors are in favor of the Standard Model. Using new lattice inputs for \(B\rightarrow K^*\) form factors and under our minimal prior assumption for the power corrections, the favor shifts toward models with chirality-flipped operators. We use the data to further constrain the hadronic form factors in \(B\rightarrow K\) and \(B\rightarrow K^*\) transitions.     

A flexible scintillation light apparatus for rare event searches

Compelling experimental evidences of neutrino oscillations and their implication that neutrinos are massive particles have given neutrinoless double beta decay ( \(\beta \beta 0\nu \) ) a central role in astroparticle physics. In fact, the discovery of this elusive decay would be a major breakthrough, unveiling that neutrino and antineutrino are the same particle and that the lepton number is not conserved. It would also impact our efforts to establish the absolute neutrino mass scale and, ultimately, understand elementary particle interaction unification. All current experimental programs to search for \(\beta \beta 0\nu \) are facing with the technical and financial challenge of increasing the experimental mass while maintaining incredibly low levels of spurious background. The new concept described in this paper could be the answer which combines all the features of an ideal experiment: energy resolution, low cost mass scalability, isotope choice flexibility and many powerful handles to make the background negligible. The proposed technology is based on the use of arrays of silicon detectors cooled to 120 K to optimize the collection of the scintillation light emitted by ultra-pure crystals. It is shown that with a 54 kg array of natural CaMoO \(_4\) scintillation detectors of this type it is possible to yield a competitive sensitivity on the half-life of the \(\beta \beta 0\nu \) of \(^{100}\) Mo as high as \(\sim \) \(10^{24}\)  years in only 1 year of data taking. The same array made of \(^{40}\) Ca \(^{\mathrm {nat}}\) MoO \(_4\) scintillation detectors (to get rid of the continuous background coming from the two neutrino double beta decay of \(^{48}\) Ca) will instead be capable of achieving the remarkable sensitivity of \(\sim \) \(10^{25}\)  years on the half-life of \(^{100}\) Mo \(\beta \beta 0\nu \) in only 1 year of measurement.     

The Distribution of the Area Under a Bessel Excursion and its Moments

A Bessel excursion is a Bessel process that begins at the origin and first returns there at some given time \(T\) . We study the distribution of the area under such an excursion, which recently found application in the context of laser cooling. The area \(A\) scales with the time as \(A \sim T^{3/2}\) , independent of the dimension, \(d\) , but the functional form of the distribution does depend on \(d\) . We demonstrate that for \(d=1\) , the distribution reduces as expected to the distribution for the area under a Brownian excursion, known as the Airy distribution, deriving a new expression for the Airy distribution in the process. We show that the distribution is symmetric in \(d-2\) , with nonanalytic behavior at \(d=2\) . We calculate the first and second moments of the distribution, as well as a particular fractional moment. We also analyze the analytic continuation from \(d<2\) to \(d>2\) . In the limit where \(d\rightarrow 4\) from below, this analytically continued distribution is described by a one-sided Lévy \(\alpha \) -stable distribution with index \(2/3\) and a scale factor proportional to \([(4-d)T]^{3/2}\) .     

Local Central Limit Theorem for Determinantal Point Processes

We prove a local central limit theorem (LCLT) for the number of points \(N(J)\) in a region \(J\) in \(\mathbb R^d\) specified by a determinantal point process with an Hermitian kernel. The only assumption is that the variance of \(N(J)\) tends to infinity as \(|J| \rightarrow \infty \) . This extends a previous result giving a weaker central limit theorem for these systems. Our result relies on the fact that the Lee–Yang zeros of the generating function for \(\{E(k;J)\}\) —the probabilities of there being exactly \(k\) points in \(J\) —all lie on the negative real \(z\) -axis. In particular, the result applies to the scaled bulk eigenvalue distribution for the Gaussian Unitary Ensemble (GUE) and that of the Ginibre ensemble. For the GUE we can also treat the properly scaled edge eigenvalue distribution. Using identities between gap probabilities, the LCLT can be extended to bulk eigenvalues of the Gaussian Symplectic Ensemble. A LCLT is also established for the probability density function of the \(k\) -th largest eigenvalue at the soft edge, and of the spacing between \(k\) -th neighbors in the bulk.     

The CMSSM and NUHM1 after LHC Run 1

We analyze the impact of data from the full Run 1 of the LHC at 7 and 8 TeV on the CMSSM with \(\mu > 0\) and \(<0\) and the NUHM1 with \(\mu > 0\) , incorporating the constraints imposed by other experiments such as precision electroweak measurements, flavour measurements, the cosmological density of cold dark matter and the direct search for the scattering of dark matter particles in the LUX experiment. We use the following results from the LHC experiments: ATLAS searches for events with \({E\!\!/}_{T}\) accompanied by jets with the full 7 and 8 TeV data, the ATLAS and CMS measurements of the mass of the Higgs boson, the CMS searches for heavy neutral Higgs bosons and a combination of the LHCb and CMS measurements of \(\mathrm{BR}(B_s \rightarrow \mu ^+\mu ^-)\) and \(\mathrm{BR}(B_d \rightarrow \mu ^+\mu ^-)\) . Our results are based on samplings of the parameter spaces of the CMSSM for both \(\mu >0\) and \(\mu <0\) and of the NUHM1 for \(\mu > 0\) with 6.8 \(\times 10^6\) , 6.2 \(\times 10^6\) and 1.6 \(\times 10^7\) points, respectively, obtained using the MultiNest tool. The impact of the Higgs-mass constraint is assessed using FeynHiggs 2.10.0, which provides an improved prediction for the masses of the MSSM Higgs bosons in the region of heavy squark masses. It yields in general larger values of \(M_h\) than previous versions of FeynHiggs, reducing the pressure on the CMSSM and NUHM1. We find that the global \(\chi ^2\) functions for the supersymmetric models vary slowly over most of the parameter spaces allowed by the Higgs-mass and the \({E\!\!/}_{T}\) searches, with best-fit values that are comparable to the \(\chi ^2/\mathrm{dof}\) for the best Standard Model fit. We provide 95 % CL lower limits on the masses of various sparticles and assess the prospects for observing them during Run 2 of the LHC.     

Quantum Measures on Finite Effect Algebras with the Riesz Decomposition Properties

One kind of generalized measures called quantum measures on finite effect algebras, which fulfil the grade-2 additive sum rule, is considered. One basis of vector space of quantum measures on a finite effect algebra with the Riesz decomposition property (RDP for short) is given. It is proved that any diagonally positive symmetric signed measure \(\lambda \) on the tensor product \(E\otimes E\) can determine a quantum measure \(\mu \) on a finite effect algebra \(E\) with the RDP such that \(\mu (x)=\lambda (x\otimes x)\) for any \(x\in E\) . Furthermore, some conditions for a grade-2 additive measure \(\mu \) on a finite effect algebra \(E\) are provided to guarantee that there exists a unique diagonally positive symmetric signed measure \(\lambda \) on \(E\otimes E\) such that \(\mu (x)=\lambda (x\otimes x)\) for any \(x\in E\) . At last, it is showed that any grade- \(t\) quantum measure on a finite effect algebra \(E\) with the RDP is essentially established by the values on a subset of \(E\) .     

Reunion Probabilities of N One-Dimensional Random Walkers with Mixed Boundary Conditions

In this work we extend the results of the reunion probability of \(N\) one-dimensional random walkers to include mixed boundary conditions between their trajectories. The level of the mixture is controlled by a parameter \(c\) , which can be varied from \(c=0\) (independent walkers) to \(c\rightarrow \infty \) (vicious walkers). The expressions are derived by using Quantum Mechanics formalism (QMf) which allows us to map this problem into a Lieb-Liniger gas (LLg) of \(N\) one-dimensional particles. We use Bethe ansatz and Gaudin’s conjecture to obtain the normalized wave-functions and use this information to construct the propagator. As it is well-known, depending on the boundary conditions imposed at the endpoints of a line segment, the statistics of the maximum heights of the reunited trajectories have some connections with different ensembles in Random Matrix Theory. Here we seek to extend those results and consider four models: absorbing, periodic, reflecting, and mixed. In all four cases, the probability that the maximum height is less or equal than \(L\) takes the form \(F_N(L)=A_N\sum _{\varvec{k}\in \Omega _{\text {B}}} \mathrm{e}^{-\sum _{j=1}^Nk_j^2}\mathcal {V}_N(\varvec{k})\) , where \(A_N\) is a normalization constant, \(\mathcal {V}_N(\varvec{k})\) contains a deformed and weighted Vandermonde determinant, and \(\Omega _{\text {B}}\) is the solution set of quasi-momenta \(\varvec{k}\) obeying the Bethe equations for that particular boundary condition.     

Modeling of HOT (111) HgCdTe MWIR detector for fast response operation

The paper reports on photoelectrical performance of the mid-wave infrared (MWIR) (111) HgCdTe high operating temperature detector for the fast response conditions. Detector structure was simulated with software APSYS by Crosslight Inc. The detailed analysis of the time response as a function of device architecture and applied voltage was performed pointing out optimal working conditions. The time response of the MWIR HgCdTe detector with 50 % cut-off wavelength of \(\lambda _{c} \approx 5.3\, \upmu \hbox {m}\) at \(T = 200\)  K was estimated at the level of \(\tau _{s} \approx \) 2,500 ps for \(V = 100\)  mV and series resistance \(R_{Series} = 510\,\Omega \) . The series resistance’s reduction enables to reach \(\tau _{s}\approx 60\!-\!500\)  ps.     

The Higgs boson inclusive decay channels H \rightarrow b\bar{b} and H \rightarrow gg up to four-loop level

The principle of maximum conformality (PMC) has been suggested to eliminate the renormalization scheme and renormalization scale uncertainties, which are unavoidable for the conventional scale setting and are usually important errors for theoretical estimations. In this paper, by applying PMC scale setting, we analyze two important inclusive Standard Model Higgs decay channels, $H\rightarrow b\bar{b}$ and $H\rightarrow gg$ , up to four-loop and three-loop levels, respectively. After PMC scale setting, it is found that the conventional scale uncertainty for these two channels can be eliminated to a high degree. There is small residual initial scale dependence for the Higgs decay widths due to unknown higher-order $\{\beta _i\}$ terms. Up to four-loop level, we obtain $\Gamma (H\rightarrow b\bar{b}) = 2.389\pm 0.073 \pm 0.041$ MeV and up to three-loop level, we obtain $\Gamma (H\rightarrow gg) = 0.373\pm 0.030$ MeV, where the first error is caused by varying $M_H=126\pm 4$ GeV and the second error for $H\rightarrow b\bar{b}$ is caused by varying the $\overline{\mathrm{MS}}$ -running mass $m_b(m_b)=4.18\pm 0.03$ GeV. Taking $H\rightarrow b\bar{b}$ as an example, we present a comparison of three BLM-based scale-setting approaches, e.g. the PMC-I approach based on the PMC–BLM correspondence, the $R_\delta $ -scheme and the seBLM approach, all of which are designed to provide effective ways to identify non-conformal $\{\beta _i\}$ -series at each perturbative order. At four-loop level, all those approaches lead to good pQCD convergence, they have almost the same pQCD series, and their predictions are almost independent on the initial renormalization scale. In this sense, those approaches are equivalent to each other.     

Double inclusive cross sections for gluon production in collision of two projectiles on two targets in the BFKL approach

We compute the normalization of the form factor entering the $B_{s}\rightarrow D_{s}\ell \nu $ decay amplitude by using numerical simulations of QCD on the lattice. From our study with $N_\mathrm{f}=2$ dynamical light quarks, and by employing the maximally twisted Wilson quark action, we obtain in the continuum limit ${\mathcal {G}}(1)= 1.052(46)$ . We also compute the scalar and tensor form factors in the region near zero recoil and find $f_0(q_0^2)/f_+(q_0^2)=0.77(2)$ , $f_T(q_0^2,m_b)/f_+(q_0^2)=1.08(7)$ , for $q_0^2=11.5\ \mathrm{GeV}^2$ . The latter results are useful for searching the effects of physics beyond the Standard Model in $B_{s}\rightarrow D_{s}\ell \nu $ decays. Our results for the similar form factors relevant to the non-strange case indicate that the method employed here can be used to achieve the precision determination of the $B\rightarrow D\ell \nu $ decay amplitude as well.     

Refractive index changes and nuclear activation upon irradiation of lithium niobate crystals with different low-mass,high-energy ions

Refractive index changes \(\Delta n\) in lithium niobate crystals upon irradiation with high-energy protons, deuterons, \(^3\) He, and \(^4\alpha \) particles (up to 14 MeV/nucleon) are created, and the accompanying, unwanted nuclear activation is investigated. The measurements give answers to the question which ion is the best choice depending on the requirements: largest values of \(\Delta n\) are achieved with \(^4\alpha \) particles, low nuclear activation with deuterons, or the best tradeoff between \(\Delta n\) and activation with \(^3\) He, respectively.     

QCD analysis of forward neutron production in DIS

Observing light-by-light scattering at the large hadron collider (LHC) has received quite some attention and it is believed to be a clean and sensitive channel to possible new physics. In this paper, we study the diphoton production at the LHC via the process \({{pp}}\rightarrow {{p}}\gamma \gamma {{p}}\rightarrow {{p}}\gamma \gamma {{p}}\) through graviton exchange in the large extra dimension (LED) model. Typically, when we do the background analysis, we also study the double Pomeron exchange of \(\gamma \gamma \) production. We compare its production in the quark–quark collision mode to the gluon–gluon collision mode and find that contributions from the gluon–gluon collision mode are comparable to the quark–quark one. Our result shows, for extra dimension \(\delta =4\) , with an integrated luminosity \(\mathcal{L} = 200\,\mathrm{fb}^{-1}\) at the 14 TeV LHC, that diphoton production through graviton exchange can probe the LED effects up to the scale \({M}_{S}=5.06 (4.51, 5.11)\,\mathrm{TeV}\) for the forward detector acceptance \(\xi _1 (\xi _2, \xi _3)\) , respectively, where \(0.0015<\xi _1<0.5\) , \(0.1<\xi _2<0.5\) , and \(0.0015<\xi _3<0.15\) .     

Antiproton low-energy collisions with Ps-atoms and true muonium atoms (μ + μ −)

Three-charge-particle collisions with participation of ultra-slow antiprotons ( \(\overline {\rm {p}}\) ) is the subject of this work. Specifically we compute the total cross sections and corresponding thermal rates of the following three-body reactions: \(\overline {\rm p}+(e^+e^-) \rightarrow \overline {\rm {H}} + e^-\) and \(\overline {\rm p}+(\mu ^+\mu ^-) \rightarrow \overline {\rm {H}}_{\mu } + \mu ^-\) , where \(e^-(\mu ^-)\) is an electron (muon) and \(e^+(\mu ^+)\) is a positron (antimuon) respectively, \(\overline {\rm {H}}=(\overline {\rm p}e^+)\) is an antihydrogen atom and \(\overline {\rm {H}}_{\mu }=(\overline {\rm p}\mu ^+)\) is a muonic antihydrogen atom, i.e. a bound state of \(\overline {\rm {p}}\) and μ ⁺. A set of two-coupled few-body Faddeev-Hahn-type (FH-type) equations is numerically solved in the framework of a modified close-coupling expansion approach.     

Maximal Distance Travelled by N Vicious Walkers Till Their Survival

We consider N Brownian particles moving on a line starting from initial positions \(\mathbf{{u}}\equiv \{u_1,u_2,\ldots u_N\}\) such that \(0 . Their motion gets stopped at time \(t_s\) when either two of them collide or when the particle closest to the origin hits the origin for the first time. For \(N=2\) , we study the probability distribution function \(p_1(m|\mathbf{{u}})\) and \(p_2(m|\mathbf{{u}})\) of the maximal distance travelled by the \(1^{\text {st}}\) and \(2^{\text {nd}}\) walker till \(t_s\) . For general N particles with identical diffusion constants \(D\) , we show that the probability distribution \(p_N(m|\mathbf{u})\) of the global maximum \(m_N\) , has a power law tail \(p_i(m|\mathbf{{u}}) \sim {N^2B_N\mathcal {F}_{N}(\mathbf{u})}/{m^{\nu _N}}\) with exponent \(\nu _N =N^2+1\) . We obtain explicit expressions of the function \(\mathcal {F}_{N}(\mathbf{u})\) and of the N dependent amplitude \(B_N\) which we also analyze for large N using techniques from random matrix theory. We verify our analytical results through direct numerical simulations.