The High Inclination Solar Mission (HISM) is a concept for an out-of-the-ecliptic mission for observing the Sun and the heliosphere. The mission profile is largely based on the Solar Polar Imager concept: initially spiraling in to a 0.48 AU ecliptic orbit, then increasing the orbital inclination at a rate of $\sim 10$ degrees per year, ultimately reaching a heliographic inclination of $>$75 degrees. The orbital profile is achieved using solar sails derived from the technology currently being developed for the Solar Cruiser mission, currently under development. HISM remote sensing instruments comprise an imaging spectropolarimeter (Doppler imager / magnetograph) and a visible light coronagraph. The in-situ instruments include a Faraday cup, an ion composition spectrometer, and magnetometers. Plasma wave measurements are made with electrical antennas and high speed magnetometers. The $7,000\,\mathrm{m}^2$ sail used in the mission assessment is a direct extension of the 4-quadrant $1,666\,\mathrm{m}^2$ Solar Cruiser design and employs the same type of high strength composite boom, deployment mechanism, and membrane technology. The sail system modelled is spun (~1 rpm) to assure required boom characteristics with margin. The spacecraft bus features a fine-pointing 3-axis stabilized instrument platform that allows full science observations as soon as the spacecraft reaches a solar distance of 0.48 AU.

Signatures of X-ray reprocessing (reflection) out of an accretion disk are commonly observed in the high-energy spectrum of accreting black holes, and can be used to probe the strong gravity region around these objects. In this paper, we extend previous work in the literature and we employ a full emission model for relativistic reflection in non-Kerr spacetime to demonstrate an approach that tests the Kerr black hole hypothesis. We analyze a composite spectrum obtained with the Proportional Counter Array in the Rossi X-ray Timing Explorer (RXTE), of the stellar-mass black hole GX 339-4 in its brightest hard state. With a remarkable sensitivity of ~0.1% and 40 million counts in the 3-45 keV band to capture the faint features in the reflection spectrum, we demonstrate that it is possible with existing data and an adequate model to place constraints on the black hole spin $a_*$ and the deformation parameter that quantifies the departure from the Kerr metric. Our measurement obtained with the best fit model, which should be regarded as principally a proof of concept, is $a_*=0.92^{+0.07}_{-0.12}$ and $\alpha_{13}=-0.76^{+0.78}_{-0.60}$ with a 90% confidence level and is consistent with the hypothesis that the compact object in GX 339-4 is a Kerr black hole. We also discuss how the physical model choice and the emissivity profile adopted could make an impact on the constraints of $\alpha_{13}$ and spin. To enable Kerr metric test using X-ray reflection spectroscopy, it is essential to improve our astrophysical understanding of accreting black holes, e.g., the natures of accretion flow and corona.

We explore the accretion properties of the black hole X-ray binary \j1550\ during its outbursts in 1998/99 and 2000. We model the disk, corona, and reflection components of X-ray spectra taken with the {\it Rossi X-ray Timing Explorer} (\rxte), using the {\tt relxill} suite of reflection models. The key result of our modeling is that the reflection spectrum in the very soft state is best explained by disk self-irradiation, i.e., photons from the inner disk are bent by the strong gravity of the black hole, and reflected off the disk surface. This is the first known detection of thermal disk radiation reflecting off the inner disk. There is also an apparent absorption line at $\sim6.9$ keV which may be evidence of an ionized disk wind. The coronal electron temperature ($kT_{\rm e}$) is, as expected, lower in the brighter outburst of 1998/99, explained qualitatively by more efficient coronal cooling due to irradiating disk photons. The disk inner radius is consistent with being within a few times the innermost stable circular orbit (ISCO) throughout the bright-hard-to-soft states (10s of $r_{\rm g}$ in gravitational units). The disk inclination is low during the hard state, disagreeing with the binary inclination value, and very close to $90^{\circ}$ in the soft state, recovering to a lower value when adopting a blackbody spectrum as the irradiating continuum.

X-ray reflection spectroscopy is a powerful tool for probing the strong gravity region of black holes and can be used for testing general relativity in the strong field regime. Simplifications of the available relativistic reflection models limit the capability of performing accurate measurements of the properties of black holes. In this paper, we present an extension of the model RELXILL_NK in which the accretion disk has a finite thickness rather than being infinitesimally thin. We employ the accretion disk geometry proposed by Taylor & Reynolds (2018) and we construct relativistic reflection models for different values of the mass accretion rate of the black hole. We apply the new model to high quality Suzaku data of the X-ray binary GRS 1915+105 to explore the impact of the thickness of the disk on tests of the Kerr metric.

We present results from the Nuclear Spectroscopic Telescope Array (NuSTAR) observations of the new black hole X-ray binary candidate MAXI J1631-479 at two epochs during its 2018-2019 outburst, which caught the source in a disk dominant state and a power-law dominant state. Strong relativistic disk reflection features are clearly detected, displaying significant variations in the shape and strength of the broad iron emission line between the two states. Spectral modeling of the reflection spectra reveals that the inner radius of the optically-thick accretion disk evolves from $<1.9$ $r_{\rm g}$ to $12\pm1$ $r_{\rm g}$ (statistical errors at 90% confidence level) from the disk dominant to the power-law dominant state. Assuming in the former case that the inner disk radius is consistent with being at the ISCO, we estimate a black hole spin of $a^*>0.94$. Given that the bolometric luminosity is similar in the two states, our results indicate that the disk truncation observed in MAXI J1631-479 in the power-law dominant state is unlikely to be driven by a global variation in the accretion rate. We propose that it may instead arise from local instabilities in the inner edge of the accretion disk at high accretion rates. In addition, we find an absorption feature in the spectra centered at $7.33\pm0.03$ keV during the disk dominant state, which is evidence for a rare case that an extremely fast disk wind ($v_{\rm out}=0.067^{+0.001}_{-0.004}~c$) is observed in a low-inclination black hole binary, with the viewing angle of $29\pm1^{\circ}$ as determined by the reflection modeling.

The Galactic black hole X-ray binary MAXI J1820+070 had a bright outburst in 2018 when it became the second brightest X-ray source in the Sky. It was too bright for X-ray CCD instruments such as XMM-Newton and Chandra, but was well observed by photon-counting instruments such as NICER and NuSTAR. We report here on the discovery of an excess emission component during the soft state. It is best modelled with a blackbody spectrum in addition to the regular disk emission, modelled either as diskbb or kerrbb. Its temperature varies from about 0.9 to 1.1 keV which is about 30 to 80 per cent higher than the inner disc temperature of diskbb. Its flux varies between 4 and 12 percent of the disc flux. Simulations of magnetised accretion discs have predicted the possibility of excess emission associated with a non-zero torque at the Innermost Stable Circular Orbit (ISCO) about the black hole, which from other NuSTAR studies lies at about 5 gravitational radii or about 60 km (for a black hole mass is 8 M). In this case the emitting region at the ISCO has a width varying between 1.3 and 4.6 km and would encompass the start of the plunge region where matter begins to fall freely into the black hole.

We explore a disc origin for the highly-blueshifted, variable absorption lines seen in the X-ray spectrum of the Narrow Line Seyfert 1 galaxy IRAS13224-3809. The blueshift corresponds to a velocity of about 0.25c. Such features in other Active Galactic Nuclei are often interpreted as UltraFast Outflows (UFOs). The velocity is of course present in the orbital motions of the inner disk. The absorption lines in IRAS13224-3809 are best seen when the flux is low and the reflection component of the disk is strong relative to the power-law continuum. The spectra are consistent with a model in which the reflection component passes through a thin, highly-ionized absorbing layer at the surface of the inner disc, the blue-shifted side of which dominates the flux due to relativistic aberration (the disc inclination is about 70 deg). No fast outflow need occur beyond the disc.

Aims. In the context of black-hole accretion disks, we aim to compute the plasma-environment effects on the atomic parameters used to model the decay of K-vacancy states in moderately charged iron ions, namely Fe ix - Fe xvi. Methods. We used the fully relativistic multiconfiguration Dirac-Fock (MCDF) method approximating the plasma electron-nucleus and electron-electron screenings with a time-averaged Debye-Huckel potential. Results. We report modified ionization potentials, K-threshold energies, wavelengths, radiative emission rates, and Auger widths for plasmas characterized by electron temperatures and densities in the ranges $10^5$ - $10^7$ K and $10^{18}$ - $10^{22}$ cm$^{-3}$. Conclusions. This study confirms that the high-resolution X-ray spectrometers onboard the future XRISM and ATHENA space missions will be capable of detecting the lowering of the K edges of these ions due to the extreme plasma conditions occurring in accretion disks around compact objects.

The aim of this study is to search for evidence of a common emission engine between radio giant pulses (GPs) and very-high-energy (VHE, E>100 GeV) gamma-rays from the Crab pulsar. 16 hours of simultaneous observations of the Crab pulsar at 1.4 GHz with the Effelsberg radio telescope and the Westerbork Synthesis Radio Telescope (WSRT), and at energies above 60 GeV with the Major Atmospheric Gamma-ray Imaging Cherenkov (MAGIC) telescopes were performed. We searched for a statistical correlation between the radio and VHE gamma-ray emission with search windows of different lengths and different time lags to the arrival times of a radio GP. A dedicated search for an enhancement in the number of VHE gamma-rays correlated with the occurrence of radio GPs was carried out separately for the P1 and P2 phase ranges respectively. 99444 radio GPs have been detected in the radio data sample. We find no significant correlation between the GPs and VHE photons in any of the search windows. Depending on phase cuts and the chosen search windows we find upper limits at 95% confidence level on an increase in VHE gamma-ray events correlated with radio GPs between 7% and 61% of the average Crab pulsar VHE flux for the P1 and P2 phase ranges respectively. This puts upper limits on the flux increase during a radio GP of 12% to 2900% (depending on search window duration and phase cuts) of the pulsed VHE flux. This is the most stringent upper limit on a correlation between gamma-ray emission and radio GPs reported so far.}

The Habitable Exoplanet Observatory, or HabEx, has been designed to be the Great Observatory of the 2030s. For the first time in human history, technologies have matured sufficiently to enable an affordable space-based telescope mission capable of discovering and characterizing Earthlike planets orbiting nearby bright sunlike stars in order to search for signs of habitability and biosignatures. Such a mission can also be equipped with instrumentation that will enable broad and exciting general astrophysics and planetary science not possible from current or planned facilities. HabEx is a space telescope with unique imaging and multi-object spectroscopic capabilities at wavelengths ranging from ultraviolet (UV) to near-IR. These capabilities allow for a broad suite of compelling science that cuts across the entire NASA astrophysics portfolio. HabEx has three primary science goals: (1) Seek out nearby worlds and explore their habitability; (2) Map out nearby planetary systems and understand the diversity of the worlds they contain; (3) Enable new explorations of astrophysical systems from our own solar system to external galaxies by extending our reach in the UV through near-IR. This Great Observatory science will be selected through a competed GO program, and will account for about 50% of the HabEx primary mission. The preferred HabEx architecture is a 4m, monolithic, off-axis telescope that is diffraction-limited at 0.4 microns and is in an L2 orbit. HabEx employs two starlight suppression systems: a coronagraph and a starshade, each with their own dedicated instrument.

Mrk501 is a very high-energy (VHE) gamma-ray blazar located at z=0.034. During a period of two weeks in July 2014, the highest X-ray activity of Mrk501 was observed in ~14 years of operation of the Neil Gehrels Swift Gamma-ray Burst Observatory. We characterize the broadband variability of Mrk501 from radio to VHE gamma rays, and evaluate whether it can be interpreted within theoretical scenarios widely used to explain the broadband emission from blazars. The temporal evolution of the most prominent and variable segments of the SED is described with a one-zone synchrotron self-Compton model with variations in the break energy of the electron energy distribution (EED), and with some adjustments in the magnetic field strength and spectral shape of the EED. A narrow feature at ~3 TeV was observed in the VHE spectrum measured on 2014 July 19 (MJD 56857.98), which is the day with the highest X-ray flux ($>0.3$ keV) measured during the entire Swift mission. This feature is inconsistent with the classical analytic functions to describe the measured VHE spectra (power law, log-parabola, and log-parabola with exponential cutoff) at more than 3$\sigma$. A fit with a log-parabola plus a narrow component is preferred over the fit with a single log-parabola at more than 4$\sigma$, and a dedicated Monte Carlo simulation estimated the significance of this extra component to be larger than 3$\sigma$. Under the assumption that this VHE spectral feature is real, we show that it can be reproduced with three distinct theoretical scenarios: a) a pileup in the EED due to stochastic acceleration; b) a structured jet with two-SSC emitting regions, with one region dominated by an extremely narrow EED; and c) an emission from an IC pair cascade.

