  A status report on the calculations of the contribution of four CP violating
operators, the $\Theta$-term, the quark EDM, the chromo EDM and the Weinberg
operator to the neutron EDM are presented. At this time, there exit precise
physical results only for the quark EDM operator by the PNDME collaboration.
First results showing signal in the contributions of the $\Theta$-term and the
connected part of the chromo EDM operator have been presented. The challenge of
divergent mixing in the chromo EDM and Weinberg operators has motivated
calculations in the gradient flow scheme. While there has been steady progress,
the challenges remaining are large. Results with $O(50\%)$ uncertainty with
control over all systematic errors can be expected for the $\Theta$-term over
the next five years. Prediction of a timeline for progress on the chromo EDM
and the Weinberg operators will depend on when the renormalization and
divergent mixing of these operators is brought under control. The most
optimistic scenario is that the gradient flow scheme provides a solution to the
numerical signal and mixing problems for both the gluonic and quark operators.
