In principle, the impact of both charge transfer and solvation effects to asymptotic parameters can be accounted for by extending the training dataset from single molecule fragments to many-body clusters. However, the increased size of training cluster renders the longrange parameter calculations excessively expensive. The chemical and conformation space of many-body clusters is also considerably larger than that of single fragments, making it much more challenging to construct the database. Considering that the remaining error showed in this work is reasonably low, we are inclined to keep the simplicity of the training dataset, which we expect to be accurate enough for nonreactive MD simulations.