NTOFPublic Web>WebPreferences>Be7npPaperDraft>Be7npPaperDraftRates (2018-05-14, AlbertoMengoni)

List of the 12th most important reactions in the BBN network

	code index	reaction	adopted	options
1	12	$^{1}H(n,\gamma)D$	`ando`	`ando skm`
2	16	$^{3}He(n,p)T$	`stlb`	`stlb de04 skm`
3	17	$^{7}Be(n,p)^{7}Li$	`ntof`	`ntof cy04 skm`
4	19	$^{7}Be(n,\alpha)^{4}He$	`ntof`	`wag ntof`
5	20	$D(p,\gamma)^{3}He$	`il16`	`il16 skm`
6	24	$^{7}Li(p,\alpha)^{4}He$	`stlb`	`stlb cf88 de04 skm`
7	26	$T(\alpha,\gamma)^{7}Li$	`stlb`	`stlb skm`
8	27	$^{3}He(\alpha,\gamma)^{7}Be$	`il16`	`il16 ncr2 skm cd08`
9	28	$D(d,n)^{3}He$	`stlb`	`stlb skm`
10	29		`stlb`	`stlb skm`
11	30	$T(d,n)^{4}He$	`stlb`	`stlb de04 cf88 skm`
12	31	$^{3}He(d,p)^{4}He$	`stlb`	`stlb de04 skm`

#wag  : Wagoner, R.V., ApJS, 18, 247 (1969)
#cf88 : Caughlan and Fowler, ADNDT 40, 283 (1988)
#skm  : Smith, Kawano and Malaney, ApjS 85, 2019 (1993)
#de04 : Descouvemont et al. ADNDT 88, 203 (2004)
#cy04 : Cyburt, PRD 70, 023505 (2004)
#ando : Ando et al. Phys. Rev. C 74, 025809 (2006)
#il16 : Iliadis et al., ApJ 831, 107 (2016)
#ncr2 : fit to ncare2 library
#stlb : fit to starlib table
#ntof : rate from n_TOF experiments

Code segments of the individual reaction rates adopted in the present BBN calculations

Reaction	12	$^{1}H(n,\gamma)D$	plot	table

    1/* H + n -> g + H2 */
    2/* rate from Ando et al. Phys. Rev. C 74, 025809 (2006) */
    3/* p + n -> d + g */
    4        a0 = 44216.0; a1 = 3.75191; a2 = 1.92934; a3 = 0.746503; a4 = 0.0197023; a5 = 3.00491e-6;
    5        b1 = 5.4678; b2 = 5.62395; b3 = 0.489312; b4 = 0.00747806;
    6        f[12]=a0*(1. + a1*T9 + a2*pow(T9,2.) + a3*pow(T9,3.) + a4*pow(T9,4.) + a5*pow(T9,5.))/(1. + b1*T9 + b2*pow(T9,2.) + b3*pow(T9,3.) + b4*pow(T9,4.));

Reaction	16	$^{3}He(n,p)T$	plot	table

    1* He3 + n -> p + H3 */
    2/* starlib rate parameters -  de04 data */
    3        a0 = 6.89423e+08;   a1 = -0.4347;      a2 = -0.199323;    a3 = 0.590373;     a4 = -0.322449;    a5 = 0.0742735;    a6 = -0.00648816;
    4        f[16] = a0*(1. + a1*pow(T9,1./2.) + a2*T9 + a3*pow(T9,3./2.) + a4*pow(T9,2.) + a5*pow(T9,5./2.) + a6*pow(T9,3.));

