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.
- nucleardatapy/astro/setup_mr.py +115 -33
- nucleardatapy/create_folder.py +2 -2
- nucleardatapy/crust/setup_crust.py +5 -5
- nucleardatapy/data/astro/HESS/J1731-347.dat +4 -0
- nucleardatapy/data/astro/NICER/J0030+0451.dat +6 -6
- nucleardatapy/data/astro/NICER/J0437-4715.dat +4 -3
- nucleardatapy/data/astro/NICER/J0614-3329.dat +4 -0
- nucleardatapy/data/astro/NICER/J0740+6620.dat +5 -5
- nucleardatapy/data/hnuclei/1991-2L-Yamamoto.csv +6 -0
- nucleardatapy/data/hnuclei/2013-2L-Ahn.csv +1 -1
- nucleardatapy/data/hnuclei/2019-2L-Ekawa.csv +7 -0
- nucleardatapy/data/matter/micro/2020-SCGF-NM-N3LO-414-E2A.dat +21 -0
- nucleardatapy/data/matter/micro/2020-SCGF-NM-N3LO-414-TD.dat +22 -0
- nucleardatapy/data/matter/micro/2020-SCGF-NM-N3LO-450-E2A.dat +20 -0
- nucleardatapy/data/matter/micro/2020-SCGF-NM-N3LO-450-TD.dat +22 -0
- nucleardatapy/data/matter/micro/2020-SCGF-NM-N3LO-500-E2A.dat +23 -0
- nucleardatapy/data/matter/micro/2020-SCGF-NM-N3LO-500-TD.dat +22 -0
- nucleardatapy/data/matter/micro/2020-SCGF-SM-N3LO-414-E2A.dat +15 -0
- nucleardatapy/data/matter/micro/2020-SCGF-SM-N3LO-414-TD.dat +21 -0
- nucleardatapy/data/matter/micro/2020-SCGF-SM-N3LO-450-E2A.dat +15 -0
- nucleardatapy/data/matter/micro/2020-SCGF-SM-N3LO-450-TD.dat +21 -0
- nucleardatapy/data/matter/micro/2020-SCGF-SM-N3LO-500-E2A.dat +20 -0
- nucleardatapy/data/matter/micro/2020-SCGF-SM-N3LO-500-TD.dat +20 -0
- nucleardatapy/data/matter/micro/2024-ABI-NM-DeltaNNLO450.dat +28 -0
- nucleardatapy/data/matter/micro/2024-ABI-NM-DeltaNNLO500.dat +28 -0
- nucleardatapy/data/matter/micro/2024-ABI-NM-DeltaNNLOgo394.dat +28 -0
- nucleardatapy/data/matter/micro/2024-ABI-NM-DeltaNNLOgo450.dat +28 -0
- nucleardatapy/data/matter/micro/2024-ABI-NM-NNLOsat.dat +28 -0
- nucleardatapy/data/matter/micro/2024-ABI-SM-DeltaNNLO450.dat +28 -0
- nucleardatapy/data/matter/micro/2024-ABI-SM-DeltaNNLO500.dat +28 -0
- nucleardatapy/data/matter/micro/2024-ABI-SM-DeltaNNLOgo394.dat +28 -0
- nucleardatapy/data/matter/micro/2024-ABI-SM-DeltaNNLOgo450.dat +28 -0
- nucleardatapy/data/matter/micro/2024-ABI-SM-NNLOsat.dat +28 -0
- nucleardatapy/fig/astro_setupGW_fig.py +4 -2
- nucleardatapy/fig/astro_setupMR_fig.py +14 -8
- nucleardatapy/fig/astro_setupMasses_fig.py +4 -2
- nucleardatapy/fig/astro_setupMtov_fig.py +6 -2
- nucleardatapy/fig/astro_setupMup_fig.py +4 -2
- nucleardatapy/fig/corr_setupEsymDen_fig.py +10 -4
- nucleardatapy/fig/corr_setupEsymLsym_fig.py +7 -3
- nucleardatapy/fig/corr_setupKsatQsat_fig.py +4 -2
- nucleardatapy/fig/crust_setupCrust_fig.py +6 -4
- nucleardatapy/fig/eos_setupAMBeq_fig.py +162 -134
- nucleardatapy/fig/eos_setupAMLeq_fig.py +35 -32
- nucleardatapy/fig/eos_setupAM_asy_lep_fig.py +39 -30
- nucleardatapy/fig/eos_setupAM_asy_nuc_fig.py +39 -30
- nucleardatapy/fig/eos_setupAM_asy_tot_fig.py +39 -30
- nucleardatapy/fig/eos_setupAM_fig.py +45 -30
- nucleardatapy/fig/eos_setupCC_fig.py +54 -24
- nucleardatapy/fig/hnuc_setupChart_fig.py +17 -14
- nucleardatapy/fig/hnuc_setupRE1LExp_fig.py +12 -2
- nucleardatapy/fig/matter_all_fig.py +55 -38
- nucleardatapy/fig/matter_setupCheck_fig.py +6 -3
- nucleardatapy/fig/matter_setupFFGLep_fig.py +13 -9
- nucleardatapy/fig/matter_setupFFGNuc_fig.py +22 -14
- nucleardatapy/fig/matter_setupHIC_fig.py +9 -9
- nucleardatapy/fig/matter_setupMicroEsym_fig.py +31 -30
- nucleardatapy/fig/matter_setupMicro_LP_fig.py +11 -5
- nucleardatapy/fig/matter_setupMicro_band_fig.py +11 -6
- nucleardatapy/fig/matter_setupMicro_effmass_fig.py +5 -4
- nucleardatapy/fig/matter_setupMicro_err_NM_fig.py +4 -2
- nucleardatapy/fig/matter_setupMicro_fig.py +30 -22
- nucleardatapy/fig/matter_setupMicro_gap_fig.py +22 -12
- nucleardatapy/fig/matter_setupNEPStats_fig.py +13 -3
- nucleardatapy/fig/matter_setupPhenoEsym_fig.py +28 -27
- nucleardatapy/fig/matter_setupPheno_fig.py +22 -17
- nucleardatapy/fig/nuc_setupBEExp_chart_fig.py +40 -31
- nucleardatapy/fig/nuc_setupBEExp_fig.py +49 -46
- nucleardatapy/fig/nuc_setupBETheo_fig.py +55 -50
- nucleardatapy/fig/nuc_setupISGMRExp_fig.py +7 -10
- nucleardatapy/fig/nuc_setupRchExp_fig.py +12 -20
- nucleardatapy/fig/nuc_setupRchTheo_fig.py +16 -20
- nucleardatapy/fig/nuc_setupRnpExp_fig.py +12 -0
- nucleardatapy/fig/nuc_setupRnpTheo_fig.py +14 -2
- nucleardatapy/hnuc/setup_re2L_exp.py +30 -7
- nucleardatapy/matter/setup_micro.py +424 -12
- nucleardatapy/matter/setup_micro_esym.py +46 -42
- nucleardatapy/matter/setup_nep.py +1 -1
- nucleardatapy/matter/setup_pheno.py +10 -6
- nucleardatapy/matter/setup_pheno_esym.py +14 -6
- nucleardatapy/nuc/setup_be_exp.py +4 -2
- nucleardatapy/nuc/setup_be_theo.py +19 -12
- {nucleardatapy-1.0.0.dist-info → nucleardatapy-1.0.1.dist-info}/METADATA +1 -1
- {nucleardatapy-1.0.0.dist-info → nucleardatapy-1.0.1.dist-info}/RECORD +87 -61
- {nucleardatapy-1.0.0.dist-info → nucleardatapy-1.0.1.dist-info}/WHEEL +0 -0
- {nucleardatapy-1.0.0.dist-info → nucleardatapy-1.0.1.dist-info}/licenses/LICENSE +0 -0
- {nucleardatapy-1.0.0.dist-info → nucleardatapy-1.0.1.dist-info}/top_level.txt +0 -0
nucleardatapy/astro/setup_mr.py
CHANGED
|
@@ -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
|
|
211
|
+
#print('ele[0]:',ele[0],' obs:',obs)
|
|
212
|
+
#print('ele[:]:',ele[:])
|
|
188
213
|
if int(ele[0]) == obs:
|
|
189
|
-
|
|
190
|
-
|
|
191
|
-
|
|
192
|
-
|
|
193
|
-
|
|
194
|
-
|
|
195
|
-
|
|
196
|
-
|
|
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
|
-
|
|
236
|
+
except ValueError:
|
|
201
237
|
self.comp = 0.0
|
|
202
|
-
self.comp_sig_up =
|
|
203
|
-
self.comp_sig_lo =
|
|
204
|
-
|
|
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
|
-
|
|
345
|
-
|
|
346
|
-
|
|
347
|
-
|
|
348
|
-
|
|
349
|
-
|
|
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
|
-
|
|
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
|
-
|
|
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
|
-
|
|
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
|
|
nucleardatapy/create_folder.py
CHANGED
|
@@ -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.
|
|
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.
|
|
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.
|
|
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.
|
|
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("
|
|
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,
|
|
2
|
-
# km, km, km, Mo, Mo, Mo, -, -, -, Hz
|
|
3
|
-
1, 13.02, 1.24, 1.06, 1.440, 0.150, 0.140, 0.
|
|
4
|
-
2, 12.71, 1.14, 1.19, 1.340, 0.150, 0.160, 0.
|
|
5
|
-
3, 14.44, 0.88, 1.05, 1.700, 0.180, 0.190, 0.
|
|
6
|
-
4, 11.71, 0.88, 0.83, 1.400, 0.130, 0.120, 0.1773, 0.0056, 0.0074, ,
|
|
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
|
|
@@ -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
|