We investigate X-ray binary (XRB) luminosity function (XLF) scaling relations for Chandra detected populations of low-mass XRBs (LMXBs) within the footprints of 24 early-type galaxies. Our sample includes Chandra and HST observed galaxies at D < 25 Mpc that have estimates of the globular cluster (GC) specific frequency (SN) reported in the literature. As such, we are able to directly classify X-ray-detected sources as being either coincident with unrelated background/foreground objects, GCs, or sources that are within the fields of the galaxy targets. We model the GC and field LMXB population XLFs for all galaxies separately, and then construct global models characterizing how the LMXB XLFs vary with galaxy stellar mass and SN. We find that our field LMXB XLF models require a component that scales with SN, and has a shape consistent with that found for the GC LMXB XLF. We take this to indicate that GCs are "seeding" the galactic field LMXB population, through the ejection of GC-LMXBs and/or the diffusion of the GCs in the galactic fields themselves. However, we also find that an important LMXB XLF component is required for all galaxies that scales with stellar mass, implying that a substantial population of LMXBs are formed "in situ," which dominates the LMXB population emission for galaxies with SN < 2. For the first time, we provide a framework quantifying how directly-associated GC LMXBs, GC-seeded LMXBs, and in-situ LMXBs contribute to LMXB XLFs in the broader early-type galaxy population.

Astrophysics spans an enormous range of questions on scales from individual planets to the entire cosmos. To address the richness of 21st century astrophysics requires a corresponding richness of telescopes spanning all bands and all messengers. Much scientific benefit comes from having the multi-wavelength capability available at the same time. Most of these bands,or measurement sensitivities, require space-based missions. Historically, NASA has addressed this need for breadth with a small number of flagship-class missions and a larger number of Explorer missions. While the Explorer program continues to flourish, there is a large gap between Explorers and strategic missions. A fortunate combination of new astrophysics technologies with new, high capacity, low dollar-per-kg to orbit launchers, and new satellite buses allow for cheaper missions with capabilities approaching strategic mission levels. NASA has recognized these developments by calling for Probe-class mission ideas for mission studies, spanning most of the electromagnetic spectrum from GeV gamma-rays to the far infrared, and the new messengers of neutrinos and ultra-high energy cosmic rays. The key insight from the Probes exercise is that order-of-magnitude advances in science performance metrics are possible across the board for initial total cost estimates in the range 500M-1B dollars.

One of the central goals of LISA is the detection of gravitational waves from the merger of supermassive black holes. Contrary to stellar-mass black hole mergers, such events are expected to be rich X-ray sources due to the accretion of material from the circumbinary disks onto the black holes. The orbital dynamics before merger is also expected to modulate the resulting X-ray emission via Doppler boosting in a quasi-periodic way, and in a simple phase relation with the gravitational wave from the inspiral of the black holes. Detecting the X-ray source would enable a precise and early localization of the binary, thus allowing many telescopes to observe the very moment of the merger. Although identifying the correct X-ray source in the relatively large LISA sky localization will be challenging due to the presence of many confounding point sources, the quasi-periodic modulation may aid in the identification. We explore the practical feasibility of such idea. We simulate populations of merging supermassive black holes, their detection with LISA and their X-ray lightcurves using a simple model. Taking the parameters of the X-ray Telescope on the proposed NASA Transient Astrophysics Probe, we then design and simulate an observation campaign that searches for the modulated X-ray source while LISA is still observing the inspiral of the black holes. Assuming a fiducial LISA detection rate of $10$ mergers per year at redshift closer than $3.5$, we expect a few detections of modulated X-ray counterparts over the nominal duration of the LISA mission.

We present the analysis of Chandra and NuSTAR spectra of NGC 4968, a local (D$\sim$44 Mpc) 12$\mu$m-selected Seyfert 2 galaxy, enshrouded within Compton-thick layers of obscuring gas. We find no evidence of variability between the Chandra and NuSTAR observations (separated by 2 years), and between the two NuSTAR observations (separated by 10 months). Using self-consistent X-ray models, we rule out the scenario where the obscuring medium is nearly spherical and uniform, contradicting the results implied by the $<$10 keV Chandra spectrum. The line-of-sight column density, from intervening matter between the source and observer that intercepts the intrinsic AGN X-ray emission, is well within the Compton-thick regime, with a minimum column density of $2\times10^{24}$ cm$^{-2}$. The average global column density is high ($> 3\times10^{23}$ cm$^{-2}$), with both Compton-thick and Compton-thin solutions permitted depending on the X-ray spectral model. The spectral models provide a range of intrinsic AGN continuum parameters and implied 2-10 keV luminosities ($L_{\rm 2-10keV,intrinsic}$), where the higher end of $L_{\rm 2-10keV,intrinsic}$ is consistent with expectations from the 12$\mu$m luminosity ($L_{\rm 2-10keV,intrinisc} \sim 7\times10^{42}$ erg s$^{-1}$). Compared with Compton-thick AGN previously observed by {\it NuSTAR}, NGC 4968 is among the most intrinsically X-ray luminous. However, despite its close proximity and relatively high intrinsic X-ray luminosity, it is undetected by the 105 month Swift-BAT survey, underscoring the importance of multi-wavelength selection for obtaining the most complete census of the most hidden black holes.

The field of astronomy has arrived at a turning point in terms of size and complexity of both datasets and scientific collaboration. Commensurately, algorithms and statistical models have begun to adapt --- e.g., via the onset of artificial intelligence --- which itself presents new challenges and opportunities for growth. This white paper aims to offer guidance and ideas for how we can evolve our technical and collaborative frameworks to promote efficient algorithmic development and take advantage of opportunities for scientific discovery in the petabyte era. We discuss challenges for discovery in large and complex data sets; challenges and requirements for the next stage of development of statistical methodologies and algorithmic tool sets; how we might change our paradigms of collaboration and education; and the ethical implications of scientists' contributions to widely applicable algorithms and computational modeling. We start with six distinct recommendations that are supported by the commentary following them. This white paper is related to a larger corpus of effort that has taken place within and around the Petabytes to Science Workshops (https://petabytestoscience.github.io/).

We present a source catalog from the first deep hard X-ray ($E>10$ keV) survey of the Small Magellanic Cloud (SMC), the NuSTAR Legacy Survey of the SMC. We observed three fields, for a total exposure time of 1 Ms, along the bar of this nearby star-forming galaxy. Fields were chosen for their young stellar and accreting binary populations. We detected 10 sources above a 3$\sigma$ significance level (4$-$25 keV) and obtained upper limits on an additional 40 sources. We reached a 3$\sigma$ limiting luminosity in the 4$-$25 keV band of $\sim$ $10^{35}$ erg s$^{-1}$, allowing us to probe fainter X-ray binary (XRB) populations than has been possible with other extragalactic NuSTAR surveys. We used hard X-ray colors and luminosities to constrain the compact-object type, exploiting the spectral differences between accreting black holes and neutron stars at $E>10$ keV. Several of our sources demonstrate variability consistent with previously observed behavior. We confirmed pulsations for seven pulsars in our 3$\sigma$ sample. We present the first detection of pulsations from a Be-XRB, SXP305 (CXO J005215.4$-$73191), with an X-ray pulse period of $305.69\pm0.16$ seconds and a likely orbital period of $\sim$1160-1180 days. Bright sources ($\gtrsim 5\times 10^{36}$ erg s$^{-1}$) in our sample have compact-object classifications consistent with their previously reported types in the literature. Lower luminosity sources ($\lesssim 5\times 10^{36}$ erg s$^{-1}$) have X-ray colors and luminosities consistent with multiple classifications. We raise questions about possible spectral differences at low luminosity between SMC pulsars and the Galactic pulsars used to create the diagnostic diagrams.

X-ray data for quasar 4C 74.26 have previously been modeled with a broad Fe K$\alpha$ emission line and reflection continuum originating in the inner part of the accretion disk around the central supermassive black hole (SMBH), i.e. the strong gravity regime. We modeled broadband X-ray spectra from $Suzaku$ and $NuSTAR$ with MYTORUS, self-consistently accounting for Fe K$\alpha$ line emission, as well as direct and reflected continuum emission, from finite column density matter. A narrow Fe K$\alpha$ emission line originating in an X-ray reprocessor with solar Fe abundance far from the central SMBH is sufficient to produce excellent fits for all spectra. For the first time, we are able to measure the global, out of the line-of-sight column density to be in the range $\sim$$1.5$ to $\sim$$2.9\times10^{24}$ cm$^{-2}$, i.e. in the Compton thick regime, while the line-of-sight column density is Compton thin in all observations. The Fe K$\alpha$ emission line is unresolved in all but one observations. The Compton scattered continuum from distant matter removes the need for relativistic broadening of the Fe K$\alpha$ emission line, which is required for SMBH spin measurements. The resolved line observation can alternatively be modeled with a relativistic model but we do not find evidence for a truncated accretion disk model. We conclude that the X-ray emission in these 4C 74.26 data is unlikely to originate in the inner accretion disk region and thus cannot be used to measure SMBH spin.

We present a simple semi-numerical model designed to explore black hole growth and galaxy evolution. This method builds on a previous model for black hole accretion that uses a semi-numerical galaxy formation model and universal Eddington ratio distribution to describe the full AGN population by independently connecting galaxy and AGN growth to the evolution of the host dark matter halos. We fit observed X-ray luminosity functions up to a redshift of z ~ 4, as well as investigate the evolution of the Eddington ratio distributions. We find that the Eddington ratio distribution evolves with redshift such that the slope of the low-Eddington accretion rate distribution increases with cosmic time, consistent with the behavior predicted in hydrodynamical simulations for galaxies with different gas fractions. We also find that the evolution of our average Eddington ratio is correlated with observed star formation histories, supporting a picture in which black holes and galaxies evolve together in a global sense. We further confirm the impact of luminosity limits on observed galaxy and halo properties by applying selection criteria to our fiducial model and comparing to surveys across a wide range of redshifts.

We present new Chandra constraints on the X-ray luminosity functions (XLFs) of X-ray binary (XRB) populations, and their scaling relations, for a sample of 38 nearby galaxies (D = 3.4-29 Mpc). Our galaxy sample is drawn primarily from the Spitzer infrared nearby galaxy survey (SINGS), and contains a wealth of Chandra (5.8 Ms total) and multiwavelength data, allowing for star-formation rates (SFRs) and stellar masses (M*) to be measured on subgalactic scales. We divided the 2478 X-ray detected sources into 21 subsamples in bins of specific-SFR (sSFR = SFR/M*) and constructed XLFs. To model the XLF dependence on sSFR, we fit a global XLF model, containing contributions from high-mass XRBs (HMXBs), low-mass XRBs (LMXBs), and background sources from the cosmic X-ray background (CXB) that respectively scale with SFR, M*, and sky area. We find an HMXB XLF that is more complex in shape than previously reported and an LMXB XLF that likely varies with sSFR, potentially due to an age dependence. When applying our global model to XLF data for each individual galaxy, we discover a few galaxy XLFs that significantly deviate from our model beyond statistical scatter. Most notably, relatively low-metallicity galaxies have an excess of HMXBs above ~10^38 erg/s and elliptical galaxies that have relatively rich populations of globular clusters (GCs) show excesses of LMXBs compared to the global model. Additional modeling of how the XRB XLF depends on stellar age, metallicity, and GC specific frequency is required to sufficiently characterize the XLFs of galaxies.

Current X-ray observations and simulations show that gravitational lensing can be used to infer the structure near the event horizons of black holes, constrain the dynamics and evolution of black-hole accretion and outflows, test general relativity in the strong-gravity regime and place constraints on the evolution of dark matter in the lensing galaxies. These science goals currently cannot be achieved in a statistically large sample of z = 0.5 - 5 lensed quasars due to the limited capabilities of current X-ray telescopes and the relatively low number (~200) of known lensed quasars. The latter limitation will be resolved with the multi-band and wide-field photometric optical survey of LSST that is expected to lead to the discovery of > 4,000 additional gravitationally lensed systems. As we show in this white paper, these science goals can be reached with an X-ray telescope having a spatial resolution of <0.5arcsec to resolve the lensed images and a collecting area of >0.5 m^2 at 1 keV.

A large fraction of the baryons and most of the metals in the Universe are unaccounted for. They likely lie in extended galaxy halos, galaxy groups, and the cosmic web, and measuring their nature is essential to understanding galaxy formation. These environments have virial temperatures >10^5.5 K, so the gas should be visible in X-rays. Here we show the breakthrough capabilities of grating spectrometers to 1) detect these reservoirs of hidden metals and mass, and 2) quantify hot gas flows, turbulence, and rotation around the Milky Way and external galaxies. Grating spectrometers are essential instruments for future X-ray missions, and existing technologies provide 50-1500-fold higher throughput compared to current orbiting instruments.

What are the most important conditions and processes governing the growth of stellar-origin compact objects? The identification of compact object type as either black hole (BH) or neutron star (NS) is fundamental to understanding their formation and evolution. To date, time-domain determination of compact object type remains a relatively untapped tool. Measurement of orbital periods, pulsations, and bursts will lead to a revolution in the study of the demographics of NS and BH populations, linking source phenomena to accretion and galaxy parameters (e.g., star formation, metallicity). To perform these measurements over sufficient parameter space, a combination of a wide-field (>5000 deg^2) transient X-ray monitor over a dynamic energy range (~1-100 keV) and an X-ray telescope for deep surveys with <5 arcsec PSF half-energy width (HEW) angular resolution are required. Synergy with multiwavelength data for characterizing the underlying stellar population will transform our understanding of the time domain properties of transient sources, helping to explain details of supernova explosions and gravitational wave event rates.

Stellar feedback -- stars regulating further star formation through the injection of energy and momentum into the interstellar medium -- operates through a complex set of processes that originate in star clusters but shape entire galaxies. A mature theory of stellar feedback is essential to a complete theory of star and galaxy formation, but the energy and momentum injected by hot gas into its surroundings remains unclear. With a next-generation X-ray observatory, breakthrough progress can be made through precision measurements of the temperature, density, velocity, and abundances of hot gas on scales of star clusters to galactic superwinds.

Much of the baryonic matter in the Universe, including the most active and luminous sources, are best studied in the X-ray band. Key advances in X-ray optics and detectors have paved the way for the Advanced X-ray Imaging Satellite (AXIS), a Probe-class mission that is a major improvement over Chandra, which has generated a steady stream of important discoveries for the past 2 decades. AXIS can be launched in the late 2020s and will transform our understanding in several major areas of astrophysics, including the growth and fueling of supermassive black holes, galaxy formation and evolution, the microphysics of cosmic plasmas, the time-variable universe, and a wide variety of cutting-edge studies. Relative to Chandra, the AXIS PSF is nearly twice as sharp on-axis; its field of view for subarcsecond imaging 70 times larger by area; its effective area at 1 keV is 10 times larger. The low-Earth orbit ensures a low and stable detector background, resulting in 50 times greater sensitivity than Chandra for extended sources. AXIS has a rapid repointing response with operations similar to Swift, but is 100 times more sensitive for time-domain science. These capabilities open up a vast discovery space and complement the next generation of astronomical observatories. A high-spectral-resolution mission (Athena) operating at the same time as a high-angular-resolution mission (AXIS) greatly increases the range of scientific discovery. AXIS will use lightweight X-ray optics made of thin single-crystal silicon mirrors developed at NASA Goddard. The detector array builds on a long legacy of X-ray CCD and provides improved photon localization, much faster readout time, and broader energy band. The estimated mission costs are consistent with the $1B Probe mission cost guideline.

Ultraluminous X-ray (ULX) pulsars are a new class of object powered by apparent super-critical accretion onto magnetized neutron stars. Three sources in this class identified so far; M82 X-2, NGC 5907 ULX-1 and NGC 7793 P13, have been found to have two properties in common; $\sim1$-s spin periods, and for NGC 5907 ULX-1 and NGC 7793 P13 periodic X-ray flux modulations on timescales of $\sim60-80$ days. M82 X-2 resides in a crowded field that includes the ULX M82 X-1 separated from X-2 by 5", as well as other bright point sources. A 60-day modulation has been observed from the region but the origin has been difficult to identify; both M82 X-1 and X-2 have been suggested as the source. In this paper we present the analysis of a systematic monitoring campaign by Chandra, the only X-ray telescope capable of resolving the crowded field. From a simple Lomb-Scargle periodogram analysis and a more sophisticated Gaussian Process analysis we find that only X-2 exhibits a periodic signal around 60 days supporting previous claims that it is the origin. We also construct a phase-averaged flux profile of the modulations from higher cadence Swift/XRT data and find that the flux variations in the Chandra data are fully consistent with the flux profile. Since the orbit of the neutron star and its companion is known to be 2.5 days, the $\sim60$-day period must be super-orbital in origin. The flux of the modulations varies by a factor of $\sim$100 from minimum to maximum, with no evidence for spectral variations, making the origin difficult to explain.

We report the analysis of simultaneous XMM-Newton+NuSTAR observations of two low-luminosity Active Galactic Nuclei (LLAGN), NGC 3998 and NGC 4579. We do not detect any significant variability in either source over the ~3 day length of the NuSTAR observations. The broad-band 0.5-60 keV spectrum of NGC 3998 is best fit with a cutoff power-law, while the one for NGC 4579 is best fit with a combination of a hot thermal plasma model, a power-law, and a blend of Gaussians to fit an Fe complex observed between 6 and 7 keV. Our main spectral results are the following: (1) neither source shows any reflection hump with a $3\sigma$ reflection fraction upper-limits $R<0.3$ and $R<0.18$ for NGC 3998 and NGC 4579, respectively; (2) the 6-7 keV line complex in NGC 4579 could either be fit with a narrow Fe K line at 6.4 keV and a moderately broad Fe XXV line, or 3 relatively narrow lines, which includes contribution from Fe XXVI; (3) NGC 4579 flux is 60% brighter than previously detected with XMM-Newton, accompanied by a hardening in the spectrum; (4) we measure a cutoff energy $E_{\rm cut}=107_{-18}^{+27}$ keV in NGC 3998, which represents the lowest and best constrained high-energy cutoff ever measured for an LLAGN; (5) NGC 3998 spectrum is consistent with a Comptonization model with either a sphere ($\tau\approx3\pm1$) or slab ($\tau\approx1.2\pm0.6$) geometry, corresponding to plasma temperatures between 20 and 150 keV. We discuss these results in the context of hard X-ray emission from bright AGN, other LLAGN, and hot accretion flow models.

We present nearly simultaneous NuSTAR and XMM-Newton observations of the nearby (832 kpc) ultraluminous X-ray source (ULX) M33 X-8. M33 X-8 has a 0.3-10 keV luminosity of LX ~ 1.4 x 10^39 erg/s, near the boundary of the "ultraluminous" classification, making it an important source for understanding the link between typical Galactic X-ray binaries and ULXs. Past studies have shown that the 0.3-10 keV spectrum of X-8 can be characterized using an advection-dominated accretion disk model. We find that when fitting to our NuSTAR and XMM-Newton observations, an additional high-energy (>10 keV) Comptonization component is required, which allows us to rule out single advection-dominated disk and classical sub-Eddington models. With our new constraints, we analyze XMM-Newton data taken over the last 17 years to show that small (~30%) variations in the 0.3-10 keV flux of M33 X-8 result in spectral changes similar to those observed for other ULXs. The two most likely phenomenological scenarios suggested by the data are degenerate in terms of constraining the nature of the accreting compact object (i.e., black hole versus neutron star). We further present a search for pulsations using our suite of data; however, no clear pulsations are detected. Future observations designed to observe M33 X-8 at different flux levels across the full 0.3-30 keV range would significantly improve our constraints on the nature of this important source.

We present a broadband X-ray spectral analysis of the M51 system, including the dual active galactic nuclei (AGN) and several off-nuclear point sources. Using a deep observation by NuSTAR, new high-resolution coverage of M51b by Chandra, and the latest X-ray torus models, we measure the intrinsic X-ray luminosities of the AGN in these galaxies. The AGN of M51a is found to be Compton thick, and both AGN have very low accretion rates ($\lambda_{\rm Edd} <10^{-4}$). The latter is surprising considering that the galaxies of M51 are in the process of merging, which is generally predicted to enhance nuclear activity. We find that the covering factor of the obscuring material in M51a is $0.26 \pm 0.03$, consistent with the local AGN obscured fraction at $L_{\rm X}\sim 10^{40}$ erg s$^{-1}$. The substantial obscuring column does not support theories that the torus, presumed responsible for the obscuration, disappears at these low accretion luminosities. However, the obscuration may have resulted from the gas infall driven by the merger rather than the accretion process. We report on several extra-nuclear sources with $L_{\rm X}>10^{39}$ erg s$^{-1}$ and find that a spectral turnover is present below 10 keV in most such sources, in line with recent results on ultraluminous X-ray sources.