Reaction	17	$^{7}Be(n,p)^{7}Li$	table

    1* Be7 + n -> p + Li7 */
    2/* rate from n_TOF  - central value */
    3        a0=   6.8048e+09; a1= -1.9706e+00; a2=  2.0419e+00; a3= -1.0687e+00; a4=  2.7172e-01; a5= -2.6699e-02; a6=  1.9610e+08; a7=  2.8899e+07; b0=  2.8114e-01;
    4        f[17] = a0*(1. + a1*pow(T9,1./2.) + a2*T9 + a3*pow(T9,3./2.) + a4*T9*T9 + a5*pow(T9,5./2.)) + a6*pow(T9/(1.+13.076*T9),3./2.)/pow(T9,3./2.) + a7/pow(T9,3./2.)*exp(-b0/T9); 
    5
    6/* rate from n_TOF  - upper */
    7        a0=   6.8181e+09; a1= -1.6226e+00; a2=  1.4457e+00; a3= -6.7629e-01; a4=  1.5737e-01; a5= -1.4373e-02; a6=  7.3135e+08; a7= -2.8413e+08; b0=  8.1347e-01; 
    8        f[17] = a0*(1. + a1*pow(T9,1./2.) + a2*T9 + a3*pow(T9,3./2.) + a4*T9*T9 + a5*pow(T9,5./2.)) + a6*pow(T9/(1.+13.076*T9),3./2.)/pow(T9,3./2.) + a7/pow(T9,3./2.)*exp(-b0/T9); 
    9
   10/* rate from n_TOF  - lower */
   11        a0=   7.5185e+09; a1= -2.4307e+00; a2=  2.9935e+00; a3= -1.7338e+00; a4=  4.6746e-01; a5= -4.7572e-02; a6= -9.7708e+08; a7= -7.0378e+09; b0=  3.1868e+00; 
   12        f[17] = a0*(1. + a1*pow(T9,1./2.) + a2*T9 + a3*pow(T9,3./2.) + a4*T9*T9 + a5*pow(T9,5./2.)) + a6*pow(T9/(1.+13.076*T9),3./2.)/pow(T9,3./2.) + a7/pow(T9,3./2.)*exp(-b0/T9);

Reaction	19	$^{7}Be(n,\alpha)^{4}He$	table

    1/* Be7 + n -> a + He4 */
    2/* n_TOF rate as in Barbagallo et al. PRL 117 (2016) 152701 */
    3        f[19]=4.81e+5 + 1.84e+6*T9 + 3.03e+6*pow(T9,1.5);

Reaction	20	$D(p,\gamma)^{3}He$	plot	table

    1/* H2 + p -> g + He3 */	
    2/* rate from Iliadis et al. (2016) fit */
    3        a0=    2.452e+03; b0=    3.732e+00; a1=  -1.632e+00; a2=   7.609e+00; a3=  -6.133e+00; a4=   6.557e+00; a5=  -1.042e+00;
    4        f[20]=a0*pow(T9,-2./3.)*exp(-b0*pow(T9,-1./3.))*(1. + a1*pow(T9,1./3.) + a2*pow(T9,2./3.) + a3*T9 + a4*pow(T9,4./3.) + a5*pow(T9,5./3.));

Reaction	24	$^{7}Li(p,\alpha)^{4}He$	plot	table

    1/* Li7 + p -> a + He4 */
    2/* starlib rate parameters -  de04 data */
    3        a0 = 1.45251e+09;    b0 = 8.44454;       a1 = -6.76368e+08;  b1 = 8.41908;       a2 = 4.0217e+09;    b2 = 25.3862;       b3 = 0.16592;       
    4        f[24] = a0/pow(T9,2./3.)*exp(-b0/pow(T9,1./3.)) + a1*pow(T9/(1.0+b3*T9),5./6.)/pow(T9,3./2.)*exp(-b1/pow(T9/(1.0+b3*T9),1./3.))+a2/pow(T9,3./2.)*exp(-b2/T9);

Reaction	26	$T(\alpha,\gamma)^{7}Li$	plot	table

    1/* H3 + a -> g + Li7  original */
    2/* starlib rate parameters -  in17 data */
    3        a0=   1.3350e+05; b0=  8.0768e+00; a1=  5.9164e+00; a2= -1.3045e+01; a3=  1.3315e+01; a4= -6.4673e+00; a5=  1.3420e+00; a6=  5.1090e+05; b1=  8.0680e+00; 
    4        f[26]=a0/pow(T9,2./3.)*exp(-b0/pow(T9,1./3.))*(1.+a1*pow(T9,1./3.)+a2*pow(T9,2./3.) + a3*T9 + a4*pow(T9,4./3.) + a5*pow(T9,5./3.)) + a6*pow(T9/(1.+0.1378*T9),5./6.)/pow(T9,3./2.)*exp(-b1/pow(T9/(1.+0.1378*T9),1./3.));

