The central engine in long gamma-ray bursts (GRBs) is thought to be a compact object produced by the core collapse of massive stars, but its exact nature (black hole or millisecond magnetar) is still debatable. Although the central engine of GRB collapsars is hidden to direct observation, its properties may be imprinted on the accompanying electromagnetic signals. We aim to decipher the generic properties of central engines that are consistent with prompt observations of long GRBs detected by the Burst Alert Telescope (BAT) on board the Neil Gehrels Swift Observatory. Adopting a generic model for the central engine, in which the engine power and activity timescale are independent of each other, we perform Monte Carlo simulations of long GRBs produced by jets that successfully breakout from the star. Our simulations consider the dependence of the jet breakout timescale on the engine luminosity and the effects of the detector's flux threshold. The two-dimensional (2D) distribution of simulated detectable bursts in the gamma-ray luminosity versus gamma-ray duration plane is consistent with the observed one for a range of parameter values describing the central engine. The intrinsic 2D distribution of simulated collapsar GRBs peaks at lower gamma-ray luminosities and longer durations than the observed one, a prediction that can be tested in the future with more sensitive detectors. Black-hole accretors, whose power and activity time are set by the large-scale magnetic flux through the progenitor star and stellar structure, respectively, are compatible with the properties of the central engine inferred by our model.

3HSP J095507.9+355101 is an extreme blazar to have been possibly associated with a high-energy neutrino, IceCube-200107A, which was detected one day before the blazar was observed and found to undergo a hard X-ray flare. Motivated by this observation, we perform a comprehensive study of the predicted multimessenger emission from 3HSP J095507.9+355101 during its recent X-ray flare, but also in the long term. We focus on single-zone leptohadronic models, where the broadband photon and high-energy neutrino emissions originate from same region of the jet, but we also explore alternative scenarios: (i) a blazar-core model, which considers neutrino production in the inner jet close to the accreting supermassive black hole, (ii) a hidden external-photon model, which considers neutrino production in the jet through interactions with photons from a possible weak broad line region, (iii) a one-zone proton synchrotron model, where high-energy protons produce $\gamma$-rays in the jet via synchrotron, and (iv) an intergalactic cascade scenario, where neutrinos are produced in the intergalactic medium by interactions of a high-energy cosmic-ray beam escaping the blazar jet. We find that the Poisson probability to detect one or more neutrinos in 8 years from 3HSP J095507.9+355101 is $\sim 2\%$ ($3\%$) when considering the most optimistic one-zone leptohadronic model (the multi-zone blazar-core model), while detection of one neutrino event during the X-ray flare is much less likely. If the association is real, then IceCube-Gen2 and other future detectors should be able to provide additional evidence for neutrino production in 3HSP J095507.9+355101 and other extreme blazars, assuming similar physical conditions.

In November 2019, MAXI detected an X-ray outburst from the known Be X-ray binary system RX J0209.6-7427 located in the outer wing of the Small Magellanic Cloud. We followed the outburst of the system with NICER which led to the discovery of X-ray pulsations with a period of 9.3 s. We analyzed simultaneous X-ray data obtained with NuSTAR and NICER allowing us to characterize the spectrum and provide an accurate estimate of its bolometric luminosity. During the outburst the maximum broadband X-ray luminosity of the system reached $1-2\times10^{39}$ erg/s, thus exceeding by about one order of magnitude the Eddington limit for a typical 1.4 $M_{\odot}$ mass neutron star (NS). Monitoring observations with Fermi/GBM and NICER allowed us to study the spin evolution of the NS and compare it with standard accretion torque models. We found that the NS magnetic field should be of the order of $3\times10^{12}$ G. We conclude that RX J0209.6-7427 exhibited one of the brightest outbursts observed from a Be X-ray binary pulsar in the Magellanic Clouds, reaching similar luminosity level to the 2016 outburst of SMC X-3. Despite the super-Eddington luminosity of RX J0209.6-7427, the NS appears to have only a moderate magnetic field strength.

We present results of the first dedicated observation of the young X-ray pulsar SXP 1062 in the broad X-ray energy band obtained during its 2019 outburst with the NuSTAR and XMM-Newton observatories. The analysis of the pulse-phase averaged and phase-resolved spectra in the energy band from 0.5 to 70 keV did not reveal any evidence for the presence of a cyclotron line. The spin period of the pulsar was found to have decreased to 979.48+/-0.06 s implying a ~10% reduction compared to the last measured period during the monitoring campaign conducted about five years ago, and is puzzling considering that the system apparently did not show major outbursts ever since. The switch of the pulsar to the spin-up regime supports the common assumption that torques acting on the accreting neutron star are nearly balanced and thus SXP 1062 likely also spins with a period close to the equilibrium value for this system. The current monitoring of the source revealed also a sharp drop of its soft X-ray flux right after the outburst, which is in drastic contrast to the behaviour during the previous outburst when the pulsar remained observable for years with only a minor flux decrease after the end of the outburst. This unexpected off state of the source lasted for at most 20 days after which SXP 1062 returned to the level observed during previous campaigns. We discuss this and other findings in context of the modern models of accretion onto strongly magnetized neutron stars.

