Next Article in Journal
Measurement-Device-Independent Quantum Key Distribution Based on Decoherence-Free Subspaces with Logical Bell State Analyzer
Next Article in Special Issue
Local Unitary Equivalence of Quantum States Based on the Tensor Decompositions of Unitary Matrices
Previous Article in Journal
A Variation of the Algorithm to Achieve the Maximum Entropy for Belief Functions
Previous Article in Special Issue
Stieltjes Transforms and R-Transforms Associated with Two-Parameter Lambert–Tsallis Functions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Double Dyson Index β Effect in Non-Hermitian Tridiagonal Matrices

by
Cleverson A. Goulart
and
Mauricio P. Pato
*
Instituto de Física, Universidade de São Paulo, Caixa Postal 66318, São Paulo 05314-970, SP, Brazil
*
Author to whom correspondence should be addressed.
Entropy 2023, 25(6), 868; https://doi.org/10.3390/e25060868
Submission received: 9 April 2023 / Revised: 24 May 2023 / Accepted: 25 May 2023 / Published: 29 May 2023
(This article belongs to the Special Issue Random Matrices: Theory and Applications)

Abstract

:
The Dyson index, β , plays an essential role in random matrix theory, as it labels the so-called “three-fold way” that refers to the symmetries satisfied by ensembles under unitary transformations. As is known, its 1 , 2, and 4 values denote the orthogonal, unitary, and symplectic classes, whose matrix elements are real, complex, and quaternion numbers, respectively. It functions, therefore, as a measure of the number of independent non-diagonal variables. On the other hand, in the case of β ensembles, which represent the tridiagonal form of the theory, it can assume any real positive value, thus losing that function. Our purpose, however, is to show that, when the Hermitian condition of the real matrices generated with a given value of β is removed, and, as a consequence, the number of non-diagonal independent variables doubles, non-Hermitian matrices exist that asymptotically behave as if they had been generated with a value 2 β . Therefore, it is as if the β index were, in this way, again operative. It is shown that this effect happens for the three tridiagonal ensembles, namely, the β –Hermite, the β –Laguerre, and the β –Jacobi ensembles.

1. Introduction

The possibility of reducing a full matrix into a tridiagonal matrix by employing unitary transformations inspired some researchers [1] to search for ensembles of tridiagonal matrices that would have the same properties of the ensembles of random matrix theory (RMT) [2]. This objective has been achieved in a seminal paper [3] in which ensembles of tridiagonal Hermitian matrices are constructed for the classes associated with the Hermite and the Laguerre polynomials. The tridiagonal matrices associated with the Jacobi polynomials, on the other hand, were obtained in [4,5,6]. In these three tridiagonal ensembles, the Dyson index, β , which in RMT has the integer values 1 , 2 , and 4 , can assume any positive real value. Therefore, the β ensemble generalized Dyson’s “three-fold way” [7] associated with invariance under the orthogonal (GOE), unitary (GUE), and symplectic (GSE) transformations of the three classes of the Gaussian ensemble.
Once the β ensembles were established, it was natural to investigate what happens if its Hermitian condition is removed. Of course, this interest parallels the introduction of the non-Hermitian random matrices by Ginibre [8] a few years after Wigner’s proposal of the Hermitian random matrices by the end of the 1950s [9]. In the tridiagonal Hermite case, this investigation was undertaken in ref. [10], in which, besides the question of introducing non-Hermiticity in β ensembles, there was also another motivation.
In fact, the main interest of this study was in the special class of non-Hermitian operators in which the operator is connected to its adjoint by a similarity transformation; namely, it satisfies the condition
A = η A η 1
where A is non-Hermitian and η is Hermitian. This condition defines a pseudo-Hermitian operator [11], and it implies that its spectrum is made of real or complex conjugate eigenvalues. This kind of non-Hermitian operator attracted a lot of interest with the discovery that Hamiltonian invariants under the combined parity and time-reversal transformation, called PT symmetry, also have real or complex conjugate spectra [12,13,14]. This led to an extension of quantum mechanics to include this special class of non-Hermitian Hamiltonians with PT symmetry [14] (see [15] for a review).
Concomitantly, there have been attempts to find, in the context of random matrix theory (RMT), ensembles of matrices satisfying the above relationship [16,17]. In [10], it was shown that in its tridiagonal form, the removal of the Hermitian condition produces matrices that satisfy the above pseudo-Hermitian constraint. We are here revisiting that work to report an important result that has since been overlooked. In fact, there, it was shown that the eigenvalues of the non-Hermitian ensemble are real and distributed with the same semi-circle law as those of the Hermitian ensemble. However, not much could be said about the fluctuations. In this respect, what we have now found is that asymptotically, that is, for arbitrary large matrices, the passage to non-Hermiticity is accompanied by the doubling of Dyson’s β index. It is important to emphasize that, in order for this doubling index effect to occur, it is necessary to take into account the presence of traces in the structure of the joint distributions of matrix elements in the case of the Hermite and Laguerre ensembles. In the Jacobi ensemble, it will be shown that the important quantity is the variance of the diagonal and the off-diagonal elements. Therefore, the effect occurs in the three β ensembles. Moreover, it affects the short-range statistics, measured by the nearest-neighbor spacing distribution (NNSD), and the long-range statistics, measured by the number variance. This is illustrated with plots of NNSD for the three ensembles and number variance for the Jacobi ensemble. The figures with similar results for the two other ensembles are omitted.

2. Preliminares

We start by recalling some results from [10]. Considering the general case of a non-Hermitian tridiagonal matrix, A, with a diagonal a = ( a n , , a 1 ) , an upper sub-diagonal b = ( b n 1 , , b 1 ) , and a lower sub-diagonal c = ( c n 1 , , c 1 ) , in which elements can assume any real value, it is easily proved that the matrix A is pseudo-Hermitian. In fact, by defining the diagonal matrix η , whose elements are given by
diag ( η ) = 1 , b n 1 c n 1 , b n 1 b n 2 c n 1 c n 2 , , b n 1 b n 2 b 1 c n 1 c n 2 c 1 ,
it is immediately verified that A and its adjoint A satisfy Equation (1), belonging, therefore, to the class of pseudo-Hermitian matrices. Moreover, by defining the diagonal matrix η whose elements are obtained by taking the square roots of the elements of η , that is,
diag ( η ) = 1 , b n 1 c n 1 , b n 1 b n 2 c n 1 c n 2 , , b n 1 b n 2 b 1 c n 1 c n 2 c 1
we find that, if the products b i c i are positive, then we can define a matrix
K = η 1 2 A η 1 2 = η 1 2 η A η 1 η 1 2 = η 1 2 A η 1 2 = K
which is Hermitian and iso-spectral with A and whose diagonal is the same of A, while the sub-diagonals are the geometrical means b n 1 c n 1 , , b 1 c 1 .
Considering now the distribution
h ν ( y ) = 2 exp ( y 2 / 2 ) y ν 1 2 ν / 2 Γ [ ν / 2 ]
with a first moment < y > = 2 Γ [ ( ν + 1 ) / 2 ] Γ ( ν / 2 ) ν and a second moment < y 2 > = ν , it is found that by constructing a new variable z = x y , where x and y are both sorted, its distribution is the K-distribution
g ν ( z = x y ) = 8 2 ν Γ 2 ν / 2 z 2 ν 1 K 0 z 2 ,
where K 0 ( x ) is a modified Bessel function of the second kind. For large values of ν , the contributions to the K-distribution come from large values of its argument. This being the case, we are allowed to replace in the expression of g ν ( z ) the Bessel function with its asymptotic expression K 0 ( z ) π 2 z exp ( z ) , which leads to
g ν ( z ) 8 π 2 ν + 1 / 2 Γ 2 ν / 2 z 2 ν 2 exp z 2 .
Using now the gamma function duplication formula
Γ μ 1 2 = 2 μ 1 2 π Γ μ 2 1 4 Γ μ 2 + 1 4
and the fact that, for large ν ,
Γ μ 2 1 4 Γ μ 2 + 1 4 Γ 2 μ 2 ,
it is found that the distribution of the geometrical mean variable can be approximated as
g ν ( z ) 2 Γ ν 1 / 2 z 2 ν 2 exp z 2 = f 2 ν 1 ( z ) f 2 ν ( z ) ,
where
f ν ( y ) = h ν ( y 2 ) 2 = 2 y ν 1 Γ ( ν / 2 ) exp ( y 2 ) .
Let us now consider the beta distribution
B ( s , t ) = 2 1 s t Γ ( s + t ) Γ ( s ) Γ ( t ) ( 1 x ) s 1 ( 1 + x ) t 1
in the asymptotic situation in which s > > 1 and t > > 1 . In this case, we can make the expansion
F ( x ) = ( s 1 ) log ( 1 x ) + ( t 1 ) log ( 1 + x ) = F ( x s ) + F ( x s ) ( x x s ) + 1 2 ! F ( x s ) ( x x s ) 2
neglecting higher-order terms. Imposing then that F ( x s ) = 0 , after approximating the gamma functions by their Stirling expressions, it is found that the beta distribution approaches the normal distribution
B ( s , t ) N t s s + t 2 , 4 ( t 1 ) ( s 1 ) ( t + s 2 ) 3 .
For normal distributions, by constructing the mean quantity
y = 1 + 1 k i = 1 k ( ± ) x i
in which x i is Gaussian distributed with an average of x ¯ i and variances of σ i , it is then found that y is distributed according to the normal distribution
N 1 + 1 k i = 1 k ( ± ) x ¯ i , 1 k 2 i = 1 k σ i 2 .

