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.

Fast radio bursts (FRBs) are bright, short-duration radio transients with very high brightness temperatures implying highly coherent emission. We suggest that the FRBs are caused by the self-focusing of an electron beam interacting with an ambient plasma right beyond the light cylinder radius of a neutron star. The magnetic field at the light cylinder radius is relatively high which can accommodate both young Crab-like systems and old millisecond pulsars addressing the diverse environments of FRBs. At the first stage, the intense pulsed-beam passing through the background plasma causes instabilities such that the trapped particles in local Buneman-type cavitons saturate the local field. The beam is then radially self-focused due to the circular electric field developed by the two-stream instability which leads to Weibel instability in the transverse direction. Finally, the non-linear saturation of the Weibel instability results in the self-modulational formation of solitons due to plasmoid instability. The resonant solitary waves are the breather-type solitons hosting relativistic particles with self-excited oscillations. The analytical solutions obtained for non-linear dispersion and solitons suggest that, near the current sheets, the relativistic bunches are accelerated/amplified by klystron-like structures due to self-excited oscillations by the induced local electric field. Boosted coherent radio emission propagates through a narrow cone with strong focusing due to radial electric field and magnetic pinching. The non-linear evolution of solitons and the stimulated emission are associated with the Buneman instability and the possibility of the presence of nanosecond shots in FRBs are investigated.

Young isolated neutron stars (INS) most commonly manifest themselves as rotationally powered pulsars (RPPs) which involve conventional radio pulsars as well as gamma-ray pulsars (GRPs) and rotating radio transients (RRATs). Some other young INS families manifest themselves as anomalous X-ray pulsars (AXPs) and soft gamma-ray repeaters (SGRs) which are commonly accepted as magnetars, i.e. magnetically powered neutron stars with decaying superstrong fields. Yet some other young INS are identified as central compact objects (CCOs) and X-ray dim isolated neutron stars (XDINSs) which are cooling objects powered by their thermal energy. Older pulsars, as a result of a previous long episode of accretion from a companion, manifest themselves as millisecond pulsars and more commonly appear in binary systems. We use Dirichlet process Gaussian mixture model (DPGMM), an unsupervised machine learning algorithm, for analyzing the distribution of these pulsar families in the parameter space of period and period derivative. We compare the average values of the characteristic age, magnetic dipole field strength, surface temperature and transverse velocity of all discovered clusters. We verify that DPGMM is robust and provides hints for inferring relations between different classes of pulsars. We discuss the implications of our findings for the magneto-thermal spin evolution models and fallback discs.

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.

We study the structure of neutron stars in scalar-tensor theories for the nonminimal coupling of the form $(1+\kappa \xi \phi^{2})\cal R$. We solve the hydrostatic equilibrium equations for two different types of scalar field potentials and three different equations of state representative of different degrees of stiffness. We obtain the mass-radius relations of the configurations and determine the allowed ranges for the term $\xi\phi^2$ at the center of the star and spatial infinity based on the measured maximum value of the mass for neutron stars and the recent constraints on the radius coming from gravitational wave observations. Thus we manage to limit the deviation of the model from general relativity. We examine the possible constraints on the parameters of the model and compare the obtained restrictions with the ones inferred from other cosmological probes that give the allowed ranges for the coupling constant only. In the case of the Higgs-like potential, we also find that the central value for the scalar field cannot be chosen arbitrarily but it depends on the vacuum expectation value of the field. Finally, we discuss the effect of the scalar field potential on the mass and the radius of the star by comparing the results obtained for the cases considered here.

We study the evolution of newborn neutron stars in high-mass X-ray binaries interacting with a wind-fed super-Eddington disk. The inner disk is regularized to a radiation-dominated quasi-spherical configuration for which we calculate the inner radius of the disk, the total luminosity of the system and the torque acting on the neutron star accordingly, following the evolution of the system through the ejector and early propeller stages. We find that the systems with $B \gtrsim 10^{13}$ G pass through a short ($\sim 20\,{\rm yr}$) ejector stage appearing as supernova impostors followed by a propeller stage lasting $\sim 10^3\,{\rm yr}$. In the super-critical propeller stage the system is still bright ($L\sim 10^{40}\,{\rm erg\, s^{-1}}$) due to the spindown power and therefore appears as an ultra-luminous X-ray source (ULX). The system evolves into pulsating ULX (PULX) when the neutron star spins down to a period ($P\sim 1$ s) allowing for accretion onto its surface to commence. Systems with lower magnetic fields, $B \sim 10^{11}$ G, pass through a long ($10^5\,{\rm yr}$) super-critical propeller stage with luminosities similar to those of the ultra-luminous super-soft sources (ULS), $L \lesssim 10^{40}\,{\rm erg\, s^{-1}}$. The equilibrium periods of these systems in the accretion stage is about $10\,{\rm ms}$, which is much smaller than the typical period range of PULX observed to date. Such systems could have a larger population, but their pulsations would be elusive due to the smaller size of the magnetosphere. Our results suggest that the ULS and some nonpulsating ULX are rapidly spinning and highly magnetized young neutron stars at the super-critical propeller stage.