The IceCube report of a $\sim 3.5\sigma$ excess of $13\pm5$ neutrino events in the direction of the blazar TXS 05056+056 in 2014-2015 and the 2017 detection of a high-energy neutrino, IceCube-170922A, during a gamma-ray flare from the same blazar, have revived the interest in scenarios for neutrino production in blazars. We perform comprehensive analyses on the long-term electromagnetic emission of TXS 05056+056 using optical, X-ray, and gamma-ray data from the All-Sky Automated Survey for Supernovae (ASAS-SN), the Neil Gehrels Swift Observatory (Swift), the Monitor of All-sky X-ray Image (MAXI), and the Fermi Large Area Telescope (Fermi-LAT). We also perform numerical modeling of the spectral energy distributions (SEDs) in four epochs prior to 2017 with contemporaneous gamma-ray and lower energy (optical and/or X-ray) data. We find that the multi-epoch SEDs are consistent with a hybrid leptonic scenario, where the gamma-rays are produced in the blazar zone via external inverse Compton scattering of accelerated electrons, and high-energy neutrinos are produced via the photomeson production process of co-accelerated protons. The multi-epoch SEDs can be satisfactorily explained with the same jet parameters and variable external photon density and electron luminosity. Using the maximal neutrino flux derived for each epoch, we put an upper limit of $\sim0.4-2$ on the muon neutrino number in ten years of IceCube observations. Our results are consistent with the IceCube-170922A detection, which can be explained as an upper fluctuation from the average neutrino rate expected from the source, but in strong tension with the 2014-2015 neutrino flare.

Recent modeling of Neutron Star Interior Composition Explorer(NICER) observations of the millisecond pulsar PSR J0030+0451 suggests that the magnetic field of the pulsar is non-dipolar. We construct a magnetic field configuration where foot points of the open field lines closely resemble the hotspot configuration from NICER observations. Using this magnetic field as input, we perform force-free simulations of the magnetosphere of PSR J0030+0451, showing the three-dimensional structure of its plasma-filled magnetosphere. Making simple and physically motivated assumptions about the emitting regions, we are able to construct the multi-wavelength lightcurves that qualitatively agree with the corresponding observations. The agreement suggests that multipole magnetic structures are the key to modeling this type of pulsars, and can be used to constrain the magnetic inclination angle and the location of radio emission.

Aims: Spectral and temporal analysis of the NuSTAR observation Galactic Be-XRB Swift J1845.7-0037. during its recent outburst. Methods: For the spectral analysis we use both phenomenological and physics-based models. We employ an often used empirical model to identify the main characteristics of the spectral shape in relation to nominal spectral characteristics of X-ray pulsars. Additionally, we used the latest version of Bulk \& Thermal comptonization model (BW), to assess the validity of the spectral components required by the empirical model and to investigate the origin of the hard X-ray emission. We also analyzed the source light-curve, studying the pulse shape at different energy ranges and tracking the spectral evolution with pulse phase by using the model independent hardness ratio (HR). Results: We find that while both the empirical and physical (BW) spectral models can produce good spectral fits, the BW model returns physically plausible best-fit values for the source parameters and does not require any additional spectral components to the non-thermal, accretion column emission. The BW model also yielded an estimation of the neutron star magnetic field placing it in the 10^12G range. Conclusions: Our results, show that the spectral and temporal characteristics of the source emission are consistent with the scattering processes expected for radiation dominated shocks within the accretion column of highly magnetized accreting neutron stars. We further indicate that physically-derived spectral models such as BW, can be used to tentatively infer fundamental source parameters, in the absence of more direct observational signatures.

