
======================================================================================
Adams, Mol. Phys. 1975, 29, 1, 307-311
======================================================================================

T = T* * 0.120032 * 1000
4 = 480
2 = 240
B = muex/kT + ln(<N>)
c " 1.063*998/NA/(1.38065e-23*480) + log(114)" = 5.00202
c " 2.828*998/NA/(1.38065e-23*480) + log(199)" = 6.00049
c " 4.116*998/NA/(1.38065e-23*480) + log(238)" = 6.50154
c "-0.736*998/NA/(1.38065e-23*240) + log(130)" = 4.49944
muex = kT(B-ln<N>)

--------------------------------------------------------------------------------------
     B     muex  T     rho     -U     p    NMC    <N>   V
a1  5.0   1.063 4.0   0.2276 1.205 1.081  500000 114  500
a2  6.0   2.828 4.0   0.3979 2.082 2.397  500000 199  500
a3  6.5   4.116 4.0   0.4754 2.462 3.285  500000 238  500
a4  4.5  -0.739 2.0   0.1301 0.893 0.238  500000 130 1000

mos

a1  5.0   1.044 4.0   0.2282 1.223        500000 114  500.8  * rand instaed of drand48
a1  5.0   0.980 4.0   0.2319 1.251        500000 116  500.8
a2  6.0   2.812 4.0   0.3988 2.112        500000 200  500.8
a3  6.5   4.081 4.0   0.4788 2.508        500000 240  500.8
a4  4.5  -0.565 2.0   0.1195 0.790        500000 119  999.2
--------------------------------------------------------------------------------------


======================================================================================
Mezei, Mol. Phys. 1987, 61, 3, 565-582
======================================================================================

for the epsilon/sigma used here the temperatures
of Mezei correspond to: T = 90.024, 90.024, 108.389 K
B values differ because same mu different V and hence <N>.
i.e.: 
c "kB*90.024*(-1.398-log(99.6))*0.001*NA" -4.49038
c "kB*90.024*(-1.23-log(101.3))*0.001*NA" -4.3773
c "kB*108.389*(2.47-log(178.0))*0.001*NA" -2.44385


--------------------------------------------------------------------------------------
     B     muex     T      rho     -U    p     NMC   <N>   V
m1 -1.398  -4.49038 0.75   0.8    5.77  0.314  3e+06  99.6 125.0
m2 -1.23   -4.3773  0.75   0.81   5.85  0.194  3e+06 101.3 125.0
m3  2.47   -2.44385 0.903  0.835  5.85  1.020  3e+06 178.0 213.2

mos

m1 -1.396  -4.49729 0.75   0.806  5.82  0.602  3e+06 100.8 125.0*
m1 -1.396  -4.50291 0.75   0.802  5.81  0.593$ 3e+06 100.3 125.0
m2 -1.23   -4.37233 0.75   0.805  5.81         3e+06 100.6 125.0
m3  2.47   -2.45445 0.903  0.845  5.93         3e+06 180.1 213.2
--------------------------------------------------------------------------------------

* from GCMD
$ P from NVT MD

