The most commonly used Lennard-Jones potential.
E = A/r12 - B/r6
A and B parameters are determined by fitting the models of liquids to experimental data (OPLS).
In the absence of such possibility they are using quantum mechanics, but do not expect the accuracy above 5-10%.