nucleardatapy 1.0.0__py3-none-any.whl → 1.0.1__py3-none-any.whl

This diff represents the content of publicly available package versions that have been released to one of the supported registries. The information contained in this diff is provided for informational purposes only and reflects changes between package versions as they appear in their respective public registries.
Files changed (87) hide show
  1. nucleardatapy/astro/setup_mr.py +115 -33
  2. nucleardatapy/create_folder.py +2 -2
  3. nucleardatapy/crust/setup_crust.py +5 -5
  4. nucleardatapy/data/astro/HESS/J1731-347.dat +4 -0
  5. nucleardatapy/data/astro/NICER/J0030+0451.dat +6 -6
  6. nucleardatapy/data/astro/NICER/J0437-4715.dat +4 -3
  7. nucleardatapy/data/astro/NICER/J0614-3329.dat +4 -0
  8. nucleardatapy/data/astro/NICER/J0740+6620.dat +5 -5
  9. nucleardatapy/data/hnuclei/1991-2L-Yamamoto.csv +6 -0
  10. nucleardatapy/data/hnuclei/2013-2L-Ahn.csv +1 -1
  11. nucleardatapy/data/hnuclei/2019-2L-Ekawa.csv +7 -0
  12. nucleardatapy/data/matter/micro/2020-SCGF-NM-N3LO-414-E2A.dat +21 -0
  13. nucleardatapy/data/matter/micro/2020-SCGF-NM-N3LO-414-TD.dat +22 -0
  14. nucleardatapy/data/matter/micro/2020-SCGF-NM-N3LO-450-E2A.dat +20 -0
  15. nucleardatapy/data/matter/micro/2020-SCGF-NM-N3LO-450-TD.dat +22 -0
  16. nucleardatapy/data/matter/micro/2020-SCGF-NM-N3LO-500-E2A.dat +23 -0
  17. nucleardatapy/data/matter/micro/2020-SCGF-NM-N3LO-500-TD.dat +22 -0
  18. nucleardatapy/data/matter/micro/2020-SCGF-SM-N3LO-414-E2A.dat +15 -0
  19. nucleardatapy/data/matter/micro/2020-SCGF-SM-N3LO-414-TD.dat +21 -0
  20. nucleardatapy/data/matter/micro/2020-SCGF-SM-N3LO-450-E2A.dat +15 -0
  21. nucleardatapy/data/matter/micro/2020-SCGF-SM-N3LO-450-TD.dat +21 -0
  22. nucleardatapy/data/matter/micro/2020-SCGF-SM-N3LO-500-E2A.dat +20 -0
  23. nucleardatapy/data/matter/micro/2020-SCGF-SM-N3LO-500-TD.dat +20 -0
  24. nucleardatapy/data/matter/micro/2024-ABI-NM-DeltaNNLO450.dat +28 -0
  25. nucleardatapy/data/matter/micro/2024-ABI-NM-DeltaNNLO500.dat +28 -0
  26. nucleardatapy/data/matter/micro/2024-ABI-NM-DeltaNNLOgo394.dat +28 -0
  27. nucleardatapy/data/matter/micro/2024-ABI-NM-DeltaNNLOgo450.dat +28 -0
  28. nucleardatapy/data/matter/micro/2024-ABI-NM-NNLOsat.dat +28 -0
  29. nucleardatapy/data/matter/micro/2024-ABI-SM-DeltaNNLO450.dat +28 -0
  30. nucleardatapy/data/matter/micro/2024-ABI-SM-DeltaNNLO500.dat +28 -0
  31. nucleardatapy/data/matter/micro/2024-ABI-SM-DeltaNNLOgo394.dat +28 -0
  32. nucleardatapy/data/matter/micro/2024-ABI-SM-DeltaNNLOgo450.dat +28 -0
  33. nucleardatapy/data/matter/micro/2024-ABI-SM-NNLOsat.dat +28 -0
  34. nucleardatapy/fig/astro_setupGW_fig.py +4 -2
  35. nucleardatapy/fig/astro_setupMR_fig.py +14 -8
  36. nucleardatapy/fig/astro_setupMasses_fig.py +4 -2
  37. nucleardatapy/fig/astro_setupMtov_fig.py +6 -2
  38. nucleardatapy/fig/astro_setupMup_fig.py +4 -2
  39. nucleardatapy/fig/corr_setupEsymDen_fig.py +10 -4
  40. nucleardatapy/fig/corr_setupEsymLsym_fig.py +7 -3
  41. nucleardatapy/fig/corr_setupKsatQsat_fig.py +4 -2
  42. nucleardatapy/fig/crust_setupCrust_fig.py +6 -4
  43. nucleardatapy/fig/eos_setupAMBeq_fig.py +162 -134
  44. nucleardatapy/fig/eos_setupAMLeq_fig.py +35 -32
  45. nucleardatapy/fig/eos_setupAM_asy_lep_fig.py +39 -30
  46. nucleardatapy/fig/eos_setupAM_asy_nuc_fig.py +39 -30
  47. nucleardatapy/fig/eos_setupAM_asy_tot_fig.py +39 -30
  48. nucleardatapy/fig/eos_setupAM_fig.py +45 -30
  49. nucleardatapy/fig/eos_setupCC_fig.py +54 -24
  50. nucleardatapy/fig/hnuc_setupChart_fig.py +17 -14
  51. nucleardatapy/fig/hnuc_setupRE1LExp_fig.py +12 -2
  52. nucleardatapy/fig/matter_all_fig.py +55 -38
  53. nucleardatapy/fig/matter_setupCheck_fig.py +6 -3
  54. nucleardatapy/fig/matter_setupFFGLep_fig.py +13 -9
  55. nucleardatapy/fig/matter_setupFFGNuc_fig.py +22 -14
  56. nucleardatapy/fig/matter_setupHIC_fig.py +9 -9
  57. nucleardatapy/fig/matter_setupMicroEsym_fig.py +31 -30
  58. nucleardatapy/fig/matter_setupMicro_LP_fig.py +11 -5
  59. nucleardatapy/fig/matter_setupMicro_band_fig.py +11 -6
  60. nucleardatapy/fig/matter_setupMicro_effmass_fig.py +5 -4
  61. nucleardatapy/fig/matter_setupMicro_err_NM_fig.py +4 -2
  62. nucleardatapy/fig/matter_setupMicro_fig.py +30 -22
  63. nucleardatapy/fig/matter_setupMicro_gap_fig.py +22 -12
  64. nucleardatapy/fig/matter_setupNEPStats_fig.py +13 -3
  65. nucleardatapy/fig/matter_setupPhenoEsym_fig.py +28 -27
  66. nucleardatapy/fig/matter_setupPheno_fig.py +22 -17
  67. nucleardatapy/fig/nuc_setupBEExp_chart_fig.py +40 -31
  68. nucleardatapy/fig/nuc_setupBEExp_fig.py +49 -46
  69. nucleardatapy/fig/nuc_setupBETheo_fig.py +55 -50
  70. nucleardatapy/fig/nuc_setupISGMRExp_fig.py +7 -10
  71. nucleardatapy/fig/nuc_setupRchExp_fig.py +12 -20
  72. nucleardatapy/fig/nuc_setupRchTheo_fig.py +16 -20
  73. nucleardatapy/fig/nuc_setupRnpExp_fig.py +12 -0
  74. nucleardatapy/fig/nuc_setupRnpTheo_fig.py +14 -2
  75. nucleardatapy/hnuc/setup_re2L_exp.py +30 -7
  76. nucleardatapy/matter/setup_micro.py +424 -12
  77. nucleardatapy/matter/setup_micro_esym.py +46 -42
  78. nucleardatapy/matter/setup_nep.py +1 -1
  79. nucleardatapy/matter/setup_pheno.py +10 -6
  80. nucleardatapy/matter/setup_pheno_esym.py +14 -6
  81. nucleardatapy/nuc/setup_be_exp.py +4 -2
  82. nucleardatapy/nuc/setup_be_theo.py +19 -12
  83. {nucleardatapy-1.0.0.dist-info → nucleardatapy-1.0.1.dist-info}/METADATA +1 -1
  84. {nucleardatapy-1.0.0.dist-info → nucleardatapy-1.0.1.dist-info}/RECORD +87 -61
  85. {nucleardatapy-1.0.0.dist-info → nucleardatapy-1.0.1.dist-info}/WHEEL +0 -0
  86. {nucleardatapy-1.0.0.dist-info → nucleardatapy-1.0.1.dist-info}/licenses/LICENSE +0 -0
  87. {nucleardatapy-1.0.0.dist-info → nucleardatapy-1.0.1.dist-info}/top_level.txt +0 -0
