[Paper Review] 3D Simulation of Convection and Spectral Line Formation in A-type Stars

M. Steffen, B. Freytag+1 more|arXiv (Cornell University)|Sep 15, 2005

Stellar, planetary, and galactic studies3 citations

TL;DR

This study presents the first 3D radiative hydrodynamics simulations of convection in A-type stars using the CO 5 BOLD code, revealing substantial convective overshoot below the He II zone and a solar-like granulation pattern with larger, less structured granules. Synthetic line profiles show a red-depressed wing inconsistent with observed inverse C-shaped bisectors in Am stars, suggesting missing physics like magnetic fields in current models.

ABSTRACT

We present first realistic numerical simulations of 3D radiative convection in the surface layers of main sequence A-type stars with Teff = 8000 K and 8500 K, log g = 4.4 and 4.0, recently performed with the CO5BOLD radiation hydrodynamics code. The resulting models are used to investigate the structure of the H+HeI and the HeII convection zones in comparison with the predictions of local and non-local convection theories, and to determine the amount of "overshoot" into the stable layers below the HeII convection zone. The simulations also predict how the topology of the photospheric granulation pattern changes from solar to A-type star convection. The influence of the photospheric temperature fluctuations and velocity fields on the shape of spectral lines is demonstrated by computing synthetic line profiles and line bisectors for some representative examples, allowing us to confront the 3D model results with observations.

Motivation & Objective

• To investigate the structure and dynamics of shallow convection zones in A-type stars using 3D hydrodynamical simulations.
• To test the validity of local and non-local convection theories against realistic 3D simulations.
• To quantify convective overshoot into stable layers beneath the He II convection zone.
• To examine how photospheric temperature and velocity fluctuations affect spectral line profiles and bisectors.
• To resolve discrepancies between observed inverse C-shaped line bisectors in slowly rotating A-type stars and synthetic profiles from 3D models.

Proposed method

• Numerical simulations performed with the CO 5 BOLD radiation hydrodynamics code, solving the Navier-Stokes equations with a Roe-type approximate Riemann solver and van Leer reconstruction.
• Radiative transfer solved using a modified Feautrier scheme on long rays with grey approximation and a Phoenix-OPAL Rosseland mean opacity table.
• Equations of state include partial ionization of H and He and H₂ formation, with strict LTE and no scattering assumed.
• Simulations use periodic lateral boundaries and closed top/bottom boundaries, with vertical grid refinement to resolve pressure scale heights.
• Synthetic line profiles and bisectors computed for Fe I λ6265.13 Å under vertical and integrated line-of-sight conditions.
• Statistical analysis via horizontal and temporal averaging of 3D velocity, entropy, and energy flux structures.

Experimental results

Research questions

• RQ1 How do 3D hydrodynamical simulations of A-type star convection compare with predictions from local mixing-length theory (MLT) and non-local convection models?
• RQ2 What is the extent and nature of convective overshoot below the He II convection zone in A-type stars?
• RQ3 How does the granulation topology in A-type stars differ from solar granulation in terms of size, structure, and filling factor?
• RQ4 To what extent do 3D model-predicted line profiles and bisectors match observed inverse C-shaped bisectors in slowly rotating A-type stars?
• RQ5 What physical mechanisms might explain the discrepancy between synthetic profiles and observations if magnetic fields or non-LTE effects are absent?

Key findings

• Local mixing-length theory (MLT) fails to reproduce the 3D hydrodynamical results, regardless of the mixing-length parameter α.
• Non-local convection theory by Kupka & Montgomery (2002) shows qualitative similarity but still exhibits significant differences in energy fluxes compared to 3D simulations.
• Overshoot below the He II convection zone is substantial, with a velocity scale height of $ H_v/H_p ightarrow 0.4 $ in model 2 and $ H_v/H_p ightarrow 0.7 $ in model 1.
• The granulation pattern in A-type stars is qualitatively similar to solar granulation but features larger granules and a lower dark area filling factor ($ f_d ightarrow 0.34 $) compared to the Sun ($ f_d ightarrow 0.53 $).
• Synthetic line bisectors from 3D models exhibit a solar-like C-shape with red-wing depression, opposite to the observed inverse C-shape in Am stars with $ T_{ m eff} ightarrow 8000 $ K.
• The observed line asymmetry (flux bisector excursion ~2 km/s in Am stars) is of the same order of magnitude as in 3D models, but the direction of asymmetry is reversed, indicating a fundamental mismatch in current models.