The thermodynamic properties of the Robin quantum well with extrapolation length are analyzed theoretically both for the canonical and two grand canonical ensembles, with special attention being paid to the situation when the energies of one or two lowest-lying states are split off from the rest of the spectrum by the large gap that is controlled by the varying . For the single split-off level, which exists for the geometry with the equal magnitudes but opposite signs of the Robin distances on the confining interfaces, the heat capacity c V of the canonical averaging is a nonmonotonic function of the temperature T with its salient maximum growing to infinity as ln 2 for the decreasing to zero extrapolation length and its position being proportional to 1/( 2 ln ). The specific heat per particle c N of the Fermi-Dirac ensemble depends nonmonotonically on the temperature too, with its pronounced extremum being foregone on the T axis by the plateau, whose value at the dying is (N − 1)/(2N)k B , with N being the number of fermions. The maximum of c N , similar to the canonical averaging, unrestrictedly increases as goes to zero and is largest for one particle. The most essential property of the Bose-Einstein ensemble is the formation, for a growing number of bosons, of the sharp asymmetric shape on the c N −T characteristics, which is more protrusive at the smaller Robin distances. This cusp-like structure is a manifestation of the phase transition to the condensate state. For two split-off orbitals, one additional maximum emerges whose position is shifted to colder temperatures with the increase of the energy gap between these two states and their higher-lying counterparts and whose magnitude approaches a -independent value. All these physical phenomena are qualitatively and quantitatively explained by the variation of the energy spectrum by the Robin distance. Parallels with other structures are drawn and similarities and differences between them are highlighted. Generalization to higher dimensions is also provided.