Up-to N³LO heavy-baryon chiral perturbation theory calculation for the M1 properties of three-nucleon systems

M1 properties, comprising magnetic moments and radiative capture of thermal neutron observables, are studied in two- and three-nucleon systems. We utilize meson exchange current derived up to N$^3$LO using heavy baryon chiral perturbation theory a la Weinberg. Calculations have been performed for several qualitatively different realistic nuclear Hamiltonians, which permits us to analyze model dependence of our results. Our results are found to be strongly correlated with the effective range parameters such as binding energies and the scattering lengths. Taking into account such correlation, the results are in good agreement with the experimental data with small model-dependence.