We use the Yang-Mills gradient flow to calculate the pseudo-scalar expansion coefficient $c_P^*(t_f)$. This quantity is a key ingredient to obtaining the chiral condensate and strange quark content of the nucleon using the Lattice QCD formulation, which can ultimately determine the spin independent (SI) elastic cross section of dark matter models involving WIMP-nucleon interactions. The goal, using the gradient flow, is to renormalize the chiral condensate and the strange content of the nucleon without a power divergent subtraction. Using Chiral symmetry and the small flow time expansion of the gradient flow, the scalar density at zero flow time can be related to the pseudo-scalar density at non zero flow time. By computing the flow time dependance of the pseudo-scalar density over multiple lattices box sizes, lattice spacings and pion masses, we can obtain the scalar density of the nucleon. Our lattice ensembles are $N_{f}=2+1$, PCAC-CS gauge field configurations, varying over $m_{\pi}\approx \{410,570,700\}$~MeV at $a=0.0907$~fm, with additional ensembles that vary $a\approx \{0.1095,0.0936,0.0684\} $~fm at $m_{\pi} \approx 700$~MeV.