Abstract A calorimetric and thermodynamic study of the (Li1-xNax)FeGe2O6 solid solution with the pyroxene structure was undertaken. The molar heat capacity at constant pressure (Cp,m) for compositions with x = (0.0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1.0) was measured using a Physical Properties Measurement System at temperatures from 2 K to 300 K and by differential scanning calorimetry between 282 and 870 K (these measurements were performed upon heating as well as upon cooling). Magnetic transitions below 20 K and structural P21/c ↔ C2/c phase transitions at higher temperatures were observed that decrease in temperature with increasing Na/(Na + Li). The P21/c ↔ C2/c transition showed a considerable temperature hysteresis that is largest for LiFeGe2O6 (∼30 K) and decreases with increasing Na/(Na + Li) in the pyroxene. The molar enthalpy and entropy changes associated with the P21/c ↔ C2/c phase transitions, ΔHtr,m and ΔStr,m, were determined. They can be described by quadratic and linear functions across the solid solution. A linear correlation between ΔHtr,m and the square of the volume strain at the transition temperature was established. The P21/c ↔ C2/c phase transition of LiFeGe2O6 is compared to that of α-β quartz, whose Cp,m was also measured. The molar enthalpy of order, ΔHord,m, was calculated from the calorimetric data. It has a maximum of ΔHord,m = (−11.6 ± 1.5)kJ·mol−1 for the P21/c ↔ C2/c phase transition of LiFeGe2O6, which is in reasonable agreement with a maximum ΔHord,m = −9.9 kJ·mol−1 calculated using density functional theory. ΔHord,m vs. T behavior was described using a tricritical Landau model with an a parameter of a = (41.8 ± 1.0)J·K−1·mol−1. All properties associated with the P21/c ↔ C2/c phase transition in LiFeGe2O6 can be calculated using this value. The Cp,m behavior at T 0.4 have the P21/c ↔ C2/c phase transition below 298.15 K, this results in a discontinuity in Scal,m298.15 behavior as function of composition at Na/(Na + Li) = 0.5. DQF was applied to model this behavior and the corresponding parameters were retrieved. The entropy of mixing behavior of the (Li,Na)FeGe2O6 solid solution at 298.15 K is characterized by positive excess molar magnetic entropies of mixing onto which entropy changes due to the structural P21/c ↔ C2/c phase transitions are superimposed in the Na-rich part of the binary system. Excess volumes of mixing for the (Li,Na)FeGe2O6 solid solution were calculated from published data and DQF parameters were derived that allow volume as function of temperature and composition to be calculated. The molar third law entropies at T = 298.15 K are Som = (188.3 ± 2)J·K−1·mol−1 and Som = (203.1 ± 2)J·K−1·mol−1 and the molar enthalpies of formation from the elements at T = 298.15 K are ΔfHom = (−1817.2 ± 9)kJ·mol−1 and ΔfHom = (−1834.0 ± 9)kJ·mol−1 for the end members LiFeGe2O6 and NaFeGe2O6, respectively. ΔfHom was calculated using density functional theory (DFT). Values for ΔfHom, Som, and molar volume, for the coefficients of Cp,m polynomials, for thermal expansion and for the bulk modulus of the end members were compiled and can be used, in combination with the derived volumetric and entropic mixing properties, for phase diagram calculations involving the (Li,Na)FeGe2O6 solid solution.