We analyze two dedicated NuSTAR observations with exposure ${\sim}190$ ks located ${\sim}10^\circ$ from the Galactic plane, one above and the other below, to search for x-ray lines from the radiative decay of sterile-neutrino dark matter. These fields were chosen to minimize astrophysical x-ray backgrounds while remaining near the densest region of the dark matter halo. We find no evidence of anomalous x-ray lines in the energy range 5--20 keV, corresponding to sterile neutrino masses 10--40 keV. Interpreted in the context of sterile neutrinos produced via neutrino mixing, these observations provide the leading constraints in the mass range 10--12 keV, improving upon previous constraints in this range by a factor ${\sim}2$. We also compare our results to Monte Carlo simulations, showing that the fluctuations in our derived limit are not dominated by systematic effects. An updated model of the instrumental background, which is currently under development, will improve NuSTAR's sensitivity to anomalous x-ray lines, particularly for energies 3--5 keV.

The Ophiuchus galaxy cluster exhibits a curious concave gas density discontinuity at the edge of its cool core. It was discovered in the Chandra X-ray image by Werner and collaborators, who considered a possibility of it being a boundary of an AGN-inflated bubble located outside the core, but discounted this possibility because it required much too powerful an AGN outburst. Using low-frequency (72-240 MHz) radio data from MWA GLEAM and GMRT, we found that the X-ray structure is, in fact, a giant cavity in the X-ray gas filled with diffuse radio emission with an extraordinarily steep radio spectrum. It thus appears to be a very aged fossil of the most powerful AGN outburst seen in any galaxy cluster ($pV\sim 5\times 10^{61}$ erg for this cavity). There is no apparent diametrically opposite counterpart either in X-ray or in the radio. It may have aged out of the observable radio band because of the cluster asymmetry. At present, the central AGN exhibits only a weak radio source, so it should have been much more powerful in the past to have produced such a bubble. The AGN is currently starved of accreting cool gas because the gas density peak is displaced by core sloshing. The sloshing itself could have been set off by this extraordinary explosion if it had occurred in an asymmetric gas core. This dinosaur may be an early example of a new class of sources to be uncovered by low-frequency surveys of galaxy clusters.

The thermal Sunyaev-Zeldovich (SZ) effect presents a relatively new tool for characterizing galaxy cluster merger shocks, traditionally studied through X-ray observations. Widely regarded as the "textbook example" of a cluster merger bow shock, the western shock front in the Bullet Cluster (1E0657-56) represents the ideal test case for such an SZ study. We aim to reconstruct a parametric model for the shock SZ signal by directly and jointly fitting deep, high-resolution interferometric data from the Atacama Large Millimeter/submillimeter Array (ALMA) and Atacama Compact Array (ACA) in Fourier space. The ALMA+ACA data are primarily sensitive to the electron pressure difference across the shock front. To estimate the shock Mach number $M$, this difference can be combined with the value for the upstream electron pressure derived from an independent Chandra X-ray analysis. In the case of instantaneous electron-ion temperature equilibration, we find $M=2.08^{+0.12}_{-0.12}$, in $\approx 2.4\sigma$ tension with the independent constraint from Chandra, $M_X=2.74\pm0.25$. The assumption of purely adiabatic electron temperature change across the shock leads to $M=2.53^{+0.33}_{-0.25}$, in better agreement with the X-ray estimate $M_X=2.57\pm0.23$ derived for the same heating scenario. We have demonstrated that interferometric observations of the SZ effect provide constraints on the properties of the shock in the Bullet Cluster that are highly complementary to X-ray observations. The combination of X-ray and SZ data yields a powerful probe of the shock properties, capable of measuring $M$ and addressing the question of electron-ion equilibration in cluster shocks. Our analysis is however limited by systematics related to the overall cluster geometry and the complexity of the post-shock gas distribution. To overcome these limitations, a joint analysis of SZ and X-ray data is needed.

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.

The diffuse hard X-ray emission that fills the Galactic center, bulge, and ridge is believed to arise from unresolved populations of X-ray binary systems. However, the identity of the dominant class of accreting objects in each region remains unclear. Recent studies of Fe line properties and the low-energy (<10 keV) X-ray continuum of the bulge indicate a major population fraction of non-magnetic cataclysmic variables (CVs), in particular quiescent dwarf novae. This is in contrast to previous high-energy (>10 keV) X-ray measurements of the bulge and ridge, which indicate a dominant population of magnetic CVs, in particular intermediate polars. In addition, NuSTAR broad-band measurements have uncovered a much heavier intermediate polar population in the central ~100 pc than previously assumed, raising the possibility that some fraction of this population extends further from the center. Here we use NuSTAR's large aperture for unfocused photons and its broad-band X-ray range to probe the diffuse continuum of the inner ~1-3$^\circ$ of the Galactic bulge. This allows us to constrain possible multi-temperature components of the spectrum, such as could indicate a mixture of soft and hard populations. Our emissivity is consistent with previous hard X-ray measurements in the bulge and ridge, with the diffuse X-ray luminosity tracing the stellar mass. The spectrum is well-described by a single-temperature thermal plasma with $kT \approx 8$ keV, with no significant emission above 20 keV. This supports that the bulge is dominated by quiescent dwarf novae; we find no evidence of a significant intermediate polar population in the hard X-ray band.

We present the results of a joint XMM-Newton and NuSTAR observation (200 ks) of the galaxy cluster Abell 523 at $z=0.104$. The peculiar morphology of the cluster radio halo and its outlier position in the radio power P(1.4 GHz) - X-ray luminosity plane make it an ideal candidate for the study of radio-X-ray correlations and for the search of inverse Compton (IC) emission. We constructed thermodynamic maps derived from the XMM observations to describe in detail the physical and dynamical state of the ICM. We performed a point-to-point comparison in terms of surface brightness between the X-ray and radio emissions, to quantify their morphological discrepancies. Making use of NuSTAR's hard X-ray focusing capability, we looked for IC emission both globally and locally, after modeling the purely thermal component with a multi-temperature description. The thermodynamic maps suggest the presence of a secondary merging process that could be responsible for the peculiar radio halo morphology. This hypothesis is supported by the comparison between the X-ray and radio surface brightnesses, which shows a broad intrinsic scatter and a series of outliers. The global NuSTAR spectrum can be explained by purely thermal gas emission, and there is no convincing evidence that an IC component is needed. The $3\sigma$ upper limit on the IC flux in the 20-80 keV band is in the range $\left[2.2 - 4.0\right] \times 10^{-13} \, \mathrm{erg} \, \mathrm{s}^{-1} \, \mathrm{cm}^{-2}$, implying a lower limit on the magnetic field strength in the range $B > [0.23 - 0.31] \, \mu G$. Locally, we looked for IC emission in the central region of the cluster radio halo finding a $3\sigma$ upper limit on the 20-80 keV non-thermal flux of $3.17 \times 10^{-14} \, \mathrm{erg} \, \mathrm{s}^{-1} \, \mathrm{cm}^{-2}$, corresponding to a lower limit on the magnetic field strength of $B \gtrsim 0.81 \, \mu G$.

