  We derive the axion-nucleon interaction Lagrangian in heavy baryon chiral
perturbation theory up to next-to-next-to-leading order. The effective
axion-nucleon coupling is calculated to a few percent accuracy.
