We improve the estimate of the axion CDM energy density by considering the new values of current quark masses, the QCD phase transition effect and a possible anharmonic effect.