We use a combined 1.2 Ms of $NuSTAR$ observations of M31 to search for X-ray lines from sterile neutrino dark matter decay. For the first time in a $NuSTAR$ analysis, we consistently take into account the signal contribution from both the focused and unfocused fields of view. We also reduce the modeling systematic uncertainty by performing spectral fits to each observation individually and statistically combining the results, instead of stacking the spectra. We find no evidence of unknown lines, and thus derive limits on the sterile neutrino parameters. Our results place stringent constraints for dark matter masses $\gtrsim 12$ keV, which reduces the available parameter space for sterile neutrino dark matter produced via neutrino mixing ($e.g.$, in the $\nu$MSM) by approximately one-third. Additional $NuSTAR$ observations, together with improved low-energy background modeling, could probe the remaining parameter space in the future. Lastly, we also report model-independent limits on generic dark matter decay rates and annihilation cross sections.

Cosmic-ray acceleration has been a long-standing mystery and despite more than a century of study, we still do not have a complete census of acceleration mechanisms. The collision of strong stellar winds in massive binary systems creates powerful shocks, which have been expected to produce high-energy cosmic-rays through Fermi acceleration at the shock interface. The accelerated particles should collide with stellar photons or ambient material, producing non-thermal emission observable in X-rays and gamma-rays. The supermassive binary star eta Carinae drives the strongest colliding wind shock in the solar neighborhood. Observations with non-focusing high-energy observatories indicate a high energy source near eta Carinae, but have been unable to conclusively identify eta Carinae as the source because of their relatively poor angular resolution. Here we present the first direct focussing observations of the non-thermal source in the extremely hard X-ray band, which is found to be spatially coincident with the star within several arc-seconds. These observations show that the source of non-thermal X-rays varies with the orbital phase of the binary, and that the photon index of the emission is similar to that derived through analysis of the gamma-ray spectrum. This is conclusive evidence that the high-energy emission indeed originates from non-thermal particles accelerated at colliding wind shocks.

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.

Galaxy clusters are massive dark matter-dominated systems filled with X-ray emitting, optically thin plasma. Their large size and relative simplicity (at least as astrophysical objects go) make them a unique laboratory to measure some of the interesting plasma properties that are inaccessible by other means but fundamentally important for understanding and modeling many astrophysical phenomena -- from solar flares to black hole accretion to galaxy formation and the emergence of the cosmological Large Scale Structure. While every cluster astrophysicist is eagerly anticipating the direct gas velocity measurements from the forthcoming microcalorimeters onboard XRISM, Athena and future missions such as Lynx, a number of those plasma properties can best be probed by high-resolution X-ray imaging of galaxy clusters. Chandra has obtained some trailblazing results, but only grazed the surface of such studies. In this white paper, we discuss why we need arcsecond-resolution, high collecting area, low relative background X-ray imagers (with modest spectral resolution), such as the proposed AXIS and the imaging detector of Lynx.

X-ray observations of supernova remnants (SNRs) in the last decade have shown that the presence of recombining plasmas is somewhat common in a certain type of object. The SNR W49B is the youngest, hottest, and most highly ionized among such objects and hence provides crucial information about how the recombination phase is reached during the early evolutionary phase of SNRs. In particular, spectral properties of radiative recombination continuum (RRC) from Fe are the key for constraining the detailed plasma conditions. Here we present imaging and spectral studies of W49B with Nuclear Spectroscopic Telescope Array (NuSTAR), utilizing the highest-ever sensitivity to the Fe RRC at > 8.8keV. We confirm that the Fe RRC is the most prominent at the western part of the SNR because of the lowest electron temperature (~ 1.2 keV) achieved there. Our spatially-resolved spectral analysis reveals a positive correlation between the electron temperature and the recombination timescale with a uniform initial temperature of ~ 4 keV, which is consistent with the rapid adiabatic cooling scenario as an origin of the overionization. This work demonstrates NuSTAR's suitability for studies of thermal emission, in addition to hard nonthermal X-rays, from young and middle-aged SNRs.

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 report on NuSTAR observations of the mixed morphology supernova remnant (SNR) W49B, focusing on its nonthermal emission. Whereas radio observations as well as recent gamma-ray observations evidenced particle acceleration in this SNR, nonthermal X-ray emission has not been reported so far. With the unprecedented sensitivity of NuSTAR in the hard X-ray band, we detect a significant power-law-like component extending up to $\sim 20~{\rm keV}$, most probably of nonthermal origin. The newly discovered component has a photon index of $\Gamma =1.4^{+1.0}_{-1.1}$ with an energy flux between 10 and 20 keV of $(3.3 \pm 0.7) \times 10^{-13}~{\rm erg}~{\rm cm}^{-2}~{\rm s}^{-1}$. The emission mechanism is discussed based on the NuSTAR data combined with those in other wavelengths in the literature. The NuSTAR data, in terms both of the spectral slope and of the flux, are best interpreted as nonthermal electron bremsstrahlung. If this scenario is the case, then the NuSTAR emission provides a new probe to sub-relativistic particles accelerated in the SNR.

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 report the evidence for the anti-correlation between pulsed fraction (PF) and luminosity of the X-ray pulsar SXP 1323, found for the first time in a luminosity range $10^{35}$--$10^{37}$ erg s$^{-1}$ from observations spanning 15 years. The phenomenon of a decrease in X-ray PF when the source flux increases has been observed in our pipeline analysis of other X-ray pulsars in the Small Magellanic Cloud (SMC). It is expected that the luminosity under a certain value decreases as the PF decreases due to the propeller effect. Above the propeller region, an anti-correlation between the PF and flux might occur either as a result of an increase in the un-pulsed component of the total emission or a decrease of the pulsed component. Additional modes of accretion may also be possible, such as spherical accretion and a change in emission geometry. At higher mass accretion rates, the accretion disk could also extend closer to the neutron star (NS) surface, where a reduced inner radius leads to hotter inner disk emission. These modes of plasma accretion may affect the change in the beam configuration to fan-beam dominant emission.

Radio relics at the peripheries of galaxy clusters are tracers of the elusive cluster merger shocks. We report the discovery of a single radio relic in the galaxy cluster PLCK G200.9-28.2 ($z=0.22$, $M_{500} = 2.7\pm0.2 \times 10^{14} M_{\odot}$) using the Giant Metrewave Radio Telescope at 235 and 610 MHz and the Karl G. Jansky Very Large Array at 1500 MHz. The relic has a size of $\sim 1 \times 0.28$ Mpc, an arc-like morphology and is located at 0.9 Mpc from the X-ray brightness peak in the cluster. The integrated spectral index of the relic is $1.21\pm0.15$. The spectral index map between 235 and 610 MHz shows steepening from the outer to the inner edge of the relic in line with the expectation from a cluster merger shock. Under the assumption of diffusive shock acceleration, the radio spectral index implies a Mach number of $3.3\pm1.8$ for the shock. The analysis of archival XMM Newton data shows that PLCK G200.9-28.2 consists of a northern brighter sub-cluster, and a southern sub-cluster in a state of merger. This cluster has the lowest mass among the clusters hosting single radio relics. The position of the Planck Sunyaev Ze'ldovich effect in this cluster is offset by 700 kpc from the X-ray peak in the direction of the radio relic, suggests a physical origin for the offset. Such large offsets in low mass clusters can be a useful tool to select disturbed clusters and to study the state of merger.

