Source author record

Zdzislaw Burda

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

ResearcherUnclaimed source record

math-ph math.MP cond-mat.stat-mech physics.soc-ph cond-mat.other hep-th math.ST nlin.CD physics.comp-ph physics.data-an q-fin.PM Quantitative Methods Statistics Theory

Catalog footprint

What is connected

11works

13topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

preprint2022arXiv

Cleaning large-dimensional covariance matrices for correlated samples

We elucidate the problem of estimating large-dimensional covariance matrices in the presence of correlations between samples. To this end, we generalize the Marcenko-Pastur equation and the Ledoit-Peche shrinkage estimator using methods of random matrix theory and free probability. We develop an efficient algorithm that implements the corresponding analytic formulas, based on the Ledoit-Wolf kernel estimation technique. We also provide an associated open-source Python library, called "shrinkage", with a user-friendly API to assist in practical tasks of estimation of large covariance matrices. We present an example of its usage for synthetic data generated according to exponentially-decaying auto-correlations.

preprint2022arXiv

Perfect cycles in the synchronous Heider dynamics in complete network

We discuss a cellular automaton simulating the process of reaching Heider balance in a fully connected network. The dynamics of the automaton is defined by a deterministic, synchronous and global update rule. The dynamics has a very rich spectrum of attractors including fixed points and limit cycles, the length and number of which change with the size of the system. In this paper we concentrate on a class of limit cycles that preserve energy spectrum of the consecutive states. We call such limit cycles perfect. Consecutive states in a perfect cycle are separated from each other by the same Hamming distance. Also the Hamming distance between any two states separated by $k$ steps in a perfect cycle is the same for all such pairs of states. The states of a perfect cycle form a very symmetric trajectory in the configuration space. We argue that the symmetry of the trajectories is rooted in the permutation symmetry of vertices of the network and a local symmetry of a certain energy function measuring the level of balance/frustration of triads.

preprint2015arXiv

Quaternionic R transform and non-hermitian random matrices

Using the Cayley-Dickson construction we rephrase and review the non-hermitian diagrammatic formalism [R. A. Janik, M. A. Nowak, G. Papp and I. Zahed, Nucl.Phys. B $\textbf{501}$, 603 (1997)], that generalizes the free probability calculus to asymptotically large non-hermitian random matrices. The main object in this generalization is a quaternionic extension of the R transform which is a generating function for planar (non-crossing) cumulants. We demonstrate that the quaternionic R transform generates all connected averages of all distinct powers of $X$ and its hermitian conjugate $X^\dagger$: $\langle\langle \frac{1}{N} \mbox{Tr} X^{a} X^{\dagger b} X^c \ldots \rangle\rangle$ for $N\rightarrow \infty$. We show that the R transform for gaussian elliptic laws is given by a simple linear quaternionic map $\mathcal{R}(z+wj) = x + σ^2 \left(μe^{2iϕ} z + w j\right)$ where $(z,w)$ is the Cayley-Dickson pair of complex numbers forming a quaternion $q=(z,w)\equiv z+ wj$. This map has five real parameters $\Re e x$, $\Im m x$, $ϕ$, $σ$ and $μ$. We use the R transform to calculate the limiting eigenvalue densities of several products of gaussian random matrices.

preprint2015arXiv

Unveiling the significance of eigenvectors in diffusing non-hermitian matrices by identifying the underlying Burgers dynamics

Following our recent letter, we study in detail an entry-wise diffusion of non-hermitian complex matrices. We obtain an exact partial differential equation (valid for any matrix size $N$ and arbitrary initial conditions) for evolution of the averaged extended characteristic polynomial. The logarithm of this polynomial has an interpretation of a potential which generates a Burgers dynamics in quaternionic space. The dynamics of the ensemble in the large $N$ is completely determined by the coevolution of the spectral density and a certain eigenvector correlation function. This coevolution is best visible in an electrostatic potential of a quaternionic argument built of two complex variables, the first of which governs standard spectral properties while the second unravels the hidden dynamics of eigenvector correlation function. We obtain general large $N$ formulas for both spectral density and 1-point eigenvector correlation function valid for any initial conditions. We exemplify our studies by solving three examples, and we verify the analytic form of our solutions with numerical simulations.

preprint2014arXiv

Universal distribution of Lyapunov exponents for products of Ginibre matrices

Starting from exact analytical results on singular values and complex eigenvalues of products of independent Gaussian complex random $N\times N$ matrices also called Ginibre ensemble we rederive the Lyapunov exponents for an infinite product. We show that for a large number $t$ of product matrices the distribution of each Lyapunov exponent is normal and compute its $t$-dependent variance as well as corrections in a $1/t$ expansion. Originally Lyapunov exponents are defined for singular values of the product matrix that represents a linear time evolution. Surprisingly a similar construction for the moduli of the complex eigenvalues yields the very same exponents and normal distributions to leading order. We discuss a general mechanism for $2\times 2$ matrices why the singular values and the radii of complex eigenvalues collapse onto the same value in the large-$t$ limit. Thereby we rederive Newman's triangular law which has a simple interpretation as the radial density of complex eigenvalues in the circular law and study the commutativity of the two limits $t\to\infty$ and $N\to\infty$ on the global and the local scale. As a mathematical byproduct we show that a particular asymptotic expansion of a Meijer G-function with large index leads to a Gaussian.

