One plus two-body random matrix ensembles with parity: Density of states and parity ratios

TL;DR

This paper introduces a random matrix ensemble with parity for fermion systems, demonstrating Gaussian-like state densities, parity ratios similar to nuclear models, and predicting ground state parity preferences.

Contribution

It generalizes the EGOE(1+2) ensemble to include parity, providing analytical formulas and numerical evidence for state densities and parity ratios in nuclear-like systems.

Findings

Gaussian form of fixed parity state densities

Parity ratios resemble Fermi-gas model predictions

Preponderance of positive parity ground states at small mixing

Abstract

One plus two-body embedded Gaussian orthogonal ensemble of random matrices with parity [EGOE(1+2)-$\pi$] generated by a random two-body interaction (modeled by GOE in two particle spaces) in the presence of a mean-field, for spinless identical fermion systems, is defined, generalizing the two-body ensemble with parity analyzed by Papenbrock and Weidenm\"{u}ller [Phys. Rev. C {\bf 78}, 054305 (2008)], in terms of two mixing parameters and a gap between the positive $(\pi=+)$ and negative $(\pi=-)$ parity single particle (sp) states. Numerical calculations are used to demonstrate, using realistic values of the mixing parameters appropriate for some nuclei, that the EGOE(1+2)-$\pi$ ensemble generates Gaussian form (with corrections) for fixed parity eigenvalue densities (i.e. state densities). The random matrix model also generates many features in parity ratios of state densities that are similar to those predicted by a method based on the Fermi-gas model for nuclei. We have also obtained, by applying the formulation due to Chang et al [Ann. Phys. (N.Y.) {\bf 66}, 137 (1971)], a simple formula for the spectral variances defined over fixed-$(m_1,m_2)$ spaces, where $m_1$ is the number of fermions in the +ve parity sp states and $m_2$ is the number of fermions in the -ve parity sp states. Similarly, using the binary correlation approximation, in the dilute limit, we have derived expressions for the lowest two shape parameters. The smoothed densities generated by the sum of fixed-$(m_1,m_2)$ Gaussians with lowest two shape corrections describe the numerical results in many situations. The model also generates preponderance of +ve parity ground states for small values of the mixing parameters and this is a feature seen in nuclear shell model results.