The ClaySor 2023 model package within the GEM-Selektor software includes an updated version of the two-site protolysis non-electrostatic surface complexation and cation exchange sorption model for illite and montmorillonite, as well as the first implementation of the generalised Caesium sorption model for illite. These models have been harmonised with the most recent PSI Chemical Thermodynamic Database 2020, resulting in the updated sorption thermodynamic database for illite and montmorillonite (STDB 2023). STDB 2023 includes surface complexation constants and cation exchange selectivity coefficients for a wide range of radionuclides and transition metals, including Cs(I), Cd(II), Co(II), Fe(II), Mn(II), Ni(II), Pb(II), Ra(II), Zn(II), Am(III), Cm(III) (for illite only), Eu(III), Pu(III, IV), Np(IV, V), Sn(IV), Th(IV), Nb(V), Pa(V), and U(IV, VI). These parameters were derived by refitting extensive experimental sorption data collected over decades at the Paul Scherrer Institute and supplemented with data from the open literature. Surface complexation constants and cation exchange selectivity coefficients are now provided with uncertainty estimates, determined using the 95% confidence intervals obtained via the Monte Carlo sampling method implemented in the GEMSFITS parameter optimization tool. To address gaps caused by missing experimental data, the surface complexation constants for Cd(II) and Fe(II) on illite, as well as those for Pu(III, IV), Np(IV), and U(IV) on both illite and montmorillonite, were estimated based on linear correlations between the stability constants of aqueous and surface complexes. For geochemical calculations, STDB 2023 is supplemented with selectivity coefficients for naturally occurring major cations, sourced from recent literature. Versions of the model and database adapted for PHREEQC and The Geochemist’s Workbench are also provided.