We present an analytical unified representation for 22 equations of state (EoS) of dense matter in neutron stars. Such analytical representations can be useful for modeling neutron star structure in modified theories of gravity with high order derivatives.

We study the structure of neutron stars in perturbative f(R) gravity models with realistic equations of state. We obtain mass-radius relations in a gravity model of the form f(R)=R+\alpha R^2. We find that deviations from the results of general relativity, comparable to the variations due to using different equations of state (EoS'), are induced for |alpha| ~ 10^9 cm^2. Some of the soft EoS' that are excluded within the framework of general relativity can be reconciled with the 2 solar mass neutron star recently observed for certain values of alpha within this range. For some of the EoS' we find that a new solution branch, which allows highly massive neutron stars, exists for values of alpha greater than a few 10^9 cm^2. We find constraints on alpha for a variety of EoS' using the recent observational constraints on the mass-radius relation. These are all 5 orders of magnitude smaller than the recent constraint obtained via Gravity Probe B for this gravity model. The associated length scale \sqrt{alpha} ~ 10^5 cm is only an order of magnitude smaller than the typical radius of a neutron star, the probe used in this test. This implies that real deviations from general relativity can be even smaller.

Accretion driven millisecond X-ray pulsars can accrete over a wide range of mass flow rates. The pulsations persist even at small accretion rates at which these objects would be expected to be in the propeller stage. We argue that the inner regions of disks around millisecond X-ray pulsars are sufficiently thick that a fraction of the inflowing matter can accrete even in the spin-down regime from regions of the disk away from the disk plane. This allows these systems to be bright throughout a wide range of mass flow rates. We model the lightcurve of SAX J1808.4-3658 during an outburst and show that the rapid decay stage can be modeled with fractional accretion in the spin-down regime.

We show that the period clustering of anomalous X-ray pulsars (AXPs) and soft gamma-ray repeaters (SGRs), their X-ray luminosities, ages and statistics can be explained with fallback disks with large initial specific angular momentum. The disk evolution models are developed by comparison to self-similar analytical models. The initial disk mass and angular momentum set the viscous timescale. An efficient torque, with (1 - w^2) dependence on the fastness parameter w leads to period clustering in the observed AXP-SGR period range under a wide range of initial conditions. The timescale t_0 for the early evolution of the fallback disk, and the final stages of fallback disk evolution, when the disk becomes passive, are the crucial determinants of the evolution. The disk becomes passive at temperatures around 100 K, which provides a natural cutoff for the X-ray luminosity and defines the end of evolution in the observable AXP and SGR phase. This low value for the minimum temperature for active disk turbulence indicates that the fallback disks are active up to a large radius greater than ~10^{12} cm. We find that transient AXPs and SGRs are likely to be older than their persistent cousins. A fallback disk with mass transfer rates corresponding to the low quiescent X-ray luminosities of the transient sources in early evolutionary phases would have a relatively lower initial mass, such that the mass-flow rate in the disk is not sufficient for the inner disk to penetrate into the light cylinder of the young neutron star, making mass accretion onto the neutron star impossible. The transient AXP phase therefore must start later. The model results imply that the transient AXP/SGRs, although older, are likely to be similar in number to persistent sources (abridged).

The recent detection of the anomalous X-ray pulsar (AXP) 4U 0142+61 in the mid-infrared with the {\it Spitzer} observatory by Z.Wang and coworkers constitutes the first instance of a disk around an AXP. We show, by analyzing earlier optical and near-IR data together with the recent data, that the overall broadband data set can be reproduced by a single model of an irradiated and viscously heated disk.

The anomalous X-ray pulsar 4U 0142+61 was recently detected in the mid infrared bands with the SPITZER Observatory (Wang, Chakrabarty & Kaplan 2006). This observation is the first instance for a disk around an AXP. From a reanalysis of optical and infrared data, we show that the observations indicate that the disk is likely to be an active disk rather than a passive dust disk beyond the light cylinder, as proposed in the discovering paper. Furthermore, we show that the irradiated accretion disk model can also account for all the optical and infrared observations of the anomalous X-ray pulsars in the persistent state.

Fallback disks are expected to form around new-born neutron stars following a supernova explosion. In almost all cases, the disk will pass through a propeller stage. If the neutron star is spinning rapidly (initial period $\sim 10$ ms) and has an ordinary magnetic moment ($\sim 10^{30}$ G cm$^3$), the rotational power transferred to the disk by the magnetic field of the neutron star will exceed the Eddington limit by many orders of magnitude, and the disk will be rapidly disrupted. Fallback disks can thus survive only around slow-born neutron stars and around black holes, assuming the latter do not torque their surrounding disks as strongly as do neutron stars. This might explain the apparent rarity of fallback disks around young compact objects.

The radiation pressure of a radio pulsar does not necessarily disrupt a surrounding disk. The position of the inner radius of a thin disk around a neutron star can be estimated by comparing the electromagnetic energy density generated by the neutron star with the kinetic energy density of the disk. Inside the light cylinder, the near zone electromagnetic field is essentially the dipole magnetic field, and the inner radius is the conventional Alfven radius. Far outside the light cylinder, in the radiation zone, $E=B$ and the electromagnetic energy density is $<S >/c \propto 1/r^2$ where $S$ is the Poynting vector. Shvartsman (1970) argued that a stable equilibrium can not be found in the radiative zone because the electromagnetic energy density dominates over the kinetic energy density, with the relative strength of the electromagnetic stresses increasing with radius. In order to check whether this is true also near the light cylinder, we employ global electromagnetic field solutions for rotating oblique magnetic dipoles (Deutsch 1955). Near the light cylinder the electromagnetic energy density increases steeply enough with decreasing $r$ to balance the kinetic energy density at a stable equilibrium. The transition from the near zone to the radiation zone is broad. The radiation pressure of the pulsar can not disrupt the disk for values of the inner radius up to about twice the light cylinder radius if the rotation axis and the magnetic axis are orthogonal. This allowed range beyond the light cylinder extends much further for small inclination angles. We discuss implications of this result for accretion driven millisecond pulsars and young neutron stars with fallback disks.

We investigate whether young neutron stars with fall-back disks can produce Anomalous X-Ray Pulsars (AXPs) within timescales indicated by the ages of associated supernova remnants. The system passes through a propeller stage before emerging as an AXP or a radio pulsar. The evolution of the disk is described by a diffusion equation which has self-similar solutions with either angular momentum or total mass of the disk conserved. We associate these two types of solutions with accretor and propeller regimes, respectively. Our numerical calculations of thin disk models with changing inner radius take into account the super-critical accretion at the early stages, and electron scattering and bound-free opacities with rich metal content. 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.

In this paper, we further develop the model for the migration of planets introduced in Del Popolo et al. (2001). We first model the protoplanetary nebula as a time-dependent accretion disc and find self-similar solutions to the equations of the accretion disc that give to us explicit formulas for the spatial structure and the temporal evolution of the nebula. These equations are then used to obtain the migration rate of the planet in the planetesimal disc and to study how the migration rate depends on the disc mass, on its time evolution and on some values of the dimensionless viscosity parameter alpha. We find that planets that are embedded in planetesimal discs, having total mass of 10^{-4}-0.1 M_{\odot}, can migrate inward a large distance for low values of alpha (e.g., alpha \simeq 10^{-3}-10^{-2}) and/or large disc mass and can survive only if the inner disc is truncated or because of tidal interaction with the star. Orbits with larger $a$ are obtained for smaller value of the disc mass and/or for larger values of alpha. This model may explain several orbital features of the recently discovered giant planets orbiting nearby stars.