Force-free pulsar magnetospheres develop a large scale poloidal electric current circuit that flows along open magnetic field lines from the neutron star to the termination shock. The electric current closes through the interior of the neutron star where it provides the torque that spins-down the star. In the present work, we study the internal electric current in an axisymmetric rotator. We evaluate the path of the electric current by requiring the minimization of internal Ohmic losses. We find that, in millisecond pulsars, the current reaches the base of the crust, while in pulsars with periods of a few seconds, the bulk of the electric current does not penetrate deeper than about $100$ m. The region of maximum spin-down torque in millisecond pulsars is the base of the crust, while in slowly spinning ones it is the outer crust. We evaluate the corresponding Maxwell stresses and find that, in typical rotation-powered radio pulsars, they are well below the critical stress that can be sustained by the crust. For magnetar-level fields, the Maxwell stresses near the surface are comparable to the critical stress and may lead to the decoupling of the crust from the rest of the stellar rotation.

We compute multiprecision solutions of the Lane-Emden equation. This differential equation arises when introducing the well-known polytropic model into the equation of hydrostatic equilibrium for a nondistorted star. Since such multiprecision computations are time-consuming, we apply to this problem parallel programming techniques and thus the execution time of the computations is drastically reduced.

In this study, we consider the problem of geometry of curved jets. We work with a simplified two-dimensional kinematical model assuming that the jet velocity is varying. We prescribe the velocity variation and study two cases concerning parsec and kilo-parsec scale jets by applying some simplified numerical simulations. For parsec scale jets, we use some general expressions and derive the shape of jets for certain numerical cases. For kilo-parsec scale jets, we study the cases of 4C53.37 and 1159+583 of Abell 2220 and Abell 1446 galaxy clusters, respectively, in an effort to test our model and compare our numerical results with observational data.