Pairs of radio relics are believed to form during cluster mergers, and are best observed when the merger occurs in the plane of the sky. Mergers can also produce radio halos, through complex processes likely linked to turbulent re-acceleration of cosmic-ray electrons. However, only some clusters with double relics also show a radio halo. Here, we present a novel method to derive upper limits on the radio halo emission, and analyse archival X-ray Chandra data, as well as galaxy velocity dispersions and lensing data, in order to understand the key parameter that switches on radio halo emission. We place upper limits on the halo power below the $P_{\rm 1.4 \, GHz}\, M_{500}$ correlation for some clusters, confirming that clusters with double relics have different radio properties. Computing X-ray morphological indicators, we find that clusters with double relics are associated with the most disturbed clusters. We also investigate the role of different mass-ratios and time-since-merger. Data do not indicate that the merger mass ratio has an impact on the presence or absence of radio halos (the null hypothesis that the clusters belong to the same group cannot be rejected). However, the data suggests that the absence of radio halos could be associated with early and late mergers, but the sample is too small to perform a statistical test. Our study is limited by the small number of clusters with double relics. Future surveys with LOFAR, ASKAP, MeerKat and SKA will provide larger samples to better address this issue.

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.

High-resolution X-ray spectroscopy with Hitomi was expected to resolve the origin of the faint unidentified E=3.5 keV emission line reported in several low-resolution studies of various massive systems, such as galaxies and clusters, including the Perseus cluster. We have analyzed the Hitomi first-light observation of the Perseus cluster. The emission line expected for Perseus based on the XMM-Newton signal from the large cluster sample under the dark matter decay scenario is too faint to be detectable in the Hitomi data. However, the previously reported 3.5 keV flux from Perseus was anomalously high compared to the sample-based prediction. We find no unidentified line at the reported high flux level. Taking into account the XMM measurement uncertainties for this region, the inconsistency with Hitomi is at a 99% significance for a broad dark-matter line and at 99.7% for a narrow line from the gas. We do not find anomalously high fluxes of the nearby faint K line or the Ar satellite line that were proposed as explanations for the earlier 3.5 keV detections. We do find a hint of a broad excess near the energies of high-n transitions of Sxvi (E=3.44 keV rest-frame) -- a possible signature of charge exchange in the molecular nebula and another proposed explanation for the unidentified line. While its energy is consistent with XMM pn detections, it is unlikely to explain the MOS signal. A confirmation of this interesting feature has to wait for a more sensitive observation with a future calorimeter experiment.

The distribution of elements produced in the inner-most layers of a supernova explosion is a key diagnostic for studying the collapse of massive stars. Here we present the results of a 2.4 Ms \textit{NuSTAR} observing campaign aimed at studying the supernova remnant Cassiopeia A (Cas A). We perform spatially-resolved spectroscopic analyses of the $^{44}$Ti ejecta which we use to determine the Doppler shift and thus the three-dimensional (3D) velocities of the $^{44}$Ti ejecta. We find an initial $^{44}$Ti mass of 1.54 $\pm$ 0.21 $\times 10^{-4}$ M$_{\odot}$ which has a present day average momentum direction of 340$^{\circ}$ $\pm$ 15$^{\circ}$ projected on to the plane of the sky (measured clockwise from Celestial North) and tilted by 58$^{\circ}$ $\pm$ 20$^{\circ}$ into the plane of the sky away from the observer, roughly opposite to the inferred direction of motion of the central compact object. We find some $^{44}$Ti ejecta that are clearly interior to the reverse shock and some that are clearly exterior to the reverse shock. Where we observe $^{44}$Ti ejecta exterior to the reverse shock we also see shock-heated iron; however, there are regions where we see iron but do not observe $^{44}$Ti. This suggests that the local conditions of the supernova shock during explosive nucleosynthesis varied enough to suppress the production of $^{44}$Ti in some regions by at least a factor of two, even in regions that are assumed to be the result of processes like $\alpha$-rich freezeout that should produce both iron and titanium.

We present the results of seven Suzaku mosaic observations (>200 ks in total) of the nearest non-cool core cluster, the Antlia Cluster (or Group), beyond its degree-scale virial radius in its eastern direction. The temperature is consistent with the scaled profiles of many other clusters. Its pressure follows the universal profile. The density slope in its outskirts is significantly steeper than that of the nearest cool core cluster (Virgo) with a similar temperature as Antlia, but shallower than those of the massive clusters. The entropy increases all the way out to R_200, which is consistent in value with the baseline model predicted by a gravity heating-only mechanism in the outskirts. Antlia is quite relaxed in this direction. However, the entropy inside ~R_500 is significantly higher than the baseline model, which is similar to many other nearby low mass clusters or groups. The enclosed gas-mass fraction does not exceed the cosmic value out to 1.3R_200. Thus, there is no evidence of significant gas clumping, electron-ion non-equipartition, or departure from the hydrostatic equilibrium approximation that are suggested to explain the entropy and gas fraction anomalies found in the outskirts of some massive clusters. We also present scaling relations for the gas fraction (f_gas,200), entropy (K_200), and temperature (T_500) using 22 groups and clusters with published data in the literature. The enclosed baryon fraction at R_200 is broadly consistent with the cosmic value. The power law slope of the K_200-T_500 relation is 0.638+/-0.205. The entropy deficit at R_200 cannot be fully accounted for by the bias or deviation in the gas fraction.

We present the 3-8 keV and 8-24 keV number counts of active galactic nuclei (AGN) identified in the NuSTAR extragalactic surveys. NuSTAR has now resolved 33-39% of the X-ray background in the 8-24 keV band, directly identifying AGN with obscuring columns up to approximately 1e25 /cm2. In the softer 3-8 keV band the number counts are in general agreement with those measured by XMM-Newton and Chandra over the flux range 5e-15 < S(3 - 8 keV)/(erg/cm2/s) < 1e-12 probed by NuSTAR. In the hard 8-24 keV band NuSTAR probes fluxes over the range 2e-14 < S(8-24 keV)/(erg/cm2/s) < 1e-12, a factor of approximately 100 fainter than previous measurements. The 8-24 keV number counts match predictions from AGN population synthesis models, directly confirming the existence of a population of obscured and/or hard X-ray sources inferred from the shape of the integrated cosmic X-ray background. The measured NuSTAR counts lie significantly above simple extrapolation with a Euclidian slope to low flux of the Swift/BAT 15-55 keV number counts measured at higher fluxes S(15-55 keV) > 1e-11 erg/cm2/s, reflecting the evolution of the AGN population between the Swift/BAT local (z<0.1) sample and NuSTAR's z~1 sample. CXB synthesis models, which account for AGN evolution, lie above the Swift/BAT measurements, suggesting that they do not fully capture the evolution of obscured AGN at low redshifts.

