# 66 neutrons
# rho        AV8'+UIX         AV8'
0.04       6.787  0.011      6.551  0.006
0.05       7.732  0.011      7.367  0.010
0.06       8.652  0.009      8.113  0.007
0.07       9.567  0.012      8.798  0.008
0.08      10.491  0.010      9.467  0.010
0.09      11.402  0.014     10.118  0.008
0.10      12.390  0.011     10.745  0.008
0.11      13.388  0.009     11.370  0.010
0.12      14.416  0.010     12.001  0.010
0.13      15.519  0.010     12.643  0.012
0.14      16.661  0.014     13.213  0.013
0.15      17.865  0.016     13.844  0.019
0.16      19.103  0.019     14.465  0.023
0.17      20.409  0.033     15.135  0.014
0.18      21.850  0.018     15.807  0.011
0.19      23.282  0.019     16.455  0.014
0.20      24.834  0.026     17.116  0.016
0.21      26.460  0.024
0.22      28.205  0.028     18.511  0.017
0.23      29.917  0.024
0.24      31.845  0.027     19.981  0.027
0.25      33.743  0.031
0.26      35.742  0.035     21.446  0.021
0.27      37.895  0.042
0.28      40.091  0.035     22.998  0.030
0.29      42.389  0.043
0.30      44.859  0.045     24.709  0.026
0.31      47.336  0.044
0.32      49.863  0.052     26.450  0.033
0.36      61.242  0.058                  
0.40      74.191  0.051     34.061  0.050
0.44      89.156  0.088                  
0.48     105.948  0.103     42.991  0.079
0.50     115.170  0.096     45.425  0.071













# 66 neutrons, AV8'+ibox + modified TNI

# E(0.16) = 17.65 +- 0.1 MeV

# FM(UIX) + VR (modified UIX)
#alpha*UIX     rho       E
0.60           0.04    6.708 0.008
0.60           0.08   10.219 0.011
0.60           0.12   13.702 0.012
0.60           0.16   17.644 0.013
0.60           0.20   22.384 0.023
0.60           0.24   27.904 0.024
0.60           0.32   41.952 0.044
0.60           0.40   60.364 0.059
0.60           0.48   83.695 0.061


# FM(UIX) + productory, v_ij*v_jk*v_ik, v(r)=exp(-2*mu*r)
#V0    mu      rho       E
550    150     0.04    6.734 0.010
550    150     0.08   10.188 0.009
550    150     0.12   13.705 0.015
550    150     0.16   17.723 0.015
550    150     0.20   22.381 0.020
550    150     0.24   27.847 0.024
550    150     0.32   41.406 0.043
550    150     0.40   58.635 0.062
550    150     0.48   79.740 0.075


# FM(UIX) + sum, v_ij*v_jk+perm., v(r)=exp(-2*mu*r), iboxv3
#V0    mu      rho       E
5.5    150     0.04    6.774 0.009
5.5    150     0.08   10.258 0.011
5.5    150     0.12   13.774 0.013
5.5    150     0.16   17.732 0.020
5.5    150     0.20   22.257 0.024
5.5    150     0.24   27.471 0.034
5.5    150     0.32   40.011 0.039
5.5    150     0.40   55.476 0.062
5.5    150     0.48   74.148 0.086 


# FM(UIX) + sum, v_ij*v_jk+perm, v(r)=exp(-2*mu*r)
#V0    mu      rho       E
1000   300     0.04    6.734 0.008
1000   300     0.06    8.516 0.008
1000   300     0.08   10.209 0.008
1000   300     0.10   11.940 0.010
1000   300     0.12   13.752 0.012
1000   300     0.14   15.681 0.014
1000   300     0.16   17.756 0.015
1000   300     0.20   22.442 0.021
1000   300     0.24   27.960 0.022
1000   300     0.32   41.689 0.043
1000   300     0.40   59.304 0.064
1000   300     0.48   81.123 0.067


# FM(UIX) + sum, v_ij*v_jk+perm, v(r)=exp(-2*mu*r)
#V0    mu      rho       E
700000 600     0.04    6.692 0.007
700000 600     0.08   10.165 0.009
700000 600     0.12   13.584 0.015
700000 600     0.16   17.564 0.016
700000 600     0.20   22.197 0.020
700000 600     0.24   27.662 0.046
700000 600     0.32   41.951 0.049
700000 600     0.40   61.248 0.085
700000 600     0.48   86.621 0.105


# sum, v_ij*v_jk+perm, v(r)=exp(-2*mu*r)
1340   300     0.08   10.165 0.017
1340   300     0.12   13.653 0.016
1340   300     0.16   17.728 0.022
1340   300     0.20   22.600 0.030
1340   300     0.24   28.501 0.027


# V3pi + V_R (as Illinois models)
#V3pi    Vr      rho       E
0.001   0.00835  0.08     10.087 0.009
0.001   0.00835  0.12     13.559 0.013
0.001   0.00835  0.16     17.698 0.016
0.001   0.00835  0.20     22.672 0.018
0.001   0.00835  0.24     28.692 0.031
0.001   0.00835  0.32     44.444 0.040
0.001   0.00835  0.40     65.972 0.074
0.001   0.00835  0.48     94.379 0.091
#V3pi    Vr      rho       E
0.0015  0.01085  0.08     10.066 0.010
0.0015  0.01085  0.12     13.558 0.015
0.0015  0.01085  0.16     17.688 0.016
0.0015  0.01085  0.20     22.686 0.020
0.0015  0.01085  0.24     28.739 0.029
0.0015  0.01085  0.32     44.601 0.054
0.0015  0.01085  0.40     66.360 0.079
0.0015  0.01085  0.48     94.755 0.098
#V3pi    Vr      rho       E
0.0022  0.0147   0.08     10.086 0.009
0.0022  0.0147   0.12     13.601 0.016
0.0022  0.0147   0.16     17.735 0.018
0.0022  0.0147   0.20     22.794 0.025
0.0022  0.0147   0.24     28.923 0.040
0.0022  0.0147   0.32     44.780 0.060
0.0022  0.0147   0.40     66.574 0.138
0.0022  0.0147   0.48     95.275 0.173



