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.
The role of the small-norm amplitudes in extended random-phase approximation (RPA) theories such as the particle-particle and hole-hole components of one-body amplitudes and the twobody amplitudes other than two-particle/two-hole components are investigated for the onedimensional Hubbard model using an extended RPA derived from the time-dependent density matrix theory.It is found that these amplitudes cannot be neglected in strongly interacting regions where the effects of ground-state correlations are significant.............
A general theory of long range correlations in the ground state of a closed shell nucleus is outlined, using a number of ideas introduced or adumbrated, but not always successfully developed, in previous work of the authors and others. These include the following: (i) Any operator can be formally expressed as a series in particle-hole (p-h) excitation and destruction operators which should converge rapidly for a suitable choice of single-particle basis at the closed shells. By means of these relations and with the help of sum rules all observables can be computed without the explicit use of wave functions. (ii) For a fixed choice of single-particle basis, the ground state energy can be computed from a variational principle in which in lowest nontrivial approximation, in consequence of (i), the parameters are amplitudes of the same type as occur in the random phase approximation (RPA). The precise relation to RPA is fully detailed. (iii) The single particle basis chosen is the natural orbital basis of Lowdin, which presents itself as the obvious choice both for formulation of the dynamics and for the algorithm of solution. (iv) Though the theory deals with pairs of fermion operators, particular attention is paid that all Pauli principle restrictions are satisfied. (v) It is verified that in the perturbation limit, known results are reproduced. The method is generally applicable to any system for which the correlations can be viewed as particle-hole excitations with respect to a reference Slater determinant.
An extended Lipkin-Meshkov-Glick model for testing the proton-neutron random-phase approximation
(pnRPA) method is developed, taking into account explicitly proton and neutron degrees of freedom. Besides the proton and neutron single-particle terms two types of residual proton-neutron interactions, one simulating a particle-particle and the other a particle-hole interaction, are included in the model Hamiltonian so that the model is exactly solvable in an isospin $\mathrm{SU}(2)\ensuremath{\bigotimes}\mathrm{SU}(2)$ basis. The behavior of the first excited (collective) state obtained by (i) exact diagonalization of the Hamiltonian matrix and (ii) with the
pnRPA is studied as a function of the model parameters and the two results are compared with each other. Furthermore, charge-changing operators simulating nuclear beta decay and their action on eigenfunctions of the model Hamiltonian are defined and transition amplitudes of them are calculated using exact, the Tamm-Dancoff, and
pnRPA eigenfunctions.
In this second work in a series of papers the coupled-cluster or
exp(S) formalism is applied to the problem of ground-state correlations in the one-component Fermi plasma. The main approximation employed is the so-called SUB2 approximation for the two-body subsystem correlation operator
S2 which provides a measure of the two-particle---two-hole component in the true ground-state wave function, and in which the coupling to larger subsystems is neglected. The SUB2 equations for Fermi systems are brought into tractable form by a "state-averaging" procedure which we develop by analogy with the comparable Bose equations, and which is shown by comparison with earlier exact results to be accurate in the metallic density regime at about the 1% level. Our final results for metallic densities include the completely integrated and self-consistent effects of the terms which by themselves generate (i) the random-phase approximation (RPA) and its associated exact long-range screening effects, (ii) the extra random-phase approximation exchange terms necessary to keep the RPA explicitly antisymmetric, (iii) the self-consistent particle-particle ladders (LAD) that describe two-particle scattering within the many-body medium and which describe the exact short-range behavior, (iv) a class of particle-hole ladder terms, and (v) the self-consistent hole-potential terms. Particular attention is paid to the important effects caused by the interference at intermediate separations of the long-range RPA and the short-range LAD effects. By comparison with recent and essentially exact stochastic simulations of the many-body Schr\"odinger equation, our results are seen to be accurate to about 1% for metallic densities, and hence to provide what is probably the currently best microscopic description available for this system. We show also that even in the low-density limit the SUB2 approximation provides a good "translationally-invariant-solid" description of both charged Fermi and Bose systems in this exact Wigner crystal regime.
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.