Probing Na in giant exoplanets with ESPRESSO and 3D NLTE stellar spectra

Neutral sodium was the first atom detected in an exoplanetary atmosphere via transmission spectroscopy and remains the most frequently detected species due to its strong doublet in the optical. However, the center-to-limb variation (CLV) of these lines in the host star can bias the Na detection.When combined with the Rossiter-McLaughlin (RM) effect, the CLV can mimic or obscure planetary absorption features. This work investigates the impact of 3D radiation hydrodynamic stellar atmospheres and non-local thermodynamic equilibrium (NLTE) radiative transfer on the modeling of the CLV+RM effect in single-line transmission spectroscopy, to improve the characterization of exoplanet atmospheres. We produced a grid of 3D NLTE synthetic spectra for Na I for FGK stars within the following parameter space: Teff=4500-6500 K, log g=4.0-5.0 and [Fe/H]=-0.5, 0, 0.5.This grid was then interpolated to match the stellar parameters of four stars hosting giant exoplanets, to correct for the CLV+RM effect in their transmission spectra. We used ESPRESSO archival observations. Our work confirms the Na detections in three systems, namely WASP-52b, WASP-76b, and WASP-127b, improving the accuracy of the measured absorption depth. Furthermore, we find that 3D NLTE stellar models can explain the spectral features in HD 209458b without the need for any planetary absorption. In the grid of synthetic spectra, we observe that the CLV effect is stronger for stars with low Teff and high log g. However, the combined effect of CLV and RM is highly dependent on the orbital geometry of the planet-star system. With the continuous improvement of instrumentation, it is crucial to use the most accurate stellar models to correct for the CLV+RM effect in high-resolution transmission spectra, to achieve the best possible characterization of exoplanet atmospheres. We make our grid of 3D NLTE spectra for Na publicly available.