This shows you the differences between two versions of the page.
Both sides previous revision Previous revision Next revision | Previous revision | ||
the_neutrino_winds-group [2014/06/09 17:33] richers |
the_neutrino_winds-group [2014/06/12 18:46] schutrumpg |
||
---|---|---|---|
Line 11: | Line 11: | ||
**Procedure** | **Procedure** | ||
+ | |||
* Calculate reaction rates given electron neutrino energy & flux, and electron anti-neutrino energy and flux | * Calculate reaction rates given electron neutrino energy & flux, and electron anti-neutrino energy and flux | ||
* modify XNet to include these reaction rates | * modify XNet to include these reaction rates | ||
Line 18: | Line 19: | ||
* Neutrino cross-sections for n and p (Burrows, Reddy, Thompson 2006) | * Neutrino cross-sections for n and p (Burrows, Reddy, Thompson 2006) | ||
* neutrino-driven wind trajectory (from Raph Hix, or arxiv.org/pdf/astro-ph/9912156.pdf, or arxiv.org/pdf/0908.1871.pdf) | * neutrino-driven wind trajectory (from Raph Hix, or arxiv.org/pdf/astro-ph/9912156.pdf, or arxiv.org/pdf/0908.1871.pdf) | ||
+ | * Git repository; https://github.com/srichers/NuWind.git | ||
- | **Adding neutrino rates to XNet* | + | **Adding neutrino rates to XNet** |
* net_preprocess.f90 | * net_preprocess.f90 | ||
* in the desc_known list, include an entry "nu11" and "au12" for neutrino and antineutrino (no reason for the numbers - arbitrary name) | * in the desc_known list, include an entry "nu11" and "au12" for neutrino and antineutrino (no reason for the numbers - arbitrary name) | ||
Line 28: | Line 30: | ||
* define formula for cross section for both neutrinos and antineutrinos in the switch statement if(iweak>0) | * define formula for cross section for both neutrinos and antineutrinos in the switch statement if(iweak>0) | ||
* rpf1 is set to inverse partition-function stuff if inverse rate, 1 if forward rate (determined by whether a "v" is in the reaction rate database entry | * rpf1 is set to inverse partition-function stuff if inverse rate, 1 if forward rate (determined by whether a "v" is in the reaction rate database entry | ||
- | * csect1 = cross section = nflux*rpf1*exp(matmul(t09,rc1)+h1) (t09,rc1,and h1 are read in from reaction database. all but h1 should be zero) | + | * csect1 = reaction rate = (number flux)*(cross section) |
- | * | + | |
+ | **Input Parameters** | ||
+ | |||
+ | Extracted input parameters from T. Fischer et al. using GSYS 2.4 (a program to extract data points from graphs) http://www.jcprg.org/gsys | ||
+ | |||
+ | **Simulations** | ||
+ | |||
+ | Ran simulations using xnet for different $Y_e$ values and different progenitor masses. We tried to see the effect of the presence of neutrinos and anti-neutrinos in the r-process abundance patern. |