@@ -13,7 +13,7 @@ def mr_sources():
13
13
  #
14
14
  if nuda.env.verb: print("\nEnter astro_mr()")
15
15
  #
16
- sources = [ 'J0030+0451', 'J0740+6620', 'J0437-4715' ]
16
+ sources = [ 'J0030+0451', 'J0740+6620', 'J0437-4715', 'J0614-3329', 'J1731-347' ]
17
17
  #
18
18
  #print('sources available in the toolkit:',sources)
19
19
  sources_lower = [ item.lower() for item in sources ]
@@ -42,6 +42,10 @@ def mr_obss( source ):
42
42
  obss = [ 1, 2, 3 ]
43
43
  elif source.lower()=='j0437-4715':
44
44
  obss = [ 1 ]
45
+ elif source.lower()=='j0614-3329':
46
+ obss = [ 1 ]
47
+ elif source.lower()=='j1731-347':
48
+ obss = [ 1 ]
45
49
  #
46
50
  #print('Observations available in the toolkit:',obss)
47
51
  #
@@ -96,7 +100,7 @@ class setupMR():
96
100
  file_in = nuda.param.path_data+'astro/NICER/J0030+0451.dat'
97
101
  if obs==1:
98
102
  #: Attribute providing the full reference to the paper to be citted.
99
- self.ref='M.C. Miller, F.K. Lamb, A.J. Dittmann, aet al., ApJL 887, L24 (2019).'
103
+ self.ref='M.C. Miller, F.K. Lamb, A.J. Dittmann, aet al., ApJL 887, L24 (2019)'
100
104
  #: Attribute providing the label the data is references for figures.
101
105
  self.label = 'J0030 Miller 2019'
102
106
  #: Attribute providing additional notes about the observation.
@@ -104,7 +108,7 @@ class setupMR():
104
108
  self.marker = 'o'
105
109
  elif obs==2:
106
110
  #: Attribute providing the full reference to the paper to be citted.
107
- self.ref='T.E. Riley, A.L. Watts, S. Bogdanov, P.S. Ray, et al., ApJ 887, L21 (2019).'
111
+ self.ref='T.E. Riley, A.L. Watts, S. Bogdanov, P.S. Ray, et al., ApJ 887, L21 (2019)'
108
112
  #: Attribute providing the label the data is references for figures.
109
113
  self.label = 'J0030 Riley 2019'
110
114
  #: Attribute providing additional notes about the observation.
@@ -112,7 +116,7 @@ class setupMR():
112
116
  self.marker = 's'
113
117
  elif obs==3:
114
118
  #: Attribute providing the full reference to the paper to be citted.
115
- self.ref='S. Vinciguerra, T. Salmi, A.L. Watts, D. Choudhury, et al., ApJ 961, 62 (2024).'
119
+ self.ref='S. Vinciguerra, T. Salmi, A.L. Watts, D. Choudhury, et al., ApJ 961, 62 (2024)'
116
120
  #: Attribute providing the label the data is references for figures.
117
121
  self.label = 'J0030 Vinciguerra 2024a'
118
122
  #: Attribute providing additional notes about the observation.
@@ -120,7 +124,7 @@ class setupMR():
120
124
  self.marker = 'x'
121
125
  elif obs==4:
122
126
  #: Attribute providing the full reference to the paper to be citted.
123
- self.ref='S. Vinciguerra, T. Salmi, A.L. Watts, D. Choudhury, et al., ApJ 961, 62 (2024).'
127
+ self.ref='S. Vinciguerra, T. Salmi, A.L. Watts, D. Choudhury, et al., ApJ 961, 62 (2024)'
124
128
  #: Attribute providing the label the data is references for figures.
125
129
  self.label = 'J0030 Vinciguerra 2024b'
126
130
  #: Attribute providing additional notes about the observation.
@@ -130,7 +134,7 @@ class setupMR():
130
134
  file_in = nuda.param.path_data+'astro/NICER/J0740+6620.dat'
131
135
  if obs==1:
132
136
  #: Attribute providing the full reference to the paper to be citted.
133
- self.ref='M.C. Miller, F.K. Lamb, A.J. Dittmann, S. Bogdanov, et al., ApJL 918, L28 (2021).'
137
+ self.ref='M.C. Miller, F.K. Lamb, A.J. Dittmann, S. Bogdanov, et al., ApJL 918, L28 (2021)'
134
138
  #: Attribute providing the label the data is references for figures.
135
139
  self.label = 'J0740 Miller 2021'
136
140
  #: Attribute providing additional notes about the observation.
@@ -138,7 +142,7 @@ class setupMR():
138
142
  self.marker = 'o'
139
143
  elif obs==2:
140
144
  #: Attribute providing the full reference to the paper to be citted.
141
- self.ref='T.E. Riley, A.L. Watts, P.S. Ray, S. Bogdanov, et al., ApJL 918, L27 (2021).'
145
+ self.ref='T.E. Riley, A.L. Watts, P.S. Ray, S. Bogdanov, et al., ApJL 918, L27 (2021)'
142
146
  #: Attribute providing the label the data is references for figures.
143
147
  self.label = 'J0740 Riley 2021'
144
148
  #: Attribute providing additional notes about the observation.
@@ -146,7 +150,7 @@ class setupMR():
146
150
  self.marker = 's'
147
151
  elif obs==3:
148
152
  #: Attribute providing the full reference to the paper to be citted.
149
- self.ref='T. Salmi, D. Choudhury, Y. Kini, T.E. Riley et al., ApJ 974, 294 (2024).'
153
+ self.ref='T. Salmi, D. Choudhury, Y. Kini, T.E. Riley et al., ApJ 974, 294 (2024)'
150
154
  #: Attribute providing the label the data is references for figures.
151
155
  self.label = 'J0740 Salmi 2024'
152
156
  #: Attribute providing additional notes about the observation.
@@ -156,12 +160,32 @@ class setupMR():
156
160
  file_in = nuda.param.path_data+'astro/NICER/J0437-4715.dat'
157
161
  if obs==1:
158
162
  #: Attribute providing the full reference to the paper to be citted.
159
- self.ref='D. Choudhury, T. Salmi, S. Vinciguerra, T.E. Riley, et al., ApJL 971, L20 (2024).'
163
+ self.ref='D. Choudhury, T. Salmi, S. Vinciguerra, T.E. Riley, et al., ApJL 971, L20 (2024)'
160
164
  #: Attribute providing the label the data is references for figures.
161
165
  self.label = 'J0437 Choudhury 2024'
162
166
  #: Attribute providing additional notes about the observation.
163
167
  self.note = "write notes about this observation."
164
168
  self.marker = 'o'
169
+ elif source.lower()=='j0614-3329':
170
+ file_in = nuda.param.path_data+'astro/NICER/J0614-3329.dat'
171
+ if obs==1:
172
+ #: Attribute providing the full reference to the paper to be citted.
173
+ self.ref='L. Mauviard, S. Guillot, T. Salmi, D. Choudhury, et al., arXiv:2506.14883'
174
+ #: Attribute providing the label the data is references for figures.
175
+ self.label = 'J0614 Mauviard 2025'
176
+ #: Attribute providing additional notes about the observation.
177
+ self.note = "write notes about this observation."
178
+ self.marker = 'o'
179
+ elif source.lower()=='j1731-347':
180
+ file_in = nuda.param.path_data+'astro/HESS/J1731-347.dat'
181
+ if obs==1:
182
+ #: Attribute providing the full reference to the paper to be citted.
183
+ self.ref='V. Doroshenko, V. Suleimanov, G. P\"uhlhofer and A. Santangelo, Nature Astronomy volume 6, pages 1444–1451 (2022)'
184
+ #: Attribute providing the label the data is references for figures.
185
+ self.label = 'J1731 HESS 2022'
186
+ #: Attribute providing additional notes about the observation.
187
+ self.note = "write notes about this observation."
188
+ self.marker = 's'
165
189
  #
166
190
  #: Attribute the observational mass of the source.
167
191
  self.mass = None
@@ -184,24 +208,49 @@ class setupMR():
184
208
  for line in file:
185
209
  if '#' in line: continue
186
210
  ele = line.split(',')
