BZPEER
More

Discover

About BZPEERSearchGraph

Research tools

HomeCollectionsResearchCollaborationPublish

Network

ExpertsCommunities

Opportunities

Listings

Account

Log inCreate account
Log inCreate account

Source author record

Romain Allez

Explore this author’s publications, collaborators and research topics. This catalog record has not been claimed by its author.

ResearcherUnclaimed source record
math.PRcond-mat.stat-mechq-fin.STcond-mat.dis-nnq-fin.MFq-fin.RM

Catalog footprint

Research overview

17works
6topics
4close collaborators

Actions

Is this your work?

Log in to claim
Open graphBrowse papers

Published work

17 publications

preprint2016arXiv

Rotational invariant estimator for general noisy matrices

We investigate the problem of estimating a given real symmetric signal matrix C\textbf{C} from a noisy observation matrix M\textbf{M} in the limit of large dimension. We consider the case where the noisy measurement M\textbf{M} comes either from an arbitrary additive or multiplicative rotational invariant perturbation. We establish, using the Replica method, the asymptotic global law estimate for three general classes of noisy matrices, significantly extending previously obtained results. We give exact results concerning the asymptotic deviations (called overlaps) of the perturbed eigenvectors away from the true ones, and we explain how to use these overlaps to "clean" the noisy eigenvalues of M\textbf{M}. We provide some numerical checks for the different estimators proposed in this paper and we also make the connection with some well known results of Bayesian statistics.

Open access
preprint2015arXiv

Eigenvector dynamics under free addition

We investigate the evolution of a given eigenvector of a symmetric (deterministic or random) matrix under the addition of a matrix in the Gaussian orthogonal ensemble. We quantify the overlap between this single vector with the eigenvectors of the initial matrix and identify precisely a "Cauchy-flight" regime. In particular, we compute the local density of this vector in the eigenvalues space of the initial matrix. Our results are obtained in a non perturbative setting and are derived using the ideas of [O. Ledoit and S. Péché, Prob. Th. Rel. Fields, {\bf 151} 233 (2011)]. Finally, we give a robust derivation of a result obtained in [R. Allez and J.-P. Bouchaud, Phys. Rev. E {\bf 86}, 046202 (2012)] to study eigenspace dynamics in a semi-perturbative regime.

Open access
preprint2015arXiv

The continuous Anderson hamiltonian in dimension two

We define the Anderson hamiltonian on the two dimensional torus R2/Z2\mathbb R^2/\mathbb Z^2. This operator is formally defined as H:=−Δ+ξ\mathscr H:= -Δ+ ξ where ΔΔ is the Laplacian operator and where ξξ belongs to a general class of singular potential which includes the Gaussian white noise distribution. We use the notion of paracontrolled distribution as introduced by Gubinelli, Imkeller and Perkowski in [14]. We are able to define the Schrödinger operator H\mathscr H as an unbounded self-adjoint operator on L2(T2)L^2(\mathbb T^2) and we prove that its real spectrum is discrete with no accumulation points for a general class of singular potential ξξ. We also establish that the spectrum is a continuous function of a sort of enhancement Ξ(ξ)Ξ(ξ) of the potential ξξ. As an application, we prove that a correctly renormalized smooth approximations Hε:=−Δ+ξε+cε\mathscr H_\varepsilon:= -Δ+ ξ_\varepsilon+c_\varepsilon (where ξεξ_\varepsilon is a smooth mollification of the Gaussian white noise ξξ and cεc_\varepsilon an explicit diverging renormalization constant) converge in the sense of the resolvent towards the singular operator H\mathscr H. In the case of a Gaussian white noise ξξ, we obtain exponential tail bounds for the minimal eigenvalue (sometimes called ground state) of the operator H\mathscr H as well as its order of magnitude log⁡L\log L when the operator is considered on a large box TL:=R2/(LZ)2\mathbb T_L:= \mathbb R^2/(L\mathbb Z)^2 with L→∞L\to \infty.

Open access
preprint2015arXiv

The eigenvectors of Gaussian matrices with an external source

