Near-Extreme Eigenvalues and the First Gap of Hermitian Random Matrices

We study the phenomenon of “crowding” near the largest eigenvalue \(\lambda _\mathrm{max}\) of random \(N \times N\) matrices belonging to the Gaussian Unitary Ensemble of random matrix theory. We focus on two distinct quantities: (i) the density of states (DOS) near \(\lambda _\mathrm{max}\) , \(\rho _\mathrm{DOS}(r,N)\) , which is the average density of eigenvalues located at a distance \(r\) from \(\lambda _\mathrm{max}\) and (ii) the probability density function of the gap between the first two largest eigenvalues, \(p_\mathrm{GAP}(r,N)\) . In the edge scaling limit where \(r = \mathcal{O}(N^{-1/6})\) , which is described by a double scaling limit of a system of unconventional orthogonal polynomials, we show that \(\rho _\mathrm{DOS}(r,N)\) and \(p_\mathrm{GAP}(r,N)\) are characterized by scaling functions which can be expressed in terms of the solution of a Lax pair associated to the Painlevé XXXIV equation. This provides an alternative and simpler expression for the gap distribution, which was recently studied by Witte et al. in Nonlinearity 26:1799, 2013. Our expressions allow to obtain precise asymptotic behaviors of these scaling functions both for small and large arguments.     

Toward a Mathematical Holographic Principle

In work started in [17] and continued in this paper our objective is to study selectors of multivalued functions which have interesting dynamical properties, such as possessing absolutely continuous invariant measures. We specify the graph of a multivalued function by means of lower and upper boundary maps \(\tau _{1}\) and \(\tau _{2}.\) On these boundary maps we define a position dependent random map \(R_{p}=\{\tau _{1},\tau _{2};p,1-p\},\) which, at each time step, moves the point \(x\) to \(\tau _{1}(x)\) with probability \(p(x)\) and to \(\tau _{2}(x)\) with probability \(1-p(x)\) . Under general conditions, for each choice of \(p\) , \(R_{p}\) possesses an absolutely continuous invariant measure with invariant density \(f_{p}.\) Let \(\varvec{\tau }\) be a selector which has invariant density function \(f.\) One of our objectives is to study conditions under which \(p(x)\) exists such that \(R_{p}\) has \(f\) as its invariant density function. When this is the case, the long term statistical dynamical behavior of a selector can be represented by the long term statistical behavior of a random map on the boundaries of \(G.\) We refer to such a result as a mathematical holographic principle. We present examples and study the relationship between the invariant densities attainable by classes of selectors and the random maps based on the boundaries and show that, under certain conditions, the extreme points of the invariant densities for selectors are achieved by bang-bang random maps, that is, random maps for which \(p(x)\in \{0,1\}.\)      

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.     

Continuous and Discrete Painlevé Equations Arising from the Gap Probability Distribution of the Finite n Gaussian Unitary Ensembles

In this paper we study the gap probability problem in the Gaussian unitary ensembles of \(n\) by \(n\) matrices : The probability that the interval \(J := (-a,a)\) is free of eigenvalues. In the works of Tracy and Widom, Adler and Van Moerbeke, and Forrester and Witte on this subject, it has been shown that two Painlevé type differential equations arise in this context. The first is the Jimbo–Miwa–Okomoto \(\sigma \) -form and the second is a particular Painlevé IV. Using the ladder operator technique of orthogonal polynomials we derive three quantities associated with the gap probability, denoted as \(\sigma _n(a)\) , \(R_n(a)\) and \(r_n(a)\) . We show that each one satisfies a second order Painlevé type differential equation as well as a discrete Painlevé type equation. In particular, in addition to providing an elementary derivation of the aforementioned \(\sigma \) -form and Painlevé IV we are able to show that the quantity \(r_n(a)\) satisfies a particular case of Chazy’s second degree second order differential equation. For the discrete equations we show that the quantity \(r_n(a)\) satisfies a particular form of the modified discrete Painlevé II equation obtained by Grammaticos and Ramani in the context of Backlund transformations. We also derive second order second degree difference equations for the quantities \(R_n(a)\) and \(\sigma _n(a)\) .     

