Resolve the origin of discrepancies in predicted and simulated column heights

Determine whether the discrepancy between the Basko–Sunyaev sinking-regime predictions and the shorter accretion-column heights obtained in relativistic radiation magnetohydrodynamic simulations arises from limitations of the simulations or from the actual physics of magnetized accretion columns.

Background

The appendix compares the Basko–Sunyaev advective sinking solution with relativistic radiation magnetohydrodynamic simulations for moderately magnetized neutron stars. The simulated columns are systematically shorter than the analytical sinking-regime predictions, and the discrepancy persists even after accounting for a reduced opacity.

The authors note that the simulations impose specific opacity limits and bottom-boundary conditions, including a non-zero downward thermal flux, and do not reach the critical radiation-pressure boundary condition assumed by the sinking solution. It therefore remains unresolved whether the disagreement is caused by these simulation limitations or reflects physical effects absent from the analytical model.