To calculate the second virial coefficient (B2), we first conducted PIMD simulation (with 32 beads) to sample the monomer geometries to account for NQE. The intramolecular potential used in the PIMD sampling was identical to the MB-pol, and the simulation was run using the MB-pol plugin in OpenMM-7.4. Then the Mayer function