  A theoretical framework is presented for the computation of rovibrational
polaritonic states of a molecule in a lossless infrared (IR) microcavity. In
the proposed approach the quantum treatment of the rotational and vibrational
motion of the molecule can be formulated using arbitrary approximations. The
cavity-induced changes in electronic structure are treated perturbatively,
which allows using the existing polished tools of standard quantum chemistry
for determining electronic molecular properties. As a case study, the
rovibrational polaritons and related thermodynamic properties of H$_2$O in an
IR microcavity are computed for varying cavity parameters and applying various
approximations to describe the molecular degrees of freedom. The self-dipole
interaction is found to be significant for nearly all light-matter coupling
strengths investigated, and the molecular polarizability proved to be important
for the correct qualitative behavior of the energy level shifts induced by the
cavity. On the other hand, the magnitude of polarization remains small,
justifying the perturbative approach for the cavity-induced changes in
electronic structure. Comparing results obtained using a high-accuracy
variational molecular model with those obtained utilizing the rigid rotor and
harmonic oscillator approximations revealed that as long as the rovibrational
model is appropriate for describing the field-free molecule, the computed
rovibropolaritonic properties can be expected to be accurate as well. Strong
light-matter coupling between the radiation mode of an IR cavity and the
rovibrational states of H$_2$O lead to minor changes in the thermodynamic
properties of the system, and these changes seem to be dominated by
non-resonant interactions between the quantum light and matter.