Ultra-luminous X-ray sources (ULXs) are off-nuclear point sources in nearby galaxies with luminosities well exceeding the Eddington limit for stellar-mass objects. It has been recognized after the discovery of pulsating ULXs (PULXs) that a fraction of these sources could be accreting neutron stars in high-mass X-ray binaries (HMXBs) though the majority of ULXs are lacking in coherent pulsations. The earliest stage of some HMXBs may harbor rapidly rotating neutron stars propelling out the matter transferred by the massive companion. The spin-down power transferred by the neutron-star magnetosphere to the accretion disk at this stage can well exceed the Eddington luminosities and the system appears as a non-pulsating ULX. In this picture, PULXs appear as super-critical mass-accreting descendants of non-pulsating ULXs. We present this evolutionary scenario within a self-consistent model of magnetosphere-disk interaction and discuss the implications of our results on the spin and magnetic field of the neutron star.

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}}$.

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.

Aql X--1 is one of the most prolific low mass X-ray binary transients (LMXBTs) showing outbursts almost annually. We present the results of our spectral analyses of RXTE/PCA observations of the 2000 and the 2011 outbursts. We investigate the spectral changes related to the changing disk-magnetosphere interaction modes of Aql X--1. The X-ray light curves of the outbursts of LMXBTs typically show phases of fast rise and exponential decay. The decay phase shows a "knee" where the flux goes from the slow decay to the rapid decay stage. We assume that the rapid decay corresponds to a weak propeller stage at which a fraction of the inflowing matter in the disk accretes onto the star. We introduce a novel method for inferring, from the light curve, the fraction of the inflowing matter in the disk that accretes onto the NS depending on the fastness parameter. We determine the fastness parameter range within which the transition from the accretion to the partial propeller stage is realized. This fastness parameter range is a measure of the scale-height of the disk in units of the inner disk radius. We applied the method to a sample of outbursts of Aql X--1 with different maximum flux and duration times. We show that different outbursts with different maximum luminosity and duration follow a similar path in the parameter space of accreted/inflowing mass flux fraction versus fastness parameter.

Using the spin and flux evolution of the accreting pulsar 4U 1626$-$67 across the 2008 torque reversal, we determine the fastness parameter dependence of the dimensionless torque acting on the pulsar. We find that the dimensionless torque is qualitatively different from the existing models: it is concave-up across the torque equilibrium whereas the existing torque models predict a concave-down (convex) relation with the fastness parameter. We show that the dimensionless torque has a cubic dependence on the fastness parameter near the torque equilibrium. We also find that the torque can not attain large values away from the equilibrium, either in the positive or the negative side, but saturates at limited values. The spin-down torque can attain a 2.5 times larger magnitude at the saturation limit than the spin-up torque. From the evolution of the frequency of quasi-periodic oscillations of 4U 1626$-$67 across the torque reversal of 1990, we determine the critical fastness parameter corresponding to torque equilibrium to be $\omega_{\rm c} \simeq 0.75$ within the framework of the beat frequency model and boundary region model for reasonable values of the model parameters. We find that the disc magnetosphere interaction becomes unstable when the inner radius approaches the corotation radius as predicted by some models, though with a longer timescale. We also find that there is an unstable regime that is triggered when the fastness parameter is 0.8 times the critical fastness parameter ($\omega=0.6$ for $\omega_{\rm c} \simeq 0.75$) possibly associated with an instability observed in numerical simulations.