Nearby galaxy surveys have long classified X-ray binaries (XRBs) by the mass category of their donor stars (high-mass and low-mass). The NuSTAR observatory, which provides imaging data at E $>10$ keV, has enabled the classification of extragalactic XRBs by their compact object type: neutron star (NS) or black hole (BH). We analyzed NuSTAR/Chandra/XMM-Newton observations from a NuSTAR-selected sample of 12 galaxies within 5 Mpc having stellar masses ($M_{\star}$) $10^{7-11}$ $M_{\odot}$ and star formation rates (SFR) $\approx0.01-15$ $M_{\odot}$ yr$^{-1}$. We detect 128 NuSTAR sources to a sensitivity of $\approx10^{38}$ erg s$^{-1}$. Using NuSTAR color-intensity and color-color diagrams we classify 43 of these sources as candidate NS and 47 as candidate BH. We further subdivide BH by accretion states (soft, intermediate, and hard) and NS by weak (Z/Atoll) and strong (accreting pulsar) magnetic field. Using 8 normal (Milky Way-type) galaxies in the sample, we confirm the relation between SFR and galaxy X-ray point source luminosity in the 4-25 and 12-25 keV energy bands. We also constrain galaxy X-ray point source luminosity using the relation $L_{\rm{X}}=\alpha M_{\star}+\beta\text{SFR}$, finding agreement with previous work. The XLF of all sources in the 4-25 and 12-25 keV energy bands matches with the $\alpha=1.6$ slope for high-mass XRBs. We find that NS XLFs suggest a decline beginning at the Eddington limit for a 1.4 $M_{\odot}$ NS, whereas the BH fraction shows an approximate monotonic increase in the 4-25 and 12-25keV energy bands. We calculate the overall ratio of BH to NS to be $\approx1$ for 4-25 keV and $\approx2$ for 12-25 keV.

We present 15 high mass X-ray binary (HMXB) candidates in the disk of M31 for which we are able to infer compact object type, spectral type of the donor star, and age using multiwavelength observations from NuSTAR, Chandra, and the Hubble Space Telescope (HST). The hard X-ray colors and luminosities from NuSTAR permit the tentative classification of accreting X-ray binary systems by compact object type, distinguishing black hole from neutron star systems. We find hard state black holes, pulsars, and non-magnetized neutron stars associated with optical point source counterparts with similar frequency. We also find nine non-magnetized neutron stars coincident with globular clusters and an equal number of pulsars with and without point source optical counterparts. We perform spectral energy distribution (SED) fitting for the most likely optical counterparts to the HMXB candidates, finding 7 likely high mass stars and 1 possible red Helium burning star. The remaining 7 HMXB optical counterparts have poor SED fits, so their companion stars remain unclassified. Using published star formation histories, we find that the majority of HMXB candidates --- X-ray sources with UV-bright point source optical counterpart candidates --- are found in regions with star formation bursts less than 50 Myr ago, with 3 associated with young stellar ages (<10 Myr). This is consistent with similar studies of HMXB populations in the Magellanic Clouds, M33, NGC 300, and NGC 2403.

We present a detailed, broadband X-ray spectral analysis of the ULX pulsar NGC 7793 P13, a known super-Eddington source, utilizing data from the $XMM$-$Newton$, $NuSTAR$ and $Chandra$ observatories. The broadband $XMM$-$Newton+NuSTAR$ spectrum of P13 is qualitatively similar to the rest of the ULX sample with broadband coverage, suggesting that additional ULXs in the known population may host neutron star accretors. Through time-averaged, phase-resolved and multi-epoch studies, we find that two non-pulsed thermal blackbody components with temperatures $\sim$0.5 and $\sim$1.5 keV are required to fit the data below 10 keV, in addition to a third continuum component which extends to higher energies and is associated with the pulsed emission from the accretion column. The characteristic radii of the thermal components appear to be similar, and are too large to be associated with the neutron star itself, so the need for two components likely indicates the accretion flow outside the magnetosphere is complex. We suggest a scenario in which the thick inner disc expected for super-Eddington accretion begins to form, but is terminated by the neutron star's magnetic field soon after its onset, implying a limit of $B \lesssim 6 \times 10^{12}$ G for the dipolar component of the central neutron star's magnetic field. Evidence of similar termination of the disc in other sources may offer a further means of identifying additional neutron star ULXs. Finally, we examine the spectrum exhibited by P13 during one of its unusual 'off' states. These data require both a hard powerlaw component, suggesting residual accretion onto the neutron star, and emission from a thermal plasma, which we argue is likely associated with the P13 system.

We present a new technique for empirically calibrating how the X-ray luminosity function (XLF) of X-ray binary (XRB) populations evolves following a star-formation event. We first utilize detailed stellar population synthesis modeling of far-UV to far-IR photometry of the nearby face-on spiral galaxy M51 to construct maps of the star-formation histories (SFHs) on subgalactic (~400 pc) scales. Next, we use the ~850 ks cumulative Chandra exposure of M51 to identify and isolate 2-7 keV detected point sources within the galaxy, and we use our SFH maps to recover the local properties of the stellar populations in which each X-ray source is located. We then divide the galaxy into various subregions based on their SFH properties (e.g., star-formation rate [SFR] per stellar mass [M*] and mass-weighted stellar age) and group the X-ray point sources according to the characteristics of the regions in which they are found. Finally, we construct and fit a parameterized XLF model that quantifies how the XLF shape and normalization evolves as a function of the XRB population age. Our best-fit model indicates the XRB XLF per unit stellar mass declines in normalization, by ~3-3.5 dex, and steepens in slope from ~10 Myr to ~10 Gyr. We find that our technique recovers results from past studies of how XRB XLFs and XRB luminosity scaling relations vary with age and provides a self-consistent picture for how the XRB XLF evolves with age.

We present Lightning, a new spectral energy distribution (SED) fitting procedure, capable of quickly and reliably recovering star formation history (SFH) and extinction parameters. The SFH is modeled as discrete steps in time. In this work, we assumed lookback times of 0-10 Myr, 10-100 Myr, 0.1-1 Gyr, 1-5 Gyr, and 5-13.6 Gyr. Lightning consists of a fully vectorized inversion algorithm to determine SFH step intensities and combines this with a grid-based approach to determine three extinction parameters. We apply our procedure to the extensive FUV-to-FIR photometric data of M51, convolved to a common spatial resolution and pixel scale, and make the resulting maps publicly available. We recover, for M51a, a peak star formation rate (SFR) between 0.1 and 5 Gyr ago, with much lower star formation activity over the last 100 Myr. For M51b, we find a declining SFR toward the present day. In the outskirt regions of M51a, which includes regions between M51a and M51b, we recover a SFR peak between 0.1 and 1 Gyr ago, which corresponds to the effects of the interaction between M51a and M51b. We utilize our results to (1) illustrate how UV+IR hybrid SFR laws vary across M51, and (2) provide first-order estimates for how the IR luminosity per unit stellar mass varies as a function of the stellar age. From the latter result, we find that IR emission from dust heated by stars is not always associated with young stars, and that the IR emission from M51b is primarily powered by stars older than 5 Gyr.

In studies of the connection between active galactic nuclei (AGN) and their host galaxies there is widespread disagreement on some key aspects stemming largely from a lack of understanding of the nature of the full underlying AGN population. Recent attempts to probe this connection utilize both observations and simulations to correct for a missed population, but presently are limited by intrinsic biases and complicated models. We take a simple simulation for galaxy evolution and add a new prescription for AGN activity to connect galaxy growth to dark matter halo properties and AGN activity to star formation. We explicitly model selection effects to produce an "observed" AGN population for comparison with observations and empirically motivated models of the local universe. This allows us to bypass the difficulties inherent in many models which attempt to infer the AGN population by inverting selection effects. We investigate the impact of selecting AGN based on thresholds in luminosity or Eddington ratio on the "observed" AGN population. By limiting our model AGN sample in luminosity, we are able to recreate the observed local AGN luminosity function and specific star formation-stellar mass distribution, and show that using an Eddington ratio threshold introduces less bias into the sample by selecting the full range of growing black holes, despite the challenge of selecting low mass black holes. We find that selecting AGN using these various thresholds yield samples with different AGN host galaxy properties.

We report the identification of a bright hard X-ray source dominating the M31 bulge above 25 keV from a simultaneous NuSTAR-Swift observation. We find that this source is the counterpart to Swift J0042.6+4112, which was previously detected in the Swift BAT All-sky Hard X-ray Survey. This Swift BAT source had been suggested to be the combined emission from a number of point sources; our new observations have identified a single X-ray source from 0.5 to 50 keV as the counterpart for the first time. In the 0.5-10 keV band, the source had been classified as an X-ray binary candidate in various Chandra and XMM studies; however, since it was not clearly associated with Swift J0042.6+4112, the previous E < 10 keV observations did not generate much attention. This source has a spectrum with a soft X-ray excess (kT~ 0.2 keV) plus a hard spectrum with a power law of Gamma ~ 1 and a cutoff around 15-20 keV, typical of the spectral characteristics of accreting pulsars. Unfortunately, any potential pulsation was undetected in the NuSTAR data, possibly due to insufficient photon statistics. The existing deep HST images exclude high-mass (>3 Msun) donors at the location of this source. The best interpretation for the nature of this source is an X-ray pulsar with an intermediate-mass (<3 Msun) companion or a symbiotic X-ray binary. We discuss other possibilities in more detail.

We present results from four new broadband X-ray observations of the extreme ultraluminous X-ray source Holmberg IX X-1 ($L_{\rm{X}} > 10^{40}$ erg/s), performed by $Suzaku$ and $NuSTAR$ in coordination. Combined with the archival data, we now have broadband observations of this remarkable source from six separate epochs. Two of these new observations probe lower fluxes than seen previously, allowing us to extend our knowledge of the broadband spectral variability exhibited. The spectra are well fit by two thermal blackbody components, which dominate the emission below 10 keV, as well as a steep ($\Gamma \sim 3.5$) powerlaw tail which dominates above $\sim$15 keV. Remarkably, while the 0.3-10.0 keV flux varies by a factor of $\sim$3 between all these epochs, the 15-40 keV flux varies by only $\sim$20%. Although the spectral variability is strongest in the $\sim$1-10 keV band, both of the thermal components are required to vary when all epochs are considered. We also re-visit the search for iron absorption features, leveraging the high-energy $NuSTAR$ data to improve our sensitivity to extreme velocity outflows in light of the ultra-fast outflow recently detected in NGC 1313 X-1. Iron absorption from a similar outflow along our line of sight can be ruled out in this case. We discuss these results in the context of super-Eddington accretion models that invoke a funnel-like geometry for the inner flow, and propose a scenario in which we have an almost face-on view of a funnel that expands to larger radii with increasing flux, resulting in an increasing degree of geometrical collimation for the emission from intermediate temperature regions.

We present the Chandra imaging and spectral analysis of NGC 4968, a nearby (z = 0.00986) Seyfert 2 galaxy. We discover extended ($\sim$1 kpc) X-ray emission in the soft band (0.5 - 2 keV) that is neither coincident with the narrow line region nor the extended radio emission. Based on spectral modeling, it is linked to on-going star formation ($\sim$2.6-4 M$_{\sun}$ yr$^{-1}$). The soft emission at circumnuclear scales (inner $\sim$400 pc) originates from hot gas, with kT $\sim$ 0.7 keV, while the most extended thermal emission is cooler (kT $\sim$ 0.3 keV). We refine previous measurements of the extreme Fe K$\alpha$ equivalent width in this source (EW = 2.5$^{+2.6}_{-1.0}$ keV), which suggests the central engine is completely embedded within Compton-thick levels of obscuration. Using physically motivated models fit to the Chandra spectrum, we derive a Compton-thick column density ($N_{\rm H} > 1.25\times10^{24}$ cm$^{-2}$) and an intrinsic hard (2-10 keV) X-ray luminosity of $\sim$3-8$\times 10^{42}$ erg s$^{-1}$ (depending on the presumed geometry of the obscurer), which is over two orders of magnitude larger than that observed. The large Fe K$\alpha$ EW suggests a spherical covering geometry, which could be confirmed with X-ray measurements above 10 keV. NGC 4968 is similar to other active galaxies that exhibit extreme Fe K$\alpha$ EWs (i.e., $>$2 keV) in that they also contain ongoing star formation. This work supports the idea that gas associated with nuclear star formation may increase the covering factor of the enshrouding gas and play a role in obscuring AGN.

We report the detection of coherent pulsations from the ultraluminous X-ray source NGC 7793 P13. The ~0.42s nearly sinusoidal pulsations were initially discovered in broadband X-ray observations using XMM-Newton and NuSTAR taken in 2016. We subsequently also found pulsations in archival XMM-Newton data taken in 2013 and 2014. The significant (>>5 sigma) detection of coherent pulsations demonstrates that the compact object in P13 is a neutron star with an observed peak luminosity of ~1e40 erg/s (assuming isotropy), well above the Eddington limit for a 1.4 M_sun accretor. This makes P13 the second ultraluminous X-ray source known to be powered by an accreting neutron star. The pulse period varies between epochs, with a slow but persistent spin up over the 2013-2016 period. This spin-up indicates a magnetic field of B ~ 1.5e12 G, typical of many accreting pulsars. The most likely explanation for the extreme luminosity is a high degree of beaming, however this is difficult to reconcile with the sinusoidal pulse profile.

We present results from sensitive, multi-epoch NuSTAR observations of the late-type star-forming galaxy M83 (d=4.6 Mpc), which is the first investigation to spatially resolve the hard (E>10 keV) X-ray emission of this galaxy. The nuclear region and ~ 20 off-nuclear point sources, including a previously discovered ultraluminous X-ray (ULX) source, are detected in our NuSTAR observations. The X-ray hardnesses and luminosities of the majority of the point sources are consistent with hard X-ray sources resolved in the starburst galaxy NGC 253. We infer that the hard X-ray emission is most likely dominated by intermediate accretion state black hole binaries and neutron star low-mass X-ray binaries (Z-sources). We construct the X-ray binary luminosity function (XLF) in the NuSTAR band for an extragalactic environment for the first time. The M83 XLF has a steeper XLF than the X-ray binary XLF in NGC 253, consistent with previous measurements by Chandra at softer X-ray energies. The NuSTAR integrated galaxy spectrum of M83 drops quickly above 10 keV, which is also seen in the starburst galaxies NGC253, NGC 3310 and NGC 3256. The NuSTAR observations constrain any AGN to be either highly obscured or to have an extremely low luminosity of $_{\sim}^<$10$^{38}$ erg/s (10-30 keV), implying it is emitting at a very low Eddington ratio. An X-ray point source consistent with the location of the nuclear star cluster with an X-ray luminosity of a few times 10$^{38}$ erg/s may be a low-luminosity AGN but is more consistent with being an X-ray binary.

We present the results of a joint Swift-NuSTAR spectroscopy campaign on M31. We focus on the five brightest globular cluster X-ray sources in our fields. Two of these had previously been argued to be black hole candidates on the basis of apparent hard-state spectra at luminosities above those for which neutron stars are in hard states. We show that these two sources are likely to be Z-sources (i.e. low magnetic field neutron stars accreting near their Eddington limits), or perhaps bright atoll sources (low magnetic field neutron stars which are just a bit fainter than this level) on the basis of simultaneous Swift and NuSTAR spectra which cover a broader range of energies. These new observations reveal spectral curvature above 6-8 keV that would be hard to detect without the broader energy coverage the NuSTAR data provide relative to Chandra and XMM-Newton. We show that the other three sources are also likely to be bright neutron star X-ray binaries, rather than black hole X-ray binaries. We discuss why it should already have been realized that it was unlikely that these objects were black holes on the basis of their being persistent sources, and we re-examine past work which suggested that tidal capture products would be persistently bright X-ray emitters. We discuss how this problem is likely due to neglecting disk winds in older work that predict which systems will be persistent and which will be transient.

We have studied high mass X-ray binary (HMXB) populations within two low-metallicity, starburst galaxies, Haro 11 and VV 114. These galaxies serve as analogs to high-redshift (z>2) Lyman break galaxies, and within the larger sample of Lyman break analogs (LBAs) are sufficiently nearby (<87 Mpc) to be spatially-resolved by Chandra. Previous studies of the X-ray emission in LBAs have found that the 2-10 keV luminosity per star formation rate (SFR) in these galaxies is elevated, potentially because of their low metallicities (12+log[O/H]= 8.3-8.4). Theoretically, the progenitors of XRBs forming in lower metallicity environments lose less mass from stellar winds over their lifetimes, producing more massive compact objects (i.e., neutron stars and black holes), and thus resulting in more numerous and luminous HMXBs per SFR. We have performed an in-depth study of the only two LBAs that have spatially-resolved 2-10 keV emission with Chandra to present the bright end of the X-ray luminosity distribution of HMXBs (L$_X>10^{39}$ erg/s; ultraluminous X-ray sources, ULXs) in these low-metallicity galaxies, based on 8 detected ULXs. Comparing with the star-forming galaxy X-ray luminosity function (XLF) presented by Mineo et al. (2012), Haro 11 and VV 114 host ~4 times more L$_X>10^{40}$ erg/s sources than expected given their SFRs. We simulate the effects of source blending from crowded lower luminosity HMXBs using the star-forming galaxy XLF and then vary the XLF shapes until we reproduce the observed point source luminosity distributions. We find that these LBAs have a shallower bright end slope than the standard XLF. If we conservatively assume that the brightest X-ray source from each galaxy is powered by an AGN rather than a HMXB and eliminate these sources from consideration, the luminosity distribution becomes poorly constrained but does appear to be consistent with a standard XLF. [Abridged]

The recent discovery by Bachetti et al. (2014) of a pulsar in M82 that can reach luminosities of up to 10^40 ergs s^-1, a factor of ~100 the Eddington luminosity for a 1.4 Msol compact object, poses a challenge for accretion physics. In order to better understand the nature of this source and its duty cycle, and in the light of several physical models that have been subsequently published, we conduct a spectral and temporal analysis of the 0.5-8 keV X-ray emission from this source from 15 years of Chandra observations. We fit the Chandra spectra of the pulsar with a power-law model and a disk black body model, subjected to interstellar absorption in M82. We carefully assess for the effect of pile-up in our observations, where 4/19 observations have a pile-up fraction >10%, which we account for during spectral modeling with a convolution model. When fitted with a power-law model, the average photon index when the source is at high luminosity (L_X>10^39 ergs s^-1) is Gamma=1.33+/-0.15. For the disk black body model, the average temperature is T=3.24+/-0.65 keV, consistent with other luminous X-ray pulsars. We also investigated the inclusion of a soft excess component and spectral break, finding that the spectra are also consistent with these features common to luminous X-ray pulsars. In addition, we present spectral analysis from NuSTAR over the 3-50 keV range where we have isolated the pulsed component. We find that the pulsed emission in this band is best fit by a power-law with a high-energy cut-off, where Gamma=0.6+/-0.3 and E_C=14^{+5}_{-3} keV. While the pulsar has previously been identified as a transient, we find from our longer-baseline study that it has been remarkably active over the 15-year period, where for 9/19 (47%) observations that we analyzed, the pulsar appears to be emitting at a luminosity in excess of 10^39 ergs s^-1, greater than 10 times its Eddington limit.

