BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Pierpaolo Vivo

Pierpaolo Vivo appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
cond-mat.stat-mechmath-phmath.MPcond-mat.mes-hallquant-phcond-mat.dis-nnmath.PRmath.STstat.OTStatistics TheoryBiological PhysicsInformation Theorymath.ITphysics.data-anphysics.soc-phQuantitative MethodsComputational ComplexityComputer Science and Game Theorycond-mat.quant-gascond-mat.softhep-thmath.COq-fin.GNq-fin.STSocial and Information Networks

Catalog footprint

What is connected

40works
25topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

40 published item(s)

preprint2022arXiv

Exact and Approximate Mean First Passage Times on Trees and other Necklace Structures: a Local Equilibrium Approach

In this work we propose a novel method to calculate mean first-passage times (MFPTs) for random walks on graphs, based on a dimensionality reduction technique for Markov State Models, known as local-equilibrium (LE). We show that for a broad class of graphs, which includes trees, LE coarse-graining preserves the MFPTs between certain nodes, upon making a suitable choice of the coarse-grained states (or clusters). We prove that this relation is exact for graphs that can be coarse-grained into a one-dimensional lattice where each cluster connects to the lattice only through a single node of the original graph. A side result of the proof generalises the well-known essential edge lemma (EEL), which is valid for reversible random walks, to irreversible walkers. Such a generalised EEL leads to explicit formulae for the MFPTs between certain nodes in this class of graphs. For graphs that do not fall in this class, the generalised EEL provides useful approximations if the graph allows a one-dimensional coarse-grained representation and the clusters are sparsely interconnected. We first demonstrate our method for the simple random walk on the $c$-ary tree, then we consider other graph structures and more general random walks, including irreversible random walks.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Graph-combinatorial approach for large deviations of Markov chains

We consider discrete-time Markov chains and study large deviations of the pair empirical occupation measure, which is useful to compute fluctuations of pure-additive and jump-type observables. We provide an exact expression for the finite-time moment generating function, which is split in cycles and paths contributions, and scaled cumulant generating function of the pair empirical occupation measure via a graph-combinatorial approach. The expression obtained allows us to give a physical interpretation of interaction and entropic terms, and of the Lagrange multipliers, and may serve as a starting point for sub-leading asymptotics. We illustrate the use of the method for a simple two-state Markov chain.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Instabilities of complex fluids with partially structured and partially random interactions

We develop a theory for thermodynamic instabilities of complex fluids composed of many interacting chemical species organised in families. This model includes partially structured and partially random interactions and can be solved exactly using tools from random matrix theory. The model exhibits three kinds of fluid instabilities: one in which the species form a condensate with a local density that depends on their family (family condensation); one in which species demix in two phases depending on their family (family demixing); and one in which species demix in a random manner irrespective of their family (random demixing). We determine the critical spinodal density of the three types of instabilities and find that the critical spinodal density is finite for both family condensation and family demixing, while for random demixing the critical spinodal density grows as the square root of the number of species. We use the developed framework to describe phase-separation instability of the cytoplasm induced by a change in pH.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

When facts fail: Bias, polarisation and truth in social networks

Online social networks provide users with unprecedented opportunities to engage with diverse opinions. At the same time, they enable confirmation bias on large scales by empowering individuals to self-select narratives they want to be exposed to. A precise understanding of such tradeoffs is still largely missing. We introduce a social learning model where most participants in a network update their beliefs unbiasedly based on new information, while a minority of participants reject information that is incongruent with their preexisting beliefs. This simple mechanism generates permanent opinion polarization and cascade dynamics, and accounts for the aforementioned tradeoff between confirmation bias and social connectivity through analytic results. We investigate the model's predictions empirically using US county-level data on the impact of Internet access on the formation of beliefs about global warming. We conclude by discussing policy implications of our model, highlighting the downsides of debunking and suggesting alternative strategies to contrast misinformation.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A shortcut through the Coulomb gas method for spectral linear statistics on random matrices

In the last decade, spectral linear statistics on large dimensional random matrices have attracted significant attention. Within the physics community, a privileged role has been played by invariant matrix ensembles for which a two dimensional Coulomb gas analogy is available. We present a critical revision of the Coulomb gas method in Random Matrix Theory (RMT) borrowing language and tools from Large Deviations Theory. This allows us to formalize an equivalent, but more effective and quicker route toward RMT free energy calculations. Moreover, we argue that this more modern viewpoint is likely to shed further light on the interesting issues of weak phase transitions and evaporation phenomena recently observed in RMT.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Correlators for the Wigner-Smith time-delay matrix of chaotic cavities