The magnetic fields of the nascent neutron stars could be submerged in the crust by rapid fallback accretion and could diffuse to the surface later in life. According to this field burial scenario young pulsars may have growing magnetic fields, a process known to result with less-than-three braking indices; larger braking indices implying longer field-growth time-scales. A nascent neutron star with a larger kick velocity would accrete less amount of matter leading to a shallower burial of its field and a more rapid field growth. Such an inverse relation between the field growth time-scale inferred from the braking indices and space velocity of pulsars was claimed in the past as a prediction of the field-burial scenario. With a braking index of $n\sim 2$ and large space velocity PSR~B0540--69 was then an outlier in the claimed relation. The object recently made a transition to a rapid spin-down state accompanied by a low braking index. This new braking index implies a much shorter time-scale for the field growth which is consistent with the high space velocity of the object, in better agreement with the claimed relation. This observation lends support to the field burial scenario and implies that the growth of the magnetic field does not proceed at a constant pace but is slowed or completely halted at times. The slow spin-down state associated with the high braking index before 2011 which lasted for at least about 30 years was then such an episode of slow field growth.

We present the X-ray light curves of the last two outbursts --2014 & 2016-- of the well known accreting millisecond X-ray pulsar (AMXP) Aquila X-1 using the monitor of all sky X-ray image (MAXI) observations in the $2-20$ keV band. After calibrating the (MAXI) count rates to the all-sky monitor (ASM) level, we report that the 2016 outburst is the most energetic event of Aql X-1, ever observed from this source. We show that 2016 outburst is a member of the long-high class according to the classification presented by G\"ung\"or et al. with $\sim68$ cnt/s maximum flux and $\sim60$ days duration time and the previous outburst, 2014, belongs to the short-low class with $\sim 25$ cnt/s maximum flux and $\sim 30$ days duration time. In order to understand differences between outbursts, we investigate the possible dependence of the peak intensity to the quiescent duration leading to the outburst and find that the outbursts following longer quiescent episodes tend to reach higher peak energetic.

The maximum magnetic moment to angular momentum conjecture, recently posed by Barrow \& Gibbons, is in tension with the generally accepted parameters of magnetars. According to the conjecture the dimensionless Schuster-Wilson-Blacket number, $c\mu/J\sqrt{G}$ where $\mu$ is the magnetic moment and $J$ is the angular momentum, should be of order unity. A resolution is to assume that not only the low-magnetic field magnetars, but all magnetars have their super-critical fields in the higher multipoles. This requires the presence of external torques other than the magnetic dipole torque, such as from fallback disks or winds, on these systems to yield the measured high spin-down rates.

The glitch activity of young pulsars arises from the exchange of angular momentum between the crust and the interior of the star. Recently, it was inferred that the moment of inertia of the crust of a neutron star is not sufficient to explain the observed glitches. Such estimates are presumed in the Einstein's general relativity in describing the hydrostatic equilibrium of neutron stars. The crust of the neutron star has a space-time curvature of 14 orders of magnitude larger than that probed in solar system tests. This makes gravity the weakest constrained physics input in the crust related processes. We calculate the ratio of the crustal to the total moment of inertia of neutron stars in the scalar-tensor theory of gravity and the non-perturbative $f({\cal R})={\cal R}+ a {\cal R}^2$ gravity. We find for the former that the crust to core ratio of the moment of inertia does not change significantly from what is inferred in general relativity. For the latter we find that the ratio increases significantly from what is inferred in general relativity in the case of high mass objects. Our results suggest that the glitch activity of pulsars may be used to probe gravity models, although the gravity models explored in this work are not appropriate candidates.

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.

General properties of neutron stars are briefly reviewed with an emphasis on the indispensability of general relativity in our understanding of these fascinating objects. In Newtonian gravity the pressure within a star merely plays the role of opposing self-gravity. In general relativity all sources of energy and momentum contribute to the gravity. As a result the pressure not only opposes gravity but also enhances it. The later role of pressure becomes more pronounced with increasing compactness, $M/R$ where $M$ and $R$ are the mass and radius of the star, and sets a critical mass beyond which collapse is inevitable. This critical mass has no Newtonian analogue; it is conceptually different than the Stoner-Landau-Chandrasekhar limit in Newtonian gravity which is attained asymptotically for ultra-relativistic fermions. For white dwarfs the general relativistic critical mass is very close to the Stoner-Landau-Chandrasekhar limit. For neutron stars the maximum mass---so called Oppenheimer-Volkoff limit---is significantly smaller than the Stoner-Landau-Chandrasekhar limit. This follows from the fact that the general relativistic correction to hydrostatic equilibrium within a neutron star is significant throughout the star, including the central part where the mass contained within radial coordinate, $m(r)$, and the Newtonian gravitational acceleration, $Gm(r)/r^2$, are small.

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.

Fallback disks are expected to form around young neutron stars. The presence of these disks can be revealed by their blackbody spectrum in the infrared, optical and UV bands. We present a re-reduction of the archival optical and infrared data of the Vela pulsar, together with the existing infrared and UV spectrum of Vela, and model their unpulsed components with the blackbody spectrum of a supernova debris disk. We invoke the quiescent disk solution of Sunyaev and Shakura for the description of the disk in the propeller stage and find the inner radius of the disk to be inside the light cylinder radius. We perform a high resolution X-ray analysis with XMM-Newton and find a narrow absorption feature at 0.57 keV which can be interpreted as K$_{\alpha}$ line of He-like oxygen (OVII). The strength of the line indicates an element over-abundance in our line of sight exceeding the amounts that would be expected from ISM. The spectral feature may originate from the pulsar wind nebula and may be partly caused by the reprocessed X-ray radiation by the fallback disk. We discuss the lower-than-3 braking index of Vela as partially due to the contribution of the propeller torques. Our results suggest that the pulsar mechanism can work simultaneously with the propeller processes and that the debris disks can survive the radiation pressure for at least $\sim 10^4$ years. As Vela is a relatively close object, and a prototypical pulsar, the presence of a disk, if confirmed, may indicate the ubiquity of debris disks around young neutron stars.

Neutron stars (NSs) are thought to be excellent laboratories for determining the equation of state (EoS) of cold dense matter. Their strong gravity suggests that they can also be used to constrain gravity models. The mass and radius (M-R) of a NS both depend on the choice of EoS and gravity, meaning that NSs cannot be simultaneously good laboratories for both of these questions. A measurement of M-R would constrain the less well known physics input. The assumption that M-R measurements can be used to constrain EoS-presumes general relativity (GR) is the ultimate model of gravity in the classical regime. We calculate the radial profile of compactness and curvature (square root of the full contraction of the Weyl tensor) within a NS and determine the domain not probed by the Solar System tests of GR. We find that, except for a tiny sphere of radius less than a millimeter at the center, the curvature is several orders of magnitude above the values present in Solar System tests. The compactness is beyond the solar surface value for r>10 m, and increases by 5 orders of magnitude towards the surface. With the density being only an order of magnitude higher than that probed by nuclear scattering experiments, our results suggest that the employment of GR as the theory of gravity describing the hydrostatic equilibrium of NSs is a rather remarkable extrapolation from the regime of tested validity, as opposed to that of EoS models. Our larger ignorance of gravity within NSs suggests that a measurement of M-R constrains gravity rather than EoS, and given that EoS has yet to be determined by nucleon scattering experiments, M-R measurements cannot tightly constrain the gravity models either. Near the surface the curvature and compactness attain their largest values, while EoS in this region is fairly well known. This renders the crust as the best site to look for deviations from GR.

We present a broad classification of all outbursts detected with the All-Sky Monitor (ASM) on the Rossi X-Ray Timing Explorer (RXTE) and the Monitor of All Sky X-Ray Image (MAXI) of Aql X-1. We identify three types of outbursts; long-high, medium-low, and short-low, based on the duration and maximum flux. We analyse the trends in the "phase-space" of flux-derivative versus flux to demonstrate the differences in the three identified outburst types. We present a spectral analysis of the observations of Aql X-1 performed by the Proportional Counter Array (PCA) onboard RXTE during the 2000 and 2011 outbursts of the long-high class and the 2010 outburst of the medium-low class. We model the source spectrum with a hybrid thermal/non-thermal hot plasma emission model (EQPAIR in XSPEC, Coppi 2000) together with a Gaussian component to model the Fe K_alpha emission line. We construct time histories of the source flux, the optical depth of the corona (tau), the seed photon temperature (kT_bb) and the hard state compactness (l_h) for these three outbursts. We show that the physical parameters of either classes reach the same values throughout the outbursts, the only difference being the maximum flux. We discuss our results in the terms of modes of interaction of the star with the disc and size of the disc kept hot by irradiation. We conclude that irradiation is the dominant physical process leading to the different classes of outbursts.

It has long been accepted that a possible mechanism for explaining the existence of magnetic white dwarfs is the merger of a binary white dwarf system, as there are viable mechanisms for producing sustainable magnetic fields within the merger product. However, the lack of rapid rotators in the magnetic white dwarf population has been always considered a problematic issue of this scenario. Smoothed Particle Hydrodynamics simulations show that in mergers in which the two white dwarfs have different masses a disc around the central compact object is formed. If the central object is magnetized it can interact with the disc through its magnetosphere. The torque applied by the disc changes the spin of the star, whereas the transferred angular momentum from the star to the disc determines the properties of the disc. In this work we build a model for the disc evolution under the effect of magnetic accretion, and for the angular momentum evolution of the star, which can be compared with the observations. Our model predicts that the magnetospheric interaction of magnetic white dwarfs with their discs results in a significant spin down, and we show that for magnetic white dwarfs with relatively strong fields (larger than 10 MG) the observed rotation periods of the magnetic white dwarf population can be reproduced. We also investigate whether turbulence can be sustained during the late phases of the evolution of the system. When a critical temperature below which turbulence is not sustained is introduced into the model, the periods of the three fast rotating, strongly magnetic, massive white dwarfs in the solar neighborhood are recovered.

A nascent neutron star may be exposed to fallback accretion soon after the proto-neutron star stage. This high accretion episode can submerge the magnetic field deep in the crust. The diffusion of the magnetic field back to the surface will take hundreds to millions of years depending on the amount of mass accreted and the consequent depth the field is buried. Neutron stars with large kick velocities will accrete less amount of fallback material leading to shallower submergence of their fields and shorter time-scales for the growth of their fields. We obtain the relation $\tau_{\rm Ohm} \propto v^{-1}$ between the space velocity of the neutron star and Ohmic time-scale for the growth of the magnetic field. We compare this with the relation between the measured transverse velocities, $v_{\perp}$ and the field growth time-scales, $\mu/\dot{\mu}$, inferred from the measured braking indices. We find that the observational data is consistent with the theoretical prediction though the small number of data precludes a strong conclusion. Measurement of the transverse velocities of pulsars B1509$-$58, J1846$-$0258, J1119$-$6127 and J1734$-$3333 would increase the number of the data and strongly contribute to understanding whether pulsar fields grow following fallback accretion.

The viscous evolution of a thin disc around a central object is considered. Such discs are described by self-similar solutions in which either all or none of the inflowing mass accretes. An approximate solution for the partial accretion regime is constructed by employing a prescription recently introduced for nonlinear heat conduction. The solution is compared with numerical simulations demonstrating that the approximate solution describes the intermediate asymptotic stage for partially accreting discs.

It has long been accepted that a possible mechanism for explaining the existence of magnetic white dwarfs is the merger of a binary white dwarf system, as there are viable mechanisms for producing sustainable magnetism within the merger product. However, the lack of rapid rotators in the magnetic white dwarf population has been always considered a problematic issue of this scenario. In order to explain this discrepancy we build a model in which the interaction between the magnetosphere of the star and the disk induces angular momentum transfer. Our model predicts that the magnetospheric interaction of magnetic white dwarfs with their disks results in a significant spin down, and we show that the observed rotation period of REJ 0317-853, which is suggested to be a product of a double degenerate merger, can be reproduced.

We study the structure of neutron stars in R+\beta\ R^{\mu \nu} R_{\mu \nu} gravity model with perturbative method. We obtain mass--radius relations for six representative equations of state (EoSs). We find that, for |\beta| ~ 10^11 cm^2, the results differ substantially from the results of general relativity. Some of the soft EoSs that are excluded within the framework of general relativity can be reconciled for certain values of \beta\ of this order with the 2 solar mass neutron star recently observed. For values of \beta\ greater than a few 10^11 cm^2 we find a new solution branch allowing highly massive neutron stars. By referring some recent observational constraints on the mass--radius relation we try to constrain the value of \beta\ for each EoS. The associated length scale \sqrt{\beta} ~ 10^6 cm is of the order of the typical radius of neutron stars implying that this is the smallest value we could find by using neutron stars as a probe. We thus conclude that the true value of \beta\ is most likely much smaller than 10^11 cm^2.

We discuss the options of the fall-back disk model of Anomalous X-Ray Pulsars (AXPs). We argue that the power-law index of the mass inflow rate during the propeller stage can be lower than those employed in earlier models. We take into account the effect of the super-critical mass inflow at the earliest stages on the inner radius of the disk and argue that the system starts as a propeller. Our results show that, assuming a fraction of the mass inflow is accreted onto the neutron star, the fall-back disk scenario can produce AXPs for acceptable parameters.