We introduce and solve exactly a family of invariant 2
xd7;2 random matrices, depending on one parameter
η, and we show that rotational invariance and real Dyson index
β are not incompatible properties. The probability density for the entries contains a weight function and a multiple trace–trace interaction term, which corresponds to the representation of the Vandermonde-squared coupling on the basis of power sums. As a result, the effective Dyson index
![View the MathML source View the MathML source](http://www.sciencedirect.com/cache/MiamiImageURL/B6TVG-4S1C2NP-1-3D/0?wchp=dGLbVlW-zSkWz)
of the ensemble can take any real value in an interval. Two weight functions (Gaussian and non-Gaussian) are explored in detail and the connections with
β-ensembles of Dumitriu–Edelman and the so-called Poisson–Wigner crossover for the level spacing are respectively highlighted. A curious spectral twinning between ensembles of different symmetry classes is unveiled: as a consequence, the identification between symmetry group (orthogonal, unitary or symplectic) and the exponent of the Vandermonde (
β=1,2,4) is shown to be potentially deceptive. The proposed technical tool more generically allows for designing actual matrix models which (i) are rotationally invariant; (ii) have a real Dyson index
c4">![View the MathML source View the MathML source](http://www.sciencedirect.com/cache/MiamiImageURL/B6TVG-4S1C2NP-1-CR/0?wchp=dGLbVlW-zSkWz)
; (iii) have a pre-assigned confining potential or alternatively level-spacing profile. The analytical results have been checked through numerical simulations with an excellent agreement. Eventually, we discuss possible generalizations and further directions of research.