“…where G andG are the gluonic field tensor and its dual, b = 1, 2, ..., 8 is the color index, g 2 /4π is the color coupling constant, N andN = N † γ 0 are the nucleon field and its Dirac adjoint, f a is the axion decay constant, and C G and C N are model-dependent dimensionless parameters. Astrophysical constraints on the axion-gluon coupling in (1) come from Big Bang nucleosynthesis [36][37][38]: m 1/4 a f a /C G 10 10 GeV 5/4 for m a 10 −16 eV and m a f a /C G 10 −9 GeV 2 for m a 10 −16 eV, assuming that axions saturate the presentday DM energy density, and from supernova energyloss bounds [35,39]: f a /C G 10 6 GeV for m a 3 × 10 7 eV. Astrophysical constraints on the axion-nucleon coupling in (1) come from supernova energy-loss bounds [39,40]: f a /C N 10 9 GeV for m a 3 × 10 7 eV, while existing laboratory constraints come from magnetometry searches for new spin-dependent forces mediated by axion exchange [41]: f a /C N 1 × 10 4 GeV for m a 10 −7 eV.…”