3. The Pseudo-Hermitian Hermite– β Ensemble

A matrix of the Hermite– β ensemble is a Hermitian tridiagonal matrix whose joint density distribution of the matrix elements is given by
P ( H β ) = 2 n 1 ( 2 π ) n / 2 j = 1 n 1 Γ [ ( N j ) β / 2 ] exp 1 2 tr H β 2 j = 1 n 1 b j ( n j ) β 1
such that the diagonal elements a i are normally distributed, namely N ( 0 , 1 ) , while the off-diagonal b j elements are distrtributed according to f ν ( y ) , as seen in Equation (11), where ν j = ( n j ) β and β is a real positive parameter. Using the lemmas of the tridiagonal form that relate the Vandermonde determinant to the off-diagonal elements, it is then found that the joint density distribution of the eigenvalues is given by [3]
P ( λ 1 , λ 2 , , λ n ) = C n exp 1 2 k = 1 n λ k 2 j > i | λ j λ i | β .
By expanding the determinant of the characteristic polynomial, it is obtained that its average satisfies the recursion relationship
< P n ( λ ) > = λ < P n 1 ( λ ) > n 1 2 β < P n 2 ( λ ) >
such that a comparison with Hermite polynomials recursion relationships leads to the identity < P n ( λ ) > = β 2 n 2 H n λ 2 β . This means that, asymptotically, the eigenvalues occupy the same compact support of the zeros of Hermite polynomials [18] defined by Wigner’s semi-circle law [2],
ρ W ( β , λ ) = 1 π β 2 n β λ 2 .
Considering the fluctuations, we assume that the nearest-neighbor distribution (NND) can be described by the NND of the 2 × 2 matrices that is found to be given by the generalized Wigner surmise
p W ( β , s ) = 2 κ β + 1 s β Γ ( β + 1 2 ) exp ( κ 2 s 2 ) ,
where κ = Γ ( β 2 + 1 ) / Γ ( β + 1 2 ) .
Let us now remove the Hermitian condition of H β by filling one of the off-diagonals by new random variables c j . In [10], this was carried out by sorting the elements of the two sub-diagonals from the same f ν ( y ) distribution of the Hermitian case. However, this choice does not keep the structure of the above joint distribution of elements in terms of the exponential of the trace. If we sort them instead from h ν ( y ) , using Equation (5), the joint distribution of the non-Hermitian matrices H ^ β can be written as
P ( H ^ β ) = 2 2 n 2 ( 2 π ) n / 2 j = 1 n 1 Γ 2 [ ( n j ) β / 2 ] exp 1 2 tr H ^ β 2 j = 1 n 1 b j c j ( n j ) β 1 .
In this case, the recursion relationship of the average characteristic polynomial becomes
< P ^ n ( λ ) > = λ < P ^ n 1 ( λ ) > ( n 1 ) β < P ^ n 2 ( λ ) > ,
which leads to the semi-circle law ρ W ( 2 β , λ ) , in which β is replaced by 2 β .
From the preliminary results stated above, H ^ β is a pseudo-Hermitian matrix whose eigenvalues are real, and, at the same time, they also are eigenvalues of a Hermitian matrix K β . The diagonal elements of K β are the same of H ^ β , while the off-diagonal elements are the geometrical mean d i = b i c i of the H ^ β elements and are therefore distributed according to Equation (6). Using, then, the asymptotic of these functions and Equation (7), the joint density distribution of the K β approaches the distribution
P ( K β ) 2 n 1 ( 2 π ) n / 2 j = 1 n 1 Γ [ ( n j ) β ] exp 1 2 tr K β 2 j = 1 n 1 , d j ( n j ) 2 β 1 ,
which is just the joint distribution of the Hermitian case with a β replaced by 2 β . The joint distribution of the eigenvalues is then given by Equation (18), with β replaced by 2 β . In Figure 1, this 2 β effect is illustrated in the case of the eigenvalues of the pseudo-Hermitian β –Hermite ensemble.

