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] (current) 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. | ||