The exponents method for calculating equilibrium concentrations of complex species in solution

The exponents method for calculating the concentrations of species in multimetal‐multiligand systems is introduced. This method uses the Newton‐Raphson method with restricted step iteration, which guarantees a monotonically decreasing objective function. Variable transformation and scaling are performed to avoid underflows and overflows during the calculations. A special linear solver using the eigenvalues and eigenvectors of the Jacobian matrix is implemented for overcoming disastrous singularity of this matrix, and the singular value decomposition method is applied for setting the initial guess. In addition, polynomial extrapolation is used for improving the performance when simulating a diagram of concentrations of species. The method was tested with 14 systems of different sizes over the whole pH range and presented robust and efficient behavior. © 1995 by John Wiley & Sons, Inc.