Modeling of color-coded III-nitride LED structures with deep quantum wells

We present, to the best of our knowledge, the first successful simulation of color-coded III-nitride light-emitting diodes (LEDs) incorporating in their active regions shallow and deep InGaN quantum wells (QWs). Dichromatic violet–aquamarine semipolar LEDs grown in Ga-polar and N-polar crystallographic orientations (Kawaguchi et al. in Appl Phys Lett 100:231110–231114, 2012) were used as an experimental benchmark. Opposite interface polarization charges in Ga-polar and N-polar LEDs provide different conditions for carrier transport and account for different shape of color-coded emission spectra. To reproduce experimentally observed trends, several effects specific for deep III-nitride QWs were essential in our modeling including \((1)\) strongly non-equilibrium character of active QW populations, \((2)\) dynamic carrier overshoot of narrow QW layers, and \((3)\) Auger-assisted QW depopulation.     

Rigorous Confidence Intervals on Critical Thresholds in 3 Dimensions

We extend the method of Balister, Bollobás and Walters (Phys. Rev. E 76:011110, 2007) for determining rigorous confidence intervals for the critical threshold of two dimensional lattices to three (and higher) dimensional lattices. We describe a method for determining a full confidence interval and apply it to show that the critical threshold for bond percolation on the simple cubic lattice is between \(0.2485\) and \(0.2490\) with \(99.9999\,\%\) confidence, and the critical threshold for site percolation on the same lattice is between \(0.3110\) and \(0.3118\) with \(99.9999\,\%\) confidence.     

Non-Equilibrium Statistical Mechanics of Turbulence

The macroscopic study of hydrodynamic turbulence is equivalent, at an abstract level, to the microscopic study of a heat flow for a suitable mechanical system (Ruelle, PNAS 109:20344–20346, 2012). Turbulent fluctuations (intermittency) then correspond to thermal fluctuations, and this allows to estimate the exponents \(\tau _p\) and \(\zeta _p\) associated with moments of dissipation fluctuations and velocity fluctuations. This approach, initiated in an earlier note (Ruelle, 2012), is pursued here more carefully. In particular we derive probability distributions at finite Reynolds number for the dissipation and velocity fluctuations, and the latter permit an interpretation of numerical experiments (Schumacher, Preprint, 2014). Specifically, if \(p(z)dz\) is the probability distribution of the radial velocity gradient we can explain why, when the Reynolds number \(\mathcal{R}\) increases, \(\ln p(z)\) passes from a concave to a linear then to a convex profile for large \(z\) as observed in (Schumacher, 2014). We show that the central limit theorem applies to the dissipation and velocity distribution functions, so that a logical relation with the lognormal theory of Kolmogorov (J. Fluid Mech. 13:82–85, 1962) and Obukhov is established. We find however that the lognormal behavior of the distribution functions fails at large value of the argument, so that a lognormal theory cannot correctly predict the exponents \(\tau _p\) and \(\zeta _p\) .     

4D Einstein equations in a 5D general Kaluza–Klein space

Based on the new point of view on space–time–matter theory developed in our paper (Bejancu, Gen Rel Grav, 2013), we obtain the $4D$ 4 D Einstein equations in a general $5D$ 5 D Kaluza–Klein space with electromagnetic potentials. In particular, we recover the $4D$ 4 D Einstein equations obtained by Wesson and Ponce de Leon (J Math Phys 33:3883, 1992) in case the electromagnetic potentials vanish identically on $\bar{M}$ M ¯ . The Riemannian horizontal connection and the $4D$ 4 D tensor calculus on $\bar{M}$ M ¯ , are the main tools in the study.     

Fluctuation–Dissipation Relation for Systems with Spatially Varying Friction