Swift J1858.6-0814 was discovered by Swift-BAT on October 25, 2018. Here we report on the first follow-up NuSTAR observation of the source, which shows variability spanning two orders of magnitude in count rate on timescales of ~10-100 s. The power-spectrum of the source does not show any quasi-periodic oscillations or periodicity, but has a large fractional rms amplitude of 147%$\pm3$%, exhibiting a number of large flares throughout the observation. The hardness ratio (defined as $R_{10-79 \rm keV}/R_{3-10 \rm keV}$) of the flares tends to be soft, while the source spans a range of hardness ratios during non-flaring periods. The X-ray spectrum of the source shows strong reflection features, which become more narrow and peaked during the non-flaring intervals. We fit an absorbed relativistic reflection model to the source spectra to place physical constraints on the system. Most notably, we find that the source exhibits a large and varying intrinsic absorbing column density ($N_{\rm H}=1.4-4.2\times10^{23}$ cm$^{-2}$). This large intrinsic absorption is further supported by the energy spectra extracted from two flares observed simultaneously by NuSTAR and NICER. We find that the inner accretion disk of the source has a low inclination, $i<29^{\circ}$ ( 3$\sigma$ upper-limit), while the iron abundance in the disk is close to solar, $A_{\rm Fe}=1.0\pm0.3$. We set a 90% confidence upper limit on the inner radius of the accretion disk of $r_{\rm in}<8 \ r_{\rm ISCO}$, and, by fixing $r_{\rm in}$ to be at $r_{\rm ISCO}$, a 90% confidence lower-limit on the spin of the black hole of $a^{*}>0.0$. Lastly, we compare the properties of Swift J1858.6$-$0814 to those of V404 Cygni and V4641 Sgr, which both show rapid flaring and a strong and variable absorption.

We present SN2018kzr, the fastest declining supernova-like transient, second only to the kilonova, AT2017gfo. SN2018kzr is characterized by a peak magnitude of $M_r = -17.98$, peak bolometric luminosity of ${\sim} 1.4 \times 10^{43}$erg s$^{\mathrm{-1}}$ and a rapid decline rate of $0.48 \pm 0.03$ mag day$^{\textrm{-1}}$ in the $r$ band. The bolometric luminosity evolves too quickly to be explained by pure $^{\mathrm{56}}$Ni heating, necessitating the inclusion of an alternative powering source. Incorporating the spin-down of a magnetized neutron star adequately describes the lightcurve and we estimate a small ejecta mass of $M_\mathrm{ej} = 0.10 \pm 0.05$ $\textrm{M}_{\odot}$. Our spectral modelling suggests the ejecta is composed of intermediate mass elements including O, Si and Mg and trace amounts of Fe-peak elements, which disfavours a binary neutron star merger. We discuss three explosion scenarios for SN2018kzr, given the low ejecta mass, intermediate mass element composition and the high likelihood of additional powering - core collapse of an ultra-stripped progenitor, the accretion induced collapse of a white dwarf and the merger of a white dwarf and neutron star. The requirement for an alternative input energy source favours either the accretion induced collapse with magnetar powering or a white dwarf - neutron star merger with energy from disk wind shocks.

The Neil Gehrels Swift Observatory carried out prompt searches for gravitational wave (GW) events detected by the LIGO/Virgo Collaboration (LVC) during the second observing run ("O2"). Swift performed extensive tiling of eight LVC triggers, two of which had very low false-alarm rates (GW 170814 and the epochal GW 170817), indicating a high confidence of being astrophysical in origin; the latter was the first GW event to have an electromagnetic counterpart detected. In this paper we describe the follow-up performed during O2 and the results of our searches. No GW electromagnetic counterparts were detected; this result is expected, as GW 170817 remained the only astrophysical event containing at least one neutron star after LVC's later retraction of some events. A number of X-ray sources were detected, with the majority of identified sources being active galactic nuclei. We discuss the detection rate of transient X-ray sources and their implications in the O2 tiling searches. Finally, we describe the lessons learned during O2, and how these are being used to improve the \swift\ follow-up of GW events. In particular, we simulate a population of GRB afterglows to evaluate our source ranking system's ability to differentiate them from unrelated and uncatalogued X-ray sources. We find that $\approx$60-70% of afterglows whose jets are oriented towards Earth will be given high rank (i.e., "interesting" designation) by the completion of our second follow-up phase (assuming their location in the sky was observed), but that this fraction can be increased to nearly 100% by performing a third follow-up observation of sources exhibiting fading behavior.

We present a high density disc reflection spectral analysis of a sample of 17 Seyfert 1 galaxies to study the inner disc densities at different black hole mass scales and accretion rates. All the available XMM-Newton observations in the archive are used. OM observations in the optical/UV band are used to estimate their accretion rates. We find that 65% of sources in our sample show a disc density significantly higher than 1e15 cm^{-3}, which was assumed in previous reflection-based spectral analyses. The best-fit disc densities show an anti-correlation with black hole mass and mass accretion rate. High density disc reflection model can successfully explain the soft excess emission and significantly reduce inferred iron abundances. We also compare our black hole spin and disc inclination angle measurements with previous analyses.

Einstein's gravity has undergone extensive tests in the weak field gravitational limit, with results in agreement with theoretical predictions. There exist theories beyond general relativity (GR) which modify gravity in the strong field regime but agree with GR in the weak field. Astrophysical black holes are believed to be described by the Kerr metric and serve as suitable candidates to test strong gravity with electromagnetic radiation. We perform such a test by fitting one Suzaku dataset of the narrow-line Seyfert 1 (NLS1) galaxy Mrk 335 with X-ray reflection spectroscopy, using the Johannsen metric to model the black hole spacetime and test for deviations from Kerr. We find the data is best modeled with a hybrid model that includes both partial covering absorption and a reflection component. This is the first time such a model has been proposed for a high-flux (low reflection) Mrk 335 dataset. We constrain the Johannsen deformation parameter $\alpha_{13}$ to $-1.5<\alpha_{13}<0.6$ with spin parameter $a_{*}>0.8$, and the $\alpha_{22}$ parameter to $-0.4<\alpha_{22}<2.1$ with $a_{*}>0.7$, both at the 99% confidence level. Although additional solutions at large deviations from the Kerr metric show statistical similarity with the ones above, further analysis suggests these solutions may be manifestations of uncertainties beyond our control and do not represent the data. Hence, our results are in agreement with the idea that the supermassive compact object at the center of Mrk 335 is described by the Kerr metric.

We present the public release version of relxill_nk, an X-ray reflection model for testing the Kerr hypothesis and general relativity. This model extends the relxill model that assumes the black hole spacetime is described by the Kerr metric. We also present relxilllp_nk, the first non-Kerr X-ray reflection model with a lamppost corona configuration, as well as all other models available in the full relxill_nk package. In all models the relevant relativistic effects are calculated through a general relativistic ray-tracing code that can be applied to any well-behaved, stationary, axisymmetric, and asymptotically flat black hole spacetime. We show that the numerical error introduced by using a ray-tracing code is not significant as compared with the observational error present in current X-ray reflection spectrum observations. In addition, we present the reflection spectrum for the Johannsen metric as calculated by relxill_nk.

We continue exploring the constraining capabilities of X-ray reflection spectroscopy to test the Kerr-nature of astrophysical black holes and we present the results of our analysis of two NuSTAR observations of Cygnus X-1 in the soft state. We find that the final measurement can strongly depend on the assumption of the intensity profile. We conclude that Cygnus X-1 is not suitable for tests of general relativity using X-ray reflection spectroscopy and we discuss the desired properties of a source to be a good candidate for our studies.

In this work, we present a spectral and temporal analysis of {\it Swift}/XRT and {\it NuSTAR} observations of GRS 1716--249 during its recent 2016--2017 outburst. This low mass X-ray binary underwent an extraordinary outburst after a long quiescence of 23 years, since its last major outburst in 1993. The source was observed over two different epochs during 2017 April, 07 and 10. The best fit joint spectral fitting in the energy range 0.5 $-$ 79.0 keV indicates that the spectrum is best described by relatively cold, weak disk blackbody emission, dominant thermal Comptonization emission, and a relativistically broadened fluorescent iron K$\alpha$ emission line. We observed a clear indication of a Compton hump around 30 keV. We also detected an excess feature of $\sim1.3$ keV. Assuming a lamp-post geometry of the corona, we constrained the inner disk radius for both observations to 11.92$^{+8.62}_{-11.92}$ R$_{ISCO}$ (i.e., an upper limit) and 10.39$^{+9.51}_{-3.02}$ R$_{ISCO}$ (where R$_{ISCO}\equiv$ radius of the innermost stable circular orbit) for the first epoch (E1) and second epoch (E2), respectively. A significant ($\sim5\sigma$) type$-$C quasi-periodic oscillation (QPO) at $1.20\pm0.04$ Hz is detected for the first time for GRS 1716--249, which drifts to $1.55\pm0.04$ Hz ($\sim6\sigma$) at the end of the second observation. The derived spectral and temporal properties show a positive correlation between the QPO frequency and the photon index.

We report results from joint NuSTAR, Swift and XMM-Newton observations of the newly discovered black hole X-ray binary candidate Swift J1658.2-4242 in the intermediate state. We observe a peculiar event in this source, with its X-ray flux rapidly decreasing by $\sim$45\% in $\sim$40~s, accompanied by only subtle changes in the shape of the broadband X-ray spectrum. In addition, we find a sudden turn-on of a transient QPO with a frequency of $6-7$~Hz around the time of the flux change, and the total fractional rms amplitude of the power spectrum increases from $\sim$2\% to $\sim$10\%. X-ray spectral and timing analyses indicate that the flux decrease is driven by intrinsic changes in the accretion flow around the black hole, rather than intervening material along the line of sight. In addition, we do not significantly detect any relativistic disk reflection component, indicating it is much weaker than previously observed while the source was in the bright hard state. We propose accretion disk instabilities triggered at a large disk radius as the origin of the fast transition in spectral and timing properties, and discuss possible causes of the unusual properties observed in Swift J1658.2-4242. The prompt flux variation detected along with the emergence of a QPO makes the event an interesting case for investigating QPO mechanisms in black hole X-ray binaries.

X-ray absorption fine structure (XAFS) in the 0.2-2 keV band is a crucial component in multi-wavelength studies of dust mineralogy, size, and shape -- parameters that are necessary for interpreting astronomical observations and building physical models across all fields, from cosmology to exoplanets. Despite its importance, many fundamental questions about dust remain open. What is the origin of the dust that suffuses the interstellar medium (ISM)? Where is the missing interstellar oxygen? How does iron, predominantly produced by Type Ia supernovae, become incorporated into dust? What is the main form of carbon in the ISM, and how does it differ from carbon in stellar winds? The next generation of X-ray observatories, employing microcalorimeter technology and $R \equiv \lambda/\Delta \lambda \geq 3000$ gratings, will provide pivotal insights for these questions by measuring XAFS in absorption and scattering. However, lab measurements of mineralogical candidates for astrophysical dust, with R > 1000, are needed to fully take advantage of the coming observations.

