Burn_nseAtDens( real(OUT) :: qbar_nse, real(OUT) :: sumyi_nse, real(OUT) :: approxtemp, real(OUT) :: edot, real(OUT) :: Yedot, real(IN) :: Ye, real(INOUT) :: dens, real(INOUT) :: emq)
This function gives the NSE final state for a constant density burn via table lookup. qbar is returned in MeV/nucleon. Only the compositional information (qbar, 1/Abar, edot, Yedot) should be used for hydrodynamics purposes. Accurate thermodynamic prperties of the NSE final state (rho, T) should be obtained by solving eint - qbar_nse = emq . This is important due to the limited accuracy of the table interpolation (i.e., the interpolated values cannot satisfy any constraints on their own.)
qbar_nse -- Average binding energy per nucleon (averaged over nse) sumyi_nse -- Sum over the abundances. (Y = X/A = X/(#of nucleons)) approxtemp -- ?? Ask Dean edot -- Neutrino and anti-neutrino loss rates [ergs/g] Yedot -- Change in electron fraction wrt time (neutronization rate) Ye -- Electron fraction dens -- Density in/out emq -- E-Q = Internal energy minus binding energy (per volume?)
Takes place after burning is complete. NSE expands after burning and the equilibrium state is a function of the density and the temperature. Been burned, flame is done, but then the hot stuff expands. This changes the composition and releases energy.