A rip in the MLI at the exit aperture of OMA, the NuSTAR optic aligned with detector focal plane module FPMA, has resulted in an increased photon flux through OMA that has manifested itself as a low energy excess. Overall, the MLI coverage has decreased by 10%, but there is an additional time-varying component, which occasionally causes the opening to increase by up to 20%. We address the problem with a calibration update, and in this paper, we describe the attributes of the problem, the implications it has on data analysis, and the solution.

We present the results from a monitoring campaign made with the Neil Gehrels Swift Observatory of the M51 galaxies, which contain several variable ultraluminous X-ray sources (ULXs). The ongoing campaign started in May 2018, and we report here on $\sim1.5$ years of observations. The campaign, which consists of 105 observations, has a typical cadence of 3--6 days, and has the goal of determining the long-term X-ray variability of the ULXs. Two of the most variable sources were ULX7 and ULX8, both of which are known to be powered by neutron stars that are exceeding their isotropic Eddington luminosities by factors of up to 100. This is further evidence that neutron star powered ULXs are the most variable. Our two main results are, first, that ULX7 exhibits a periodic flux modulation with a period of 38 days varying over a magnitude and a half in flux from peak to trough. Since the orbital period of the system is known to be 2 days, the modulation is super-orbital, which is a near-ubiquitous property of ULX pulsars. Secondly we identify a new transient ULX, M51 XT-1, the onset of which occurred during our campaign, reaching a peak luminosity of $\sim10^{40}$ erg s$^{-1}$, before gradually fading over the next $\sim200$ days until it slipped below the detection limit of our observations. Combined with the high-quality Swift/XRT lightcurve of the transient, serendipitous observations made with Chandra and XMM-Newton provide insights into the onset and evolution of a likely super-Eddington event.

Two recent observations of the nearby galaxy NGC 6946 with NuSTAR, one simultaneous with an XMM-Newton observation, provide an opportunity to examine its population of bright accreting sources from a broadband perspective. We study the three known ultraluminous X-ray sources (ULXs) in the galaxy, and find that ULX-1 and ULX-2 have very steep power-law spectra with $\Gamma=3.6^{+0.4}_{-0.3}$ in both cases. Their properties are consistent with being super-Eddington accreting sources with the majority of their hard emission obscured and down-scattered. ULX-3 (NGC 6946 X-1) is significantly detected by both XMM-Newton and NuSTAR at $L_{\rm X}=(6.5\pm0.1)\times10^{39}$ erg s$^{-1}$, and has a power-law spectrum with $\Gamma=2.51\pm0.05$. We are unable to identify a high-energy break in its spectrum like that found in other ULXs, but the soft spectrum likely hinders our ability to detect one. We also characterise the new source, ULX-4, which is only detected in the joint XMM-Newton and NuSTAR observation, at $L_{\rm X}=(2.27\pm0.07)\times10^{39}$ erg s$^{-1}$, and is absent in a Chandra observation ten days later. It has a very hard cut-off power-law spectrum with $\Gamma=0.7\pm0.1$ and $E_{\rm cut}=11^{+9}_{-4}$ keV. We do not detect pulsations from ULX-4, but its transient nature can be explained either as a neutron star ULX briefly leaving the propeller regime or as a micro-tidal disruption event induced by a stellar-mass compact object.

We report on NuSTAR observations of transient pulsations in the neutron star X-ray binary SMC X-1. The transition from non-pulsing to pulsing states was not accompanied by a large change in flux. Instead, both pulsing and non-pulsing states were observed in a single observation during the low-flux super-orbital state. During the high-state, we measure a pulse period of $P = 0.70117(9)\,\mathrm{s}$ at $T_{ref} = 56145\,\mathrm{MJD}$. Spectral analysis during non-pulsing and pulsing states reveals that the observations can be consistently modeled by an absorbed power law with a phenomenological cutoff resembling a Fermi-Dirac distribution, or by a partially obscured cutoff power law. The shapes of the underlying continua show little variability between epochs, while the covering fraction and column density vary between super-orbital states. The strength of pulsations also varies, leading us to infer that the absence and reemergence of pulsations are related to changing obscuration, such as by a warped accretion disk. SMC X-1 is accreting near or above its Eddington limit, reaching an unabsorbed X-ray luminosity of $L_{\rm X}({\rm 2-10~keV}) \approx 5 \times 10^{38}\, {\rm erg}\, {\rm s}^{-1}$. This suggests that SMC X-1 may be a useful local analog to ultraluminous X-ray pulsars (ULXPs), which likewise exhibit strong variability in their pulsed fractions, as well as flux variability on similar timescales. In particular, the gradual pulse turn-on which has been observed in M82 X-2 is similar to the behavior we observe in SMC X-1. Thus we propose that pulse fraction variability of ULXPs may also be due to variable obscuration.

