We present results for the neutron electric dipole moment due to the dimension-4 and 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.