In this paper we discuss two approximate methods previously suggested for modeling hyperfine spectral line emission for molecules whose collisional transitions rates between hyperfine levels are unknown. Hyperfine structure is seen in the rotational spectra of many commonly observed molecules such as HCN, HNC, NH3, N2H+, and C17O. The intensities of these spectral lines can be modeled by numerical techniques such as Lambda-iteration that alternately solve the equations of statistical equilibrium and the equation of radiative transfer. However, these calculations require knowledge of both the radiative and collisional rates for all transitions. For most commonly observed radio frequency spectral lines, only the net collisional rates between rotational levels are known. For such cases, two approximate methods have been suggested. The first method, hyperfine statistical equilibrium (HSE), distributes the hyperfine level populations according to their statistical weight, but allows the population of the rotational states to depart from local thermodynamic equilibrium (LTE). The second method, the proportional method approximates the collision rates between the hyperfine levels as fractions of the net rotational rate apportioned according to the statistical degeneracy of the final hyperfine levels. The second method is able to model non-LTE hyperfine emission. We compare simulations of N2H+ hyperfine lines made with approximate and more exact rates and find that satisfactory results are obtained.

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.

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.

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.

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.

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).

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.

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.

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.

A large number of neutron stars (NSs), ~10^9, populate the Galaxy, but only a tiny fraction of them is observable during the short radio pulsar lifetime. The majority of these isolated NSs, too cold to be detectable by their own thermal emission, should be visible in X-rays as a result of accretion from the interstellar medium. The ROSAT all sky survey has however shown that such accreting isolated NSs are very elusive: only a few tentative candidates have been identified, contrary to theoretical predictions that up to several thousands should be seen. We suggest that the fundamental reason for this discrepancy lies in the use of the standard Bondi formula to estimate the accretion rates. We compute the expected source counts using updated estimates of the pulsar velocity distribution, realistic hydrogen atmosphere spectra, and a modified expression for the Bondi accretion rate as suggested by recent MHD simulations, and supported by direct observations in the case of accretion around supermassive black holes in nearby galaxies and in our own. We find that, whereas the inclusion of atmospheric spectra partly compensates for the reduction in the counts due to the higher mean velocities of the new distribution, the modified Bondi formula dramatically suppresses the source counts. The new predictions are consistent with a null detection at the ROSAT sensitivity.

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.

We investigate the structure of self-gravitating polytropic stellar systems. We present a method which allows to obtain approximate analytical solutions, $\psi_{n+\epsilon}({\bf x})$, of the nonlinear Poisson equation with the polytropic index $n+\epsilon$, given the solution $\psi_n({\bf x})$ with the polytropic index n, for any positive or negative $\epsilon$ such that $|\epsilon|\ll1$. Application of this method to the spherically symmetric stellar polytropes with $n\simeq5$ yields the solutions which describe spatially bound systems if n<5 and the formation of a second core if n>5. A heuristic approximate expressions for the radial profiles are also presented. Due to the duality between stellar and gas polytropes, our results are valid for gaseous, self-gravitating polytropic systems (e.g., molecular clouds) with the index $\gamma\simeq 6/5$. Stability of such systems and observational consequences for both stellar and gaseous systems are discussed.

The spectra of the first galaxies and quasars in the Universe should be strongly absorbed shortward of their rest-frame Lyman-alpha wavelength by neutral hydrogen (HI) in the intervening intergalactic medium. However, the Lyman-alpha line photons emitted by these sources are not eliminated but rather scatter until they redshift out of resonance and escape due to the Hubble expansion of the surrounding intergalactic HI. We calculate the resulting brightness distribution and the spectral shape of the diffuse Lyman-alpha line emission around high redshift sources, before the intergalactic medium was reionized. Typically, the Lyman-alpha photons emitted by a source at z=10 scatter over a characteristic angular radius of order 15 arcseconds around the source and compose a line which is broadened and redshifted by about a thousand km/s relative to the source. The scattered photons are highly polarized. Detection of the diffuse Lyman-alpha halos around high redshift sources would provide a unique tool for probing the neutral intergalactic medium before the epoch of reionization. On sufficiently large scales where the Hubble flow is smooth and the gas is neutral, the Lyman-alpha brightness distribution can be used to determine the cosmological mass densities of baryons and matter.

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.

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.

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.

The deprojection of axisymmetric density distributions is generally indeterminate to within the addition of certain axisymmetric distributions (konus densities) that are invisible in projection. The known class of konus densities is expanded considerably here through the introduction of {\em semikonus functions}. These functions are closed with respect to multiplication in ordinary space, and the real part of an arbitrary polynomial of semikonus functions is a konus function. This property facilitates the construction of semikonus (and konus) functions with tailored properties, such as asymptotic forms. We also develop a simple technique for constructing several classes of konus distributions with arbitrary density profiles in the equatorial plane.

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.

Since the mass of the electron is very small relative to atomic masses, Thomson scattering of low-energy photons ($h\nu \ll m_ec^2$) produces thermal Doppler frequency shifts that are much larger than atomic Doppler widths. A method is developed here to evaluate the electron scattering emissivity from a given radiation field which is considerably faster than previous methods based on straightforward evaluation of the scattering integral. This procedure is implemented in our multilevel radiative code (MALI), which now takes full account of the effects of noncoherent electron scattering on level populations, as well as on the emergent spectrum. Calculations using model atmospheres of hot, low-gravity stars display not only the expected broad wings of strong emission lines but also effects arising from the scattering of photons across continuum edges. In extreme cases this leads to significant shifts of the ionization equilibrium of helium.

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.

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.

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.