We investigate the long-term variability of the K$\alpha$ line of iron in the spectra of two Ultra Compact X-ray Sources (UCXBs) with C/O-rich donors. We revisit archival observations by five different X-ray telescopes, over a ~twenty year period. Adopting physically motivated models for the spectral continuum, we probe the long-term evolution of the source emission in a self-consistent manner enabling physical interpretation of potential variability in the primary X-ray emission continuum and/or any emission lines from reflection off the accretion disk. We find that the spectral shape and flux of the source emission (for both objects) has remained almost constant throughout all the observations, displaying only minor variability in some spectral parameters and the source flux (largest variation is a ~25% drop in the flux of Swift J1756.9-2508). We note a striking variability of the Fe K$\alpha$ line which fluctuates from a notable equivalent width of ~66-100 eV in 4U 1543-624 and ~170 eV in Swift J1756.9-2508, to non-detections with upper limits of 2-8 eV. We argue that the disappearance of the iron line is due to the screening of the Fe K$\alpha$ line by the overabundant oxygen in the C/O-rich UCXBs. This effect is cancelled when oxygen becomes fully ionized in the inner disk region, resulting in the variability of the Fe K$\alpha$ line in an otherwise unaltered spectral shape. This finding supports earlier predictions on the consequences of H-poor, C/O-rich accretion disk on reflection induced fluorescent lines in the spectra of UCXBs.

We have analyzed the X-ray spectra of all known Ultra Compact X-ray Binaries (UCXBs), with the purpose of constraining the chemical composition of their accretion disk and donor star. Our investigation was focused on the presence (or absence) of the Fe K${\alpha}$ emission line, which was used as the probe of chemical composition of the disk, based on previously established theoretical predictions for the reflection of X-ray radiation off the surface of C/O-rich or He-rich accretion disks in UCXBs. We have contrasted the results of our spectral analysis to the history of type I X-ray bursts from these systems, which can also indicate donor star composition. We found that UCXBs with prominent and persistent iron K${\alpha}$ emission also featured repeat bursting activity. On the other hand, the UCXBs for which no iron line was detected, appear to have few or no type I X-ray bursts detected over more than a decade of monitoring. Based on Monte Carlo simulations, demonstrating a strong correlation between the Fe K${\alpha}$ line strength and the abundance of C and O in the accretion disk material and given the expected correlation between the H/He abundance and the recurrence rate of type I X-ray bursts, we propose that there is a considerable likelihood that UCXBs with persistent iron emission have He-rich donors, while those that do not, likely have C/O or O/Ne/Mg-rich donors. Our results strongly advocate for the development of more sophisticated simulations of X-ray reflection from hydrogen-poor accretion disks.

In the current work we explore the applicability of standard theoretical models of accretion to the observed properties of M51 ULX-7. The spin-up rate and observed X-ray luminosity are evidence of a neutron star with a surface magnetic field of $2-7\times10^{13}$ G, rotating near equilibrium. Analysis of the X-ray light-curve of the system (Swift/XRT data) reveals the presence of a $\sim$39 d super-orbital period. We argue that the super-orbital periodicity is due to disc precession, and that material is accreted onto the neutron star at a constant rate throughout it. Moreover, by attributing this modulation to the free precession of the neutron star we estimate a surface magnetic field strength of $3-4\times10^{13}$ G. The agreement of these two independent estimates provide strong constraints on the surface polar magnetic field strength of the neutron star.

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

NGC300 ULX1 is an ultra-luminous X-ray pulsar, showing an unprecedented spin evolution, from about 126 s to less than 20 s in only 4 years, consistent with steady mass accretion rate. Following its discovery we have been monitoring the system with Swift/XRT and NICER to further study its properties. We found that even though the observed flux of the system dropped by a factor of $\gtrsim$20, the spin-up rate remained almost constant. A possible explanation is that the decrease in the observed flux is a result of increased absorption of obscuring material due to outflows or a precessing accretion disk.

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.

NGC 300 ULX1 is a newly identified ultra-luminous X-ray pulsar. The system is associated with the supernova impostor SN 2010da that was later classified as a possible supergiant Be X-ray binary. In this work we report on the spin period evolution of the neutron star based on all the currently available X-ray observations of the system. We argue that the X-ray luminosity of the system has remained almost constant since 2010, at a level above ten times the Eddington limit. Moreover, we find evidence that the spin period of the neutron star evolved from ~126 s down to ~18 s within a period of about 4 years. We explain this unprecedented spin evolution in terms of the standard accretion torque theory. An intriguing consequence for NGC 300 ULX1 is that a neutron star spin reversal should have occurred a few years after the SN 2010da event.

