Spacetime Harmonic Functions and A Hawking Mass Monotonicity Formula for Initial Data Sets