We consider a diffusive matrix process (Xt)t≥0(X_t)_{t\ge 0} defined as Xt:=A+HtX_t:=A+H_t where AA is a given deterministic Hermitian matrix and (Ht)t≥0(H_t)_{t\ge 0} is a Hermitian Brownian motion. The matrix AA is the "external source" that one would like to estimate from the noisy observation XtX_t at some time t>0t>0. We investigate the relationship between the non-perturbed eigenvectors of the matrix AA and the perturbed eigenstates at some time tt for the three relevant scaling relations between the time tt and the dimension NN of the matrix XtX_t. We determine the asymptotic (mean-squared) projections of any given non-perturbed eigenvector ∣ψj0⟩|ψ_j^0\rangle, associated to an eigenvalue aja_j of AA which may lie inside the bulk of the spectrum or be isolated (spike) from the other eigenvalues, on the orthonormal basis of the perturbed eigenvectors ∣ψit⟩,i≠j|ψ_i^t\rangle,i\neq j. We derive a Burgers type evolution equation for the local resolvent (z−Xt)ii−1(z-X_t)_{ii}^{-1}, describing the evolution of the local density of a given initial state ∣ψj0⟩|ψ_j ^0\rangle. We are able to solve this equation explicitly in the large NN limit, for any initial matrix AA. In the case of one isolated eigenvector ∣ψj0⟩|ψ_j^0\rangle, we prove a central limit Theorem for the overlap ⟨ψj0∣ψjt⟩\langle ψ_j^0|ψ_j^t\rangle. When properly centered and rescaled by a factor N\sqrt{N}, this overlap converges in law towards a centered Gaussian distribution with an explicit variance depending on tt. Our method is based on analyzing the eigenvector flow under the Dyson Brownian motion.

Open access
preprint2014arXiv

From Sine kernel to Poisson statistics

We study the Sineβ_β process introduced in [B. Valkó and B. Virág. Invent. math. (2009)] when the inverse temperature ββ tends to 0. This point process has been shown to be the scaling limit of the eigenvalues point process in the bulk of ββ-ensembles and its law is characterized in terms of the winding numbers of the Brownian carrousel at different angular speeds. After a careful analysis of this family of coupled diffusion processes, we prove that the Sineβ_β point process converges weakly to a Poisson point process on R\mathbb{R}. Thus, the Sineβ_β point processes establish a smooth crossover between the rigid clock (or picket fence) process (corresponding to β=∞β=\infty) and the Poisson process.

Open access
preprint2014arXiv

Index Distribution of the Ginibre Ensemble

Complex systems, and in particular random neural networks, are often described by randomly interacting dynamical systems with no specific symmetry. In that context, characterizing the number of relevant directions necessitates fine estimates on the Ginibre ensemble. In this Letter, we compute analytically the probability distribution of the number of eigenvalues NRN_R with modulus greater than RR (the index) of a large N×NN\times N random matrix in the real or complex Ginibre ensemble. We show that the fraction NR/N=pN_R/N=p has a distribution scaling as exp⁡(−βN2ψR(p))\exp(-βN^2 ψ_R(p)) with β=1β=1 (respectively β=1/2β=1/2) for the complex (resp. real) Ginibre ensemble. For any p∈[0,1]p\in[0,1], the equilibrium spectral densities as well as the rate function ψR(p)ψ_R(p) are explicitly derived. This function displays a third order phase transition at the critical (minimum) value pR∗=1−R2p^*_R=1-R^2, associated to a phase transition of the Coulomb gas. We deduce that, in the central regime, the fluctuations of the index NRN_R around its typical value pR∗Np^*_R N scale as N1/3N^{1/3}.

Open access
preprint2014arXiv

Random matrices in non-confining potentials

