Fig6a.dat: Data for Figure. 6a; temperature in kelvin (column 1) and molar susceptibility (column 2).



Fig6b.dat: Data for Figure. 6b; The file is made up of 4 pairs of columns corresponding to applied field in tesla and magnetisation in Bohr magnetons per formula unit for measurements at 0.63 K (Columns 1 and 2), 1.52 K (Columns 3 and 4), 4.0 K (columns 5 and 6) and 16.1 K (columns 7 and 8).



Fig6binsert.dat: Data for insert to Figure. 6b; applied field in tesla (column 1) and differential susceptibility in arbitrary units (column 2)



Fig6c.dat: Data for Figure. 6c; Temperature in kelvin (column 1) and total heat capacity divided by temperature (column 2).



Fig6d.dat: Data for Figure. 6d; The file is made up of 10 pairs of columns corresponding to temperature in kelvin and magnetic heat capacity for measurements at 0 T (columns 1 and 2), 1 T (columns 3 and 4), 2 T (columns 5 and 6), 3 T (columns 7 and 8), 4 T (columns 9 and 10), 5 T (columns 11 and 12), 6 T (columns 13 and 14), 7 T (columns 15 and 16), 8 T (columns 17 and 18) and 9 T (columns 19 and 20).



Fig7a.dat: Data for Figure. 7a; The file is made up of 7 pairs of columns corresponding to applied field in tesla and ESR transmission intensity in arbitrary units for measurements at 628.8 GHz (columns 1 and 2), 609.6 GHz (columns 3 and 4), 428.8 GHz (columns 5 and 6), 412.8 GHz (columns 7 and 8), 324 GHz (columns 9 and 10), 168.0 GHz (columns 11 and 12) and 156 GHz (columns 13 and 14).



Fig7b.dat: Data for Figure. 7b; The file is made up of 4 pairs of columns corresponding to applied field in tesla and ESR transmission intensity in arbitrary units for simulated measurement at 5 K and 412.8 GHz (columns 1 and 2) and 412.8 GHz data taken at 5 K (columns 3 and 4), 30 K (columns 5 and 6) and 80 K (columns 7 and 8)



Fig7c.dat: Data for Figure. 7c; The file is made up of 5 pairs of columns corresponding to field in tesla and frequency in GHz for observed resonances observed in ESR data and fitted to a Hamiltonian containing D and E terms and the zeeman term for applied field parallel to the x-axis (columns 1 and 2), y-axis (columns 3 and 4), z-axis (columns 5 and 6) and off-axis resonance (columns 7 and 8). Columns 9 and 10 contain field in tesla and frequency in GHz for off-axis resonances that weren't fitted to above Hamiltonian. Further columns show the resultant simulation from the fit to above Hamiltonian is contained in columns 11, 12, 13 (field in tesla) and 14 (frequency in GHz) for field parallel to the x-axis;  columns 15, 16, 17, 18 (field in tesla) and 19 (frequency in GHz) for field parallel to the y-axis; columns 20, 21, 22, 23 (field in tesla) and 24 (frequency in GHz) for field parallel to the z-axis; and columns 25 (field in tesla) and 26 (frequency in GHz) for fields parallel to the off-axis direction.



Fig8a.dat: Data for Figure. 8a; Temperature in kelvin (column 1) and molar susceptibility (column 2).



Fig8b.dat: Data for Figure. 8b; The file is made up of 4 pairs of columns corresponding to field in tesla and magnetisation in Bohr magnetons per formula unit for field sweeps at 0.6 K (columns 1 and 2), 2 K (columns 3 and 4), 4 K (columns 5 and 6) and 10 K (columns 7 and 8).



Fig8c.dat: Data for Figure. 8c; Field in tesla (column 1) and dM/dH data in arbitrary units (column 2).



Fig8d.dat: Data for Figure. 8d; Temperature in kelvin (column 1) and total heat capacity divided by temperature (column 2).



Fig8e.dat: Data for Figure. 8e; The file is made up of 11 pairs of columns corresponding to temperature in kelvin and magnetic heat capacity for measurements at 0 T (columns 1 and 2), 1 T (columns 3 and 4), 1.5 T (columns 5 and 6), 2 T (columns 7 and 8), 3 T (columns 9 and 10), 4 T (columns 11 and 12), 5 T (columns 13 and 14), 6 T (columns 15 and 16), 7 T (columns 17 and 18) and 8 T (columns 19 and 20), 9 T (columns 21 and 22).



