Evan edited untitled.tex  over 8 years ago

Commit id: 1b9f42904d30039cf4d868dbbf9640900e29f6af

deletions | additions      

       

The goal is to create a simple yet complete set of conditions that can be easily simulated by many core-collapse groups. We want to strictly adhere to the same procedure for everything to minimize differences. Below we outline in detail the conditions for our tests, please comment if these are too strict, not strict enough, or if other conditions need to be considered.  \section{Initial Conditions}  \subsection{Non-neutrino Physics}  \begin{enumerate}  \item Progenitor: $20\,M_\odot$ model from \cite{2015arXiv151004643S}. This model is publicly available from \url{http://www.ucolick.org/~sukhbold/#iii}, and as stated on this website, will be available from ApJ after publication. This openly available, up to date progenitor model is ideal for this study.  \item For mapping the progenitor, we will use density, temperature, and ye. Careful to note the definition of the radial coordinate in the initial model (radial coordinate is the location of the outer edge of the zone, the velocity is also defined at this radius, the remaining required quantities (rho, temperature, ye) are defined as cell averages).  \item  Equation of state: The SFHo nuclear equation of state fromMattias  Hempel et al. available from \url{http://phys-merger.physik.unibas.ch/~hempel/eos.html} or \url{http://www.stellarcollapse.org/equationofstate}. This EOS extends down to densities of $1\,\mathrm{g}\,\mathrm{cm}^{-3}$. In this EOS, NSE is assumed down to these densities (which is incorrect for supernovae). However, to eliminate (for now) issues related to low density equations of state and nuclear reaction networks, this study will use only the SFHo equation of state for all densities, temperatures, and ye's. \item Boundary and Boundary Conditions: The SFHo EOS only goes down to $0.1\,$MeV, therefore the outer boundary must be closer than $1.2\times 10^9\,\mathrm{cm}$. For this comparison, we would like the outer boundary to be taken as $10^9\,\mathrm{cm}$. For the outer  boundary conditions, fix the density and velocity so as to maintain a constant mass accretion rate. This is not the most physical boundary condition, but ensures the same condition is used by different groups. \item For multidimensional simulations, do collapse in 1D up to 15\,ms after bounce, then transition to 2D. Add perturbations via a clearly defined manner.  \item Perform simulations using both Newtonian gravity and some form of general relativistic gravity (either effective potential, true GR).  \end{enumerate}  \subsection{Neutrino Physics}  \begin{enumerate}  \item The SFHo EOS contains light clusters. For this study, all neutrino interactions (both absorption and scattering) on light clusters is ignored.  \item For scattering on free nucleons, use the Breunn (1985) rates. Include weak magnetism and recoil corrections.  \item For scattering on heavy nuclei, use the Breunn (1985) rate, include corrections a la Burrows, Reddy, and Thompson (2006): Ion-ion correlations, electron polarization correction, and a correction for the nuclear form factor.  \item Include inelastic neutrino-electron scattering only in 1D simulations. This will include inelastic scattering during the collapse phase of 2D simulations as up to 15ms postbounce will be done in 1D.  \item For pair processes, include electron-positron annihilation to neutrino pairs (and the reverse). Ignore all other pair processes, nucleon-nucleon Bremsstrahlung and neutrino pair conversion to other neutrino pairs.  \item Use Bruenn (1985) rates for absorption rates (on nucleons and nuclei). Include weak magnetism and recoil corrections for nucleon rates. Do not include any nucleon potentials.  \end{enumerate}  \section{Results}