Energy minimization methods Shell model

In most electronegativity equalization models, if the energy is quadratic in the charges (as in Eq. [36]), the minimization condition (Eq. [41]) leads to a coupled set of linear equations for the charges. As with the polarizable point dipole and shell models, solving for the charges can be done by matrix inversion, iteration, or extended Lagrangian methods. 

---