Optimal convergence rates in L2 for a first order system least squares finite element method

We analyze a divergence based first order system least squares method applied to a second order elliptic model problem with homogeneous boundary conditions. We prove optimal convergence in the L2(Ω) norm for the scalar variable. Numerical results confirm our findings.