We consider invariant matrix processes diffusing in non-confining cubic potentials of the form Va(x)=x3/3−ax,a∈RV_a(x)= x^3/3 - a x, a\in \mathbb{R}. We construct the trajectories of such processes for all time by restarting them whenever an explosion occurs, from a new (well chosen) initial condition, insuring continuity of the eigenvectors and of the non exploding eigenvalues. We characterize the dynamics of the spectrum in the limit of large dimension and analyze the stationary state of this evolution explicitly. We exhibit a sharp phase transition for the limiting spectral density ρaρ_a at a critical value a=a∗a=a^*. If a≥a∗a\geq a^*, then the potential VaV_a presents a well near x=ax=\sqrt{a} deep enough to confine all the particles inside, and the spectral density ρaρ_a is supported on a compact interval. If a<a∗a<a^* however, the steady state is in fact dynamical with a macroscopic stationary flux of particles flowing across the system. In this regime, the eigenvalues allocate according to a stationary density profile ρaρ_{a} with full support in R\mathbb{R}, flanked with heavy tails such that ρa(x)∼Ca/x2ρ_{a}(x)\sim C_a /x^2 as x→±∞x\to \pm \infty. Our method applies to other non-confining potentials and we further investigate a family of quartic potentials, which were already studied in Brézin et al. to count planar diagrams.

Open access
preprint2014arXiv

Tracy-Widom at high temperature

