(1), Zn (2) and Cd (3)] have been synthesized and characterized by powder and single-crystal X-ray diffraction, vibrational spectroscopy (IR and Raman) and simultaneous thermal analysis (TG-DTA-MS).Crystal structure analyses of compounds 1-3 have revealed that both phosphonate groups of (pyridin-1-ium-3,5-diyl)diphosphonate ligands serve to extend the metal cations into a strongly undulated grid structure parallel to the (101) plane. The adjacent monolayers are held together via multiple hydrogen bonds and offset face-to-face p-p interactions. The interlayer interaction energy values, estimated using the CRYSTAL code, have revealed that the stabilization energies of compounds 1-3 are quite similar, where 1 is least and 3 most advantageous. Computational studies have been further related to the results of thermal analyses. In the case of compound 1, which has unpaired electrons at the metal centre, the investigation has been supplemented by magnetic measurements. The magnetic properties of 1 (1.7-300 K) have been analyzed assuming that any possible exchange interactions are very weak.The obtained best-fit parameters are: D (axial splitting of 4 T 1g term) ¼ À478 cm À1 , l (spin-orbit coupling) ¼ À103 cm À1 , a (orbital reduction factor) ¼ 1.38, and zJ 0 (mean-field exchange parameter) ¼ À0.11 cm À1 . The negative value of D implies that the ground state is the orbital doublet, which cannot be treated within the spin-Hamiltonian approach.