We study the Wigner-Smith time-delay matrix $Q$ of a ballistic quantum dot supporting $N$ scattering channels. We compute the $v$-point correlators of the power traces $\mathrm{Tr} Q^κ$ for arbitrary $v\geq1$ at leading order for large $N$ using techniques from the random matrix theory approach to quantum chromodynamics. We conjecture that the cumulants of the $\mathrm{Tr} Q^κ$'s are integer-valued at leading order in $N$ and include a MATHEMATICA code that computes their generating functions recursively.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

From non-ergodic eigenvectors to local resolvent statistics and back: a random matrix perspective

We study the statistics of the local resolvent and non-ergodic properties of eigenvectors for a generalised Rosenzweig-Porter $N\times N$ random matrix model, undergoing two transitions separated by a delocalised non-ergodic phase. Interpreting the model as the combination of on-site random energies $\{a_i\}$ and a structurally disordered hopping, we found that each eigenstate is delocalised over $N^{2-γ}$ sites close in energy $|a_j-a_i|\leq N^{1-γ}$ in agreement with Kravtsov \emph{et al}, arXiv:1508.01714. Our other main result, obtained combining a recurrence relation for the resolvent matrix with insights from Dyson's Brownian motion, is to show that the properties of the non-ergodic delocalised phase can be probed studying the statistics of the local resolvent in a non-standard scaling limit.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Large deviations of radial statistics in the two-dimensional one-component plasma

The two-dimensional one-component plasma is an ubiquitous model for several vortex systems. For special values of the coupling constant $βq^2$ (where $q$ is the particles charge and $β$ the inverse temperature), the model also corresponds to the eigenvalues distribution of normal matrix models. Several features of the system are discussed in the limit of large number $N$ of particles for generic values of the coupling constant. We show that the statistics of a class of radial observables produces a rich phase diagram, and their asymptotic behaviour in terms of large deviation functions is calculated explicitly, including next-to-leading terms up to order 1/N. We demonstrate a split-off phenomenon associated to atypical fluctuations of the edge density profile. We also show explicitly that a failure of the fluid phase assumption of the plasma can break a genuine $1/N$-expansion of the free energy. Our findings are corroborated by numerical comparisons with exact finite-N formulae valid for $βq^2=2$.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Large deviations of spread measures for Gaussian matrices

For a large $n\times m$ Gaussian matrix, we compute the joint statistics, including large deviation tails, of generalized and total variance - the scaled log-determinant $H$ and trace $T$ of the corresponding $n\times n$ covariance matrix. Using a Coulomb gas technique, we find that the Laplace transform of their joint distribution $\mathcal{P}_n(h,t)$ decays for large $n,m$ (with $c=m/n\geq 1$ fixed) as $\hat{\mathcal{P}}_n(s,w)\approx \exp\left(-βn^2 J(s,w)\right)$, where $β$ is the Dyson index of the ensemble and $J(s,w)$ is a $β$-independent large deviation function, which we compute exactly for any $c$. The corresponding large deviation functions in real space are worked out and checked with extensive numerical simulations. The results are complemented with a finite $n,m$ treatment based on the Laguerre-Selberg integral. The statistics of atypically small log-determinants is shown to be driven by the split-off of the smallest eigenvalue, leading to an abrupt change in the large deviation speed.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Large-$N$ expansion for the time-delay matrix of ballistic chaotic cavities

We consider the $1/N$-expansion of the moments of the proper delay times for a ballistic chaotic cavity supporting $N$ scattering channels. In the random matrix approach, these moments correspond to traces of negative powers of Wishart matrices. For systems with and without broken time reversal symmetry (Dyson indices $β=1$ and $β=2$) we obtain a recursion relation, which efficiently generates the coefficients of the $1/N$-expansion of the moments. The integrality of these coefficients and their possible diagrammatic interpretation is discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Number statistics for $β$-ensembles of random matrices: applications to trapped fermions at zero temperature

Let $\mathcal{P}_β^{(V)} (N_{\cal I})$ be the probability that a $N\times N$ $β$-ensemble of random matrices with confining potential $V(x)$ has $N_{\cal I}$ eigenvalues inside an interval ${\cal I}=[a,b]$ of the real line. We introduce a general formalism, based on the Coulomb gas technique and the resolvent method, to compute analytically $\mathcal{P}_β^{(V)} (N_{\cal I})$ for large $N$. We show that this probability scales for large $N$ as $\mathcal{P}_β^{(V)} (N_{\cal I})\approx \exp\left(-βN^2 ψ^{(V)}(N_{\cal I} /N)\right)$, where $β$ is the Dyson index of the ensemble. The rate function $ψ^{(V)}(k_{\cal I})$, independent of $β$, is computed in terms of single integrals that can be easily evaluated numerically. The general formalism is then applied to the classical $β$-Gaussian (${\cal I}=[-L,L]$), $β$-Wishart (${\cal I}=[1,L]$) and $β$-Cauchy (${\cal I}=[-L,L]$) ensembles. Expanding the rate function around its minimum, we find that generically the number variance ${\rm Var}(N_{\cal I})$ exhibits a non-monotonic behavior as a function of the size of the interval, with a maximum that can be precisely characterized. These analytical results, corroborated by numerical simulations, provide the full counting statistics of many systems where random matrix models apply. In particular, we present results for the full counting statistics of zero temperature one-dimensional spinless fermions in a harmonic trap.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Random pure states: quantifying bipartite entanglement beyond the linear statistics

We analyze the properties of entangled random pure states of a quantum system partitioned into two smaller subsystems of dimensions $N$ and $M$. Framing the problem in terms of random matrices with a fixed-trace constraint, we establish, for arbitrary $N \leq M$, a general relation between the $n$-point densities and the cross-moments of the eigenvalues of the reduced density matrix, i.e. the so-called Schmidt eigenvalues, and the analogous functionals of the eigenvalues of the Wishart-Laguerre ensemble of the random matrix theory. This allows us to derive explicit expressions for two-level densities, and also an exact expression for the variance of von Neumann entropy at finite $N,M$. Then we focus on the moments $\mathbb{E}\{K^a\}$ of the Schmidt number $K$, the reciprocal of the purity. This is a random variable supported on $[1,N]$, which quantifies the number of degrees of freedom effectively contributing to the entanglement. We derive a wealth of analytical results for $\mathbb{E}\{K^a\}$ for $N = 2$ and $N=3$ and arbitrary $M$, and also for square $N = M$ systems by spotting for the latter a connection with the probability $P(x_{min}^{GUE} \geq \sqrt{2N}ξ)$ that the smallest eigenvalue $x_{min}^{GUE}$ of a $N\times N$ matrix belonging to the Gaussian Unitary Ensemble is larger than $\sqrt{2N}ξ$. As a byproduct, we present an exact asymptotic expansion for $P(x_{min}^{GUE} \geq \sqrt{2N}ξ)$ for finite $N$ as $ξ\to \infty$. Our results are corroborated by numerical simulations whenever possible, with excellent agreement.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Universal covariance formula for linear statistics on random matrices

We derive an analytical formula for the covariance $\mathrm{Cov}(A,B)$ of two smooth linear statistics $A=\sum_i a(λ_i)$ and $B=\sum_i b(λ_i)$ to leading order for $N\to\infty$, where $\{λ_i\}$ are the $N$ real eigenvalues of a general one-cut random-matrix model with Dyson index $β$. The formula, carrying the universal $1/β$ prefactor, depends on the random-matrix ensemble only through the edge points $[λ_-,λ_+]$ of the limiting spectral density. For $A=B$, we recover in some special cases the classical variance formulas by Beenakker and Dyson-Mehta, clarifying the respective ranges of applicability. Some choices of $a(x)$ and $b(x)$ lead to a striking \emph{decorrelation} of the corresponding linear statistics. We provide two applications - the joint statistics of conductance and shot noise in ideal chaotic cavities, and some new fluctuation relations for traces of powers of random matrices.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A unified fluctuation formula for one-cut $β$-ensembles of random matrices

Using a Coulomb gas approach, we compute the generating function of the covariances of power traces for one-cut $β$-ensembles of random matrices in the limit of large matrix size. This formula depends only on the support of the spectral density, and is therefore universal for a large class of models. This allows us to derive a closed-form expression for the limiting covariances of an arbitrary one-cut $β$-ensemble. As particular cases of the main result we consider the classical $β$-Gaussian, $β$-Wishart and $β$-Jacobi ensembles, for which we derive previously available results as well as new ones within a unified simple framework. We also discuss the connections between the problem of trace fluctuations for the Gaussian Unitary Ensemble and the enumeration of planar maps.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Invariant sums of random matrices and the onset of level repulsion

We compute analytically the joint probability density of eigenvalues and the level spacing statistics for an ensemble of random matrices with interesting features. It is invariant under the standard symmetry groups (orthogonal and unitary) and yet the interaction between eigenvalues is not Vandermondian. The ensemble contains real symmetric or complex hermitian matrices $\mathbf{S}$ of the form $\mathbf{S}=\sum_{i=1}^M \langle \mathbf{O}_i \mathbf{D}_i\mathbf{O}_i^{\mathrm{T}}\rangle$ or $\mathbf{S}=\sum_{i=1}^M \langle \mathbf{U}_i \mathbf{D}_i\mathbf{U}_i^\dagger\rangle$ respectively. The diagonal matrices $\mathbf{D}_i=\mathrm{diag}\{λ_1^{(i)},\ldots,λ_N^{(i)}\}$ are constructed from real eigenvalues drawn \emph{independently} from distributions $p^{(i)}(x)$, while the matrices $\mathbf{O}_i$ and $\mathbf{U}_i$ are all orthogonal or unitary. The average $\langle\cdot\rangle$ is simultaneously performed over the symmetry group and the joint distribution of $\{λ_j^{(i)}\}$. We focus on the limits i.) $N\to\infty$ and ii.) $M\to\infty$, with $N=2$. In the limit i.), the resulting sum $\mathbf{S}$ develops level repulsion even though the original matrices do not feature it, and classical RMT universality is restored asymptotically. In the limit ii.) the spacing distribution attains scaling forms that are computed exactly: for the orthogonal case, we recover the $β=1$ Wigner's surmise, while for the unitary case an entirely new universal distribution is obtained. Our results allow to probe analytically the microscopic statistics of the sum of random matrices that become asymptotically free. We also give an interpretation of this model in terms of radial random walks in a matrix space. The analytical results are corroborated by numerical simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Joint statistics of quantum transport in chaotic cavities

