Precision calculation of the axion-nucleon coupling in chiral perturbation theory