This can be solved by using the (quasi-)Newton method with respect to R. For the full optimization problem consisting of Eqs. ( )-( ), we have two coupled optimization problems with respect to participation coefficients a and nuclear coordinates R. The former seeks a that maximizes the functionality p for a given fixed R. The latter searches for the equilibrium geometry with local energy minimization for a given a.