Next-generation spectro-polarimetric broadband surveys will probe cosmic magnetic fields in unprecedented detail, using the magneto-optical effect known as Faraday rotation. However, non-parametric methods such as RMCLEAN can introduce non-observable linearly polarised flux into a fitted model at negative wavelengths squared. This leads to Faraday rotation structures that are consistent with the observed data, but would be impossible or difficult to measure. We construct a convex non-parametric QU-fitting algorithm to constrain the flux at negative wavelengths squared to be zero. This allows the algorithm to recover structures that are limited in complexity to the observable region in wavelength squared. We verify this approach on simulated broadband data sets where we show that it has a lower root mean square error and that it can change the scientific conclusions for real observations. We advise using this prior in next-generation broadband surveys that aim to uncover complex Faraday depth structures. We provide a public Python implementation of the algorithm at https://github.com/Luke-Pratley/Faraday-Dreams.

Broadband X-ray data of five magnetars show that their hard X-ray pulses suffer periodic phase modulations, at a period ∼ 104 times their pulse period. The phenomenon is interpreted as a result of free precession of neutron stars that are prolately deformed to an asphericity of ∼ 10−4, by the magnetic stress of toroidal fields reaching ∼ 1016 G. The behavior is absent in their soft X-ray pulses, probably due to a higher emission symmetry. The ultra-high toroidal fields, considered common to magnetars, may persist longer than their dipole fields.