Gradient estimates for the Lagrangian mean curvature equation with critical and supercritical phase