Random matrix theory for transition strength densities in finite quantum systems: Results from embedded unitary ensembles

Embedded random matrix ensembles are generic models for describing statistical properties of finite isolated interacting quantum many-particle systems. For the simplest spinless fermion (or boson) systems, with say m$m$ fermions (or bosons) in N$N$ single particle states and interacting via k$k$-body interactions, we have EGUE(k$k$) embedded GUE of k$k$-body interactions] with GUE embedding and the embedding algebra is U(N)$U\left(N\right)$. A finite quantum system, induced by a transition operator, makes transitions from its states to the states of the same system or to those of another system. Examples are electromagnetic transitions (then the initial and final systems are same), nuclear beta and double beta decay (then the initial and final systems are different), particle addition to or removal from a given system and so on. Towards developing a complete statistical theory for transition strength densities (transition strengths multiplied by the density of states at the initial and final energies), we have derived formulas for the lower order bivariate moments of the strength densities generated by a variety of transition operators. Firstly, for a spinless fermion system, using EGUE(k$k$) representation for a Hamiltonian that is k$k$-body and an independent EGUE(t$t$) representation for a transition operator that is t$t$-body and employing the embedding U(N)$U\left(N\right)$ algebra, finite-N$N$ formulas for moments up to order four are derived, for the first time, for the transition strength densities. Secondly, formulas for the moments up to order four are also derived for systems with two types of spinless fermions and a transition operator similar to beta decay and neutrinoless beta decay operators. In addition, moments formulas are also derived for a transition operator that removes k₀$k_0$ number of particles from a system of m$m$ spinless fermions. In the dilute limit, these formulas are shown to reduce to those for the EGOE version derived using the asymptotic limit theory of Mon and French (1975). Numerical results obtained using the exact formulas for two-body (k=2$k=2$) Hamiltonians (in some examples for k=3$k=3$ and 4$4$) and the asymptotic formulas clearly establish that in general the smoothed (with respect to energy) form of the bivariate transition strength densities take bivariate Gaussian form for isolated finite quantum systems. Extensions of these results to bosonic systems and EGUE ensembles with further symmetries are discussed.