4. The Pseudo-Hermitian Laguerre– β Ensemble

Let B be the bidiagonal matrix with diagonal elements z = ( z n , z n 1 , , z 1 ) and sub-diagonal elements x = ( x n 1 , x n 2 , , x 1 ) , whose joint distribution of the elements is given by
P ( B β ) i = 1 n z n i + 1 2 α ( i 1 ) β 1 exp 1 2 z n i + 1 2 j = 1 n 1 x n j ( n j ) β 1 exp 1 2 x n j 2 ,
where α = m β / 2 with m n . Therefore, the elements obey the distribution h ν ( y ) , where ν = 2 α ( i 1 ) β for the diagonal elements and ν = ( n j ) β for the sub-diagonal elements. The Laguerre β ensemble is then defined as the tridiagonal matrix
L β = B B T = z n 2 z n x n 1 z n x n 1 z n 1 2 + x n 1 2 z n 1 x n 2 . . . . . . . . . z 3 x 2 z 2 2 + x 2 2 z 2 x 1 z 2 x 1 z 1 2 + x 1 2 .
As shown in Ref. [3], in order to change from a bidiagonal to a tridiagonal matrix, the Jacobian
J B L = 2 n z 1 i = 2 n z i 2 1
is used to give
P ( L β ) 2 n z 1 2 α ( n 1 ) β 2 exp 1 2 tr L β i = 1 n 1 z i + 1 2 α ( n 1 i ) β 3 j = 1 n 1 x j j β 1
such that, using lemmas of the tridiagonal form, it is found that the joint distribution of the eigenvalues is given by [3]
P ( λ 1 , λ 2 , , λ n ) = C n , m exp 1 2 i = 1 n λ i i = 1 n λ i ( m n + 1 ) β 2 1 i j | λ j λ i | β .
Let us start by obtaining the average characteristic polynomial whose derivation is also valid in the pseudo-Hermitian case. By expanding the determinant polynomial and taking the average of elements, we have the recursion relationship
< P n ( λ ) > = m β λ < R n 1 ( λ ) > m ( n 1 ) β 2 < R n 2 ( λ ) > ,
where the R-polynomials satisfy the recursion relationship
< R n 1 ( λ ) > = ( m + n 2 ) β λ < R n 2 ( λ ) > ( m 1 ) ( n 2 ) β 2 < R n 3 ( λ ) > .
These above relationships, respectively, are solved by the two generalized Laguerre polynomials < R n 1 ( λ ) > = ( n 1 ) ! β n 1 L n 1 m n + 1 ( λ β ) and < P n ( λ ) > = ( n ) ! β n L n m n ( λ β ) . The density of the eigenvalues, which is also that of the zeros of the polynomials, is given by the Marchenko–Pastur expression
ρ ( λ ) = 1 2 β π λ ( λ + λ ) ( λ λ ) ,
where λ ± = n β ( m n ± 1 ) 2 , and from this, we have the cumulative function
N ( λ ) = 1 4 π β 4 λ + λ arctan λ + ( λ λ ) λ ( λ + λ ) + ( λ + + λ ) arccos λ + + λ 2 λ λ + λ + 2 ( λ + λ ) ( λ λ ) .
To construct the pseudo-Hermitian–Laguerre ensemble, we introduce a new bidiagonal matrix C with a diagonal w = ( w n , w n 1 , , w 1 ) , and a lower sub-diagonal y = ( y n 1 , y n 2 , , y 1 ) , whose elements are sorted from the same distribution of the elements of B. Once this is carried out, a non-Hermitian β –Laguerre matrix is obtained by taking the product L ^ β = B C T that produces the matrix
L ^ β = z n w n z n y n 1 w n x n 1 z n 1 w n 1 + x n 1 y n 1 z n 1 y n 2 w m 1 x m 2 . . . . . . . . w 3 x 2 z 2 w 2 + x 2 y 2 z 2 y 1 w 2 x 1 z 1 w 1 + x 1 y 1
Immediately, the pseudo-Hermitian nature of this matrix is proved, and, furthermore, it can be seen that it is iso-spectral with the Hermitian matrix M β = D D T in which D is the bidiagonal matrix
D = z n w n x n 1 y n 1 z n 1 w n 1 x n 2 y n 2 z n 2 w n 2 . . . . . . . x 2 y 2 z 2 w 2 x 1 y 1 z 1 w 1 .
Denoting the diagonal elements as a i = z i w i and the sub-diagonals as b i = x i y i , of course, they are distributed according to g ν ( y ) , as per Equation (6), with ν = 2 α ( j 1 ) β for the diagonals and ν = ( n j ) β for the sub-diagonals. Replacing g ν ( y ) with its asymptotic expression, we obtain
P ( M β ) a 1 4 α 2 ( n 1 ) β 2 exp tr M β i = 1 n 1 a i + 1 4 α 2 ( n 1 j ) β 3 j = 1 n 1 b j 2 j β 1 ,
and for the eigenvalues, we obtain
P ( λ 1 , λ 2 , , λ n ) = C n , m exp i = 1 n λ i i = 1 n λ i ( m n + 1 ) β 1 i j | λ j λ i | 2 β .
In Figure 1, the 2 β effect is illustrated in the case of the pseudo-Hermitian β –Laguerre ensemble.

5. The Pseudo-Hermitian Jacobi– β Ensemble

A matrix of the Jacobi β ensemble is a Hermitian tridiagonal matrix [4,5,6,19]
J β = a 1 b 1 b 1 a 2 b 2 . . . . . . . . . b n 2 a n 1 b n 1 b n 1 a n
where the diagonal elements are given by
a k + 1 = ( 1 α 2 k 1 ) α 2 k ( 1 + α 2 k 1 ) α 2 k 2
and the off-diagonal elements are given by
b k + 1 = ( 1 α 2 k 1 ) ( 1 α 2 k 2 ) ( 1 + α 2 k + 1 ) ,
where, with 0 k 2 n 1 , the α k are β -distributed, following Equation (12), with variables given by
B 2 n k 2 4 β + a + 1 , 2 n k 2 4 β + b + 1 for   k   even
and
B 2 n k 3 4 β + a + b + 2 , 2 n k 1 4 β for   k   odd ,
aside from the special cases α 2 n 1 = α 1 = 1 . It has been proven that the joint distribution of the eigenvalues of the these matrices is given by
P ( λ 1 , λ 2 , , λ n ) = C n i = 1 n ( 2 λ i ) a 1 ( 2 + λ i ) b 1 j > i | λ j λ i | β .
The density of the eigenvalues satisfies an inverted semi-circle law
ρ ( λ ) = n π 4 λ 2 + c ,
where the constant c is a small parameter necessary to adjust the results of the numerical simulations. As a consequence, the unfolding of the eigenvalues is performed with the cumulative function
N ( λ ) = n π arcsin ( λ 2 ) + c λ .
As the above parameters of the beta distributions increase for large matrix sizes, it is justified to assume that when n > > 1 , they can be replaced by the Gaussian approximations
N 2 ( b a ) ( 2 n k 2 ) β + 2 a + 2 b , 2 [ ( 2 n k 2 ) β + 4 a ) ] [ ( 2 n k 2 ) β + 4 b ] [ ( 2 n k 2 ) β + 2 a + 2 b ] 3 for   k   even
and
N β 2 ( a + b ) 4 ( 2 n k 2 ) β + 2 a + 2 b , 2 [ ( 2 n k 3 ) β + 4 a + 4 b + 4 ) ] [ ( 2 n k 1 ) β 4 ] [ ( 2 n k 2 ) β + 2 a + 2 b ] 3 for   k   odd ,
which were derived in the preliminary section.
The important point to observe in these equations is that the average positions and the variances of the above Gaussians approach zero when the size n of the matrices becomes arbitrarily large. As a consequence, the above expressions connecting the matrix elements to the variables α k can be linearized as a k + 1 α 2 k α 2 k 2 for the diagonal elements and b k + 1 1 ( α 2 k 1 α 2 k + 1 ) / 2 for the off-diagonal ones. These approximated expressions show that, asymptotically, the diagonal and the off-diagonal elements decouple, and, moreover, all elements become Gaussian distributed, such that we have for the diagonal elements
P ( a k + 1 ) N α ¯ 2 k α ¯ 2 k 2 , σ 2 k 2 + σ 2 k 2 2 ,
and for the off-diagonal elements, we have
P ( b k + 1 ) N 1 α ¯ 2 k 1 α ¯ 2 k + 1 2 , σ 2 k 1 2 + σ 2 k + 1 2 4 .
Therefore, the off-diagonal elements fluctuate around one, while the diagonal ones fluctuate around zero. Thus, it is reasonable to expect that the off-diagonal elements have a determinant role in the behavior of the eigenvalue properties inasmuch as, for tridiagonal matrices, the eigenvalues are directly connected to the off-diagonal elements. The precision of this approximated description as compared to the exact expressions is attested to in Figure 2.
Therefore, the matrix elements are constructed using the 2n−1 random components of the vector α k , which are shared by both the diagonal and the sub-diagonal elements. Thus, it does not seem to have a clear way to remove the Hermitian condition. On the other hand, the basic idea of the 2 β Dyson effect is that it is caused by doubling the number of independent variables, especially those in the off-diagonal matrix elements. Considering now the removal of the Hermitian condition of the matrices, this can be executed by just sorting a new vector α ˜ k from which only the odd elements are used to fill the lower sub-diagonal. Therefore, only the number of odd elements are doubled. This defines a lower sub-diagonal whose elements c k + 1 differ from b k + 1 by odd terms. As the elements of these two off-diagonal are positive, these non-Hermitian matrices, besides being pseudo-Hermitian, are iso-spectral with the Hermitian matrix whose off-diagonal elements are given by b k + 1 c k + 1 .
Turning now to the asymptotic description, consistent with our level of approximations, the geometrical mean of the two off-diagonal elements can be linearized as
b k + 1 c k + 1 1 [ ( α 2 k 1 α 2 k + 1 ) + ( α 2 k 1 α 2 k + 1 ) ] / 4 ,
where the primes denote α variables belonging to the c k + 1 sub-diagonal. Using the result (16) from Section 2, we can write
P ( b k + 1 c k + 1 ) N 1 α ¯ 2 k 1 α ¯ 2 k + 1 2 , 1 2 σ 2 k 1 2 + σ 2 k + 1 2 4
such that the factor 1 / 2 affects the variances of the odd k terms as
2 [ ( 2 n k 3 ) 2 β + 8 ( a + b + 1 ) ] [ ( 2 n k 1 ) 2 β 8 ] [ ( 2 n k 2 ) 2 β + 2 ( a + b ) ] 3 .
The above expression shows the occurrence of a 2 β effect, which affects only the odd terms of the off-diagonal elements. This means that we are dealing with a hybrid pseudo-Hermitian ensemble that combines the even elements generated with β and the odd ones sorted with 2 β . The evidence from numerical simulations (see Figure 3) suggests that the spectral statistics, that is, the NNSD, can be fitted with a parameter 2 β . It is interesting to observe that 2 β is the geometric mean of β and 2 β . To ’fix’ this, what we can do is, in the pseudo-Hermitian case, generate the even elements with a parameter 2 β while sorting the odd elements with β . The results of the numerical simulations following this procedure demonstrate a clear 2 β effect, as is shown in Figure 4 for the NNSD and in Figure 5 for the number variance.