We study the joint statistics of conductance $G$ and shot noise $P$ in chaotic cavities supporting a large number $N$ of open electronic channels in the two attached leads. We determine the full phase diagram in the $(G,P)$ plane, employing a Coulomb gas technique on the joint density of transmission eigenvalues, as dictated by Random Matrix Theory. We find that in the region of typical fluctuations, conductance and shot noise are uncorrelated and jointly Gaussian, and away from it they fluctuate according to a different joint rate function in each phase of the $(G,P)$ plane. Different functional forms of the rate function in different regions emerge as a direct consequence of third order phase transitions in the associated Coulomb gas problem.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Large deviations of the maximum of independent and identically distributed random variables

A pedagogical account of some aspects of Extreme Value Statistics (EVS) is presented from the somewhat non-standard viewpoint of Large Deviation Theory. We address the following problem: given a set of $N$ i.i.d. random variables $\{X_1,\ldots,X_N\}$ drawn from a parent probability density function (pdf) $p(x)$, what is the probability that the maximum value of the set $X_{\mathrm{max}}=\max_i X_i$ is "atypically larger" than expected? The cases of exponential and Gaussian distributed variables are worked out in detail, and the right rate function for a general pdf in the Gumbel basin of attraction is derived. The Gaussian case convincingly demonstrates that the full rate function cannot be determined from the knowledge of the limiting distribution (Gumbel) alone, thus implying that it indeed carries additional information. Given the simplicity and richness of the result and its derivation, its absence from textbooks, tutorials and lecture notes on EVS for physicists appears inexplicable.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A singular-potential random matrix model arising in mean-field glassy systems

We consider an invariant random matrix model where the standard Gaussian potential is distorted by an additional single pole of order $m$. We compute the average or macroscopic spectral density in the limit of large matrix size, solving the loop equation with the additional constraint of vanishing trace on average. The density is generally supported on two disconnected intervals lying on the two sides of the pole. In the limit of having no pole, we recover the standard semicircle. Obtained in the planar limit, our results apply to matrices with orthogonal, unitary or symplectic symmetry alike. The orthogonal case with $m=2$ is motivated by an application to spin glass physics. In the Sherrington-Kirkpatrick mean-field model, in the paramagnetic phase and for sufficiently large systems the spin glass susceptibility is a random variable, depending on the realization of disorder. It is essentially given by a linear statistics on the eigenvalues of the coupling matrix. As such its large deviation function can be computed using standard Coulomb fluid techniques. The resulting free energy of the associated fluid precisely corresponds to the partition function of our random matrix model. Numerical simulations provide an excellent confirmation of our analytical results.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Phase transitions and edge scaling of number variance in Gaussian random matrices

We consider $N\times N$ Gaussian random matrices, whose average density of eigenvalues has the Wigner semi-circle form over $[-\sqrt{2},\sqrt{2}]$. For such matrices, using a Coulomb gas technique, we compute the large $N$ behavior of the probability $\mathcal{P}_{\scriptscriptstyle N,L}(N_L)$ that $N_L$ eigenvalues lie within the box $[-L,L]$. This probability scales as $\mathcal{P}_{\scriptscriptstyle N,L}(N_L=κ_L N)\approx\exp\left(-β N^2 ψ_L(κ_L)\right)$, where $β$ is the Dyson index of the ensemble and $ψ_L(κ_L)$ is a $β$-independent rate function that we compute exactly. We identify three regimes as $L$ is varied: (i) $\, N^{-1}\ll L<\sqrt{2}$ (bulk), (ii) $\ L\sim\sqrt{2}$ on a scale of $\mathcal{O}(N^{-{2}/{3}})$ (edge) and (iii) $\ L > \sqrt{2}$ (tail). We find a dramatic non-monotonic behavior of the number variance $V_N(L)$ as a function of $L$: after a logarithmic growth $\propto \ln (N L)$ in the bulk (when $L \sim {\cal O}(1/N)$), $V_N(L)$ decreases abruptly as $L$ approaches the edge of the semi-circle before it decays as a stretched exponential for $L > \sqrt{2}$. This "drop-off" of $V_N(L)$ at the edge is described by a scaling function $\tilde V_β$ which smoothly interpolates between the bulk (i) and the tail (iii). For $β= 2$ we compute $\tilde V_2$ explicitly in terms of the Airy kernel. These analytical results, verified by numerical simulations, directly provide for $β=2$ the full statistics of particle-number fluctuations at zero temperature of 1d spinless fermions in a harmonic trap.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Phase transitions in the condition number distribution of Gaussian random matrices

