A Landau-Devonshire thermodynamic energy density function for ferroelectric wurtzite aluminum scandium nitride (Al1−xScxN) solid solution is developed. It is parametrized using available experimental and theoretical data, enabling the accurate reproduction of composition-dependent ferroelectric properties, such as spontaneous polarization, dielectric permittivity, and piezoelectric constants, for both bulk and thin films. The maximum concentration of Sc for the wurtzite structure to remain ferroelectric is found to be 61 at. %. A detailed analysis of Al1−xScxN thin films reveals that the ferroelectric phase transition and properties are insensitive to substrate strain. This study lays the foundation for quantitative modeling of novel ferroelectric wurtzite solid solutions.