Astrophysical black hole systems are the ideal laboratories for testing Einstein's theory of gravity in the strong field regime. We have recently developed a framework which uses the reflection spectrum of black hole systems to perform precision tests of general relativity by testing the Kerr black hole hypothesis. In this paper, we analyze XMM-Newton and NuSTAR observations of the supermassive black hole in the Seyfert 1 galaxy MCG-06-30-15 with our disk reflection model. We consider the Johannsen metric with the deformation parameters $\alpha_{13}$ and $\alpha_{22}$, which quantify deviations from the Kerr metric. For $\alpha_{22} = 0$, we obtain the black hole spin $0.928 < a_* < 0.983$ and $-0.44 < \alpha_{13} < 0.15$. For $\alpha_{13} = 0$, we obtain $0.885 < a_* < 0.987$ and $-0.12 < \alpha_{22} < 1.05$. The Kerr solution is recovered for $\alpha_{13} = \alpha_{22} = 0$. Thus, our results include the Kerr solution within statistical uncertainties. Systematic uncertainties are difficult to account for, and we discuss some issues in this regard.

We employ the accretion disk reflection model RELXILL_NK to test the spacetime geometry around the stellar-mass black hole in GRS 1915+105. We adopt the Johannsen metric with the deformation parameters $\alpha_{13}$ and $\alpha_{22}$, for which the Kerr solution is recovered when $\alpha_{13} = \alpha_{22} = 0$. We analyze a NuSTAR observation of 2012, obtaining vanishing and non-vanishing values of the deformation parameters depending on the astrophysical model adopted. Similar difficulties were not found in our previous tests with other sources. The results of this work can shed light on the choice of sources suitable for testing the Kerr metric using X-ray reflection spectroscopy and on the parts of our reflection models that more urgently require improvement.

We present the study of a sample of seven "bare" active galactic nuclei (AGN) observed with Suzaku. We interpret the spectrum of these sources with a relativistic reflection component and we employ our model RELXILL_NK to test the Kerr nature of their supermassive black holes. We constrain the Johannsen deformation parameters $\alpha_{13}$ and $\alpha_{22}$, in which the Kerr metric is recovered when $\alpha_{13} = \alpha_{22} = 0$. All our measurements are consistent with the hypothesis that the spacetime geometry around these supermassive objects is described by the Kerr solution. For some sources, we obtain quite strong constraints on $\alpha_{13}$ and $\alpha_{22}$ when compared to those found in our previous studies. We discuss the systematic uncertainties in our tests and the implications of our results.

X-ray reflection spectroscopy is a promising technique for testing general relativity in the strong field regime, as it can be used to test the Kerr black hole hypothesis. In this context, the parametrically deformed black hole metrics proposed by Konoplya, Rezzolla \& Zhidenko (Phys. Rev. D93, 064015, 2016) form an important class of non-Kerr black holes. We implement this class of black hole metrics in \textsc{relxill\_nk}, which is a framework we have developed for testing for non-Kerr black holes using X-ray reflection spectroscopy. We compare the results from this new model in the Kerr limit with existing codes, illustrate the effect of the leading deformation parameters on the X-ray reflection spectrum and present constraints on some KRZ parameters using data from the supermassive black hole in Ark 564.

Thanks to high-resolution and non-dispersive spectrometers onboard future X-ray missions such as XRISM and Athena, we are finally poised to answer important questions about the formation and evolution of galaxies and large-scale structure. However, we currently lack an adequate understanding of many atomic processes behind the spectral features we will soon observe. Large error bars on parameters as critical as transition energies and atomic cross sections can lead to unacceptable uncertainties in the calculations of e.g., elemental abundance, velocity, and temperature. Unless we address these issues, we risk limiting the full scientific potential of these missions. Laboratory astrophysics, which comprises theoretical and experimental studies of the underlying physics behind observable astrophysical processes, is therefore central to the success of these missions.

X-ray reverberation mapping has emerged as a powerful probe of microparsec scales around AGN, and with high sensitivity detectors, its full potential in echo-mapping the otherwise inaccessible disk-corona at the black hole horizon scale will be revealed.

We present the Spectroscopic Time-Resolving Observatory for Broadband Energy X-rays (STROBE-X), a probe-class mission concept selected for study by NASA. It combines huge collecting area, high throughput, broad energy coverage, and excellent spectral and temporal resolution in a single facility. STROBE-X offers an enormous increase in sensitivity for X-ray spectral timing, extending these techniques to extragalactic targets for the first time. It is also an agile mission capable of rapid response to transient events, making it an essential X-ray partner facility in the era of time-domain, multi-wavelength, and multi-messenger astronomy. Optimized for study of the most extreme conditions found in the Universe, its key science objectives include: (1) Robustly measuring mass and spin and mapping inner accretion flows across the black hole mass spectrum, from compact stars to intermediate-mass objects to active galactic nuclei. (2) Mapping out the full mass-radius relation of neutron stars using an ensemble of nearly two dozen rotation-powered pulsars and accreting neutron stars, and hence measuring the equation of state for ultradense matter over a much wider range of densities than explored by NICER. (3) Identifying and studying X-ray counterparts (in the post-Swift era) for multiwavelength and multi-messenger transients in the dynamic sky through cross-correlation with gravitational wave interferometers, neutrino observatories, and high-cadence time-domain surveys in other electromagnetic bands. (4) Continuously surveying the dynamic X-ray sky with a large duty cycle and high time resolution to characterize the behavior of X-ray sources over an unprecedentedly vast range of time scales. STROBE-X's formidable capabilities will also enable a broad portfolio of additional science.

Atoll sources are accreting neutron star (NS) low-mass X-ray binaries. We present a spectral analysis of four persistent atoll sources (GX 3+1, 4U 1702$-$429, 4U 0614+091, and 4U 1746$-$371) observed for $\sim20$ ks each with NuSTAR to determine the extent of the inner accretion disk. These sources range from an apparent luminosity of $0.006-0.11$ of the Eddington limit (assuming the empirical limit of $3.8\times10^{38}$ ergs s$^{-1}$). Broad Fe emission features shaped by Doppler and relativistic effects close to the NS were firmly detected in three of these sources. The position of the disk appears to be close to the innermost stable circular orbit (ISCO) in each case. For GX 3+1, we determine $R_{in}=1.8^{+0.2}_{-0.6}\ R_{\mathrm{ISCO}}$ (90% confidence level) and an inclination of $27^{\circ}-31^{\circ}$. For 4U 1702$-$429, we find a $R_{in}=1.5_{-0.4}^{+1.6}\ R_{\mathrm{ISCO}}$ and inclination of $53^{\circ}-64^{\circ}$. For 4U 0614+091, the disk has a position of $R_{in}=1.3_{-0.2}^{+5.4}\ R_{\mathrm{ISCO}}$ and inclination of $50^{\circ}-62^{\circ}$. If the disk does not extend to the innermost stable circular orbit, we can place conservative limits on the magnetic field strength in these systems in the event that the disk is truncated at the Alfv\'{e}n radius. This provides the limit at the poles of $B\leq6.7\times10^{8}$ G, $3.3\times10^{8}$ G, and $14.5\times10^{8}$ G for GX 3+1, 4U 1702$-$429, and 4U 0614+091, respectively. For 4U 1746$-$371, we argue that the most plausible explanation for the lack of reflection features is a combination of source geometry and strong Comptonization. We place these sources among the larger sample of NSs that have been observed with NuSTAR.

A number of works point to the presence of narrow emission features at unusual energies in the X-ray spectra of active galactic nuclei (AGN) or to potentially low iron abundances as possible evidence for the spallation of iron. With the imminence of high-resolution calorimeter spectroscopy, the potential to test spallation models in astrophysical sources will soon be possible. Previously determined nuclear spallation reactions of Fe are recalculated making use of improved total inelastic and partial reaction cross sections that result in different absolute and relative abundances of the main spallation elements Mn, Cr, V, and Ti. The effects of ionisation and dynamics near the black hole are examined in simulated spectra with CCD and calorimeter (i.e. Hitomi) resolution. With high-resolution, differences in relative abundances and ionisation should be detectable if spallation is originating at large distances from the black hole (e.g. torus or disc wind) where blurring is not significant. If spallation were occurring in the inner accretion disc, it would likely be undetected as blurring effects would cause significant blending of spectral features.

We present results from a deep, coordinated $XMM$-$Newton$+$NuSTAR$ observation of the Seyfert 2 galaxy IRAS 00521-7054. The $NuSTAR$ data provide the first detection of this source in high-energy X-rays ($E > 10$ keV), and the broadband data show this to be a highly complex source which exhibits relativistic reflection from the inner accretion disc, further reprocessing by more distant material, neutral absorption, and evidence for ionised absorption in an extreme, ultrafast outflow ($v_{\rm{out}} \sim 0.4c$). Based on lamppost disc reflection models, we find evidence that the central supermassive black hole is rapidly rotating ($a > 0.77$), consistent with previous estimates from the profile of the relativistic iron line, and that the accretion disc is viewed at a fairly high inclination ($i \sim 59^{\circ}$). Based on extensive simulations, we find the ultrafast outflow is detected at $\sim$4$\sigma$ significance (or greater). We also estimate that the extreme outflow should be sufficient to power galaxy-scale feedback, and may even dominate the energetics of the total output from the system.

A population of globular clusters (GCs) has been recently established by the Fermi-LAT telescope as a new class of GeV $\gamma$-ray sources. Leptons accelerated to TeV energies, in the inner magnetospheres of MSPs or in their wind regions, should produce $\gamma$-rays through the inverse Compton scattering in the dense radiation field from the huge population of stars. We have conducted deep observations of the globular cluster M15 with the MAGIC telescopes and used 165 hrs in order to search for $\gamma$-ray emission. A strong upper limit on the TeV $\gamma$-ray flux $<3.2\times 10^{-13}\mathrm{cm^{-2}s^{-1}}$ above 300 GeV ($<0.26\%$ of the Crab nebula flux) has been obtained. We interpret this limit as a constraint on the efficiency of the acceleration of leptons in the magnetospheres of the MSPs. We constrain the injection rate of relativistic leptons, $\eta_{\rm e}$, from the MSPs magnetospheres and their surrounding. We conclude that $\eta_{\rm e}$ must be lower than expected from the modelling of high energy processes in MSP inner magnetospheres. For leptons accelerated with the power law spectrum in the MSP wind regions, $\eta_{\rm e}$ is constrained to be much lower than derived for the wind regions around classical pulsars. These constraints are valid for the expected range of magnetic field strengths within the GC and for the range of likely energies of leptons injected from the inner magnetospheres, provided that the leptons are not removed from the globular cluster very efficiently due to advection process. We discuss consequences of these constraints for the models of radiation processes around millisecond pulsars.

The HAWC Collaboration released the 2HWC catalog of TeV sources, in which 19 show no association with any known high-energy (HE; E > 10 GeV) or very-high-energy (VHE; E > 300 GeV) sources. This catalog motivated follow-up studies by both the MAGIC and Fermi-LAT observatories with the aim of investigating gamma-ray emission over a broad energy band. In this paper, we report the results from the first joint work between HAWC, MAGIC and Fermi-LAT on three unassociated HAWC sources: 2HWC J2006+341, 2HWC J1907+084* and 2HWC J1852+013*. Although no significant detection was found in the HE and VHE regimes, this investigation shows that a minimum 1 degree extension (at 95% confidence level) and harder spectrum in the GeV than the one extrapolated from HAWC results are required in the case of 2HWC J1852+013*, while a simply minimum extension of 0.16 degrees (at 95% confidence level) can already explain the scenario proposed by HAWC for the remaining sources. Moreover, the hypothesis that these sources are pulsar wind nebulae is also investigated in detail.

