Constrained density functional theory calculation with iterative optimization

An iterative optimization approach that simultaneously minimizes the energy and optimizes the Lagrange multipliers enforcing desired constraints is presented. The method is tested on previously established benchmark systems and it is proved to be efficient and accurate. The approach can also be efficiently used when the constraint is not a scalar quantity but a spatially varying function such as the charge density distribution.