We show that, by constructing a magnetic-field density functional theory for the response of many-electron systems to the magnetic field generated by two magnetic dipoles, one can calculate the orbital portion of the indirect spin-spin coupling, nuclear magnetic shielding tenser, and diamagnetic susceptibility of any such system.