To realize multimodal optimization, we used the Differentiable Molecular Force Field (DMFF) program to perform trajectory reweighting, which avoids expensive backpropagation through MD trajectory. During the optimization, we targeted a large number of systems at different temperatures and formulations simultaneously. We constructed an objective function encompassing multiple properties to obtain parameters suitable for various systems. Specifically, the loss function is constructed as following: