We report on the core-collapse supernova simulation we conducted for a 11.2M progenitor model in three-dimensional space up to 20 ms after bounce, using a radiation hydrodynamics code with full Boltzmann neutrino transport. We solve the six-dimensional Boltzmann equations for three neutrino species and the three-dimensional compressible Euler equations with Furusawa and Togashi's nuclear equation of state. We focus on the prompt convection at ~10 ms after bounce and investigate how neutrinos are transported in the convective matter. We apply a new analysis based on the eigenvalues and eigenvectors of the Eddington tensor and make a comparison between the Boltzmann transport results and the M1 closure approximation in the transition regime between the optically thick and thin limits. We visualize the eigenvalues and eigenvectors using an ellipsoid, in which each principal axis is parallel to one of the eigenvectors and has a length proportional to the corresponding eigenvalue. This approach enables us to understand the difference between the Eddington tensor derived directly from the Boltzmann simulation and the one given by the M1 prescription from a new perspective. We find that the longest principal axis of the ellipsoid is almost always nearly parallel to the energy flux in the M1 closure approximation whereas in the Boltzmann simulation it becomes perpendicular in some transition regions, where the mean free path is ~0.1 times the radius. In three spatial dimensions, the convective motions make it difficult to predict where this happens and possibly improves the closure relation there.

Using the Boltzmann-radiation-hydrodynamics code which solves the Boltzmann equation for the neutrino transport, we present the results of the simulations with the nuclear equations of state (EOSs) of Lattimer-Swesty (LS) and Furusawa-Shen (FS). We extend the simulation time of the LS model and conducted thorough investigations, though the previous paper reported some of the results briefly. Only the LS model shows the shock revival. This seems to originate from the nuclear composition: the different nuclear composition results in the different energy loss by the photodissociation and hence the different strength of the prompt convection and the later neutrino-driven convection. The proto-neutron star seen in the FS model is more compact than that in the LS model since the existence of multi-nuclear species softens the EOS. For the behavior of neutrinos, we examined the flux and the Eddington tensor of neutrinos. In the optically thick region, the diffusion of neutrinos and the dragging by the matter motion determine the flux. In the optically thin region, the free-streaming determines it. The Eddington tensor is compared with that obtained from the M1-closure relation. The M1-closure scheme overestimates the contribution from the velocity-dependent terms in the semi-transparent region.

Neutrinos are believed to have a key role in the explosion mechanism of core-collapse supernovae as they carry most of the energy released by the gravitational collapse of a massive star. If their flavor is converted fast inside the neutrino sphere, the supernova explosion may be influenced. This paper is reporting the results of the extended work of our previous paper. We perform a thorough survey of the ELN crossing in one of our self-consistent, realistic Boltzmann simulations in two spatial dimensions under axisymmetry for the existence of the crossings between $\nu_e$ and $\bar\nu_e$ angular distributions, or the electron lepton number (ELN) crossing. We report for the first time the positive detections deep inside the core of the massive star in the vicinity of neutrino sphere at $r$ $\approx$ 16 - 21 km. We find that low values of the electron fraction $Y_e$ produced by convective motions together with the appearance of light elements are critically important to give rise to the ELN crossing by enhancing the chemical potential difference between proton and neutron, and hence by mitigating the Fermi-degeneracy of $\nu_e$. Since the region of positive detection are sustained and, in fact, expanding with time, it may have an impact on the explosion of core-collapse supernovae, observational neutrino astronomy and nucleosynthesis of heavy nuclei.

We present a much improved equation of state for neutron star matter, QHC19, with a smooth crossover from the hadronic regime at lower densities to the quark regime at higher densities. We now use the Togashi et al.~equation of state (Togashi:2017), a generalization of the Akmal-Pandharipande-Ravenhall equation of state of uniform nuclear matter, in the entire hadronic regime; the Togashi equation of state consistently describes non-uniform as well as uniform matter, and matter at beta equilibrium without the need for an interpolation between pure neutron and symmetric nuclear matter. We describe the quark matter regime at higher densities with the Nambu--Jona--Lasinio model, now identifying tight constraints on the phenomenological universal vector repulsion between quarks and the pairing interaction between quarks arising from the requirements of thermodynamic stability and causal propagation of sound. The resultant neutron star properties agree very well with the inferences of the LIGO/Virgo collaboration, from GW170817, of the pressure vs. baryon density, neutron star radii, and tidal deformabilities. The maximum neutron star mass allowed by QHC19 is 2.35 $M_\odot$, consistent with all neutron star mass determinations.