187
- #print('ele[0]:',ele[0],' obs:',obs,' ele[:]:',ele[:])
211
+ #print('ele[0]:',ele[0],' obs:',obs)
212
+ #print('ele[:]:',ele[:])
188
213
  if int(ele[0]) == obs:
189
- self.rad = float(ele[1])
190
- self.rad_sig_up = float(ele[2])
191
- self.rad_sig_lo = float(ele[3])
192
- self.mass = float(ele[4])
193
- self.mass_sig_up = float(ele[5])
194
- self.mass_sig_lo = float(ele[6])
195
- #print('ele[7]:',ele[7],' len:',len(ele[7]))
196
- if len(ele[7])>1:
214
+ try:
215
+ self.rad = float(ele[1])
216
+ self.rad_sig_up = float(ele[2])
217
+ self.rad_sig_lo = float(ele[3])
218
+ except ValueError:
219
+ self.rad = 0.0
220
+ self.rad_sig_up = 0.0
221
+ self.rad_sig_lo = 0.0
222
+ try:
223
+ self.mass = float(ele[4])
224
+ self.mass_sig_up = float(ele[5])
225
+ self.mass_sig_lo = float(ele[6])
226
+ except ValueError:
227
+ self.mass = 0.0
228
+ self.mass_sig_up = 0.0
229
+ self.mass_sig_lo = 0.0
230
+ #print('ele[7]:',ele[7],isinstance(ele[7].replace(' ',''), (int, float)),' len:',len(ele[7]))
231
+ #print('ele[10]:',ele[10],isinstance(ele[10].replace(' ',''), (int, float)))
232
+ try:
197
233
  self.comp = float(ele[7])
198
234
  self.comp_sig_up = float(ele[8])
199
235
  self.comp_sig_lo = float(ele[9])
200
- else:
236
+ except ValueError:
201
237
  self.comp = 0.0
202
- self.comp_sig_up = 1.0
203
- self.comp_sig_lo = 1.0
204
- self.latexCite = ele[10].replace('\n','').replace(' ','')
238
+ self.comp_sig_up = 0.0
239
+ self.comp_sig_lo = 0.0
240
+ try:
241
+ self.spin = float(ele[10])
242
+ except ValueError:
243
+ self.spin = 0.0
244
+ # if len(ele[7])>1:
245
+ # self.comp = float(ele[7])
246
+ # self.comp_sig_up = float(ele[8])
247
+ # self.comp_sig_lo = float(ele[9])
248
+ # else:
249
+ # self.comp = 0.0
250
+ # self.comp_sig_up = 1.0
251
+ # self.comp_sig_lo = 1.0
252
+ self.latexCite = ele[11].replace('\n','').replace(' ','')
253
+ #print('outputs:',self.rad,self.mass,self.comp,self.spin,self.latexCite)
205
254
  #
206
255
  # compute compactness
207
256
  #
@@ -318,8 +367,8 @@ class setupMRAverage():
318
367
  # construct the distribution of observations in ay
319
368
  ar = np.linspace(rmin,rmax,300); ar=np.array( ar )
320
369
  am = np.linspace(mmin,mmax,300); am=np.array( am )
321
- print('cmax:',cmax)
322
- print('cmin:',cmin)
370
+ #print('cmax:',cmax)
371
+ #print('cmin:',cmin)
323
372
  ac = np.linspace(cmin,cmax,300); ac=np.array( ac )
324
373
  #ac = 0.5 * nuda.cst.rshsol_si / 1.e3 * am / ar
325
374
  ayr = np.zeros(300); ayr=np.array( ayr )
@@ -341,12 +390,36 @@ class setupMRAverage():
341
390
  stdr = sum ( ayr*ar**2 )
342
391
  stdm = sum ( aym*am**2 )
343
392
  stdc = sum ( ayc*ac**2 )
344
- self.rad_cen = cenr / nor
345
- self.mass_cen = cenm / nom
346
- self.comp_cen = cenc / noc
347
- self.rad_sig_std = round( math.sqrt( stdr/nor - self.rad_cen**2 ), 3 )
348
- self.mass_sig_std = round( math.sqrt( stdm/nom - self.mass_cen**2 ), 3 )
349
- self.comp_sig_std = round( math.sqrt( stdc/noc - self.comp_cen**2 ), 3 )
393
+ try:
394
+ if nor != 0.0:
395
+ self.rad_cen = cenr / nor
396
+ self.rad_sig_std = round( math.sqrt( stdr/nor - self.rad_cen**2 ), 3 )
397
+ else:
398
+ self.rad_cen = 0.0
399
+ self.rad_sig_std = 0.0
400
+ except ZeroDivisionError:
401
+ self.rad_cen = 0.0
402
+ self.rad_sig_std = 0.0
403
+ try:
404
+ if nom != 0.0:
405
+ self.mass_cen = cenm / nom
406
+ self.mass_sig_std = round( math.sqrt( stdm/nom - self.mass_cen**2 ), 3 )
407
+ else:
408
+ self.mass_cen = 0.0
409
+ self.mass_sig_std = 0.0
410
+ except ZeroDivisionError:
411
+ self.mass_cen = 0.0
412
+ self.mass_sig_std = 0.0
413
+ try:
414
+ if noc != 0.0:
415
+ self.comp_cen = cenc / noc
416
+ self.comp_sig_std = round( math.sqrt( stdc/noc - self.comp_cen**2 ), 3 )
417
+ else:
418
+ self.comp_cen = 0.0
419
+ self.comp_sig_std = 0.0
420
+ except ZeroDivisionError:
421
+ self.comp_cen = 0.0
422
+ self.comp_sig_std = 0.0
350
423
  self.rad_cen = round( self.rad_cen, 3)
351
424
  self.mass_cen = round( self.mass_cen, 3)
352
425
  self.comp_cen = round( self.comp_cen, 3)
@@ -397,11 +470,20 @@ def gauss( ax, cent, sig_up, sig_lo ):
397
470
  gauss = []
398
471
  for x in ax:
399
472
  if x < cent:
400
- z = ( x - cent ) / sig_lo
473
+ if sig_lo != 0.0:
474
+ z = ( x - cent ) / sig_lo
475
+ else:
476
+ z= x * 0.0
401
477
  norm = sig_lo * fac
402
478
  else:
403
- z = ( x - cent ) / sig_up
479
+ if sig_up != 0.0:
480
+ z = ( x - cent ) / sig_up
481
+ else:
482
+ z= x * 0.0
404
483
  norm = sig_up * fac
405
- gauss.append( math.exp( -0.5*z**2 ) / norm )
484
+ if norm != 0.0:
485
+ gauss.append( math.exp( -0.5*z**2 ) / norm )
486
+ else:
487
+ gauss.append( 0.0 )
406
488
  return gauss
407
489
 
@@ -1,7 +1,7 @@
1
1
 
2
2
  import os
3
3
 
4
- def create_folder_fig():
4
+ def create_folder_fig( folder = 'figs'):
5
5
  #
6
- os.system('mkdir -p figs/')
6
+ os.system('mkdir -p '+folder+'/')
7
7
  #
@@ -176,7 +176,7 @@ class setupCrust():
176
176
  self.label = 'PCPFDDG-BSK22-2018'
177
177
  # NEP
178
178
  nep = nuda.matter.setupNEP( model = 'ESkyrme', param = 'BSk22' )
179
- self.nsat = nep.sat
179
+ self.nsat = nep.nsat
180
180
  self.Esym = nep.Esym
181
181
  self.Lsym = nep.Lsym
182
182
  self.Ksym = nep.Ksym
@@ -193,7 +193,7 @@ class setupCrust():
193
193
  self.label = 'PCPFDDG-BSK24-2018'
194
194
  # NEP
195
195
  nep = nuda.matter.setupNEP( model = 'ESkyrme', param = 'BSk24' )
196
- self.nsat = nep.sat
196
+ self.nsat = nep.nsat
197
197
  self.Esym = nep.Esym
198
198
  self.Lsym = nep.Lsym
199
199
  self.Ksym = nep.Ksym
@@ -210,7 +210,7 @@ class setupCrust():
210
210
  self.label = 'PCPFDDG-BSK25-2018'
211
211
  # NEP
212
212
  nep = nuda.matter.setupNEP( model = 'ESkyrme', param = 'BSk25' )
213
- self.nsat = nep.sat
213
+ self.nsat = nep.nsat
214
214
  self.Esym = nep.Esym
215
215
  self.Lsym = nep.Lsym
216
216
  self.Ksym = nep.Ksym
