Rhabdophane is a hydrous phosphate that commonly replaces monazite as a weathering product in critical mineral deposits during the alteration of rare earth elements (REE) bearing carbonatites and alkaline igneous complexes. It is an important host to the light (L)REE (i.e., La to Gd) but the stability and structure of binary solid solutions between the Ce and the other LREE endmembers have not yet been determined experimentally. Here we present room temperature calorimetric experiments that were used to measure the enthalpy of precipitation of rhabdophane (Ce_1−xREE_xPO₄·nH₂O; REE = La, Pr, Nd, Sm, Eu, and Gd). The solids were characterized using X-ray diffraction, scanning electron microscopy, Raman spectroscopy, and the role of water in the rhabdophane structure was further determined using thermogravimetric analysis coupled with differential scanning calorimetry. The calorimetric experiments indicate a non-ideal behavior for all of the binary solid solutions investigated with an excess enthalpy of mixing (ΔH^ex) described by a 2- to 3-term Guggenheim parameters equation. The solid solutions were categorized into three groups: (1) binary Ce-La and Ce-Pr which display positive ΔH^ex values with a slight asymmetry; (2) binary Ce-Nd and Ce-Sm which display negative ΔH^ex values with a nearly symmetric shape; (3) Ce-Eu and Ce-Gd which display both negative and positive ΔH^ex values with nearly symmetric shape. The excess Gibbs energy (ΔG^ex) of the solid solutions was further investigated using a thermodynamic analysis approach of aqueous-solid solution equilibria and the optimization programs GEMS and GEMSFITS. The resulting ΔG^ex values combined with the calorimetric ΔH^ex values indicate that there is likely an excess entropy contribution implying important short-range structural modifications in the solid solutions dependent on the deviation of the REE ionic radii from the size of Ce³⁺. These observations corroborate with the trends in the Raman v₁ stretching bands of the PO₄-site. The excess molar volumes determined from X-ray diffraction analysis further indicate an overall asymmetric behavior in all of the studied binary solid solutions, which becomes increasingly important from La to Gd. The pronounced short-range order–disorder occurring in groups 2 and 3 solid solutions mimics some of the behavior observed from previous studies in anhydrous monazite solid solutions. This study highlights the potential to use the chemistry and the structural modifications of rhabdophane as potential indicators of formation conditions in geologic systems and permits improving our modeling capabilities of REE partitioning in critical minerals systems.