The mechanisms producing fast variability of the $\gamma$-ray emission in active galactic nuclei are under debate. The MAGIC telescopes detected a fast very high energy (VHE, E$>100$ GeV) $\gamma$-ray flare from BL Lacertae on 2015 June 15. The flare had a maximum flux of $(1.5\pm 0.3)\times 10^{-10}$ photons cm$^{-2}$ s$^{-1}$ and halving time of $26\pm8$ minutes. The MAGIC observations were triggered by a high state in the optical and high energy (HE, E$>100$ MeV) $\gamma$-ray bands. In this paper we present the MAGIC VHE $\gamma$-ray data together with multiwavelength data from radio, optical, X-rays, and HE $\gamma$ rays from 2015 May 1 to July 31. Well-sampled multiwavelength data allow us to study the variability in detail and compare it to the other epochs when fast VHE $\gamma$-ray flares have been detected from this source. Interestingly, we find that the behaviour in radio, optical, X-rays and HE $\gamma$-rays is very similar to two other observed VHE $\gamma$-ray flares. In particular, also during this flare there was an indication of rotation of the optical polarization angle and of activity at the 43\,GHz core. These repeating patterns indicate a connection between the three events. We also test modelling of the spectral energy distribution, based on constraints from the light curves and VLBA observations, with two different geometrical setups of two-zone inverse Compton models. In addition we model the $\gamma$-ray data with the star-jet interaction model. We find that all of the tested emission models are compatible with the fast VHE $\gamma$-ray flare, but all have some tension with the multiwavelength observations.

SNR G24.7+0.6 is a 9.5 kyrs radio and $\gamma$-ray supernova remnant evolving in a dense medium. In the GeV regime, SNR G24.7+0.6 (3FHL\,J1834.1--0706e/FGES\,J1834.1--0706) shows a hard spectral index ($\Gamma$$\sim$2) up to $200$\,GeV, which makes it a good candidate to be observed with Cherenkov telescopes such as MAGIC. We observed the field of view of \snr\ with the MAGIC telescopes for a total of 31 hours. We detect very high energy $\gamma$-ray emission from an extended source located 0.34\degr\ away from the center of the radio SNR. The new source, named \mgc\ is detected up to 5\,TeV, and its spectrum is well-represented by a power-law function with spectral index of $2.74 \pm 0.08$. The complexity of the region makes the identification of the origin of the very-high energy emission difficult, however the spectral agreement with the LAT source and overlapping position at less than 1.5$\sigma$ point to a common origin. We analysed 8 years of \fermi-LAT data to extend the spectrum of the source down to 60\,MeV. \fermi-LAT and MAGIC spectra overlap within errors and the global broad band spectrum is described by a power-law with exponential cutoff at $1.9\pm0.5$\,TeV. The detected $\gamma$-ray emission can be interpreted as the results of proton-proton interaction between the supernova and the CO-rich surrounding.

We present the analysis of the first NuSTAR observations ($\sim 220$ ks), simultaneous with the last SUZAKU observations ($\sim 50$ ks), of the active galactic nucleus of the bright Seyfert 1 galaxy Mrk 509. The time-averaged spectrum in the $1-79$ keV X-ray band is dominated by a power-law continuum ($\Gamma\sim 1.8-1.9$), a strong soft excess around 1 keV, and signatures of X-ray reflection in the form of Fe K emission ($\sim 6.4$ keV), an Fe K absorption edge ($\sim 7.1$ keV), and a Compton hump due to electron scattering ($\sim 20-30$ keV). We show that these data can be described by two very different prescriptions for the soft excess: a warm ($kT\sim 0.5-1$ keV) and optically thick ($\tau\sim10-20$) Comptonizing corona, or a relativistically blurred ionized reflection spectrum from the inner regions of the accretion disk. While these two scenarios cannot be distinguished based on their fit statistics, we argue that the parameters required by the warm corona model are physically incompatible with the conditions of standard coronae. Detailed photoionization calculations show that even in the most favorable conditions, the warm corona should produce strong absorption in the observed spectrum. On the other hand, while the relativistic reflection model provides a satisfactory description of the data, it also requires extreme parameters, such as maximum black hole spin, a very low and compact hot corona, and a very high density for the inner accretion disk. Deeper observations of this source are thus necessary to confirm the presence of relativistic reflection, and to further understand the nature of its soft excess.

Aims: To investigate the extension of the very-high-energy spectral tail of the Crab pulsar at energies above 400 GeV. Methods: We analyzed $\sim$320 hours of good quality data of Crab with the MAGIC telescope, obtained from February 2007 until April 2014. Results: We report the most energetic pulsed emission ever detected from the Crab pulsar reaching up to 1.5 TeV. The pulse profile shows two narrow peaks synchronized with the ones measured in the GeV energy range. The spectra of the two peaks follow two different power-law functions from 70 GeV up to 1.5 TeV and connect smoothly with the spectra measured above 10 GeV by the Large Area Telescope (LAT) on board of the Fermi satellite. When making a joint fit of the LAT and MAGIC data, above 10 GeV, the photon indices of the spectra differ by 0.5$\pm$0.1. Conclusions: We measured with the MAGIC telescopes the most energetic pulsed photons from a pulsar to date. Such TeV pulsed photons require a parent population of electrons with a Lorentz factor of at least $5\times 10^6$. These results strongly suggest IC scattering off low energy photons as the emission mechanism and a gamma-ray production region in the vicinity of the light cylinder.

We report on observations of the pulsar / Be star binary system PSR J2032+4127 / MT91 213 in the energy range between 100 GeV and 20 TeV with the VERITAS and MAGIC imaging atmospheric Cherenkov telescope arrays. The binary orbit has a period of approximately 50 years, with the most recent periastron occurring on 2017 November 13. Our observations span from 18 months prior to periastron to one month after. A new, point-like, gamma-ray source is detected, coincident with the location of PSR J2032+4127 / MT91 213. The gamma-ray light curve and spectrum are well-characterized over the periastron passage. The flux is variable over at least an order of magnitude, peaking at periastron, thus providing a firm association of the TeV source with the pulsar / Be star system. Observations prior to periastron show a cutoff in the spectrum at an energy around 0.5 TeV. This result adds a new member to the small population of known TeV binaries, and it identifies only the second source of this class in which the nature and properties of the compact object are firmly established. We compare the gamma-ray results with the light curve measured with the X-ray Telescope (XRT) on board the Neil Gehrels \textit{Swift} Observatory and with the predictions of recent theoretical models of the system. We conclude that significant revision of the models is required to explain the details of the emission we have observed, and we discuss the relationship between the binary system and the overlapping steady extended source, TeV J2032+4130.

It is thought that the spacetime metric around astrophysical black holes is well described by the Kerr solution of Einstein's gravity. However, a robust observational evidence of the Kerr nature of these objects is still lacking. Here we fit the X-ray spectrum of the stellar-mass black hole in GS 1354-645 with a disk reflection model beyond Einstein's gravity in order test the Kerr black hole hypothesis. We consider the Johannsen metric with the deformation parameters $\alpha_{13}$ and $\alpha_{22}$. The Kerr metric is recovered for $\alpha_{13} = \alpha_{22} = 0$. For $\alpha_{22} = 0$, our measurements of the black hole spin and of the deformation parameter $\alpha_{13}$ are $a_* > 0.975$ and $-0.34 < \alpha_{13} < 0.16$, respectively. For $\alpha_{13} = 0$, we find $a_* > 0.975$ and $-0.09 < \alpha_{22} < 0.42$. All the reported uncertainties are at 99% of confidence level for two relevant parameters.

For the first time in human history, technologies have matured sufficiently to enable a mission capable of discovering and characterizing habitable planets like Earth orbiting sunlike stars other than the Sun. At the same time, such a platform would enable unique science not possible from ground-based facilities. This science is broad and exciting, ranging from new investigations of our own solar system to a full range of astrophysics disciplines. The Habitable Exoplanet Observatory, or HabEx, is one of four studies currently being undertaken by NASA in preparation for the 2020 Astrophysics Decadal Survey. HabEx has been designed to be the Great Observatory of the 2030s, with community involvement through a competed and funded Guest Observer (GO) program. This interim report describes the HabEx baseline concept, which is a space-based 4-meter diameter telescope mission concept with ultraviolet (UV), optical, and near-infrared (near-IR) imaging and spectroscopy capabilities. More information on HabEx can be found at https://www.jpl.nasa.gov/habex

PKS 1510-089 is a flat spectrum radio quasar strongly variable in the optical and GeV range. We search for low-state VHE gamma-ray emission from PKS 1510-089. We aim to characterize and model the source in a broad-band context, which would provide a baseline over which high states and flares could be better understood. We use daily binned Fermi-LAT flux measurements of PKS 1510-089 to characterize the GeV emission and select the observation periods of MAGIC during low state of activity. For the selected times we compute the average radio, IR, optical, UV, X-ray and gamma-ray emission to construct a low-state spectral energy distribution of the source. The broadband emission is modelled within an External Compton scenario with a stationary emission region through which plasma and magnetic field are flowing. We perform also the emission-model-independent calculations of the maximum absorption in the broad line region (BLR) using two different models. Results. The MAGIC telescopes collected 75 hrs of data during times when the Fermi-LAT flux measured above 1 GeV was below 3x10-8cm-2s-1, which is the threshold adopted for the definition of a low gamma-ray activity state. The data show a strongly significant (9.5{\sigma}) VHE gamma-ray emission at the level of (4.27+-0.61stat)x10-12cm-2s-1 above 150GeV, a factor 80 smaller than the highest flare observed so far from this object. Despite the lower flux, the spectral shape is consistent with earlier detections in the VHE band. The broad-band emission is compatible with the EC scenario assuming a large emission region located beyond the BLR. For the first time the gamma-ray data allow us to place a limit on the location of the emission region during a low gamma-ray state of a FSRQ. For the used model of the BLR, the 95% C.L. on the location of the emission region allows us to place it at the distance >74% of the outer radius of the BLR.

Fast radio bursts (FRBs) are bright flashes observed typically at GHz frequencies with millisecond duration, whose origin is likely extragalactic. Their nature remains mysterious, motivating searches for counterparts at other wavelengths. FRB 121102 is so far the only source known to repeatedly emit FRBs and is associated with a host galaxy at redshift $z \simeq 0.193$. We conducted simultaneous observations of FRB 121102 with the Arecibo and MAGIC telescopes during several epochs in 2016--2017. This allowed searches for millisecond-timescale burst emission in very-high-energy (VHE) gamma rays as well as the optical band. While a total of five FRBs were detected during these observations, no VHE emission was detected, neither of a persistent nature nor burst-like associated with the FRBs. The average integral flux upper limits above 100 GeV at 95% confidence level are $6.6 \times 10^{-12}~\mathrm{photons\ cm^{-2}\ s^{-1}}$ (corresponding to luminosity $L_{\rm VHE} \lesssim 10^{45}~\mathrm{erg\ s^{-1}}$) over the entire observation period, and $1.2 \times 10^{-7}~ \mathrm{photons\ cm^{-2}\ s^{-1}}$ ($L_{\rm VHE} \lesssim 10^{49}~\mathrm{erg\ s^{-1}}$) over the total duration of the five FRBs. We constrain the optical U-band flux to be below 8.6 mJy at 5-$\sigma$ level for 1-ms intervals around the FRB arrival times. A bright burst with U-band flux $29~\mathrm{mJy}$ and duration $\sim 12$ ms was detected 4.3 s before the arrival of one FRB. However, the probability of spuriously detecting such a signal within the sampled time space is 1.5% (2.2 $\sigma$, post-trial), i.e. consistent with the expected background. We discuss the implications of the obtained upper limits for constraining FRB models.

Clusters of galaxies are the largest known gravitationally bound structures in the Universe, with masses around $10^{15}\ M_\odot$, most of it in the form of dark matter. The ground-based Imaging Atmospheric Cherenkov Telescope MAGIC made a deep survey of the Perseus cluster of galaxies using almost 400 h of data recorded between 2009 and 2017. This is the deepest observational campaign so far on a cluster of galaxies in the very high energy range. We search for gamma-ray signals from dark matter particles in the mass range between 200 GeV and 200 TeV decaying into standard model pairs. We apply an analysis optimized for the spectral and morphological features expected from dark matter decays and find no evidence of decaying dark matter. From this, we conclude that dark matter particles have a decay lifetime longer than $\sim10^{26}$~s in all considered channels. Our results improve previous lower limits found by MAGIC and represent the strongest limits on decaying dark matter particles from ground-based gamma-ray instruments.

