This is the first of a series of papers on a systematic application of the
NRQED bound state theory of Caswell and Lepage to higher-order radiative
corrections to the hyperfine structure of the muonium ground state. This paper
describes the calculation of the α2(Zα) radiative correction. Our
result for the complete α2(Zα) correction is 0.424(4) kHz, which
reduces the theoretical uncertainty significantly. The remaining uncertainty is
dominated by that of the numerical evaluation of the nonlogarithmic part of the
α(Zα)2 term and logarithmic terms of order α4.Comment: 56 pages, Rev.tex V3.0 and epsf.tex. 12 postscript files are called
in the text. Version accepted by Phys. Rev. D. A new table is adde