Methods. We analyzed broadband XMM-Newton and NuSTAR observations of NGC 300 ULX1, performing phase-averaged and phase-resolved spectroscopy. We compared two physically motivated models for the source spectrum: Non-thermal accretion column emission modeled by a power law with high-energy exponential roll-off (AC model) vs multicolor thermal emission from an optically thick accretion envelope plus a hard power-law tail (MCAE model). We combine the findings of our Bayesian analysis with qualitative physical considerations to evaluate the suitability of each model. Results. The low-energy part (<2 keV) of the source spectrum is dominated by non-pulsating, multicolor thermal emission. The (pulsating) high energy continuum is more ambiguous. If modelled with the AC model a residual structure is detected, that can be modeled using a broad Gaussian absorption line centered at ~12 keV. However, the same residuals can be successfully modeled using the MCAE model, without the need for the absorption-like feature. Model comparison, using the Bayesian approach strongly indicates that the MCAE model -- without the absorption line -- is the preferred model. Conclusions. All models considered strongly indicate the presence of an accretion disk truncated at a large distance from the central object, as has been recently suggested for a large fraction of both pulsating and non-pulsating ULXs. The hard, pulsed emission is not described by a smooth spectral continuum. If modelled by a broad Gaussian absorption line, the fit residuals can be interpreted as a cyclotron scattering feature (CRSF) compatible with a ~10^12 G magnetic field. However, the MCAE model describes the spectral and temporal characteristics of the source emission, without the need for an absorption feature and yields physically meaningful parameter values. Therefore strong doubts are cast on the presence of a CRSF in NGC 300 ULX1.

NGC 300 ULX1 is the fourth to be discovered in the class of the ultra-luminous X-ray pulsars. Pulsations from NGC 300 ULX1 were discovered during simultaneous XMM-Newton / NuSTAR observations in Dec. 2016. The period decreased from 31.71 s to 31.54 s within a few days, with a spin-up rate of -5.56 x 10^{-7} s s^{-1}, likely one of the largest ever observed from an accreting neutron star. Archival Swift and NICER observations revealed that the period decreased exponentially from ~45 s to ~17.5 s over 2.3 years. The pulses are highly modulated with a pulsed fraction strongly increasing with energy and reaching nearly 80% at energies above 10keV. The X-ray spectrum is described by a power-law and a disk black-body model, leading to a 0.3-30 keV unabsorbed luminosity of 4.7 x 10^{39} erg s^{-1}. The spectrum from an archival XMM-Newton observation of 2010 can be explained by the same model, however, with much higher absorption. This suggests, that the intrinsic luminosity did not change much since that epoch. NGC 300 ULX1 shares many properties with supergiant high mass X-ray binaries, however, at an extreme accretion rate.

B[e] Supergiants are luminous evolved massive stars. The mass-loss during this phase creates a complex circumstellar environment with atomic, molecular, and dusty regions usually found in rings or disk-like structures. For a better comprehension of the mechanisms behind the formation of these rings, detailed knowledge about their structure and dynamics is essential. To address that, we obtained high-resolution optical and near-infrared spectra for 8 selected Galactic B[e] Supergiants, for which CO emission has been detected. Assuming Keplerian rotation for the disk, we combine the kinematics obtained from the CO bands in the near-IR with those obtained by fitting the forbidden emission [OI] $\lambda$5577, [OI] $\lambda\lambda$6300,6363, and [CaII] $\lambda\lambda$7291,7323 lines in the optical to probe the disk structure. We find that the emission originates from multiple ring structures around all B[e] Supergiants, with each one of them displaying a unique combination of rings regardless of whether the object is part of a binary system. The confirmed binaries display spectroscopic variations of their line intensities and profiles as well as photometric variability, whereas the ring structures around the single stars are stable.

We report the results of AstroSat and NuSTAR observations of the Be/X-ray binary pulsar SXP 15.3 in the Small Magellanic Cloud during its outburst in late 2017, when the source reached a luminosity level of ~ 10^{38} erg s^{-1}, close to the Eddington limit. The unprecedented broadband coverage of the source allowed us to perform timing and spectral analysis between 3 and 80 keV. The pulse profile exhibits a significant energy dependence, and morphs from a double peaked profile to a single broad pulse at energies >15 keV. This can be explained by a spectral hardening during an intensity dip seen between the two peaks of the pulse profile. We detect a Cyclotron Resonance Scattering Feature (CRSF) at ~5 keV in the X-ray spectrum, independent of the choice of the continuum model. This indicates a magnetic field strength of 6x10^{11} G for the neutron star.

We study the temporal and spectral characteristics of SMC X-3 during its 2016 outburst. To probe accretion onto highly magnetized neutron stars (NSs), at the Eddington limit. We obtained XMM-Newton observations of SMC X-3 and combined them with long-term observations by Swift. We studied the temporal and spectral behavior of the source, and its short- and long-term evolution. We constructed a simple toy-model to gain insight into the complex emission pattern of SMC X-3. We confirm the pulse period of the system, derived by previous works and note that the pulse has a complex, three peaked shape. We find that the pulsed emission is dominated by hard photons, while at energies below ~1 keV, the emission is virtually non-pulsating. We further find that the shape of the pulse profile and the short and long-term evolution of the source light-curve can be explained by invoking a combination of a "fan" and a "polar" beam. Our spectroscopic analysis, reveals a two-component emission: a hard power law and a soft thermal component. We find that the latter produces the bulk of the non-pulsating emission and is most likely the result of reprocessing of the primary hard emission, by optically thick material that partially obscures the central source. We also detect strong emission lines from highly ionized metals. The strength of the emission lines are strongly phase-depended. The energy and temporal evolution and the shape of the pulse profile and the long-term spectra evolution of the source are consistent with the expected emission pattern of the accretion column in the super-critical regime, while the presence of a large reprocessing region is consistent with the analysis of previously studied X-ray pulsars observed at high accretion rates. The presence of this region is consistent with recently proposed works that suggest that high-B NSs occupy a considerable fraction of ULXs.