We present a NuSTAR, Chandra, and XMM--Newton survey of nine of the nearest ultraluminous infrared galaxies (ULIRGs). The unprecedented sensitivity of NuSTAR at energies above 10 keV enables spectral modeling with far better precision than was previously possible. Six of the nine sources observed were detected sufficiently well by NuSTAR to model in detail their broadband X-ray spectra, and recover the levels of obscuration and intrinsic X-ray luminosities. Only one source (IRAS 13120--5453) has a spectrum consistent with a Compton--thick AGN, but we cannot rule out that a second source (Arp 220) harbors an extremely highly obscured AGN as well. Variability in column density (reduction by a factor of a few compared to older observations) is seen in IRAS 05189--2524 and Mrk 273, altering the classification of these border-line sources from Compton-thick to Compton-thin. The ULIRGs in our sample have surprisingly low observed fluxes in high energy (>10 keV) X-rays, especially compared to their bolometric luminosities. They have lower ratios of unabsorbed 2--10 keV to bolometric luminosity, and unabsorbed 2--10 keV to mid-IR [O IV] line luminosity than do Seyfert 1 galaxies. We identify IRAS 08572+3915 as another candidate intrinsically X-ray weak source, similar to Mrk 231. We speculate that the X-ray weakness of IRAS 08572+3915 is related to its powerful outflow observed at other wavelengths.

We present nearly simultaneous Chandra and NuSTAR observations of two actively star-forming galaxies within 50 Mpc: NGC 3256 and NGC 3310. Both galaxies are detected by both Chandra and NuSTAR, which together provide the first-ever spectra of these two galaxies spanning 0.3-30 keV. The X-ray emission from both galaxies is spatially resolved by Chandra; we find that hot gas dominates the E < 1-3 keV emission while ultraluminous X-ray sources (ULXs) dominate at E > 1-3 keV. The NuSTAR galaxy-wide spectra of both galaxies follow steep power-law distributions with Gamma ~ 2.6 at E > 5-7 keV, similar to the spectra of bright individual ULXs and other galaxies that have been studied by NuSTAR. We find that both NGC 3256 and NGC 3310 have X-ray detected sources coincident with nuclear regions; however, the steep NuSTAR spectra of both galaxies restricts these sources to be either low luminosity AGN or non-AGN in nature (e.g., ULXs or crowded X-ray sources that reach L2-10 keV ~ 10^40 erg/s cannot be ruled out). Combining our constraints on the 0.3-30 keV spectra of NGC 3256 and NGC 3310 with equivalent measurements for nearby star-forming galaxies M83 and NGC 253, we analyze the SFR-normalized spectra of these starburst galaxies. The spectra of all four galaxies show sharply declining power-law slopes above 3-6 keV due to ULX populations. Our observations therefore constrain the average spectra of luminous accreting binaries (i.e., ULXs). This result is similar to the super-Eddington accreting ULXs that have been studied individually in a targeted NuSTAR ULX program. We also find that NGC 3310 exhibits a factor of ~3-10 elevation of X-ray emission over the other star-forming galaxies. We argue that the excess is most likely explained by the relatively low metallicity of the young stellar population in NGC 3310.

We report on simultaneous observations of the local starburst system Arp 299 with NuSTAR and Chandra, which provides the first resolved images of this galaxy up to energies of ~ 45 keV. Fitting the 3-40 keV spectrum reveals a column density of $N_{\rm H}$ ~ 4 x10^{24} cm^{-2}, characteristic of a Compton-thick AGN, and a 10-30 keV luminosity of 1.2x 10^{43} ergs s^{-1}. The hard X-rays detected by NuSTAR above 10 keV are centered on the western nucleus, Arp 299-B, which previous X-ray observations have shown to be the primary source of neutral Fe-K emission. Other X-ray sources, including Arp 299-A, the eastern nucleus which is also thought to harbor an AGN, as well as X-ray binaries, contribute $\lesssim 10%$ to the 10-20 keV emission from the Arp 299 system. The lack of significant emission above 10 keV other than that attributed to Arp 299-B suggests that: a) any AGN in Arp 299-A must be heavily obscured ($N_{\rm H}$ > 10^{24} cm^{-2}) or have a much lower luminosity than Arp 299-B and b) the extranuclear X-ray binaries have spectra that cut-off above ~10 keV. Such soft spectra are characteristic of ultraluminous X-ray (ULX) sources observed to date by NuSTAR.

We present results for two Ultraluminous X-ray Sources (ULXs), IC 342 X-1 and IC 342 X-2, using two epochs of XMM-Newton and NuSTAR observations separated by $\sim$7 days. We observe little spectral or flux variability above 1 keV between epochs, with unabsorbed 0.3--30 keV luminosities being $1.04^{+0.08}_{-0.06} \times 10^{40}$ erg s$^{-1}$ for IC 342 X-1 and $7.40\pm0.20 \times 10^{39}$ erg s$^{-1}$ for IC 342 X-2, so that both were observed in a similar, luminous state. Both sources have a high absorbing column in excess of the Galactic value. Neither source has a spectrum consistent with a black hole binary in low/hard state, and both ULXs exhibit strong curvature in their broadband X-ray spectra. This curvature rules out models that invoke a simple reflection-dominated spectrum with a broadened iron line and no cutoff in the illuminating power-law continuum. X-ray spectrum of IC 342 X-1 can be characterized by a soft disk-like black body component at low energies and a cool, optically thick Comptonization continuum at high energies, but unique physical interpretation of the spectral components remains challenging. The broadband spectrum of IC 342 X-2 can be fit by either a hot (3.8 keV) accretion disk, or a Comptonized continuum with no indication of a seed photon population. Although the seed photon component may be masked by soft excess emission unlikely to be associated with the binary system, combined with the high absorption column, it is more plausible that the broadband X-ray emission arises from a simple thin blackbody disk component. Secure identification of the origin of the spectral components in these sources will likely require broadband spectral variability studies.

Prior to the launch of NuSTAR, it was not feasible to spatially resolve the hard (E > 10 keV) emission from galaxies beyond the Local Group. The combined NuSTAR dataset, comprised of three ~165 ks observations, allows spatial characterization of the hard X-ray emission in the galaxy NGC 253 for the first time. As a follow up to our initial study of its nuclear region, we present the first results concerning the full galaxy from simultaneous NuSTAR, Chandra, and VLBA monitoring of the local starburst galaxy NGC 253. Above ~10 keV, nearly all the emission is concentrated within 100" of the galactic center, produced almost exclusively by three nuclear sources, an off-nuclear ultraluminous X-ray source (ULX), and a pulsar candidate that we identify for the first time in these observations. We detect 21 distinct sources in energy bands up to 25 keV, mostly consisting of intermediate state black hole X-ray binaries. The global X-ray emission of the galaxy - dominated by the off-nuclear ULX and nuclear sources, which are also likely ULXs - falls steeply (photon index >~ 3) above 10 keV, consistent with other NuSTAR-observed ULXs, and no significant excess above the background is detected at E > 40 keV. We report upper limits on diffuse inverse Compton emission for a range of spatial models. For the most extended morphologies considered, these hard X-ray constraints disfavor a dominant inverse Compton component to explain the {\gamma}-ray emission detected with Fermi and H.E.S.S. If NGC 253 is typical of starburst galaxies at higher redshift, their contribution to the E > 10 keV cosmic X-ray background is < 1%.

We present results from the coordinated broadband X-ray observations of the extreme ultraluminous X-ray source Holmberg IX X-1 performed by NuSTAR, XMM-Newton and Suzaku in late 2012. These observations provide the first high-quality spectra of Holmberg IX X-1 above 10 keV to date, extending the X-ray coverage of this remarkable source up to ~30 keV. Broadband observations were undertaken at two epochs, between which Holmberg IX X-1 exhibited both flux and strong spectral variability, increasing in luminosity from Lx = (1.90+/-0.03)e40 erg/s to Lx = (3.35+/-0.03)e40 erg/s. Neither epoch exhibits a spectrum consistent with emission from the standard low/hard accretion state seen in Galactic black hole binaries, that would have been expected if Holmberg IX X-1 harbors a truly massive black hole accreting at substantially sub-Eddington accretion rates. The NuSTAR data confirm that the curvature observed previously in the 3-10 keV bandpass does represent a true spectral cutoff. During each epoch, the spectrum appears to be dominated by two optically thick thermal components, likely associated with an accretion disk. The spectrum also shows some evidence for a non-thermal tail at the highest energies, which may further support this scenario. The available data allow for either of the two thermal components to dominate the spectral evolution, although both scenarios require highly non-standard behavior for thermal accretion disk emission.

We present Nustar 3-40 keV observations of the optically selected Type 2 quasar (QSO2) SDSS J1034+6001 or Mrk 34. The high-quality hard X-ray spectrum and archival XMM-Newton data can be fitted self-consistently with a reflection-dominated continuum and strong Fe Kalpha fluorescence line with equivalent-width >1 keV. Prior X-ray spectral fitting below 10 keV showed the source to be consistent with being obscured by Compton-thin column densities of gas along the line-of-sight, despite evidence for much higher columns from multiwavelength data. NuSTAR now enables a direct measurement of this column, and shows that Nh lies in the Compton-thick (CT) regime. The new data also show a high intrinsic 2-10 keV luminosity of L_{2-10}~10^{44} erg/s, in contrast to previous low-energy X-ray measurements for which L_{2-10}<~10^{43} erg/s (i.e. X-ray selection below 10 keV does not pick up this source as an intrinsically luminous obscured quasar). Both the obscuring column and the intrinsic power are about an order of magnitude (or more) larger than inferred from pre-NuSTAR X-ray spectral fitting. Mrk34 is thus a 'gold standard' CT QSO2 and is the nearest non-merging system in this class, in contrast to the other local CT quasar NGC6240 which is currently undergoing a major merger coupled with strong star-formation. For typical X-ray bolometric correction factors, the accretion luminosity of Mrk34 is high enough to potentially power the total infrared luminosity. X-ray spectral fitting also shows that thermal emission related to star-formation is unlikely to drive the observed bright soft component below ~3 keV, favoring photionization instead.

We present direct constraints on how the formation of low-mass X-ray binary (LMXB) populations in galactic fields depends on stellar age. In this pilot study, we utilize Chandra and Hubble Space Telescope (HST) data to detect and characterize the X-ray point source populations of three nearby early-type galaxies: NGC 3115, 3379, and 3384. The luminosity-weighted stellar ages of our sample span 3-10 Gyr. X-ray binary population synthesis models predict that the field LMXBs associated with younger stellar populations should be more numerous and luminous per unit stellar mass than older populations due to the evolution of LMXB donor star masses. Crucially, the combination of deep Chandra and HST observations allows us to test directly this prediction by identifying and removing counterparts to X-ray point sources that are unrelated to the field LMXB populations, including LMXBs that are formed dynamically in globular clusters, Galactic stars, and background AGN/galaxies. We find that the "young" early-type galaxy NGC 3384 (~2-5 Gyr) has an excess of luminous field LMXBs (L_X > (5-10) x 10^37 erg/s) per unit K-band luminosity (L_K; a proxy for stellar mass) than the "old" early-type galaxies NGC 3115 and 3379 (~8-10 Gyr), which results in a factor of ~2-3 excess of LX/LK for NGC 3384. This result is consistent with the X-ray binary population synthesis model predictions; however, our small galaxy sample size does not allow us to draw definitive conclusions on the evolution field LMXBs in general. We discuss how future surveys of larger galaxy samples that combine deep Chandra and HST data could provide a powerful new benchmark for calibrating X-ray binary population synthesis models.

Using self-consistent, physically motivated models, we investigate the X-ray obscuration in 19 Type 2 [OIII] 5007 \AA\ selected AGN, 9 of which are local Seyfert 2 galaxies and 10 of which are Type 2 quasar candidates. We derive reliable line-of-sight and global column densities for these objects, which is the first time this has been reported for an AGN sample; 4 AGN have significantly different global and line-of-sight column densities. Five sources are heavily obscured to Compton-thick. We comment on interesting sources revealed by our spectral modeling, including a candidate ``naked'' Sy2. After correcting for absorption, we find that the ratio of the rest-frame, 2-10 keV luminosity (L$_{\rm 2-10keV,in}$) to L$_{\rm [OIII]}$ is 1.54 $\pm$ 0.49 dex which is essentially identical to the mean Type 1 AGN value. The Fe K$\alpha$ luminosity is significantly correlated with L$_{\rm [OIII]}$, but with substantial scatter. Finally, we do not find a trend between L$_{\rm 2-10keV,in}$ and global or line-of-sight column density, between column density and redshift, between column density and scattering fraction or between scattering fraction and redshift.

We present high-energy (3--30 keV) {\it NuSTAR} observations of the nearest quasar, the ultraluminous infrared galaxy (ULIRG) Markarian 231 (Mrk 231), supplemented with new and simultaneous low-energy (0.5--8 keV) data from {\it Chandra}. The source was detected, though at much fainter levels than previously reported, likely due to contamination in the large apertures of previous non-focusing hard X-ray telescopes. The full band (0.5--30 keV) X-ray spectrum suggests the active galactic nucleus (AGN) in Mrk 231 is absorbed by a patchy and Compton-thin (N$_{\rm H} \sim1.2^{+0.3}_{-0.3}\times10^{23}$ cm$^{-2}$) column. The intrinsic X-ray luminosity (L$_{\rm 0.5-30 keV}\sim1.0\times10^{43}$ erg s$^{1}$) is extremely weak relative to the bolometric luminosity where the 2--10 keV to bolometric luminosity ratio is $\sim$0.03% compared to the typical values of 2--15%. Additionally, Mrk 231 has a low X-ray-to-optical power law slope ($\alpha_{\rm OX}\sim-1.7$). It is a local example of a low-ionization broad absorption line (LoBAL) quasar that is intrinsically X-ray weak. The weak ionizing continuum may explain the lack of mid-infrared [O IV], [Ne V], and [Ne VI] fine-structure emission lines which are present in sources with otherwise similar AGN properties. We argue that the intrinsic X-ray weakness may be a result of the super-Eddington accretion occurring in the nucleus of this ULIRG, and may also be naturally related to the powerful wind event seen in Mrk 231, a merger remnant escaping from its dusty cocoon.

We have discovered a persistent, but highly variable X-ray source in the nearby starburst galaxy NGC 253. The source varies at the level of a factor of about 5 in count rate on timescales of a few hours. Two long observations of the source with Chandra and XMM-Newton show suggestive evidence for the source having a period of about 14-15 hours, but the time sampling in existing data is insufficient to allow a firm determination that the source is periodic. Given the amplitude of variation and the location in a nuclear starburst, the source is likely to be a Wolf-Rayet X-ray binary, with the tentative period being the orbital period of the system. In light of the fact that we have demonstrated that careful examination of the variability of moderately bright X-ray sources in nearby galaxies can turn up candidate Wolf-Rayet X-ray binaries, we discuss the implications of Wolf-Rayet X-ray binaries for predictions of the gravitational wave source event rate, and, potentially, interpretations of the events.

We present the results of NuSTAR and XMM-Newton observations of the two ultraluminous X-ray sources (ULX) NGC 1313 X-1 and X-2. The combined spectral bandpass of the two satellites enables us to produce the first spectrum of X-1 between 0.3 and 30 keV, while X-2 is not significantly detected by NuSTAR above 10 keV. The NuSTAR data demonstrate that X-1 has a clear cutoff above 10 keV, whose presence was only marginally detectable with previous X-ray observations. This cutoff rules out the interpretation of X-1 as a black hole in a standard low/hard state, and it is deeper than predicted for the downturn of a broadened iron line in a reflection-dominated regime. The cutoff differs from the prediction of a single-temperature Comptonization model. Further, a cold disk-like black body component at ~0.3 keV is required by the data, confirming previous measurements by XMM-Newton only. We observe a spectral transition in X-2, from a state with high luminosity and strong variability to a lower-luminosity state with no detectable variability, and we link this behavior to a transition from a super-Eddington to a sub-Eddington regime.

Following a serendipitous detection with the NuSTAR observatory, we present a multi-epoch spectral and temporal analysis of an extreme ultraluminous X-ray source (ULX) located in the outskirts of the Circinus galaxy, hereafter Circinus ULX5, including coordinated follow-up observations with XMM-Newton and NuSTAR. The NuSTAR data presented here represent one of the first instances of a ULX reliably detected at hard (E>10 keV) X-rays. Circinus ULX5 is variable on long timescales by at least a factor of ~5 in flux, and was caught in a historically bright state during our 2013 observations (0.3-30.0 keV luminosity of 1.6e40 erg/s). During this epoch, the source displayed a curved 3-10 keV spectrum, broadly similar to other bright ULXs. While in many cases the interpretation of spectral components in ULXs is uncertain, the spectral and temporal properties of the all the high quality datasets currently available strongly support a simple disk-corona model reminiscent of that invoked for Galactic binaries, with the accretion disk becoming more prominent as the luminosity increases. However, although the disk temperature and luminosity are remarkably well correlated across all timescales currently probed, the observed luminosity follows T^(1.70+/-0.17), flatter than that expected for simple blackbody radiation. The spectral variability displayed by Circinus ULX5 is highly reminiscent of that observed from some known Galactic black hole binaries (BHBs) at high luminosities. This comparison would imply a black hole mass of ~90 Msun for Circinus ULX5. However, given the diverse behavior observed from Galactic BHB accretion disks, this mass estimate is still uncertain. Finally, the limits placed on any undetected iron absorption features with the 2013 dataset imply that we are not viewing the central regions of Circinus ULX5 through any extreme super-Eddington outflow. [abridged]

(Abridged) We analyzed the {\it Chandra} and {\it XMM-Newton} archival observations for 72 type 2 quasars at $z<1$. These objects were was selected based on the [O III]$\lambda$5007 optical emission line which we assume to be an approximate indicator of the intrinsic AGN luminosity. We find that the means of the column density and photon index of our sample are $\log N_{\rm H}=23.0$ cm$^{-2}$ and $\Gamma=1.87$ respectively, which are consistent with results from deep X-ray surveys. The observed ratios of hard X-ray and [O III] line luminosities imply that the majority of our sample suffer significant amounts of obscuration in the hard X-ray band. A more physically realistic model which accounts for both Compton scattering and a potential partial covering of the central X-ray source was used to estimate the true absorbing column density. We find that the absorbing column density estimates based on simple power-law models significantly underestimate the actual absorption in approximately half of the sources. Eleven sources show a prominent Fe K$\alpha$ emission line, and we detect this line in the other sources through a joint fit (spectral stacking). The correlation between the Fe K$\alpha$ and [O III] fluxes and the inverse correlation of the equivalent width of Fe K$\alpha$ line with the ratio of hard X-ray and [O III] fluxes is consistent with previous results for lower luminosity Seyfert 2 galaxies. We conclude that obscuration is the cause of the weak hard X-ray emission rather than intrinsically low X-ray luminosities. We find that about half of the population of optically-selected type 2 quasars are likely to be Compton-thick. We also find no evidence that the amount of X-ray obscuration depends on the AGN luminosity (over a range of more than three orders-of-magnitude in luminosity).

We present the first public version (v0.2) of the open-source and community-developed Python package, Astropy. This package provides core astronomy-related functionality to the community, including support for domain-specific file formats such as Flexible Image Transport System (FITS) files, Virtual Observatory (VO) tables, and common ASCII table formats, unit and physical quantity conversions, physical constants specific to astronomy, celestial coordinate and time transformations, world coordinate system (WCS) support, generalized containers for representing gridded as well as tabular data, and a framework for cosmological transformations and conversions. Significant functionality is under active development, such as a model fitting framework, VO client and server tools, and aperture and point spread function (PSF) photometry tools. The core development team is actively making additions and enhancements to the current code base, and we encourage anyone interested to participate in the development of future Astropy versions.

Advanced LIGO and Advanced Virgo will be all-sky monitors for merging compact objects within a few hundred Mpc. Finding the electromagnetic counterparts to these events will require an understanding of the transient sky at low red-shift (z<0.1). We performed a systematic search for extragalactic, low red-shift, transient events in the XMM-Newton Slew Survey. In a flux limited sample, we found that highly-variable objects comprised 10% of the sample, and that of these, 10% were spatially coincident with cataloged optical galaxies. This led to 4x10^-4 transients per square degree above a flux threshold of 3x10^-12 erg cm-2 s-1 [0.2-2 keV] which might be confused with LIGO/Virgo counterparts. This represents the first extragalactic measurement of the soft X-ray transient rate within the Advanced LIGO/Virgo horizon. Our search revealed six objects that were spatially coincident with previously cataloged galaxies, lacked evidence for optical AGNs, displayed high luminosities around 10^43 erg s-1, and varied in flux by more than a factor of ten when compared with the ROSAT All-Sky Survey. At least four of these displayed properties consistent with previously observed tidal disruption events.

