Better model-independent theoretical predictions and input to nucleon structure functions, including parton distribution functions are vital to improving our understanding of how the nucleon is formed by quarks and gluons. Better determination here is important for experiments probing Hadron structure, including at the LHC and future experiments at JLab, or the EIC. We show how to examine the gluon form factors using gradient flow at non-zero Q^2. By using a second order extension of the Feynman-Hellmann theorem we also examine the nucleon Compton structure functions T_1(omega,Q^2) without the need for calculations of four-point correlators. The calculation is performed for multiple lattice spacings, volumes and a large range of Q^2 up to 15 GeV^2, allowing examination of lattice artefacts and higher twist contributions.