Minimization methods for the one-particle dirac equation

Taking into account relativistic effects in quantum chemistry is crucial for accurate computations involving heavy atoms. Standard numerical methods can deal with the problem of variational collapse and the appearance of spurious roots only in special cases. The goal of this Letter is to provide a general and robust method to compute particle bound states of the Dirac equation.