As mentioned in Section 1.2.1, the LJ parameters were optimized first, and the resulting force field from step 1 was used as input for valence parameter training. With the optimization procedure described in Section 2.3, the ForceBalance run satisfied the convergence criteria after 14 steps, and the drop in objective function value is shown in SI Figure .