Search

Abridged) We model the X-ray properties of millisecond pulsars (MSPs) by considering hot spot emission from a weakly magnetized rotating neutron star (NS) covered by an optically-thick hydrogen atmosphere. We investigate the limitations of using the thermal X-ray pulse profiles of MSPs to constrain the mass-to-radius ($M/R$) ratio of the underlying NS. The accuracy is strongly dependent on the viewing angle and magnetic inclination. For certain systems, the accuracy is ultimately limited only by photon statistics implying that future X-ray observatories could, in principle, achieve constraints on $M/R$ and hence the NS equation of state to better than $\sim$5%. We demonstrate that valuable information regarding the basic properties of the NS can be extracted even from X-ray data of fairly limited photon statistics through modeling of archival spectroscopic and timing observations of the nearby isolated PSRs J0030+0451 and J2124--3358. The X-ray emission from these pulsars is consistent with the presence of a hydrogen atmosphere and a dipolar magnetic field configuration, in agreement with previous findings for PSR J0437--4715. For both MSPs, the favorable geometry allows us to place interesting limits on the allowed $M/R$ of NSs. Assuming 1.4 M$_{\odot}$, the stellar radius is constrained to be $R > 9.4$ km and $R > 7.8$ km (68% confidence) for PSRs J0030+0451 and J2124--3358, respectively. We explore the prospects of using future observatories such as \textit{Constellation-X} and \textit{XEUS} to conduct blind X-ray timing searches for MSPs not detectable at radio wavelengths due to unfavorable viewing geometry. Using the observational constraints on the pulsar obliquities we are also able to place strong constraints on the magnetic field evolution model proposed by Ruderman., Comment: 9 pages, 7 figures, published in the Astrophysical Journal (Volume 689, Issue 1, pp. 407-415)

We present a model of thermal X-ray emission from hot spots on the surface of a rotating compact star with an unmagnetized light-element atmosphere. An application to ROSAT, Chandra, and XMM-Newton X-ray observations of the nearest known rotation-powered millisecond pulsar (MSP) PSR J0437--4715 reveals that the thermal emission from this pulsar is fully consistent with such a model, enabling constraints on important properties of the underlying neutron star. We confirm that the observed thermal X-ray pulsations from J0437--4715 are incompatible with blackbody emission and require the presence of an optically thick, light element (most likely hydrogen) atmosphere on the neutron star surface. The morphology of the X-ray pulse profile is consistent with a global dipole configuration of the pulsar magnetic field but suggests an off-center magnetic axis, with a displacement of 0.8-3 km from the stellar center. For an assumed mass of 1.4 M$_{\odot}$, the model restricts the allowed stellar radii to R=6.8-13.8 km (90% confidence) and R>6.7 km (99.9% confidence), which is consistent with standard NS equations of state and rules out an ultracompact star smaller than its photon sphere. Deeper spectroscopic and timing observations of this and other nearby radio MSPs with current and future X-ray facilities (Constellation-X and XEUS) can provide further insight into the fundamental properties of neutron stars., Comment: 10 pages, 8 figures; accepted for publication in the Astrophysical Journal

The majority of X-ray-detected rotation-powered millisecond pulsars (MSPs) appear to exhibit predominantly thermal emission, believed to originate from the heated magnetic polar caps of the pulsar. In the nearest MSP, J0437--4715 a faint PL is also observed at >3 keV, usually associated with magnetospheric emission processes. However, the hard emission in this and other similar MSPs may instead be due to weak Comptonization of the thermal polar cap emission by energetic electrons/positrons of small optical depth most likely in the pulsar magnetosphere. This spectral model implies that all soft X-rays are of purely thermal origin, which has important implications in the study of neutron stars., Comment: Proceedings of the 363. WE-Heraeus Seminar on: Neutron Stars and Pulsars (Posters and contributed talks) Physikzentrum Bad Honnef, Germany, May.14-19, 2006, eds. W.Becker, H.H.Huang, MPE Report 291, pp.5-8

X-ray emission from many rotation-powered millisecond pulsars (MSPs) is observed to be of predominantly thermal nature. In PSR J0437--4715, the nearest MSP known, an additional faint power-law tail is observed above 2.5 keV, commonly attributed to non-thermal magnetospheric radiation. We propose that the hard emission in this and other similar MSPs is instead due to weak Comptonization of the thermal (blackbody or hydrogen atmosphere) polar cap emission by energetic electrons/positrons of small optical depth in the pulsar magnetosphere. This spectral model implies that all soft X-rays are of purely thermal origin, which has profound implications in the study of neutron star structure and fundamental pulsar physics., Comment: 4 pages, 3 figures, submitted for publication in the Astrophysical Journal Letters

