A comparison between existing nuclear magnetic resonance measurements and calculations of the scalar spin-spin interaction (J-coupling) in deuterated molecular hydrogen (HD) yields stringent constraints on anomalous spin-dependent potentials between nucleons at the atomic scale (∼ 1Å). The dimensionless coupling constant g p P g N P /4π associated with exchange of pseudoscalar (axion-like) bosons between nucleons is constrained to be less than 5 × 10 −7 for boson masses in the range of 5 keV, representing improvement by a factor of 100 over previous constraints. The dimensionless coupling constant g −19 for bosons of mass < ∼ 1000 eV, improving constraints at this distance scale by a factor of 100 for proton-proton couplings and more than 8 orders of magnitude for neutron-proton couplings.PACS numbers: 14.80. Va, 11.40.Ha, Over the past few decades, searches for anomalous spin-dependent forces have drawn considerable interest as a signature of the axion [1], a hypothetical pseudoscalar Goldstone boson [2,3] arising out of the Peccei-Quinn solution to the strong-CP problem [4,5]. Axions are also appealing as a candidate for dark matter [6]. Other exotic spin-dependent interactions are predicted by a variety of novel theories such as those involving para-photons [7] and unparticles [8].The possible spin-dependent forces that could arise from exchange of scalar/pseudoscalar or vector gauge bosons are enumerated in Ref. [9]. Another theoretical framework for considering anomalous spin-dependent interactions comes from introducing nonzero torsion into general relativity, which causes gravity to acquire scalar and vector components that manifest as new spin-mass and spin-spin couplings [10][11][12][13]. Spontaneous Lorentz violation may also generate exotic spindependent interactions [14].Recent experiments have significantly improved constraints on long-range (tens of cm) dipole-dipole interactions between neutrons [15, 16] and monopole-dipole interactions between electrons and nucleons [17]. Earlier work constrained monopole-dipole interactions between nuclei [18] and anomalous dipole-dipole interactions between electrons [19,20] and electrons and nuclei [21]. Constraints on new pseudoscalar interactions can also be obtained from spin-independent tests of the inverse square law [22]. Other laboratory experiments constrain anomalous monopole-dipole couplings on length scales in the range of µm to mm [23][24][25][26]. Spin relaxation studies of hyperpolarized 3 He gas have been used to limit anomalous dipole-dipole interactions between neutrons at length scales of about 100 nm [27]. Because of the possibility that new force-mediating bosons may be massive and have limited range, it is of considerable interest to find experimental techniques to search for anomalous spin-dependent interactions at even shorter distances.Here we discuss constraints on the existence of anomalous dipole-dipole forces on angstrom length scales. Our constraints are obtained by comparison of nuclear magnetic resonance (NMR) measurements and th...