Reaction	27	$^{3}He(\alpha,\gamma)^{7}Be$	plot	table

    1/* He3 + a -> g + Be7 */
    2/* rate from Iliadis et al. (2016) fit */
    3        a0=   9.8079e+06; b0=   1.4192e+01; a1= -1.3272e+00; a2=  1.9911e-01; a3=  1.0387e+00; a4= -6.7320e-01; a5=  1.2511e-01; a6=  6.1164e+06; b1=  1.2829e+01;
    4        T9F=T9/(1.+0.1071*T9);
    5        f[27]= a0*pow(T9,-2./3.)*exp(-b0/pow(T9,1./3.))*(1.+a1*pow(T9,1./3.) + a2*pow(T9,2./3.) + a3*T9 + a4*pow(T9,4./3.) + a5*pow(T9,5./3.)) + a6*pow(T9F,5./6.)*pow(T9,-3./2.)*exp(-b1/pow(T9F,1./3.));

Reaction	28	$D(d,n)^{3}He$	plot	table

    1/* H2 + H2 -> n + He3 */
    2/* starlib rate parameters -  in17 data */
    3        a0 = 3.95e+8; b0 = 4.259; a1 = 0.098; a2 = 1.00239; a3 = 1.08419; a4 =-0.743092; a5 = 0.139404;       
    4        f[28]=a0*pow(T9,-2./3.)*exp(-b0/pow(T9,1./3.))*(1.+a1*pow(T9,1./3.)+a2*pow(T9,2./3.)+a3*T9+a4*pow(T9,4./3.)+a5*pow(T9,5./3.));

Reaction	29		plot	table

    1/* H2 + H2 -> p + H3 */
    2/* starlib rate parameters -  in17 data */
    3        a0 = 4.12649e+08; b0 = 4.258; a1 = 0.304711; a2 = -0.446767; a3 = 2.23021; a4 = -1.19482; a5 = 0.210592;   
    4        f[29]= a0*pow(T9,-2./3.)*exp(-b0/pow(T9,1./3.))*(1.+a1*pow(T9,1./3.)+a2*pow(T9,2./3.)+a3*T9+a4*pow(T9,4./3.)+a5*pow(T9,5./3.));

Reaction	30	$T(d,n)^{4}He$	plot	table

    1/* H3 + H2 -> n + He4 */
    2/* starlib rate parameters -  de04 data */
    3        a0 = 6.32853e+10; b0 = 4.52388; b1 = 0.072028;     a1 = 6.47865;      a2 = -29.3575;     a3 = 4.35456;      a4 = 280.001;      a5 = -253.094;     a6 = 7.56418e+08;  b2 = 0.472317;     
    4        f[30]=a0*pow(T9,-2./3.)*exp(-b0/pow(T9,1./3.)-pow(T9/b1,2.))*(1. + a1*pow(T9,1./3.) + a2*pow(T9,2./3.) + a3*T9 + a4*pow(T9,4./3.) + a5*pow(T9,5./3.)) + a6*pow(T9,-2./3.)*exp(-b2/T9);

Reaction	31	$^{3}He(d,p)^{4}He$	plot	table

    1/* He3 + H2 -> p + He4 */
    2/* starlib rate parameters -  de04 data */
    3        a0 = 2.46165e+10; b0 = 7.15127; b1 = 0.241809; a1 = 8.69492; a2 = -22.7967; a3 = 21.7527; a4 = 6.97; a5 = 34.9149; a6 = 4.82217e+08; b2 = 1.71366;
    4        f[31]=a0*pow(T9,-2./3.)*exp(-b0/pow(T9,1./3.)-pow(T9/b1,2.))*(1.+ a1*pow(T9,1./3.)+ a2*pow(T9,2./3.) + a3*T9 + a4*pow(T9,4./3.) + a5*pow(T9,5./3.)) + a6/pow(T9,1./2.)*exp(-b2/T9);

-- AlbertoMengoni - 2018-02-20

Topic revision: r10 - 2018-05-14 - AlbertoMengoni

NTOFPublic

Public webs

- Cern Search
- TWiki Search
- Google Search
NTOFPublic All webs

Copyright &© 2008-2023 by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
or Ideas, requests, problems regarding TWiki? use Discourse or Send feedback