With the Boltzmann-radiation-hydrodynamics code, which we have developed to solve numerically the Boltzmann equations for neutrino transfer, the Newtonian hydrodynamics equations, and the Newtonian self-gravity simultaneously and consistently, we simulate the collapse of a rotating core of the progenitor with a zero-age-main-sequence mass of $11.2\,M_\odot$ and a shelluler rotation of $1\,{\rm rad\,s^{-1}}$ at the center. We pay particular attention in this paper to the neutrino distribution in phase space, which is affected by the rotation. By solving the Boltzmann equations directly, we can assess the rotation-induced distortion of the angular distribution in momentum space, which gives rise to the rotational component of the neutrino flux. We compare the Eddington tensors calculated both from the raw data and from the M1-closure approximation. We demonstrate that the Eddington tensor is determined by complicated interplays of the fluid velocity and the neutrino interactions and that the M1-closure, which assumes that the Eddington factor is determined by the flux factor, fails to fully capture this aspect, especially in the vicinity of the shock. We find that the error in the Eddington factor reaches $\sim 20\%$ in our simulation. This is due not to the resolution but to the different dependence of the Eddington and flux factors on the angular profile of the neutrino distribution function, and hence modification to the closure relation is needed.

We perform an extensive study of the influence of nuclear weak interactions on core-collapse supernovae (CCSNe), paying particular attention to consistency between nuclear abundances in the equation of state (EOS) and nuclear weak interactions. We compute properties of uniform matter based on the variational method. For inhomogeneous nuclear matter, we take a full ensemble of nuclei into account with various finite-density and thermal effects and directly use the nuclear abundances to compute nuclear weak interaction rates. To quantify the impact of a consistent treatment of nuclear abundances on CCSN dynamics, we carry out spherically symmetric CCSN simulations with full Boltzmann neutrino transport, systematically changing the treatment of weak interactions, EOSs, and progenitor models. We find that the inconsistent treatment of nuclear abundances between the EOS and weak interaction rates weakens the EOS dependence of both the dynamics and neutrino signals. We also test the validity of two artificial prescriptions for weak interactions of light nuclei and find that both prescriptions affect the dynamics. Furthermore, there are differences in neutrino luminosities by ~10% and in average neutrino energies by 0.25-1 MeV from those of the fiducial model. We also find that the neutronization burst neutrino signal depends on the progenitor more strongly than on the EOS, preventing a detection of this signal from constraining the EOS.

I have systematically investigated the equations of state (EOSs) in nuclear statistical equilibrium under thermodynamic conditions relevant for core collapse of massive stars by varying the bulk properties of nuclear matter, the mass data for neutron-rich nuclei, and the finite-temperature modifications of the nuclear model. It is found that the temperature dependence of the nuclear free energies has a significant impact on the entropy and nuclear composition, which affect the dynamics of core-collapse supernovae. There is a little influence from the bulk properties and the mass data. For all models, common nuclei that are likely to contribute to core-deleptonization are those near $Z\approx30$ and $N\approx50$. A model with a semi-empirical expression for internal degrees of freedom, however, overestimates the number densities of magic nuclei with $N\approx50$ and $82$, while a model, in which nuclear shell effects are not considered, underestimates the number densities of heavy nuclei, and especially of the magic nuclei. Other models, which include the temperature dependence of shell effects in the internal degrees of freedom and/or in the nuclear internal free energy, indicate that the difference in population between magic nuclei and non-magic nuclei disappears as the temperature increases. The construction of complete statistical EOS will require further theoretical and experimental studies of medium-mass, neutron-rich nuclei with proton numbers 25-45 and neutron numbers 40-85.