The supernova impostor SN 2010da located in the nearby galaxy NGC 300, later identified as a likely supergiant B[e] high-mass X-ray binary, was simultaneously observed by NuSTAR and XMM-Newton between 2016 December 16 and 20, over a total time span of 310 ks. We report the discovery of a strong periodic modulation in the X-ray flux with a pulse period of 31.6 s and a very rapid spin-up, and confirm therefore that the compact object is a neutron star. We find that the spin period is changing from 31.71 s to 31.54 s over that period, with a spin-up rate of -5.56 x 10-7 s s-1, likely the largest ever observed from an accreting neutron star. The spectrum is described by a power-law and a disk black-body model, leading to a 0.3-30 keV unabsorbed luminosity of 4.7 x 10^39 erg s-1. Applying our best-fit model successfully to the spectra of an XMM-Newton observation from 2010, suggests that the lower fluxes of NGC 300 ULX1 reported from observations around that time are caused by a large amount of absorption, while the intrinsic luminosity was similar as seen in 2016. A more constant luminosity level is also consistent with the long-term pulse period evolution approaching an equilibrium value asymptotically. We conclude that the source is another candidate for the new class of ultraluminous X-ray pulsars.

The Large Magellanic Cloud (LMC) currently hosts around 23 high mass X-ray binaries (HMXBs) of which most are Be/X-ray binaries. The LMC XMM-Newton survey provided follow-up observations of previously known X-ray sources that were likely HMXBs, as well as identifying new HMXB candidates. In total 19 candidate HMXBs were selected based on their X-ray hardness ratios. In this paper we present red and blue optical spectroscopy, obtained with SALT and the SAAO 1.9-m telescope, plus a timing analysis of the long term optical light curves from OGLE to confirm the nature of these candidates. We find that 9 of the candidates are new Be/X-ray Binaries, substantially increasing the LMC Be/X-ray binary population. Furthermore, we present the optical properties of these new systems, both individually and as a group of all the BeXBs identified by the XMM-Newton survey of the LMC.

We report on the X-ray and optical properties of two high-mass X-ray binary systems located in the Large Magellanic Cloud (LMC). Based on the obtained optical spectra, we classify the massive companion as a supergiant star in both systems. Timing analysis of the X-ray events collected by XMM-Newton revealed the presence of coherent pulsations (spin period $\sim$2013 s) for XMMU J053108.3-690923 and fast flaring behaviour for XMMU J053320.8-684122. The X-ray spectra of both systems can be modelled sufficiently well by an absorbed power-law, yielding hard spectra and high intrinsic absorption from the environment of the systems. Due to their combined X-ray and optical properties we classify both systems as SgXRBs: the 19$^{\rm th}$ confirmed X-ray pulsar and a probable supergiant fast X-ray transient in the LMC, the second such candidate outside our Galaxy.

Aims: In light of recent discoveries of pulsating ULXs and recently introduced models, that propose neutron stars (NSs) as the central engines of ULXs, we revisit the (XMM-Newton and NuSTAR) spectra of eighteen well known ULXs, in search of indications that favor this new hypothesis. Results: We confirm the, previously noted, presence of the low energy (<6 keV) spectral rollover and argue that it could be interpreted as thermal emission. The spectra are well described by a double thermal model consisting of a "hot" (>1 keV) and a "cool" (<0.7 keV) multicolor black body (MCB). Under the assumption that the "cool" MCD emission originates in a disk truncated at the NS magnetosphere, we find that all ULXs in our sample are consistent with accretion onto a highly magnetized (B>10^12G) NS. We note a strong correlation between the strength of the magnetic field, the temperature of the "hot" thermal component and the total unabsorbed luminosity. Examination of the NuSTAR data supported this interpretation and confirmed the presence of a weak, high energy (>15 keV) tail, most likely the result of modification of the MCB emission by inverse Compton scattering. We also note that the apparent high energy tail, may simply be the result of mis-modeling of MCB emission with an atypical temperature (T) vs radius (r) gradient, using a standard MCD model with a fixed gradient of T~r^-0.75. Conclusions: We have offered a new physical interpretation for the dual-thermal spectra of ULXs. We find that the best-fit derived parameters, are in excellent agreement with recent predictions that favor super-critically accreting NSs as the engines of a large fraction of ULXs. Nevertheless, the considerable degeneracy between models and the lack of unequivocal evidence cannot rule out other equally plausible interpretations. Our results warrant further exploration of this new framework.

