Neutron electric dipole moment from isovector quark chromo-electric dipole moment

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.