Spectra of Random Hermitian Matrices with a Small-Rank External Source: The Supercritical and Subcritical Regimes

Random Hermitian matrices with a source term arise, for instance, in the study of non-intersecting Brownian walkers and sample covariance matrices. We consider the case when the n×n external source matrix has two distinct real eigenvalues: a with multiplicity r and zero with multiplicity n?r. The source is small in the sense that r is finite or $r=\mathcal{O}(n^{\gamma})$ , for 0<γ<1. For a Gaussian potential, Péché (Probab. Theory Relat. Fields 134:127–173, 2006) showed that for |a| sufficiently small (the subcritical regime) the external source has no leading-order effect on the eigenvalues, while for |a| sufficiently large (the supercritical regime) r eigenvalues exit the bulk of the spectrum and behave as the eigenvalues of the r×r Gaussian unitary ensemble (GUE). We establish the universality of these results for a general class of analytic potentials in the supercritical and subcritical regimes.     

Large n Limit of Gaussian Random Matrices with External Source, Part I

We consider the random matrix ensemble with an external sourcedefined on n×n Hermitian matrices, where A is a diagonal matrix with only two eigenvalues ±a of equal multiplicity. For the case a>1, we establish the universal behavior of local eigenvalue correlations in the limit n, which is known from unitarily invariant random matrix models. Thus, local eigenvalue correlations are expressed in terms of the sine kernel in the bulk and in terms of the Airy kernel at the edge of the spectrum. We use a characterization of the associated multiple Hermite polynomials by a 3×3-matrix Riemann-Hilbert problem, and the Deift/Zhou steepest descent method to analyze the Riemann-Hilbert problem in the large n limit.Dedicated to Freeman Dyson on his eightieth birthdayThe first author was supported in part by NSF Grants DMS-9970625 and DMS-0354962.The second author was supported in part by projects G.0176.02 and G.0455.04 of FWO-Flanders, by K.U.Leuven research grant OT/04/24, and by INTAS Research Network NeCCA 03-51-6637.     

The Largest Eigenvalue of Real Symmetric,Hermitian and Hermitian Self-dual Random Matrix Models with Rank One External Source,Part I

We consider the limiting location and limiting distribution of the largest eigenvalue in real symmetric (β=1), Hermitian (β=2), and Hermitian self-dual (β=4) random matrix models with rank 1 external source. They are analyzed in a uniform way by a contour integral representation of the joint probability density function of eigenvalues. Assuming the “one-band” condition and certain regularities of the potential function, we obtain the limiting location of the largest eigenvalue when the nonzero eigenvalue of the external source matrix is not the critical value, and further obtain the limiting distribution of the largest eigenvalue when the nonzero eigenvalue of the external source matrix is greater than the critical value. When the nonzero eigenvalue of the external source matrix is less than or equal to the critical value, the limiting distribution of the largest eigenvalue will be analyzed in a subsequent paper. In this paper we also give a definition of the external source model for all β>0.     

Large <Emphasis Type="Italic">n</Emphasis> Limit of Gaussian Random Matrices with External Source,Part III: Double Scaling Limit

We consider the double scaling limit in the random matrix ensemble with an external source

Large <Emphasis Type="Italic">n</Emphasis> Limit of Gaussian Random Matrices with External Source,Part II

We continue the study of the Hermitian random matrix ensemble with external source where A has two distinct eigenvalues ±a of equal multiplicity. This model exhibits a phase transition for the value a=1, since the eigenvalues of M accumulate on two intervals for a>1, and on one interval for 0<a<1. The case a>1 was treated in Part I, where it was proved that local eigenvalue correlations have the universal limiting behavior which is known for unitarily invariant random matrices, that is, limiting eigenvalue correlations are expressed in terms of the sine kernel in the bulk of the spectrum, and in terms of the Airy kernel at the edge. In this paper we establish the same results for the case 0<a<1. As in Part I we apply the Deift/Zhou steepest descent analysis to a 3×3-matrix Riemann-Hilbert problem. Due to the different structure of an underlying Riemann surface, the analysis includes an additional step involving a global opening of lenses, which is a new phenomenon in the steepest descent analysis of Riemann-Hilbert problems.The first and third author are supported in part by INTAS Research Network NeCCA 03-51-6637 and by NATO Collaborative Linkage Grant PST.CLG.979738. The first author is supported in part by RFBR 05-01-00522 and the program “Modern problems of theoretical mathematics” RAS(DMS). The second author is supported in part by the National Science Foundation (NSF) Grant DMS-0354962. The third author is supported in part by FWO-Flanders projects G.0176.02 and G.0455.04 and by K.U.Leuven research grant OT/04/24 and by the European Science Foundation Program Methods of Integrable Systems, Geometry, Applied Mathematics (MISGAM) and the European Network in Geometry, Mathematical Physics and Applications (ENIGMA)     