A multiwavelength campaign was organized to take place between March and July of 2012. Excellent temporal coverage was obtained with more than 25 instruments, including the MAGIC, FACT and VERITAS Cherenkov telescopes, the instruments on board the Swift and Fermi spacecraft, and the telescopes operated by the GASP-WEBT collaboration. Mrk 501 showed a very high energy (VHE) gamma-ray flux above 0.2 TeV of $\sim$0.5 times the Crab Nebula flux (CU) for most of the campaign. The highest activity occurred on 2012 June 9, when the VHE flux was $\sim$3 CU, and the peak of the high-energy spectral component was found to be at $\sim$2 TeV. This study reports very hard X-ray spectra, and the hardest VHE spectra measured to date for Mrk 501. The fractional variability was found to increase with energy, with the highest variability occurring at VHE, and a significant correlation between the X-ray and VHE bands. The unprecedentedly hard X-ray and VHE spectra measured imply that their low- and high-energy components peaked above 5 keV and 0.5 TeV, respectively, during a large fraction of the observing campaign, and hence that Mrk 501 behaved like an extreme high-frequency- peaked blazar (EHBL) throughout the 2012 observing season. This suggests that being an EHBL may not be a permanent characteristic of a blazar, but rather a state which may change over time. The one-zone synchrotron self-Compton (SSC) scenario can successfully describe the segments of the SED where most energy is emitted, with a significant correlation between the electron energy density and the VHE gamma-ray activity, suggesting that most of the variability may be explained by the injection of high-energy electrons. The one-zone SSC scenario used reproduces the behaviour seen between the measured X-ray and VHE gamma-ray fluxes, and predicts that the correlation becomes stronger with increasing energy of the X-rays.

We present a spectral and timing analysis of the newly reported Galactic X-ray transient Swift J1658.2-4242 observed by NuSTAR and Swift. The broad-band X-ray continuum is typical of a black hole binary in the bright hard state, with a photon index of $\Gamma=1.63\pm0.02$ and a low coronal temperature of $kT_{\rm e}=22\pm1$ keV, corresponding to a low spectral cutoff well constrained by NuSTAR. Spectral modeling of the relativistic disk reflection features, consisting of a broad Fe K$\alpha$ line and the Compton reflection hump, reveals that the black hole is rapidly spinning with the spin parameter of $a^{*}>0.96$, and the inner accretion disk is viewed at a high inclination angle of $i=64^{+2}_{-3}{^\circ}$ (statistical errors, 90% confidence). The high inclination is independently confirmed by dips in the light curves, which can be explained by absorbing material located near the disk plane temporarily obscuring the central region. In addition, we detect an absorption line in the NuSTAR spectra centered at $7.03^{+0.04}_{-0.03}$ keV. If associated with ionized Fe K absorption lines, this provides evidence for the presence of outflowing material in the low/hard state of a black hole binary candidate. A timing analysis shows the presence of a type-C QPO in the power spectrum, with the frequency increasing from $\sim0.14$ Hz to $\sim0.21$ Hz during the single NuSTAR exposure. Our analysis reveals that Swift J1658.2-4242 displays characteristics typical for a black hole binary that is viewed at a high inclination angle, making it a good system for studying the accretion geometry in black hole binaries.

A neutrino with energy of $\sim$290 TeV, IceCube-170922A, was detected in coincidence with the BL Lac object TXS~0506+056 during enhanced gamma-ray activity, with chance coincidence being rejected at $\sim 3\sigma$ level. We monitored the object in the very-high-energy (VHE) band with the MAGIC telescopes for $\sim$41 hours from 1.3 to 40.4 days after the neutrino detection. Day-timescale variability is clearly resolved. We interpret the quasi-simultaneous neutrino and broadband electromagnetic observations with a novel one-zone lepto-hadronic model, based on interactions of electrons and protons co-accelerated in the jet with external photons originating from a slow-moving plasma sheath surrounding the faster jet spine. We can reproduce the multiwavelength spectra of TXS 0506+056 with neutrino rate and energy compatible with IceCube-170922A, and with plausible values for the jet power of $\sim 10^{45} - 4 \times 10^{46} {\rm erg \ s^{-1}}$. The steep spectrum observed by MAGIC is concordant with internal $\gamma\gamma$ absorption above a few tens of GeV entailed by photohadronic production of a $\sim$290 TeV neutrino, corroborating a genuine connection between the multi-messenger signals. In contrast to previous predictions of predominantly hadronic emission from neutrino sources, the gamma-rays can be mostly ascribed to inverse Compton up-scattering of external photons by accelerated electrons. The X-ray and VHE bands provide crucial constraints on the emission from both accelerated electrons and protons. We infer that the maximum energy of protons in the jet co-moving frame can be in the range $\sim 10^{14}$ to $10^{18}$ eV.

The MAGIC telescopes observed S2 0109+22 in 2015 July during its flaring activity in high energy gamma rays observed by Fermi-LAT. We analyse the MAGIC data to characterise the very high energy (VHE) gamma-ray emission of S2 0109+22, which belongs to the subclass of intermediate synchrotron peak (ISP) BL Lac objects. We study the multi-frequency emission in order to investigate the source classification. Finally, we compare the source long-term behaviour to other VHE gamma-ray emitting (TeV) blazars. We performed a temporal and spectral analysis of the data centred around the MAGIC interval of observation (MJD 57225-57231). Long-term radio and optical data have also been investigated using the discrete correlation function. The redshift of the source is estimated through optical host-galaxy imaging and also using the amount of VHE gamma-ray absorption. The quasi-simultaneous multi-frequency spectral energy distribution (SED) is modelled with the conventional one-zone synchrotron self-Compton (SSC) model. MAGIC observations resulted in the detection of the source at a significance level of $5.3\,\sigma$. The VHE gamma-ray emission of S2 0109+22 is variable on a daily time scale. VHE gamma-ray luminosity of the source is lower than the average of TeV BL Lacs. The optical polarization, and long-term optical/radio behaviour of the source are different from the general population of TeV blazars. All these findings agree with the classification of the source as an ISP BL Lac object. We estimate the source redshift as $z = 0.36 \pm 0.07$. The SSC parameters describing the SED are rather typical for blazars.

The BL Lac object S5~0716+714, a highly variable blazar, underwent an impressive outburst in January 2015 (Phase A), followed by minor activity in February (Phase B). The MAGIC observations were triggered by the optical flux observed in Phase A, corresponding to the brightest ever reported state of the source in the R-band. The comprehensive dataset collected is investigated in order to shed light on the mechanism of the broadband emission. Multi-wavelength light curves have been studied together with the broadband Spectral Energy Distributions (SEDs). The data set collected spans from radio, optical photometry and polarimetry, X-ray, high-energy (HE, 0.1 GeV < E < 100 GeV) with \textit{Fermi}-LAT to the very-high-energy (VHE, E>100 GeV) with MAGIC. The flaring state of Phase A was detected in all the energy bands, providing for the first time a multi-wavelength sample of simultaneous data from the radio band to the VHE. In the constructed SED the \textit{Swift}-XRT+\textit{NuSTAR} data constrain the transition between the synchrotron and inverse Compton components very accurately, while the second peak is constrained from 0.1~GeV to 600~GeV by \textit{Fermi}+MAGIC data. The broadband SED cannot be described with a one-zone synchrotron self-Compton model as it severely underestimates the optical flux in order to reproduce the X-ray to $\gamma$-ray data. Instead we use a two-zone model. The EVPA shows an unprecedented fast rotation. An estimation of the redshift of the source by combined HE and VHE data provides a value of $z = 0.31 \pm 0.02_{stats} \pm 0.05_{sys}$, confirming the literature value. The data show the VHE emission originating in the entrance and exit of a superluminal knot in and out a recollimation shock in the inner jet. A shock-shock interaction in the jet seems responsible for the observed flares and EVPA swing. This scenario is also consistent with the SED modelling.

We report on the detection of flaring activity from the Fanaroff-Riley I radio galaxy NGC 1275 in very-high-energy (VHE, E $>$ 100 GeV) gamma rays with the MAGIC telescopes. Observations were performed between 2016 September and 2017 February as part of a monitoring program. The brightest outburst with $\sim1.5$ times the Crab Nebula flux above 100 GeV (C.U.) was observed during the night between 2016 December 31 and 2017 January 1 (fifty times higher than the mean previously measured in two observational campaigns between 2009 and 2011). Significant variability of the day-by-day light curve was measured, the shortest flux-doubling time-scales was found to be of $(611\pm101)$ min. The combined spectrum of the MAGIC data during the strongest flare state and simultaneous data from the Fermi-LAT around 2017 January 1 follows a power-law with an exponential cutoff at the energy $(492\pm35)$ GeV. Simultaneous optical flux density measurements in the R-band obtained with the KVA telescope are also presented and the correlation between the optical and gamma-ray emission is investigated. Due to possible internal pair-production, the fast flux variability constrains the Doppler factor to values which are inconsistent with a large viewing angle as observed in the radio band. We investigate different scenarios for the explanation of fast gamma-ray variability, namely emission from: magnetospheric gaps, relativistic blobs propagating in the jet (mini-jets) or external cloud (or star) entering the jet. We find that the only plausible model to account for the luminosities here observed would be the production of gamma rays in a magnetospheric gap around the central black hole only in the eventuality of an enhancement of the magnetic field threading the hole from its equipartition value with the gas pressure in the accretion flow.

A search for tau neutrino induced showers with the MAGIC telescopes is presented. The MAGIC telescopes located at an altitude of 2200 m a.s.l. in the Canary Island of La Palma, can point towards the horizon or a few degrees below across an azimuthal range of about 80 degrees. This provides a possibility to search for air showers induced by tau leptons arising from interactions of tau neutrinos in the Earth crust or the surrounding ocean. In this paper we show how such air showers can be discriminated from the background of very inclined hadronic showers by using Monte Carlo simulations. Taking into account the orography of the site, the point source acceptance and the event rates expected have been calculated for a sample of generic neutrino fluxes from photo-hadronic interactions in AGNs. The analysis of about 30 hours of data taken towards the sealeads to a 90\% C.L. point source limit for tau neutrinos in the energy range from $1.0 \times 10^{15}$ eV to $3.0 \times 10^{18}$ eV of about $E_{\nu_{\tau}}^{2}\times \phi (E_{\nu_{\tau}}) < 2.0 \times 10^{-4}$ GeV cm$^{-2}$ s$^{-1}$ for an assumed power-law neutrino spectrum with spectral index $\gamma$=-2. However, with 300 hours and in case of an optimistic neutrino flare model, limits of the level down to $E_{\nu_{\tau}}^{2}\times \phi (E_{\nu_{\tau}}) < 8.4 \times 10^{-6}$ GeV cm$^{-2}$ s$^{-1}$ can be expected.

We present Neutron Star Interior Composition Explorer (NICER) observations of the neutron star low-mass X-ray binary Serpens X-1 during the early mission phase in 2017. With the high spectral sensitivity and low-energy X-ray passband of NICER, we are able to detect the Fe L line complex in addition to the signature broad, asymmetric Fe K line. We confirm the presence of these lines by comparing the NICER data to archival observations with XMM-Newton/RGS and NuSTAR. Both features originate close to the innermost stable circular orbit (ISCO). When modeling the lines with the relativistic line model RELLINE, we find the Fe L blend requires an inner disk radius of $1.4_{-0.1}^{+0.2}$ $R_{\mathrm{ISCO}}$ and Fe K is at $1.03_{-0.03}^{+0.13}$ $R_{\mathrm{ISCO}}$ (errors quoted at 90%). This corresponds to a position of $17.3_{-1.2}^{+2.5}$ km and $12.7_{-0.4}^{+1.6}$ km for a canonical neutron star mass ($M_{\mathrm{NS}}=1.4\ M_{\odot}$) and dimensionless spin value of $a=0$. Additionally, we employ a new version of the RELXILL model tailored for neutron stars and determine that these features arise from a dense disk and supersolar Fe abundance.