@@ -227,7 +227,7 @@ class setupCrust():
227
227
  self.label = 'PCPFDDG-BSK26-2018'
228
228
  # NEP
229
229
  nep = nuda.matter.setupNEP( model = 'ESkyrme', param = 'BSk26' )
230
- self.nsat = nep.sat
230
+ self.nsat = nep.nsat
231
231
  self.Esym = nep.Esym
232
232
  self.Lsym = nep.Lsym
233
233
  self.Ksym = nep.Ksym
@@ -659,7 +659,7 @@ class setupCrust():
659
659
  if nuda.env.verb: print("Enter print_outputs()")
660
660
  #
661
661
  print("- Print output:")
662
- print(" modcrust:",self.modcrust)
662
+ print(" model:",self.model)
663
663
  print(" ref: ",self.ref)
664
664
  print(" label: ",self.label)
665
665
  print(" note: ",self.note)
@@ -0,0 +1,4 @@
1
+ # Victor Doroshenko, Valery Suleimanov, Gerd Pühlhofer & Andrea Santangelo, Nature Astronomy volume 6, pages 1444–1451 (2022)
2
+ # obs, R_e, sig+, sig-, M, sig+, sig-, C, sig+, sig-, spin, latexCite
3
+ # km, km, km, Mo, Mo, Mo, -, -, -, Hz
4
+ 1, 10.40, 0.86, 0.78, 0.77, 0.20, 0.17, , , , , VDoroshenko:2022
@@ -1,6 +1,6 @@
1
- # obs, R_e, sig+, sig-, M, sig+, sig-, C, sig+, sig-, spin, dist, latexCite
2
- # km, km, km, Mo, Mo, Mo, -, -, -, Hz, pc
3
- 1, 13.02, 1.24, 1.06, 1.440, 0.150, 0.140, 0.163, 0.008, 0.009, , , MCMiller:2019
4
- 2, 12.71, 1.14, 1.19, 1.340, 0.150, 0.160, 0.156, 0.008, 0.010, , , TERiley:2019
5
- 3, 14.44, 0.88, 1.05, 1.700, 0.180, 0.190, 0.179, 0.011, 0.022, , , SVinciguerra:2024
6
- 4, 11.71, 0.88, 0.83, 1.400, 0.130, 0.120, 0.1773, 0.0056, 0.0074, , , SVinciguerra:2024
1
+ # obs, R_e, sig+, sig-, M, sig+, sig-, C, sig+, sig-, spin, latexCite
2
+ # km, km, km, Mo, Mo, Mo, -, -, -, Hz
3
+ 1, 13.02, 1.24, 1.06, 1.440, 0.150, 0.140, 0.1630, 0.008, 0.009, , MCMiller:2019
4
+ 2, 12.71, 1.14, 1.19, 1.340, 0.150, 0.160, 0.1560, 0.008, 0.010, , TERiley:2019
5
+ 3, 14.44, 0.88, 1.05, 1.700, 0.180, 0.190, 0.1790, 0.011, 0.022, , SVinciguerra:2024
6
+ 4, 11.71, 0.88, 0.83, 1.400, 0.130, 0.120, 0.1773, 0.0056, 0.0074, , SVinciguerra:2024
@@ -1,3 +1,4 @@
1
- # obs, R_e, sig+, sig-, M, sig+, sig-, C, sig+, sig-, latexCite
2
- #
3
- 1, 11.360, 0.945, 0.629, 1.418, 0.037, 0.037, 0.1847, 0.0097, 0.0143, DChoudhury:2024
1
+ # obs, R_e, sig+, sig-, M, sig+, sig-, C, sig+, sig-, spin, latexCite
2
+ # km, km, km, Mo, Mo, Mo, -, -, -, Hz
3
+ 1, 11.360, 0.945, 0.629, 1.418, 0.037, 0.037, 0.1847, 0.0097, 0.0143, , DChoudhury:2024
4
+ #2, 12.800, 0.800, 0.800, PStammler:2025
@@ -0,0 +1,4 @@
1
+ #Lucien Mauviard, Sebastien Guillot, Tuomo Salmi, Devarshi Choudhury, et al., https://doi.org/10.48550/arXiv.2506.14883
2
+ # obs, R_e, sig+, sig-, M, sig+, sig-, C, sig+, sig-, spin, latexCite
3
+ # km, km, km, Mo, Mo, Mo, -, -, -, Hz
4
+ 1, 10.29, 1.01, 0.86, 1.44, 0.06, 0.07, 0.21, 0.02, 0.02, 318, LMauviard:2025
@@ -1,5 +1,5 @@
1
- # obs, R_e, sig+, sig-, M, sig+, sig-, C, sig+, sig-, latexCite
2
- #
3
- 1, 13.70, 2.60, 1.50, 2.080, 0.070, 0.070, , , , MCMiller:2021
4
- 2, 12.39, 1.30, 0.98, 2.072, 0.067, 0.066, , , , TERiley:2021
5
- 3, 12.49, 1.28, 0.88, 2.073, 0.069, 0.069, , , , TSalmi:2024
1
+ # obs, R_e, sig+, sig-, M, sig+, sig-, C, sig+, sig-, spin, latexCite
2
+ # km, km, km, Mo, Mo, Mo, -, -, -, Hz
3
+ 1, 13.70, 2.60, 1.50, 2.080, 0.070, 0.070, , , , , MCMiller:2021
4
+ 2, 12.39, 1.30, 0.98, 2.072, 0.067, 0.066, , , , , TERiley:2021
5
+ 3, 12.49, 1.28, 0.88, 2.073, 0.069, 0.069, , , , , TSalmi:2024
@@ -0,0 +1,6 @@
1
+ # Binding and bound energies for LL-hypernuclei
2
+ # Ref. Yasuo Yamamoto, Hideo Takaki, Kiyomi Ikeda, Progress of Theoretical Physics, Volume 86, Issue 4, October 1991, Pages 867–875, https://doi.org/10.1143/ptp/86.4.867
3
+ # Z, name, N, LLBE, err, LLDBE, err, probe
4
+ # (MeV), (MeV), (MeV), (MeV),
5
+ 4, Be, 4, 8.50, 0.70, -4.90, 0.70, Emul KEK-E176 1991
6
+ 3, B, 8, 27.60, 0.70, 4.80, 0.70, Emul KEK-E176 1991
@@ -2,4 +2,4 @@
2
2
  # Ref. J. K. Ahn, H. Akikawa, S. Aoki, K. Arai, Phys. Rev. C 88, 014003 (2013)
3
3
  # Z, name, N, LLBE, err, LLDBE, err, probe
4
4
  # (MeV), (MeV), (MeV), (MeV),
