Initial GAFF parameters were obtained using Ambertools, and the initial structures in different compositions were constructed using Packmol. All non-differentiable simulations were conducted using OpenMM. The entire optimization process was conducted on a single V100 card, undergoing over 60 iterations and taking approximately 10 days to complete. All procedures were fully automated, requiring no manual intervention.