6. Conclusions

In conclusion, we have extended the universality of RMT Wigner–Dyson statistics to a class of non-Hermitian matrices with real eigenvalues. Matrices of this class are connected to their adjoints by a similarity transformation, a discrete symmetry condition that defines a pseudo-Hermitian matrix. In our case, the matrices are constructed by removing the Hermitian condition of the tridiagonal matrices of the three so-called β ensembles. By doing this, the number of off-diagonal independent elements doubles, and, as also occurs in standard random matrix theory, the Dyson β index also doubles. One obvious question is to investigate if this 2 β effect can be observed in pseudo-Hermitian full random matrices. However, it is not clear how, in this case, a pseudo-Hermitian model can be constructed in which the number of independent variables doubles. Another question is if the effect can be observed in a Hamiltonian with real eigenvalues, as this can lead to a link between PT symmetry and quantum chaos. As a final remark, we emphasize that the origin of the 2 β effect stems from the asymptotic behavior of the geometric mean calculated with the distribution functions of the matrix elements of the three ensembles. This suggests another universality associated with random matrix theories.

Author Contributions

Conceptualization, M.P.P.; Formal analysis, C.A.G.; Investigation, C.A.G.; Writing—review & editing, M.P.P. All authors have read and agreed to the published version of the manuscript.