We study the statistics of the condition number $κ=λ_{\mathrm{max}}/λ_{\mathrm{min}}$ (the ratio between largest and smallest squared singular values) of $N\times M$ Gaussian random matrices. Using a Coulomb fluid technique, we derive analytically and for large $N$ the cumulative $\mathcal{P}[κ<x]$ and tail-cumulative $\mathcal{P}[κ>x]$ distributions of $κ$. We find that these distributions decay as $\mathcal{P}[κ<x]\approx\exp\left(-βN^2 Φ_{-}(x)\right)$ and $\mathcal{P}[κ>x]\approx\exp\left(-βN Φ_{+}(x)\right)$, where $β$ is the Dyson index of the ensemble. The left and right rate functions $Φ_{\pm}(x)$ are independent of $β$ and calculated exactly for any choice of the rectangularity parameter $α=M/N-1>0$. Interestingly, they show a weak non-analytic behavior at their minimum $\langleκ\rangle$ (corresponding to the average condition number), a direct consequence of a phase transition in the associated Coulomb fluid problem. Matching the behavior of the rate functions around $\langleκ\rangle$, we determine exactly the scale of typical fluctuations $\sim\mathcal{O}(N^{-2/3})$ and the tails of the limiting distribution of $κ$. The analytical results are in excellent agreement with numerical simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Entanglement production in non-ideal cavities and optimal opacity

We compute analytically the distributions of concurrence $\bm{\mathcal{C}}$ and squared norm $\bm{\mathcal{N}}$ for the production of electronic entanglement in a chaotic quantum dot. The dot is connected to the external world via one ideal and one partially transparent lead, characterized by the opacity $γ$. The average concurrence increases with $γ$ while the average squared norm of the entangled state decreases, making it less likely to be detected. When a minimal detectable norm $\bm{\mathcal{N}}_0$ is required, the average concurrence is maximal for an optimal value of the opacity $γ^\star(\bm{\mathcal{N}}_0)$ which is explicitly computed as a function of $\bm{\mathcal{N}}_0$. If $\bm{\mathcal{N}}_0$ is larger than the critical value $\bm{\mathcal{N}}_0^\star\simeq 0.3693\dots$, the average entanglement production is maximal for the completely ideal case, a direct consequence of an interesting bifurcation effect.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Index Distribution of Cauchy Random Matrices

Using a Coulomb gas technique, we compute analytically the probability $\mathcal{P}_β^{(C)}(N_+,N)$ that a large $N\times N$ Cauchy random matrix has $N_+$ positive eigenvalues, where $N_+$ is called the index of the ensemble. We show that this probability scales for large $N$ as $\mathcal{P}_β^{(C)}(N_+,N)\approx \exp\left[-βN^2 ψ_C(N_+/N)\right]$, where $β$ is the Dyson index of the ensemble. The rate function $ψ_C(κ)$ is computed in terms of single integrals that are easily evaluated numerically and amenable to an asymptotic analysis. We find that the rate function, around its minimum at $κ=1/2$, has a quadratic behavior modulated by a logarithmic singularity. As a consequence, the variance of the index scales for large $N$ as $\mathrm{Var}(N_+)\sim σ_C\ln N$, where $σ_C=2/(βπ^2)$ is twice as large as the corresponding prefactor in the Gaussian and Wishart cases. The analytical results are checked by numerical simulations and against an exact finite $N$ formula which, for $β=2$, can be derived using orthogonal polynomials.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Outage Capacity for the Optical MIMO Channel

MIMO processing techniques in fiber optical communications have been proposed as a promising approach to meet increasing demand for information throughput. In this context, the multiple channels correspond to the multiple modes and/or multiple cores in the fiber. In this paper we characterize the distribution of the mutual information with Gaussian input in a simple channel model for this system. Assuming significant cross talk between cores, negligible backscattering and near-lossless propagation in the fiber, we model the transmission channel as a random complex unitary matrix. The loss in the transmission may be parameterized by a number of unutilized channels in the fiber. We analyze the system in a dual fashion. First, we evaluate a closed-form expression for the outage probability, which is handy for small matrices. We also apply the asymptotic approach, in particular the Coulomb gas method from statistical mechanics, to obtain closed-form results for the ergodic mutual information, its variance as well as the outage probability for Gaussian input in the limit of large number of cores/modes. By comparing our analytic results to simulations, we see that, despite the fact that this method is nominally valid for large number of modes, our method is quite accurate even for small to modest number of channels.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Statistics of quantum transport in weakly non-ideal chaotic cavities

