Radiative transfer with scattering for domain-decomposed 3D MHD simulations of cool stellar atmospheres

We present the implementation of a radiative transfer solver with coherent scattering in the new BIFROST code for radiative magneto-hydrodynamical (MHD) simulations of stellar surface convection. The code is fully parallelized using MPI domain decomposition, which allows for large grid sizes and improved resolution of hydrodynamical structures. We apply the code to simulate the surface granulation in a solar-type star, ignoring magnetic fields, and investigate the importance of coherent scattering for the atmospheric structure. A scattering term is added to the radiative transfer equation, requiring an iterative computation of the radiation field. We use a short-characteristics-based Gauss-Seidel acceleration scheme to compute radiative flux divergences for the energy equation. The effects of coherent scattering are tested by comparing the temperature stratification of three 3D time-dependent hydrodynamical atmosphere models of a solar-type star: without scattering, with continuum scattering only, and with both continuum and line scattering. We show that continuum scattering does not have a significant impact on the photospheric temperature structure for a star like the Sun. Including scattering in line-blanketing, however, leads to a decrease of temperatures by about 350,K below log tau < -4. The effect is opposite to that of 1D hydrostatic models in radiative equilibrium, where scattering reduces the cooling effect of strong LTE lines in the higher layers of the photosphere. Coherent line scattering also changes the temperature distribution in the high atmosphere, where we observe stronger fluctuations compared to a treatment of lines as true absorbers.

Chemical similarities between Galactic bulge and local thick disk red giant stars

The evolution of the Milky Way bulge and its relationship with the other Galactic populations is still poorly understood. The bulge has been suggested to be either a merger-driven classical bulge or the product of a dynamical instability of the inner disk. To probe the star formation history, the initial mass function and stellar nucleosynthesis of the bulge, we performed an elemental abundance analysis of bulge red giant stars. We also completed an identical study of local thin disk, thick disk and halo giants to establish the chemical differences and similarities between the various populations. High-resolution infrared spectra of 19 bulge giants and 49 comparison giants in the solar neighborhood were acquired with Gemini/Phoenix. All stars have similar stellar parameters but cover a broad range in metallicity. A standard 1D local thermodynamic equilibrium analysis yielded the abundances of C, N, O and Fe. A homogeneous and differential analysis of the bulge, halo, thin disk and thick disk stars ensured that systematic errors were minimized. We confirm the well-established differences for [O/Fe] (at a given metallicity) between the local thin and thick disks. For the elements investigated, we find no chemical distinction between the bulge and the local thick disk, which is in contrast to previous studies relying on literature values for disk dwarf stars in the solar neighborhood. Our findings suggest that the bulge and local thick disk experienced similar, but not necessarily shared, chemical evolution histories. We argue that their formation timescales, star formation rates and initial mass functions were similar.