On Consistency of the Quantum-Like Representation Algorithm

In this paper we continue to study so-called “inverse Born’s rule problem”: to construct a representation of probabilistic data of any origin by a complex probability amplitude which matches Born’s rule. The corresponding algorithm—quantum-like representation algorithm (QLRA)—was recently proposed by A. Khrennikov (Found. Phys. 35(10):1655–1693, 2005; Physica E 29:226–236, 2005; Dokl. Akad. Nauk 404(1):33–36, 2005; J. Math. Phys. 46(6):062111–062124, 2005; Europhys. Lett. 69(5):678–684, 2005). Formally QLRA depends on the order of conditioning. For two observables (of any origin, e.g., physical or biological) a and b, b|a- and a|b conditional probabilities produce two representations, say in Hilbert spaces H ^b|a and H ^a|b . In this paper we prove that under “natural assumptions” (which hold, e.g., for quantum observables represented by operators with nondegenerate spectra) these two representations are unitary equivalent. This result proves the consistency of QLRA.     

On SLE Martingales in Boundary WZW Models

Following Bettelheim et al. (Phys Rev Lett 95:251601, 2005), we consider the boundary WZW model on a half-plane with a cut growing according to the Schramm–Loewner stochastic evolution and the boundary fields inserted at the tip of the cut and at infinity. We study necessary and sufficient conditions for boundary correlation functions to be SLE martingales. Necessary conditions come from the requirement for the boundary field at the tip of the cut to have a depth two null vector. Sufficient conditions are established using Knizhnik–Zamolodchikov equations for boundary correlators. Combining these two approaches, we show that in the case of G = SU(2) the boundary correlator is an SLE martingale if and only if the boundary field carries spin 1/2. In the case of G = SU(n) and the level k = 1, there are several situations when boundary one-point correlators are SLE _κ -martingales. If the boundary field is labelled by the defining n-dimensional representation of SU(n), we obtain \varkappa = 2{\varkappa=2} . For n even, by choosing the boundary field labelled by the (unique) self-adjoint fundamental representation, we get \varkappa = 8/(n + 2){\varkappa=8/(n {+} 2)} . We also study the situation when the distance between the two boundary fields is finite, and we show that in this case the SLE_\varkappa{{\rm SLE}_\varkappa} evolution is replaced by SLE_\varkappa,r{{\rm SLE}_{\varkappa,\rho}} with r = \varkappa -6{\rho=\varkappa -6} .     

Noncrossing Normal Ordering for Functions of Boson Operators

Normally ordered forms of functions of boson operators are important in many contexts in particular concerning Quantum Field Theory and Quantum Optics. Beginning with the seminal work of Katriel (Lett. Nuovo Cimento 10(13):565–567, 1974), in the last few years, normally ordered forms have been shown to have a rich combinatorial structure, mainly in virtue of a link with the theory of partitions. In this paper, we attempt to enrich this link. By considering linear representations of noncrossing partitions, we define the notion of noncrossing normal ordering. Given the growing interest in noncrossing partitions, because of their many unexpected connections (like, for example, with free probability), noncrossing normal ordering appears to be an intriguing notion. We explicitly give the noncrossing normally ordered form of the functions (a ^r (a ^† ) ^s ) ⁿ ) and (a ^r +(a ^† ) ^s ) ⁿ , plus various special cases. We are able to establish for the first time bijections between noncrossing contractions of these functions, k-ary trees and sets of lattice paths.     

Non-Intersecting Squared Bessel Paths: Critical Time and Double Scaling Limit

We consider the double scaling limit for a model of n non-intersecting squared Bessel processes in the confluent case: all paths start at time t = 0 at the same positive value x = a, remain positive, and are conditioned to end at time t = 1 at x = 0. After appropriate rescaling, the paths fill a region in the tx–plane as n → ∞ that intersects the hard edge at x = 0 at a critical time t = t ^*. In a previous paper, the scaling limits for the positions of the paths at time t ≠ t ^* were shown to be the usual scaling limits from random matrix theory. Here, we describe the limit as n → ∞ of the correlation kernel at critical time t ^* and in the double scaling regime. We derive an integral representation for the limit kernel which bears some connections with the Pearcey kernel. The analysis is based on the study of a 3 × 3 matrix valued Riemann-Hilbert problem by the Deift-Zhou steepest descent method. The main ingredient is the construction of a local parametrix at the origin, out of the solutions of a particular third-order linear differential equation, and its matching with a global parametrix.     