We consider statistics of electronic transport in chaotic cavities where time-reversal symmetry is broken and one of the leads is weakly non-ideal, i.e. it contains tunnel barriers characterized by tunneling probabilities $Γ_i$. Using symmetric function expansions and a generalized Selberg integral, we develop a systematic perturbation theory in $1-Γ_i$ valid for arbitrary number of channels, and obtain explicit formulas up to second order for the average and variance of the conductance, and for the average shot-noise. Higher moments of the conductance are considered to leading order.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Equilibrium strategy and population-size effects in lowest unique bid auctions

In lowest unique bid auctions, $N$ players bid for an item. The winner is whoever places the \emph{lowest} bid, provided that it is also unique. We use a grand canonical approach to derive an analytical expression for the equilibrium distribution of strategies. We then study the properties of the solution as a function of the mean number of players, and compare them with a large dataset of internet auctions. The theory agrees with the data with striking accuracy for small population size $N$, while for larger $N$ a qualitatively different distribution is observed. We interpret this result as the emergence of two different regimes, one in which adaptation is feasible and one in which it is not. Our results question the actual possibility of a large population to adapt and find the optimal strategy when participating in a collective game.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Invariant $β$-Wishart ensembles, crossover densities and asymptotic corrections to the Marchenko-Pastur law

We construct a diffusive matrix model for the $β$-Wishart (or Laguerre) ensemble for general continuous $β\in [0,2]$, which preserves invariance under the orthogonal/unitary group transformation. Scaling the Dyson index $β$ with the largest size $M$ of the data matrix as $β=2c/M$ (with $c$ a fixed positive constant), we obtain a family of spectral densities parametrized by $c$. As $c$ is varied, this density interpolates continuously between the Mar\vcenko-Pastur ($c\to \infty$ limit) and the Gamma law ($c\to 0$ limit). Analyzing the full Stieltjes transform (resolvent) equation, we obtain as a byproduct the correction to the Mar\vcenko-Pastur density in the bulk up to order 1/M for all $β$ and up to order $1/M^2$ for the particular cases $β=1,2$.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Large deviations of the top eigenvalue of large Cauchy random matrices

We compute analytically the probability density function (pdf) of the largest eigenvalue $λ_{\max}$ in rotationally invariant Cauchy ensembles of $N\times N$ matrices. We consider unitary ($β= 2$), orthogonal ($β=1$) and symplectic ($β=4$) ensembles of such heavy-tailed random matrices. We show that a central non-Gaussian regime for $λ_{\max} \sim \mathcal{O}(N)$ is flanked by large deviation tails on both sides which we compute here exactly for any value of $β$. By matching these tails with the central regime, we obtain the exact leading asymptotic behaviors of the pdf in the central regime, which generalizes the Tracy-Widom distribution known for Gaussian ensembles, both at small and large arguments and for any $β$. Our analytical results are confirmed by numerical simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Compact smallest eigenvalue expressions in Wishart-Laguerre ensembles with or without fixed-trace

The degree of entanglement of random pure states in bipartite quantum systems can be estimated from the distribution of the extreme Schmidt eigenvalues. For a bipartition of size M\geq N, these are distributed according to a Wishart-Laguerre ensemble (WL) of random matrices of size N x M, with a fixed-trace constraint. We first compute the distribution and moments of the smallest eigenvalue in the fixed trace orthogonal WL ensemble for arbitrary M\geq N. Our method is based on a Laplace inversion of the recursive results for the corresponding orthogonal WL ensemble by Edelman. Explicit examples are given for fixed N and M, generalizing and simplifying earlier results. In the microscopic large-N limit with M-N fixed, the orthogonal and unitary WL distributions exhibit universality after a suitable rescaling and are therefore independent of the constraint. We prove that very recent results given in terms of hypergeometric functions of matrix argument are equivalent to more explicit expressions in terms of a Pfaffian or determinant of Bessel functions. While the latter were mostly known from the random matrix literature on the QCD Dirac operator spectrum, we also derive some new results in the orthogonal symmetry class.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Moments of Wishart-Laguerre and Jacobi ensembles of random matrices: application to the quantum transport problem in chaotic cavities

We collect explicit and user-friendly expressions for one-point densities of the real eigenvalues $\{λ_i\}$ of $N\times N$ Wishart-Laguerre and Jacobi random matrices with orthogonal, unitary and symplectic symmetry. Using these formulae, we compute integer moments $τ_n=<\sum_{i=1}^Nλ_i^n>$ for all symmetry classes without any large $N$ approximation. In particular, our results provide exact expressions for moments of transmission eigenvalues in chaotic cavities with time-reversal or spin-flip symmetry and supporting a finite and arbitrary number of electronic channels in the two incoming leads.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Number of relevant directions in Principal Component Analysis and Wishart random matrices