The dwarf spheroidal galaxy Ursa Major II (UMaII) is believed to be one of the most dark-matter dominated systems among the Milky Way satellites and represents a suitable target for indirect dark matter (DM) searches. The MAGIC telescopes carried out a deep observation campaign on UMaII between 2014 and 2016, collecting almost one hundred hours of good-quality data. This campaign enlarges the pool of DM targets observed at very high energy (E$\gtrsim$50GeV) in search for signatures of dark matter annihilation in the wide mass range between $\sim$100 GeV and $\sim$100 TeV. To this end, the data are analyzed with the full likelihood analysis, a method based on the exploitation of the spectral information of the recorded events for an optimal sensitivity to the explored dark matter models. We obtain constraints on the annihilation cross-section for different channels that are among the most robust and stringent achieved so far at the TeV mass scale from observations of dwarf satellite galaxies.

Recently, we have extended the X-ray reflection model relxill to test the spacetime metric in the strong gravitational field of astrophysical black holes. In the present Letter, we employ this extended model to analyze XMM-Newton, NuSTAR, and Swift data of the supermassive black hole in 1H0707-495 and test deviations from a Kerr metric parametrized by the Johannsen deformation parameter $\alpha_{13}$. Our results are consistent with the hypothesis that the spacetime metric around the black hole in 1H0707-495 is described by the Kerr solution.

Here we study a 1-200 keV energy spectrum of the black hole binary Cygnus X-1 taken with NuSTAR and Suzaku. This is the first report of a NuSTAR observation of Cyg X-1 in the intermediate state, and the observation was taken during the part of the binary orbit where absorption due to the companion's stellar wind is minimal. The spectrum includes a multi-temperature thermal disk component, a cutoff power-law component, and relativistic and non-relativistic reflection components. Our initial fits with publicly available constant density reflection models (relxill and reflionx) lead to extremely high iron abundances (>9.96 and 10.6(+1.6)(-0.9) times solar, respectively). Although supersolar iron abundances have been reported previously for Cyg X-1, our measurements are much higher and such variability is almost certainly unphysical. Using a new version of reflionx that we modified to make the electron density a free parameter, we obtain better fits to the spectrum even with solar iron abundances. We report on how the higher density (n_e = (3.98(+0.12)(-0.25))E20 cm-3) impacts other parameters such as the inner radius and inclination of the disk.

The very-high-energy (VHE, $\gtrsim 100$ GeV) $\gamma$-ray MAGIC observations of the blazar S4 0954+65, were triggered by an exceptionally high flux state of emission in the optical. This blazar has a disputed redshift of z=0.368 or z$\geqslant$0.45 and an uncertain classification among blazar subclasses. The exceptional source state described here makes for an excellent opportunity to understand physical processes in the jet of S4 0954+65 and thus contribute to its classification. We investigate the multiwavelength (MWL) light curve and spectral energy distribution (SED) of the S4 0954+65 blazar during an enhanced state in February 2015 and put it in context with possible emission scenarios. We collect photometric data in radio, optical, X-ray, and $\gamma$ ray. We study both the optical polarization and the inner parsec-scale jet behavior with 43 GHz data. Observations with the MAGIC telescopes led to the first detection of S4 0954+65 at VHE. Simultaneous data with Fermi-LAT at high energy $\gamma$ ray\ (HE, 100 MeV < E < 100 GeV) also show a period of increased activity. Imaging at 43 GHz reveals the emergence of a new feature in the radio jet in coincidence with the VHE flare. Simultaneous monitoring of the optical polarization angle reveals a rotation of approximately 100$^\circ$. (...) The broadband spectrum can be modeled with an emission mechanism commonly invoked for flat spectrum radio quasars, i.e. inverse Compton scattering on an external soft photon field from the dust torus, also known as external Compton. The light curve and SED phenomenology is consistent with an interpretation of a blob propagating through a helical structured magnetic field and eventually crossing a standing shock in the jet, a scenario typically applied to flat spectrum radio quasars (FSRQs) and low-frequency peaked BL Lac objects (LBL).

The flat-spectrum radio quasar PKS 1441+25 at a redshift of z = 0.940 is detected between 40 and 250 GeV with a significance of 25.5 {\sigma} using the MAGIC telescopes. Together with the gravitationally lensed blazar QSO B0218+357 (z = 0.944), PKS 1441+25 is the most distant very high energy (VHE) blazar detected to date. The observations were triggered by an outburst in 2015 April seen at GeV energies with the Large Area Telescope on board Fermi. Multi-wavelength observations suggest a subdivision of the high state into two distinct flux states. In the band covered by MAGIC, the variability time scale is estimated to be 6.4 +/- 1.9 days. Modeling the broadband spectral energy distribution with an external Compton model, the location of the emitting region is understood as originating in the jet outside the broad line region (BLR) during the period of high activity, while being partially within the BLR during the period of low (typical) activity. The observed VHE spectrum during the highest activity is used to probe the extragalactic background light at an unprecedented distance scale for ground-based gamma-ray astronomy.

We report on a NuSTAR observation of the recently discovered bright black hole candidate MAXI J1535-571. NuSTAR observed the source on MJD 58003 (five days after the outburst was reported). The spectrum is characteristic of a black hole binary in the hard state. We observe clear disk reflection features, including a broad Fe K$\alpha$ line and a Compton hump peaking around 30 keV. Detailed spectral modeling reveals narrow Fe K$\alpha$ line complex centered around 6.5 keV on top of the strong relativistically broadened Fe K$\alpha$ line. The narrow component is consistent with distant reflection from moderately ionized material. The spectral continuum is well described by a combination of cool thermal disk photons and a Comptonized plasma with the electron temperature $kT_{\rm e}=19.7\pm{0.4}$ keV. An adequate fit can be achieved for the disk reflection features with a self-consistent relativistic reflection model that assumes a lamp-post geometry for the coronal illuminating source. The spectral fitting measures a black hole spin $a>0.84$, inner disk radius $R_{\rm in}<2.01~r_{\rm ISCO}$, and a lamp-post height $h=7.2^{+0.8}_{-2.0} r_{\rm g}$ (statistical errors, 90% confidence), indicating no significant disk truncation and a compact corona. Although the distance and mass of this source are not currently known, this suggests the source was likely in the brighter phases of the hard state during this NuSTAR observation.

We present a spectral and timing study of the NuSTAR and Swift observations of the black hole candidate IGR J17091-3624 in the hard state during its outburst in 2016. Disk reflection is detected in each of the NuSTAR spectra taken in three epochs. Fitting with relativistic reflection models reveals that the accretion disk is truncated during all epochs with $R_{\rm in}>10~r_{\rm g}$, with the data favoring a low disk inclination of $\sim 30^{\circ}-40^{\circ}$. The steepening of the continuum spectra between epochs is accompanied by a decrease in the high energy cut-off: the electron temperature $kT_{\rm e}$ drops from $\sim 64$ keV to $\sim 26$ keV, changing systematically with the source flux. We detect type-C QPOs in the power spectra with frequency varying between 0.131 Hz and 0.327 Hz. In addition, a secondary peak is found in the power spectra centered at about 2.3 times the QPO frequency during all three epochs. The nature of this secondary frequency is uncertain, however a non-harmonic origin is favored. We investigate the evolution of the timing and spectral properties during the rising phase of the outburst and discuss their physical implications.

The reflection spectroscopic model RELXILL is commonly implemented in studying relativistic X-ray reflection from accretion disks around black holes. We present a systematic study of the model's capability to constrain the dimensionless spin and ionization parameters from $\sim$6,000 NuSTAR simulations of a bright X-ray source employing the lamppost geometry. We employ high count spectra to show the limitations in the model without being confused with limitations in signal-to-noise. We find that both parameters are well-recovered at 90% confidence with improving constraints at higher reflection fraction, high spin, and low source height. We test spectra across a broad range - first at 10$^6-$10$^7$ and then $\sim$10$^5$ total source counts across the effective 3-79 keV band of NuSTAR, and discover a strong dependence of the results on how fits are performed around the starting parameters, owing to the complexity of the model itself. A blind fit chosen over an approach that carries some estimates of the actual parameter values can lead to significantly worse recovery of model parameters. We further stress on the importance to span the space of nonlinear-behaving parameters like $log~\xi$ carefully and thoroughly for the model to avoid misleading results. In light of selecting fitting procedures, we recall the necessity to pay attention to the choice of data binning and fit statistics used to test the goodness of fit by demonstrating the effect on the photon index $\Gamma$. We re-emphasize and implore the need to account for the detector resolution while binning X-ray data and using Poisson fit statistics instead while analyzing Poissonian data.

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

MAGIC (Major Atmospheric Gamma Imaging Cherenkov) is a system of two 17 m diameter, F/1.03 Imaging Atmospheric Cherenkov Telescopes (IACT). They are dedicated to the observation of gamma rays from galactic and extragalactic sources in the very high energy range (VHE, 30 GeV to 100 TeV). This submission contains links to the proceedings for the 35th International Cosmic Ray Conference (ICRC2017), held in Bexco, Busan, Korea from the 12th to the 17th of July, 2017.

Spontaneous breaking of Lorentz symmetry at energies on the order of the Planck energy or lower is predicted by many quantum gravity theories, implying non-trivial dispersion relations for the photon in vacuum. Consequently, gamma-rays of different energies, emitted simultaneously from astrophysical sources, could accumulate measurable differences in their time of flight until they reach the Earth. Such tests have been carried out in the past using fast variations of gamma-ray flux from pulsars, and more recently from active galactic nuclei and gamma-ray bursts. We present new constraints studying the gamma-ray emission of the galactic Crab Pulsar, recently observed up to TeV energies by the MAGIC collaboration. A profile likelihood analysis of pulsar events reconstructed for energies above 400GeV finds no significant variation in arrival time as their energy increases. Ninety-five percent~CL limits are obtained on the effective Lorentz invariance violating energy scale at the level of $E_{\mathrm{QG}_1} > 5.5\cdot 10^{17}$GeV ($4.5\cdot 10^{17}$GeV) for a linear, and $E_{\mathrm{QG}_2} > 5.9\cdot 10^{10}$GeV ($5.3\cdot 10^{10}$GeV) for a quadratic scenario, for the subluminal and the superluminal cases, respectively. A substantial part of this study is dedicated to calibration of the test statistic, with respect to bias and coverage properties. Moreover, the limits take into account systematic uncertainties, found to worsen the statistical limits by about 36--42\%. Our constraints would have resulted much more competitive if the intrinsic pulse shape of the pulsar between 200GeV and 400GeV was understood in sufficient detail and allowed inclusion of events well below 400GeV.

The microquasar Cygnus X-1 displays the two typical soft and hard X-ray states of a black-hole transient. During the latter, Cygnus X-1 shows a one-sided relativistic radio-jet. Recent detection of the system in the high energy (HE; $E\gtrsim60$ MeV) gamma-ray range with \textit{Fermi}-LAT associates this emission with the outflow. Former MAGIC observations revealed a hint of flaring activity in the very high-energy (VHE; $E\gtrsim100$ GeV) regime during this X-ray state. We analyze $\sim97$ hr of Cygnus X-1 data taken with the MAGIC telescopes between July 2007 and October 2014. To shed light on the correlation between hard X-ray and VHE gamma rays as previously suggested, we study each main X-ray state separately. We perform an orbital phase-folded analysis to look for variability in the VHE band. Additionally, to place this variability behavior in a multiwavelength context, we compare our results with \textit{Fermi}-LAT, \textit{AGILE}, \textit{Swift}-BAT, \textit{MAXI}, \textit{RXTE}-ASM, AMI and RATAN-600 data. We do not detect Cygnus X-1 in the VHE regime. We establish upper limits for each X-ray state, assuming a power-law distribution with photon index $\Gamma=3.2$. For steady emission in the hard and soft X-ray states, we set integral upper limits at 95\% confidence level for energies above 200 GeV at $2.6\times10^{-12}$~photons cm$^{-2}$s$^{-1}$ and $1.0\times10^{-11}$~photons cm$^{-2}$s$^{-1}$, respectively. We rule out steady VHE gamma-ray emission above this energy range, at the level of the MAGIC sensitivity, originating in the interaction between the relativistic jet and the surrounding medium, while the emission above this flux level produced inside the binary still remains a valid possibility.

