VMR of our best fit atmospheric model against the dataset from Benneke et al. (2019) at nominal irradiation (175 (Z/H)_⊙, K_zz = 10⁶ cm² s⁻¹). For clarity, only absorbing species and N₂ are represented. FeH, TiO, and VO are not represented as their respective abundances are, respectively, < 10⁻¹², < 10⁻³⁵, and < 10⁻³⁶ within our pressure grid. The bump of the PH₃ VMR at 200 kPa is due to the chemical equilibrium between PH₃ and P₂ peaking in favour of P₂ at this pressure.