We compute analytically, for large $N$, the probability $\mathcal{P}(N_+,N)$ that a $N\times N$ Wishart random matrix has $N_+$ eigenvalues exceeding a threshold $Nζ$, including its large deviation tails. This probability plays a benchmark role when performing the Principal Component Analysis of a large empirical dataset. We find that $\mathcal{P}(N_+,N)\approx\exp(-βN^2 ψ_ζ(N_+/N))$, where $β$ is the Dyson index of the ensemble and $ψ_ζ(κ)$ is a rate function that we compute explicitly in the full range $0\leq κ\leq 1$ and for any $ζ$. The rate function $ψ_ζ(κ)$ displays a quadratic behavior modulated by a logarithmic singularity close to its minimum $κ^\star(ζ)$. This is shown to be a consequence of a phase transition in an associated Coulomb gas problem. The variance $Δ(N)$ of the number of relevant components is also shown to grow universally (independent of $ζ)$ as $Δ(N)\sim (βπ^2)^{-1}\ln N$ for large $N$.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Phase transitions in the distribution of the Andreev conductance of superconductor-metal junctions with multiple transverse modes

We compute analytically the full distribution of Andreev conductance $G_{\mathrm{NS}}$ of a metal-superconductor interface with a large number $N_c$ of transverse modes, using a random matrix approach. The probability distribution $\mathcal{P}(G_{\mathrm{NS}},N_c)$ in the limit of large $N_c$ displays a Gaussian behavior near the average value $<G_{\mathrm{NS}}>= (2-\sqrt{2}) N_c$ and asymmetric power-law tails in the two limits of very small and very large $G_{\mathrm{NS}}$. In addition, we find a novel third regime sandwiched between the central Gaussian peak and the power law tail for large $G_{\mathrm{NS}}$. Weakly non-analytic points separate these four regimes---these are shown to be consequences of three phase transitions in an associated Coulomb gas problem.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Asymptotics of Selberg-like integrals: The unitary case and Newton's interpolation formula

We investigate the asymptotic behavior of the Selberg-like integral $$ \frac1{N!}\int_{[0,1]^N}x_1^p\prod_{i<j}(x_i-x_j)^2\prod_ix_i^{a-1}(1-x_i)^{b-1}dx_i$$, as $N\to\infty$ for different scalings of the parameters $a$ and $b$ with $N$. Integrals of this type arise in the random matrix theory of electronic scattering in chaotic cavities supporting $N$ channels in the two attached leads. Making use of Newton's interpolation formula, we show that an asymptotic limit exists and we compute it explicitly.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Entangled random pure states with orthogonal symmetry: exact results

We compute analytically the density $\varrho_{N,M}(λ)$ of Schmidt eigenvalues, distributed according to a fixed-trace Wishart-Laguerre measure, and the average Rényi entropy $\langle\mathcal{S}_q\rangle$ for reduced density matrices of entangled random pure states with orthogonal symmetry $(β=1)$. The results are valid for arbitrary dimensions $N=2k,M$ of the corresponding Hilbert space partitions, and are in excellent agreement with numerical simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

How many eigenvalues of a Gaussian random matrix are positive?

We study the probability distribution of the index ${\mathcal N}_+$, i.e., the number of positive eigenvalues of an $N\times N$ Gaussian random matrix. We show analytically that, for large $N$ and large $\mathcal{N}_+$ with the fraction $0\le c=\mathcal{N}_+/N\le 1$ of positive eigenvalues fixed, the index distribution $\mathcal{P}({\mathcal N}_+=cN,N)\sim\exp[-βN^2 Φ(c)]$ where $β$ is the Dyson index characterizing the Gaussian ensemble. The associated large deviation rate function $Φ(c)$ is computed explicitly for all $0\leq c \leq 1$. It is independent of $β$ and displays a quadratic form modulated by a logarithmic singularity around $c=1/2$. As a consequence, the distribution of the index has a Gaussian form near the peak, but with a variance $Δ(N)$ of index fluctuations growing as $Δ(N)\sim \log N/βπ^2$ for large $N$. For $β=2$, this result is independently confirmed against an exact finite $N$ formula, yielding $Δ(N)= \log N/2π^2 +C+\mathcal{O}(N^{-1})$ for large $N$, where the constant $C$ has the nontrivial value $C=(γ+1+3\log 2)/2π^2\simeq 0.185248...$ and $γ=0.5772...$ is the Euler constant. We also determine for large $N$ the probability that the interval $[ζ_1,ζ_2]$ is free of eigenvalues. Part of these results have been announced in a recent letter [\textit{Phys. Rev. Lett.} {\bf 103}, 220603 (2009)].

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Largest Schmidt eigenvalue of entangled random pure states and conductance distribution in chaotic cavities

