Kawano コードはSarkarのページのBig Bang NucleosynthesisのToolsというコーナーにある。Kawanoの "Let's go: Early universe 2. Primordial nucleosynthesis the computer way"というmanualもそこからspire経由で探せるようになっている。Arxivに置いてないので、Fermi Labのserverとかからとってこないといけない。

上の図のRefs.

\bibitem{Kawanocode}L. H. Kawano, FERMILAB-Pub-92/04-A (1992).

\bibitem{NACRE} C. Angulo {\it et al.}, Nucl. Phys. A. {\bf 656}, 3 (1999).

%\cite{Cyburt:2001pp}
\bibitem{Cyburt:2001pp}
R.~H.~Cyburt, B.~D.~Fields and K.~A.~Olive,
%``The NACRE Thermonuclear Reaction Compilation and Big Bang Nucleosynthesis,''
New Astron.\ {\bf 6}, 215 (2001)
%[arXiv:astro-ph/0102179].
%%CITATION = ASTRO-PH 0102179;%%

\bibitem{Fields1998} B. D. Fields and K. A. Olive, Astrophys. J. {\bf 506}, 177 (1998)

%\cite{Izotov:2003xn}
\bibitem{Izotov:2003xn}
Y.~I.~Izotov and T.~X.~Thuan,
%``Systematic effects and a new determination of the primordial abundance of
%4He and dY/dZ from observations of blue compact galaxies,''
arXiv:astro-ph/0310421.
%%CITATION = ASTRO-PH 0310421;%%

%\cite{Kirkman:2003uv}
\bibitem{Kirkman:2003uv}
D.~Kirkman, D.~Tytler, N.~Suzuki, J.~M.~O'Meara and D.~Lubin,
%``The cosmological baryon density from the deuterium to hydrogen ratio towards
%QSO absorption systems: D/H towards Q1243+3047,''
Astrophys. J., Suppl. Ser. {\bf 149}, 1 (2003)
%[arXiv:astro-ph/0302006].
%%CITATION = ASTRO-PH 0302006;%%

\bibitem{Ryan2000} S. G. Ryan, T. C. Beers, K. A. Olive, B. D. Fields, and J. E. Norris, Astrophys. J. Lett., {\bf 530}, L57 (2000)



上の絵にはのっていないが、

%\cite{Bonifacio:2002yx}
\bibitem{Bonifacio:2002yx}
P.~Bonifacio {\it et al.},
%``The lithium content of the globular cluster NGC 6397,''
Astron. Astrophys. {\bf 390}, 91 (2002)
%[arXiv:astro-ph/0204332].
%%CITATION = ASTRO-PH 0204332;%%



という論文でLi7のabundanceがreportされている。Abstract にはA(Li)=2.34 \pm 0.056の様にしかかいてないが、本文中のConclusion sectionに、「0.056はstatisticalでこれに0.06のsystematicを足さないといけない」となにげに書いてあるので注意。

追加したい計算：

DとLi7の理論エラーをLogでも出す（観測のエラーがそのようにreportされるから）。

Input files:

computationpara.dat（以下の値はKawanoのデフォルトのはず。めったに変更する必要なし）

0.3d0
0.03d0
1.0d-04
1.0d+02　開始温度（10^9 K単位)
1.0d-02　終了温度（10^9 K単位)
1.0d-25
30


modelpara.dat　（元のコードから結構変えた。今は整数値に限らないN_nuで計算できるようにしたため、すべてのニュートリノのxiが同じ場合しか計算できないようになっているはずなので、注意。）

1.0d0　　　重力定数（の普通の値との比）
885.7d0　　中性子の寿命
3.0d0　　　ニュートリノ世代数（runpara.datで制御しているので、ここでの値は使わない）
0.0d0　　eta（これもrunpara.datで制御しているので、ここでの値は使わない）
0.0d0　cosmological constant
1.0d-8 xi_e （これもrunpara.datで制御しているので、ここでの値は使わない）
0.0d0 xi_mu
0.0d0 xi_tau


runpara.dat　（各種ループやモンテカルロ用。標準的な値にしたい時は、minをその値にして、回数を1にしておく。xi_eがzeroの時は、コードの内部で1.d-8を代入するようになっている。このため、thermのところでweak interactionのfitting公式ではなく、数値積分が計算されるので、xi_e でループしたときnon-zeroのところとなめらかにつながるようになっている。）

3 BBNネットーワークの大きさ。今は一番小さい３でいい。
'MCend' 1000　　 MCする回数。１にすると、乱数をふらない。
1.d-10　　　　　eta_min
1.99526d-9 eta_max
0.0d0 対数目盛りでのetaの増分。ここがzeroなら下の等間隔eta幅で計算。
'eta' 4.d-10 4.d-12 101　etaのmin, inc, 回数　（max = min+inc*(回数-1)になるはず）
'xi_e' 0.d0 0.005d0 1　　xi（どのニュートリノ種についても共通）のmin, inc, 回数
'#_nu' 3.d0 1.d0 1　　　ニュートリノ種類数のmin, inc, 回数。（非整数も可）
'rAlpha' 1.0d0 0.01d0 1　fine structure constantの標準値との比のmin, inc, 回数。



etaを1^-10から1^-9までlogで20等分計算したいときは、３〜５行目を

1.d-10　　
1.d-9
0.05d0

とする。