We investigate the high-temperature effect on the nuclear matter that consists of mixture of nucleons and all nuclei in the dense and hot stellar environment. The individual nuclei are described within the compressible liquid-drop model that is based on Skyrme interactions for bulk energies and that takes into account modifications of the surface and Coulomb energies at finite temperatures and densities. The free-energy density is minimized with respect to the individual equilibrium densities of all heavy nuclei and the nuclear composition. We find that their optimized equilibrium densities become smaller and smaller at high temperatures because of the increase of thermal contributions to bulk free energies and the reduction of surface energies. The neutron-rich nuclei become unstable and disappear one after another at some temperatures. The calculations are performed for two sets of model parameters leading to different values of the slope parameter in the nuclear symmetry energy. It is found that the larger slope parameter reduces the equilibrium densities and the melting temperatures. We also compare the new model with some other approaches and find that the mass fractions of heavy nuclei in the previous calculations that omit vaporization are underestimated at $T \lesssim 10$MeV and overestimated at $T\gtrsim10$~MeV. The further sophistication of calculations of nuclear vaporization and of light clusters would be required to construct the equation of state for explosive astrophysical phenomena.

We present the first results of our spatially axisymmetric core-collapse supernova simulations with full Boltzmann neutrino transport, which amount to a time-dependent 5-dimensional (2 in space and 3 in momentum space) problem in fact. Special relativistic effects are fully taken into account with a two-energy-grid technique. We performed two simulations for a progenitor of 11.2M, employing different nuclear equations-of-state (EOS's): Lattimer and Swesty's EOS with the incompressibility of K = 220MeV (LS EOS) and Furusawa's EOS based on the relativistic mean field theory with the TM1 parameter set (FS EOS). In the LS EOS the shock wave reaches ~700km at 300ms after bounce and is still expanding whereas in the FS EOS it stalled at ~200km and has started to recede by the same time. This seems to be due to more vigorous turbulent motions in the former during the entire post-bounce phase, which leads to higher neutrino-heating efficiency in the neutrino-driven convection. We also look into the neutrino distributions in momentum space, which is the advantage of the Boltzmann transport over other approximate methods. We find non-axisymmetric angular distributions with respect to the local radial direction, which also generate off-diagonal components of the Eddington tensor. We find that the r {\theta}-component reaches ~10% of the dominant rr-component and, more importantly, it dictates the evolution of lateral neutrino fluxes, dominating over the {\theta}{\theta}-component, in the semi-transparent region. These data will be useful to further test and possibly improve the prescriptions used in the approximate methods.

This paper is a sequel to our previous one (Kato et al.2015), which calculated the luminosities and spectra of electron-type anti-neutrinos ($\bar{\nu}_e$'s) from the progenitors of core-collapse supernovae. Expecting that a capability to detect electron-type neutrinos ($\nu_e$'s) will increase dramatically with the emergence of liquid-argon detectors such as DUNE, we broaden the scope in this study to include all-flavors of neutrinos emitted from the pre-bounce phase. We pick up three progenitor models of an electron capture supernova (ECSN) and iron-core collapse supernovae (FeCCSNe). We find that the number luminosities reach $\sim10^{57} \mathrm{s^{-1}}$ and $\sim10^{53} \mathrm{s^{-1}}$ at maximum for $\nu_e$ and $\bar{\nu}_e$, respectively. We also estimate the numbers of detection events at terrestrial neutrino detectors including DUNE, taking flavor oscillations into account and assuming the distance to the progenitors to be 200 pc. It is demonstrated that $\bar{\nu}_e$'s from the ECSN-progenitor will be undetected at almost all detectors, whereas we will be able to observe $\gtrsim$15900 $\nu_e$'s at DUNE for the inverted mass hierarchy. From the FeCCSN-progenitors, the number of $\bar{\nu}_e$ events will be largest for JUNO, 200-900 $\bar{\nu}_e$'s, depending on the mass hierarchy whereas the number of $\nu_e$ events at DUNE is $\gtrsim$2100 for the inverted mass hierarchy. These results imply that the detection of $\bar{\nu}_e$'s is useful to distinguish FeCCSN- from ECSN-progenitors, while $\nu_e$'s will provide us with detailed information on the collapse phase regardless of the type and mass of progenitor.

We construct new equations of state for baryons at sub-nuclear densities for the use in core-collapse supernova simulations. The abundance of various nuclei is obtained together with thermodynamic quantities. The formulation is an extension of the previous model, in which we adopted the relativistic mean field theory with the TM1 parameter set for nucleons, the quantum approach for $d$, $t$, $h$ and $\alpha$ as well as the liquid drop model for the other nuclei under the nuclear statistical equilibrium. We reformulate the model of the light nuclei other than $d$, $t$, $h$ and $\alpha$ based on the quasi-particle description. Furthermore, we modify the model so that the temperature dependences of surface and shell energies of heavy nuclei could be taken into account. The pasta phases for heavy nuclei and the Pauli- and self-energy shifts for $d$, $t$, $h$ and $\alpha$ are taken into account in the same way as in the previous model. We find that nuclear composition is considerably affected by the modifications in this work, whereas thermodynamical quantities are not changed much. In particular, the washout of shell effect has a great impact on the mass distribution above $T \sim 1$ MeV. This improvement may have an important effect on the rates of electron captures and coherent neutrino scatterings on nuclei in supernova cores.

We have constructed a nuclear equation of state (EOS) that includes a full nuclear ensemble for use in core-collapse supernova simulations. It is based on the EOS for uniform nuclear matter that two of the authors derived recently, applying a variational method to realistic two- and there-body nuclear forces. We utilize an extended liquid drop model of heavy nuclei and a quantum-theoretical mass evaluation for light nuclei. In addition to realistic nuclear forces, the inclusion of in-medium effects on the full ensemble of nuclei makes the new EOS one of the most realistic EOS's for supernova simulations. We make comparisons with the FYSS EOS, which is based on the same formulation for the nuclear ensemble but adopts the relativistic mean field (RMF) theory with the TM1 parameter set for uniform nuclear matter. The new EOS is softer than the FYSS EOS around and above nuclear saturation densities. We find that neutron-rich nuclei with small mass numbers are more abundant in the new EOS than in the FYSS EOS because of the larger saturation densities and smaller symmetry energy of nuclei in the former. We apply the two EOS's to 1D supernova simulations and find that the new EOS gives lower electron fractions and higher temperatures in the collapse phase owing to the smaller symmetry energy. As a result, the inner core has smaller masses for the new EOS. It is more compact, on the other hand, due to the softness of the new EOS and bounces at higher densities. The ensuing outward propagations of the shock wave in the outer core are very similar in the two simulations, which may be an artifact, though, caused by the use of the same tabulated electron capture rates for heavy nuclei ignoring differences in the nuclear composition between the two EOS's in these computations.

We investigate the influences of the nuclear composition on the weak interaction rates of heavy nuclei during the core collapse of massive stars. The nuclear abundances in nuclear statistical equilibrium (NSE) are calculated by some equation of state (EOS) models including in-medium effects on nuclear masses. We systematically examine the sensitivities of electron capture and neutrino-nucleus scattering on heavy nuclei to the nuclear shell effects and the single nucleus approximation. We find that the washout of shell effects at high temperatures brings significant change to weak rates by smoothing the nuclear abundance distribution: the electron capture rate decreases by $\sim$20$\%$ in the early phase and increases by $\sim$40$\%$ in the late phase at most, while the cross section for neutrino-nucleus scattering is reduced by $\sim$15$\%$. This is because the open-shell nuclei become abundant instead of those with closed neutron shells as the shell effects disappear. We also find that the single-nucleus description based on the average values leads to underestimations of weak rates. Electron captures and neutrino coherent scattering on heavy nuclei are reduced by $\sim$80$\%$ in the early phase and by $\sim$5$\%$ in the late phase, respectively. These results indicate that NSE like EOS accounting for shell washout is indispensable for the reliable estimation of weak interaction rates in simulations of core-collapse supernovae.

We investigate the mass fractions and in-medium properties of heavy nuclei in stellar matter at characteristic densities and temperatures for supernova (SN) explosions. The individual nuclei are described within the compressible liquid-drop model taking into account modifications of bulk, surface and Coulomb energies. The equilibrium properties of nuclei and the full ensemble of heavy nuclei are calculated self-consistently. It is found that heavy nuclei in the ensemble are either compressed or decompressed depending on the isospin asymmetry of the system. The compression or decompression has a little influence on the binding energies, total mass fractions and average mass numbers of heavy nuclei, although the equilibrium densities of individual nuclei themselves are changed appreciably above one hundredth of normal nuclear density. We find that nuclear structure in single nucleus approximation deviates from the actual one obtained in the multi-nucleus description, since the density of free nucleons is different between these two descriptions. This study indicates that a multi-nucleus description is required to realistically account for in-medium effects on the nuclear structure in supernova matter.

We present a newly developed moving-mesh technique for the multi-dimensional Boltzmann-Hydro code for the simulation of core-collapse supernovae (CCSNe). What makes this technique different from others is the fact that it treats not only hydrodynamics but also neutrino transfer in the language of the 3+1 formalism of general relativity (GR), making use of the shift vector to specify the time evolution of the coordinate system. This means that the transport part of our code is essentially general relativistic although in this paper it is applied only to the moving curvilinear coordinates in the flat Minknowski spacetime, since the gravity part is still Newtonian. The numerical aspect of the implementation is also described in detail. Employing the axisymmetric two-dimensional version of the code, we conduct two test computations: oscillations and runaways of proto-neutron star (PNS). We show that our new method works fine, tracking the motions of PNS correctly. We believe that this is a major advancement toward the realistic simulation of CCSNe.

We study transitions of hadronic matter (HM) to 3-flavor quark matter (3QM), regarding the conversion processes as combustion and describing them hydrodynamically. Under the assumption that HM is metastable with their free energies being larger than those of 3QM but smaller than those of 2-flavor quark matter (2QM), we consider in this paper the conversion induced by diffusions of seed 3QM. This is a sequel to our previous paper, in which the shock-induced conversion was studied in the same frame work. We not only pay attention to the jump condition on both sides of the conversion front but the structures inside the front are also considered by taking into account what happens during the conversion processes on the time scale of weak interactions. We employ for HM the Shen's EOS, which is based on the relativistic mean field theory, and the bag model-based EOS for QM just as in the previous paper. We demonstrated in that paper that in this combination of EOS's the combustion will occur for a wide range of the bag constant and strong coupling constant in the so-called endothermic regime, in which the Hugoniot curve for combustion runs below the initial state. We find that weak deflagration nearly always occurs and that weak detonation is possible only when the diffusion constant is (unrealistically) large and the critical strange fraction is small. The velocities of the conversion front are ~ $10^3-10^7$ cm/s depending on the initial temperature and density as well as the parameters in the QM EOS and become particularly small when the final state is in the mixed phase. Finally we study linear stability of the laminar weak-deflagration front and find that it is unstable in the exothermic regime (Darrius-Landau instability) but stable in the endothermic regime, which is quite contrary to the ordinary combustions.

We study transitions of hadronic matter (HM) to 3-flavor quark matter (3QM) locally, regarding the conversion processes as combustion and describing them hydrodynamically. Not only the jump condition on both sides of the conversion front but the structures inside the front are also considered by taking into account what happens during the conversion processes on the time scale of weak interactions as well as equations of state (EOS's) in the mixed phase. Under the assumption that HM is metastable with their free energies being larger than those of 3QM but smaller than those of 2-flavor quark matter (2QM), we consider the transition via 2QM triggered by a rapid density rise in a shock wave. Based on the results, we discuss which combustion modes (strong/weak detonation) may be realized. HM is described by an EOS based on the relativistic mean field theory and 2, 3QM's are approximated by the MIT bag model. We demonstrate for a wide range of bag constant and strong coupling constant in this combination of EOS's that the combustion may occur in the so-called endothermic regime, in which the Hugoniot curve for combustion runs below the one for the shock wave in P-V plane, and which has no terrestrial counter part. We find that strong detonation always occurs. Depending on the EOS of quark matter (QM) as well as the density of HM and the Mach number of the detonation front, deconfinement from HM to 2QM is either completed or not completed in the shock wave. In the latter case, which is more likely if the EOS of QM ensures that deconfinement occurs above the nuclear saturation density and that the maximum mass of cold quark stars is larger than two solar mass, the conversion continues further via the mixing state of HM and 3QM on the time scale of weak interactions.

We perform numerical experiments to investigate the influence of inelastic neutrino reactions with light nuclei on the standing accretion shock instability (SASI). The time evolution of shock waves is calculated with a simple light-bulb approximation for the neutrino transport and a multi-nuclei equation of state. The neutrino absorptions and inelastic interactions with deuterons, tritons, helions and alpha particles are taken into account in the hydrodynamical simulations in addition to the ordinary charged-current interactions with nucleons. Axial symmetry is assumed but no equatorial symmetry is imposed. We show that the heating rates of deuterons reach as high as 10% of those of nucleons around the bottom of the gain region. On the other hand, alpha particles are heated near the shock wave, which is important when the shock wave expands and the density and temperature of matter become low. It is also found that the models with heating by light nuclei have different evolutions from those without it in the non-linear phase of SASI. This results is because matter in the gain region has a varying density and temperature and there appear sub-regions that are locally rich in deuterons and alpha particles. Although the light nuclei are never dominant heating sources and they work favorably for shock revival in some cases and unfavorably in other cases, they are non-negligible and warrant further investigation.

We construct new equations of state for baryons at sub-nuclear densities for the use in core-collapse simulations of massive stars. The abundance of various nuclei is obtained together with thermodynamic quantities. A model free energy is constructed, based on the relativistic mean field theory for nucleons and the mass formula for nuclei with the proton number up to ~ 1000. The formulation is an extension of the previous model, in which we adopted the liquid drop model to all nuclei under the nuclear statistical equilibrium. We reformulate the new liquid drop model so that the temperature dependences of bulk energies could be taken into account. Furthermore, we extend the region in the nuclear chart, in which shell affects are included, by using theoretical mass data in addition to experimental ones. We also adopt a quantum theoretical mass evaluation of light nuclei, which incorporates the Pauli- and self-energy shifts that are not included in the ordinary liquid drop model. The pasta phases for heavy nuclei are taken into account in the same way as in the previous model. We find that the abundances of heavy nuclei are modified by the shell effects of nuclei and temperature dependence of bulk energies. These changes may have an important effect on the rates of electron captures and coherent neutrino scatterings on nuclei in supernova cores. The abundances of light nuclei are also modified by the new mass evaluation, which may affect the heating and cooling rates of supernova cores and shocked envelopes.

We compare three different statistical models for the equation of state (EOS) of stellar matter at subnuclear densities and temperatures (0.5-10 MeV) expected to occur during the collapse of massive stars and supernova explosions. The models introduce the distributions of various nuclear species in nuclear statistical equilibrium, but use somewhat different nuclear physics inputs. It is demonstrated that the basic thermodynamical quantities of stellar matter under these conditions are similar, except in the region of high densities and low temperatures. We demonstrate that mass and isotopic distributions have considerable differences related to the different assumptions of the models on properties of nuclei at these stellar conditions. Overall, the three models give similar trends, but the details reflect the uncertainties related to the modeling of medium effects, such as the temperature and density dependence of surface and bulk energies of heavy nuclei, and the nuclear shell structure effects. We discuss importance of new physics inputs for astrophysical calculations from experimental data obtained in intermediate energy heavy-ion collisions, in particular, the similarities of the conditions reached during supernova explosions and multifragmentation reactions.

We investigate the emergence of hyperons in black-hole-forming failed supernovae, which are caused by the dynamical collapse of nonrotating massive stars. We perform neutrino-radiation hydrodynamical simulations in general relativity adopting realistic hyperonic equation-of-state (EOS). Attractive and repulsive cases are examined for the potential of \Sigma\ hyperons. Since hyperons soften the EOS, they shorten the time interval from the bounce to black hole formation, which corresponds to the duration of neutrino emission. This effect is larger for the attractive case than the repulsive case because \Sigma\ hyperons appear more easily. In addition, we investigate the impacts of pions to find that they also promotes the recollapse towards the black hole formation.

We calculate a new equation of state for baryons at sub-nuclear densities meant for the use in core-collapse simulations of massive stars. The abundance of various nuclei is obtained together with the thermodynamic quantities. The formulation is the NSE description and the liquid drop approximation of nuclei. The model free energy to minimize is calculated by relativistic mean field theory for nucleons and the mass formula for nuclei with the atomic number up to ~ 1000. We have also taken into account the pasta phase, thanks to which the transition to uniform nuclear matter in our EOS occurs in the conventional manner: nuclei are not dissociated to nucleons but survive right up to the transition to uniform nuclear matter. We find that the free energy and other thermodynamical quantities are not very different from those given in the Shen's EOS, one of the standard EOS's that adopt the single nucleus approximation. The average mass is systematically different, on the other hand, which may have an important ramification to the rates of electron captures and coherent neutrino scatterings on nuclei in supernova cores. It is also interesting that the root mean square of the mass number is not very different from the average mass number, since the former is important for the evaluation of coherent scattering rates on nuclei but has been unavailable so far. The EOS table is currently under construction, which will include the weak interaction rates.