AbstractWe aim to improve the accuracy of radiative energy transport in three-dimensional radiation hydrodynamical simulations in ANTARES (A Numerical Tool for Astrophysical RESearch). We implement in the ANTARES short-characteristics numerical schemes a modification of the Bézier interpolant solver. This method yields a smoother surface structure in simulations of solar convection and reduces the artifacts appearing due to the limited number of rays along which the integration is done. Reducing such artifacts leads to increased stability of the code. We show that our new implementation achieves a better agreement of the temperature structure and its gradient with a semi-empirical model derived from observations, as well as of synthetic spectral-line profiles with the observed solar spectrum.