Fig10a.dat: Data for Figure. 10a; Temperature in kelvin (column 1) and molar susceptibility (column 2).



Fig10ainsert.dat: Data for insert to Figure. 10a; Temperature in kelvin (column 1) and inverse molar susceptibility (column 2).



Fig10b.dat: Data for Figure. 10b; Temperature in kelvin (column 1) and d(Chi*T)/dT data (column 2).


Fig10c.dat: Data for Figure. 10c; The file is made up of 4 pairs of columns corresponding to field in tesla and magnetisation in Bohr magnetons per formula unit for field sweeps at 0.6 K (columns 1 and 2), 1.6 K (columns 3 and 4), 5 K (columns 5 and 6) and 12 K (columns 7 and 8).


Fig10d.dat: Data for Figure. 10d; The file is made up of 4 pairs of columns corresponding to field in tesla and dM/dH in arbitrary units for field sweeps at 0.6 K (columns 1 and 2), 1.6 K (columns 3 and 4), 5 K (columns 5 and 6) and 12 K (columns 7 and 8).


Fig10dinsert.dat: Data for insert to Figure. 10d; Temperature in kelvin (column 1) and d2MdH2 data (column 2).


Fig11a.dat: Data for Figure. 11a; Temperature in kelvin (column 1) and total heat capacity divided by temperature (column 2).



Fig11b.dat: Data for Figure. 11b; Temperature in kelvin (column 1), entropy (column 2), upper bound on entropy (column 3) and lower bound on entropy (column 4).



Fig11c.dat: Data for Figure. 11c; The file is made up of 10 pairs of columns corresponding to temperature in kelvin and magnetic heat capacity for measurements at 0 T (columns 1 and 2), 1 T (columns 3 and 4), 2 T (columns 5 and 6), 3 T (columns 7 and 8), 4 T (columns 9 and 10), 5 T (columns 11 and 12), 6 T (columns 13 and 14), 7 T (columns 15 and 16), 8 T (columns 17 and 18) and 9 T (columns 19 and 20).



Fig11cLeftInsert.dat: Data for the left insert of Figure. 11c; The file is made up of 6 pairs of columns corresponding to temperature in kelvin and magnetic heat capacity for measurements at 6.5 T (columns 1 and 2), 7 T (columns 3 and 4), 7.5 T (columns 5 and 6),

8 T (columns 7 and 8), 8.5 T (columns 9 and 10) and 9 T (columns 11 and 12)



Fig11cRightInsert.dat: Data for the right insert of Figure. 11c; Applied field normalised to the g-factor in units of Kelvin (column 1) and temperature in units of kelvin (column 2) and temperature error in units of kelvin (column 3)


Fig12a: Muon-spin relaxation spectra and fits are provided for data shown, measured at T=2 K (muon_2k.dat), and at 3.5 K for model 1 (muon_3_5k_v1.dat) and model 2 (muon_3_5k_v2.dat). The columns are explained in the file headers.

Fig12b_data.dat: Neutron powder diffraction data shown as red circles in the main pane of figure 12b. Left column gives d spacing in units Angstrom, and right column gives intensity, in arbitrary units.

Fig12b_fit.dat: Neutron powder diffraction fit shown as a black line in the main pane of figure 12b. Left column gives d spacing in units Angstrom, and right column gives intensity, in arbitrary units.

Fig12b_difference.dat: Neutron powder diffraction difference curve between data and fit shown as a blue line in the main pane of figure 12b. Left column gives d spacing in units Angstrom, and right column gives intensity, in arbitrary units.

Fig12b_left_inset_calculation.dat: Neutron powder diffraction calculation shown as a black line in the left inset of figure 12b. Left column gives d spacing in units Angstrom, and right column gives intensity, in arbitrary units.

Fig12b_right_inset_calculation.dat: Neutron powder diffraction calculation shown as a black line in the right inset of figure 12b. Left column gives d spacing in units Angstrom, and right column gives intensity, in arbitrary units.

Fig12c_Moment_temp_dep.dat: Data for moment size (column 2, arbitrary units) and associated error (column 3) vs temperature (column 1, kelvin).