Abstract. The nuclear reaction code EMPIRE is presented as a useful tool for nuclear astrophysics. EMPIRE combines a variety of the reaction models with a comprehensive library of input parameters providing a diversity of options for the user. With exclusion of the directsemidirect capture all reaction mechanisms relevant to the nuclear astrophysics energy range of interest are implemented in the code. Comparison to experimental data show consistent agreement for all relevant channels.
IntroductionFor the majority of elements heavier than iron, nucleosynthesis occurs mainly via neutron capture reactions in the r-and s-processes. The fission channel needs to be considered in the r-process for heavy nuclei (Z > 80) [1] in which fission barriers are low. Of particular importance are the superheavy elements (SHE) since they may determine the upper end of the nucleosynthesis flow [2] which may only be reached via certain nuclear physics inputs [3]. Although the calculated and measured fission cross sections globally differ by less than a factor of three [4], current extrapolations for r-process simulation models are still unreliable. The neutron spectrum typical of s-process sites can be described by a Maxwell-Boltzmann distribution (MBD). Neutron captures have their relevant energies around the maximum of the Maxwell-Boltzmann distribution (MBD), E o ≈ kT [5]. Stellar neutron capture cross sections in the mass region of the weak s-process are often not available with the required accuracy for abundance determination. The main mechanism of synthesis for the p-nuclei are photodisintegration reactions. A reprocessing of initial s-seed distribution [6] occurs via (γ,n) reactions with subsequent alpha (γ,α) or proton emission (γ,p) branchings. Simulation models have had problems reproducing the abundance distribution for the p-nuclei for Type II supernova [7]. The situation is ameliorated when considering Type 1a supernovae as a site of synthesis [8].The heavier nuclei are formed via reactions at a high enough level density such that the Hauser Feshbach (HF) formulation applies. EMPIRE [9] is a nuclear reaction model code that accepts nucleons, γ's, and several light ions (including α-particles and deuterons) in the incoming and outgoing channels. In the presence of direct reactions, the Engelbrecht-Weidenmüller transformation can be invoked to account for the interference between compound nucleus (CN) and direct. Resonance treatment (r-matrix formalism) is beyond the scope of the code and will therefore be ignored (neutron resonances [10] can be invoked from the resonance module). The bulk of the incoming flux is treated in terms of the HF formulation with the width fluctuation correction accounted for by the HRTW [11] model. The incident channel can be treated in terms