Realization of the Spectrum Generating Algebra for the Generalized Kratzer Potentials

The dynamical symmetries of the Kratzer-type molecular potentials (generalized Kratzer molecular potentials) are studied by using the factorization method. The creation and annihilation (ladder) operators for the radial eigenfunctions satisfying quantum dynamical algebra SU(1,1) are established. Factorization method is a very simple method of calculating the matrix elements from these ladder operators. The matrix elements of different functions of r, r\fracddrr\frac{d}{dr}, their sum Γ₁ and difference Γ₂ are evaluated in a closed form. The exact bound state energy eigenvalues E _n,ℓ and matrix elements of r, r\fracddrr\frac{d}{dr}, their sum Γ₁ and difference Γ₂ are calculated for various values of n and ℓ quantum numbers for CO and NO diatomic molecules for the two potentials. The results obtained are in very good agreement with those obtained by other methods.     

Inverse resonance scattering for Jacobi operators

The Jacobi operator (Jf) _n = a _n−1 f _n−1 +a _n f _n+1 + b _n f _n on ℤ with real finitely supported sequences (a _n − 1) _n∈ℤ and (b _n ) _n∈ℤ is considered. The inverse problem for two mappings (including their characterization): (a _n , b _n , n ∈ ℤ) → {the zeros of the reflection coefficient} and (a _n , b _n , n ∈ ℤ) → {the eigenvalues and the resonances} is solved. All Jacobi operators with the same eigenvalues and resonances are also described.     

Edge Universality for Orthogonal Ensembles of Random Matrices

We prove edge universality of local eigenvalue statistics for orthogonal invariant matrix models with real analytic potentials and one interval limiting spectrum. Our starting point is the result of Shcherbina (Commun. Math. Phys. 285, 957–974, 2009) on the representation of the reproducing matrix kernels of orthogonal ensembles in terms of scalar reproducing kernel of corresponding unitary ensemble.     

Integrable Fredholm Operators and Dual Isomonodromic Deformations

The Fredholm determinants of a special class of integrable integral operators K supported on the union of m curve segments in the complex λ-plane are shown to be the τ-functions of an isomonodromic family of meromorphic covariant derivative operators , having regular singular points at the 2m endpoints of the curve segments, and a singular point of Poincaré index 1 at infinity. The rank r of the corresponding vector bundle over the Riemann sphere equals the number of distinct terms in the exponential sum defining the numerator of the integral kernel. The matrix Riemann–Hilbert problem method is used to deduce an identification of the Fredholm determinant as a τ-function in the sense of Segal–Wilson and Sato, i.e., in terms of abelian group actions on the determinant line bundle over a loop space Grassmannian. An associated dual isomonodromic family of covariant derivative operators , having rank n= 2m, and r finite regular singular points located at the values of the exponents defining the kernel of K is derived. The deformation equations for this family are shown to follow from an associated dual set of Riemann–Hilbert data, in which the r?les of the r exponential factors in the kernel and the 2m endpoints of its support are interchanged. The operators are analogously associated to an integral operator whose Fredholm determinant is equal to that of K. Received: 10 June 1997 / Received revised: 16 February 2001 / Accepted: 27 November 2001     

Quantum Gates and Quantum Algorithms with Clifford Algebra Technique

We use the Clifford algebra technique (J. Math. Phys. 43:5782, 2002; J. Math. Phys. 44:4817, 2003), that is nilpotents and projectors which are binomials of the Clifford algebra objects γ ^a with the property {γ ^a ,γ ^b }₊=2η ^ab , for representing quantum gates and quantum algorithms needed in quantum computers in a simple and an elegant way. We identify n-qubits with the spinor representations of the group SO(1,3) for a system of n spinors. Representations are expressed in terms of products of projectors and nilpotents; we pay attention also on the nonrelativistic limit. An algorithm for extracting a particular information out of a general superposition of 2 ⁿ qubit states is presented. It reproduces for a particular choice of the initial state the Grover’s algorithm (Proc. 28th Annual ACM Symp. Theory Comput. 212, 1996).     