5
- 2, He, 2, 6.91, 0.16, 0.67, 0.17, Emul
5
+ 2, He, 2, 6.91, 0.16, 0.67, 0.17, Emul NAGARA KEK-E373 2013
@@ -0,0 +1,7 @@
1
+ # Binding and bound energies for LL-hypernuclei
2
+ # Ref. H. Ekawa et al., PTEP 2019, 021D02.
3
+ # Z, name, N, LLBE, err, LLDBE, err, probe
4
+ # (MeV), (MeV), (MeV), (MeV),
5
+ 4, Be, 4, 15.05, 0.11, 1.63, 0.11, (K-,K+) JPARC-E07 2019
6
+ 4, Be, 5, 19.07, 0.11, 1.87, 0.37, (K-,K+) JPARC-E07 2019
7
+ 4, Be, 6, 13.68, 0.11, -2.70, 1.00, (K-,K+) JPARC-E07 2019
@@ -0,0 +1,21 @@
1
+ # rho [fm-3] T [MeV] E_TOT/A [MeV] E_KIN/A [MeV] E_POT/A [MeV]
2
+ 0.20000000E-01 0.00000000E+00 0.41913002E+01 0.96698074E+01 -0.54785072E+01 0.16261796E+00 0.22002090E+00
3
+ 0.30000000E-01 0.00000000E+00 0.50804628E+01 0.12412858E+02 -0.73323955E+01 0.34925549E+00 0.27600771E+00
4
+ 0.32000000E-01 0.00000000E+00 0.54590146E+01 0.12978051E+02 -0.75190368E+01 0.14931802E+00 0.22040522E+00
5
+ 0.40000000E-01 0.00000000E+00 0.62078998E+01 0.14976306E+02 -0.87684063E+01 0.84105694E-01 0.17034625E+00
6
+ 0.60000000E-01 0.00000000E+00 0.73788084E+01 0.19313598E+02 -0.11934789E+02 0.43566048E+00 0.12482157E+00
7
+ 0.64000000E-01 0.00000000E+00 0.79906531E+01 0.20129133E+02 -0.12138480E+02 0.19335889E+00 0.24265078E+00
8
+ 0.80000000E-01 0.00000000E+00 0.92486022E+01 0.23336970E+02 -0.14088367E+02 0.54483190E-01 0.92736911E-04
9
+ 0.96000000E-01 0.00000000E+00 0.10493915E+02 0.26324610E+02 -0.15830695E+02 0.64101603E+00 0.94174971E+00
10
+ 0.10000000E+00 0.00000000E+00 0.10624562E+02 0.27026693E+02 -0.16402131E+02 0.12922241E+00 0.51281243E-01
11
+ 0.12000000E+00 0.00000000E+00 0.11564460E+02 0.30460376E+02 -0.18895916E+02 0.75365783E+00 0.21570968E-01
12
+ 0.12800000E+00 0.00000000E+00 0.12258865E+02 0.31877289E+02 -0.19618424E+02 0.70786033E+00 0.74959546E-01
13
+ 0.14000000E+00 0.00000000E+00 0.13861600E+02 0.33635175E+02 -0.19773575E+02 0.10466647E+00 0.68182628E-01
14
+ 0.16000000E+00 0.00000000E+00 0.16039735E+02 0.36699867E+02 -0.20660131E+02 0.16905664E+00 0.14588913E+00
15
+ 0.18000000E+00 0.00000000E+00 0.17825303E+02 0.39716876E+02 -0.21891573E+02 0.15743319E+00 0.31835316E-01
16
+ 0.20000000E+00 0.00000000E+00 0.20174044E+02 0.42621598E+02 -0.22447554E+02 0.11762932E+00 0.31570418E-02
17
+ 0.24000000E+00 0.00000000E+00 0.25723600E+02 0.48558482E+02 -0.22834882E+02 0.16525810E+00 0.15048843E+00
18
+ 0.28000000E+00 0.00000000E+00 0.31679628E+02 0.54393180E+02 -0.22713552E+02 0.62328270E+00 0.66160029E+00
19
+ 0.32000000E+00 0.00000000E+00 0.37163514E+02 0.59968721E+02 -0.22805207E+02 0.28225974E+00 0.59007041E+00
20
+ 0.36000000E+00 0.00000000E+00 0.42023488E+02 0.63665078E+02 -0.21641591E+02 0.36277403E+00 0.32459982E-01
21
+ 0.40000000E+00 0.00000000E+00 0.52185556E+02 0.73350862E+02 -0.21165306E+02 0.26235105E+01 0.23113296E+01
@@ -0,0 +1,22 @@
1
+ # rho [fm-3] T [MeV] mu [MeV] E/A [MeV] S_1/A F/A [MeV] P [MeV fm-3]
2
+ 0.10000000E-01 0.00000000E+00 0.43668569E+01 0.25987421E+01 0.00000000E+00 0.25987421E+01 0.17681148E-01
3
+ 0.15000000E-01 0.00000000E+00 0.53870389E+01 0.35346218E+01 0.00000000E+00 0.35346218E+01 0.27786257E-01
4
+ 0.20000000E-01 0.00000000E+00 0.62516898E+01 0.41439137E+01 0.00000000E+00 0.41439137E+01 0.42155523E-01
5
+ 0.30000000E-01 0.00000000E+00 0.80928725E+01 0.49814102E+01 0.00000000E+00 0.49814102E+01 0.93343869E-01
6
+ 0.32000000E-01 0.00000000E+00 0.84155745E+01 0.55026567E+01 0.00000000E+00 0.55026567E+01 0.93213368E-01
7
+ 0.40000000E-01 0.00000000E+00 0.96185225E+01 0.60009747E+01 0.00000000E+00 0.60009747E+01 0.14470191E+00
8
+ 0.60000000E-01 0.00000000E+00 0.12373721E+02 0.73117973E+01 0.00000000E+00 0.73117973E+01 0.30371545E+00
9
+ 0.64000000E-01 0.00000000E+00 0.12915035E+02 0.81341882E+01 0.00000000E+00 0.81341882E+01 0.30597421E+00
10
+ 0.80000000E-01 0.00000000E+00 0.15446115E+02 0.92267598E+01 0.00000000E+00 0.92267598E+01 0.49754845E+00
11
+ 0.96000000E-01 0.00000000E+00 0.18222028E+02 0.10431012E+02 0.00000000E+00 0.10431012E+02 0.74793758E+00
12
+ 0.10000000E+00 0.00000000E+00 0.19097157E+02 0.10592070E+02 0.00000000E+00 0.10592070E+02 0.85050872E+00
13
+ 0.12000000E+00 0.00000000E+00 0.23361307E+02 0.11557871E+02 0.00000000E+00 0.11557871E+02 0.14164123E+01
14
+ 0.12800000E+00 0.00000000E+00 0.25287934E+02 0.12254644E+02 0.00000000E+00 0.12254644E+02 0.16682611E+01
15
+ 0.14000000E+00 0.00000000E+00 0.28002552E+02 0.13861600E+02 0.00000000E+00 0.13861600E+02 0.19797333E+01
16
+ 0.16000000E+00 0.00000000E+00 0.33316634E+02 0.16039735E+02 0.00000000E+00 0.16039735E+02 0.27643038E+01
17
+ 0.18000000E+00 0.00000000E+00 0.39373368E+02 0.17825303E+02 0.00000000E+00 0.17825303E+02 0.38786517E+01
18
+ 0.20000000E+00 0.00000000E+00 0.45614955E+02 0.20174044E+02 0.00000000E+00 0.20174044E+02 0.50881822E+01
19
+ 0.24000000E+00 0.00000000E+00 0.58459676E+02 0.25723600E+02 0.00000000E+00 0.25723600E+02 0.78566584E+01
20
+ 0.28000000E+00 0.00000000E+00 0.71623658E+02 0.31679628E+02 0.00000000E+00 0.31679628E+02 0.11184328E+02
21
+ 0.32000000E+00 0.00000000E+00 0.83217235E+02 0.37163514E+02 0.00000000E+00 0.37163514E+02 0.14737191E+02
22
+ 0.36000000E+00 0.00000000E+00 0.97520831E+02 0.42023488E+02 0.00000000E+00 0.42023488E+02 0.19979044E+02
@@ -0,0 +1,20 @@
1
+ # rho [fm-3] T [MeV] E_TOT/A [MeV] E_KIN/A [MeV] E_POT/A [MeV]
2
+ 0.20000000E-01 0.00000000E+00 0.41657809E+01 0.96825377E+01 -0.55167568E+01 0.17118476E+00 0.23202281E+00
3
+ 0.30000000E-01 0.00000000E+00 0.50216893E+01 0.12447897E+02 -0.74262080E+01 0.38330311E+00 0.29402864E+00
4
+ 0.32000000E-01 0.00000000E+00 0.53530691E+01 0.13029864E+02 -0.76767950E+01 0.22966292E+00 0.22770640E+00
5
+ 0.40000000E-01 0.00000000E+00 0.64216991E+01 0.15043028E+02 -0.86213289E+01 0.15829508E+00 0.16649694E+00
6
+ 0.60000000E-01 0.00000000E+00 0.73335035E+01 0.19341325E+02 -0.12007822E+02 0.31392981E+00 0.22579378E-01
7
+ 0.64000000E-01 0.00000000E+00 0.77906989E+01 0.20164171E+02 -0.12373472E+02 0.36318291E+00 0.26526393E+00
8
+ 0.80000000E-01 0.00000000E+00 0.92962529E+01 0.23372601E+02 -0.14076348E+02 0.12002559E+00 0.37018091E+00
9
+ 0.96000000E-01 0.00000000E+00 0.10404343E+02 0.26271600E+02 -0.15867258E+02 0.30503521E-01 0.93634841E-01
10
+ 0.10000000E+00 0.00000000E+00 0.10879062E+02 0.27055519E+02 -0.16176457E+02 0.76113627E-01 0.22271940E+00
11
+ 0.12000000E+00 0.00000000E+00 0.12282901E+02 0.30529203E+02 -0.18246303E+02 0.86787541E-01 0.21396893E+00
12
+ 0.12800000E+00 0.00000000E+00 0.12165178E+02 0.31903869E+02 -0.19738691E+02 0.66569231E+00 0.68935687E-01
13
+ 0.14000000E+00 0.00000000E+00 0.13039971E+02 0.33895409E+02 -0.20855438E+02 0.78253098E+00 0.11951527E+00
14
+ 0.16000000E+00 0.00000000E+00 0.14912759E+02 0.36918290E+02 -0.22005531E+02 0.82559068E+00 0.16839848E+00
15
+ 0.18000000E+00 0.00000000E+00 0.17308818E+02 0.39915624E+02 -0.22606805E+02 0.21365357E+00 0.35554315E-01
16
+ 0.20000000E+00 0.00000000E+00 0.19416739E+02 0.42880607E+02 -0.23463869E+02 0.19032138E+00 0.55451658E-02
17
+ 0.24000000E+00 0.00000000E+00 0.22366972E+02 0.48367291E+02 -0.26000319E+02 0.14946694E+01 0.19756319E+00
18
+ 0.28000000E+00 0.00000000E+00 0.27982135E+02 0.53905147E+02 -0.25923012E+02 0.40852900E+00 0.12524215E+00
19
+ 0.32000000E+00 0.00000000E+00 0.34157475E+02 0.59182012E+02 -0.25024537E+02 0.58623950E+00 0.78688057E+00
20
+ 0.36000000E+00 0.00000000E+00 0.38107297E+02 0.64443022E+02 -0.26335724E+02 0.35768921E+00 0.94990033E+00
@@ -0,0 +1,22 @@
1
+ # rho [fm-3] T [MeV] mu [MeV] E/A [MeV] S_1/A F/A [MeV] P [MeV fm-3]
2
+ 0.10000000E-01 0.00000000E+00 0.43459235E+01 0.25850924E+01 0.00000000E+00 0.25850924E+01 0.17608311E-01
3
+ 0.15000000E-01 0.00000000E+00 0.53604558E+01 0.35262487E+01 0.00000000E+00 0.35262487E+01 0.27513106E-01
4
+ 0.20000000E-01 0.00000000E+00 0.66968125E+01 0.41362726E+01 0.00000000E+00 0.41362726E+01 0.51210798E-01
5
+ 0.30000000E-01 0.00000000E+00 0.80123887E+01 0.49398604E+01 0.00000000E+00 0.49398604E+01 0.92175849E-01
6
+ 0.32000000E-01 0.00000000E+00 0.83936762E+01 0.53204847E+01 0.00000000E+00 0.53204847E+01 0.98342128E-01
7
+ 0.40000000E-01 0.00000000E+00 0.95646916E+01 0.62913252E+01 0.00000000E+00 0.62913252E+01 0.13093466E+00
8
+ 0.60000000E-01 0.00000000E+00 0.12324462E+02 0.72527659E+01 0.00000000E+00 0.72527659E+01 0.30430176E+00
9
+ 0.64000000E-01 0.00000000E+00 0.12807441E+02 0.77478480E+01 0.00000000E+00 0.77478480E+01 0.32381398E+00
10
+ 0.80000000E-01 0.00000000E+00 0.15272577E+02 0.93173938E+01 0.00000000E+00 0.93173938E+01 0.47641465E+00
11
+ 0.96000000E-01 0.00000000E+00 0.17976245E+02 0.10339071E+02 0.00000000E+00 0.10339071E+02 0.73316870E+00
12
+ 0.10000000E+00 0.00000000E+00 0.18737917E+02 0.10778360E+02 0.00000000E+00 0.10778360E+02 0.79595564E+00
13
+ 0.12000000E+00 0.00000000E+00 0.22701164E+02 0.12293341E+02 0.00000000E+00 0.12293341E+02 0.12489388E+01
14
+ 0.12800000E+00 0.00000000E+00 0.24416021E+02 0.12159705E+02 0.00000000E+00 0.12159705E+02 0.15688084E+01
15
+ 0.14000000E+00 0.00000000E+00 0.27095860E+02 0.13039971E+02 0.00000000E+00 0.13039971E+02 0.19678244E+01
16
+ 0.16000000E+00 0.00000000E+00 0.31489385E+02 0.15487453E+02 0.00000000E+00 0.15487453E+02 0.25603090E+01
17
+ 0.18000000E+00 0.00000000E+00 0.36535720E+02 0.17308818E+02 0.00000000E+00 0.17308818E+02 0.34608423E+01
18
+ 0.20000000E+00 0.00000000E+00 0.41852326E+02 0.19416739E+02 0.00000000E+00 0.19416739E+02 0.44871175E+01
19
+ 0.24000000E+00 0.00000000E+00 0.52198617E+02 0.22366972E+02 0.00000000E+00 0.22366972E+02 0.71595948E+01
20
+ 0.28000000E+00 0.00000000E+00 0.63272413E+02 0.27982135E+02 0.00000000E+00 0.27982135E+02 0.98812777E+01
21
+ 0.32000000E+00 0.00000000E+00 0.73666724E+02 0.34157475E+02 0.00000000E+00 0.34157475E+02 0.12642960E+02
22
+ 0.36000000E+00 0.00000000E+00 0.84131744E+02 0.38107297E+02 0.00000000E+00 0.38107297E+02 0.16568801E+02
@@ -0,0 +1,23 @@
1
+ # rho [fm-3] T [MeV] E_TOT/A [MeV] E_KIN/A [MeV] E_POT/A [MeV]
2
+ 0.10000000E-01 0.00000000E+00 0.25087242E+01 0.66124314E+01 -0.41037072E+01 0.23736791E+00 0.92243562E-01
3
+ 0.20000000E-01 0.00000000E+00 0.41556490E+01 0.10054047E+02 -0.58983975E+01 0.16224846E+00 0.21612753E+00
4
+ 0.30000000E-01 0.00000000E+00 0.50265052E+01 0.13007956E+02 -0.79814504E+01 0.34867092E+00 0.29921062E+00
5
+ 0.32000000E-01 0.00000000E+00 0.54452468E+01 0.13649187E+02 -0.82039404E+01 0.10835225E+00 0.22499772E+00
6
+ 0.40000000E-01 0.00000000E+00 0.60505479E+01 0.15846296E+02 -0.97957481E+01 0.19368868E+00 0.16529118E+00
7
+ 0.60000000E-01 0.00000000E+00 0.73629856E+01 0.20586053E+02 -0.13223068E+02 0.39727356E+00 0.11784186E+00
8
+ 0.64000000E-01 0.00000000E+00 0.81127666E+01 0.21487528E+02 -0.13374761E+02 0.67967444E-01 0.10654290E+00
9
+ 0.80000000E-01 0.00000000E+00 0.91101936E+01 0.25206723E+02 -0.16096529E+02 0.22995657E+00 0.28580600E+00
10
+ 0.96000000E-01 0.00000000E+00 0.10242757E+02 0.28537261E+02 -0.18294504E+02 0.19277562E+00 0.24915610E+00
11
+ 0.10000000E+00 0.00000000E+00 0.10620913E+02 0.29335758E+02 -0.18714845E+02 0.44294932E-01 0.10626409E+00
12
+ 0.12000000E+00 0.00000000E+00 0.11364693E+02 0.33306837E+02 -0.21942144E+02 0.70214354E+00 0.13572460E+00
13
+ 0.12800000E+00 0.00000000E+00 0.13100656E+02 0.34712372E+02 -0.21611716E+02 0.41582521E+00 0.28424456E-01
14
+ 0.14000000E+00 0.00000000E+00 0.13673425E+02 0.36637612E+02 -0.22964186E+02 0.25021458E-02 0.84550741E-01
15
+ 0.16000000E+00 0.00000000E+00 0.14779873E+02 0.40349123E+02 -0.25569250E+02 0.85323381E+00 0.22573325E+00
16
+ 0.18000000E+00 0.00000000E+00 0.17606202E+02 0.43529405E+02 -0.25923203E+02 0.58869062E-02 0.62700487E-01
17
+ 0.20000000E+00 0.00000000E+00 0.19979610E+02 0.47058889E+02 -0.27079279E+02 0.22097598E+00 0.22379158E-02
18
+ 0.24000000E+00 0.00000000E+00 0.24550382E+02 0.53342274E+02 -0.28791892E+02 0.68409435E+00 0.41573452E-01
19
+ 0.28000000E+00 0.00000000E+00 0.30803410E+02 0.59447523E+02 -0.28644112E+02 0.92409806E-01 0.53842228E-01
20
+ 0.32000000E+00 0.00000000E+00 0.36732724E+02 0.65253967E+02 -0.28521243E+02 0.53291433E+00 0.22593595E-01
21
+ 0.36000000E+00 0.00000000E+00 0.42254296E+02 0.70974378E+02 -0.28720083E+02 0.19632918E+01 0.26776891E+00
22
+ 0.40000000E+00 0.00000000E+00 0.52332618E+02 0.76851901E+02 -0.24519283E+02 0.56914386E+00 0.10200168E+00
23
+ 0.44000000E+00 0.00000000E+00 0.59287999E+02 0.80240591E+02 -0.20952591E+02 0.30346902E+00 0.11790481E+01
@@ -0,0 +1,22 @@
1
+ # rho [fm-3] T [MeV] mu [MeV] E/A [MeV] S_1/A F/A [MeV] P [MeV fm-3]
2
+ 0.10000000E-01 0.00000000E+00 0.42737590E+01 0.25087242E+01 0.00000000E+00 0.25087242E+01 0.17650347E-01
3
+ 0.15000000E-01 0.00000000E+00 0.53582801E+01 0.35367341E+01 0.00000000E+00 0.35367341E+01 0.27323189E-01
4
+ 0.20000000E-01 0.00000000E+00 0.60879971E+01 0.42698546E+01 0.00000000E+00 0.42698546E+01 0.36362852E-01
5
+ 0.30000000E-01 0.00000000E+00 0.79670145E+01 0.49348308E+01 0.00000000E+00 0.49348308E+01 0.90965514E-01
6
+ 0.32000000E-01 0.00000000E+00 0.83283411E+01 0.56088532E+01 0.00000000E+00 0.56088532E+01 0.87023611E-01
7
+ 0.40000000E-01 0.00000000E+00 0.95715885E+01 0.58961997E+01 0.00000000E+00 0.58961997E+01 0.14701555E+00
8
+ 0.60000000E-01 0.00000000E+00 0.12356439E+02 0.72937461E+01 0.00000000E+00 0.72937461E+01 0.30376157E+00
9
+ 0.64000000E-01 0.00000000E+00 0.12928223E+02 0.82540299E+01 0.00000000E+00 0.82540299E+01 0.29914836E+00
10
+ 0.80000000E-01 0.00000000E+00 0.15455832E+02 0.89670334E+01 0.00000000E+00 0.89670334E+01 0.51910390E+00
11
+ 0.96000000E-01 0.00000000E+00 0.18129853E+02 0.10239001E+02 0.00000000E+00 0.10239001E+02 0.75752182E+00
12
+ 0.10000000E+00 0.00000000E+00 0.19123871E+02 0.10531601E+02 0.00000000E+00 0.10531601E+02 0.85922700E+00
13
+ 0.12000000E+00 0.00000000E+00 0.23320545E+02 0.11359904E+02 0.00000000E+00 0.11359904E+02 0.14352770E+01
14
+ 0.12800000E+00 0.00000000E+00 0.25032316E+02 0.13102998E+02 0.00000000E+00 0.13102998E+02 0.15269526E+01
15
+ 0.14000000E+00 0.00000000E+00 0.27902356E+02 0.13663117E+02 0.00000000E+00 0.13663117E+02 0.19934935E+01
16
+ 0.16000000E+00 0.00000000E+00 0.32816378E+02 0.14779873E+02 0.00000000E+00 0.14779873E+02 0.28858408E+01
17
+ 0.18000000E+00 0.00000000E+00 0.39281034E+02 0.17606202E+02 0.00000000E+00 0.17606202E+02 0.39014697E+01
18
+ 0.20000000E+00 0.00000000E+00 0.46593462E+02 0.19979610E+02 0.00000000E+00 0.19979610E+02 0.53227705E+01
19
+ 0.24000000E+00 0.00000000E+00 0.61205726E+02 0.24550382E+02 0.00000000E+00 0.24550382E+02 0.87972828E+01
20
+ 0.28000000E+00 0.00000000E+00 0.77496254E+02 0.30803410E+02 0.00000000E+00 0.30803410E+02 0.13073996E+02
21
+ 0.32000000E+00 0.00000000E+00 0.94918693E+02 0.36732724E+02 0.00000000E+00 0.36732724E+02 0.18619510E+02
22
+ 0.36000000E+00 0.00000000E+00 0.11346746E+03 0.42254296E+02 0.00000000E+00 0.42254296E+02 0.25636738E+02
@@ -0,0 +1,15 @@
1
+ # rho [fm-3] T [MeV] E_TOT/A [MeV] E_KIN/A [MeV] E_POT/A [MeV]
2
+ 0.30000000E-01 0.00000000E+00 -0.63894153E+01 0.12334033E+02 -0.18723448E+02 0.21364439E+00 0.31957126E+00
3
+ 0.32000000E-01 0.00000000E+00 -0.65486758E+01 0.12636287E+02 -0.19184963E+02 0.23953226E+00 0.55058904E+00
4
+ 0.40000000E-01 0.00000000E+00 -0.74992084E+01 0.13747550E+02 -0.21246759E+02 0.86093955E-01 0.45037576E+00
5
+ 0.60000000E-01 0.00000000E+00 -0.97569757E+01 0.16062580E+02 -0.25819555E+02 0.35760003E-01 0.49802524E+00
6
+ 0.64000000E-01 0.00000000E+00 -0.10241275E+02 0.16548300E+02 -0.26789575E+02 0.16952119E+00 0.51299416E+00
7
+ 0.80000000E-01 0.00000000E+00 -0.11677112E+02 0.18277687E+02 -0.29954800E+02 0.47661659E+00 0.95259017E+00
8
+ 0.96000000E-01 0.00000000E+00 -0.12619164E+02 0.20130160E+02 -0.32749324E+02 0.20311108E+00 0.22460433E+00
9
+ 0.10000000E+00 0.00000000E+00 -0.12990146E+02 0.20500599E+02 -0.33490745E+02 0.16958006E-01 0.43596073E+00
10
+ 0.12000000E+00 0.00000000E+00 -0.14115592E+02 0.22581546E+02 -0.36697138E+02 0.97231288E-01 0.28350964E+00
11
+ 0.12800000E+00 0.00000000E+00 -0.14616082E+02 0.23428053E+02 -0.38044135E+02 0.83735932E-01 0.35709105E+00
12
+ 0.14000000E+00 0.00000000E+00 -0.15261233E+02 0.24736887E+02 -0.39998120E+02 0.94053543E-01 0.79111467E-01
13
+ 0.16000000E+00 0.00000000E+00 -0.15463196E+02 0.26601463E+02 -0.42064659E+02 0.41579297E-02 0.34295856E+00
14
+ 0.18000000E+00 0.00000000E+00 -0.15815713E+02 0.28605993E+02 -0.44421706E+02 0.14610417E+00 0.13399562E+00
15
+ 0.20000000E+00 0.00000000E+00 -0.15778297E+02 0.30192407E+02 -0.45970705E+02 0.96986608E-01 0.28316980E+00
@@ -0,0 +1,21 @@
1
+ # rho [fm-3] T [MeV] mu [MeV] E/A [MeV] S_1/A F/A [MeV] P [MeV fm-3]
2
+ 0.20000000E-01 0.00000000E+00 -0.93024043E+01 -0.46926017E+01 0.00000000E+00 -0.46926017E+01 -0.92196052E-01
3
+ 0.30000000E-01 0.00000000E+00 -0.10858707E+02 -0.63894153E+01 0.00000000E+00 -0.63894153E+01 -0.13407876E+00
4
+ 0.32000000E-01 0.00000000E+00 -0.11230298E+02 -0.65486758E+01 0.00000000E+00 -0.65486758E+01 -0.14981192E+00
5
+ 0.40000000E-01 0.00000000E+00 -0.12377231E+02 -0.74652548E+01 0.00000000E+00 -0.74652548E+01 -0.19647903E+00
6
+ 0.60000000E-01 0.00000000E+00 -0.15495961E+02 -0.97783856E+01 0.00000000E+00 -0.97783856E+01 -0.34305451E+00
7
+ 0.64000000E-01 0.00000000E+00 -0.15942145E+02 -0.10241275E+02 0.00000000E+00 -0.10241275E+02 -0.36485574E+00
8
+ 0.80000000E-01 0.00000000E+00 -0.17792006E+02 -0.11677112E+02 0.00000000E+00 -0.11677112E+02 -0.48919153E+00
9
+ 0.96000000E-01 0.00000000E+00 -0.18719010E+02 -0.12619164E+02 0.00000000E+00 -0.12619164E+02 -0.58558517E+00
10
+ 0.10000000E+00 0.00000000E+00 -0.18918612E+02 -0.12990146E+02 0.00000000E+00 -0.12990146E+02 -0.59284664E+00
11
+ 0.12000000E+00 0.00000000E+00 -0.19178101E+02 -0.14115592E+02 0.00000000E+00 -0.14115592E+02 -0.60750107E+00
12
+ 0.12800000E+00 0.00000000E+00 -0.18995106E+02 -0.14791775E+02 0.00000000E+00 -0.14791775E+02 -0.53802636E+00
13
+ 0.14000000E+00 0.00000000E+00 -0.18332826E+02 -0.15261233E+02 0.00000000E+00 -0.15261233E+02 -0.43002306E+00
14
+ 0.16000000E+00 0.00000000E+00 -0.16522628E+02 -0.15463196E+02 0.00000000E+00 -0.15463196E+02 -0.16950904E+00
15
+ 0.18000000E+00 0.00000000E+00 -0.13756101E+02 -0.15815713E+02 0.00000000E+00 -0.15815713E+02 0.37073022E+00
16
+ 0.20000000E+00 0.00000000E+00 -0.10631430E+02 -0.15778297E+02 0.00000000E+00 -0.15778297E+02 0.10293735E+01
17
+ 0.24000000E+00 0.00000000E+00 -0.40867975E+00 -0.14709416E+02 0.00000000E+00 -0.14709416E+02 0.34321766E+01
18
+ 0.28000000E+00 0.00000000E+00 0.12834465E+02 -0.12496853E+02 0.00000000E+00 -0.12496853E+02 0.70927689E+01
19
+ 0.32000000E+00 0.00000000E+00 0.28550502E+02 -0.93970777E+01 0.00000000E+00 -0.93970777E+01 0.12143225E+02
20
+ 0.36000000E+00 0.00000000E+00 0.46965244E+02 -0.45623782E+01 0.00000000E+00 -0.45623782E+01 0.18549944E+02
21
+ 0.40000000E+00 0.00000000E+00 0.66036683E+02 0.90440161E+00 0.00000000E+00 0.90440161E+00 0.26052912E+02
@@ -0,0 +1,15 @@
1
+ # rho [fm-3] T [MeV] E_TOT/A [MeV] E_KIN/A [MeV] E_POT/A [MeV]
2
+ 0.30000000E-01 0.00000000E+00 -0.61554390E+01 0.12270185E+02 -0.18425624E+02 0.22012863E+00 0.36262631E+00
3
+ 0.32000000E-01 0.00000000E+00 -0.63383429E+01 0.12840780E+02 -0.19179123E+02 0.28227333E+00 0.41822138E+00
4
+ 0.40000000E-01 0.00000000E+00 -0.74056456E+01 0.13489849E+02 -0.20895495E+02 0.24237224E+00 0.54892751E+00
5
+ 0.60000000E-01 0.00000000E+00 -0.95092348E+01 0.17576348E+02 -0.27085583E+02 0.25843061E-02 0.47064058E+00
6
+ 0.64000000E-01 0.00000000E+00 -0.10038964E+02 0.18307448E+02 -0.28346412E+02 0.18162415E+00 0.45968955E+00
7
+ 0.80000000E-01 0.00000000E+00 -0.11200838E+02 0.20394028E+02 -0.31594866E+02 0.94796151E-01 0.44386266E+00
8
+ 0.96000000E-01 0.00000000E+00 -0.12446458E+02 0.22124216E+02 -0.34570673E+02 0.15313601E+00 0.29405495E+00
9
+ 0.10000000E+00 0.00000000E+00 -0.12339861E+02 0.22447738E+02 -0.34787599E+02 0.15590837E+00 0.33826082E+00
10
+ 0.12000000E+00 0.00000000E+00 -0.13552936E+02 0.24669178E+02 -0.38222115E+02 0.96008434E-01 0.12364333E+00
11
+ 0.12800000E+00 0.00000000E+00 -0.13804676E+02 0.25446200E+02 -0.39250876E+02 0.60844602E-01 0.22899078E+00
12
+ 0.14000000E+00 0.00000000E+00 -0.14193752E+02 0.26729830E+02 -0.40923583E+02 0.11089745E+00 0.14869100E+00
13
+ 0.16000000E+00 0.00000000E+00 -0.14314444E+02 0.28739481E+02 -0.43053925E+02 0.22839281E+00 0.23091399E+00
14
+ 0.18000000E+00 0.00000000E+00 -0.14764287E+02 0.30761830E+02 -0.45526116E+02 0.98801311E-01 0.13092863E+00
15
+ 0.20000000E+00 0.00000000E+00 -0.14514000E+02 0.32656401E+02 -0.47170401E+02 0.11630126E+00 0.12260756E+00
@@ -0,0 +1,21 @@
1
+ # rho [fm-3] T [MeV] mu [MeV] E/A [MeV] S_1/A F/A [MeV] P [MeV fm-3]
2
+ 2.00000000E-02 0.00000000E+00 -9.09551272E+00 -4.75059173E+00 0.00000000E+00 -4.75059173E+00 -8.68984199E-02
3
+ 3.00000000E-02 0.00000000E+00 -1.06038116E+01 -6.18827908E+00 0.00000000E+00 -6.18827908E+00 -1.32465976E-01
4
+ 3.20000000E-02 0.00000000E+00 -1.08837002E+01 -6.33179290E+00 0.00000000E+00 -6.33179290E+00 -1.45661034E-01
5
+ 4.00000000E-02 0.00000000E+00 -1.20941177E+01 -7.45275322E+00 0.00000000E+00 -7.45275322E+00 -1.85654577E-01
6
+ 6.00000000E-02 0.00000000E+00 -1.50784008E+01 -9.54709529E+00 0.00000000E+00 -9.54709529E+00 -3.31878333E-01
7
+ 6.40000000E-02 0.00000000E+00 -1.54055026E+01 -1.00445255E+01 0.00000000E+00 -1.00445255E+01 -3.43102540E-01
8
+ 8.00000000E-02 0.00000000E+00 -1.67417643E+01 -1.11523177E+01 0.00000000E+00 -1.11523177E+01 -4.47155729E-01
9
+ 9.60000000E-02 0.00000000E+00 -1.74821049E+01 -1.24014462E+01 0.00000000E+00 -1.24014462E+01 -4.87743239E-01
10
+ 1.00000000E-01 0.00000000E+00 -1.76440833E+01 -1.23085835E+01 0.00000000E+00 -1.23085835E+01 -5.33549975E-01
11
+ 1.20000000E-01 0.00000000E+00 -1.72308912E+01 -1.34921954E+01 0.00000000E+00 -1.34921954E+01 -4.48643494E-01
12
+ 1.28000000E-01 0.00000000E+00 -1.68165974E+01 -1.37555281E+01 0.00000000E+00 -1.37555281E+01 -3.91816864E-01
13
+ 1.40000000E-01 0.00000000E+00 -1.58770992E+01 -1.41427189E+01 0.00000000E+00 -1.41427189E+01 -2.42813240E-01
14
+ 1.60000000E-01 0.00000000E+00 -1.34489381E+01 -1.42713196E+01 0.00000000E+00 -1.42713196E+01 1.31581042E-01
15
+ 1.80000000E-01 0.00000000E+00 -1.00410752E+01 -1.47247802E+01 0.00000000E+00 -1.47247802E+01 8.43066900E-01
16
+ 2.00000000E-01 0.00000000E+00 -6.03304636E+00 -1.45310037E+01 0.00000000E+00 -1.45310037E+01 1.69959146E+00
17
+ 2.40000000E-01 0.00000000E+00 5.34131193E+00 -1.31356352E+01 0.00000000E+00 -1.31356352E+01 4.43446730E+00
18
+ 2.80000000E-01 0.00000000E+00 1.94422181E+01 -1.06195309E+01 0.00000000E+00 -1.06195309E+01 8.41728973E+00
19
+ 3.20000000E-01 0.00000000E+00 3.62030849E+01 -6.22439064E+00 0.00000000E+00 -6.22439064E+00 1.35767922E+01
20
+ 3.60000000E-01 0.00000000E+00 5.40294307E+01 -1.50586342E+00 0.00000000E+00 -1.50586342E+00 1.99927059E+01
21
+ 4.00000000E-01 0.00000000E+00 7.32345740E+01 4.19037193E+00 0.00000000E+00 4.19037193E+00 2.76176808E+01