Funding

M.P.P. is a member of the Brazilian National Institute of Science and Technology Quantum Information (INCT-IQ). C.A.G. is supported by the Brazilian agency CAPES under grant number 88887.616887/2021-00.

Institutional Review Board Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Baker, T.; Forrester, P. The Calogero-Sutherland Model and Generalized Classical Polynomials. Commun. Math. Phys. 1997, 188, 175. [Google Scholar] [CrossRef]
  2. Mehta, M.L. Random Matrices, 3rd ed.; Academic Press: London, UK, 2004. [Google Scholar]
  3. Dumitriu, I.; Edelman, A. Matrix models for beta ensembles. J. Math. Phys. 2002, 43, 5830. [Google Scholar] [CrossRef]
  4. Lippert, R.A. A matrix model for the β-Jacobi ensemble. J. Math. Phys. 2003, 44, 4807–4816. [Google Scholar] [CrossRef]
  5. Killip, T.; Nenciu, I. Matrix models for circular ensembles. Int. Math. Res. Not. 2004, 50, 2665. [Google Scholar] [CrossRef]
  6. Edelman, A.; Sutton, B. The Beta-Jacobi Matrix Model, the CS Decomposition, and Generalized Singular Value Problems. Found. Comput. Math. 2008, 8, 259–285. [Google Scholar] [CrossRef]
  7. Dyson, F.J. The Threefold Way. Algebraic Structure of Symmetry Groups and Ensembles in Quantum Mechanics. J. Math. Phys. 1962, 3, 1199. [Google Scholar] [CrossRef]
  8. Ginibre, J. Statistical Ensembles of Complex, Quaternion, and Real Matrices. J. Math. Phys. 1965, 6, 440. [Google Scholar] [CrossRef]
  9. Wigner, E.P. On the Distribution of the Roots of Certain Symmetric Matrices. Ann. Math. 1958, 68, 325. [Google Scholar] [CrossRef]
  10. Bohigas, O.; Pato, M.P. Non-Hermitian β-ensemble with real eigenvalues. AIP Adv. 2013, 3, 032130. [Google Scholar] [CrossRef]
  11. Mostafazadeh, A. Pseudo-Hermiticity versus PT symmetry: The necessary condition for the reality of the spectrum of a non-Hermitian Hamiltonian. J. Math. Phys. 2002, 43, 205–214. [Google Scholar] [CrossRef]
  12. Bender, C.M.; Boettcher, S. Real Spectra in Non-Hermitian Hamiltonians Having PT-Symmetry. Phys. Rev. Lett. 1998, 80, 5243. [Google Scholar] [CrossRef]
  13. Bender, C.M.; Boettcher, S.; Meisinger, P.N. PT-symmetric quantum mechanics. J. Math. Phys. 1999, 40, 2201. [Google Scholar] [CrossRef]
  14. Bender, C.M.; Brody, D.C.; Jones, H.F. Complex Extension of Quantum Mechanics. Phys. Rev. Lett. 2002, 89, 270401. [Google Scholar] [CrossRef] [PubMed]
  15. Bender, C.M. Making sense of non-Hermitian Hamiltonians. Rep. Prog. Phys. 2007, 70, 947. [Google Scholar] [CrossRef]
  16. Jain, S.R. Random matrix theories and exactly solvable models. Czechoslov. J. Phys. 2006, 56, 1021–1032. [Google Scholar] [CrossRef]
  17. Srivastava, S.C.L.; Jain, S.R. Pseudo-Hermitian random matrix theory. Fortschr. Phys. 2013, 61, 276. [Google Scholar] [CrossRef]
  18. Szego, G. Orthogonal Polynomials; American Mathematical Society: New York, NY, USA, 1939. [Google Scholar]
  19. Forrester, P.J. Log-Gases and Random Matrices; Princeton University Press: Princeton, NJ, USA, 2010. [Google Scholar]