We present the largest-scale comparison to date between observed extragalactic X-ray binary (XRB) populations and theoretical models of their production. We construct observational X-ray luminosity functions (oXLFs) using Chandra observations of 12 late-type galaxies from the Spitzer Infrared Nearby Galaxy Survey (SINGS). For each galaxy, we obtain theoretical XLFs (tXLFs) by combining XRB synthetic models, constructed with the population synthesis code StarTrack, with observational star formation histories (SFHs). We identify highest-likelihood models both for individual galaxies and globally, averaged over the full galaxy sample. Individual tXLFs successfully reproduce about half of oXLFs, but for some galaxies we are unable to find underlying source populations, indicating that galaxy SFHs and metallicities are not well matched and/or XRB modeling requires calibration on larger observational samples. Given these limitations, we find that best models are consistent with a product of common envelope ejection efficiency and central donor concentration ~=0.1, and a 50% uniform -- 50% "twins" initial mass-ratio distribution. We present and discuss constituent subpopulations of tXLFs according to donor, accretor and stellar population characteristics. The galaxy-wide X-ray luminosity due to low-mass and high-mass XRBs, estimated via our best global model tXLF, follows the general trend expected from the L_x - star formation rate and L_x - stellar mass relations of Lehmer et al 2010. Our best models are also in agreement with modeling of the evolution both of XRBs over cosmic time and of the galaxy X-ray luminosity with redshift.

In this paper, we study the relationship between the 2-10 keV X-ray luminosity (L_X), assumed to originate from X-ray binaries (XRBs), and star formation rate (SFR) in UV-selected z<0.1 Lyman break analogs (LBAs). We present Chandra observations for four new GALEX-selected LBAs. Including previously studied LBAs, Haro 11 and VV 114, we find that LBAs demonstrate L_X/SFR ratios that are elevated by ~1.5sigma compared to local galaxies, similar to the ratios found for stacked LBGs in the early Universe (z>2). We show that these LBAs are unlikely to harbor AGN, based on their optical and X-ray spectra and the spatial distribution of the X-rays in three spatially extended cases. We expect that high-mass X-ray binaries (HMXBs) dominate the X-ray emission in these galaxies, based on their high specific SFRs (sSFRs=SFR/M* > 10^{-9}/yr), which suggest the prevalence of young stellar populations. Since both LBAs and LBGs have lower dust attenuations and metallicities compared to similar samples of more typical local galaxies, we investigate the effects of dust extinction and metallicity on the L_X/SFR for the broader population of galaxies with high sSFRs (>10^{-10}/yr). The estimated dust extinctions (corresponding to column densities of N_H<10^{22}/cm^2) are expected to have insignificant effects on observed L_X/SFR ratio for the majority of galaxy samples. We find that the observed relationship between L_X/SFR and metallicity appears consistent with theoretical expectations from X-ray binary population synthesis models. Therefore, we conclude that lower metallicities, related to more luminous HMXBs such as ultraluminous X-ray sources (ULXs), drive the elevated L_X/SFR observed in our sample of z<0.1 LBAs. The relatively metal-poor, active mode of star formation in LBAs and distant z>2 LBGs may yield higher total HMXB luminosity than found in typical galaxies in the local Universe.

We present results from three nearly simultaneous NuSTAR and Chandra monitoring observations between 2012 Sep 2 and 2012 Nov 16 of local star-forming galaxy NGC 253. The 3-40 keV NuSTAR intensity of the inner 20 arcsec (~400 pc) nuclear region varied by a factor of ~2 across the three monitoring observations. The Chandra data reveal that the nuclear region contains three bright X-ray sources, including a luminous (L2-10 keV ~ few x 10^39 erg/s) point source ~1 arcsec from the dynamical center of the galaxy (within the 3sigma positional uncertainty of the dynamical center); this source drives the overall variability of the nuclear region at energies >3 keV. We make use of the variability to measure the spectra of this single hard X-ray source when it was in bright states. The spectra are well described by an absorbed (NH ~ 1.6 x 10^23 cm^-2) broken power-law model with spectral slopes and break energies that are typical of ultraluminous X-ray sources (ULXs), but not AGN. A previous Chandra observation in 2003 showed a hard X-ray point source of similar luminosity to the 2012 source that was also near the dynamical center (~0.4 arcsec); however, this source was offset from the 2012 source position by ~1 arcsec. We show that the probability of the 2003 and 2012 hard X-ray sources being unrelated is >>99.99% based on the Chandra spatial localizations. Interestingly, the Chandra spectrum of the 2003 source (3-8 keV) is shallower in slope than that of the 2012 hard X-ray source. Its proximity to the dynamical center and harder Chandra spectrum indicate that the 2003 source is a better AGN candidate than any of the sources detected in our 2012 campaign; however, we were unable to rule out a ULX nature for this source. Future NuSTAR and Chandra monitoring would be well equipped to break the degeneracy between the AGN and ULX nature of the 2003 source, if again caught in a high state.

Using the spectra of a sample of $\sim$28,000 nearby obscured active galaxies from Data Release 7 of the Sloan Digital Sky Survey (SDSS), we probe the connection between AGN activity and star formation over a range of radial scales in the host galaxy. We use the extinction-corrected luminosity of the [OIII] 5007 \AA\ line as a proxy of intrinsic AGN power and supermassive black hole (SMBH) accretion rate. The star formation rates (SFRs) are taken from the MPA-JHU value-added catalog and are measured through the 3$^{\prime\prime}$ SDSS aperture. We construct matched samples of galaxies covering a range in redshifts. With increasing redshift, the projected aperture size encompasses increasing amounts of the host galaxy. This allows us to trace the radial distribution of star-formation as a function of AGN luminosity. We find that the star formation becomes more centrally concentrated with increasing AGN luminosity and Eddington ratio. This implies that such circumnuclear star formation is associated with AGN activity, and that it increasingly dominates over omnipresent disk star formation at higher AGN luminosities, placing critical constraints on theoretical models that link host galaxy star formation and SMBH fueling. We parametrize this relationship and find that the star formation on radial scales $<$1.7 kpc, when including a constant disk component, has a sub-linear dependence on SMBH accretion rate: $SFR \propto \dot{M}^{0.36}$, suggesting that angular momentum transfer through the disk limits accretion efficiency rather than the supply from stellar mass loss.

Emission from X-ray binaries (XRBs) is a major component of the total X-ray luminosity of normal galaxies, so X-ray studies of high redshift galaxies allow us to probe the formation and evolution of X-ray binaries on very long timescales. In this paper, we present results from large-scale population synthesis models of binary populations in galaxies from z = 0 to 20. We use as input into our modeling the Millennium II Cosmological Simulation and the updated semi-analytic galaxy catalog by Guo et al. (2011) to self-consistently account for the star formation history (SFH) and metallicity evolution of each galaxy. We run a grid of 192 models, varying all the parameters known from previous studies to affect the evolution of XRBs. We use our models and observationally derived prescriptions for hot gas emission to create theoretical galaxy X-ray luminosity functions (XLFs) for several redshift bins. Models with low CE efficiencies, a 50% twins mass ratio distribution, a steeper IMF exponent, and high stellar wind mass loss rates best match observational results from Tzanavaris & Georgantopoulos (2008), though they significantly underproduce bright early-type and very bright (Lx > 10d41) late-type galaxies. These discrepancies are likely caused by uncertainties in hot gas emission and SFHs, AGN contamination, and a lack of dynamically formed Low-mass XRBs. In our highest likelihood models, we find that hot gas emission dominates the emission for most bright galaxies. We also find that the evolution of the normal galaxy X-ray luminosity density out to z = 4 is driven largely by XRBs in galaxies with X-ray luminosities between 10d40 and 10d41 erg/s.

The Nuclear Spectroscopic Telescope Array (NuSTAR) mission, launched on 13 June 2012, is the first focusing high-energy X-ray telescope in orbit. NuSTAR operates in the band from 3 -- 79 keV, extending the sensitivity of focusing far beyond the ~10 keV high-energy cutoff achieved by all previous X-ray satellites. The inherently low-background associated with concentrating the X-ray light enables NuSTAR to probe the hard X-ray sky with a more than one-hundred-fold improvement in sensitivity over the collimated or coded-mask instruments that have operated in this bandpass. Using its unprecedented combination of sensitivity, spatial and spectral resolution, NuSTAR will pursue five primary scientific objectives, and will also undertake a broad program of targeted observations. The observatory consists of two co-aligned grazing-incidence X-ray telescopes pointed at celestial targets by a three-axis stabilized spacecraft. Deployed into a 600 km, near-circular, 6degree inclination orbit, the Observatory has now completed commissioning, and is performing consistent with pre-launch expectations. NuSTAR is now executing its primary science mission, and with an expected orbit lifetime of ten years, we anticipate proposing a guest investigator program, to begin in Fall 2014.

High redshift galaxies permit the study of the formation and evolution of X-ray binary populations on cosmological timescales, probing a wide range of metallicities and star-formation rates. In this paper, we present results from a large scale population synthesis study that models the X-ray binary populations from the first galaxies of the universe until today. We use as input to our modeling the Millennium II Cosmological Simulation and the updated semi-analytic galaxy catalog by Guo et al. (2011) to self-consistently account for the star formation history and metallicity evolution of the universe. Our modeling, which is constrained by the observed X-ray properties of local galaxies, gives predictions about the global scaling of emission from X-ray binary populations with properties such as star-formation rate and stellar mass, and the evolution of these relations with redshift. Our simulations show that the X-ray luminosity density (X-ray luminosity per unit volume) from X-ray binaries in our Universe today is dominated by low-mass X-ray binaries, and it is only at z>2.5 that high-mass X-ray binaries become dominant. We also find that there is a delay of ~1.1 Gyr between the peak of X-ray emissivity from low-mass Xray binaries (at z~2.1) and the peak of star-formation rate density (at z~3.1). The peak of the X-ray luminosity from high-mass X-ray binaries (at z~3.9), happens ~0.8 Gyr before the peak of the star-formation rate density, which is due to the metallicity evolution of the Universe.

We present the results of the Suzaku observation of the Seyfert 2 galaxy NGC 4507. This source is one of the X-ray brightest Compton-thin Seyfert 2s and a candidate for a variable absorber. Suzaku caught NGC 4507 in a highly absorbed state characterised by a high column density (NH \sim8 x10^23 cm^-2), a strong reflected component (R\sim 1.9) and a high equivalent width Fe K alpha emission line (EW\sim 500 eV). The Fe K alpha emission line is unresolved at the resolution of the Suzaku CCDs (sigma < 30 eV or FWHM < 3000 km s^-1) and most likely originates in a distant absorber. The Fe K beta emission line is also clearly detected and its intensity is marginally higher than the theoretical value for low ionisation Fe. A comparison with previous observations performed with XMM-Newton and BeppoSAX reveals that the X-ray spectral curvature changes on a timescale of a few months. We analysed all these historical observations, with standard models as well as with a most recent model for a toroidal reprocessor and found that the main driver of the observed 2-10 keV spectral variability is a change of the line-of-sight obscuration, varying from \sim4x10^23 cm^-2 to \sim9 x 10^23 cm^-2. The primary continuum is also variable, although its photon index does not appear to vary, while the Fe K alpha line and reflection component are consistent with being constant across the observations. This suggests the presence of a rather constant reprocessor and that the observed line of sight NH variability is either due to a certain degree of clumpiness of the putative torus or due to the presence of a second clumpy absorber.

We have explored the interplay of star formation and AGN activity in soft X-rays (0.5-2 keV) in two samples of Seyfert 2 galaxies (Sy2s). Using a combination of low resolution CCD spectra from Chandra and XMM-Newton, we modeled the soft emission of 34 Sy2s using power law and thermal models. For the 11 sources with high signal-to-noise Chandra imaging of the diffuse host galaxy emission, we estimate the luminosity due to star formation by removing the AGN, fitting the residual emission. The AGN and star formation contributions to the soft X-ray luminosity (i.e. L$_{x,AGN}$ and L$_{x,SF}$) for the remaining 24 Sy2s were estimated from the power law and thermal luminosities derived from spectral fitting. These luminosities were scaled based on a template derived from XSINGS analysis of normal star forming galaxies. To account for errors in the luminosities derived from spectral fitting and the spread in the scaling factor, we estimated L$_{x,AGN}$ and L$_{x,SF}$ from Monte Carlo simulations. These simulated luminosities agree with L$_{x,AGN}$ and L$_{x,SF}$ derived from Chandra imaging analysis within a 3\sigma\ confidence level. Using the infrared [NeII]12.8\mu m and [OIV]26\mu m lines as a proxy of star formation and AGN activity, respectively, we independently disentangle the contributions of these two processes to the total soft X-ray emission. This decomposition generally agrees with L$_{x,SF}$ and L$_{x,AGN}$ at the 3\sigma\ level. In the absence of resolvable nuclear emission, our decomposition method provides a reasonable estimate of emission due to star formation in galaxies hosting type 2 AGN.

We study a combined sample of 264 star-forming, 51 composite, and 73 active galaxies using optical spectra from SDSS and mid-infrared (mid-IR) spectra from the Spitzer Infrared Spectrograph. We examine optical and mid-IR spectroscopic diagnostics that probe the amount of star formation and relative energetic contributions from star formation and an active galactic nucleus (AGN). Overall we find good agreement between optical and mid-IR diagnostics. Misclassifications of galaxies based on the SDSS spectra are rare despite the presence of dust obscuration. The luminosity of the [NeII] 12.8 \mu m emission-line is well correlated with the star formation rate (SFR) measured from the SDSS spectra, and this holds for the star forming, composite, and AGN-dominated systems. AGN show a clear excess of [NeIII] 15.6 \mu m emission relative to star forming and composite systems. We find good qualitative agreement between various parameters that probe the relative contributions of the AGN and star formation, including: the mid-IR spectral slope, the ratio of the [NeV] 14.3 \mu m to [NeII] \mu m 12.8 fluxes, the equivalent widths of the 7.7, 11.3, and 17 $\mu m$ PAH features, and the optical "D" parameter which measures the distance a source lies from the locus of star forming galaxies in the optical BPT emission-line diagnostic diagram. We also consider the behavior of the three individual PAH features by examining how their flux ratios depend upon the degree of AGN-dominance. We find that the PAH 11.3 \mu m feature is significantly suppressed in the most AGN-dominated systems.

We present a {\it Chandra} observation of IRAS 19254--7245, a nearby ULIRG also known as {\it the Superantennae}. The high spatial resolution of {\it Chandra} allows us to disentangle for the first time the diffuse starburst emission from the embedded Compton-thick AGN. The 2-10 keV spectrum of the AGN emission is fitted by a flat power-law $\Gamma=1.3$) and a He-like Fe K$\alpha$ line with EW$\sim$1.5 keV, consistent with previous observations. The Fe K$\alpha$ line profile could be resolved as a blend of a neutral 6.4 keV line and an ionized 6.7 keV (He-like) or 6.9 keV (H-like) line. Variability is detected compared with the previous {\it XMM-Newton} and {\it suzaku} observations, demonstrating the compact size of the iron line emission. We fit the spectrum of the galaxy-scale extended emission excluding the AGN and other bright point sources with a soft thermal component with kT~0.8 keV. The luminosity of the extended emission is about one order of magnitude lower than that of the AGN. The basic physical and structural properties of the extended emission are fully consistent with a galactic wind being driven by the starburst (no contribution to the feedback by the AGN is required). A candidate ultra-luminous X-ray source is detected 8\arcsec\ south of the southern nucleus. The 0.3-10 keV luminosity of this off-nuclear point source is ~$6\times 10^{40}$ erg s$^{-1}$ if the emission is isotropic and the source is associated with the Superantennae.

We present 0.5-2 keV, 2-8 keV, 4-8 keV, and 0.5-8 keV cumulative and differential number counts (logN-logS) measurements for the recently completed ~4 Ms Chandra Deep Field-South (CDF-S) survey, the deepest X-ray survey to date. We implement a new Bayesian approach, which allows reliable calculation of number counts down to flux limits that are factors of ~1.9-4.3 times fainter than the previously deepest number-counts investigations. In the soft band, the most sensitive bandpass in our analysis, the ~4 Ms CDF-S reaches a maximum source density of ~27,800 deg-2. By virtue of the exquisite X-ray and multiwavelength data available in the CDF-S, we are able to measure the number counts from a variety of source populations (active galactic nuclei [AGNs], normal galaxies, and Galactic stars) and subpopulations (as a function of redshift, AGN absorption, luminosity, and galaxy morphology), and test models that describe their evolution. We find that AGNs still dominate the X-ray number counts down to the faintest flux levels for all bands and reach a limiting soft-band source density of ~14,900 deg-2, the highest reliable AGN source density measured at any wavelength. We find that the normal-galaxy counts rise rapidly near the flux limits, and at the limiting soft-band flux, reach source densities of ~12,700 deg-2 and make up 46+/-5% of the total number counts. The rapid rise of the galaxy counts toward faint fluxes, and significant normal-galaxy contributions to the overall number counts, indicate that normal galaxies will overtake AGNs just below the ~4 Ms soft-band flux limit and will provide a numerically significant new X-ray source population in future surveys that reach below the ~4 Ms sensitivity limit. We show that a future ~10 Ms CDF-S would allow for a significant increase in X-ray detected sources, with many of the new sources being cosmologically distant (z > 0.6) normal galaxies.

We have used the European VLBI Network (EVN) to observe a sample of Lyman Break Analogs (LBAs), nearby (z < 0.3) galaxies with properties similar to the more distant Lyman Break Galaxies (LBGs). The study of LBGs may help define the feed-back relationship between black holes (BHs) and their host galaxies. Previous VLA observations have shown that the kpc-scale radio emission from LBAs is dominated by starbursts. The main targets of this VLBI experiment were selected because they possessed emission-line properties between starbursts and Type 2 (obscured) AGN. Eight targets (three star-forming LBAs, four composite LBAs, and one Type 1 AGN) were observed at 5 GHz, four of which were also observed at 1.7 GHz. One star-forming LBA and one composite LBA were detected above 5 \sigma at 1.7 GHz (only), while the AGN was detected at 5 GHz. In both LBAs, the radio luminosity (LR) exceeds that expected from supernovae (remnants) based on a comparison with Arp220, Arp229A and Mrk273, by factors of 2 - 8. The composite LBA exhibits a compact core emitting around 10% of the VLA flux density. The high Tb of 3.5E7 K and excess core L_R with respect to the L_R/L_X relation of radio-quiet AGN indicate that this LBA possesses an obscured AGN (MBH ~ E5-E7 M_sun). While weak AGN may co-exist with the starbursts as shown in at least one of the LBAs, their contribution to the total radio flux is fairly minimal. Our results show that the detection of such weak AGN presents a challenge at radio, X-ray and optical emission-line wavelengths at z ~ 0.2, indicating the great difficulties that need to be overcome in order to study similar processes at high redshift when these types of galaxies were common.

A non-proprietary Mission Concept available for presentation to NASA, providing High Area, High Resolution Imaging Spectroscopy and Timing with Arcmin Angular Resolution Submitted in response to NASA 2011 RFI NNH11ZDA018L 'Concepts for the Next NASA X-ray Astronomy Mission'

We report the discovery of a Compton-thick Active Galactic Nucleus (AGN) at z=4.76 in the 4 Ms Chandra Deep Field South. This object was selected as a V-band dropout in HST/ACS images and previously recognized as an AGN from optical spectroscopy. The 4 Ms Chandra observations show a significant (~4.2sigma) X-ray detection at the V-band dropout position. The X-ray source displays a hardness ratio of HR=0.23+-0.24, which, for a source at z~5, is highly suggestive of Compton-thick absorption. The source X-ray spectrum is seen above the background level in the energy range of ~0.9-4 keV, i.e., in the rest-frame energy range of ~5-23 keV. When fixing the photon index to Gamma=1.8, the measured column density is N_H=1.4^{+0.9}_{-0.5} x 10^{24} cm^{-2}, which is Compton-thick. To our knowledge, this is the most distant heavily obscured AGN, confirmed by X-ray spectral analysis, discovered so far. The intrinsic (de-absorbed), rest-frame luminosity in the 2-10 keV band is ~2.5 x 10^{44} erg s^{-1}, which places this object among type-2 quasars. The Spectral Energy Distribution shows that massive star formation is associated with obscured black hole accretion. This system may have then been caught during a major co-eval episode of black hole and stellar mass assembly at early times. The measure of the number density of heavily obscured AGN at high redshifts will be crucial to reconstruct the BH/galaxy evolution history from the beginning.

