  We present results from our lattice QCD study of the contribution of the
isovector quark cEDM (qcEDM) operator to the neutron EDM. The calculation was
carried out on four 2+1+1-flavor highly improved staggered quark ensembles
(provided to us by the MILC collaboration) using Wilson-clover quarks to
construct correlation functions. We use the nonsinglet axial Ward identity
including corrections up to O(a) to show how to control the power-divergent
mixing of the isovector qcEDM operator with the lower dimensional pseudoscalar
operator. Results for the nEDM are presented after conversion to the MS scheme
at the leading-log order.