PSR J2032+4127 is a young and rapidly rotating pulsar on a highly eccentric orbit around the high-mass Be star MT91 213. X-ray monitoring of the binary system over a $\sim$4000 day period with Swift has revealed an increase of the X-ray luminosity which we attribute to the synchrotron emission of the shocked pulsar wind. We use Swift X-ray observations to infer a clumpy stellar wind with $r^{-2}$ density profile and constrain the Lorentz factor of the pulsar wind to $10^5<\gamma_w<10^6$. We investigate the effects of an axisymmetric stellar wind with polar gradient on the X-ray emission. Comparison of the X-ray light curve hundreds of days before and after the periastron can be used to explore the polar structure of the wind.

Type IIn supernovae (SNe), a rare subclass of core collapse SNe, explode in dense circumstellar media that have been modified by the SNe progenitors at their last evolutionary stages. The interaction of the freely expanding SN ejecta with the circumstellar medium gives rise to a shock wave propagating in the dense SN environment, which may accelerate protons to multi-PeV energies. Inelastic proton-proton collisions between the shock-accelerated protons and those of the circumstellar medium lead to multi-messenger signatures. Here, we evaluate the possible neutrino signal of type IIn SNe and compare with IceCube observations. We employ a Monte Carlo method for the calculation of the diffuse neutrino emission from the SN IIn class to account for the spread in their properties. The cumulative neutrino emission is found to be ~ 10 per cent of the observed IceCube neutrino flux above 60 TeV. Type IIn SNe would be the dominant component of the diffuse astrophysical flux, only if 4 per cent of all core collapse SNe were of this type and 20 to 30 per cent of the shock energy was channeled to accelerated protons. Lower values of the acceleration efficiency are accessible by the observation of a single type IIn SN as a neutrino point source with IceCube using up-going muon neutrinos. Such an identification is possible in the first year following the SN shock breakout for sources within 20 Mpc.

The EXTraS project to explore the X-ray Transient and variable Sky searches for coherent signals in the X-ray archival data of XMM-Newton. XMM-Newton performed more than 400 pointed observations in the region of the Large Magellanic Cloud (LMC). We inspected the results of the EXTraS period search to systematically look for new X-ray pulsators in our neighbour galaxy. We analysed the XMM-Newton observations of two sources from the 3XMM catalogue which show significant signals for coherent pulsations. 3XMM J051259.8-682640 was detected as source with hard X-ray spectrum in two XMM-Newton observations, revealing a periodic modulation of the X-ray flux with 956~s. As optical counterpart we identify an early-type star with Halpha emission. The OGLE I-band light curve exhibits a regular pattern with three brightness dips which mark a period of ~1350 d. The X-ray spectrum of 3XMM J051034.6-670353 is dominated by a super-soft blackbody-like emission component (kT ~ 70 eV) which is modulated by nearly 100% with a period of ~1418 s. From GROND observations we suggest a star with r' = 20.9 mag as possible counterpart of the X-ray source. 3XMM J051259.8-682640 is confirmed as a new Be/X-ray binary pulsar in the Large Magellanic Cloud. We discuss the long-term optical period as likely orbital period which would be the longest known from a high-mass X-ray binary. The spectral and temporal properties of the super-soft source 3XMM J051034.6-670353 are very similar to those of RX J0806.3+1527 and RX J1914.4+2456 suggesting that it belongs to the class of double-degenerate polars and is located in our Galaxy rather than in the LMC.

Ap Librae is one out of a handful of low-frequency peaked blazars to be detected at TeV $\gamma$-rays and the only one with an identified X-ray jet. Combined observations of Fermi-LAT at high energies (HE) and of H.E.S.S at very high energies (VHE) revealed a striking spectral property of Ap Librae; the presence of a broad high-energy component that extends more than nine orders of magnitude in energy and is, therefore, hard to be explained by the usual single-zone synchrotron self-Compton model. We show that the superposition of different emission components related to photohadronic interactions can explain the $\gamma$-ray emission of Ap Librae without invoking external radiation fields. We present two indicative model fits to the spectral energy distribution of Ap Librae where the VHE emission is assumed to originate from a compact, sub-pc scale region of the jet. A robust prediction of our model is VHE flux variability on timescales similar to those observed at X-rays and HE $\gamma$-rays, which can be further used to distinguish between a sub-pc or kpc scale origin of the TeV emission. We thus calculate the expected variability signatures at X-rays, HE and VHE $\gamma$-rays and show that quasi-simultaneous flares are expected, with larger amplitude flares appearing at $\gamma$-rays. We assess the detectability of VHE variability from Ap Librae with CTA, next generation of IACTs. We show that $\sim$hr timescale variability at $E_{\gamma}>0.1$ TeV could be detectable at high significance with shorter exposure times than current Cherenkov telescopes.

