Precision calculation of γd↦πnn within chiral perturbation theory