When a particle diffuses in a medium with spatially dependent friction coefficient \(\alpha (r)\) at constant temperature \(T\) , it drifts toward the low friction end of the system even in the absence of any real physical force \(f\) . This phenomenon, which has been previously studied in the context of non-inertial Brownian dynamics, is termed “spurious drift”, although the drift is real and stems from an inertial effect taking place at the short temporal scales. Here, we study the diffusion of particles in inhomogeneous media within the framework of the inertial Langevin equation. We demonstrate that the quantity which characterizes the dynamics with non-uniform \(\alpha (r)\) is not the displacement of the particle \(\Delta r=r-r^0\) (where \(r^0\) is the initial position), but rather \(\Delta A(r)=A(r)-A(r^0)\) , where \(A(r)\) is the primitive function of \(\alpha (r)\) . We derive expressions relating the mean and variance of \(\Delta A\) to \(f\) , \(T\) , and the duration of the dynamics \(\Delta t\) . For a constant friction coefficient \(\alpha (r)=\alpha \) , these expressions reduce to the well known forms of the force-drift and fluctuation–dissipation relations. We introduce a very accurate method for Langevin dynamics simulations in systems with spatially varying \(\alpha (r)\) , and use the method to validate the newly derived expressions.     

Microcanonical Analysis of the Curie–Weiss Anisotropic Quantum Heisenberg Model in a Magnetic Field

The anisotropic quantum Heisenberg model with Curie-Weiss-type interactions is studied analytically in several variants of the microcanonical ensemble. (Non)equivalence of microcanonical and canonical ensembles is investigated by studying the concavity properties of entropies. The microcanonical entropy \(s(e,\varvec{m})\) is obtained as a function of the energy \(e\) and the magnetization vector \({\varvec{m}}\) in the thermodynamic limit. Since, for this model, \(e\) is uniquely determined by \({\varvec{m}}\) , the same information can be encoded either in \(s(\varvec{m})\) or \(s(e,m_1,m_2)\) . Although these two entropies correspond to the same physical setting of fixed \(e\) and \({\varvec{m}}\) , their concavity properties differ. The entropy \(s_{{\varvec{h}}}(u)\) , describing the model at fixed total energy \(u\) and in a homogeneous external magnetic field \({\varvec{h}}\) of arbitrary direction, is obtained by reduction from the nonconcave entropy \(s(e,m_1,m_2)\) . In doing so, concavity, and therefore equivalence of ensembles, is restored. \(s_{{\varvec{h}}}(u)\) has nonanalyticities on surfaces of co-dimension 1 in the \((u,\varvec{h})\) -space. Projecting these surfaces into lower-dimensional phase diagrams, we observe that the resulting phase transition lines are situated in the positive-temperature region for some parameter values, and in the negative-temperature region for others. In the canonical setting of a system coupled to a heat bath of positive temperatures, the nonanalyticities in the microcanonical negative-temperature region cannot be observed, and this leads to a situation of effective nonequivalence even when formal equivalence holds.     

Frame-dragging fields and spin 1 gravitomagnetic radiation

Experimental results published in 2004 (Ciufolini and Pavlis in Nature 431:958–960, 2004) and 2011 (Everitt et al. in Phys Rev Lett 106:221101, 1–5, 2011) have confirmed the frame-dragging phenomenon for a spinning earth predicted by Einstein’s field equations. Since this is observed as a precession caused by the gravitomagnetic (GM) field of the rotating body, these experiments may be viewed as measurements of a GM field. The effect is encapsulated in the classic steady state solution for the vector potential field $\zeta $ of a spinning sphere–a solution applying to a sphere with angular momentum J and describing a field filling space for all time (Weinberg in Gravitation and Cosmology, Wiley, New York, 1972). In a laboratory setting one may visualise the case of a sphere at rest $(\zeta =0, \text{ t}<0)$ , being spun up by an external torque at $\text{ t}=0$ to the angular momentum J: the $\zeta $ field of the textbook solution cannot establish itself instantaneously over all space at $\text{ t}=0$ , but must propagate with the velocity c, implying the existence of a travelling GM wave field yielding the textbook $\zeta $ field for large enough t (Tolstoy in Int J Theor Phys 40(5):1021–1031, 2001). The linearized GM field equations of the post-Newtonian approximation being isomorphic with Maxwell’s equations (Braginsky et al. in Phys Rev D 15(6):2047–2060, 1977), such GM waves are dipole waves of spin 1. It is well known that in purely gravitating systems conservation of angular momentum forbids the existence of dipole radiation (Misner et al. in Gravitation, Freeman & Co., New York, 1997); but this rule does not prohibit the insertion of angular momentum into the system from an external source–e.g., by applying a torque to our laboratory sphere.     