Clusters of galaxies are the most massive gravitationally-bound objects in the Universe and are still forming. They are thus important probes of cosmological parameters and a host of astrophysical processes. Knowledge of the dynamics of the pervasive hot gas, which dominates in mass over stars in a cluster, is a crucial missing ingredient. It can enable new insights into mechanical energy injection by the central supermassive black hole and the use of hydrostatic equilibrium for the determination of cluster masses. X-rays from the core of the Perseus cluster are emitted by the 50 million K diffuse hot plasma filling its gravitational potential well. The Active Galactic Nucleus of the central galaxy NGC1275 is pumping jetted energy into the surrounding intracluster medium, creating buoyant bubbles filled with relativistic plasma. These likely induce motions in the intracluster medium and heat the inner gas preventing runaway radiative cooling; a process known as Active Galactic Nucleus Feedback. Here we report on Hitomi X-ray observations of the Perseus cluster core, which reveal a remarkably quiescent atmosphere where the gas has a line-of-sight velocity dispersion of 164+/-10 km/s in a region 30-60 kpc from the central nucleus. A gradient in the line-of-sight velocity of 150+/-70 km/s is found across the 60 kpc image of the cluster core. Turbulent pressure support in the gas is 4% or less of the thermodynamic pressure, with large scale shear at most doubling that estimate. We infer that total cluster masses determined from hydrostatic equilibrium in the central regions need little correction for turbulent pressure.

The results of long XMM-Newton X-ray observations of the NW radio relic of Abell 3667 are presented. A shock is detected at the sharp outer edge of the radio relic, both in the X-ray surface brightness and the temperature profiles. The Mach number is M = 2.54^+0.80_-0.43. The temperature jump at the shock is larger than expected from the density jump, which may indicate that a dynamically important magnetic field aligned primarily parallel to the shock front is present. The gas temperature rises gradually over several arc minutes within the shock region. This could indicate that the shock energy is initially dissipated into some mix of thermal and nonthermal (e.g., turbulence) components, and that the nonthermal energy decays into heat in the post-shock region. The observed radio relic can be powered if ~0.2% of the energy dissipated in the shock goes into the (re)acceleration of relativistic electrons. We show that the observed steepening of the radio spectrum with distance behind the shock is consistent with radiative losses by the radio-emitting electrons. However, the radio spectrum immediately behind the shock is flatter than expected for linear diffusive shock acceleration of thermal electrons. This suggests that the shock re-accelerates a pre-existing population of relativistic electrons. We also detect a bright, cool region (the "Mushroom") to the south of the radio relic, which we propose is the remnant cool core of a merging subcluster, and that this subcluster was the driver for the observed NW shock. In this model, the properties of Abell 3667 are mainly the result of an offset binary merger, and the cluster is being observed about 1 Gyr after first core passage. We predict that deeper X-ray or SZ observations of the SE radio relic will reveal a second merger shock at the outer edge.

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.

Deep (103 ks) \chandra\ observations of Abell 665 have revealed rich structures in this merging galaxy cluster, including a strong shock and two cold fronts. The newly discovered shock has a Mach number of $M$ = 3.0 $\pm$ 0.6, propagating in front of a cold disrupted cloud. This makes Abell~665 the second cluster where a strong merger shock of $M \approx$ 3 has been detected, after the Bullet cluster. The shock velocity from jump conditions is consistent with (2.7 $\pm$ 0.7) $\times$ 10$^3$ km sec$^{-1}$. The new data also reveal a prominent southern cold front, with potentially heated gas ahead of it. Abell 665 also hosts a giant radio halo. There is a hint of diffuse radio emission extending to the shock at the north, which needs to be examined with better radio data. This new strong shock provides a great opportunity to study the re-acceleration model with the X-ray and radio data combined.

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 present initial results and the source catalog from the NuSTAR survey of the Extended Chandra Deep Field South (hereafter, ECDFS) - currently the deepest contiguous component of the NuSTAR extragalactic survey program. The survey covers the full ~30 arcmin x 30 arcmin area of this field to a maximum depth of ~360 ks (~220 ks when corrected for vignetting at 3-24 keV), reaching sensitivity limits of ~1.3 x 10^-14 erg/cm2/s (3-8 keV), ~3.4 x 10^-14 erg/cm2/s (8-24 keV) and ~3.0 x 10^-14 erg/cm2/s (3-24 keV). Fifty four (54) sources are detected over the full field, although five of these are found to lie below our significance threshold once contaminating flux from neighboring (i.e., blended) sources is taken into account. Of the remaining 49 that are significant, 19 are detected in the 8-24 keV band. The 8-24 keV to 3-8 keV band ratios of the twelve sources that are detected in both bands span the range 0.39-1.7, corresponding to a photon index range of Gamma ~ 0.5-2.3, with a median photon index of 1.70 +/- 0.52. The redshifts of the 49 sources in our main sample span the range z = 0.21-2.7, and their rest-frame 10-40 keV luminosities (derived from the observed 8-24 keV fluxes) span the range L(10-40 keV) ~ (0.7-300) x 10^43 erg/s, sampling below the "knee" of the X-ray luminosity function out to z ~ 0.8-1. Finally, we identify one NuSTAR source that has neither a Chandra nor an XMM-Newton counterpart, but that shows evidence of nuclear activity at infrared wavelengths, and thus may represent a genuine, new X-ray source detected by NuSTAR in the ECDFS.

To provide the census of the sources contributing to the X-ray background peak above 10 keV, NuSTAR is performing extragalactic surveys using a three-tier "wedding cake" approach. We present the NuSTAR survey of the COSMOS field, the medium sensitivity and medium area tier, covering 1.7 deg2 and overlapping with both Chandra and XMM-Newton data. This survey consists of 121 observations for a total exposure of ~3 Ms. To fully exploit these data, we developed a new detection strategy, carefully tested through extensive simulations. The survey sensitivity at 20% completeness is 5.9, 2.9 and 6.4 x 10^-14 erg/cm2/s in the 3-24 keV, 3-8 keV and 8-24 keV bands, respectively. By combining detections in 3 bands, we have a sample of 91 NuSTAR sources with 10^42 -10^45.5 erg/s luminosities and redshift z=0.04-2.5. Thirty two sources are detected in the 8-24 keV band with fluxes ~100 times fainter than sources detected by Swift-BAT. Of the 91 detections, all but four are associated with a Chandra and/or XMM-Newton point-like counterpart. One source is associated with an extended lower energy X-ray source. We present the X-ray (hardness ratio and luminosity) and optical-to-X-ray properties. The observed fraction of candidate Compton-thick AGN measured from the hardness ratio is between 13-20%. We discuss the spectral properties of NuSTAR J100259+0220.6 (ID 330) at z=0.044, with the highest hardness ratio in the entire sample. The measured column density exceeds 10^24 /cm2, implying the source is Compton-thick. This source was not previously recognized as such without the >10 keV data.

