Modelling Barium Isotopes in Metal-Poor Stars
Abstract
The principal theory concerning the origin of the elements heavier than the Fe-peak,
such as Ba, strongly suggest that for old, metal-poor environments, the rapid (r-)
process is the most likely path taken in their synthesis, while the slow (s-) process
becomes more substantial in younger, more metal-rich stellar populations.
In this work I test this theory by evaluating the isotope ratios of Ba. It is understood
that Ba consists of seven stable isotopes, five of which are synthesised by the
two neutron-capture processes. The two odd isotopes, 135,137Ba, as well as 138Ba are
synthesised via both the r- and s-processes while two of the even isotopes, 134,136Ba are
synthesised via the s-process only. The relative contribution of the r- and s-process
to these isotopes can be understood via nucleosynthesis calculations and is described
using the parameter fodd, where fodd = [N (135Ba) + N (137Ba)] /N (Ba). Low values of
fodd (~0.11) indicate an s-process regime, while high values of fodd (~0.46) indicate
an r-process regime.
In the Ba II 4554 A line the even isotopes lie close to the line centre, while the
odd isotopes, which are hyperfine split because of their non-zero nuclear spin, lie in
the wings of the line. From an analysis of the line profile shape, one can determine
whether Ba has been synthesised primarily through the r-process or s-process; a broad,
asymmetric line would indicate a high r-process contribution, while a line with a deeper
core and shallower wings would indicate a high s-process contribution.
Using the radiative transfer code ATLAS, which assumes local thermodynamic
equilibrium (LTE) and employs 1-dimensional (1D) KURUCZ06 model atmospheres,
I synthesised line profiles for six metal-poor stars: HD140283, HD122563, HD88609,
HD84937, BD-04 3208 and BD+26 3578 - for a range of isotope ratios. All six
are of sufficiently low metallicity that Ba was expected to have an r-process origin.
These were fit to high resolution (R\equiv \lamda/\Delta\lamda = 90 000 - 95 000), high signal-to-noise to the Ba II 4554 A line which has multiple components. In the first test, synthetic
spectra were computed using the non local thermodynamic equilibrium (NLTE) radiative
transfer code MULTI. The synthetic line profiles were fit to a number of lines in
HD140283. Although this technique might have improved the fit in the line core, it
was found that such a treatment did not improve upon fitting errors associated with
the best fit 1D LTE synthetic profiles.
The second test used a 3-dimensional (3D) radiative transfer code (LINFOR3D) that
employed 3D, time-dependent atmospheres produced with CO5BOLD. The 3D synthetic
pro les were fit to a selection of Fe lines and improvements over the poor fits produced
by the 1D LTE synthesis were seen. It was found that the 3D synthesis could almost
completely reproduce the line asymmetries seen in the observed stellar spectrum. This
result suggests that further work to refine the 3D calculations and synthesis code would
be valuable.