For
m number of bosons, carrying spin ($\cs=\spin$) degree of freedom, in
Ω number of single particle orbitals, each doubly degenerate, we introduce and analyze embedded Gaussian orthogonal ensemble of random matrices generated by random two-body interactions that are spin (
S) scalar [BEGOE(2)-$\cs$]. Embedding algebra for the BEGOE(2)-$\cs$ ensemble and also for BEGOE(1+2)-$\cs$ that includes the mean-field one-body part is
U(2Ω)⊃U(Ω)⊗SU(2) with SU(2) generating spin. A method for constructing the ensembles in fixed-(
m,S) spaces has been developed. Numerical calculations show that for BEGOE(2)-$\cs$, the fixed-
(m,S) density of states is close to Gaussian and level fluctuations follow GOE in the dense limit. For BEGOE(1+2)-$\cs$, generically there is Poisson to GOE transition in level fluctuations as the interaction strength (measured in the units of the average spacing of the single particle levels defining the mean-field) is increased. The interaction strength needed for the onset of the transition is found to decrease with increasing
S. Covariances in energy centroids and spectral variances are analyzed. Propagation formula is derived for the variance propagator for the fixed-
(m,S) ensemble averaged spectral variances. Variance propagator clearly shows, by applying the Jacquod and Stone prescription, that the BEGOE(2)-$\cs$ ensemble generates ground states with spin
S=Smax. This is further corroborated by analyzing the structure of the ground states in the presence of the exchange interaction
S^2 in BEGOE(1+2)-$\cs$. Natural spin ordering is also observed with random interactions. Going beyond these, we also introduce pairing symmetry in the space defined by BEGOE(2)-$\cs$. Expectation values of the pairing Hamiltonian show that random interactions exhibit pairing correlations in the ground state region.
In this paper we discuss in detail the
P(I)′s, angular momentum I probabilities in the ground states, of many-body systems interacting via a two-body random ensemble (TBRE). In particular, we extensively apply an approach introduced in an earlier paper and compare the predicted
P(I)′s with those obtained by diagonalizing a TBRE Hamiltonian. We begin with a few solvable cases, such as fermions in a small single-
j shell and d boson systems, where elegant agreements between the predicted
P(I)′s and those obtained by diagonalizing a TBRE Hamiltonian are achieved. We find that d boson systems systematically present counterexamples of angular momentum 0 ground state dominance when the number of d bosons is $6\ensuremath{\kappa}\ifmmode\pm\else\textpm\fi{}1$ with $\ensuremath{\kappa}$ a natural number, which suggests that certain fundamental symmetry (say, time reversal invariance) of the Hamiltonian cannot ensure the occurrence of angular momentum 0 ground state dominance. Next, we apply the same approach to more complicated cases, such as even or odd number of fermions in a large single-
j shell or a many-
j shell,
sd-boson or
sdg-boson systems, etc. We find that the simple approach proposed in an earlier paper is also well applicable, and thus it is a universal approach. The numerical experiments provide a guideline to tell which interactions are essential to produce a sizable
P(I) in a many-body system. This disproves a popular idea that the angular momentum 0 ground state (0 g.s.) dominance may be independent of two-body interactions. Some matrix elements which are useful to understand the observed regularities are given or addressed in detail. In this paper we also report a synchronous staggering between the 0 g.s. probabilities of even numbers of fermions in a single-
j shells and j g.s. probabilities of odd numbers of fermions in a single-
j shell when j is small. The low seniority chain of 0 g.s. using the same set of two-body interactions is confirmed, but it is noted that contribution to the total 0 g.s. probability beyond this chain may be more important for even numbers of fermions in a single-
j shell. Some interesting results taking a displaced two-body random ensemble are presented for the I ground state probabilities.
An ensemble with random n-body interactions is investigated in the presence of symmetries. A striking emergence of regularities in the spectra, ground state spins, and isospins is discovered. In even-particle systems the strong preponderance of spin zero ground states is observed with three, four, and higher-order random many-body interactions, and similar trends are discovered in systems that possess both rotational and isospin symmetry. Various types of correlations from pairing to spectral sequences and correlations across different masses are explored. A search for an interpretation is presented.
Random matrix ensembles for a system of m number of fermions or bosons in Ω number of single particle levels each r-fold degenerate and interacting with two-body forces are considered. The spectrum generating algebra for these systems is U(rΩ) and a subalgebra of interest is U(rΩ) ⊃ U(Ω) ⊗ SU(r) algebra. Now, for random two-body interactions preserving SU(r) symmetry, one can introduce embedded Gaussian unitary ensemble of random matrices with U(Ω)⊗SU(r) embedding and this class of ensembles are denoted by EGUE(2)-SU(r). Ensembles with r = 1,2 and 4 for fermions correspond to spinless fermions, fermions with spin and fermions with Wigner's spin-isospin SU(4) symmetry respectively. Similarly, for bosons r = 1, 2 and 3 correspond to spinless bosons, two species boson systems and bosons with spin one respectively. The distinction between fermions and bosons is in the U(Ω) irreducible representations. General formulation based on Wigner-Racah algebra for lower order moments of the one- and two-point functions in eigenvalues generated by EGUE(2)- SU(r) is briefly reviewed. The final formulas for the moments involve only SU(Ω) Racah coefficients. For the fourth moment of the one-point function for r > 1 and for the higher order (> 4) bivariate moments of the two-point function for r ≥ 1, formulas are not available for the SU(Ω) Racah coefficients that are needed. It is necessary to derive analytical formulas for these or develop methods that give asymptotic results (an example for this is given in the paper) or develop methods that allow for their numerical evaluation. This important open problem is discussed in some detail.
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.
Configuration interaction calculations for the Sm I
J=0+ and
4+ with matrix dimensions 1351 and 7325 are carried out. The eigenvalues and eigenfunctions are analyzed for density of states, strength functions, and information entropy in wave functions and compared with the Gaussian forms given by two-body random matrix ensembles (TBRE). Signatures of TBRE are clearly seen in the Sm I example.