Virasoro conformal blocks are universal ingredients of correlation functions of two-dimensional conformal field theories (2d CFTs) with Virasoro symmetry. It is acknowledged that in the (classical) limit of large central charge of the Virasoro algebra and large external, and intermediate conformal weights with fixed ratios of these parameters Virasoro blocks exponentiate to functions known as Zamolodchikovs’ classical blocks. The latter are special functions which have awesome mathematical and physical applications. Uniformization, monodromy problems, black holes physics, quantum gravity, entanglement, quantum chaos, holography, $$ \mathcal{N} $$
N
= 2 gauge theory and quantum integrable systems (QIS) are just some of contexts, where classical Virasoro blocks are in use. In this paper, exploiting known connections between power series and integral representations of (quantum) Virasoro blocks, we propose new finite closed formulae for certain multi-point classical Virasoro blocks on the sphere. Indeed, combining classical limit of Virasoro blocks expansions with a saddle point asymptotics of Dotsenko-Fateev (DF) integrals one can relate classical Virasoro blocks with a critical value of the “Dotsenko-Fateev matrix model action”. The latter is the “DF action” evaluated on a solution of saddle point equations which take the form of Bethe equations for certain QIS (Gaudin spin models). A link with integrable models is our main motivation for this research line. For instance, analogous quantities as the “DF on-shell action” appear in 2d CFT realization of the so-called Richardson’s solution of the reduced BCS model describing physics of ultra-small superconducting grains. Precisely, the Richardson solution and its particular limit characterizing the rational Gaudin model have known implementation in 2d CFT in terms of a free field representation of certain (perturbed and unperturbed) WZW blocks. The WZW analogue of the “DF on-shell action” appears here and plays a crucial role, it is a generating function for eigenvalues of quantum integrals of motion. An exploration of relationships between power expansions and Coulomb gas representations of the aforementioned WZW blocks could pave a way for new analytical tools useful in the study of models of the Richardson-Gaudin type.