Millisecond magnetars are often invoked as the central engine of some gamma-ray bursts (GRBs), specifically the ones showing a plateau phase. We argue that an apparent plateau phase may not be realized if the magnetic field of the nascent magnetar is in a transient rapid decay stage. Some GRBs that lack a clear plateau phase may also be hosting millisecond magnetars. We present an approximate analytical solution of the coupled set of equations describing the evolution of the angular velocity and the inclination angle between rotation and magnetic axis of a neutron star in the presence of a co-rotating plasma. We also show how the solution can be generalized to the case of evolving magnetic fields. We determine the evolution of the spin period, inclination angle, magnetic dipole moment and braking index of six putative magnetars associated with GRB 091018, GRB 070318, GRB 080430, GRB 090618, GRB 110715A, GRB 140206A through fitting, via Bayesian analysis, the X-ray afterglow light curves by using our recent model [\c{S}a\c{s}maz Mu\c{s} et al. 2019]. We find that within the first day following the formation of the millisecond magnetar, the inclination angle aligns rapidly, the magnetic dipole field may decay by a few times and the braking index varies by an order of magnitude.

The afterglow emission of some gamma-ray bursts (GRBs) show a shallow decay (plateau) phase implying continuous injection of energy. The source of this energy is very commonly attributed to the spin-down power of a nascent millisecond magnetar. The magnetic dipole radiation torque is considered to be the mechanism causing the spin-down of the neutron star. This torque has a component working for the alignment of the angle between rotation and magnetic axis, i.e., inclination angle, which has been neglected in modelling GRB afterglow light curves. Here, we demonstrate the evolution of the inclination angle and magnetic dipole moment of nascent magnetars associated with GRBs. We constrain the initial inclination angle, magnetic dipole moment and rotation period of seven magnetars by modelling the seven long-GRB afterglow light curves. We find that, in its first day, the inclination angle of a magnetar decreases rapidly. The rapid alignment of the magnetic and rotation axis may address the lack of persistent radio emission from mature magnetars. We also find that in three cases the magnetic dipole moments of magnetars decrease exponentially to a value a few times smaller than the initial value. The braking index of nascent magnetars, as a result of the alignment and magnetic dipole moment decline, is variable during the afterglow phase and always greater than three.

Deviations from the predictions of general relativity due to energy-momentum squared gravity (EMSG) are expected to become pronounced in the high density cores of neutron stars. We derive the hydrostatic equilibrium equations in EMSG and solve them numerically to obtain the neutron star mass-radius relations for four different realistic equations of state. We use the existing observational measurements of the masses and radii of neutron stars to constrain the free parameter, $\alpha ,$ that characterizes the coupling between matter and spacetime in EMSG. We show that $-10^{-38}\,\mathrm{cm^{3}/erg}<\alpha <+10^{-37}\,\mathrm{cm^{3}/erg}$. Under this constraint, we discuss what contributions EMSG can provide to the physics of neutron stars, in particular, their relevance to the so called \textit{hyperon puzzle} in neutron stars. We also discuss how EMSG alters the dynamics of the early universe from the predictions of the standard cosmological model. We show that EMSG leaves the standard cosmology safely unaltered back to $t\sim 10^{-4}$ seconds at which the energy density of the universe is $\sim 10^{34}\,\mathrm{erg\,cm^{-3}}$.

The vacuum solutions around a spherically symmetric and static object in the Starobinsky model are studied with a perturbative approach. The differential equations for the components of the metric and the Ricci scalar are obtained and solved by using the method of matched asymptotic expansions. The presence of higher order terms in this gravity model leads to the formation of a boundary layer near the surface of the star allowing the accommodation of the extra boundary conditions on the Ricci scalar. Accordingly, the metric can be different from the Schwarzschild solution near the star depending on the value of the Ricci scalar at the surface of the star while matching the Schwarzschild metric far from the star.

We study the structure of relativistic stars in $\mathcal{R}+\alpha \mathcal{R}^{2}$ theory using the method of matched asymptotic expansion to handle the higher order derivatives in field equations arising from the higher order curvature term. We find solutions, parametrized by $\alpha$, for uniform density stars. We obtain the mass-radius relations and study the dependence of maximum mass on $\alpha$. We find that $M_{\max}$ is almost linearly proportional to $\alpha$. For each $\alpha$ the maximum mass configuration has the biggest compactness parameter ($\eta = GM/Rc^2$), and we argue that the general relativistic stellar configuration corresponding to $\alpha=0$ is the least compact among these.

The recently discovered rotationally powered pulsar PSR J1640-4631 is the first to have a braking index measured, with high enough precision, that is greater than three. An inclined magnetic rotator in vacuum or plasma would be subject not only to spin-down but also to an alignment torque. The vacuum model can address the braking index only for an almost orthogonal rotator that is incompatible with the single peaked pulse profile. The magnetic dipole model with the corotating plasma predicts braking indices between $3-3.25$. We find that the braking index of $3.15$ is consistent with two different inclination angles, $18.5\pm 3$ degrees and $56 \pm 4$ degrees. The smaller angle is preferred given the pulse profile has a single peak and the radio output of the source is weak. We infer the change in the inclination angle to be at the rate $-0.23$ degrees per century, three times smaller in absolute value than the rate recently observed for the Crab pulsar.

The recent detection of pulsations from the ultra luminous X-ray source (ULX) NuSTAR J095551+6940.8 in M82 by Bachetti et al. indicates that the object is an accreting neutron star in a high mass X-ray binary (HMXB) system. The super-Eddington luminosity of the object implies that the magnetic field is sufficiently strong to suppress the scattering cross-section unless its beam is viewed at a favourable angle. We show that the torque equilibrium condition for the pulsar indicates the dipole magnetic field of the neutron star is $6.7 \times 10^{13}$ G, two orders of magnitude higher than that estimated by Bachetti et al., and further point to the possibility that even stronger magnetic fields could well be in the higher multipoles. This supports the recent view that magnetars descent from HMXBs if the magnetic field decays an order of magnitude during the process of transition.