A new Fokker-Planck equation is developed for treating resonance line scattering, especially relevant to the treatment of Lyman alpha in the early universe. It is a "corrected" form of the equation of Rybicki & Dell'Antonio that now obeys detailed balance, so that the approach to thermal equilibrium is properly described. The new equation takes into account the energy changes due to scattering off moving particles, the recoil term of Basko, and stimulated scattering. One result is a surprising unification of the equation for resonance line scattering and the Kompaneets equation. An improved energy exchange formula due to resonance line scattering is derived. This formula is compared to previous formulas of Madau, Meikson, & Rees (1997) and Chen & Miralda-Escud\'e (2004)., Comment: 10 pages, no figures, submitted to ApJ

Current X-ray missions are providing high-quality X-ray spectra from neutron stars (NSs) in quiescent low-mass X-ray binaries (qLMXBs). This has motivated us to calculate new hydrogen-atmosphere models, including opacity due to free-free absorption and Thomson scattering, thermal electron conduction, and self-irradiation by photons from the compact object. We have constructed a self-consistent grid of neutron star models covering a wide range of surface gravities as well as effective temperatures, which we make available to the scientific community. We present multi-epoch Chandra X-ray observations of the qLMXB X7 in the globular cluster 47 Tuc, which is remarkably nonvariable on timescales from minutes to years. Its high-quality X-ray spectrum is adequately fit by our hydrogen-atmosphere model without any hard power-law component or narrow spectral features. If a mass of 1.4 Msol is assumed, our spectral fits require that its radius be in the range R=14.5^{+1.8}_{-1.6} km (90% confidence), larger than expected from currently preferred models of NS interiors. If its radius is assumed to be 10 km, then a mass of M=2.20^{+0.03}_{-0.16} Msol is required. Using models with the appropriate surface gravity for each value of the mass and radius becomes important for interpretation of the highest quality data., Comment: 14 pages, 6 figures, ApJ in press (replaced with accepted version)

We develop a method of analyzing radio frequency spectral line observations to derive data on the temperature, density, velocity, and molecular abundance of the emitting gas. The method incorporates a radiative transfer code with a new technique for handling overlapping hyperfine emission lines within the accelerated lambda iteration algorithm and a heuristic search algorithm based on simulated annnealing. We apply this method to new observations of N_2H^+ in three Lynds clouds thought to be starless cores in the first stages of star formation and determine their density structure. A comparison of the gas densities derived from the molecular line emission and the millimeter dust emission suggests that the required dust mass opacity is about kappa_{1.3mm}=0.04 cm^2/g, consistent with models of dust grains that have opacities enhanced by ice mantles and fluffy aggregrates., Comment: 42 pages, 17 figures, to appear in ApJ

A soft component of thermal emission is very commonly observed from the surfaces of quiescent, accreting neutron stars. We searched with Chandra for such a surface component of emission from the dynamical black-hole candidate XTE J1118+480 (= J1118), which has a primary mass of about 8 solar masses. None was found, as one would expect if the compact X-ray source is a bona fide black hole that possesses an event horizon. The X-ray spectrum of J1118 is well-fitted by a simple power-law model. In our search for a thermal component, we fitted our Chandra data to a power-law model (with slope and column density fixed) plus a series of nine hydrogen-atmosphere models with radii ranging from 9/8 to 2.8 Schwarzschild radii. For the more compact models, we included the important effect of self-irradiation of the atmosphere. Because of the extremely low column density to J1118, we obtained very strong limits on a hypothetical thermal source: kT(infinity) < 0.011 keV and L(infinity) < 9.4E30 erg/s (99% confidence level). In analogy with neutron stars, we examine two possible sources of thermal radiation from a hypothetical surface of J1118: deep crustal heating and accretion in quiescence. We conclude that if J1118 possesses a surface composed of normal nuclear matter, we would have detected thermal radiation from this surface. The absence of thermal radiation indicates that J1118 is a true black hole with an event horizon., Comment: 41 pages, 8 figures; ApJ, in press; 20% longer including a new section on exotic alternatives to black holes, a new figure, and further discussion of SAX J1808.4-3658