We have used XMM-Newton to observe six Lyman Break Analogs (LBAs): members of the rare population of local galaxies that have properties that are very similar to distant Lyman Break Galaxies. Our six targets were specifically selected because they have optical emission-line properties that are intermediate between starbursts and Type 2 (obscured) AGN. Our new X-ray data provide an important diagnostic of the presence of an AGN. We find X-ray luminosities of order 10^{42} erg/s and ratios of X-ray to far-IR luminosities that are higher than values in pure starburst galaxies by factors ranging from ~ 3 to 30. This strongly suggests the presence of an AGN in at least some of the galaxies. The ratios of the luminosities of the hard (2-10 keV) X-ray to [O III]\lambda 5007 emission-line are low by about an order-of-magnitude compared to Type 1 AGN, but are consistent with the broad range seen in Type 2 AGN. Either the AGN hard X-rays are significantly obscured or the [O III] emission is dominated by the starburst. We searched for an iron emission line at ~ 6.4 keV, which is a key feature of obscured AGN, but only detected emission at the ~ 2\sigma level. Finally, we find that the ratios of the mid-infrared (24\mu m) continuum to [O III]\lambda 5007 luminosities in these LBAs are higher than the values for Type 2 AGN by an average of 0.8 dex. Combining all these clues, we conclude that an AGN is likely to be present, but that the bolometric luminosity is produced primarily by an intense starburst. If these black holes are radiating at the Eddington limit, their masses would lie in the range of 10^5 to 10^6 M_{sun}. These objects may offer ideal local laboratories to investigate the processes by which black holes grew in the early universe.

We have analyzed archival \textit{Chandra} and \textit{XMM-Newton} data for two nearly complete homogeneously selected samples of type 2 Seyfert galaxies (Sy2s). These samples were selected based on intrinsic Active Galactic Nuclei (AGN) flux proxies: a mid-infrared (MIR) sample from the original IRAS 12$\mu$m survey and an optical ([OIII]$\lambda$ 5007 \AA\ flux limited) sample from the Sloan Digital Sky Survey (SDSS), providing a total of 45 Sy2s. As the MIR and [OIII] fluxes are largely unaffected by AGN obscuration, these samples can present an unbiased estimate of the Compton-thick (column density N$_H > 10^{24}$ cm$^{-2}$) subpopulation. We find that the majority of this combined sample is likely heavily obscured, as evidenced by the 2-10 keV X-ray attenuation (normalized by intrinsic flux diagnostics) and the large Fe K$\alpha$ equivalent widths (several hundred eV to over 1 keV). A wide range of these obscuration diagnostics is present, showing a continuum of column densities, rather than a clear segregation into Compton-thick and Compton-thin sub-populations. We find that in several instances, the fitted column densities severely under-represent the attenuation implied by these obscuration diagnostics, indicating that simple X-ray models may not always recover the intrinsic absorption. We compared AGN and host galaxy properties, such as intrinsic luminosity, central black hole mass, accretion rate, and star formation rate with obscuration diagnostics. No convincing evidence exists to link obscured sources with unique host galaxy populations from their less absorbed counterparts. Finally, we estimate that a majority of these Seyfert 2s will be detectable in the 10-40 keV range by the future NuSTAR mission, which would confirm whether these heavily absorbed sources are indeed Compton-thick.

The Wide Field X-Ray Telescope (WFXT) is a medium-class mission designed to be 2-orders-of-magnitude more sensitive than any previous or planned X-ray mission for large area surveys and to match in sensitivity the next generation of wide-area optical, IR and radio surveys. Using an innovative wide-field X-ray optics design, WFXT provides a field of view of 1 square degree (10 times Chandra) with an angular resolution of 5" (Half Energy Width, HEW) nearly constant over the entire field of view, and a large collecting area (up to 1 m^2 at 1 keV, > 10x Chandra) over the 0.1-7 keV band. WFXTs low-Earth orbit also minimizes the particle background. In five years of operation, WFXT will carry out three extragalactic surveys at unprecedented depth and address outstanding questions in astrophysics, cosmology and fundamental physics. In this article, we illustrate the mission concept and the connection between science requirements and mission parameters.

We investigate the scientific impact of the Wide Field X-ray Telescope mission. We present simulated images and spectra of X-ray sources as observed from the three surveys planned for the nominal 5-year WFXT lifetime. The goal of these simulations is to provide WFXT images of the extragalactic sky in different energy bands based on accurate description of AGN populations, normal and star forming galaxies, groups and clusters of galaxies. The images are realized using a detailed PSF model, instrumental and physical backgrounds/foregrounds, accurate model of the effective area and the related vignetting effect. Thanks to this comprehensive modelization of the WFXT properties, the simulated images can be used to evaluate the flux limits for detection of point and extended sources, the effect of source confusion at very faint fluxes, and in general the efficiency of detection algorithms. We also simulate the spectra of the detected sources, in order to address specific science topics which are unique to WFXT. Among them, we focus on the characterization of the Intra Cluster Medium (ICM) of high-z clusters, and in particular on the measurement of the redshift from the ICM spectrum in order to build a cosmological sample of galaxy clusters. The end-to-end simulation procedure presented here, is a valuable tool in optimizing the mission design. Therefore, these simulations can be used to reliably characterize the WFXT discovery space and to verify the connection between mission requirements and scientific goals. Thanks to this effort, we can conclude on firm basis that an X-ray mission optimized for surveys like WFXT is necessary to bring X-ray astronomy at the level of the optical, IR, submm and radio wavebands as foreseen in the coming decade.

We present new Chandra observations that complete a sample of seventeen (17) luminous infrared galaxies (LIRGs) with D < 60 Mpc and low Galactic column densities of N_H < 5 X 10^20 cm^-2. The LIRGs in our sample have total infrared (8-1000um) luminosities in the range of L_IR ~ (1-8) X 10^11 L_sol. The high-resolution imaging and X-ray spectral information from our Chandra observations allow us to measure separately X-ray contributions from active galactic nuclei (AGNs) and normal galaxy processes (e.g., X-ray binaries and hot gas). We utilized total infrared plus UV luminosities to estimate star-formation rates (SFRs) and K-band luminosities and optical colors to estimate stellar masses (M*) for the sample. Under the assumption that the galaxy-wide 2-10 keV luminosity (LX) traces the combined emission from high mass X-ray binaries (HMXBs) and low mass X-ray binaries (LMXBs), and that the power output from these components are linearly correlated with SFR and M*, respectively, we constrain the relation LX = alpha M* + beta SFR. To achieve this, we construct a Chandra-based data set composed of our new LIRG sample combined with additional samples of less actively star-forming normal galaxies and more powerful LIRGs and ultraluminous infrared galaxies (ULIRGs) from the literature. Using these data, we measure best-fit values of alpha = (9.05 +/- 0.37) X 10^28 ergs s^-1 Msol^-1 and beta = (1.62 +/- 0.22) X 10^39 ergs s^-1 (Msol yr^-1)^-1. This scaling provides a more physically meaningful estimate of LX, with ~0.1-0.2 dex less scatter, than a direct linear scaling with SFR (abridged).

We use WIRC, IR images of the Antennae (NGC 4038/4039) together with the extensive catalogue of 120 X-ray point sources (Zezas et al. 2006) to search for counterpart candidates. Using our proven frame-tie technique, we find 38 X-ray sources with IR counterparts, almost doubling the number of IR counterparts to X-ray sources first identified in Clark et al. (2007). In our photometric analysis, we consider the 35 IR counterparts that are confirmed star clusters. We show that the clusters with X-ray sources tend to be brighter, K_s ~16 mag, with (J-K_s) = 1.1 mag. We then use archival HST images of the Antennae to search for optical counterparts to the X-ray point sources. We employ our previous IR-to-X-ray frame-tie as an intermediary to establish a precise optical-to-X-ray frame-tie with <0.6 arcsec rms positional uncertainty. Due to the high optical source density near the X-ray sources, we determine that we cannot reliably identify counterparts. Comparing the HST positions to the 35 identified IR star cluster counterparts, we find optical matches for 27 of these sources. Using Bruzual-Charlot spectral evolutionary models, we find that most clusters associated with an X-ray source are massive, ~10^6 M_sun, young, ~10^6 yr, with moderate metallicities, Z=0.05.

We consider five indicators for intrinsic AGN luminosity: the luminosities of the [OIII]$\lambda$5007 line, the [OIV]25.89$\mu$m line, the mid-infrared (MIR) continuum emission by the torus, and the radio and hard X-ray (E $>$ 10keV) continuum emission. We compare these different proxies using two complete samples of low-redshift type 2 AGN selected in a homogeneous way based on different indicators: an optically selected [OIII] sample and a mid-infrared selected 12$\mu$m sample. We examine the correlations between all five different proxies, and find better agreement for the [OIV], MIR, and [OIII] luminosities than for the hard X-ray and radio luminosities. Next, we compare the ratios of the fluxes of the different proxies to their values in unobscured Type 1 AGN. The agreement is best for the ratio of the [OIV] and MIR fluxes, while the ratios of the hard X-ray to [OIII], [OIV], and MIR fluxes are systematically low by about an order-of-magnitude in the Type 2 AGN, indicating that hard X-ray selected samples do not represent the full Type 2 AGN population. In a similar spirit, we compare different optical and MIR diagnostics of the relative energetic contributions of AGN and star formation processes in our samples of Type 2 AGN. We find good agreement between the various diagnostic parameters, such as the equivalent width of the MIR polycyclic aromatic hydrocarbon features, the ratio of the MIR [OIV]/[NeII] emission-lines, the spectral index of the MIR continuum, and the commonly used optical emission-line ratios. Finally, we test whether the presence of cold gas associated with star-formation leads to an enhanced conversion efficiency of AGN ionizing radiation into [OIII] or [OIV] emission. We find that no compelling evidence exists for this scenario for the luminosities represented in this sample (L$_{bol}$ $\approx$ 10$^{9}$ - 8 $\times$ 10$^{11}$ L$_{\sun}$). (abridged)

We present HST UV/optical imaging, Spitzer mid-IR photometry, and optical spectroscopy of a sample of 30 low-redshift (z=0.1-0.3) galaxies chosen from SDSS/GALEX to be accurate local analogs of the high-z Lyman Break Galaxies. The Lyman Break Analogs (LBAs) are similar in mass, metallicity, dust, SFR, size and gas velocity dispersion, thus enabling a detailed investigation of processes that are important at high-z. The optical emission line properties of LBAs are also similar to those of LBGs, indicating comparable conditions in their ISM. In the UV, LBAs are characterized by complexes of massive star-forming "clumps", while in the optical they most often show evidence for (post-)mergers/interactions. In 6 cases, we find an extremely massive (>10^9 Msun) compact (R~100 pc) dominant central object (DCO). The DCOs are preferentially found in LBAs with the highest mid-IR luminosities and correspondingly high SFRs (15-100 Msun/yr). We show that the massive SF clumps (including the DCOs) have masses much larger than the nuclear super star clusters seen in normal late type galaxies. However, the DCOs have masses, sizes, and densities similar to the excess-light/central-cusps seen in typical elliptical galaxies with masses similar to the LBA galaxies. We suggest that the DCOs form in present-day examples of the dissipative mergers at high redshift that are believed to have produced the central-cusps in local ellipticals. More generally, the properties of the LBAs are consistent with the idea that instabilities in a gas-rich disk lead to very massive star-forming clumps that eventually coalesce to form a spheroid. We speculate that the DCOs are too young at present to be growing a supermassive black hole because they are still in a supernova-dominated outflow phase.

(abridged)The majority of Active Galactic Nuclei (AGN) suffer from significant obscuration by surrounding dust and gas. X-ray surveys in the 2-10 keV band will miss the most heavily-obscured AGN in which the absorbing column density exceeds $\sim10^{24}$cm$^{-2}$ (the Compton-thick AGN). It is therefore vital to know the fraction of AGN that are missed in such X-rays surveys and to determine if these AGN represent some distinct population in terms of the fundamental properties of AGN and/or their host galaxies. In this paper we present the analysis of \textit{XMM-Newton} X-ray data for a complete sample of 17 low-redshift Type 2 Seyfert galaxies chosen from the Sloan Digital Sky Survey based solely on the high observed flux of the [OIII]$\lambda$5007 emission-line. This line is formed in the Narrow Line Region hundreds of parsecs away from the central engine. Thus, unlike the X-ray emission, it is not affected by obscuration due to the torus surrounding the black hole. It therefore provides a useful isotropic indicator of the AGN luminosity. As additional indicators of the intrinsic AGN luminosity, we use the Spitzer Space Telescope to measure the luminosities of the mid-infrared continuum and the [OIV]25.89$\mu$m narrow emission-line. We then use the ratio of the 2-10 keV X-ray luminosity to the [OIII], [OIV], and mid-infrared luminosities to assess the amount of X-ray obscuration and to distinguish between Compton-thick and Compton-thin objects. We find that the majority of the sources suffer significant amounts of obscuration: the observed 2-10 keV emission is depressed by more than an order-of-magnitude in 11 of the 17 cases (as expected for Compton-thick sources).

We discuss here a long Suzaku observation of IRAS 19254-7245 (also known as the Superantennae), one of the brightest and well studied Ultra Luminous Infrared Galaxies in the local Universe. This long observation provided the first detection of IRAS 19254-7245 above 10 keV, and measured a 15-30 keV flux of ~5x10^(-12) erg cm^-2 s^-1. The detection above 10 keV has allowed us to unveil, for the first time, the intrinsic luminosity of the AGN hosted in IRAS 19254-7245, which is strongly absorbed (Nh ~ 3x10^(24) cm^-2) and has an intrinsic luminosity in the QSO regime (L(2-10 keV) ~ 3 x 10^(44) erg s^-1). The 2-10 keV spectrum of IRAS 19254-7245 is remarkably hard (Gamma~1.2), and presents a strong iron line (EW ~0.7 keV), clearly suggesting that below 10 keV we are seeing only reprocessed radiation. Since the energy of the Fe K emission is found to be at ~6.7 keV, consistent with He-like Fe, and its EW is too high to be explained in a starburst dominated scenario, we suggest that the 2--10 keV emission of IRAS 19254-7245 is dominated by reflection/scattering from highly ionized matter. Indeed, within this latter scenario we found that the photon index of the illuminating source is Gamma=1.87 (+0.11,-0.28), in excellent agreement with the mean value found for radio quiet unobscured AGN.

We discuss the central role played by X-ray studies to reconstruct the past history of formation and evolution of supermassive Black Holes (BHs), and the role they played in shaping the properties of their host galaxies. We shortly review the progress in this field contributed by the current X-ray and multiwavelength surveys. Then, we focus on the outstanding scientific questions that have been opened by observations carried out in the last years and that represent the legacy of Chandra and XMM, as for X-ray observations, and the legacy of the SDSS, as for wide area surveys: 1) When and how did the first supermassive black holes form? 2) How does cosmic environment regulate nuclear activity (and star formation) across cosmic time? 3) What is the history of nuclear activity in a galaxy lifetime? We show that the most efficient observational strategy to address these questions is to carry out a large-area X-ray survey, reaching a sensitivity comparable to that of deep Chandra and XMM pointings, but extending over several thousands of square degrees. Such a survey can only be carried out with a Wide-Field X-ray Telescope (WFXT) with a high survey speed, due to the combination of large field of view and large effective area, i.e., grasp, and sharp PSF. We emphasize the important synergies that WFXT will have with a number of future groundbased and space telescopes, covering from the radio to the X-ray bands and discuss the immense legacy value that such a mission will have for extragalactic astronomy at large.

Sensitive, wide-area X-ray surveys which would be possible with the WFXT will detect huge samples of virialized objects spanning the mass range from sub-groups to the most massive clusters, and extending in redshift to beyond z=2. These samples will be an excellent dataset for carrying out many traditional cosmological tests using the cluster mass function and power spectrum. Uniquely, WFXT will be able not only to detect clusters but also to make detailed X-ray measurements for a large number of clusters and groups right from the survey data. Very high quality measurements of the cluster mass function and spatial correlation over a very wide range of masses, spatial scales, and redshifts, will be useful for expanding the cosmological discovery space, and in particular, in searching for departures from the "concordant" Lambda-CDM cosmological model. Finding such departures would have far-reaching implications on our understanding of the fundamental physics which governs the Universe.

We discuss the central role played by the X-ray study of hot baryons within galaxy clusters to reconstruct the assembly of cosmic structures and to trace the past history of star formation and accretion onto supermassive Black Holes (BHs). We shortly review the progress in this field contributed by the current generation of X-ray telescopes. Then, we focus on the outstanding scientific questions that have been opened by observations carried out in the last years and that represent the legacy of Chandra and XMM: (a) When and how is entropy injected into the inter-galactic medium (IGM)? (b) What is the history of metal enrichment of the IGM? (c) What physical mechanisms determine the presence of cool cores in galaxy clusters? (d) How is the appearance of proto-clusters at z~2 related to the peak of star formation activity and BH accretion? We show that a highly efficient observational strategy to address these questions is to carry out a large-area X-ray survey, reaching a sensitivity comparable to that of deep Chandra and XMM pointings, but extending over several thousands of square degrees. A similar survey can only be carried out with a Wide-Field X-ray Telescope (WFXT), which combines a high survey speed with a sharp PSF across the entire FoV. We emphasize the important synergies that WFXT will have with a number of future ground-based and space telescopes, covering from the radio to the X-ray bands. Finally, we discuss the immense legacy value that such a mission will have for extragalactic astronomy at large.

We explore the multidimensional, multiwavelength selection of quasars from mid-IR (MIR) plus optical data, specifically from Spitzer-IRAC and the Sloan Digital Sky Survey (SDSS). We apply modern statistical techniques to combined Spitzer MIR and SDSS optical data, allowing up to 8-D color selection of quasars. Using a Bayesian selection method, we catalog 5546 quasar candidates to an 8.0 um depth of 56 uJy over an area of ~24 sq. deg; ~70% of these candidates are not identified by applying the same Bayesian algorithm to 4-color SDSS optical data alone. Our selection recovers 97.7% of known type 1 quasars in this area and greatly improves the effectiveness of identifying 3.5<z<5 quasars. Even using only the two shortest wavelength IRAC bandpasses, it is possible to use our Bayesian techniques to select quasars with 97% completeness and as little as 10% contamination. This sample has a photometric redshift accuracy of 93.6% (Delta Z +/-0.3), remaining roughly constant when the two reddest MIR bands are excluded. While our methods are designed to find type 1 (unobscured) quasars, as many as 1200 of the objects are type 2 (obscured) quasar candidates. Coupling deep optical imaging data with deep mid-IR data could enable selection of quasars in significant numbers past the peak of the quasar luminosity function (QLF) to at least z~4. Such a sample would constrain the shape of the QLF and enable quasar clustering studies over the largest range of redshift and luminosity to date, yielding significant gains in our understanding of quasars and the evolution of galaxies.

While the exceptional sensitivity of Chandra and XMM-Newton has resulted in revolutionary studies of the Galactic neighborhood in the soft (<10 keV) X-ray band, there are many open questions. We discuss these issues and how they would be addressed by very wide-area (> 100 sq. deg.) X-ray surveys.

What is the contribution of mass, metals and energy from starburst galaxies to the Intergalactic Medium? Starburst galaxies drive galactic-scale outflows or "superwinds" that may be responsible for removing metals from galaxies and polluting the Intergalactic Medium (IGM). In the last decade tremendous progress was made in mapping cool entrained gas in superwinds through UV/optical imaging and absorption line spectroscopy. These studies demonstrated that superwinds are ubiquitous in galaxies forming stars at high surface densities and that the most powerful starbursts can drive outflows near escape velocity. Theoretical models of galaxy evolution have begun to incorporate superwinds, using various ad-hoc prescriptions based on our knowledge of the cool gas. However, these efforts are fundamentally impeded by our lack of information about the hot phase of these outflows. The hot X-ray emitting phase of a superwind contains the majority of its energy and newly-synthesized metals, and given its high specific energy and inefficient cooling it is also the component most likely escape from the galaxy's gravitational potential well. Knowledge of the chemical composition and velocity of the hot gas are crucial to assess the energy and chemical feedback from a starburst. A high priority for the next decade is to enable direct measurements to be made of the rates at which starburst galaxies of all masses eject gas, metals, and energy into the IGM. This will require a high sensitivity X-ray imaging spectrometer capable of measuring velocities in faint diffuse X-ray emission with a velocity accuracy of ~ 100 km/s. Such spectral resolution automatically allows detailed line-based plasma diagnostics, and thus composition, energetics and flow rates can be derived.

