Reproduction of exact solutions of Lipkin model by nonlinear higher random-phase approximation

It is shown that the random-phase approximation (RPA) method with its nonlinear higher generalization, which was previously considered as approximation except for a very limited case, reproduces the exact solutions of the Lipkin model. The nonlinear higher RPA is based on an equation nonlinear on eigenvectors and includes many particle-many hole components in the creation operator of the excited states. We demonstrate the exact character of solutions analytically for the particle number [Formula: see text] and, numerically, for [Formula: see text]. This finding indicates that the nonlinear higher RPA is equivalent to the exact Schrödinger equation, which opens up new possibilities for realistic calculations in many-body problems.