User Tools

Site Tools


the_neutrino_winds-group

Differences

This shows you the differences between two versions of the page.

Link to this comparison view

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 databaseall but h1 should be zero+    * csect1 = reaction rate = (number flux)*(cross section) 
-    + 
 +**Input Parameters** 
 + 
 +Extracted input parameters ​from TFischer et al. using GSYS 2.4 (a program to extract data points from graphshttp://​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.
the_neutrino_winds-group.txt · Last modified: 2014/06/12 18:46 by schutrumpg