While often the point sources in X-ray surveys are dominated by AGN, with the high sensitivity of modern X-ray telescopes such as Chandra and XMM-Newton normal/starburst galaxies are also being detected in large numbers. We have made use of Bayesian statistics for both the selection of galaxies from deep X-ray surveys and in the analysis of the luminosity functions for galaxies. These techniques can be used to similarly select galaxies from wide-area X-ray surveys and to analyze their luminosity function. The prospects for detecting galaxies and AGN from a proposed ``wide-deep'' XMM-Newton survey and from future wide-area X-ray survey missions (such as WFXT and eRosita) are also discussed.

We use the previously-identified 15 infrared star-cluster counterparts to X-ray point sources in the interacting galaxies NGC 4038/4039 (the Antennae) to study the relationship between total cluster mass and X-ray binary number. This significant population of X-Ray/IR associations allows us to perform, for the first time, a statistical study of X-ray point sources and their environments. We define a quantity, \eta, relating the fraction of X-ray sources per unit mass as a function of cluster mass in the Antennae. We compute cluster mass by fitting spectral evolutionary models to K_s luminosity. Considering that this method depends on cluster age, we use four different age distributions to explore the effects of cluster age on the value of \eta and find it varies by less than a factor of four. We find a mean value of \eta for these different distributions of \eta = 1.7 x 10^-8 M_\sun^-1 with \sigma_\eta = 1.2 x 10^-8 M_\sun^-1. Performing a \chi^2 test, we demonstrate \eta could exhibit a positive slope, but that it depends on the assumed distribution in cluster ages. While the estimated uncertainties in \eta are factors of a few, we believe this is the first estimate made of this quantity to ``order of magnitude'' accuracy. We also compare our findings to theoretical models of open and globular cluster evolution, incorporating the X-ray binary fraction per cluster.

We present the results of the simultaneous deep XMM and Chandra observations of the bright Seyfert 1.9 galaxy MCG-5-23-16, which is thought to have one of the best known examples of a relativistically broadened iron K-alpha line. The time averaged spectral analysis shows that the iron K-shell complex is best modeled with an unresolved narrow emission component (FWHM < 5000 km/s, EW ~ 60 eV) plus a broad component. This latter component has FWHM ~ 44000 km/s and EW ~ 50 eV. Its profile is well described by an emission line originating from an accretion disk viewed with an inclination angle ~ 40^\circ and with the emission arising from within a few tens of gravitational radii of the central black hole. The time-resolved spectral analysis of the XMM EPIC-pn spectrum shows that both the narrow and broad components of the Fe K emission line appear to be constant in time within the errors. We detected a narrow sporadic absorption line at 7.7 keV which appears to be variable on a time-scale of 20 ksec. If associated with Fe XXVI Ly-alpha this absorption is indicative of a possibly variable, high ionization, high velocity outflow. The variability of this absorption feature appears to rule out a local (z=0) origin. The analysis of the XMM RGS spectrum reveals that the soft X-ray emission of MCG-5-23-16 is likely dominated by several emission lines superimposed on an unabsorbed scattered power-law continuum. The lack of strong Fe L shell emission together with the detection of a strong forbidden line in the O VII triplet is consistent with a scenario where the soft X-ray emission lines are produced in a plasma photoionized by the nuclear emission.

We obtained a wide-band spectrum of the Compton-thick Seyfert 2 galaxy Mrk 3 with Suzaku. The observed spectrum was clearly resolved into weak, soft power-law emission, a heavily absorbed power-law component, cold reflection, and many emission lines. The heavily absorbed component, absorbed by gas with a column density of 1.1x10^24 cm^-2, has an intrinsic 2--10 keV luminosity of ~1.6x10^43 erg s^-1, and is considered to be direct emission from the Mrk 3 nucleus. The reflection component was interpreted as reflection of the direct light off cold, thick material; the reflection fraction $R$ was 1.36+/-0.20. The cold material is inferred to be located > 1 pc from the central black hole of Mrk 3 due to the low ionization parameter of iron (xi < 1 erg cm s^-1) and the narrow iron line width (s < 22 eV). A Compton shoulder to the iron line was detected, but the intensity of the shoulder component was less than that expected from spherically distributed Compton-thick material. The weak, soft power-law emission is considered to be scattered light by ionized gas. The existence of many highly-ionized lines from O, Ne, Mg, Si, S, and Fe in the observed spectrum indicates that the ionized gas has a broad ionized structure, with xi=10--1000. The scattering fraction with respect to the direct light was estimated to be 0.9+/-0.2%, which indicates that the column density of the scattering region is about 3.6x10^22 cm^-2. This high-quality spectrum obtained by Suzaku can be considered a template for studies of Seyfert 2 galaxies.