Dilaton gravity, (quasi-) black holes,and scalar charge

Electrically charged dust is considered in the framework of Einstein–Maxwell–dilaton gravity with a Lagrangian containing the interaction term \(P(\chi )F_{\mu \nu }F^{\mu \nu }\) , where \(P(\chi )\) is an arbitrary function of the dilaton scalar field \(\chi \) , which can be normal or phantom. Without assumption of spatial symmetry, we show that static configurations exist for arbitrary functions \(g_{00} = \exp (2\gamma (x^{i}))\) ( \(i=1,2,3\) ) and \(\chi =\chi (\gamma )\) . If \(\chi = \mathrm{const}\) , the classical Majumdar–Papapetrou (MP) system is restored. We discuss solutions that represent black holes (BHs) and quasi-black holes (QBHs), deduce some general results and confirm them by examples. In particular, we analyze configurations with spherical and cylindrical symmetries. It turns out that cylindrical BHs and QBHs cannot exist without negative energy density somewhere in space. However, in general, BHs and QBHs can be phantom-free, that is, can exist with everywhere nonnegative energy densities of matter, scalar and electromagnetic fields.     

The Colored Hofstadter Butterfly for the Honeycomb Lattice

We rely on a recent method for determining edge spectra and we use it to compute the Chern numbers for Hofstadter models on the honeycomb lattice having rational magnetic flux per unit cell. Based on the bulk-edge correspondence, the Chern number \(\sigma _\mathrm{H}\) is given as the winding number of an eigenvector of a \(2 \times 2\) transfer matrix, as a function of the quasi-momentum \(k\in (0,2\pi )\) . This method is computationally efficient (of order \(\mathcal {O}(n^4)\) in the resolution of the desired image). It also shows that for the honeycomb lattice the solution for \(\sigma _\mathrm{H}\) for flux \(p/q\) in the \(r\) -th gap conforms with the Diophantine equation \(r=\sigma _\mathrm{H}\cdot p+ s\cdot q\) , which determines \(\sigma _\mathrm{H}\mod q\) . A window such as \(\sigma _\mathrm{H}\in (-q/2,q/2)\) , or possibly shifted, provides a natural further condition for \(\sigma _\mathrm{H}\) , which however turns out not to be met. Based on extensive numerical calculations, we conjecture that the solution conforms with the relaxed condition \(\sigma _\mathrm{H}\in (-q,q)\) .     

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.     

NMR in Magnetic Field of the Earth: Pre-Polarization of Nuclei with Alternating Magnetic Field

The polarization of nuclei in the low static magnetic field \(B_0\) with an alternating magnetic field \(B^{*} (B^{*} \gg B_0)\) at a very low frequency \(f_m\) (but \(f_m\gg 1\) / \({T_1}\) , where \(T_1\) is the spin-lattice relaxation time) has been investigated. The question of the optimization of the energy consumption during the pre-polarization is also considered. The possibilities of the method are illustrated by the observation of nuclear magnetic resonance signals from a few liquids.     

Semiclassical strings in supergravity PFT

