The most fundamental problem in statistics is the inference of an unknown probability distribution from a finite number of samples. For a specific observed data set, answers to the following questions would be desirable: (1) Estimation: Which candidate distribution provides the best fit to the observed data?, (2) Goodness-of-fit: How concordant is this distribution with the observed data?, and (3) Uncertainty: How concordant are other candidate distributions with the observed data? A simple unified approach for univariate data that addresses these traditionally distinct statistical notions is presented called "maximum fidelity". Maximum fidelity is a strict frequentist approach that is fundamentally based on model concordance with the observed data. The fidelity statistic is a general information measure based on the coordinate-independent cumulative distribution and critical yet previously neglected symmetry considerations. An approximation for the null distribution of the fidelity allows its direct conversion to absolute model concordance (p value). Fidelity maximization allows identification of the most concordant model distribution, generating a method for parameter estimation, with neighboring, less concordant distributions providing the "uncertainty" in this estimate. Maximum fidelity provides an optimal approach for parameter estimation (superior to maximum likelihood) and a generally optimal approach for goodness-of-fit assessment of arbitrary models applied to univariate data. Extensions to binary data, binned data, multidimensional data, and classical parametric and nonparametric statistical tests are described. Maximum fidelity provides a philosophically consistent, robust, and seemingly optimal foundation for statistical inference. All findings are presented in an elementary way to be immediately accessible to all researchers utilizing statistical analysis.

XMM-Newton RGS spectra of MCG-6-30-15 and Mrk 766 exhibit complex discrete structure, which was interpreted in a paper by Branduardi-Raymont et al. (2001) as evidence for the existence of relativistically broadened Lyman alpha emission from carbon, nitrogen, and oxygen, produced in the inner-most regions of an accretion disk around a Kerr black hole. This suggestion was subsequently criticized in a paper by Lee et al. (2001), who argued that for MCG-6-30-15, the Chandra HETG spectrum, which is partially overlapping the RGS in spectral coverage, is adequately fit by a dusty warm absorber model, with no relativistic line emission. We present a reanalysis of the original RGS data sets in terms of the Lee et al. (2001) model, and demonstrate that spectral models consisting of a smooth continuum with ionized and dust absorption alone cannot reproduce the RGS spectra of both objects. The original relativistic line model with warm absorption proposed by Branduardi-Raymont et al. (2001) provides a superior fit to the RGS data, both in the overall shape of the spectrum and in the discrete absorption lines. Limits on the amount of X-ray absorption by dust particles are discussed. We also discuss a possible theoretical interpretation for the putative relativistic Lyman alpha line emission in terms of the photoionized surface layers of the inner regions of an accretion disk.

We present a detailed model of the discrete X-ray spectroscopic features expected from steady-state, low-density photoionized plasmas. We apply the Flexible Atomic Code (FAC) to calculate all of the necessary atomic data for the full range of ions relevant for the X-ray regime. These calculations have been incorporated into a simple model of a cone of ions irradiated by a point source located at its tip (now available as the XSPEC model PHOTOION). For each ionic species in the cone, photoionization is balanced by recombination and ensuing radiative cascades, and photoexcitation of resonance transitions is balanced by radiative decay. This simple model is useful for diagnosing X-ray emission mechanisms, determining photoionization/photoexcitation/recombination rates, fitting temperatures and ionic emission measures, and probing geometrical properties (covering factor/column densities/radial filling factor/velocity distributions) of absorbing/reemitting regions in photoionized plasmas. Such plasmas have already been observed in diverse astrophysical X-ray sources, including active galactic nuclei, X-ray binaries, cataclysmic variables, and stellar winds of early-type stars, and may also provide a significant contribution to the X-ray spectra of gamma-ray-burst afterglows and the intergalactic medium.

We present Chandra HETG spectra of seven cataclysmic variables. We find that they divide unambiguously into two distinct types. Spectra of the first type are remarkably well fit by a simple cooling flow model, which assumes only steady-state isobaric radiative cooling. The maximum temperature and the normalization, which provides a highly precise measurement of the accretion rate, are the only free parameters of this model. Spectra of the second type are grossly inconsistent with a cooling flow model. They instead exhibit a hard continuum, and show strong H-like and He-like ion emission but little Fe L-shell emission, which is consistent with expectations for line emission from a photoionized plasma. Using a simple photoionization model, we argue that the observed line emission for these sources can be driven entirely by the hard continuum. The physical significance of these two distinct types of X-ray spectra is also explored.