# E(0.16) = 16.05 +- 0.1 MeV

# FM(UIX) + VR (modified UIX)
#alpha*UIX     rho       E
0.19           0.08    9.872 0.008
0.19           0.12   12.864 0.011
0.19           0.16   16.102 0.016
0.19           0.20   19.640 0.019
0.19           0.24   23.676 0.021
0.19           0.32   33.279 0.034
0.19           0.40   45.280 0.041
0.19           0.48   59.746 0.066


# FM(UIX) + sum, v_ij*v_jk+perm, v(r)=exp(-2*mu*r)
#V0    mu      rho       E
1.7    150     0.08    9.895 0.012
1.7    150     0.12   12.904 0.012
1.7    150     0.16   16.081 0.013
1.7    150     0.20   19.627 0.019
1.7    150     0.24   23.528 0.026
1.7    150     0.32   32.677 0.036
1.7    150     0.40   43.702 0.048
1.7    150     0.48   56.603 0.079


# FM(UIX) + sum, v_ij*v_jk+perm, v(r)=exp(-2*mu*r)
#V0    mu      rho       E
280    300     0.04    6.652 0.012
280    300     0.08    9.865 0.014
280    300     0.12   12.818 0.014
280    300     0.16   16.010 0.016
280    300     0.20   19.495 0.021
280    300     0.24   23.455 0.026
280    300     0.32   32.691 0.039
280    300     0.40   44.131 0.053
280    300     0.48   57.736 0.077


# FM(UIX) + sum, v_ij*v_jk+perm, v(r)=exp(-2*mu*r)
#V0    mu      rho       E
200000 600     0.04    6.646 0.010
200000 600     0.08    9.836 0.010
200000 600     0.12   12.826 0.021
200000 600     0.16   16.048 0.021
200000 600     0.20   19.605 0.018
200000 600     0.24   23.666 0.021
200000 600     0.32   33.391 0.041
200000 600     0.40   45.726 0.066
200000 600     0.48   61.129 0.073


# V3pi + V_R (as Illinois models)
#V3pi    Vr      rho       E
0.0022  0.0125   0.08      9.712 0.013
0.0022  0.0125   0.12     12.713 0.012
0.0022  0.0125   0.16     16.012 0.016
0.0022  0.0125   0.20     19.903 0.020
0.0022  0.0125   0.24     24.436 0.032
0.0022  0.0125   0.32     35.713 0.067
0.0022  0.0125   0.40     50.902 0.062
0.0022  0.0125   0.48     70.506 0.098






# fit to several EOS
# E(rho) = a*rho**b + c*rho**d   
# E(rho) is the energy per neutron in MeV, and rho is the number density in fm**-3

#Hamiltonian                  E_sym  L                 a          b            c          d
  FM(UIX)+VR(modified UIX)      33.7   52.9            34.1 (7)   0.512 (7)    350 (3)    2.39 (2)
5 FM(UIX)+sum,mu=150            33.7   51.5            30.0 (6)   0.475 (6)    251 (2)    2.12 (1)
  FM(UIX)+sum,mu=300            33.8   53.0            32.2 (2)   0.495 (4)    315 (2)    2.29 (1)
  FM(UIX)+sum,mu=600            33.6   52.3            36.3 (5)   0.531 (4)    399 (4)    2.54 (2). 
  V3pi+VR(no FM),V3pi=0.001     33.7   55.6            32.9 (9)   0.50 (1)     450 (4)    2.51 (2)
  V3pi+VR(no FM),V3pi=0.0015    33.7   55.8            32.5 (3)   0.498 (4)    452 (1)    2.499 (7)
6 V3pi+VR(no FM),V3pi=0.0022    33.8   56.2            32.4 (9)   0.50 (1)     452 (5)    2.49 (2)

  FM(UIX)+VR(modified UIX)      32.1   41.3            30.4 (6)   0.476 (7)    195 (1)    2.22 (2)
  FM(UIX)+sum,mu=150            32.1   40.8            27.4 (7)   0.443 (9)    164.2 (9)  2.04 (2)
2 FM(UIX)+sum,mu=150 (fit2)     32.1   40.8            30.5 (1)   0.48  (1)    167 (1)    2.12 (1)
3 FM(UIX)+sum,mu=300            32.0   40.6            31.3 (7)   0.488 (7)    180 (3)    2.20 (2)
  FM(UIX)+sum,mu=600            32.0   41.3            33.0 (3)   0.503 (3)    213 (1)    2.34 (1)
4 V3pi+VR(no FM),V3pi=0.0022    32.0   44.0            32 (1)     0.49 (1)     297 (4)    2.47 (3)

  FM(UIX)                       31.3   35.7            24.3 (6)   0.404 (8)    120.8 (5)  1.90 (2)

1 AV8'                          30.5   31.3            31.0 (3)   0.486 (3)    112 (1)    2.26 (2)
7 AV8'+UIX                      35.1   63.6            34.5 (2)   0.514 (2)    492 (1)    2.436 (4)