We present the first direct measurements of the rest-frame 10-40 keV X-ray luminosity function (XLF) of Active Galactic Nuclei (AGNs) based on a sample of 94 sources at 0.1 < z <3, selected at 8-24 keV energies from sources in the NuSTAR extragalactic survey program. Our results are consistent with the strong evolution of the AGN population seen in prior, lower-energy studies of the XLF. However, different models of the intrinsic distribution of absorption, which are used to correct for selection biases, give significantly different predictions for the total number of sources in our sample, leading to small, systematic differences in our binned estimates of the XLF. Adopting a model with a lower intrinsic fraction of Compton-thick sources and a larger population of sources with column densities N_H ~ 10^{23-24} /cm2 or a model with stronger Compton reflection component (with a relative normalization of R ~ 2 at all luminosities) can bring extrapolations of the XLF from 2-10 keV into agreement with our NuSTAR sample. Ultimately, X-ray spectral analysis of the NuSTAR sources is required to break this degeneracy between the distribution of absorbing column densities and the strength of the Compton reflection component and thus refine our measurements of the XLF. Furthermore, the models that successfully describe the high-redshift population seen by NuSTAR tend to over-predict previous, high-energy measurements of the local XLF, indicating that there is evolution of the AGN population that is not fully captured by the current models.

Line emission from dark matter is well motivated for some candidates e.g. sterile neutrinos. We present the first search for dark matter line emission in the 3-80keV range in a pointed observation of the Bullet Cluster with NuSTAR. We do not detect any significant line emission and instead we derive upper limits (95% CL) on the flux, and interpret these constraints in the context of sterile neutrinos and more generic dark matter candidates. NuSTAR does not have the sensitivity to constrain the recently claimed line detection at 3.5keV, but improves on the constraints for energies of 10-25keV.

We present NuSTAR observations of the powerful radio galaxy Cygnus A, focusing on the central absorbed active galactic nucleus (AGN). Cygnus A is embedded in a cool-core galaxy cluster, and hence we also examine archival XMM-Newton data to facilitate the decomposition of the spectrum into the AGN and intracluster medium (ICM) components. NuSTAR gives a source-dominated spectrum of the AGN out to >70keV. In gross terms, the NuSTAR spectrum of the AGN has the form of a power law (Gamma~1.6-1.7) absorbed by a neutral column density of N_H~1.6x10^23 cm^-2. However, we also detect curvature in the hard (>10keV) spectrum resulting from reflection by Compton-thick matter out of our line-of-sight to the X-ray source. Compton reflection, possibly from the outer accretion disk or obscuring torus, is required even permitting a high-energy cutoff in the continuum source; the limit on the cutoff energy is E_cut>111keV (90% confidence). Interestingly, the absorbed power-law plus reflection model leaves residuals suggesting the absorption/emission from a fast (15,000-26,000km/s), high column-density (N_W>3x10^23 cm^-2), highly ionized (xi~2,500 erg cm/s) wind. A second, even faster ionized wind component is also suggested by these data. We show that the ionized wind likely carries a significant mass and momentum flux, and may carry sufficient kinetic energy to exercise feedback on the host galaxy. If confirmed, the simultaneous presence of a strong wind and powerful jets in Cygnus A demonstrates that feedback from radio-jets and sub-relativistic winds are not mutually exclusive phases of AGN activity but can occur simultaneously.

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 present deep ($>$2.4 Ms) observations of the Cassiopeia A supernova remnant with {\it NuSTAR}, which operates in the 3--79 keV bandpass and is the first instrument capable of spatially resolving the remnant above 15 keV. We find that the emission is not entirely dominated by the forward shock nor by a smooth "bright ring" at the reverse shock. Instead we find that the $>$15 keV emission is dominated by knots near the center of the remnant and dimmer filaments near the remnant's outer rim. These regions are fit with unbroken power-laws in the 15--50 keV bandpass, though the central knots have a steeper ($\Gamma \sim -3.35$) spectrum than the outer filaments ($\Gamma \sim -3.06$). We argue this difference implies that the central knots are located in the 3-D interior of the remnant rather than at the outer rim of the remnant and seen in the center due to projection effects. The morphology of $>$15 keV emission does not follow that of the radio emission nor that of the low energy ($<$12 keV) X-rays, leaving the origin of the $>$15 keV emission as an open mystery. Even at the forward shock front we find less steepening of the spectrum than expected from an exponentially cut off electron distribution with a single cutoff energy. Finally, we find that the GeV emission is not associated with the bright features in the {\it NuSTAR} band while the TeV emission may be, suggesting that both hadronic and leptonic emission mechanisms may be at work.

We present the results of a 55ks NuSTAR observation of the core of the Coma Cluster. The global spectrum can be explained by thermal gas emission, with a conservative 90% upper limit to non-thermal inverse Compton (IC) emission of $5.1 \times 10^{-12}$ erg cm$^{-2}$ s$^{-1}$ in a 12 arcmin $\times$ 12 arcmin field of view . The brightness of the thermal component in this central region does not allow more stringent upper limits on the IC component when compared with non-imaging instruments with much larger fields of view where claims of detections have been made. Future mosaic NuSTAR observations of Coma will further address this issue. The temperature map shows a relatively uniform temperature distribution with a gradient from the hot northwest side to the cooler southeast, in agreement with previous measurements. The temperature determination is robust given the flat effective area and low background in the 3-20 keV band, making NuSTAR an ideal instrument to measure high temperatures in the intracluster medium.

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 use a 43-ks XMM-Newton observation to investigate the nature of sources first distinguished by a follow-up Chandra observation of the field surrounding INTEGRAL source IGR J17448-3232, which includes extended emission and a bright point source previously classified as a blazar. We establish that the extended emission is a heretofore unknown massive galaxy cluster hidden behind the Galactic bulge. The emission-weighted temperature of the cluster within the field of view is 8.8 keV, with parts of the cluster reaching temperatures of up to 12 keV; no cool core is evident. At a redshift of 0.055, the cluster is somewhat under-luminous relative to the X-ray luminosity-temperature relation, which may be attributable to its dynamical state. We present a preliminary analysis of its properties in this paper. We also confirm that the bright point source is a blazar, and we propose that it is either a flat spectrum radio quasar or a low-frequency peaked BL Lac object. We find four other fainter sources in the field, which we study and tentatively identify. Only one, which we propose is a foreground Galactic X-ray binary, is hard enough to contribute to IGR J17448-3232, but it is too faint to be significant. We thus determine that IGR J17448-3232 is in fact the galaxy cluster up to $\approx$45 keV and the blazar beyond.

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

The low-mass X-ray binary Cen X-4 is the brightest and closest (<1.2 kpc) quiescent neutron star transient. Previous 0.5-10 keV X-ray observations of Cen X-4 in quiescence identified two spectral components: soft thermal emission from the neutron star atmosphere and a hard power-law tail of unknown origin. We report here on a simultaneous observation of Cen X-4 with NuSTAR (3-79 keV) and XMM-Newton (0.3-10 keV) in 2013 January, providing the first sensitive hard X-ray spectrum of a quiescent neutron star transient. The 0.3-79 keV luminosity was 1.1 x 10^(33) erg/s (for D=1kpc), with around 60 percent in the thermal component. We clearly detect a cutoff of the hard spectral tail above 10 keV, the first time such a feature has been detected in this source class. We show that thermal Comptonization and synchrotron shock origins for the hard X-ray emission are ruled out on physical grounds. However, the hard X-ray spectrum is well fit by a thermal bremsstrahlung model with an 18 keV electron temperature, which can be understood as arising either in a hot layer above the neutron star atmosphere or in a radiatively-inefficient accretion flow (RIAF). The power-law cutoff energy may be set by the degree of Compton cooling of the bremsstrahlung electrons by thermal seed photons from the neutron star surface. Lower thermal luminosities should lead to higher (possibly undetectable) cutoff energies. We compare Cen~X-4's behavior with PSR J1023+0038, IGR J18245-2452, and XSS J12270-4859, which have shown transitions between LMXB and radio pulsar modes at a similar X-ray luminosity.

We present NuSTAR high energy X-ray observations of the pulsar wind nebula (PWN)/supernova remnant G21.5-0.9. We detect integrated emission from the nebula up to ~40 keV, and resolve individual spatial features over a broad X-ray band for the first time. The morphology seen by NuSTAR agrees well with that seen by XMM-Newton and Chandra below 10 keV. At high energies NuSTAR clearly detects non-thermal emission up to ~20 keV that extends along the eastern and northern rim of the supernova shell. The broadband images clearly demonstrate that X-ray emission from the North Spur and Eastern Limb results predominantly from non-thermal processes. We detect a break in the spatially integrated X-ray spectrum at ~9 keV that cannot be reproduced by current SED models, implying either a more complex electron injection spectrum or an additional process such as diffusion compared to what has been considered in previous work. We use spatially resolved maps to derive an energy-dependent cooling length scale, $L(E) \propto E^{m}$ with $m = -0.21 \pm 0.01$. We find this to be inconsistent with the model for the morphological evolution with energy described by Kennel & Coroniti (1984). This value, along with the observed steepening in power-law index between radio and X-ray, can be quantitatively explained as an energy-loss spectral break in the simple scaling model of Reynolds (2009), assuming particle advection dominates over diffusion. This interpretation requires a substantial departure from spherical magnetohydrodynamic (MHD), magnetic-flux-conserving outflow, most plausibly in the form of turbulent magnetic-field amplification.

Asymmetry is required by most numerical simulations of stellar core-collapse explosions, but the form it takes differs significantly among models. The spatial distribution of radioactive 44Ti, synthesized in an exploding star near the boundary between material falling back onto the collapsing core and that ejected into the surrounding medium, directly probes the explosion asymmetries. CassiopeiaA is a young, nearby, core-collapse remnant from which 44Ti emission has previously been detected but not imaged. Asymmetries in the explosion have been indirectly inferred from a high ratio of observed 44Ti emission to estimated 56Ni emission, from optical light echoes, and from jet-like features seen in the X-ray and optical ejecta. Here we report spatial maps and spectral properties of the 44Ti in Cassiopeia A. This may explain the unexpected lack of correlation between the 44Ti and iron X-ray emission, the latter being visible only in shock-heated material. The observed spatial distribution rules out symmetric explosions even with a high level of convective mixing, as well as highly asymmetric bipolar explosions resulting from a fast-rotating progenitor. Instead, these observations provide strong evidence for the development of low-mode convective instabilities in core-collapse supernovae.

The search for diffuse non-thermal inverse Compton (IC) emission from galaxy clusters at hard X-ray energies has been undertaken with many instruments, with most detections being either of low significance or controversial. Background and contamination uncertainties present in the data of non-focusing observatories result in lower sensitivity to IC emission and a greater chance of false detection. We present 266ks NuSTAR observations of the Bullet cluster, detected from 3-30 keV. NuSTAR's unprecedented hard X-ray focusing capability largely eliminates confusion between diffuse IC and point sources; however, at the highest energies the background still dominates and must be well understood. To this end, we have developed a complete background model constructed of physically inspired components constrained by extragalactic survey field observations, the specific parameters of which are derived locally from data in non-source regions of target observations. Applying the background model to the Bullet cluster data, we find that the spectrum is well - but not perfectly - described as an isothermal plasma with kT=14.2+/-0.2 keV. To slightly improve the fit, a second temperature component is added, which appears to account for lower temperature emission from the cool core, pushing the primary component to kT~15.3 keV. We see no convincing need to invoke an IC component to describe the spectrum of the Bullet cluster, and instead argue that it is dominated at all energies by emission from purely thermal gas. The conservatively derived 90% upper limit on the IC flux of 1.1e-12 erg/s/cm^2 (50-100 keV), implying a lower limit on B>0.2{\mu}G, is barely consistent with detected fluxes previously reported. In addition to discussing the possible origin of this discrepancy, we remark on the potential implications of this analysis for the prospects for detecting IC in galaxy clusters in the future.

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 report the discovery of a giant radio halo in a new, hot, X-ray luminous galaxy cluster recently found by Planck, PLCKG171.9-40.7. The radio halo was found using Giant Metrewave Radio Telescope observations at 235 MHz and 610 MHz, and in the 1.4 GHz data from a NRAO Very Large Array Sky Survey pointing that we have reanalyzed. The diffuse radio emission is coincident with the cluster X-ray emission, has an extent of ~1 Mpc and a radio power of ~5x 10^24 W/Hz at 1.4 GHz. Its integrated radio spectrum has a slope of alpha~1.8 between 235 MHz and 1.4 GHz, steeper than that of a typical giant halo. The analysis of the archival XMM-Newton X-ray data shows that the cluster is hot (~10 keV) and disturbed, consistent with X-ray selected clusters hosting radio halos. This is the first giant radio halo discovered in one of the new clusters found by Planck.