preprint2014arXiv

Universal microscopic correlation functions for products of truncated unitary matrices

We investigate the spectral properties of the product of $M$ complex non-Hermitian random matrices that are obtained by removing $L$ rows and columns of larger unitary random matrices uniformly distributed on the group ${\rm U}(N+L)$. Such matrices are called truncated unitary matrices or random contractions. We first derive the joint probability distribution for the eigenvalues of the product matrix for fixed $N,\ L$, and $M$, given by a standard determinantal point process in the complex plane. The weight however is non-standard and can be expressed in terms of the Meijer G-function. The explicit knowledge of all eigenvalue correlation functions and the corresponding kernel allows us to take various large $N$ (and $L$) limits at fixed $M$. At strong non-unitarity, with $L/N$ finite, the eigenvalues condense on a domain inside the unit circle. At the edge and in the bulk we find the same universal microscopic kernel as for a single complex non-Hermitian matrix from the Ginibre ensemble. At the origin we find the same new universality classes labelled by $M$ as for the product of $M$ matrices from the Ginibre ensemble. Keeping a fixed size of truncation, $L$, when $N$ goes to infinity leads to weak non-unitarity, with most eigenvalues on the unit circle as for unitary matrices. Here we find a new microscopic edge kernel that generalizes the known results for M=1. We briefly comment on the case when each product matrix results from a truncation of different size $L_j$.

preprint2013arXiv

Collective Correlations of Brodmann Areas fMRI Study with RMT-Denoising

We study collective behavior of Brodmann regions of human cerebral cortex using functional Magnetic Resonance Imaging (fMRI) and Random Matrix Theory (RMT). The raw fMRI data is mapped onto the cortex regions corresponding to the Brodmann areas with the aid of the Talairach coordinates. Principal Component Analysis (PCA) of the Pearson correlation matrix for 41 different Brodmann regions is carried out to determine their collective activity in the idle state and in the active state stimulated by tapping. The collective brain activity is identified through the statistical analysis of the eigenvectors to the largest eigenvalues of the Pearson correlation matrix. The leading eigenvectors have a large participation ratio. This indicates that several Broadmann regions collectively give rise to the brain activity associated with these eigenvectors. We apply random matrix theory to interpret the underlying multivariate data.

preprint2013arXiv

Free products of large random matrices - a short review of recent developments

We review methods to calculate eigenvalue distributions of products of large random matrices. We discuss a generalization of the law of free multiplication to non-Hermitian matrices and give a couple of examples illustrating how to use these methods in practice. In particular we calculate eigenvalue densities of products of Gaussian Hermitian and non-Hermitian matrices including combinations of GUE and Ginibre matrices.

preprint2007arXiv

Random Levy Matrices Revisited

We compare eigenvalue densities of Wigner random matrices whose elements are independent identically distributed (iid) random numbers with a Levy distribution and maximally random matrices with a rotationally invariant measure exhibiting a power law spectrum given by stable laws of free random variables. We compute the eigenvalue density of Wigner-Levy (WL) matrices using (and correcting) the method by Bouchaud and Cizeau (BC), and of free random Levy (FRL) rotationally invariant matrices by adapting results of free probability calculus. We compare the two types of eigenvalue spectra. Both ensembles are spectrally stable with respect to the matrix addition. The discussed ensemble of FRL matrices is maximally random in the sense that it maximizes Shannon's entropy. We find a perfect agreement between the numerically sampled spectra and the analytical results already for matrices of dimension N=100. The numerical spectra show very weak dependence on the matrix size N as can be noticed by comparing spectra for N=400. After a pertinent rescaling spectra of Wigner-Levy matrices and of symmetric FRL matrices have the same tail behavior. As we discuss towards the end of the paper the correlations of large eigenvalues in the two ensembles are however different. We illustrate the relation between the two types of stability and show that the addition of many randomly rotated Wigner-Levy matrices leads by a matrix central limit theorem to FRL spectra, providing an explicit realization of the maximal randomness principle.

preprint2006arXiv

Spectral properties of empirical covariance matrices for data with power-law tails

We present an analytic method for calculating spectral densities of empirical covariance matrices for correlated data. In this approach the data is represented as a rectangular random matrix whose columns correspond to sampled states of the system. The method is applicable to a class of random matrices with radial measures including those with heavy (power-law) tails in the probability distribution. As an example we apply it to a multivariate Student distribution.

preprint2005arXiv

Correlated Wishart Matrices and Critical Horizons

We discuss a practical method to determine the eigenvalue spectrum of the empirical correlation matrix. The method is based on the analysis of the behavior of a conformal map at a critical horizon which is defined as a border line of the physical Riemann sheet of this map. The map is a convenient representation of the Marčenko-Pastur equation.