We report on the results of a $\sim$40 d multi-wavelength monitoring of the Be X-ray binary system IGR J05007-7047 (LXP 38.55). During that period the system was monitored in the X-rays using the Swift telescope and in the optical with multiple instruments. When the X-ray luminosity exceeded $10^{36}$ erg/s we triggered an XMM-Newton ToO observation. Timing analysis of the photon events collected during the XMM-Newton observation reveals coherent X-ray pulsations with a period of 38.551(3) s (1 {\sigma}), making it the 17$^{th}$ known high-mass X-ray binary pulsar in the LMC. During the outburst, the X-ray spectrum is fitted best with a model composed of an absorbed power law ($\Gamma =0.63$) plus a high-temperature black-body (kT $\sim$ 2 keV) component. By analysing $\sim$12 yr of available OGLE optical data we derived a 30.776(5) d optical period, confirming the previously reported X-ray period of the system as its orbital period. During our X-ray monitoring the system showed limited optical variability while its IR flux varied in phase with the X-ray luminosity, which implies the presence of a disk-like component adding cooler light to the spectral energy distribution of the system.

We report on the first detection of X-ray dust scattered rings from the Galactic low mass X-ray binary V404 Cyg. The observation of the system with Swift/XRT on June 30 2015 revealed the presence of five concentric ring-like structures centred at the position of V404 Cyg. Follow-up Swift/XRT observations allowed a time-dependent study of the X-ray rings. Assuming that these are the result of small-angle, single X-ray scattering by dust grains along the line of sight, we find that their angular size scales as $\theta \propto\sqrt{t}$ in agreement with theoretical predictions. The dust grains are concentrated in five dust layers located at about 2.12, 2.05, 1.63, 1.50 and 1.18 kpc from the observer. These coincide roughly with locations of enhanced extinction as determined by infrared photometry. Assuming that the grain size distribution is described by a generalized Mathis-Rumpl-Nordsieck model, we find that the power-law index of the most distant cloud is $q\sim 4.4$, while $q \sim 3.5-3.7$ in all other clouds. We constrain at a $3\sigma$ level the maximum grain size of the intermediate dust layers in the range $0.16-0.20\,\mu$m and set a lower limit of $\sim 0.2\,\mu$m in the other clouds. Hints of an exponential cutoff at the angular intensity profile of the outermost X-ray ring suggest that the smallest grains have sizes $0.01 \mu{\rm m}\le \alpha_{\min} \lesssim 0.03\,\mu$m. Based on the relative ratios of dust column densities we find the highest dust concentration at $\sim 1.6$ kpc. Our results indicate a gradient in the dust properties within 1 kpc from V404 Cyg.

We present a comprehensive X-ray study of the population of supernova remnants (SNRs) in the LMC. Using primarily XMM-Newton, we conduct a systematic spectral analysis of LMC SNRs to gain new insights on their evolution and the interplay with their host galaxy. We combined all the archival XMM observations of the LMC with those of our Very Large Programme survey. We produced X-ray images and spectra of 51 SNRs, out of a list of 59. Using a careful modelling of the background, we consistently analysed all the X-ray spectra and measure temperatures, luminosities, and chemical compositions. We investigated the spatial distribution of SNRs in the LMC and the connection with their environment, characterised by various SFHs. We tentatively typed all LMC SNRs to constrain the ratio of core-collapse to type Ia SN rates in the LMC. We compared the X-ray-derived column densities to HI maps to probe the three-dimensional structure of the LMC. This work provides the first homogeneous catalogue of X-ray spectral properties of LMC SNRs. It offers a complete census of LMC SNRs exhibiting Fe K lines (13% of the sample), or revealing contribution from hot SN ejecta (39%). Abundances in the LMC ISM are found to be 0.2-0.5 solar, with a lower [$\alpha$/Fe] than in the Milky Way. The ratio of CC/type Ia SN in the LMC is $N_{\mathrm{CC}}/N_{\mathrm{Ia}} = 1.35(_{-0.24}^{+0.11})$, lower than in local SN surveys and galaxy clusters. Comparison of X-ray luminosity functions of SNRs in Local Group galaxies reveals an intriguing excess of bright objects in the LMC. We confirm that 30 Doradus and the LMC Bar are offset from the main disc of the LMC, to the far and near sides, respectively. (abridged)

