  We present results for the neutron electric dipole moment due to the
dimension 4 and dimension 6 gluonic CP violation, and the isovector quark
chromoelectric dipole moment using clover valence quarks on HISQ dynamical
ensembles generated by the MILC Collaboration. For the gluonic operators, we
use the gradient flow scheme to obtain divergence-free continuum results. For
the chromoelectric dipole moment operator, we use the unflowed local operator
but discuss how the quadratically divergent mixing with the pseudoscalar
operator can be controlled nonperturbatively.