Using the combined spectral and spatial resolving power of the Low Energy Transmission Grating (LETGS) on board Chandra, we obtain separate spectra from the bright central source of NGC 1068 (Primary region), and from a fainter bright spot 4" to the NE (Secondary region). Both spectra are dominated by line emission from H- and He-like ions of C through S, and from Fe L-shell ions, but also include narrow radiative recombination continua, indicating that most of the soft X-ray emission arises in low-temperature (kT few eV) photoionized plasma. We confirm the conclusions of Kinkhabwala et al. (2002), based on XMM-Newton RGS observations, that the entire nuclear spectrum can be explained by recombination/radiative cascade following photoionization, and radiative decay following photoexcitation, with no evidence for hot, collisionally ionized plasma. In addition, this model also provides an excellent fit to the spectrum of the Secondary region, albeit with radial column densities a factor of three lower, as would be expected given its distance from the source of the ionizing continuum. The remarkable overlap and kinematical agreement of the optical and X-ray line emission, coupled with the need for a distribution of ionization parameter to explain the X-ray spectra, collectively imply the presence of a distribution of densities (over a few orders of magnitude) at each radius in the ionization cone. Relative abundances of all elements are consistent with Solar abundance, except for N, which is 2-3 times Solar. The long wavelength spectrum beyond 30 A is rich of L-shell transitions of Mg, Si, S, and Ar, and M-shell transitions of Fe. The velocity dispersion decreases with increasing ionization parameter, as deduced from these long wavelength lines and the Fe-L shell lines.

High resolution spectroscopic observations of Seyfert galaxies with Chandra and XMM-Newton allows us to study the detailed ionization and thermal structures of the X-ray absorbing/emitting material in the circumnuclear environment. The vast improvement in the spectral resolving power by more than an order of magnitude enables us, for the first time, to unambiguously distinguish the dominant line emission mechanisms and to measure its dynamical properties as well. The X-ray band harbors spectral transitions from a wide range of ionization states, including valence-shell transitions in K-shell and L-shell ions from most cosmically abundant elements, as well as inner-shell transitions of iron and other mid-Z elements, which can be probed through absorption measurements. The X-ray spectrum, therefore, provides simultaneous velocity and column density constraints of highly ionized to only slightly ionized gas harbored in many of these systems. We summarize recent results that have emerged from observations of Seyfert 2 galaxies with the grating spectrometers onboard Chandra and XMM-Newton. We give particular emphasis to an empirical physical model that we have developed based on the observed spectra, and how it can be used for comparative studies with Seyfert 1 galaxies to test the AGN unification scenarios.

We present the first high-resolution, soft-X-ray spectrum of the prototypical Seyfert 2 galaxy, NGC 1068. This spectrum was obtained with the XMM-Newton Reflection Grating Spectrometer. Emission lines from H-like and He-like low-Z ions (from C to Si) and Fe-L-shell ions dominate the spectrum. Strong, narrow radiative recombination continua (RRC) for several ions are also present, implying that most of the observed soft-X-ray emission arises in low-temperature (few eV) plasma. This plasma is photoionized by the inferred nuclear continuum (obscured along our line of sight), as in the unified model of active galactic nuclei (AGN). We find excess emission (compared with pure recombination) in all resonance lines (np to 1s) up to the photoelectric edge, demonstrating the importance of photoexcitation as well. We introduce a simple model of a cone of plasma irradiated by the nuclear continuum; the line emission we observe along our line of sight perpendicular to the cone is produced through recombination/radiative cascade following photoionization and radiative decay following photoexcitation. A remarkably good fit is obtained to the H-like/He-like ionic line series, with inferred radial ionic column densities consistent with recent observations of warm absorbers in Seyfert 1 galaxies. Previous Chandra imaging revealed a large (extending out to 500 pc) ionization cone containing most of the X-ray flux, implying that the warm absorber in NGC 1068 is a large-scale outflow. To explain the ionic column densities, a broad, flat distribution in the logarithm of the ionization parameter ($\xi=L_X/n_e r^2$) is necessary, spanning $\log\xi=0$--3. This suggests either radially-stratified ionization zones or the existence of a broad density distribution (spanning a few orders of magnitude) at each radius.

We present high-resolution soft-X-ray spectra of the prototypical Seyfert 2 galaxy, NGC 1068, taken with XMM-Newton RGS and Chandra LETGS. Its rich emission-line spectrum is dominated by recombination in a warm plasma (bright, narrow radiative recombination continua provide the ``smoking gun''), which is photoionized by the inferred nuclear power-law continuum. Radiative decay following photoexcitation of resonant transitions is also significant. A self-consistent model of an irradiated cone of gas is capable of reproducing the hydrogenic/heliumlike ionic line series in detail. The radial ionic column densities we infer are consistent with absorption measurements (the warm absorber) in Seyfert 1 galaxies. This strongly suggests that the emission spectrum we observe from NGC 1068 emanates from its warm absorber. The observed extent of the ionization-cone/warm absorber in NGC 1068 of about 300 pc implies that a large fraction of the gas associated with generic warm absorbers may typically exist on the hundreds-of-parsec scale rather than much closer to the nucleus (e.g., less than a parsec). Spatially-resolved spectroscopy using the LETGS of two distinct emission regions yields two noticeably different spectra. We show that these differences are solely due to differing radial column densities. A fairly flat distribution in ionization parameter is necessary to explain the inferred radial ionic column densities of all spectra. This must primarily be due to a broad density distribution at each radius, spanning roughly 0.1-100 cm$^{-3}$. (Abridged)

