Computation of forces arising from the polarizable continuum model within the domain-decomposition paradigm