M1 properties, comprising magnetic moments and radiative capture of thermal neutron observables, are studied in two-and three-nucleon systems. We use meson exchange current derived up to N 3 LO using heavy baryon chiral perturbation theoryà 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 correlations, the results are in good agreement with the experimental data with small model dependence. YOUNG-HO SONG, RIMANTAS LAZAUSKAS, AND TAE-SUN PARK PHYSICAL REVIEW C 79, 064002 (2009) F JM (r) = j J (qr)Y JM (r), (20) F M JL (r) = j L (qr)Y M J L1 (r), where j L (qr) is the spherical Bessel function; Y JM and Y M J L1are spherical and vector-spherical harmonics, respectively; r is a vector describing the particle (nucleon or meson), which interacts with EM field. Then the electric and magnetic multipoles read064002-4 UP TO N 3 LO HEAVY-BARYON CHIRAL . . . PHYSICAL REVIEW C 79, 064002 (2009) E JM = i J + 1 2J + 1 1/2 F M J J −1 (r) − J 2J + 1 1/2 F M J J +1 (r) · j c + qF M J J (r) · j µ .