Free Energies of Dilute Bose Gases: Upper Bound

We derive an upper bound on the free energy of a Bose gas at density ϱ and temperature T. In combination with the lower bound derived previously by Seiringer (Commun. Math. Phys. 279(3): 595–636, 2008), our result proves that in the low density limit, i.e., when a ³ ϱ≪1, where a denotes the scattering length of the pair-interaction potential, the leading term of Δf, the free energy difference per volume between interacting and ideal Bose gases, is equal to 4pa(2r²-[r-r_c]²₊)4\pi a(2\varrho^{2}-[\varrho-\varrho_{c}]^{2}_{+}). Here, ϱ _c (T) denotes the critical density for Bose–Einstein condensation (for the ideal Bose gas), and [⋅]₊=max {⋅,0} denotes the positive part.     

Effect of the velocity-dependent potentials on the energy eigenvalues of the Morse potential

We investigate the effect of the isotropic velocity-dependent potentials on the bound state energy eigenvalues of the Morse potential for any quantum states. When the velocity-dependent term is used as a constant parameter, ρ(r) = ρ ₀, the energy eigenvalues can be obtained analytically by using the Pekeris approximation. When the velocity-dependent term is considered as an harmonic oscillator type, ρ(r) = ρ ₀ r ², we show how to obtain the energy eigenvalues of the Morse potential without any approximation for any n and ℓ quantum states by using numerical calculations. The calculations have been performed for different energy eigenvalues and different numerical values of ρ ₀, in order to show the contribution of the velocity-dependent potential on the energy eigenvalues of the Morse potential.     

CORAL: QSPR models for solubility of [C60] and [C70] fullerene derivatives

Quantitative structure–property relationships (QSPRs) between the molecular structure of [C₆₀] and [C₇₀] fullerene derivatives and their solubility in chlorobenzene (mg/mL) have been established by means of CORAL (CORrelations And Logic) freeware. The CORAL models are based on representation of the molecular structure by simplified molecular input line entry system (SMILES). Three random splits into the training and the external validation sets have been examined. The ranges of statistical characteristics of these models are as follows: n = 18, r ² = 0.748–0.815, s = 15.1 –17.5 (mg/mL), F = 47–71 (training set); n = 9, r ² = 0.806–0.936, s = 12.5–17.5 (mg/mL), F = 29–103 (validation set).     

On Connected Diagrams and Cumulants of Erdős-Rényi Matrix Models

Regarding the adjacency matrices of n-vertex graphs and related graph Laplacian we introduce two families of discrete matrix models constructed both with the help of the Erdős-Rényi ensemble of random graphs. Corresponding matrix sums represent the characteristic functions of the average number of walks and closed walks over the random graph. These sums can be considered as discrete analogues of the matrix integrals of random matrix theory. We study the diagram structure of the cumulant expansions of logarithms of these matrix sums and analyze the limiting expressions as n → ∞ in the cases of constant and vanishing edge probabilities.     

Local Semicircle Law at the Spectral Edge for Gaussian <Emphasis Type="Italic">β</Emphasis>-Ensembles

We study the local semicircle law for Gaussian β-ensembles at the edge of the spectrum. We prove that at the almost optimal level of n^-2/3+e{n^{-2/3+\epsilon}}, the local semicircle law holds for all β ≥ 1 at the edge. The proof of the main theorem relies on the calculation of the moments of the tridiagonal model of Gaussian β-ensembles up to the p _n -moment, where p_n = O(n^2/3-e){p_n = O(n^{2/3-\epsilon})}. The result is analogous to the result of Sinai and Soshnikov (Funct Anal Appl 32(2), 1998) for Wigner matrices, but the combinatorics involved in the calculations are different.     

Simple Collaboration Eavesdropping on the Improved Multiparty Quantum Secret Sharing Protocol

We propose a new attack strategy for the improvement n-party (n≥4) case [S. Lin, F. Gao, Q.Y. Wen, F.C. Zhu in Opt. Commun. 281:4553, 2008] of the multiparty quantum secret sharing protocol [Z.J. Zhang, G. Gao, X. Wang, L.F. Han, S.H. Shi in Opt. Commun. 269:418, 2007]. Our attack strategy is an interesting collaboration eavesdropping and much simpler than that in the paper [T.Y. Wang, Q.Y. Wen, F. Gao, S. Lin, F.C. Zhu in Phys. Lett. A 373:65, 2008].