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
Next revision Both sides next revision
the_neutrino_winds-group [2014/06/09 16:27]
richers
the_neutrino_winds-group [2014/06/09 17:52]
richers
Line 7: Line 7:
  
 Simulate neutrino interactions changing the composition of a neutrino-driven wind. Simulate neutrino interactions changing the composition of a neutrino-driven wind.
 +  * $\nu_e$ + n --> e + p
 +  * $\bar{\nu}_e$ + p --> e<​sup>​+</​sup>​ + n
  
 **Procedure** **Procedure**
Line 16: Line 18:
   * 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**
 +  * 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)
 +    * increase ndesc by 2 (152-->​154) to account for additional entries
 +    * introduce new case in idesc switch statement (we called it case 153 and 154). set iec to 11 and 12 (matching the names above). Farther down, iwk1(n) is set to iec, where n is the number of the reaction.
 +  * common.f90
 +    * iweak>0 means include all (including weak) interactions. iweak<0 means ONLY weak interactions
 +    * 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
 +    * 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)
 +    * 
the_neutrino_winds-group.txt ยท Last modified: 2014/06/12 18:46 by schutrumpg