We have re-analyzed FUSE data and obtained new Chandra observations of Haro 11, a local (D_L=88 Mpc) UV luminous galaxy. Haro 11 has a similar far-UV luminosity (10^10.3 L_\odot), UV surface brightness (10^9.4 L_\odot kpc^-2), SFR, and metallicity to that observed in Lyman Break Galaxies (LBGs). We show that Haro 11 has extended, soft thermal (kT~0.68 keV) X-ray emission with a luminosity and size which scales with the physical properties (e.g. SFR, stellar mass) of the host galaxy. An enhanced alpha/Fe, ratio of ~4 relative to solar abundance suggests significant supernovae enrichment. These results are consistent with the X-ray emission being produced in a shock between a supernovae driven outflow and the ambient material. The FUV spectra show strong absorption lines similar to those observed in LBG spectra. A blueshifted absorption component is identified as a wind outflowing at ~200-280 km/s. OVI\lambda\lambda1032,1038 emission, the dominant cooling mechanism for coronal gas at T~10^5.5 K is also observed. If associated with the outflow, the luminosity of the OVI emission suggests that <20% of the total mechanical energy from the supernovae and solar winds is being radiated away. This implies that radiative cooling through OVI is not significantly inhibiting the growth of the outflowing gas. In contradiction to the findings of Bergvall et al 2006, we find no convincing evidence of Lyman continuum leakage in Haro 11. We conclude that the wind has not created a `tunnel' allowing the escape of a significant fraction of Lyman continuum photons and place a limit on the escape fraction of f_{esc}<2%. Overall, both Haro 11 and a previously observed LBG analogue VV 114, provide an invaluable insight into the X-ray and FUV properties of high redshift LBGs.

A Suzaku observation of the nucleus of the radio-loud AGN Centaurus A in 2005 has yielded a broadband spectrum spanning 0.3 to 250 keV. The net exposure times after screening were: 70 ks per X-ray Imaging Spectrometer (XIS) camera, 60.8 ks for the Hard X-ray Detector (HXD) PIN, and 17.1 ks for the HXD-GSO. The hard X-rays are fit by two power-laws of the same slope, absorbed by columns of 1.5 and 7 * 10^{23} cm^{-2} respectively. The spectrum is consistent with previous suggestions that the power-law components are X-ray emission from the sub-pc VLBI jet and from Bondi accretion at the core, but it is also consistent with a partial covering interpretation. The soft band is dominated by thermal emission from the diffuse plasma and is fit well by a two-temperature VAPEC model, plus a third power-law component to account for scattered nuclear emission, jet emission, and emission from X-ray Binaries and other point sources. Narrow fluorescent emission lines from Fe, Si, S, Ar, Ca and Ni are detected. The Fe K alpha line width yields a 200 light-day lower limit on the distance from the black hole to the line-emitting gas. Fe, Ca, and S K-shell absorption edges are detected. Elemental abundances are constrained via absorption edge depths and strengths of the fluorescent and diffuse plasma emission lines. The high metallicity ([Fe/H]=+0.1) of the circumnuclear material suggests that it could not have originated in the relatively metal-poor outer halo unless enrichment by local star formation has occurred. Relative abundances are consistent with enrichment from Type II and Ia supernovae.

We present soft (0.5-2 keV) X-ray luminosity functions (XLFs) in the Great Observatories Origins Deep Survey (GOODS) fields, derived for galaxies at z~0.25 and 0.75. SED fitting was used to estimate photometric redshifts and separate galaxy types, resulting in a sample of 40 early-type galaxies and 46 late-type galaxies. We estimate k-corrections for both the X-ray/optical and X-ray/NIR flux ratios, which facilitates the separation of AGN from the normal/starburst galaxies. We fit the XLFs with a power-law model using both traditional and Markov-Chain Monte Carlo (MCMC) procedures. The XLFs differ between z<0.5 and z>0.5, at >99% significance levels for early-type, late-type and all (early and late-type) galaxies.We also fit Schechter and log-normal models to the XLFs, fitting the low and high redshift XLFs for a given sample simultaneously assuming only pure luminosity evolution. In the case of log-normal fits, the results of MCMC fitting of the local FIR luminosity function were used as priors for the faint and bright-end slopes (similar to ``fixing'' these parameters at the FIR values except here the FIR uncertainty is included). The best-fit values of the change in log L* with redshift were dlogL* = 0.23 +/- 0.16 dex (for early-type galaxies) and 0.34 +/- 0.12 dex (for late-type galaxies), corresponding to (1+z)^1.6 and (1+z)^2.3. These results were insensitive to whether the Schechter or log-normal function was adopted.

We use deep J and Ks images of the Antennae (NGC 4038/9) obtained with WIRC on the Palomar 200-inch telescope, together with the Chandra X-ray source list of Zezas et al. (2002a), to search for IR counterparts to X-ray point sources. We establish an X-ray/IR astrometric frame tie with 0.5" rms residuals over a \~4.3' field. We find 13 ``strong'' IR counterparts brighter than Ks = 17.8 mag and < 1.0" from X-ray sources, and an additional 6 ``possible'' IR counterparts between 1.0" and 1.5" from X-ray sources. The surface density of IR sources near the X-ray sources suggests only ~2 of the ``strong'' counterparts and ~3 of the ``possible'' counterparts are chance superpositions of unrelated objects. Comparing both strong and possible IR counterparts to our photometric study of ~220 Antennae, IR clusters, we find the IR counterparts to X-ray sources are \~1.2 mag more luminous in Ks than average non-X-ray clusters. We also note that the X-ray/IR matches are concentrated in the spiral arms and ``overlap'' regions of the Antennae. This implies that these X-ray sources lie in the most ``super'' of the Antennae's Super Star Clusters, and thus trace the recent massive star formation history here. Based on the N_H inferred from the X-ray sources without IR counterparts, we determine that the absence of most of the ``missing'' IR counterparts is because they are intrinsically less luminous in the IR, implying that they trace a different (possibly older) stellar population.

We report on a 100 ks Suzaku observation of the bright, nearby (z=0.008486) Seyfert 1.9 galaxy MCG -5-23-16. The broad-band (0.4-100 keV) X-ray spectrum allows us to determine the nature of the high energy emission with little ambiguity. The X-ray continuum consists of a cutoff power-law of photon index $\Gamma=1.9$, absorbed through Compton-thin matter of column density $N_{\rm H}=1.6\times10^{22}$ cm$^{-2}$. A soft excess is observed below 1 keV and is likely a combination of emission from scattered continuum photons and distant photoionized gas. The iron K line profile is complex, showing narrow neutral iron K$\alpha$ and K$\beta$ emission, as well as a broad line which can be modeled by a moderately inclined accretion disk. The line profile shows either the disk is truncated at a few tens of gravitational radii, or the disk emissivity profile is relatively flat. A strong Compton reflection component is detected above 10 keV, which is best modeled by a combination of reflection off distant matter and the accretion disk. The reflection component does not appear to vary. The overall picture is that this Seyfert 1.9 galaxy is viewed at moderate (50 degrees) inclination through Compton-thin matter at the edge of a Compton-thick torus covering $2\pi$ steradians, consistent with unified models.

We have developed a framework for the Monte-Carlo simulation of the X-Ray Telescopes (XRT) and the X-ray Imaging Spectrometers (XIS) onboard Suzaku, mainly for the scientific analysis of spatially and spectroscopically complex celestial sources. A photon-by-photon instrumental simulator is built on the ANL platform, which has been successfully used in ASCA data analysis. The simulator has a modular structure, in which the XRT simulation is based on a ray-tracing library, while the XIS simulation utilizes a spectral "Redistribution Matrix File" (RMF), generated separately by other tools. Instrumental characteristics and calibration results, e.g., XRT geometry, reflectivity, mutual alignments, thermal shield transmission, build-up of the contamination on the XIS optical blocking filters (OBF), are incorporated as completely as possible. Most of this information is available in the form of the FITS (Flexible Image Transport System) files in the standard calibration database (CALDB). This simulator can also be utilized to generate an "Ancillary Response File" (ARF), which describes the XRT response and the amount of OBF contamination. The ARF is dependent on the spatial distribution of the celestial target and the photon accumulation region on the detector, as well as observing conditions such as the observation date and satellite attitude. We describe principles of the simulator and the ARF generator, and demonstrate their performance in comparison with in-flight data.

The extended X-ray emission from 'the Cap' region located 11' (11.6 kpc) above the disk of the starburst galaxy M82 has been observed with Suzaku and XMM-Newton. Owing to the good energy resolution and the large collecting area of the XIS on Suzaku, combined with similar properties of the EPIC instrument on XMM-Newton, we have clearly detected K-shell emission lines from O VII, O VIII, Ne X, Mg XI, Mg XII and the Fe-L complex. Two optically-thin thermal plasma components are required to fit the observed X-ray spectra. We have determined the metal abundances of O, Ne, Mg, Si and Fe in this region for the first time. Their metal abundance ratios agree well with those of metal-poor stars and the model prediction of metals synthesized by type-II supernovae, but they are not consistent with the metallicities of type-Ia supernovae. This result is support for the idea that the origin of the metals in the Cap is type-II supernovae explosions occurring in the starburst regions in the M82 galaxy. We discuss the possible contribution from sputtered dust grains to the metals in the Cap. An emission line consistent with the C VI transition of n=4 to 1 at 0.459 keV is marginally detected, although it is not statistically significant at the 99% confidence level; the presence of this line would suggest charge-exchange processes in the Cap.

We present detailed time-averaged X-ray spectroscopy in the 0.5--10 keV band of the Seyfert~1.9 galaxy NGC 2992 with the Suzaku X-ray Imaging Spectrometers (XIS). We model the complex continuum in detail. There is an Fe K line emission complex that we model with broad and narrow lines and we show that the intensities of the two components are decoupled at a confidence level >3sigma. The broad Fe K line has an EW of 118 (+32,-61) eV and could originate in an accretion disk (with inclination angle greater than ~30 degrees). The narrow Fe Kalpha line has an EW of 163 (+47,-26) eV and is unresolved FWHM <4090 km/s) and likely originates in distant matter. The absolute flux in the narrow line implies that the column density out of the line-of-sight could be much higher than measured in the line-of-sight, and that the mean (historically-averaged) continuum luminosity responsible for forming the line could be a factor of several higher than that measured from the data. We also detect the narrow Fe Kbeta line with a high signal-to-noise ratio and describe a new robust method to constrain the ionization state of Fe responsible for the Fe Kalpha and Fe Kbeta lines that does not require any knowledge of possible gravitational and Doppler energy shifts affecting the line energies. For the distant line-emitting matter (e.g. the putative obscuring torus) we deduce that the predominant ionization state is lower than Fe VIII (at 99% confidence), conservatively taking into account residual calibration uncertainties in the XIS energy scale and theoretical and experimental uncertainties in the Fe K fluorescent line energies. From the limits on a possible Compton-reflection continuum it is likely that the narrow Fe Kalpha and Fe Kbeta lines originate in a Compton-thin structure.

We present a systematic analysis of archival HST WFPC2 ``Association'' data sets that correlate with the Chandra positions of a set of 44 ultra-luminous X-ray sources (ULXs) of nearby galaxies. We have improved the Chandra-HST relative astrometry whenever possible. Disparate numbers of potential ULX counterparts are found, and in some cases none are found. The lack of or low number of counterparts in some cases may be due to insufficient depth in the WFPC2 images. Particularly in late-type galaxies, the HST image in the ULX region was often complex or crowded. We therefore address various scenarios for the nature of the ULX since it is not known which, if any, of the sources found are true counterparts. The optical luminosities of the sources are typically in the range 10^4-6 L_sun. In several cases color information is available, with the colors roughly tending to be more red in early-type galaxies. This suggests that, in general, the (potential) counterparts found in early-type galaxies are likely to be older stellar populations, and are probably globular clusters. Several early-type galaxy counterparts have blue colors, which may be due to younger stellar populations in the host galaxies, however these could also be background sources. In spiral galaxies the sources may also be due to localized structure in the disks rather than bound stellar systems. Alternatively some of the counterparts in late-type galaxies may be isolated supergiant stars. The observed X-ray/optical flux ratio is diluted by the optical emission of the cluster in cases where the system is an X-ray binary in a cluster, particularly in the case of a low-mass X-ray binaries in old cluster. (abridged)

We have analyzed FUSE, XMM, and Chandra observations of VV 114, a local galaxy merger with strong similarities to typical high-redshift Lyman Break Galaxies (LBGs). Diffuse thermal X-ray emission encompassing VV114 has been observed by Chandra and XMM. This region of hot (kT~0.59 keV) gas has an enhanced alpha to iron element ratio relative to solar abundances and follows the same relation as typical starbursts between its properties (luminosity, size, and temperature) and those of the starburst galaxy (star formation rate, dust temperature, galaxy mass). These results are consistent with the X-ray gas having been produced by shocks driven by a galactic superwind. The FUSE observations of VV 114 show strong, broad interstellar absorption lines with a pronounced blueshifted component(similar to what is seen in LBGs). This implies an outflow of material moving at 300-400 km/s relative to VV 114. The properties of the strong OVI absorption line are consistent with radiative cooling at the interface between the hot outrushing gas seen in X-rays and the cooler material seen in the other outflowing ions in the FUSE data. We show that the wind in VV114 has not created a ``tunnel'' that enables more than a small fraction (< few percent) of the ionizing photons from VV114 to escape into the IGM. Taken together, these data provide a more complete physical basis for understanding the outflows that seem to be generic in LBGs. This will lead to improved insight into the role that such outflows play in the evolution of galaxies and the inter-galactic medium.

We are carrying out sensitive X-ray observations with Chandra and XMM of type II quasars selected from the Sloan Digital Sky Survey based on their optical emission line properties. We present observations of four objects at redshifts 0.4 < z < 0.8 and an analysis of the archival data for four additional objects in the same redshift range. Six of the eight were detected in X-rays; five of them have sufficient signal to derive spectral information. All of the detected sources have intrinsic luminosities L(2-10 keV) > 5 x 10^43 erg s^-1. The five with sufficient counts for spectral fitting show evidence for significant absorption (N_H >~ a few x 10^22 cm^-2). At least three of the objects likely have N_H > 10^23 cm^-2; some may be Compton-thick (N_H > 10^{24} cm^-2). In the five objects for which we could fit spectra, the slopes tend to be significantly flatter than is typically observed in AGN; it is possible that this is due either to reprocessing of the nuclear emission or to a line of sight that passes through patchy absorption.

In this paper we report that a bright, X-ray source in the Antennae galaxies (NGC 4038/9), previously identified as an ultra-luminous X-ray source, is in fact a background quasar. We identify an isolated infrared and optical counterpart within 0.3+/-0.5 arcseconds the X-ray source X-37. After acquiring an optical spectrum of its counterpart, we use the narrow [OIII] and broad H_alpha emission lines to identify X-37 as a quasar at a redshift of z=0.26. Through a U, V, and K_s photometric analysis, we demonstrate that most of the observable light along this line of sight is from the quasar. We discuss the implications of this discovery and the importance of acquiring spectra for optical and IR counterparts to ULXs.

In this paper we assess the relationship of the population of Active Galactic Nuclei (AGN) selected by hard X-rays to the traditional population of AGN with strong optical emission lines. First, we study the emission-line properties of a new hard X-ray selected sample of 47 local AGN (classified optically as both Type 1 and Type 2 AGN). We find that the hard X- ray (3-20 keV) and [OIII]$\lambda$5007 optical emission-line luminosities are well-correlated over a range of about four orders-of-magnitude in luminosity (mean luminosity ratio 2.15 dex with a standard deviation of $\sigma$ = 0.51 dex). Second, we study the hard X-ray properties of a sample of 55 local AGN selected from the literature on the basis of the flux in the [OIII] line. The correlation between the hard X-ray (2-10 keV) and [OIII] luminosity for the Type 1 AGN is consistent with what is seen in the hard X-ray selected sample. However, the Type 2 AGN have a much larger range in the luminosity ratio, and many are very weak in hard X-rays (as expected for heavily absorbed AGN). We then compare the hard X-ray (3-20 keV) and [OIII] luminosity functions of AGN in the local universe. These have similar faint-end slopes with a luminosity ratio of 1.60 dex (0.55 dex smaller than the mean value for individual hard X-ray selected AGN). We conclude that at low redshift, selection by narrow optical emission- lines will recover most AGN selected by hard X-rays (with the exception of BL Lac objects). However, selection by hard X-rays misses a significant fraction of the local AGN population with strong emission lines.

We have analyzed observations with the Chandra X-ray Observatory of the diffuse emission by hot gas in 7 dwarf starburst galaxies, 6 edge-on starburst galaxies, and 9 Ultra Luminous Infrared Galaxies. These systems cover ranges of \~ 10^4 in X-ray luminosity and several thousand in star formation rate and K-band luminosity (a proxy for stellar mass). Despite this range in fundamental parameters, we find that the properties of the diffuse X-ray emission are very similar in all three classes of starburst galaxies. The spectrum of the diffuse emission is well fit by thermal emission from gas with kT ~ 0.25 to 0.8 keV and with several-times-solar abundance ratios of alpha elements to Fe. The ratio of the thermal X-ray to far-infrared luminosity is roughly constant, as is the characteristic surface brightness of the diffuse X-ray emission. The size of the diffuse X-ray source increases systematically with both far-infrared and K-Band luminosity. All three classes show strong morphological relationships between the regions of hot gas probed by the diffuse X-ray emission and the warm gas probed by optical line emission. These findings suggest that the same physical mechanism is producing the diffuse X-ray emission in the three types of starbursts, and are consistent with that mechanism being shocks driven by a galactic "superwind'' powered by the kinetic energy collectively supplied by stellar winds and supernovae in the starburst.

We have cross-correlated X-ray catalogs derived from archival Chandra ACIS observations with a Sloan Digital Sky Survey (SDSS) Data Release 2 (DR2) galaxy catalog to form a sample of 42 serendipitously X-ray detected galaxies over the redshift interval 0.03 < z < 0.25. This pilot study will help fill in the "redshift gap" between local X-ray-studied samples of normal galaxies and those in the deepest X-ray surveys. Our chief purpose is to compare optical spectroscopic diagnostics of activity (both star-formation and accretion) with X-ray properties of galaxies. Our work supports a normalization value of the X-ray-star-formation-rate (X-ray-SFR) correlation consistent with the lower values published in the literature. The difference is in the allocation of X-ray emission to high-mass X-ray binaries relative to other components such as hot gas, low-mass X-ray binaries, and/or AGN. We are able to quantify a few pitfalls in the use of lower-resolution, lower signal-to-noise optical spectroscopy to identify X-ray sources (as has necessarily been employed for many X-ray surveys). Notably, we find a few AGN that likely would have been misidentified as non-AGN sources in higher-redshift studies. However, we do not find any X-ray hard, highly X-ray-luminous galaxies lacking optical spectroscopic diagnostics of AGN activity. Such sources are members of the "X-ray Bright, Optically Normal Galaxy" (XBONG) class of AGN.

The cosmological star formation rate in the combined Chandra Deep Fields North and South is derived from our X-Ray Luminosity Function for Galaxies in these Deep Fields. Mild evolution is seen up to redshift order unity with SFR ~ (1 + z)^{2.7}. This is the first directly observed normal star-forming galaxy X-ray luminosity function (XLF) at cosmologically interesting redshifts (z>0). This provides the most direct measure yet of the X-ray derived cosmic star-formation history of the Universe. We make use of Bayesian statistical methods to classify the galaxies and the two types of AGN, finding the most useful discriminators to be the X-ray luminosity, X-ray hardness ratio, and X-ray to optical flux ratio. There is some residual AGN contamination in the sample at the bright end of the luminosity function. Incompleteness slightly flattens the XLF at the faint end of the luminosity function. The XLF has a lognormal distribution and agrees well with the radio and infrared luminosity functions. However, the XLF does not agree with the Schechter luminosity function for the H-alpha LF indicating that additional and different physical processes may be involved in the establishment of the lognormal form of the XLF. The agreement of our star formation history points with the other star formation determinations in different wavebands (IR, Radio, H-alpha) gives an interesting constraint on the IMF, and X-rays may be measuring directly the binary star formation history of the Universe. X-ray studies will continue to be useful for probing the star formation history of the universe by avoiding problems of obscuration. Star formation may therefore be measured in more detail by deep surveys with future x-ray missions.

We present the results of a 10 ks XMM-Newton observation of NGC 3998, a ``type-I'' LINER galaxy. Our goal is to test the extent to which radiatively-inefficient accretion flow (RIAF) models and/or scaled-down AGN models are consistent with the observed properties of NGC 3998. The upper-limit for narrow Fe-K emission derived from a combined fit of the XMM-Newton and BeppoSAX spectra is 25 eV, which is one of the strictest limits to date for any AGN. This significantly rules out Fe-K emission as is expected to be observed in typical Seyfert 1 galaxies. The lack of any reflection features suggests that any optically-thick, geometrically-thin accretion disk must be truncated, probably at a radius of order 100-300 R_s. RIAF models fit the UV to X-ray spectral energy distribution of NGC 3998 reasonably well. In these models the mid-IR flux also constrains the emission from any outer thin disk component that might be present. The UV to X-ray SED is also consistent with a Comptonized thin disk with a very low accretion rate, in which case the lack of Fe-K emission may be due to an ionized accretion disk. Accretion models in general do not account for the observed radio flux of NGC 3998, and the radio flux may be due to a jet. Recent jet models may also be consistent with the nuclear fluxes of NGC 3998 in general, including the X-ray, optical/UV and mid-IR bands. We also derive nuclear fluxes using archival HST WFPC2 data to constrain the SED of NGC 3998. We discuss a possible OM U band and USNO-B detection of the NGC 3998 ULX.

We present a statistical analysis of the properties of galaxies containing ultraluminous X-ray objects (ULXs). Our primary goal is to establish the fraction of galaxies containing a ULX as a function of ULX luminosity. Our sample is based on ROSAT HRI observations of galaxies. We find that ~ 12% of galaxies contain at least one ULX with L_X > 10^39 erg/s and ~ 1% of galaxies contain at least one ULX with L_X > 10^40 erg/s. These ULX frequencies are lower limits since ROSAT HRI observation would miss absorbed ULXs (i.e., with N_H >~ 10^21 cm^-2) and those within ~ 10" of the nucleus (due to the positional error circle of the ROSAT HRI). The Hubble type distribution of galaxies with a ULX differs significantly from the distribution of types for nearby RC3 galaxies, but does not differ significantly from the galaxy type distribution of galaxies observed by the HRI in general. We find no increase in the mean FIR luminosity or FIR / K band luminosity ratio for galaxies with a ULX relative to galaxies observed by the HRI in general, however this result is also most likely biased by the soft bandpass of the HRI and the relatively low number of high SFR galaxies observed by the HRI with enough sensitivity to detect a ULX.

We analyzed 1441 Chandra X-ray point sources in 32 nearby galaxies. The total point-source X-ray luminosity L_XP is well correlated with B, K, and FIR+UV luminosities of spiral host galaxies, and with the B and K luminosities for ellipticals. This suggests an intimate connection between L_XP and both the old and young stellar populations, for which K and FIR+UV luminosities are proxies for the galaxy mass M and star-formation rate SFR. We derive proportionality constants 1.3E29 erg/s/Msol and 0.7E39 erg/s/(Msol/yr), which can be used to estimate the old and young components from M and SFR, respectively. The cumulative X-ray luminosity functions for the point sources have quite different slopes for the spirals (gamma ~= 0.5-0.8) and ellipticals (gamma ~= 1.4), implying *the most luminous point sources dominate L_XP* for the spirals. Most of the point sources have X-ray colors that are consistent with either LMXBs or Ultraluminous X-ray sources (ULXs a.k.a. IXOs) and we rule out classical HMXBs (e.g. neutron-star X-ray pulsars) as contributing much to L_XP. Thus, for spirals, the ULXs dominate L_XP. We estimate that >~20% of all ULXs found in spirals originate from the older (pop II) stellar populations, indicating that many of the ULXs that have been found in spiral galaxies are in fact pop II ULXs, like those in elliptical galaxies. The linear dependence of L_XP on the SFR argues for either a steepening in the X-ray luminosity function of the young (pop I) X-ray source population at L_X >~10^(38.5-39) erg/s, or a decreasing efficiency for producing all types of young X-ray point sources as the galaxy SFR increases.

We present initial results from a Chandra survey of a complete sample of the 8 nearest (z <= 0.04) ultraluminous IR galaxies (ULIRGs), and also include the IR-luminous galaxy NGC 6240 for comparison. In this paper we use the hard X-rays (2-8 kev) to search for the possible presence of an obscured AGN. In every case, a hard X-ray source is detected in the nuclear region. If we divide the sample according to the optical/IR spectroscopic classification (starburst vs. AGN), we find that the 5 ``starburst'' ULIRGs have hard X-ray luminosities about an order-of-magnitude smaller than the 3 ``AGN'' ULIRGs. NGC 6240 has an anomalously high hard X-ray luminosity compared to the ``starburst'' ULIRGs. The Fe Kalpha line is convincingly detected in only two ULIRGs. The weakness of the Fe-K emission in these ULIRGs generally suggests that the hard X-ray spectrum is not dominated by reflection from high N_H neutral material. The hard X-ray continuum flux ranges from a few X 10^3 to a few X 10^-5 of the far-IR flux, similar to values in pure starbursts, and several orders-of-magnitude smaller than in Compton-thin AGN. The upper limits on the ratio of the Fe Kalpha to far-IR flux are below the values measured in Compton-thick type 2 Seyfert galaxies. While very large column densities of molecular gas are observed in the nuclei of these galaxies, we find no evidence that the observed X-ray sources are obscured by Compton-thick material. Thus, our new hard X-ray data do not provide direct evidence that powerful ``buried quasars'' dominate the overall energetics of most ultraluminous infrared galaxies.

XAssist is a NASA AISR-funded project for the automation of X-ray astrophysics, with emphasis on galaxies. It is nearing completion of its initially funded effort, and is working well for Chandra and ROSAT HRI data. Initial support for XMM-Newton data is present as well. It is capable of data reprocessing, source detection, and preliminary spatial, temporal and spectral analysis for each source with sufficient counts. The bulk of the system is written in Python, which in turn drives underlying software (CIAO for Chandra data, etc.). Future work will include a GUI (mainly for beginners and status monitoring) and the exposure of at least some functionality as web services. The latter will help XAssist to eventually become part of the VO, making advanced queries possible, such as determining the X-ray fluxes of counterparts to HST or SDSS sources (including the use of unpublished X-ray data), and add the ability of ``on-the-fly'' X-ray processing. Pipelines are running on ROSAT, Chandra and now XMM-Newton observations of galaxies to demonstrate XAssist's capabilities, and the results are available online (in real time) at http://www.xassist.org . XAssist itself as well as various associated projects are available for download.

Chandra X-Ray Observatory (CXO) ACIS-S spectra of three quasars which lie behind three foreground damped Lyman alpha (DLA) absorbers are analyzed in order to attempt to determine the amount of photoelectric absorption due to metals present in their x-ray spectra. These absorbers are the three largest neutral hydrogen column density absorption-line systems known at low redshift (0.313 \le z_abs \le 0.524). They have HI column densities which lie between 3E21 and 5E21 atoms/cm^2. At these redshifts the amount of photoelectric absorption at x-ray energies is primarily an indicator of the oxygen abundance. Since the column densities of these systems are so high, one would expect accurate metallicity measurements of them to yield a robust estimate of the column-density-weighted cosmic neutral-gas-phase metallicity at z \approx 0.4. We consider cases where the DLA gas has solar element abundance ratios and ones with the alpha-group element abundance ratios enhanced. For the adopted assumptions, the column-density-weighted cosmic neutral-gas-phase metallicity of the non-enhanced elements (e.g. Zn) at z \approx 0.4 likely lies in the range \approx 0.04-0.38 Z_solar.

We present the SHEEP survey for serendipitously-detected hard X-ray sources in ASCA GIS images. In a survey area of $\sim 40$ deg$^{2}$, 69 sources were detected in the 5-10 keV band to a limiting flux of $\sim 10^{-13}$ erg cm$^{-2}$ s$^{-1}$. The number counts agree with those obtained by the similar BeppoSAX HELLAS survey, and both are in close agreement with ASCA and BeppoSAX 2-10 keV surveys. Spectral analysis of the SHEEP sample reveals that the 2-10 and 5-10 keV surveys do not sample the same populations, however, as we find considerably harder spectra, with an average $\Gamma\sim1.0$ assuming no absorption. The implication is that the agreement in the number counts is coincidental, with the 5-10 keV surveys gaining approximately as many hard sources as they lose soft ones, when compared to the 2-10 keV surveys. This is hard to reconcile with standard AGN ``population synthesis'' models for the X-ray background, which posit the existence of a large population of absorbed sources. We find no evidence of the population hardening at faint fluxes, with the exception that the few very brightest objects are anomalously soft. 53 of the SHEEP sources have been covered by ROSAT in the pointed phase. Of these 32 were detected. An additional 3 were detected in the RASS. As expected the sources detected with ROSAT are systematically softer than those detected with ASCA alone, and of the sample as a whole (truncated).

ROSAT, and now Chandra, X-ray images allow studies of extranuclear X-ray point sources in galaxies other than our own. X-ray observations of normal galaxies with ROSAT and Chandra have revealed that off-nuclear, compact, Intermediate-luminosity (Lx[2-10 keV] >= 1e39 erg/s) X-ray Objects (IXOs, a.k.a. ULXs [Ultraluminous X-ray sources]) are quite common. Here we present a catalog and finding charts for 87 IXOs in 54 galaxies, derived from all of the ROSAT HRI imaging data for galaxies with cz <= 5000 km/s from the Third Reference Catalog of Bright Galaxies (RC3). We have defined the cutoff Lx for IXOs so that it is well above the Eddington luminosity of a 1.4 Msun black hole (10^38.3 erg/s), so as not to confuse IXOs with ``normal'' black hole X-ray binaries. This catalog is intended to provide a baseline for follow-up work with Chandra and XMM, and with space- and ground-based survey work at wavelengths other than X-ray. We demonstrate that elliptical galaxies with IXOs have a larger number of IXOs per galaxy than non-elliptical galaxies with IXOs, and note that they are not likely to be merely high-mass X-ray binaries with beamed X-ray emission, as may be the case for IXOs in starburst galaxies. Approximately half of the IXOs with multiple observations show X-ray variability, and many (19) of the IXOs have faint optical counterparts in DSS optical B-band images. Follow-up observations of these objects should be helpful in identifying their nature.

This paper presents a comprehensive study of the X-ray properties of low-ionization nuclear emission-line regions (LINERs) and low-luminosity Seyfert galaxies based on observations obtained with ASCA. We analyzed data of 53 observations of 21 LINERs and 17 low-luminosity Seyferts. X-ray emission has been detected in all but one object. The X-ray luminosities in the 2-10 keV band range from 4e39 ergs/s to 5e41 ergs/s, which are 1-3 orders of magnitude smaller than in classical Seyfert galaxies. The X-ray spectra of most objects are well described by a canonical model which consists of (1) a soft component from a thermal plasma with kT<1 keV and (2) a hard component represented by a power law with a photon index~1.8 or thermal bremsstrahlung emission with kT~10 keV. Several objects do not require the soft thermal component, and their continua are well fitted by a single power-law model. Some objects show heavy absorption with column densities in excess of 1e23 cm^-2. We detect in several objects Fe K line emission with equivalent widths ranging from 50 eV to 2 keV. Variability on timescales less than a day is uncommon in our sample. By comparing multiple observations made with ASCA or with published observations from other satellites, we show that at least eight objects are variable on timescales of a week to several years. We find that the morphologies of many objects, both in the soft and hard bands, are consistent with being pointlike relative to the telescope PSF; a few are clearly extended in either or both energy bands.

We present the results of a systematic analysis of Fe K emission lines in low-luminosity AGNs (LLAGNs) observed with ASCA. We used a sample of LLAGNs with small intrinsic absorption to compare with luminous AGNs and accretion models. Some objects show an Fe line at 6.4 keV or 6.7 keV, and the rest of the objects show no significant Fe K emission line. We made a composite spectrum of sources without detected Fe emission lines to search for weak lines and obtained an upper limit on the equivalent width (EW) of a narrow line at 6.4 keV of 145 eV, which is marginally consistent with EWs seen in Seyfert 1s. Even in those objects with detected Fe lines, a skewed broad component, typically observed in Seyfert 1s, is not seen. The weakness of the Fe line is consistent with the absence of an optically thick disk in the vicinity of a central black hole. Alternatively, it might be due to the ionization of a standard disk. The lack of big blue bump favors the former.

We present a detailed analysis of the hard X-ray (>2.5 keV) EPIC spectra from the first observations of NGC 4151 made by XMM-Newton. We fit the spectra with a model consisting of a power-law continuum modified by line-of-sight absorption (arising in both partially photoionized and neutral gas) plus additional iron-K emission and absorption features. This model provides an excellent overall fit to the EPIC spectra. The iron K-alpha line is well modelled as a narrow Gaussian component. In contrast to several earlier studies based on data from ASCA, a relativistically broadened iron K-alpha emission feature is not required by the XMM-Newton data. The upper limit on the flux contained in any additional broad line is ~8% of that in the narrow line. The measured intrinsic line width (sigma =32+_7 eV) may be ascribed to (i) the doublet nature of the iron K-alpha line and (ii) emission from low ionization states of iron, ranging from neutral up to ~FeXVII. The additional iron absorption edge arises in cool material and implies factor ~2 overabundance of iron in this component.

We discuss the observable consequences for the detection of galaxies in the X-ray bandpass resulting from a peak in the cosmic star-formation rate at a redshift > 1. Following White & Ghosh, we assume a large evolution in the X-ray/B luminosity ratio at z ~ 0.5-1.5 resulting from the X-ray binaries that have evolved from stars formed at z > 1-2. Using the HDF-N redshift survey data and the locally observed X-ray/B luminosity ratio as a guide, we estimate a median X-ray flux (2-10 keV) on the order of 8 x 10^-18 ergs/s/cm^2 for galaxies in the HDF-N, which is consistent with a signal derived from a stacking analysis of the HDF-N Chandra data by Brandt et al. (2001). We also predict the number counts in deep X-ray surveys expected from normal galaxies at high redshift.

Using 100ks of XMM-Newton data on the Lockman Hole, we show how the X-ray background 'spectral paradox' is being resolved. We find the summed spectra of the Type I AGN, Type II AGN, and the unidentified objects. We conclude that the hard slope of the XRB is caused predominantly by latter sources, i.e. the Type II AGN and unidentified objects.

We present the first high resolution X-ray image of the dwarf elliptical galaxy NGC 3226. The data were obtained during an observation of the nearby Seyfert Galaxy NGC 3227 using the Chandra X-ray Observatory. We detect a point X-ray source spatially consistent with the optical nucleus of NGC 3226 and a recently-detected, compact, flat-spectrum, radio source. The X-ray spectrum can be measured up to ~10 keV and is consistent with a power law with a photon index 1.7 <~ Gamma <~ 2.2, or thermal bremmstrahlung emission with 4 <~ kT <~ 10 keV. In both cases the luminosity in the 2--10 keV band ~10^{40} h_{75}^{-1} erg/s. We find marginal evidence that the nucleus varies within the observation. These characteristics support evidence from other wavebands that NGC 3226 harbors a low-luminosity, active nucleus. We also comment on two previously-unknown, fainter X-ray sources <~ 15 arcsec from the nucleus of NGC 3226. Their proximity to the nucleus (with projected distances <~ 1.3/h_{75} kpc) suggests both are within NGC 3226, and thus have luminosities (~few x 10^{38} -- few x 10^{39} erg/s) consistent with black-hole binary systems.

We report preliminary results from an arcsecond-resolution X-ray survey of nearby galaxies using ACIS on board Chandra. The total sample consists of 41 low-luminosity AGNs, including Seyferts, LINERs, and LINER/H II transition objects. In the initial subsample of 24 objects observed thus far, we detect in 62% of the objects a compact, pointlike source astrometrically coincident with either the optical or radio position of the nucleus. The high detection rate strongly suggests that the majority of the objects do contain weakly active, AGN-like cores, presumably powered by central massive black holes. The 2-10 keV luminosities of the nuclear sources range from 10^38 to 10^41 erg/s, with a median value of 2x10^38 erg/s. Our detection limit corresponds to 8x10^37 erg/s for the typical sample distance of 12 Mpc; this limit is two orders of magnitude fainter than the weakest sources of this kind previously studied using ASCA or BeppoSAX. The new data extend toward lower luminosities the known linear correlation between hard X-ray and H-alpha luminosity for broad-line AGNs. Many narrow-line objects do contain X-ray cores, consistent with either weak AGNs or X-ray binary systems, but they have X-ray luminosities a factor of 10 below the L_X-L(Halpha) relation of the broad-line sources. Their distributions of photon energies show no indication of exceptionally high absorption. The optical line emission in these nuclei is likely powered, at least in part, by stellar processes.

We present the first results of the hard (5-10 keV) ASCA GIS survey SHEEP (Search for the High Energy Extragalactic Population). We have analysed 149 fields covering an area of 39 sq. deg detecting 69 sources. Several of these appear to be associated with QSOs and Seyfert-1 galaxies but with hard X-ray spectra, probably due to high absorption. Indeed, the hardness ratio analysis shows that the spectra of the majority of our sources can be represented with a ``scatterer'' model similar to obscured Seyfert galaxies locally. According to this model, our sources present high intrinsic absorption (logN_H~23) but also significant amounts of soft X-ray emission coming from scattered light.

We report on the first deep X-ray survey with the XMM-Newton observatory during the performance verification phase. The field of the Lockman Hole, one of the best studied sky areas over a very wide range of wavelengths, has been observed. A total of ~100 ksec good exposure time has been accumulated. Combining the images of the European Photon Imaging Camera (EPIC) detectors we reach a flux limit of 0.31, 1.4 and 2.4 X 10^{-15} erg cm^{-2} s^{-1}, respectively in the 0.5-2, 2-10, and 5-10 keV band. Within an off-axis angle of 10 arcmin we detect 148, 112 and 61 sources, respectively. The log(N)-log(S) relation in the three bands is compared with previous results. In particular in the 5-10 keV band these observations present the deepest X-ray survey ever, about a factor 20 more sensitive than the previous BeppoSAX observations. Using X-ray spectral diagnostics and the set of previously known, spectroscopically identified ROSAT sources in the field, the new sources can be classified. XMM-Newton detects a significant number (~40%) of X-ray sources with hard, probably intrinsically absorbed X-ray spectra, confirming a prediction of the population synthesis models for the X-ray background.

The high XMM-Newton throughput allows a detailed investigation of the spatial, spectral and variability properties of the extended and point source emission of the starburst galaxy NGC 253 simultaneously. We characterize the brightest sources by their hardness ratios, detect a bright X-ray transient, and show the spectrum and light curve of the brightest point source, most likely a black-hole X-ray binary. The unresolved emission of two disk regions can be modeled by two thin thermal plasma components of 0.13 and 0.4 keV plus residual harder emission, with the lower temperature component originating from above the disk, the nuclear spectrum by a three temperature plasma (~0.6, 0.9, and 6 keV) with the higher temperatures increasingly absorbed. The high temperature component most likely originates from the starburst nucleus. No non-thermal component, that would point at a significant contribution from an active nucleus (AGN), is needed. Assuming that type IIa supernova remnants (SNRs) are mostly responsible for the E>4 keV emission, the detection with EPIC of the 6.7 keV line allows us to estimate a supernova rate within the nuclear starburst of 0.2 yr^-1. RGS spectra and EPIC images reveal that the limb-brightening of the plume is mostly seen in higher ionization emission lines, while in the lower ionization lines, and below 0.5 keV, the plume is more homogeneously structured. (abridged)

Low-luminosity AGN (with X-ray luminosity < 1e42 ergs/s) far outnumber ordinary AGN, and are therefore perhaps more relevant to our understanding of AGN phenomena and the relationship between AGN and host galaxies. Many normal galaxies harbor LINER and starburst nuclei, which, together with LLAGN, are a class of ``low-activity'' galaxies that have a number of surprisingly similar X-ray characteristics, despite their heterogenous optical classification. This strongly supports the hypothesis of an AGN-starburst connection. Further, X-ray observations of normal galaxies without starburst or AGN-like activity in their nuclei offer opportunities to study populations of X-ray binaries, HII regions, and warm or hot ISM under different conditions than is often the case in the Milky Way. The results of recent X-ray observations of these types of galaxies are reviewed, and what we hope to learn about both nearby and high redshift galaxies of each type from observations with forthcoming and planned satellites is discussed.

We report X-ray fluxes in the 2--10 keV band from LINERs (low-ionization nuclear emission-line regions) and low-luminosity Seyfert galaxies obtained with the ASCA satellite. Observed X-ray luminosities are in the range between 4e39 and 5e41 ergs/s, which are significantly smaller than that of the ``classical'' low-luminosity Seyfert 1 galaxy NGC 4051. We found that X-ray luminosities in 2--10 keV of LINERs with broad Halpha emission in their optical spectra (LINER 1s) are proportional to their Halpha luminosities. This correlation strongly supports the hypothesis that the dominant ionizing source in LINER 1s is photoionization by hard photons from low-luminosity AGNs. On the other hand, the X-ray luminosities of most LINERs without broad Halpha emission (LINER 2s) in our sample are lower than LINER 1s at a given Halpha luminosity. The observed X-ray luminosities in these objects are insufficient to power their Halpha luminosities, suggesting that their primary ionizing source is other than an AGN, or that an AGN, if present, is obscured even at energies above 2 keV.

We present results of new ASCA observations of the low-luminosity AGN (LLAGN) NGC 4579 obtained on 1998 December 18 and 28, and we report on detection of variability of an iron K emission line. The X-ray luminosities in the 2--10 keV band for the two observations are nearly identical (LX $\approx$ 2$\times10^{41}$ ergs/s), but they are $\sim$35% larger than that measured in 1995 July by Terashima et al. An Fe K emission line is detected at $6.39\pm0.09$ keV (source rest frame) which is lower than the line energy $6.73^{+0.13}_{-0.12}$ keV in the 1995 observation. If we fit the Fe lines with a blend of two Gaussians centered at 6.39 keV and 6.73 keV, the intensity of the 6.7 keV line decreases, while the intensity of the 6.4 keV line increases, within an interval of 3.5 yr. This variability rules out thermal plasmas in the host galaxy as the origin of the ionized Fe line in this LLAGN. The detection and variability of the 6.4 keV line indicates that cold matter subtends a large solid angle viewed from the nucleus and that it is located within $\sim1$ pc from the nucleus. It could be identified with an optically thick standard accretion disk. If this is the case, a standard accretion disk is present at the Eddington ratio of $L_{\rm Bol}/L_{\rm Eddington} \sim 2\times10^{-3}$. A broad disk-line profile is not clearly seen and the structure of the innermost part of accretion disk remains unclear.

We present ASCA observations of LINERs without broad H$\alpha$ emission in their optical spectra. The sample of "type 2" LINERs consists of NGC 404, 4111, 4192, 4457, and 4569. We have detected X-ray emission from all the objects except for NGC 404; among the detected objects are two so-called transition objects (NGC 4192 and NGC 4569), which have been postulated to be composite nuclei having both an HII region and a LINER component. The images of NGC 4111 and NGC 4569 in the soft (0.5-2 keV) and hard (2-7 keV) X-ray bands are extended on scales of several kpc. The X-ray spectra of NGC 4111, NGC 4457 and NGC 4569 are well fitted by a two-component model that consists of soft thermal emission with $kT\sim0.65$ keV and a hard component represented by a power law (photon index $\sim$ 2) or by thermal bremsstrahlung emission ($kT\sim$ several keV). The extended hard X-rays probably come from discrete sources, while the soft emission most likely originates from hot gas produced by active star formation in the host galaxy. We have found no clear evidence for the presence of active galactic nuclei (AGNs) in the sample. If an AGN component is the primary ionization source of the optical emission lines, then it must be heavily obscured with a column density significantly larger than $10^{23}$ cm$^{-2}$. Alternatively, the optical emission could be ionized by a population of exceptionally hot stars.

We obtained observations of NGC 4639 with ASCA in order to investigate its mildly active Seyfert 1 nucleus at hard X-ray energies. Koratkar et al. (1995) have previously shown that the nucleus is a pointlike source in the ROSAT soft X-ray band. We detected in the 2-10 keV band a compact central source with a luminosity of 8.3E+40 erg/s. Comparison of the ASCA data with archival data taken with the Einstein and ROSAT satellites shows that the nucleus varies on timescales of months to years. The variability could be intrinsic, or it could be caused by variable absorption. More rapid variability, on a timescale of \~10^4 s, may be present in the ASCA data. The spectrum from 0.5 to 10 keV is well described by a model consisting of a lightly absorbed (N_H = 7.3E+20 cm^-2) power law with a photon index of 1.68. We find no evidence for significant emission from a thermal plasma; if present, it can account for no more than 25% of the flux in the 0.5-2.0 keV band. The limited photon statistics of our data do not allow us to place significant limits on the presence of iron K emission. (abridged)

We report on the detection of hard (2-10 keV) X-ray variability in the starburst galaxy M82 over the course of 9 ASCA observations. Long-term variability occurred on a time scale of days, with a change in flux by a factor of up to 4, corresponding to a point-source luminosity of L(2-10 keV) ~ 6e40 ergs/s. Short-term variability with an amplitude of ~ 1.4 on a time scale of hours was observed during the longest observation. This demonstrates that a large fraction of the hard X-ray emission of M82 (depending on the flux state) is from a compact region and is probably due to an accreting source. The 2-10 keV luminosity of the source is a lower limit to its Eddington luminosity, implying a blackhole mass of at least ~ 460 solar masses, or a mass intermediate to that of normal AGN and stellar-mass blackhole candidates.

The results of a 0.4-10.0 keV ASCA spectral analysis of a sample of low-luminosity AGN (LLAGN; M51, NGC 3147, NGC 4258), low-ionization nuclear emission line regions (LINERs; NGC 3079, NGC 3310, NGC 3998, NGC 4579, NGC 4594) and starburst galaxies (M82, NGC 253, NGC 3628 and NGC 6946) are presented. In spite of the heterogeneous optical classifications of these galaxies, the X-ray spectra are fit well by a ``canonical'' model consisting of an optically-thin Raymond-Smith plasma ``soft'' component with T ~ 7 x 10^6 K and a ``hard'' component that can be modeled by either a power-law with a photon index ~ 1.7 or a thermal bremsstrahlung with T ~ 6 x 10^7 K. The soft-component 0.4-10 keV instrinsic luminosities tend to be on the order 10^39-40 ergs/s while the hard-component luminosities tend to be on the order of 10^40-41 ergs/s. The detection of line emission is discussed. An analysis of the short-term variability properties was given in Ptak et al. (1998) and detailed interpretation of these results will be given in Paper II. (abridged)

We report the results from an ASCA observation of the high-luminosity, radio-loud quasar PKS 0637-752 (redshift 0.654), covering the 0.8-15 keV band in the quasar-frame. We find the source to have a luminosity ~10^46 erg/s in the 2-10 keV band, a factor of ~3 lower than during a previous GINGA observation. The continuum appears to be well modeled by a simple power-law with a photon index of 1.64 +/- 0.07, with no evidence for absorption by material intrinsic to the quasar, or Fe-K emission (with an equivalent width < 80 eV at 90% confidence). However we do find evidence for a narrow emission line at an energy 1.60 +/- 0.07 keV and equivalent width 59 [+38,-34] eV (both in the quasar frame). Line emission at these energies has not been observed in any other active galaxy or quasar to date. We reject the possibility that this line is the result of instrumental artifacts, and briefly explore possible identifications.

(Abridged) We present the results from {\it ASCA} observations of NGC 3227 along with those from a {\it ROSAT} observation. The 0.6--10 keV continuum is consistent with a powerlaw with $\Gamma \sim 1.6$, flatter than that typically observed in Seyfert galaxies. Fe $K$-shell emission is observed with an EW and profile typical of Seyfert 1 galaxies. The {\it ASCA} observations in 1993 reveal absorption by a screen $N_{H,z}^{ion} \simeq 3\times10^{21} cm^{-2}}$ of ionized material with $U_X \simeq 0.01$. Joint analysis of the {\it ASCA} and {\it ROSAT} data show an additional screen of neutral material instrinsic to NGC 3227 with $N_{H,z}^{neu} \sim $few$\times10^{20} cm^{-2}}$. NGC 3227 exhibits spectral variability both within and between the observations, due to short-term variability in the continuum emission and longer-term variability in the absorption. Time-resolved spectroscopy show the continuum steepened by $\Delta \Gamma \simeq 0.1$ during a flare of duration $\sim 10^4$s. The column increased by a factor of $\sim10$ by the 1995 epoch, while the continuum is consistent with that observed in 1993. The 1995 data also show evidence that the warm absorber allows $\sim$10% of the nuclear emission to escape without attenuation.

As a class, LINERs and Low-Luminosity AGN tend to show little or no significant short-term variability (i.e., with time-scales less than a day). This is a marked break for the trend of increased variability in Seyfert 1 galaxies with decreased luminosity. We propose that this difference is due to the lower accretion rate in LINERs and LLAGN which is probably causing the accretion flow to be advection-dominated. This results in a larger characteristic size for the X-ray producing region than is the case in ``normal'' AGN. Short-term variability may be caused by a localized instability or occultation events, but we note that such events would likely be accompanied by broad-band spectral changes. Since the ADAF is more compact in a Kerr metric, it is possible that the X-ray emission from ADAFs around rotating blackholes would be more variable than X-ray emission from ADAFs in a Schwarzchild metric. Similar variability arguments also apply to other wavelengths, and accordingly multiwavelength monitoring of LLAGN could serve to ``map'' the ADAF regions.

We present the results of an ASCA observation of the LINER NGC 4579. A point-like X-ray source is detected at the nucleus with a 2-10 keV luminosity of 1.5x10^41 ergs/s assuming a distance of 16.8 Mpc. The X-ray spectrum is represented by a combination of a power-law with a photon index of ~1.7 and soft thermal component with kT~0.9 keV. An iron K emission line is detected at 6.73+/-0.13 keV (rest frame) with an equivalent width of 490 +180/-190 eV and is statistically significant at more than 99.9 % confidence. The line center energy is consistent with Helium-like iron and is significantly higher than 6.4 keV which is expected from fluorescence by "cold" (or a lower ionization state of) iron. The iron line profile shows no significant red tail in contrast to Seyfert 1 galaxies although the statistics are limited. The line center energy, equivalent width, and profile are consistent with an origin in an ionized accretion disk. However the large mass accretion rate necessary to ionize the accretion disk is not consistent with the observed luminosity and normal accretion models.

The ASCA results of the starburst galaxy M82 are presented. The X-rays in the 0.5--10 keV band exhibit a thin thermal spectrum with emission lines from highly ionized magnesium, silicon, and sulfur, as well as a hard tail extending to higher than 10keV energy. The soft X-rays are spatially extended, while the hard X-rays show an unresolved point-like structure with possible a long-term flux variability. The flux ratio of the emission lines and the spatially extended structure in the low-energy band indicate that at least two-temperature thin thermal plasmas are present. The abundances of the oxygen, neon, magnesium, silicon, sulfur, and iron in the thin thermal plasmas are found to be significantly lower than the cosmic value. Neither type-Ia nor type-II supernova explosions can reproduce the observed abundance ratio. The origin of the unresolved hard X-rays is uncertain, but is probably an obscured low-luminosity AGN.

We present a new, simple method for analysing gapped time-series which is particularly suited to probe the X-ray power spectra of active galactic nuclei (AGN).The method is applied to a four-day observation of the Seyfert 1 galaxy MCG -6-30-15 with ASCA. MCG -6-30-15 is well-known for rapid, large amplitude X-ray variability. In this source, our method probes the power spectrum in the range 0.001 Hz to 0.05 Hz, thus extending the high-frequency coverage by a factor of over 50 compared to previous studies. The ASCA data out a cessation of variability up to 0.016 Hz at greater than 90% confidence. The power-law index of the power spectrum is consistent with -1.4.At 0.021 Hz we can place an upper limit of 7% on the intrinsic fluctuations of the source luminosity at greater than 99% confidence. Above this frequency Poisson noise becomes dominant but we can place an upper limit of 10% on intrinsic fluctuations up to 0.05 Hz at 98.9% confidence. The data argue against orbital modulation as the origin of the variability since it would imply upper limits on the mass of the putative black hole of less than half a million solar solar masses and therefore super-Eddington luminosities.

We present the results of the first imaging X-ray observations in the 0.4-10.0 keV bandpass of the nearby starburst galaxies M82 and NGC 253. The ASCA spectra of both M82 and NGC 253 are complex with strong line emission from O, Ne, Fe, Mg, S, and Si, allowing elemental abundances to be estimated in the X-ray band for the first time in these sources. Two components are required to fit the spectra of both galaxies, with a "soft" component well described by a thermal model with a temperature of ~ 10^6-7 K and a "hard" component well described by either a thermal model (T ~ 10^8 K) or a power-law model (energy index ~ 0.8-1.0). We find that different models (with different continua) yield absolute abundances that differ by more than an order of magnitude, while relative abundances are more robust and suggest an underabundance of Fe (inferred from the Fe-L complex) relative to alpha-burning elements, possibly as a result of dust depletion. Most of the soft flux (which originates mostly within the central kpc of M82 and NGC 253) is consistent with starburst models of supernovae-heated ISM and, to a lesser extent, starburst-driven superwind emission and the direct emission from supernova (SN). The hard component in both galaxies may have some contribution from ~ 10^8 K superwind emission or individual SN, although most of the emission probably originates in point sources (most likely blackhole candidates or mini-AGN with M_Edd > 3-20 solar masses) and, possibly, inverse-Compton scattering of IR photons. The similarity of the spectral X-ray characteristics of NGC 253 and M82 to some LINERs and low-luminosity AGN suggests a link between AGN and starbursts (e.g., both may contain an accretion-driven emission component).

We present the results of ASCA observations of a heterogenous sample of 15 spiral galaxies. 8 are LINERs or low-luminosity AGN (LLAGN), 5 are starburst galaxies and 2 are normal spiral galaxies. We find that in all cases the ASCA spectra can be described by a canonical model consisting of a power-law with a photon index ~ 1.7-2.0, plus a soft optically thin emission component with kT ~ 0.6-0.8 keV. The implied element abundances are often sub-solar. The soft component is usually extended and the nuclear, point-like emission is sometimes absorbed by column densities in the range log nH ~ 21-23 cm-2. The relative luminosities of the soft and hard components vary from galaxy to galaxy. For the LINERs, the 2-10 keV luminosity of the hard component is typically log L ~ 40-41 ergs/s whereas the 0.5-2.0 keV luminosity of the soft component is typically log L ~ 39-40 ergs/s. For starbursts, the 2-10 keV luminosity of is log L ~ 39-40 ergs/s, somewhat lower than the corresponding luminosity of most of the LINERs in our sample. The hard component is similar to the observed X-ray spectra of quasars and also to the intrinsic X-ray spectra of classical Seyfert galaxies. Most of the galaxies in our sample exhibit no significant (dI/I > 20%) short-term variability (with timescales of a day or less) whereas long-term variability is common. We present a case study of the LINER M81 in detail where there is evidence of large-amplitude (dI/I ~ 70%) variability over several weeks. There is also clear evidence for a broad, complex Fe-K emission line which is compatible with an origin in an accretion disk viewed at ~ 40 deg. These results suggest a strong connection between classical AGN, LINERs, and starburst galaxies.