A strategy to evaluate the distribution of the largest Schmidt eigenvalue for entangled random pure states of bipartite systems is proposed. We point out that the multiple integral defining the sought quantity for a bipartition of sizes N, M is formally identical (upon simple algebraic manipulations) to the one providing the probability density of Landauer conductance in open chaotic cavities supporting N and M electronic channels in the two leads. Known results about the latter can then be straightforwardly employed in the former problem for both systems with broken (β = 2) and preserved (β = 1) time reversal symmetry. The analytical results, yielding a continuous but not everywhere analytic distribution, are in excellent agreement with numerical simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Probability distributions of Linear Statistics in Chaotic Cavities and associated phase transitions

We establish large deviation formulas for linear statistics on the $N$ transmission eigenvalues $\{T_i\}$ of a chaotic cavity, in the framework of Random Matrix Theory. Given any linear statistics of interest $A=\sum_{i=1}^N a(T_i)$, the probability distribution $\mathcal{P}_A(A,N)$ of $A$ generically satisfies the large deviation formula $\lim_{N\to\infty}[-2\log\mathcal{P}_A(Nx,N)/βN^2]=Ψ_A(x)$, where $Ψ_A(x)$ is a rate function that we compute explicitly in many cases (conductance, shot noise, moments) and $β$ corresponds to different symmetry classes. Using these large deviation expressions, it is possible to recover easily known results and to produce new formulas, such as a closed form expression for $v(n)=\lim_{N\to\infty}\mathrm{var}(\mathcal{T}_n)$ (where $\mathcal{T}_n=\sum_{i}T_i^n$) for arbitrary integer $n$. The universal limit $v^\star=\lim_{n\to\infty} v(n)=1/2πβ$ is also computed exactly. The distributions display a central Gaussian region flanked on both sides by non-Gaussian tails. At the junction of the two regimes, weakly non-analytical points appear, a direct consequence of phase transitions in an associated Coulomb gas problem. Numerical checks are also provided, which are in full agreement with our asymptotic results in both real and Laplace space even for moderately small $N$. Part of the results have been announced in [P. Vivo, S.N. Majumdar and O. Bohigas, {\it Phys. Rev. Lett.} {\bf 101}, 216809 (2008)].

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Random Matrix approach to collective behavior and bulk universality in protein dynamics

Covariance matrices of amino acid displacements, commonly used to characterize the large-scale movements of proteins, are investigated through the prism of Random Matrix Theory. Bulk universality is detected in the local spacing statistics of noise-dressed eigenmodes, which is well described by a Brody distribution with parameter $β\simeq 0.8$. This finding, supported by other consistent indicators, implies a novel quantitative criterion to single out the collective degrees of freedom of the protein from the majority of high-energy, localized vibrations.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Nonlinear Random Matrix Statistics, symmetric functions and hyperdeterminants

Nonlinear statistics (i.e. statistics of permanents) on the eigenvalues of invariant random matrix models are considered for the three Dyson's symmetry classes $β=1,2,4$. General formulas in terms of hyperdeterminants are found for $β=2$. For specific cases and all $β$s, more computationally efficient results are obtained, based on symmetric functions expansions. As an application, we consider the case of quantum transport in chaotic cavities extending results from [D.V. Savin, H.-J. Sommers and W. Wieczorek, {\it Phys. Rev. B} {\bf 77}, 125332 (2008)].

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The Index Distribution of Gaussian Random Matrices

We compute analytically, for large N, the probability distribution of the number of positive eigenvalues (the index N_{+}) of a random NxN matrix belonging to Gaussian orthogonal (β=1), unitary (β=2) or symplectic (β=4) ensembles. The distribution of the fraction of positive eigenvalues c=N_{+}/N scales, for large N, as Prob(c,N)\simeq\exp[-βN^2 Φ(c)] where the rate function Φ(c), symmetric around c=1/2 and universal (independent of $β$), is calculated exactly. The distribution has non-Gaussian tails, but even near its peak at c=1/2 it is not strictly Gaussian due to an unusual logarithmic singularity in the rate function.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Universal Correlations and Power-Law Tails in Financial Covariance Matrices

Signatures of universality are detected by comparing individual eigenvalue distributions and level spacings from financial covariance matrices to random matrix predictions. A chopping procedure is devised in order to produce a statistical ensemble of asset-price covariances from a single instance of financial data sets. Local results for the smallest eigenvalue and individual spacings are very stable upon reshuffling the time windows and assets. They are in good agreement with the universal Tracy-Widom distribution and Wigner surmise, respectively. This suggests a strong degree of robustness especially in the low-lying sector of the spectra, most relevant for portfolio selections. Conversely, the global spectral density of a single covariance matrix as well as the average over all unfolded nearest-neighbour spacing distributions deviate from standard Gaussian random matrix predictions. The data are in fair agreement with a recently introduced generalised random matrix model, with correlations showing a power-law decay.

0 citations0 reviewsNo code linkedOpen access