From optical spectroscopy of X-ray sources observed as part of ChaMP, we present redshifts and classifications for a total of 1569 Chandra sources from our targeted spectroscopic follow up using the FLWO, SAAO, WIYN, CTIO, KPNO, Magellan, MMT and Gemini telescopes, and from archival SDSS spectroscopy. We classify the optical counterparts as 50% BLAGN, 16% NELG, 14% ALG, and 20% stars. We detect QSOs out to z~5.5 and galaxies out to z~3. We have compiled extensive photometry from X-ray to radio bands. Together with our spectroscopic information, this enables us to derive detailed SEDs for our extragalactic sources. We fit a variety of templates to determine bolometric luminosities, and to constrain AGN and starburst components where both are present. While ~58% of X-ray Seyferts require a starburst event to fit observed photometry only 26% of the X-ray QSO population appear to have some kind of star formation contribution. This is significantly lower than for the Seyferts, especially if we take into account torus contamination at z>1 where the majority of our X-ray QSOs lie. In addition, we observe a rapid drop of the percentage of starburst contribution as X-ray luminosity increases. This is consistent with the quenching of star formation by powerful QSOs, as predicted by the merger model, or with a time lag between the peak of star formation and QSO activity. We have tested the hypothesis that there should be a strong connection between X-ray obscuration and star-formation but we do not find any association between X-ray column density and star formation rate both in the general population or the star-forming X-ray Seyferts. Our large compilation also allows us to report here the identification of 81 XBONG, 78 z>3 X-ray sources and 8 Type-2 QSO candidates. Also we have identified the highest redshift (z=5.4135) X-ray selected QSO with optical spectroscopy.

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.

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.

We present new XMM-Newton observations of the northwest (NW) radio relic region in the cluster Abell 3667. We detect a jump in the X-ray surface brightness and X-ray temperature at the sharp outer edge of the radio relic which indicate that this is the location of a merger shock with a Mach number of about 2. Comparing the radio emission to the shock properties implies that approximately 0.2% of the dissipated shock kinetic energy goes into accelerating relativistic electrons. This is an order of magnitude smaller than the efficiency of shock acceleration in many Galactic supernova remnants, which may be due to the lower Mach numbers of cluster merger shocks. The X-ray and radio properties indicate that the magnetic field strength in the radio relic is >= 3 muG, which is a very large field at a projected distance of ~2.2 Mpc from the center of a cluster. The radio spectrum is relatively flat at the shock, and steepens dramatically with distance behind the shock. This is consistent with radiative losses by the electrons and the post-shock speed determined from the X-ray properties. The Cygnus A radio source is located in a merging cluster of galaxies. This appears to be an early-stage merger. Our recent Suzaku observation confirm the presence of a hot region between the two subclusters which agrees with the predicted shocked region. The high spectral resolution of the CCDs on Suzaku allowed us to measure the radial component of the merger velocity, Delta v_r \approx 2650 km/s.

The search for diffuse non-thermal, inverse Compton (IC) emission from galaxy clusters at hard X-ray energies has been underway for many years, with most detections being either of low significance or controversial. In this work, we investigate 14-195 keV spectra from the Swift Burst Alert Telescope (BAT) all-sky survey for evidence of non-thermal excess emission above the exponentially decreasing tail of thermal emission in the flux-limited HIFLUGCS sample. To account for the thermal contribution at BAT energies, XMM-Newton EPIC spectra are extracted from coincident spatial regions so that both thermal and non-thermal spectral components can be determined simultaneously. We find marginally significant IC components in six clusters, though after closer inspection and consideration of systematic errors we are unable to claim a clear detection in any of them. The spectra of all clusters are also summed to enhance a cumulative non-thermal signal not quite detectable in individual clusters. After constructing a model based on single-temperature fits to the XMM-Newton data alone, we see no significant excess emission above that predicted by the thermal model determined at soft energies. This result also holds for the summed spectra of various subgroups, except for the subsample of clusters with diffuse radio emission. For clusters hosting a diffuse radio halo, a relic, or a mini-halo, non-thermal emission is initially detected at the \sim5-sigma confidence level - driven by clusters with mini-halos - but modeling and systematic uncertainties ultimately degrade this significance. In individual clusters, the non-thermal pressure of relativistic electrons is limited to \sim10% of the thermal electron pressure, with stricter limits for the more massive clusters, indicating that these electrons are likely not dynamically important in the central regions of clusters.

The Coma cluster of galaxies hosts the brightest radio halo known and has therefore been the target of numerous searches for associated inverse Compton (IC) emission, particularly at hard X-ray energies where the IC signal must eventually dominate over thermal emission. The most recent search with the Suzaku Hard X-ray Detector (HXD) failed to confirm previous IC detections with RXTE and BeppoSAX, instead setting an upper limit 2.5 times below their nonthermal flux. However, this discrepancy can be resolved if the IC emission is very extended, beyond the scale of the cluster radio halo. Using reconstructed sky images from the 58-month Swift BAT all sky survey, the feasibility of such a solution is investigated. Building on Renaud et al., we test and implement a method for extracting the fluxes of extended sources, assuming specified spatial distributions. BAT spectra are jointly fit with an XMM-Newton EPIC-pn spectrum derived from mosaic observations. We find no evidence for large-scale IC emission at the level expected from the previously detected nonthermal fluxes. For all nonthermal spatial distributions considered, which span the gamut of physically reasonable IC models, we determine upper limits for which the largest (most conservative) limit is <4.2x10^{-12} erg/s/cm^2 (20-80 keV), which corresponds to a lower limit on the magnetic field B>0.2uG. A nominal flux upper limit of <2.7x10^{-12} erg/s/cm^2, with corresponding B>0.25uG, is derived for the most probable IC distribution given the size of the radio halo and likely magnetic field radial profile.

The collisionless accretion shock at the outer boundary of a galaxy cluster should primarily heat the ions instead of electrons since they carry most of the kinetic energy of the infalling gas. Near the accretion shock, the density of the intracluster medium is very low and the Coulomb collisional timescale is longer than the accretion timescale. Electrons and ions may not achieve equipartition in these regions. Numerical simulations have shown that the Sunyaev-Zel'dovich observables (e.g., the integrated Comptonization parameter Y) for relaxed clusters can be biased by a few percent. The Y-mass relation can be biased if non-equipartition effects are not properly taken into account. Using a set of hydrodynamical simulations, we have calculated three potential systematic biases in the Y-mass relations introduced by non-equipartition effects during the cross-calibration or self-calibration when using the galaxy cluster abundance technique to constraint cosmological parameters. We then use a semi-analytic technique to estimate the non-equipartition effects on the distribution functions of Y (Y functions) determined from the extended Press-Schechter theory. Depending on the calibration method, we find that non-equipartition effects can induce systematic biases on the Y functions, and the values of the cosmological parameters Omega_8, sigma_8, and the dark energy equation of state parameter w can be biased by a few percent. In particular, non-equipartition effects can introduce an apparent evolution in w of a few percent in all of the systematic cases we considered. Techniques are suggested to take into account the non-equipartition effect empirically when using the cluster abundance technique to study precision cosmology. We conclude that systematic uncertainties in the Y-mass relation of even a few percent can introduce a comparable level of biases in cosmological parameter measurements.

