The intramolecular terms were initialized using GAFF 3 with AM1-BCC charge. The optimization targets are the experimental density and the center-of-mass (COM) RDF from neutron diffraction at 30 • C (303 K), respectively. During the optimization, we keep monitoring the effective sample size of the estimator, and resampling will be run if the sample size is smaller than a predetermined threshold.