Thus, the derivative of the wavefunction parameters are not explicitly evaluated because the multipliers have been determined so that the Lagrangian satisfies Eq. ( ). The last three terms in Eq. ( ) account for the displacement of tesserae. In this study, the position of tesserae is allowed to move, and these terms are evaluated in each geometry optimization step. With C-PCM, no approximations are required for the evaluation.