MAGIC, a system of two imaging atmospheric Cherenkov telescopes, achieves its best performance under dark conditions, i.e. in absence of moonlight or twilight. Since operating the telescopes only during dark time would severely limit the duty cycle, observations are also performed when the Moon is present in the sky. Here we develop a dedicated Moon-adapted analysis to characterize the performance of MAGIC under moonlight. We evaluate energy threshold, angular resolution and sensitivity of MAGIC under different background light levels, based on Crab Nebula observations and tuned Monte Carlo simulations. This study includes observations taken under non-standard hardware configurations, such as reducing the camera photomultiplier tubes gain by a factor ~1.7 (Reduced HV settings) with respect to standard settings (Nominal HV) or using UV-pass filters to strongly reduce the amount of moonlight reaching the cameras of the telescopes. The Crab Nebula spectrum is correctly reconstructed in all the studied illumination levels, that reach up to 30 times brighter than under dark conditions. The main effect of moonlight is an increase in the analysis energy threshold and in the systematic uncertainties on the flux normalization. The sensitivity degradation is constrained to be below 10%, within 15-30% and between 60 and 80% for Nominal HV, Reduced HV and UV-pass filter observations, respectively. No worsening of the angular resolution was found. Thanks to observations during moonlight, the maximal duty cycle of MAGIC can be increased from ~18%, under dark nights only, to up to ~40% in total with only moderate performance degradation.

We perform an X-ray spectral analysis of the brightest and cleanest bare AGN known so far, Ark 120, in order to determine the process(es) at work in the vicinity of the SMBH. We present spectral analysis of data from an extensive campaign observing Ark 120 in X-rays with XMM-Newton (4$\times$120 ks, 2014 March 18-24), and NuSTAR (65.5 ks, 2014 March 22). During this very deep X-ray campaign, the source was caught in a high flux state similar to the earlier 2003 XMM-Newton observation, and about twice as bright as the lower-flux observation in 2013. The spectral analysis confirms the "softer when brighter" behaviour of Ark 120. The four XMM-Newton/pn spectra are characterized by the presence of a prominent soft X-ray excess and a significant FeK$\alpha$ complex. The continuum is very similar above about 3 keV, while significant variability is present for the soft X-ray excess. We find that relativistic reflection from a constant-density, flat accretion disk cannot simultaneously produce the soft excess, broad FeK$\alpha$ complex, and hard X-ray excess. Instead, Comptonization reproduces the broadband (0.3-79 keV) continuum well, together with a contribution from a mildly relativistic disk reflection spectrum. During this 2014 observational campaign, the soft X-ray spectrum of Ark 120 below $\sim$0.5 keV was found to be dominated by Comptonization of seed photons from the disk by a warm ($kT_{\rm e}$$\sim$0.5 keV), optically-thick corona ($\tau$$\sim$9). Above this energy, the X-ray spectrum becomes dominated by Comptonization from electrons in a hot optically thin corona, while the broad FeK$\alpha$ line and the mild Compton hump result from reflection off the disk at several tens of gravitational radii.

The large jet kinetic power and non-thermal processes occurring in the microquasar SS 433 make this source a good candidate for a very high-energy (VHE) gamma-ray emitter. Gamma-ray fluxes have been predicted for both the central binary and the interaction regions between jets and surrounding nebula. Also, non-thermal emission at lower energies has been previously reported. We explore the capability of SS 433 to emit VHE gamma rays during periods in which the expected flux attenuation due to periodic eclipses and precession of the circumstellar disk periodically covering the central binary system is expected to be at its minimum. The eastern and western SS433/W50 interaction regions are also examined. We aim to constrain some theoretical models previously developed for this system. We made use of dedicated observations from MAGIC and H.E.S.S. from 2006 to 2011 which were combined for the first time and accounted for a total effective observation time of 16.5 h. Gamma-ray attenuation does not affect the jet/medium interaction regions. The analysis of a larger data set amounting to 40-80 h, depending on the region, was employed. No evidence of VHE gamma-ray emission was found. Upper limits were computed for the combined data set. We place constraints on the particle acceleration fraction at the inner jet regions and on the physics of the jet/medium interactions. Our findings suggest that the fraction of the jet kinetic power transferred to relativistic protons must be relatively small to explain the lack of TeV and neutrino emission from the central system. At the SS433/W50 interface, the presence of magnetic fields greater 10$\mu$G is derived assuming a synchrotron origin for the observed X-ray emission. This also implies the presence of high-energy electrons with energies up to 50 TeV, preventing an efficient production of gamma-ray fluxes in these interaction regions.

It is widely believed that the bulk of the Galactic cosmic rays are accelerated in supernova remnants (SNRs). However, no observational evidence of the presence of particles of PeV energies in SNRs has yet been found. The young historical SNR Cassiopeia A (Cas A) appears as one of the best candidates to study acceleration processes. Between December 2014 and October 2016 we observed Cas A with the MAGIC telescopes, accumulating 158 hours of good-quality data. We derived the spectrum of the source from 100 GeV to 10 TeV. We also analysed $\sim$8 years of $Fermi$-LAT to obtain the spectral shape between 60 MeV and 500 GeV. The spectra measured by the LAT and MAGIC telescopes are compatible within the errors and show a clear turn off (4.6 $\sigma$) at the highest energies, which can be described with an exponential cut-off at $E_c = 3.5\left(^{+1.6}_{-1.0}\right)_{\textit{stat}} \left(^{+0.8}_{-0.9}\right)_{\textit{sys}}$ TeV. The gamma-ray emission from 60 MeV to 10 TeV can be attributed to a population of high-energy protons with spectral index $\sim$2.2 and energy cut-off at $\sim$10 TeV. This result indicates that Cas A is not contributing to the high energy ($\sim$PeV) cosmic-ray sea in a significant manner at the present moment. A one-zone leptonic model fails to reproduce by itself the multi-wavelength spectral energy distribution. Besides, if a non-negligible fraction of the flux seen by MAGIC is produced by leptons, the radiation should be emitted in a region with a low magnetic field (B$\lessapprox$100$\mu$G) like in the reverse shock.)

The microquasar V404 Cygni underwent a series of outbursts in 2015, June 15-31, during which its flux in hard X-rays (20-40 keV) reached about 40 times the Crab Nebula flux. Because of the exceptional interest of the flaring activity from this source, observations at several wavelengths were conducted. The MAGIC telescopes, triggered by the INTEGRAL alerts, followed-up the flaring source for several nights during the period June 18-27, for more than 10 hours. One hour of observation was conducted simultaneously to a giant 22 GHz radio flare and a hint of signal at GeV energies seen by Fermi-LAT. The MAGIC observations did not show significant emission in any of the analysed time intervals. The derived flux upper limit, in the energy range 200--1250 GeV, is 4.8$\times 10^{-12}$ ph cm$^{-2}$ s$^{-1}$. We estimate the gamma-ray opacity during the flaring period, which along with our non-detection, points to an inefficient acceleration in the V404\,Cyg jets if VHE emitter is located further than $1\times 10^{10}$ cm from the compact object.

B1957+20 is a millisecond pulsar located in a black widow type compact binary system with a low mass stellar companion. The interaction of the pulsar wind with the companion star wind and/or the interstellar plasma is expected to create plausible conditions for acceleration of electrons to TeV energies and subsequent production of very high energy {\gamma} rays in the inverse Compton process. We performed extensive observations with the MAGIC telescopes of B1957+20. We interpret results in the framework of a few different models, namely emission from the vicinity of the millisecond pulsar, the interaction of the pulsar and stellar companion wind region, or bow shock nebula. No significant steady very high energy {\gamma}-ray emission was found. We derived a 95% confidence level upper limit of 3.0 x 10 -12 cm -2 s -1 on the average {\gamma}-ray emission from the binary system above 200 GeV. The upper limits obtained with MAGIC constrain, for the first time, different models of the high-energy emission in B1957+20. In particular, in the inner mixed wind nebula model with mono-energetic injection of electrons, the acceleration efficiency of electrons is constrained to be below ~(2-10)% of the pulsar spin down power. For the pulsar emission, the obtained upper limits for each emission peak are well above the exponential cut-off fits to the Fermi-LAT data, extrapolated to energies above 50 GeV. The MAGIC upper limits can rule out a simple power-law tail extension through the sub-TeV energy range for the main peak seen at radio frequencies.

We present the first X-ray reflection model for testing the assumption that the metric of astrophysical black holes is described by the Kerr solution. We employ the formalism of the transfer function proposed by Cunningham. The calculations of the reflection spectrum of a thin accretion disk are split into two parts: the calculation of the transfer function and the calculation of the local spectrum at any emission point in the disk. The transfer function only depends on the background metric and takes into account all the relativistic effects (gravitational redshift, Doppler boosting, and light bending). Our code computes the transfer function for a spacetime described by the Johannsen metric and can easily be extended to any stationary, axisymmetric, and asymptotically flat spacetime. Transfer functions and single line shapes in the Kerr metric are compared with those calculated from existing codes to check that we reach the necessary accuracy. We also simulate some observations with NuSTAR and LAD/eXTP and fit the data with our new model to show the potential capabilities of current and future observations to constrain possible deviations from the Kerr metric.

The relativistic jets created by some active galactic nuclei are important agents of AGN feedback. In spite of this, our understanding of what produces these jets is still incomplete. X-ray observations, which can probe the processes operating in the central regions in immediate vicinity of the supermassive black hole, the presumed jet launching point, are potentially particularly valuable in illuminating the jet formation process. Here, we present the hard X-ray NuSTAR observations of the radio-loud quasar 4C 74.26 in a joint analysis with quasi-simultaneous, soft X-ray Swift observations. Our spectral analysis reveals a high-energy cut-off of 183$_{-35}^{+51}$ keV and confirms the presence of ionized reflection in the source. From the average spectrum we detect that the accretion disk is mildly recessed with an inner radius of $R_\mathrm{in}=4-180\,R_\mathrm{g}$. However, no significant evolution of the inner radius is seen during the three months covered by our NuSTAR campaign. This lack of variation could mean that the jet formation in this radio-loud quasar differs from what is observed in broad-line radio galaxies.

Ark 564 is an archetypal Narrow line Seyfert 1 that has been well observed in soft X-rays from 0.3-10 keV, revealing a steep spectrum, strong soft excess, iron K emission line and dramatic variability on the order of hours. Because of its very steep spectrum, observations of the source above 10 keV have been sparse. We report here on the first NuSTAR observation of Ark 564. The source was observed for 200 ks with NuSTAR, 50 ks of which were concurrent with Suzaku observations. NuSTAR and Suzaku observed a dramatic flare, in which the hard emission is clearly delayed with respect to the soft emission, consistent with previous detections of a low-frequency hard lag found in XMM-Newton data. The NuSTAR spectrum is well described by a low-temperature Comptonisation continuum (with an electron temperature of 15 +/- 2 keV), which irradiates a highly ionised disc. No further relativistic broadening or ionized absorption is required. These spectral results show that Ark 564 has one of the lowest temperature coronae observed by NuSTAR to date. We discuss possible reasons for low-temperature coronae in high-Eddington sources.

The extragalactic VHE gamma-ray sky is rich in blazars. These are jetted active galactic nuclei viewed at a small angle to the line-of-sight. Only a handful of objects viewed at a larger angle are known so far to emit above 100 GeV. Multi-wavelength studies of such objects up to the highest energies provide new insights into the particle and radiation processes of active galactic nuclei. We report the results from the first multi-wavelength campaign observing the TeV detected nucleus of the active galaxy IC 310, whose jet is observed at a moderate viewing angle of 10 deg - 20 deg. The multi-instrument campaign was conducted between 2012 Nov. and 2013 Jan., and involved observations with MAGIC, Fermi, INTEGRAL, Swift, OVRO, MOJAVE and EVN. These observations were complemented with archival data from the AllWISE and 2MASS catalogs. A one-zone synchrotron self-Compton model was applied to describe the broad-band spectral energy distribution. IC 310 showed an extraordinary TeV flare at the beginning of the campaign, followed by a low, but still detectable TeV flux. Compared to previous measurements, the spectral shape was found to be steeper during the low emission state. Simultaneous observations in the soft X-ray band showed an enhanced energy flux state and a harder-when-brighter spectral shape behaviour. No strong correlated flux variability was found in other frequency regimes. The broad-band spectral energy distribution obtained from these observations supports the hypothesis of a double-hump structure. The harder-when-brighter trend in the X-ray and VHE emission is consistent with the behaviour expected from a synchrotron self-Compton scenario. The contemporaneous broad-band spectral energy distribution is well described with a one-zone synchrotron self-Compton model using parameters that are comparable to those found for other gamma-ray-emitting misaligned blazars.

