This paper is concerned with the new error analysis of a Hodge-decomposition based
finite element method for the time-dependent Ginzburg-Landau equations in superconductivity. In this approach, the original equation of magnetic potential $\boldsymbol{A}$ is replaced by a new
system consisting of four scalar variables. As a result, the conventional Lagrange finite
element method (FEM) can be applied to problems defined on non-smooth domains. It is
known that due to the low regularity of $\boldsymbol{A},$ conventional FEM, if applied to the original
Ginzburg-Landau system directly, may converge to the unphysical solution. The main
purpose of this paper is to establish an optimal error estimate for the order parameter in
spatial direction, as previous analysis only gave a sub-optimal convergence rate analysis
for all three variables due to coupling of variables. The analysis is based on a nonstandard
quasi-projection for ψ and the corresponding negative-norm estimate for the classical Ritz
projection. Our numerical experiments confirm the optimal convergence of $ψ_h.$