We start by briefly summarizing the Hamiltonian hybrid particle-field molecular dynamics approach, for which parameters will be optimized. Consider a system of M molecules, with the i th molecule containing N i particles at positions R i = {r j } N i j=1 with conjugate momenta P i = {p j } N i j=1 , subject to the Hamiltonian 14