Abell 3667 is the archetype of a merging cluster with radio relics. The NW radio relic is the brightest cluster relic or halo known, and is believed to be due to a strong merger shock. We have observed the NW relic for 40 ksec of net XMM time. We observe a global decline of temperature across the relic from 6 to 1 keV, similar to the Suzaku results. Our new observations reveal a sharp change of both temperature and surface brightness near the position of the relic. The increased X-ray emission on the relic can be equivalently well described by either a thermal or nonthermal spectral model. The parameters of the thermal model are consistent with a Mach number M~2 shock and a shock speed of ~1200 km s^-1. The energy content of the relativistic particles in the radio relic can be explained if they are (re)-accelerated by the shock with an efficiency of ~0.2%. Comparing the limit on the inverse Compton X-ray emission with the measured radio synchrotron emission, we set a lower limit to the magnetic field in the relic of 3 muG. If the emission from the relic is non-thermal, this lower limit is in fact the required magnetic field.

The brightest cluster radio halo known resides in the Coma cluster of galaxies. The relativistic electrons producing this diffuse synchrotron emission should also produce inverse Compton emission that becomes competitive with thermal emission from the ICM at hard X-ray energies. Thus far, claimed detections of this emission in Coma are controversial (Fusco-Femiano et al. 2004; Rossetti & Molendi 2004). We present a Suzaku HXD-PIN observation of the Coma cluster in order to nail down its non-thermal hard X-ray content. The contribution of thermal emission to the HXD-PIN spectrum is constrained by simultaneously fitting thermal and non-thermal models to it and a spatially equivalent spectrum derived from an XMM-Newton mosaic of the Coma field (Schuecker et al. 2004). We fail to find statistically significant evidence for non-thermal emission in the spectra, which are better described by only a single or multi-temperature model for the ICM. Including systematic uncertainties, we derive a 90% upper limit on the flux of non-thermal emission of 6.0x10^-12 erg/s/cm^2 (20-80 keV, for photon index of 2.0), which implies a lower limit on the cluster-averaged magnetic field of B>0.15 microG. Our flux upper limit is 2.5x lower than the detected non-thermal flux from RXTE (Rephaeli & Gruber 2002) and BeppoSAX (Fusco-Femiano et al. 2004). However, if the non-thermal hard X-ray emission in Coma is more spatially extended than the observed radio halo, the Suzaku HXD-PIN may miss some fraction of the emission. A detailed investigation indicates that ~50-67% of the emission might go undetected, which could make our limit consistent with these detections. The thermal interpretation of the hard Coma spectrum is consistent with recent analyses of INTEGRAL (Eckert et al. 2007) and Swift (Ajello et al. 2009) data.

Wide-band Suzaku data on the merging cluster Abell 3667 were examined for hard X-ray emission in excess to the known thermal component. Suzaku detected X-ray signals in the wide energy band from 0.5 to 40 keV. The hard X-ray (> 10 keV) flux observed by the HXD around the cluster center cannot be explained by a simple extension of the thermal emission with average temperature of ~7 keV. The emission is most likely an emission from a very hot (kT > 13.2 keV) thermal component around the cluster center, produced via a strong heating process in the merger. In the north-west radio relic, no signature of non-thermal emission was observed. Using the HXD, the overall upper-limit flux within a 34'x34' field-of-view around the relic is derived to be 5.3e-12 erg s-1 cm-2 in the 10-40 keV band, after subtracting the ICM contribution estimated using the XIS or the XMM-Newton spectra. Directly on the relic region, the upper limit is further tightened by the XIS data to be less than 7.3e-13 erg s-1 cm-2, when converted into the 10--40 keV band. The latter value suggest that the average magnetic field within the relic is higher than 1.6 uG. The non-thermal pressure due to magnetic fields and relativistic electrons may be as large as ~20% of the thermal pressure in the region.

Sensitive surveys of the Cosmic Microwave Background will detect thousands of galaxy clusters via the Sunyaev-Zel'dovich (SZ) effect. Two SZ observables, the central or maximum and integrated Comptonization parameters y_max and Y, relate in a simple way to the total cluster mass, which allow the construction of mass functions (MFs) that can be used to estimate cosmological parameters such as Omega_M, sigma_8, and the dark energy parameter w. However, clusters form from the mergers of smaller structures, events that can disrupt the equilibrium of intracluster gas upon which SZ-M relations rely. From a set of N-body/hydrodynamical simulations of binary cluster mergers, we calculate the evolution of Y and y_max over the course of merger events and find that both parameters are transiently "boosted," primarily during the first core passage. We then use a semi-analytic technique developed by Randall et al. (2002) to estimate the effect of merger boosts on the distribution functions YF and yF of Y and y_max, respectively, via cluster merger histories determined from extended Press-Schechter (PS) merger trees. We find that boosts do not induce an overall systematic effect on YFs, and the values of Omega_M, sigma_8, and w were returned to within 2% of values expected from the nonboosted YFs. The boosted yFs are significantly biased, however, causing Omega_M to be underestimated by 15-45%, sigma_8 to be overestimated by 10-25%, and w to be pushed to more negative values by 25-45%. We confirm that the integrated SZ effect, Y, is far more robust to mergers than y_max, as previously reported by Motl et al. (2005) and similarly found for the X-ray equivalent Y_X, and we conclude that Y is the superior choice for constraining cosmological parameters.

We present X-ray and optical analysis of 188 AGN identified from 497 hard X-ray (f (2.0-8.0 keV) > 2.7x10^-15 erg cm^-2 s^-1) sources in 20 Chandra fields (1.5 deg^2) forming part of the Chandra Multi-wavelength Project. These medium depth X-ray observations enable us to detect a representative subset of those sources responsible for the bulk of the 2-8 keV Cosmic X-ray Background. Brighter than our optical spectroscopic limit, we achieve a reasonable degree of completeness (77% of X-ray sources with counter-parts r'< 22.5 have been classified): broad emission line AGN (62%), narrow emission line galaxies (24%), absorption line galaxies (7%), stars (5%) or clusters (2%). We find that most X-ray unabsorbed AGN (NH<10^22 cm^-2) have optical properties characterized by broad emission lines and blue colors, similiar to optically-selected quasars from the Sloan Digital Sky Survey but with a slighly broader color distribution. However, we also find a significant population of redder (g'-i'>1.0) AGN with broad optical emission lines. Most of the X-ray absorbed AGN (10^22<NH<10^24 cm^-2) are associated with narrow emission line galaxies, with red optical colors characteristically dominated by luminous, early type galaxy hosts rather than from dust reddening of an AGN. We also find a number of atypical AGN; for instance, several luminous AGN show both strong X-ray absorption (NH>10^22 cm^-2) and broad emission lines. Overall, we find that 81% of X-ray selected AGN can be easily interpreted in the context of current AGN unification models. Most of the deviations seem to be due to an optical contribution from the host galaxies of the low luminosity AGN.