In this paper, we introduce the bulk viscosity in the formalism of modified gravity theory in which the gravitational action contains a general function \(f(R,T)\) , where \(R\) and \(T\) denote the curvature scalar and the trace of the energy–momentum tensor, respectively, within the framework of a flat Friedmann–Robertson–Walker model. As an equation of state for a prefect fluid, we take \(p=(\gamma -1)\rho \) , where \(0 \le \gamma \le 2\) and a viscous term as a bulk viscosity due to the isotropic model, of the form \(\zeta =\zeta _{0}+\zeta _{1}H\) , where \(\zeta _{0}\) and \(\zeta _{1}\) are constants, and \(H\) is the Hubble parameter. The exact non-singular solutions to the corresponding field equations are obtained with non-viscous and viscous fluids, respectively, by assuming a simplest particular model of the form of \(f(R,T) = R+2f(T)\) , where \(f(T)=\alpha T\) ( \(\alpha \) is a constant). A big-rip singularity is also observed for \(\gamma <0\) at a finite value of cosmic time under certain constraints. We study all possible scenarios with the possible positive and negative ranges of \(\alpha \) to analyze the expansion history of the universe. It is observed that the universe accelerates or exhibits a transition from a decelerated phase to an accelerated phase under certain constraints of \(\zeta _0\) and \(\zeta _1\) . We compare the viscous models with the non-viscous one through the graph plotted between the scale factor and cosmic time and find that the bulk viscosity plays a major role in the expansion of the universe. A similar graph is plotted for the deceleration parameter with non-viscous and viscous fluids and we find a transition from decelerated to accelerated phase with some form of bulk viscosity.     

DGP cosmological model with generalized Ricci dark energy

The Higgs-boson decay \(h \rightarrow \gamma \ell ^+ \ell ^-\) for various lepton states \(\ell = (e, \, \mu , \, \tau )\) is analyzed. The differential decay width and forward–backward asymmetry are calculated as functions of the dilepton invariant mass in a model where the Higgs boson interacts with leptons and quarks via a mixture of scalar and pseudoscalar couplings. These couplings are partly constrained from data on the decays to leptons, \(h \rightarrow \ell ^+ \ell ^-\) , and quarks \(h \rightarrow q \bar{q} \) (where \(q = (c, \, b)\) ), while the Higgs couplings to the top quark are chosen from the two-photon and two-gluon decay rates. Nonzero values of the forward–backward asymmetry will manifest effects of new physics in the Higgs sector. The decay width and asymmetry integrated over the dilepton invariant mass are also presented.     

Phase Transitions in Layered Systems

We consider the Ising model on \(\mathbb Z\times \mathbb Z\) where on each horizontal line \(\{(x,i), x\in \mathbb Z\}\) , called “layer”, the interaction is given by a ferromagnetic Kac potential with coupling strength \(J_{ \gamma }(x,y)={ \gamma }J({ \gamma }(x-y))\) , where \(J(\cdot )\) is smooth and has compact support; we then add a nearest neighbor ferromagnetic vertical interaction of strength \({ \gamma }^{A}\) , where \(A\ge 2\) is fixed, and prove that for any \(\beta \) larger than the mean field critical value there is a phase transition for all \({ \gamma }\) small enough.     

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 }.\)      

Heterogeneous Memorized Continuous Time Random Walks in an External Force Fields

In this paper, we study the anomalous diffusion of a particle in an external force field whose motion is governed by nonrenewal continuous time random walks with correlated memorized waiting times, which involves Reimann–Liouville fractional derivative or Reimann–Liouville fractional integral. We show that the mean squared displacement of the test particle \(X_{x}\) which is dependent on its location \(x\) of the form (El-Wakil and Zahran, Chaos Solitons Fractals, 12, 1929–1935, 2001) 1 $$\begin{aligned} \langle \mathbb {X}_x^2\rangle (t)=\langle (\Delta X_x(t))^2\rangle _0\sim |x|^{-\theta }t^{\gamma }, \quad 0<\gamma <1, \quad \theta =d_w-2, \end{aligned}$$ where \(d_w>2\) is the anomalous exponent, the diffusion exponent \(\gamma \) is dependent on the model parameters. We obtain the Fokker–Planck-type dynamic equations, and their stationary solutions are of the Boltzmann–Gibbs form. These processes obey a generalized Einstein–Stokes–Smoluchowski relation and the second Einstein relation. We observe that the asymptotic behavior of waiting times and subordinations are of stretched Gaussian distributions. We also discuss the time averaged in the case of an harmonic potential, and show that the process exhibits aging and ergodicity breaking.