Soft X-ray spectroscopy of Seyfert 2 galaxies offers perhaps the best method to probe the possible connection between AGN activity and star formation. Obscuration of powerful radiation from the inferred nucleus allows for detailed study of circumnuclear emission regions. And soft X-ray spectroscopy of these regions allows for robust discrimination between warm gas radiatively driven by the AGN and hot collisionally-driven gas possibly associated with star formation. A simple model of a (bi-)cone of gas photoionized and photoexcited by a nuclear power-law continuum is sufficient to explain the soft X-ray spectra of all Seyfert 2 galaxies so far observed by the XMM-Newton and Chandra satellites. An upper limit of around 10 percent to an additional hot, collisionally-driven gas contribution to the soft X-ray regime appears to hold for five different Seyfert 2 galaxies, placing interesting constraints on circumnuclear star formation.

X-ray spectroscopy of Seyfert 2 galaxies provides an excellent probe of the circumnuclear environment in active galactic nuclei. The grating experiments on both Chandra and XMM-Newton have now provided the first high resolution spectra of several of the brightest Seyfert 2's. We present Chandra HETG data on Markarian 3 and XMM-Newton RGS data on NGC 1068. In both cases, the spectra are dominated by emission lines due to radiative recombination following photoionization, photoexcitation, and fluorescence. There is no evidence for any significant contribution from collisionally-heated gas.

Preliminary analysis of the X-ray spectrum of NGC 1068 obtained by the RGS spectrometer on board XMM-Newton is presented. A physically consistent model is developed in order to quantitatively describe the reprocessing of the central AGN continuum source into the discrete X-ray emission observed in Seyfert 2 galaxies. All the important atomic processes are taken into account, including photoexcitation, which has been neglected in some previous models. The model fits the high resolution NGC 1068 data very well, which implies that the contribution of hot collisional gas to the X-ray spectrum of NGC 1068 is negligible.

Giant pulses are short, intense outbursts of radio emission with a power-law intensity distribution that have been observed from the Crab Pulsar and PSR B1937+21. We have undertaken a systematic study of giant pulses from PSR B1937+21 using the Arecibo telescope at 430, 1420, and 2380 MHz. At 430 MHz, interstellar scattering broadens giant pulses to durations of $\sim50 \mu$secs, but at higher frequencies the pulses are very short, typically lasting only $\sim1$-$2 \mu$secs. At each frequency, giant pulses are emitted only in narrow ($\lsim10 \mu$s) windows of pulse phase located $\sim 55$-$70 \mu$sec after the main and interpulse peaks. Although some pulse-to-pulse jitter in arrival times is observed, the mean arrival phase appears stable; a timing analysis of the giant pulses yields precision competitive with the best average profile timing studies. We have measured the intensity distribution of the giant pulses, confirming a roughly power-law distribution with approximate index of -1.8, contributing $\gsim0.1%$ to the total flux at each frequency. We also find that the intensity of giant pulses falls off with a slightly steeper power of frequency than the ordinary radio emission.

We calculate the large-angle cross-correlation between the cosmic-microwave-background (CMB) temperature and the x-ray-background (XRB) intensity expected in an open Universe with cold dark matter (CDM) and a nearly scale-invariant spectrum of adiabatic density perturbations. Results are presented as a function of the nonrelativistic-matter density $\Omega_0$ (in units of the critical density) and the x-ray bias $b_x$ (evaluated at a redshift $z\simeq1$ in evolving-bias models) for both an open Universe and a flat cosmological-constant Universe. Recent experimental upper limits to the amplitude of this cross-correlation provide a new constraint to the $\Omega_0$-$b_x$ parameter space that open-CDM models (and the open-inflation models that produce them) must satisfy.

Rates for detection of weakly-interacting massive-particle (WIMP) dark matter are usually carried out assuming the Milky Way halo is an isothermal sphere. However, it is possible that our halo is not precisely spherical; it may have some bulk rotation; and the radial profile may differ from that of an isothermal sphere. In this paper, we calculate detection rates in observationally consistent alternative halo models that produce the same halo contributions to the local and asymptotic rotation speeds to investigate the effects of theoretical uncertainty of the WIMP spatial and velocity distribution. We use self-consistent models to take into account the effects of various mass distributions on the local velocity distribution. The local halo density may be increased up to a factor of 2 by flattening or by an alternative radial profile (which may also decrease the density slightly). However, changes in the WIMP velocity distribution in these models produce only negligible changes in the WIMP detection rate. Reasonable bulk rotations lead to only an $O(10%)$ effect on event rates. We also show how the nuclear recoil spectrum in a direct-detection experiment could provide information on the shape and rotation of the halo.