Gauge invariance, infrared/collinear singularitiesand tree level matrix element for $e^ + e^- \to \nu_e \bar{\nu}_e \gamma \gamma$