We present first results from a series of $NuSTAR$ observations of the black hole X-ray binary V404 Cyg obtained during its summer 2015 outburst, primarily focusing on observations during the height of this outburst activity. The $NuSTAR$ data show extreme variability in both the flux and spectral properties of the source. This is partly driven by strong and variable line-of-sight absorption, similar to previous outbursts. The latter stages of this observation are dominated by strong flares, reaching luminosities close to Eddington. During these flares, the central source appears to be relatively unobscured and the data show clear evidence for a strong contribution from relativistic reflection, providing a means to probe the geometry of the innermost accretion flow. Based on the flare properties, analogy with other Galactic black hole binaries, and also the simultaneous onset of radio activity, we argue that this intense X-ray flaring is related to transient jet activity during which the ejected plasma is the primary source of illumination for the accretion disk. If this is the case, then our reflection modelling implies that these jets are launched in close proximity to the black hole (as close as a few gravitational radii), consistent with expectations for jet launching models that tap either the spin of the central black hole, or the very innermost accretion disk. Our analysis also allows us to place the first constraints on the black hole spin for this source, which we find to be $a^* > 0.92$ (99% statistical uncertainty, based on an idealized lamppost geometry).

Context. PKS 1510-089 is one of only a few flat spectrum radio quasars detected in the VHE (very-high-energy, > 100 GeV) gamma-ray band. Aims. We study the broadband spectral and temporal properties of the PKS 1510-089 emission during a high gamma-ray state. Methods. We performed VHE gamma-ray observations of PKS 1510-089 with the MAGIC telescopes during a long high gamma-ray state in May 2015. In order to perform broadband modelling of the source, we have also gathered contemporaneous multiwavelength data in radio, IR, optical photometry and polarization, UV, X-ray and GeV gamma-ray ranges. We construct a broadband spectral energy distribution (SED) in two periods, selected according to VHE gamma-ray state. Results. PKS 1510-089 has been detected by MAGIC during a few day-long observations performed in the middle of a long, high optical and gamma-ray state, showing for the first time a significant VHE gamma-ray variability. Similarly to the optical and gamma-ray high state of the source detected in 2012, it was accompanied by a rotation of the optical polarization angle and the emission of a new jet component observed in radio. However, due to large uncertainty on the knot separation time, the association with the VHE gamma-ray emission cannot be firmly established. The spectral shape in the VHE band during the flare is similar to the ones obtained during previous measurements of the source. The observed flux variability sets for the first time constraints on the size of the region from which VHE gamma rays are emitted. We model the broadband SED in the framework of the external Compton scenario and discuss the possible emission site in view of multiwavelength data and alternative emission models.

Active galactic nucleus (AGN) feedback is the process by which supermassive black holes in the centres of galaxies may moderate the growth of their hosts. Gas outflows from supermassive black holes release huge quantities of energy into the interstellar medium, clearing the surrounding gas. The most extreme of these, the ultra-fast outflows (UFOs), are the subset of X-ray detected outflows with velocities higher than 10,000 km/s, believed to originate in relativistic disc winds, a few hundred gravitational radii from the black hole. The absorption features produced by these outflows are variable, but no clear link has been found between the behaviour of the X-ray continuum and the energy or equivalent width of the outflow features due to the long time-scales of quasar variability. Here, we present the detection of multiple absorption lines from an extreme ultra-fast gas flow in the X-ray spectrum of the active galactic nucleus IRAS 13224-3809, at 0.236+/-0.006 times the speed of light (71,000 km/s), where the absorption is strongly anti-correlated with the emission from the inner regions of the accretion disk. If the gas flow is identified as a genuine outflow then it is in the fastest 5 per cent of such winds, and its variability is hundreds of times faster than in other variable winds, allowing us to observe in hours what would take months in a quasar. We find signatures of the wind simultaneously in both low and high energy detectors, which are consistent with a single ionized outflow, linking the two phenomena. The detection of the wind responding to the emission from the inner disk demonstrates a connection between accretion processes occurring on very different scales, with the X-rays from within a few gravitational radii of the black hole ionizing the relativistically outflowing gas as the flux rises.

In this work we present data from observations with the MAGIC telescopes of SN 2014J detected in January 21 2014, the closest Type Ia supernova since Imaging Air Cherenkov Telescopes started to operate. We probe the possibility of very-high-energy (VHE; $E\geq100$ GeV) gamma rays produced in the early stages of Type Ia supernova explosions. We performed follow-up observations after this supernova explosion for 5 days, between January 27 and February 2 in 2014. We search for gamma-ray signal in the energy range between 100 GeV and several TeV from the location of SN 2014J using data from a total of $\sim5.5$ hours of observations. Prospects for observing gamma-rays of hadronic origin from SN 2014J in the near future are also being addressed. No significant excess was detected from the direction of SN 2014J. Upper limits at 95$\%$ confidence level on the integral flux, assuming a power-law spectrum, d$F/$d$E\propto E^{-\Gamma}$, with a spectral index of $\Gamma=2.6$, for energies higher than 300 GeV and 700 GeV, are established at $1.3\times10^{-12}$ and $4.1\times10^{-13}$ photons~cm$^{-2}$s$^{-1}$, respectively. For the first time, upper limits on the VHE emission of a Type Ia supernova are established. The energy fraction isotropically emitted into TeV gamma rays during the first $\sim10$ days after the supernova explosion for energies greater than 300 GeV is limited to $10^{-6}$ of the total available energy budget ($\sim 10^{51}$ erg). Within the assumed theoretical scenario, the MAGIC upper limits on the VHE emission suggest that SN 2014J will not be detectable in the future by any current or planned generation of Imaging Atmospheric Cherenkov Telescopes.

We present the first detection of the nearby (z=0.084) low-luminosity BL Lac object 1ES 1741+196 in the very high energy (VHE: E$>$100 GeV) band. This object lies in a triplet of interacting galaxies. Early predictions had suggested 1ES 1741+196 to be, along with several other high-frequency BL Lac sources, within the reach of MAGIC detectability. Its detection by MAGIC, later confirmed by VERITAS, helps to expand the small population of known TeV BL Lacs. The source was observed with the MAGIC telescopes between 2010 April and 2011 May, collecting 46 h of good quality data. These observations led to the detection of the source at 6.0 $\sigma$ confidence level, with a steady flux $\mathrm{F}(> 100 {\rm GeV}) = (6.4 \pm 1.7_{\mathrm{stat}}\pm 2.6_{\mathrm{syst}}) \cdot 10^{-12}$ ph cm$^{-2}$ s$^{-1}$ and a differential spectral photon index $\Gamma = 2.4 \pm 0.2_{\mathrm{stat}} \pm 0.2_{\mathrm{syst}}$ in the range of $\sim$80 GeV - 3 TeV. To study the broad-band spectral energy distribution (SED) simultaneous with MAGIC observations, we use KVA, Swift/UVOT and XRT, and Fermi/LAT data. One-zone synchrotron-self-Compton (SSC) modeling of the SED of 1ES 1741+196 suggests values for the SSC parameters that are quite common among known TeV BL Lacs except for a relatively low Doppler factor and slope of electron energy distribution. A thermal feature seen in the SED is well matched by a giant elliptical's template. This appears to be the signature of thermal emission from the host galaxy, which is clearly resolved in optical observations.

47 Tuc X9 is a low mass X-ray binary (LMXB) in the globular cluster 47 Tucanae, and was previously thought to be a cataclysmic variable. However, Miller-Jones et al. (2015) recently identified a radio counterpart to X9 (inferring a radio/X-ray luminosity ratio consistent with black hole LMXBs), and suggested that the donor star might be a white dwarf. We report simultaneous observations of X9 performed by Chandra, Nustar and Australia Telescope Compact Array. We find a clear 28.18$\pm$0.02 min periodic modulation in the Chandra data, which we identify as the orbital period, confirming this system as an ultracompact X-ray binary. Our X-ray spectral fitting provides evidence for photoionized gas having a high oxygen abundance in this system, which indicates a C/O white dwarf donor. We also identify reflection features in the hard X-ray spectrum, making X9 the faintest LMXB to show X-ray reflection. We detect a $\sim$ 6.8 day modulation in the X-ray brightness by a factor of 10, in archival Chandra, Swift, and Rosat data. The simultaneous radio/X-ray flux ratio is consistent with either a black hole primary or a neutron star primary, if the neutron star is a transitional millisecond pulsar. Considering the measured orbital period (with other evidence of a white dwarf donor), and the lack of transitional millisecond pulsar features in the X-ray light curve, we suggest that this could be the first ultracompact black hole X-ray binary identified in our Galaxy.

We present a spectral analysis of the NuSTAR and XMM-Newton observations of the Seyfert 1h galaxy/ULIRG IRAS 05189-2524 taken in 2013. We find evidence for relativistic disk reflection in the broadband X-ray spectrum: a highly asymmetric broad Fe K$\alpha$ emission line extending down to 3 keV and a Compton scattering component above 10 keV. Physical modeling with a self-consistent disk reflection model suggests the accretion disk is viewed at an intermediate angle with a super-solar iron abundance, and a mild constraint can be put on the high-energy cutoff of the power-law continuum. We test the disk reflection modeling under different absorption scenarios. A rapid black hole (BH) spin is favored, however we cannot place a model-independent tight constraint on the value. The high reflection fraction ($R_{\rm ref} \simeq$ 2.0-3.2) suggests the coronal illuminating source is compact and close to the BH (lying within 8.7 $R_{\rm g}$ above the central BH), where light-bending effects are important.

Continuum and reflection spectral models have each been widely employed in measuring the spins of accreting black holes. However, the two approaches have not been implemented together in a photon-conserving, self-consistent framework. We develop such a framework using the black-hole X-ray binary GX 339-4 as a touchstone source, and we demonstrate three important ramifications: (1) Compton scattering of reflection emission in the corona is routinely ignored, but is an essential consideration given that reflection is linked to the regimes with strongest Comptonization. Properly accounting for this causes the inferred reflection fraction to increase substantially, especially for the hard state. Another important impact of the Comptonization of reflection emission by the corona is the downscattered tail. Downscattering has the potential to mimic the relativistically-broadened red wing of the Fe line associated with a spinning black hole. (2) Recent evidence for a reflection component with a harder spectral index than the power-law continuum is naturally explained as Compton-scattered reflection emission. (3) Photon conservation provides an important constraint on the hard state's accretion rate. For bright hard states, we show that disk truncation to large scales (R much larger than the ISCO radius) is unlikely as this would require accretion rates far in excess of the brightest soft states. Our principal conclusion is that when modeling relativistically-broadened reflection, spectral models should allow for coronal Compton scattering of the reflection features, and when possible, take advantage of the additional constraining power from linking to the thermal disk component.

MCG-5-23-16 was targeted in early 2015 with a half mega-seconds observing campaign using NuSTAR. Here we present the spectral analysis of these datasets along with an earlier observation and study the relativistic reflection and the primary coronal source. The data show strong reflection features in the form of both narrow and broad iron lines plus a Compton reflection hump. A cutoff energy is significantly detected in all exposures. The shape of the reflection spectrum does not change in the two years spanned by the observations, suggesting a stable geometry. A strong positive correlation is found between the cutoff energy and both the hard X-ray flux and spectral index. The measurements imply that the coronal plasma is not at the runaway electron-positron pair limit, and instead contains mostly electrons. The observed variability in the coronal properties is driven by a variable optical depth. A constant heating to cooling ratio is measured implying that there is a feedback mechanism in which a significant fraction of the photons cooling the corona are due to reprocessed hard X-rays.