Figure 1. NNSD for pseudo-Hermitian matrices of size n = 128 of the Hermite (blue circles) and Laguerre ensembles (red squares) for the indicated values of β = 0.5 and β = 2 . The black full lines were calculated with Wigner surmise and Equation (21), with β = 1 and β = 4 , respectively.
Figure 1. NNSD for pseudo-Hermitian matrices of size n = 128 of the Hermite (blue circles) and Laguerre ensembles (red squares) for the indicated values of β = 0.5 and β = 2 . The black full lines were calculated with Wigner surmise and Equation (21), with β = 1 and β = 4 , respectively.
Entropy 25 00868 g001
Figure 2. Comparison of the Gaussian approximation with the exact calculation of NNSD for β = 1 and β = 2 , as indicated in the figure.The calculations were performed with 50 matrices of size n = 450.
Figure 2. Comparison of the Gaussian approximation with the exact calculation of NNSD for β = 1 and β = 2 , as indicated in the figure.The calculations were performed with 50 matrices of size n = 450.
Entropy 25 00868 g002
Figure 3. The hybrid pseudo-Hermitian β –Jacobi ensemble in which even elements are sorted with β , while the odd ones, asymptotically, become sorted with 2 β . The blue line is calculated with Equation (21) using as parameters the geometric mean of β and 2 β , that is, 2 β . The calculations were performed with 50 matrices of size n = 360.
Figure 3. The hybrid pseudo-Hermitian β –Jacobi ensemble in which even elements are sorted with β , while the odd ones, asymptotically, become sorted with 2 β . The blue line is calculated with Equation (21) using as parameters the geometric mean of β and 2 β , that is, 2 β . The calculations were performed with 50 matrices of size n = 360.
Entropy 25 00868 g003
Figure 4. NNSD for 100 matrices of size 120 in which the even elements of the pseudo-Hermitian β –Jacobi ensemble are sorted with 2 β (black and blue dots) for the indicated values of beta. The full lines were calculated with Wigner surmise, as seen in Equation (21), with β = 1 and β = 4 , respectively.
Figure 4. NNSD for 100 matrices of size 120 in which the even elements of the pseudo-Hermitian β –Jacobi ensemble are sorted with 2 β (black and blue dots) for the indicated values of beta. The full lines were calculated with Wigner surmise, as seen in Equation (21), with β = 1 and β = 4 , respectively.
Entropy 25 00868 g004
Figure 5. The verification of the 2 β effect in the number variance statistics: black dots were calculated with 600 matrices of size n = 120 of the pseudo-Hermitian β –Jacobi ensemble sorted with β = 0.5 , and the red full line corresponds to the GOE number variance.
Figure 5. The verification of the 2 β effect in the number variance statistics: black dots were calculated with 600 matrices of size n = 120 of the pseudo-Hermitian β –Jacobi ensemble sorted with β = 0.5 , and the red full line corresponds to the GOE number variance.
Entropy 25 00868 g005
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Goulart, C.A.; Pato, M.P. The Double Dyson Index β Effect in Non-Hermitian Tridiagonal Matrices. Entropy 2023, 25, 868. https://doi.org/10.3390/e25060868

AMA Style

Goulart CA, Pato MP. The Double Dyson Index β Effect in Non-Hermitian Tridiagonal Matrices. Entropy. 2023; 25(6):868. https://doi.org/10.3390/e25060868

Chicago/Turabian Style

Goulart, Cleverson A., and Mauricio P. Pato. 2023. "The Double Dyson Index β Effect in Non-Hermitian Tridiagonal Matrices" Entropy 25, no. 6: 868. https://doi.org/10.3390/e25060868

APA Style

Goulart, C. A., & Pato, M. P. (2023). The Double Dyson Index β Effect in Non-Hermitian Tridiagonal Matrices. Entropy, 25(6), 868. https://doi.org/10.3390/e25060868

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop