Random matrix ensembles defined by a mean-field one body plus a chaos generating random two-body interaction [called embedded ensembles of (1+2)-body interactions] predict for wave functions, in the chaotic domain, an essentially one-parameter Gaussian forms for the energy dependence of the number of principal components (NPC) and the localization length l(H) (defined by information entropy), which are two important measures of chaos in finite interacting many-particle systems. Numerical embedded ensemble calculations and nuclear shell-model results, for NPC and l(H), are compared with the theory. These analyses clearly point out that for realistic finite interacting many-particle systems, in the chaotic domain, wave-function structure is given by (1+2)-body embedded random matrix ensembles.
Finite interacting Fermi systems with a mean-field and a chaos generating two-body interaction are modeled by one plus two-body embedded Gaussian orthogonal ensemble of random matrices with spin degree of freedom [called EGOE(1+2)-s]. Numerical calculations are used to demonstrate that, as lambda , the strength of the interaction (measured in the units of the average spacing of the single-particle levels defining the mean-field), increases, generically there is Poisson to GOE transition in level fluctuations, Breit-Wigner to Gaussian transition in strength functions (also called local density of states) and also a duality region where information entropy will be the same in both the mean-field and interaction defined basis. Spin dependence of the transition points lambda_{c} , lambdaF, and lambdad , respectively, is described using the propagator for the spectral variances and the formula for the propagator is derived. We further establish that the duality region corresponds to a region of thermalization. For this purpose we compared the single-particle entropy defined by the occupancies of the single-particle orbitals with thermodynamic entropy and information entropy for various lambda values and they are very close to each other at lambda=lambdad.
Embedded random matrix ensembles generated by random interactions (of low body rank and usually two-body) in the presence of a one-body mean field, introduced in nuclear structure physics, are now established to be indispensable in describing statistical properties of a large number of isolated finite quantum many-particle systems. Lie algebra symmetries of the interactions, as identified from nuclear shell model and the interacting boson model, led to the introduction of a variety of embedded ensembles (EEs). These ensembles with a mean field and chaos generating two-body interaction generate in three different stages, delocalization of wave functions in the Fock space of the mean-field basis states. The last stage corresponds to what one may call thermalization and complex nuclei, as seen from many shell model calculations, lie in this region. Besides briefly describing them, their recent applications to nuclear structure are presented and they are (i) nuclear level densities with interactions; (ii) orbit occupancies; (iii) neutrinoless double beta decay nuclear transition matrix elements as transition strengths. In addition, their applications are also presented briefly that go beyond nuclear structure and they are (i) fidelity, decoherence, entanglement and thermalization in isolated finite quantum systems with interactions; (ii) quantum transport in disordered networks connected by many-body interactions with centrosymmetry; (iii) semicircle to Gaussian transition in eigenvalue densities with [Formula: see text]-body random interactions and its relation to the Sachdev–Ye–Kitaev (SYK) model for majorana fermions.
Employing hamiltonians defined by two-body embedded Gaussian orthogonal ensemble of random matrices(EGOE(2)) plus a mean-field producing one-body part, strength functions (for states defined by the one-body part) are constructed for various values of the strength of the chaos generating two-body part. Numerical calculations for six and seven fermion systems clearly demonstrate Breit-Wigner to Gaussian transition, in the chaotic domain, in strength functions as found earlier in nuclear shell model and Lipkin-Meshkov-Glick model calculations.
There is a newly emerging understanding that in the chaotic domain of isolated finite interacting many particle systems smoothed densities define the statistical description of these systems and these densities follow from embedded (two-body) random matrix ensembles and their various deformations. These ensembles predict that the smoothed form of matrix elements of a transition operator between the chaotic eigenstates weighted by the densities at the two ends (i.e. the bivariate strength density) will be a bivariate Gaussian with the bivariate correlation coefficient arising out of the non-commutability of the hamiltonian and the transition operator involved. The ensemble theory extends to systems with a mean-field and a chaos generating two-body interaction (as in nuclei, atoms and diffusive quantum dots). These developments in many-body quantum chaos are described with special reference to one-body transition operators.
Embedded random matrix ensembles operating in nuclear shell model spaces, with nucleons occupying a finite set of single particle orbits and interacting via a two-body interaction, form the basis for statistical shell model. With sufficiently strong interaction, the level densities in shell model spaces take close to a Gaussian form and transition strength distributions close to a bivariate Gaussian form. In practice, partitioning via spherical configurations (
m~) and angular momentum
J (also isospin where appropriate) are essential. The resulting statistical spectroscopy or statistical shell model was applied successfully in the past in some studies of nuclear level densities, orbit occupancies,
β-decay matrix elements and so on. Going beyond these, recently it is recognized that embedded ensembles, in a better approximation, generate in-fact
q-normal form (
q=1 gives Gaussian and
q=0 Wigner's semi-circle) for density of eigenvalues, bivariate
q-normal form for transition strengths and conditional
q-normal form for strength functions. These then allow us to develop statistical shell model with
q-normal forms. These new developments in embedded ensembles and statistical shell model are briefly reviewed in this paper. Also described, using some examples, is the role of the
q parameter in generating statistical properties of general quantum many-particle systems.