A kinetic equation for Compton scattering is given that differs from the Kompaneets equation in several significant ways. By using an inverse differential operator this equation allows treatment of problems for which the radiation field varies rapidly on the scale of the width of the Compton kernel. This inverse operator method describes, among other effects, the thermal Doppler broadening of spectral lines and continuum edges, and automatically incorporates the process of Compton heating/cooling. It is well adapted for inclusion into a numerical iterative solution of radiative transfer problems. The equivalent kernel of the new method is shown to be a positive function and with reasonable accuracy near the intitial frequency, unlike the Kompaneets kernel, which is singular and not wholly positive. It is shown that iterates of the inverse operator kernel can be easily calculated numerically, and a simple summation formula over these iterates is derived that can be efficiently used to compute Comptonized spectra. It is shown that the new method can be used for initial value and other problems with no more numerical effort than the Kompaneets equation, and that it more correctly describes the solution over times comparable to the mean scattering time., Comment: 27 pages, 5 figures, to be published in ApJ. Minor changes, including one reference corrected

In Loeb & Rybicki (1999; paper I) it was shown that before reionization, the scattering of \lya photons from a cosmological source might lead to a fairly compact ($\sim 15''$) \lya halo around the source. Observations of such halos could constrain the properties of the neutral intergalactic medium (IGM), and in particular yield the cosmological density parameters of baryons and matter on scales where the Hubble flow is unperturbed. Paper I did not treat the polarization of this scattered radiation, but did suggest that the degree of such polarization might be large. In this Letter we report on improved calculations for these \lya halos, now accounting for the polarization of the radiation field. The polarization is linear and is oriented tangentially to the projected displacement from the center of the source. The degree of polarization is found to be 14% at the core radius, where the intensity has fallen to half of the central value. It rises to 32% and 45% at the radii where the intensity has fallen to one-tenth and one-hundreth of the central intensity, respectively. At larger radii the degree of polarization rises further, asymptotically to 60%. Such high values of polarization should be easily observable and provide a clear signature of the phenomenon of \lya halos surrounding sources prior to reionization., Comment: 8 pages, 2 Postscript figures, accepted by Astrophysical Journal Letters; some typos corrected; added two paragraphs at the end of section 3 concerning detectability of Lyman alpha halos

Radio observations establish the B/A magnification ratio of gravitational lens 0957+561 at about 0.75. Yet, for more than 15 years, the optical magnfication ratio has been between 0.9 and 1.12. The accepted explanation is microlensing of the optical source. However, this explanation is mildly discordant with (i) the relative constancy of the optical ratio, and (ii) recent data indicating possible non-achromaticity in the ratio. To study these issues, we develop a statistical formalism for separately measuring, in a unified manner, the magnification ratio of the fluctuating and constant parts of the light curve. Applying the formalism to the published data of Kundi\'c et al. (1997), we find that the magnification ratios of fluctuating parts in both the g and r colors agrees with the magnification ratio of the constant part in g-band, and tends to disagree with the r-band value. One explanation could be about 0.1 mag of consistently unsubtracted r light from the lensing galaxy G1, which seems unlikely. Another could be that 0957+561 is approaching a caustic in the microlensing pattern., Comment: 12 pages including 1 PostScript figure

The observed profile of spectral lines from a relativistic accretion disk can constrain parameters such as the disk geometry and the rotation of the central black hole. The formation of the spectral line in a disk generally has been modeled with simple assumptions such as local isotropy of emission. Here we consider line formation in the presence of velocity gradients induced by the differential flow in the disk. In this case the emission can have anisotropy in the form of an azimuthal dependence relative to the local principle axes of shear. Since the physical conditions in a disk are uncertain in detail, we investigate this effect with simple parameterized models based on Sobolev theory to highlight the overall character of the changes in the line profile. We find that velocity gradients generally cause a relative increase of flux in the red wing, hence the inner radius of the disk would be underestimated if the effect were not taken into consideration. If the inner radius is used as a signature of black hole rotation, as when the disk is not emissive within the marginally stable circular orbit, then the inferred rotation would be overestimated in cases where the emissivity of the disk has fairly shallow fall-off with radius. If the disk were emissive even within the marginally stable orbit, then the local azimuthal anisotropy of emission produces features in the line profile which distinguish rotating from nonrotating black holes., Comment: 16 pages, LaTeX, including 5 Postscript figures