Two \textit{NuSTAR} observations of the luminous X-ray pulsar LMC X-4 in October and November 2015 captured several bright accretion flares from this source, which has a long history of stable pulse and superorbital behavior. We present a timing analysis of these data in which we detect a rapid pulse turn-on in association with the accretion flares, during which the source reaches super-Eddington luminosities. Pulsations, which are normally seen from this source, are found to only occur for approximately one hour before and during the bright flares. Beyond one hour before and after the flares, we find pulsations to be weak or nonexistent, with fractional RMS amplitudes of less than 0.05. At the onset of the flare, the pulse profiles exhibit a phase shift of 0.25 cycles that could be associated with a change in the emission geometry. This increase in pulse strength occurring well before the flare cannot be explained by the propeller effect, and potentially offers a connection between the magnetic properties of pulsars that accrete close to their Eddington limits and ULX pulsars.

With the first direct detection of merging black holes in 2015, the era of gravitational wave (GW) astrophysics began. A complete picture of compact object mergers, however, requires the detection of an electromagnetic (EM) counterpart. We report ultraviolet (UV) and X-ray observations by Swift and the Nuclear Spectroscopic Telescope ARray (NuSTAR) of the EM counterpart of the binary neutron star merger GW170817. The bright, rapidly fading ultraviolet emission indicates a high mass ($\approx0.03$ solar masses) wind-driven outflow with moderate electron fraction ($Y_{e}\approx0.27$). Combined with the X-ray limits, we favor an observer viewing angle of $\approx 30^{\circ}$ away from the orbital rotation axis, which avoids both obscuration from the heaviest elements in the orbital plane and a direct view of any ultra-relativistic, highly collimated ejecta (a gamma-ray burst afterglow).

We apply three data science techniques, Nonnegative Matrix Factorization (NMF), Principal Component Analysis (PCA) and Independent Component Analysis (ICA), to simulated X-ray energy spectra of a particular class of super-massive black holes. Two competing physical models, one whose variable components are additive and the other whose variable components are multiplicative, are known to successfully describe X-ray spectral variation of these super-massive black holes, within accuracy of the contemporary observation. We hope to utilize these techniques to compare the viability of the models by probing the mathematical structure of the observed spectra, while comparing advantages and disadvantages of each technique. We find that PCA is best to determine the dimensionality of a dataset, while NMF is better suited for interpreting spectral components and comparing them in terms of the physical models in question. ICA is able to reconstruct the parameters responsible for spectral variation. In addition, we find that the results of these techniques are sufficiently different that applying them to observed data may be a useful test in comparing the accuracy of the two spectral models.

We performed time lag analysis on the X-ray light curves of Centaurus A (Cen A) obtained by the Gas Slit Camera (GSC) aboard the Monitor of All-sky X-ray Image ({\it MAXI}) in three energy bands (2--4 keV, 4--10 keV, and 10--20 keV). We discovered a soft X-ray lag relative to higher energies (soft lag) on a timescale of days in a flaring episode by employing the discrete correlation function (DCF) and the z-transformed discrete correlation function (ZDCF) method. In the episode, a peak and a centroid in the DCF and the ZDCF was observed at a soft lag of $\sim 5$ days in 2--4 keV versus 4--10 keV and in 4--10 keV versus 10--20 keV, and $\sim 10$ days in 2--4 keV versus 10--20 keV. We found it difficult to explain the observed X-ray variation by a single energy injection with the one-zone synchrotron self-Compton (SSC) model, in which the soft lags in these three energy bands reflect the different cooling times of the relativistic electrons, by assuming the magnetic field and minimum Lorentz factor estimated from a broad-band SED. Alternatively, if the phenomenon is interpreted as cooling of Comptonizing electrons in a corona covering the accretion disk, the temperature of the corona producing the variable X-rays should be $\sim$10 keV for reconciliation with the soft lag in the energy range of 2--20 keV.

