  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.
