  We present a novel algorithm which can overcome the drawbacks of the
conventional linear scaling method with minimal computational overhead. This is
achieved by introducing additional constraints, thus eliminating the redundancy
of the orbitals. The performance of our algorithm is evaluated in ab initio
molecular-dynamics simulations as well as in a model system.
