The internal topology, strength and dynamics of magnetic fields in strongly magnetised neutron stars play a crucial role in the secular and dynamical evolution of these exotic objects. I will present both analytic calculations and nonlinear simulations of strongly magnetised neutron stars, with particular emphasis on determining stable magnetic field equilibrium states. The existence of various new branches of solutions in the context of magnetar physics and gravitational wave observations will be presented.