We identify a new candidate for a Be/X-ray binary in the XMM-Newton slew survey and archival Swift observations that is located in the transition region of the Wing of the Small Magellanic Cloud and the Magellanic Bridge. We investigated and classified this source with follow-up XMM-Newton and optical observations. We model the X-ray spectra and search for periodicities and variability in the X-ray observations and the OGLE I-band light curve. The optical counterpart has been classified spectroscopically, with data obtained at the SAAO 1.9 m telescope, and photometrically, with data obtained using GROND at the MPG 2.2 m telescope. The X-ray spectrum is typical of a high-mass X-ray binary with an accreting neutron star. We detect X-ray pulsations, which reveal a neutron-star spin period of P = (264.516+-0.014) s. The source likely shows a persistent X-ray luminosity of a few 10^35 erg/s and in addition type-I outbursts that indicate an orbital period of ~146 d. A periodicity of 0.867 d, found in the optical light curve, can be explained by non-radial pulsations of the Be star. We identify the optical counterpart and classify it as a B1-2II-IVe star. This confirms SXP 265 as a new Be/X-ray binary pulsar originating in the tidal structure between the Magellanic Clouds.

We present new optical spectroscopy of 20 candidate counterparts of 17 X-ray sources in the direction of the M31 disc. By comparing the X-ray catalogue from the XMM-Newton survey of M31 with star catalogues from the Local Group Galaxy Survey, we chose counterpart candidates based on optical colour and X-ray hardness. We have discovered 17 counterpart candidates with spectra containing stellar features. Eight of these are early-type stars of O or B type in M31, with hard X-ray spectra, making them good HMXB candidates. Three of these eight exhibit emission lines, which we consider to be the strongest HMXB candidates. In addition, our spectra reveal two likely Galactic cataclysmic variables, one foreground M star, two probable LMXBs related to M31 globular clusters, one emission line region with an embedded Wolf-Rayet star, and one newly-discovered supernova remnant. Finally, two of the sources have stellar spectra with no features indicative of association with an X-ray source.

We observed RX J0520.5-6932 in the X-rays and studied the optical light curve of its counterpart to verify it as a Be/X-ray binary. We performed an XMM-Newton anticipated target of opportunity observation in January 2013 during an X-ray outburst of the source in order to search for pulsations and derive its spectral properties. We monitored the source with Swift to follow the evolution of the outburst and to look for further outbursts to verify the regular pattern seen in the optical light curve with a period of ~24.4 d. The XMM-Newton EPIC light curves show coherent X-ray pulsations with a period of 8.035331(15) s (1 sigma). The X-ray spectrum can be modelled by an absorbed power law with photon index of ~0.8, an additional black-body component with temperature of ~0.25 keV and an Fe K line. Phase-resolved X-ray spectroscopy reveals that the spectrum varies with pulse phase. We confirm the identification of the optical counterpart within the error circle of XMM-Newton at an angular distance of ~0.8 arcsec, which is an O9Ve star with known Halpha emission. By analyzing the combined data from three OGLE phases we derived an optical period of 24.43 d.The X-ray pulsations and long-term variability, as well as the properties of the optical counterpart, confirm that RX J0520.5-6932 is a Be/X-ray binary pulsar in the Large Magellanic Cloud. Based on the X-ray monitoring of the source we conclude that the event in January 2013 was a moderately bright type-I X-ray outburst, with a peak luminosity of 1.79e36 erg/s.

We observed the newly discovered X-ray source Swift J053041.9-665426 in the X-ray and optical regime to confirm its proposed nature as a high mass X-ray binary. We obtained XMM-Newton and Swift X-ray data, along with optical observations with the ESO Faint Object Spectrograph, to investigate the spectral and temporal characteristics of Swift J053041.9-665426. The XMM-Newton data show coherent X-ray pulsations with a period of 28.77521(10) s (1 sigma). The X-ray spectrum can be modelled by an absorbed power law with photon index within the range 0.76 to 0.87. The addition of a black body component increases the quality of the fit but also leads to strong dependences of the photon index, black-body temperature and absorption column density. We identified the only optical counterpart within the error circle of XMM-Newton at an angular distance of ~0.8 arcsec, which is 2MASS J05304215-6654303. We performed optical spectroscopy from which we classify the companion as a B0-1.5Ve star. The X-ray pulsations and long-term variability, as well as the properties of the optical counterpart, confirm that Swift J053041.9-665426 is a new Be/X-ray binary pulsar in the Large Magellanic Cloud.