We present the light curve data of a remarkable blazer 3C 454.3 (z=0.859) in optical, X-ray, and gamma-ray bands. Since January 2008, we have been monitoring this object using the 50 cm MITSuME, a optical telescope, and detected several flares including extraordinary and simultaneous flares in the $\gamma$-ray and optical bands in November 2010. Additionally, the Monitor of All-sky Image (MAXI) has been observing 3C 454.3 continuously since August 2009. Using these data and gamma-ray flux observed with Fermi-LAT, we discuss features and correlations of flux variations between the energy bands.

Multi-wavelength observations of the black-widow binary system 2FGL J2339.6-0532 are reported. The Fermi gamma-ray source 2FGL J2339.6-0532 was recently categorized as a black widow in which a recycled millisecond pulsar (MSP) is evaporating up the companion star with its powerful pulsar wind. Our optical observations show clear sinusoidal light curves due to the asymmetric temperature distribution of the companion star. Assuming a simple geometry, we constrained the range of the inclination angle of the binary system to 52$^{\circ}$ < i < 59$^{\circ}$, which enables us to discuss the interaction between the pulsar wind and the companion in detail. The X-ray spectrum consists of two components: a soft, steady component that seems to originate from the surface of the MSP, and a hard variable component from the wind-termination shock near the companion star. The measured X-ray luminosity is comparable to the bolometric luminosity of the companion, meaning that the heating efficiency is less than 0.5. In the companion orbit, 10$^{11}$ cm from the pulsar, the pulsar wind is already in particle dominant-stage, with a magnetization parameter of $\sigma$ < 0.1. In addition, we precisely investigated the time variations of the X-ray periodograms and detected a weakening of orbital modulation. The observed phenomenon may be related to an unstable pulsar-wind activity or a weak mass accretion, both of which can result in the temporal extinction of radio-pulse.

We present results from a new set of 30 cosmological simulations of galaxy clusters, including the effects of radiative cooling, star formation, supernova feedback, black hole growth and AGN feedback. We first demonstrate that our AGN model is capable of reproducing the observed cluster pressure profile at redshift, z~0, once the AGN heating temperature of the targeted particles is made to scale with the final virial temperature of the halo. This allows the ejected gas to reach larger radii in higher-mass clusters than would be possible had a fixed heating temperature been used. Such a model also successfully reduces the star formation rate in brightest cluster galaxies and broadly reproduces a number of other observational properties at low redshift, including baryon, gas and star fractions; entropy profiles outside the core; and the X-ray luminosity-mass relation. Our results are consistent with the notion that the excess entropy is generated via selective removal of the densest material through radiative cooling; supernova and AGN feedback largely serve as regulation mechanisms, moving heated gas out of galaxies and away from cluster cores. However, our simulations fail to address a number of serious issues; for example, they are incapable of reproducing the shape and diversity of the observed entropy profiles within the core region. We also show that the stellar and black hole masses are sensitive to numerical resolution, particularly the gravitational softening length; a smaller value leads to more efficient black hole growth at early times and a smaller central galaxy.

We develop a statistical method to measure the interaction cross-section of Dark Matter, exploiting the continuous minor merger events in which small substructures fall into galaxy clusters. We find that by taking the ratio of the distances between the galaxies and Dark Matter, and galaxies and gas in accreting sub-halos, we form a quantity that can be statistically averaged over a large sample of systems whilst removing any inherent line-of-sight projections. In order to interpret this ratio as a cross-section of Dark Matter we derive an analytical description of sub-halo infall which encompasses; the force of the main cluster potential, the drag on a gas sub-halo, a model for Dark Matter self-interactions and the resulting sub-halo drag, the force on the gas and galaxies due to the Dark Matter sub-halo potential, and finally the buoyancy on the gas and Dark Matter. We create mock observations from cosmological simulations of structure formation and find that collisionless Dark Matter becomes physically separated from X-ray gas by up to 20h^-1 kpc. Adding realistic levels of noise, we are able to predict achievable constraints from observational data. Current archival data should be able to detect a difference in the dynamical behaviour of Dark Matter and standard model particles at 6 sigma, and measure the total interaction cross-section sigma/m with 68% confidence limits of +/- 1cm2g^-1. We note that this method is not restricted by the limited number of major merging events and is easily extended to large samples of clusters from future surveys which could potentially push statistical errors to 0.1cm^2g^-1.