Implicit solvent models aim to represent the effect a solvent exerts on a solute by an electrostatic continuum. This approximation not only reduces the computational burden due to a much smaller number of explicit particles in the system but also significantly accelerates sampling efficiency due to the instantaneous averaging of solvent configurations and reduced viscosity . Many implicit solvent models, especially those that are used in conjunction with QM methods, consider the solute's electron density when solving the Poisson-Boltzmann equation . Examples include the conductor-like screening (COSMO) , the conductor-like polarizable continuum (CPCM) , and the solvation model based on density (SMD) models. A more complex solvent model that aims to include explicit-solvent effects was coined COSMO-RS . Starting from perfectly screened molecules in the conductor-like approach, the COSMO-RS model uses methods from statistical thermodynamics to pair charged surface segments. Due to the relatively high computational burden associated with solving the Poisson-Boltzmann equation, force-field methods are typically paired with implicit solvent models that are based on the generalized Born (GB) model, which does not consider the solute's electron density but rather atomic monopoles. Examples of these models are GB-HCT , GB-OBC , GB-Neck , and GB-Neck2 .