Abstract. In this talk, we motivate the calculation of the matrix elements of novel CP violating operators, the quark EDM and the quark chromo EDM operators, within the nucleon state using lattice QCD. These matrix elements, combined with the bound on the neutron EDM, would provide stringent constraints on beyond the standard model physics, especially as the next generation of neutron EDM experiments reduce the current bound. We then present our lattice strategy for the calculation of these matrix elements, in particular we describe the use of the Schrodinger source method to reduce the calculation of the 4-point to 3-point functions needed to evaluate the quark chromo EDM contribution. We end with a status report on the quality of the signal obtained in the lattice calculations of the connected contributions to the quark chromo EDM operator and the pseudoscalar operator it mixes with under renormalization.