With the ansatz that a data set's correlation matrix has a certain parametrized form (one general enough, however, to allow the arbitrary specification of a slowly-varying decorrelation distance and population variance) the general machinery of Wiener or optimal filtering can be reduced from $O(n^3)$ to $O(n)$ operations, where $n$ is the size of the data set. The implied vast increases in computational speed can allow many common sub-optimal or heuristic data analysis methods to be replaced by fast, relatively sophisticated, statistical algorithms. Three examples are given: data rectification, high- or low- pass filtering, and linear least squares fitting to a model with unaligned data points., Comment: 7 pages, LaTeX with REVTeX 3.0 macros, no figures. A toolkit with implementations (in Fortran 90) of the algorithms is available by anonymous ftp to cfata4.harvard.edu

The time-dependent spectral profile of a resonance line in a homogeneous expanding medium is studied by numerically solving an improved Fokker-Planck diffusion equation. The solutions are used to determine the time required to reach a quasi-static solution near the line center. A simple scaling law for this relaxation time is derived and is fitted to the numerical results. The results are applied to the case of Lyman alpha scattering during primordial recombination of hydrogen. For a wide range of cosmological models it is found that the relaxation times are smaller than the recombination timescale, although not by a very large factor. Thus the standard assumption of a quasi-static solution in cosmological recombination calculations is reasonably valid, and should not cause substantial errors in the solutions., Comment: 20 pages text and 10 figures, in 30 pages of uuencoded, compressed postscript. CFA preprint no. 3753

Curve of growth analysis, applied to the Lyman series absorption ratios deduced in our previous paper, yields a measurement of the logarithmic slope of distribution of \Lya\ clouds in column density $N$. The observed exponential distribution of the clouds' equivalent widths $W$ is then shown to require a broad distribution of velocity parameters $b$, extending up to 80 km s$^{-1}$. We show how the exponential itself emerges in a natural way. An absolute normalization for the differential distribution of cloud numbers in $z$, $N$, and $b$ is obtained. By detailed analysis of absorption fluctuations along the line of sight we are able to put upper limits on the cloud-cloud correlation function $\xi$ on several megaparsec length scales. We show that observed $b$ values, if thermal, are incompatible, in several different ways, with the hypothesis of equilibrium heating and ionization by a background UV flux. Either a significant component of $b$ is due to bulk motion (which we argue against on several grounds), or else the clouds are out of equilibrium, and hotter than is implied by their ionization state, a situation which could be indicative of recent adiabatic collapse., Comment: 32 pages, LaTeX using aastex30 macros, submitted to Ap.J

Techniques for the statistical analysis of the \Lya\ forest in high redshift quasars are developed, and applied to the low resolution (25 \AA) spectra of 29 of the 33 quasars in the Schneider-Schmidt-Gunn (SSG) sample.We find that the mean absorption increases with $z$ approximately as a power law $(1+z)^{\gamma+1}$ with $\gamma = 2.46\pm 0.37$. The mean ratio of \Lya\ to Lyman $\beta$ absorption in the clouds is $0.476\pm 0.054$. We also detect, and obtain ratios, for Lyman $\beta$, $\gamma$, and possibly $\epsilon$. We are also able to quantify the fluctuations of the absorption around its mean, and find that these are comparable to, or perhaps slightly larger than, that expected from an uncorrelated distribution of clouds. The techniques in this paper, which include the use of bootstrap resampling of the quasar sample to obtain estimated errors and error covariances, and a mathematical treatment of absorption from a (possibly non-uniform) stochastic distribution of lines, should be applicable to future, more extensive, data sets., Comment: 29 pages, LaTeX using aastex30 macros, forthcoming as CfA preprint

Earlier analyses of the linear instability of line-driven stellar winds are extended to the case, relevant to Wolf-Rayet stars, in which the continuum remains optically thick well above the sonic point. It is found that an optically thick flow driven by pure scattering lines is stabilized by the drag effect of the diffuse, scattered radiation. However, even a relatively small photon destruction probability can cause a flow with continuum optical thickness much greater than 1 to remain unstable, with a given growth rate. The implications of these results for the variability characteristics of winds from Wolf-Rayet stars are briefly discussed.

Numerical radiation-hydrodynamics simulations of the nonlinear evolution of instabilities in radiatively driven stellar winds have been performed. The results show a strong tendency for the unstable flow to form rather sharp rarefactions in which the highest speed material has very low density. The qualitative features of the model agree well with the reqirements of displaced narrow absorption components in UV lines.