Monte-Carlo radiation hydrodynamic simulations of line-driven disc winds: relaxing the isothermal approximation

Disc winds play a crucial role in many accreting astrophysical systems across all scales. In accreting white dwarfs (AWDs) and active galactic nuclei (AGN), radiation pressure on spectral lines is a promising wind-driving mechanism. However, the efficiency of line driving is extremely sensitive to the ionization state of the flow, making it difficult to construct a reliable physical picture of these winds. Recently, we presented the first radiation-hydrodynamic (RHD) simulations for AWDs that incorporated detailed, multi-dimensional ionization calculations via fully frequency-dependent radiative transfer, using the Sirocco code coupled to PLUTO. These simulations produced much weaker line-driven winds (Mdot_wind / Mdot_acc < 1e-5 for our adopted parameters) than earlier studies using more approximate treatments of ionization and radiative transfer (which yielded Mdot_wind / Mdot_acc ~ 1e-4). One remaining limitation of our work was the assumption of an isothermal outflow. Here, we relax this by adopting an ideal gas equation of state and explicitly solving for the multi-dimensional temperature structure of the flow. In the AWD setting, accounting for the thermal state of the wind does not change the overall conclusions drawn from the isothermal approximation. Our new simulations confirm the line-driving efficiency problem: the predicted outflows are too highly ionized, meaning they neither create optimal driving conditions nor reproduce the observed ultraviolet wind signatures. Possible solutions include wind clumping on sub-grid scales, a softer-than-expected spectral energy distribution, or additional driving mechanisms. With the physics now built into our simulations, we are well-equipped to also explore line-driven disc winds in AGN.