Force field-based models are a Newtonian mechanics approximation of reality and are inherently noisy. Coupling models from different molecular scale domains (including single, gas-phase molecules up to multimolecule, condensed phase ensembles) is difficult, which is also the case for finding solutions that transfer well between the scales. In this contribution, we introduce a surrogate-assisted algorithm to optimize Lennard-Jones parameters for target data from different scale domains to overcome the difficulties named above. Specifically, our approach combines a surrogate-assisted global evolutionary optimization method with a presampling phase that takes advantage of one scale domain being less computationally expensive to evaluate. The algorithm's components were evaluated individually, elucidating their individual merits. Our findings show that the process of parametrizing force fields can significantly benefit from both the presampling method, which alleviates the need to have a good initial guess for the parameters, and the surrogate model, which improves efficiency.