  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.