We investigate the marginal distribution of the bottom eigenvalues of the stochastic Airy operator when the inverse temperature ββ tends to 00. We prove that the minimal eigenvalue, whose fluctuations are governed by the Tracy-Widom ββ 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 kk. 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 [L. Dumaz and B. 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 β→0β\to 0. As an application, we investigate the maximal eigenvalues statistics of βNβ_N-ensembles when the repulsion parameter βN→0β_N\to 0 when N→+∞N\to +\infty. We study the double scaling limit N→+∞,βN→0N\to +\infty, β_N \to 0 and argue with heuristic and numerical arguments that the statistics of the marginal distributions can be deduced following the ideas of [A. Edelman and B. D. Sutton. J. Stat. Phys. 127, 6, 1121-1165 (2007)] and [J. A. Ramírez, B. Rider and B. Virág. J. Amer. Math. Soc. 24, 919-944 (2011)] from our later study of the stochastic Airy operator.

Open access
preprint2012arXiv

A diffusive matrix model for invariant ββ-ensembles

We define a new diffusive matrix model converging towards the ββ-Dyson Brownian motion for all β∈[0,2]β\in [0,2] that provides an explicit construction of ββ-ensembles of random matrices that is invariant under the orthogonal/unitary group. We also describe the eigenvector dynamics of the limiting matrix process; we show that when β<1β< 1 and that two eigenvalues collide, the eigenvectors of these two colliding eigenvalues fluctuate very fast and take the uniform measure on the orthocomplement of the eigenvectors of the remaining eigenvalues.

Open access
preprint2012arXiv

Eigenvector dynamics: general theory and some applications

We propose a general framework to study the stability of the subspace spanned by PP consecutive eigenvectors of a generic symmetric matrix H0{\bf H}_0, when a small perturbation is added. This problem is relevant in various contexts, including quantum dissipation (H0{\bf H}_0 is then the Hamiltonian) and financial risk control (in which case H0{\bf H}_0 is the assets return covariance matrix). We argue that the problem can be formulated in terms of the singular values of an overlap matrix, that allows one to define a "fidelity" distance. We specialize our results for the case of a Gaussian Orthogonal H0{\bf H}_0, for which the full spectrum of singular values can be explicitly computed. We also consider the case when H0{\bf H}_0 is a covariance matrix and illustrate the usefulness of our results using financial data. The special case where the top eigenvalue is much larger than all the other ones can be investigated in full detail. In particular, the dynamics of the angle made by the top eigenvector and its true direction defines an interesting new class of random processes.

Open access
preprint2012arXiv

Invariant ββ-ensembles and the Gauss-Wigner crossover

We define a new diffusive matrix model converging towards the ββ -Dyson Brownian motion for all β∈[0,2]β\in [0,2] that provides an explicit construction of ββ-ensembles of random matrices that is invariant under the orthogonal/unitary group. For small values of ββ, our process allows one to interpolate smoothly between the Gaussian distribution and the Wigner semi-circle. The interpolating limit distributions form a one parameter family that can be explicitly computed. This also allows us to compute the finite-size corrections to the semi-circle.

Open 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 β∈[0,2]β\in [0,2], which preserves invariance under the orthogonal/unitary group transformation. Scaling the Dyson index ββ with the largest size MM of the data matrix as β=2c/Mβ=2c/M (with cc a fixed positive constant), we obtain a family of spectral densities parametrized by cc. As cc is varied, this density interpolates continuously between the Mar\vcenko-Pastur (c→∞c\to \infty limit) and the Gamma law (c→0c\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/M21/M^2 for the particular cases β=1,2β=1,2.

Open access
preprint2012arXiv

Lognormal scale invariant random measures

In this article, we consider the continuous analog of the celebrated Mandelbrot star equation with lognormal weights. Mandelbrot introduced this equation to characterize the law of multiplicative cascades. We show existence and uniqueness of measures satisfying the aforementioned continuous equation; these measures fall under the scope of the Gaussian multiplicative chaos theory developed by J.P. Kahane in 1985 (or possibly extensions of this theory). As a by product, we also obtain an explicit characterization of the covariance structure of these measures. We also prove that qualitative properties such as long-range independence or isotropy can be read off the equation.

Open access
preprint2012arXiv

Marchenko Pastur type theorem for independent MRW processes: convergence of the empirical spectral measure

We study the asymptotic of the spectral distribution for large empirical covariance matrices composed of independent Multifractal Random Walk processes. The asymptotic is taken as the observation lag shrinks to 0. In this setting, we show that there exists a limiting spectral distribution whose Stieltjes transform is uniquely characterized by equations which we specify. We also illustrate our results by numerical simulations.

Open access
preprint2011arXiv

Eigenvector dynamics: theory and some applications

We propose a general framework to study the stability of the subspace spanned by PP consecutive eigenvectors of a generic symmetric matrix H0{\bf H}_0, when a small perturbation is added. This problem is relevant in various contexts, including quantum dissipation (H0{\bf H}_0 is then the Hamiltonian) and risk control (in which case H0{\bf H}_0 is the assets return correlation matrix). We specialize our results for the case of a Gaussian Orthogonal H0{\bf H}_0, or when H0{\bf H}_0 is a correlation matrix. We illustrate the usefulness of our framework using financial data.

Open access
preprint2011arXiv

Principal Regression Analysis and the index leverage effect

We revisit the index leverage effect, that can be decomposed into a volatility effect and a correlation effect. We investigate the latter using a matrix regression analysis, that we call `Principal Regression Analysis' (PRA) and for which we provide some analytical (using Random Matrix Theory) and numerical benchmarks. We find that downward index trends increase the average correlation between stocks (as measured by the most negative eigenvalue of the conditional correlation matrix), and makes the market mode more uniform. Upward trends, on the other hand, also increase the average correlation between stocks but rotates the corresponding market mode {\it away} from uniformity. There are two time scales associated to these effects, a short one on the order of a month (20 trading days), and a longer time scale on the order of a year. We also find indications of a leverage effect for sectorial correlations as well, which reveals itself in the second and third mode of the PRA.

Open access
preprint2010arXiv

Individual and collective stock dynamics: intra-day seasonalities

We establish several new stylised facts concerning the intra-day seasonalities of stock dynamics. Beyond the well known U-shaped pattern of the volatility, we find that the average correlation between stocks increases throughout the day, leading to a smaller relative dispersion between stocks. Somewhat paradoxically, the kurtosis (a measure of volatility surprises) reaches a minimum at the open of the market, when the volatility is at its peak. We confirm that the dispersion kurtosis is a markedly decreasing function of the index return. This means that during large market swings, the idiosyncratic component of the stock dynamics becomes sub-dominant. In a nutshell, early hours of trading are dominated by idiosyncratic or sector specific effects with little surprises, whereas the influence of the market factor increases throughout the day, and surprises become more frequent.

Open access

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.