We study the deferred correction method for the magnetohydrodynamics (MHD) system written in Elsasser variables. The proposed algorithm is based on the penalty projection with grad-div stabilized Taylor Hood solutions of the Elsasser formulation. In this way, second order accuracy of the method in time is obtained through the deferred correction method with excellent mass conservation properties. For the proposed method, stability , is rigorously proven and numerical experiments are presented to verify the proposed scheme and theory. (C) 2019 Elsevier B.V. All rights reserved.