Commit d6146555 authored by Chu, Ran's avatar Chu, Ran
Browse files

Updated README.md: adding more info.

parent a376ff4c
# WeakLib Tables
- WeakLib source code at GitHub:<a href="https://github.com/starkiller-astro/weaklib"> https://github.com/starkiller-astro/weaklib</a>
- **WeakLib source code**:<a href="https://github.com/starkiller-astro/weaklib"> https://github.com/starkiller-astro/weaklib</a>
- The provenance of these tables are as follows:
- The provenance of these tables:
( Origin/Type of Table ) - ( Table Identifier ) - ( Resolution of Table in rho-T-Ye ) - (Opacity Type If Applied)
......@@ -11,7 +11,7 @@
wl-EOS : This is a Weaklib generated EOS Table.
wl-Op : This is a Weaklib generated Opacity Table.
- Table Identifier:
- Model Identifier:
LS220, SFHo, SFHx
......@@ -20,88 +20,145 @@
EOS example: 25-40-100 is 20 pts/decade in density, 40 pts/decade in temperature, and 100 pts/1.00 in electron fraction
Op example: 25-40-100-E40 is 20 pts/decade in density, 40 pts/decade in temperature, 100 pts/1.00 in electron fraction, 40 pts from 0.1 MeV - 300 MeV in energy
- Opacity Type and Tabulated Rates
- Opacity Type:
- B85 : based on <a href="http://adsabs.harvard.edu/pdf/1985ApJS...58..771B">Bruenn (1985) DOI: 10.1086/191056</a>
- B85 : based on Bruenn (1985)
- AbEm : neutrino absorptivity 1/&lambda;(e) and emissivity j(e) on nucleons and nuclei.
- AbEm : neutrino absorptivity 1/&lambda;(e) and emissivity j(e) on nucleons and nuclei
Tabulated(e,rho,T,ye) = 1/&lambda;(e,rho,T,ye) + j(e,rho,T,ye)
Tabulated(e,rho,T,ye) = 1/λ(e,rho,T,ye) + j(e,rho,T,ye)
for electron-type neutrino and electron-type antineutrino
in unit: [cm-1]
- Iso : neutrino isoenergetic scattering on nucleons and nuclei.
- Iso : neutrino isoenergetic scattering on nucleons and nuclei
Tabulated(e,1,rho,T,ye) = 1/2 * 4&pi;/(cvel*(hc)**3) * e2 * Phi_0
Tabulated(e,2,rho,T,ye) = 3/2 * 4&pi;/(cvel*(hc)**3) * e2 * Phi_1
Tabulated(e,1,rho,T,ye) = 1/2 * 4π/(cvel*(hc)**3) * e2 * Phi_0
Tabulated(e,2,rho,T,ye) = 3/2 * 4π/(cvel*(hc)**3) * e2 * Phi_1
for neutrino and antineutrino
in unit: [cm-1]
- NES : neutrino-electron scattering.
Tabulated(ep,e,1,T,eta) = H0i * 2&pi;/(cvel*(hc)**3) = zero moment of the "i" NES outgoing scattering kernal
Tabulated(ep,e,2,T,eta) = H0ii * 2&pi;/(cvel*(hc)**3) = zero moment of the "ii" NES outgoing scattering kernal
Tabulated(ep,e,3,T,eta) = H1i * 2&pi;/(cvel*(hc)**3) = first moment of the "i" NES outgoing scattering kernal
Tabulated(ep,e,4,T,eta) = H1ii * 2&pi;/(cvel*(hc)**3) = first moment of the "ii" NES outgoing scattering kernal
- Pair : neutrino production from pair process.
- NES : neutrino-electron scattering
Tabulated(ep,e,1,T,eta) = J0i * 2&pi;/(cvel*(hc)**3) = zero moment of the "i" pair annihilation kernal
Tabulated(ep,e,2,T,eta) = J0ii * 2&pi;/(cvel*(hc)**3) = zero moment of the "ii" pair annihilation kernal
Tabulated(ep,e,3,T,eta) = J1i * 2&pi;/(cvel*(hc)**3) = first moment of the "i" pair annihilation kernal
Tabulated(ep,e,4,T,eta) = J1ii * 2&pi;/(cvel*(hc)**3) = first moment of the "ii" pair annihilation kernal
Tabulated(ep,e,1,T,eta) = H0i * 2π/(cvel*(hc)**3) = zero moment of the "i" NES outgoing scattering kernal
Tabulated(ep,e,2,T,eta) = H0ii * 2π/(cvel*(hc)**3) = zero moment of the "ii" NES outgoing scattering kernal
Tabulated(ep,e,3,T,eta) = H1i * 2π/(cvel*(hc)**3) = first moment of the "i" NES outgoing scattering kernal
Tabulated(ep,e,4,T,eta) = H1ii * 2π/(cvel*(hc)**3) = first moment of the "ii" NES outgoing scattering kernal
in unit: [MeV-3 cm-1]
- Ref: <a href="https://ui.adsabs.harvard.edu/abs/1985ApJS...58..771B/abstract"> Bruenn(1985), 10.1086/191056</a>
- Pair : neutrino production from pair process
Tabulated(ep,e,1,T,eta) = J0i * 2π/(cvel*(hc)**3) = zero moment of the "i" pair annihilation kernal
Tabulated(ep,e,2,T,eta) = J0ii * 2π/(cvel*(hc)**3) = zero moment of the "ii" pair annihilation kernal
Tabulated(ep,e,3,T,eta) = J1i * 2π/(cvel*(hc)**3) = first moment of the "i" pair annihilation kernal
Tabulated(ep,e,4,T,eta) = J1ii * 2π/(cvel*(hc)**3) = first moment of the "ii" pair annihilation kernal
in unit: [MeV-3 cm-1]
## Implement Example
- For interpolating with a given rho-t-ye profile, see source code:
weaklib/Distributions/UnitTests/Opacities
- <a href="https://github.com/starkiller-astro/weaklib/tree/master/Distributions/UnitTests/Opacities">weaklib/Distributions/UnitTests/Opacities</a>
- <a href="https://github.com/starkiller-astro/weaklib/tree/master/Distributions/Tools">weaklib/Distributions/Tools</a>
- For using tables in a real simulation
- <a href="https://github.com/endeve/thornado/tree/master/Modules/EquationOfState">Thornado/Modules/EquationOfState</a>
- <a href="https://github.com/endeve/thornado/tree/master/Modules/Opacities">Thornado/Modules/Opacities</a>
- FLASH5( ongoing )
## Avaliable Table:
#### SFHo
HighRes:
Size | FileName
------|-------------------------------
352M | wl-EOS-SFHo-25-50-100-Standard.h5
411M | wl-EOS-SFHo-25-50-100.h5 ( extended ye to 0.7 )
Size | FileName | Cover Range | IDTag
------ |----------------------------------------|------------------------------|------------------------
| | (1.661e4, 3.164e15) g cm-3 | wl-EOS-SFHo-25-50-100
352M | wl-EOS-SFHo-25-50-100-Standard.h5 | (1.161e9, 1.839e12) K | 7-24-17, no high Ye bck
| | (0.010, 0.600) |
------ |----------------------------------------|------------------------------|------------------------
| | (1.661e4, 3.164e15) g cm-3 | wl-EOS-SFHo-25-50-100
411M | wl-EOS-SFHo-25-50-100.h5 | (1.161e9, 1.839e12) K | 11-13-17, with clusters,
| | (0.010, 0.700) | renorm
LowRes:
Size | FileName
------|-------------------------------
53M | wl-EOS-SFHo-15-25-50.h5
274M | wl-Op-SFHo-15-25-50-E40-B85-AbEm.h5
549M | wl-Op-SFHo-15-25-50-E40-B85-Iso.h5
237M | wl-Op-SFHo-15-25-50-E40-B85-NES.h5
237M | wl-Op-SFHo-15-25-50-E40-B85-Pair.h5
Size | FileName | Cover Range | IDTag
------ |----------------------------------------|------------------------------|------------------------
| | (1.661e4, 3.164e15) g cm-3 | wl-EOS-SFHo-25-50-100
53M | wl-EOS-SFHo-15-25-50.h5 | (1.161e9, 1.839e12) K | 7-24-17, no high Ye bck
| | (0.010, 0.600) |
------ |----------------------------------------|------------------------------|------------------------
| | same as above, and |
274M | wl-Op-SFHo-15-25-50-E40-B85-AbEm.h5 | (1.000e-1, 3.000e2) MeV | NA
| | nue, nuebar |
------ |----------------------------------------|------------------------------|------------------------
| | same as above, and |
549M | wl-Op-SFHo-15-25-50-E40-B85-Iso.h5 | (1.000e-1, 3.000e2) MeV | NA
| | nue, nuebar |
| | (same for mu, tau) |
------ |----------------------------------------|------------------------------|------------------------
| | same as above, and |
237M | wl-Op-SFHo-15-25-50-E40-B85-NES.h5 | (1.000e-1, 3.000e2) MeV | NA
| | (1.000e-3, 2.500e3) in eta |
------ |----------------------------------------|------------------------------|------------------------
| | same as above, and |
237M | wl-Op-SFHo-15-25-50-E40-B85-Pair.h5 | (1.000e-1, 3.000e2) MeV | NA
| | (1.000e-3, 2.500e3) in eta |
#### LS220
HighRes:
Size | FileName
------|-------------------------------
246M | wl-EOS-LS220-25-50-100.h5 ( extended ye to 1.0 )
195M | wl-EOS-LSwBCK-20-40-100.h5 ( extended ye to 1.0 )
101M | wl-EOS-LS220-20-40-100-Lower-T.h5
87M | wl-EOS-LS220-20-40-100.h5
Size | FileName | Cover Range | IDTag
------ |----------------------------------------|------------------------------|---------------------------
| | (1.000e7, 1.000e15) g cm-3 | wl-EOS-LS220BCK-25-50-100,
259M | wl-EOS-LS220-25-50-100.h5 | (1.996e9, 1.000e12) K | 5-25-17
| | (0.060, 1.000) |
------ |----------------------------------------|------------------------------|---------------------------
| | (1.000e7, 1.000e15) g cm-3 | wl-EOS-LS220BCK-20-40-100,
205M | wl-EOS-LSwBCK-20-40-100.h5 | (1.996e9, 1.000e12) K | 4-17-17
| | (0.060, 1.000) |
------ |----------------------------------------|------------------------------|---------------------------
| | (1.000e7, 1.000e15) g cm-3 | wl-EOS-LS220-20-40-100,
106M | wl-EOS-LS220-20-40-100.h5 | (1.996e9, 1.000e12) K | 2-20-17
| | (0.060, 0.540) |
LowRes:
Size | FileName
------|-------------------------------
24M | wl-EOS-LS220-15-25-50-Lower-T.h5
126M | wl-Op-LS220-15-25-50-Lower-T-E40-B85-AbEm.h5
251M | wl-Op-LS220-15-25-50-Lower-T-E40-B85-Iso.h5
199M | wl-Op-LS220-15-25-50-Lower-T-E40-B85-NES.h5
199M | wl-Op-LS220-15-25-50-Lower-T-E40-B85-Pair.h5
Size | FileName | Cover Range | IDTag
------ |------------------------------------------------|------------------------------|------------------------
| | (1.000e7, 1.000e15) g cm-3 |
24M | wl-EOS-LS220-15-25-50-Lower-T.h5 | (1.996e9, 1.000e12) K | NA
| | (0.060, 0.540) |
------ |------------------------------------------------|------------------------------|------------------------
| | same as above, and |
126M | wl-Op-LS220-15-25-50-Lower-T-E40-B85-AbEm.h5 | (1.000e-1, 3.000e2) MeV | NA
| | nue, nuebar |
------ |------------------------------------------------|------------------------------|------------------------
| | same as above, and |
251M | wl-Op-LS220-15-25-50-Lower-T-E40-B85-Iso.h5 | (1.000e-1, 3.000e2) MeV | NA
| | nue, nuebar |
| | (same for mu, tau) |
------ |------------------------------------------------|------------------------------|------------------------
| | same as above, and |
199M | wl-Op-LS220-15-25-50-Lower-T-E40-B85-NES.h5 | (1.000e-1, 3.000e2) MeV | NA
| | (1.000e-3, 2.500e3) in eta |
------ |------------------------------------------------|------------------------------|------------------------
| | same as above, and |
199M | wl-Op-LS220-15-25-50-Lower-T-E40-B85-Pair.h5 | (1.000e-1, 3.000e2) MeV | NA
| | (1.000e-3, 2.500e3) in eta |
#### SFHx
HighRes:
Size | FileName
------|-------------------------------
435M | wl-EOS-SFHx-25-50-100.h5 ( extended ye to 0.7 )
Size | FileName | Cover Range | IDTag
------ |----------------------------------------|------------------------------|------------------------
| | (1.661e4, 3.164e15) g cm-3 | wl-EOS-NL3-25-50-100
435M | wl-EOS-SFHx-25-50-100.h5 | (1.161e9, 1.839e12) K | 12-04-17, with clusters
| | (0.010, 0.700) | renorm, 0.7
## Ask For Help
- R. Chu : rchu@vols.utk.edu
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment