Radiative Relativistic Magnetohydrodynamic Simulations Of Neutron Star Column Accretion