Nuclear axial current operators to fourth order in chiral effective field theory

We present the complete derivation of the nuclear axial charge and current operators as well as the pseudoscalar operators to fourth order in the chiral expansion relative to the dominant one-body contribution using the method of unitary transformation. We demonstrate that the unitary ambiguity in the resulting operators can be eliminated by the requirement of renormalizability and by matching of the pion-pole contributions to the nuclear forces. We give expressions for the renormalized single-, two- and three-nucleon contributions to the charge and current operators and pseudoscalar operators including the relevant relativistic corrections. We also verify explicitly the validity of the continuity equation.