The All-sky Medium Energy Gamma-ray Observatory (AMEGO) is a probe class mission concept that will provide essential contributions to multimessenger astrophysics in the late 2020s and beyond. AMEGO combines high sensitivity in the 200 keV to 10 GeV energy range with a wide field of view, good spectral resolution, and polarization sensitivity. Therefore, AMEGO is key in the study of multimessenger astrophysical objects that have unique signatures in the gamma-ray regime, such as neutron star mergers, supernovae, and flaring active galactic nuclei. The order-of-magnitude improvement compared to previous MeV missions also enables discoveries of a wide range of phenomena whose energy output peaks in the relatively unexplored medium-energy gamma-ray band.

Supergiant X-ray binaries usually comprise a neutron star accreting from the wind of a OB supergiant companion. They are classified as classical systems and the supergiant fast X-ray transients (SFXTs). The different behavior of these sub-classes of sources in X-rays, with SFXTs displaying much more pronounced variability, is usually (at least) partly ascribed to different physical properties of the massive star clumpy stellar wind. In case of SFXTs, a systematic investigation of the effects of clumps on flares/outbursts of these sources has been reported by Bozzo et al. (2017) exploiting the capabilities of the instruments on-board XMM-Newton to perform a hardness-resolved spectral analysis on timescales as short as a few hundreds of seconds. In this paper, we use six XMM-Newton observations of IGR J18027-2016 to extend the above study to a classical supergiant X-ray binary and compare the findings with those derived in the case of SFXTs. As these observations of IGR J18027-2016 span different orbital phases, we also study its X-ray spectral variability on longer timescales and compare our results with previous publications. Although obtaining measurements of the clump physical properties from X-ray observations of accreting supergiant X-ray binaries was already proven to be challenging, our study shows that similar imprints of clumps are found in the X-ray observations of the supergiant fast X-ray transients and at least one classical system, i.e. IGR J18027-2016. This provides interesting perspectives to further extend this study to many XMM-Newton observations already performed in the direction of other classical supergiant X-ray binaries.

Vela X-1 is an eclipsing high mass X-ray binary (HMXB) consisting of a 283s accreting X-ray pulsar in a close orbit of 8.964 days around the B0.5Ib supergiant HD77581 at a distance of just 2.4 kpc. The system is considered a prototype of wind-accreting HMXB and it has been used as a baseline in different theoretical or modelling studies. We discuss the observational properties of the system and the use of the observational data as laboratory to test recent developments in modelling the accretion process in High-Mass X-ray Binaries (e.g., Sander et al. 2018; El Mellah et al. 2018), which range from detailed descriptions of the wind acceleration to modelling of the structure of the flow of matter close to the neutron star and its variations.

Enabled by the Fermi Large Area Telescope, we now know young and recycled pulsars fill the gamma-ray sky, and we are beginning to understand their emission mechanism and their distribution throughout the Galaxy. However, key questions remain: Is there a large population of pulsars near the Galactic center? Why do the most energetic pulsars shine so brightly in MeV gamma rays but not always at GeV energies? What is the source and nature of the pair plasma in pulsar magnetospheres, and what role does the polar cap accelerator play? Addressing these questions calls for a sensitive, wide-field MeV telescope, which can detect the population of MeV-peaked pulsars hinted at by Fermi and hard X-ray telescopes and characterize their spectral shape and polarization.

In this paper we present the science potential of the enhanced X-ray Timing and Polarimetry (eXTP) mission for studies of strongly magnetized objects. We will focus on the physics and astrophysics of strongly magnetized objects, namely magnetars, accreting X-ray pulsars, and rotation powered pulsars. We also discuss the science potential of eXTP for QED studies. Developed by an international Consortium led by the Institute of High Energy Physics of the Chinese Academy of Sciences, the eXTP mission is expected to be launched in the mid 2020s.

In this White Paper we present the potential of the enhanced X-ray Timing and Polarimetry (eXTP) mission for studies related to Observatory Science targets. These include flaring stars, supernova remnants, accreting white dwarfs, low and high mass X-ray binaries, radio quiet and radio loud active galactic nuclei, tidal disruption events, and gamma-ray bursts. eXTP will be excellently suited to study one common aspect of these objects: their often transient nature. Developed by an international Consortium led by the Institute of High Energy Physics of the Chinese Academy of Science, the eXTP mission is expected to be launched in the mid 2020s.

In this paper we describe the potential of the enhanced X-ray Timing and Polarimetry (eXTP) mission for studies related to accretion flows in the strong field gravity regime around both stellar-mass and supermassive black-holes. eXTP has the unique capability of using advanced 'spectral-timing-polarimetry' techniques to analyze the rapid variations with three orthogonal diagnostics of the flow and its geometry, yielding unprecedented insight into the inner accreting regions, the effects of strong field gravity on the material within them and the powerful outflows which are driven by the accretion process.

We have modeled hydrodynamical tori in the Klu{\'z}niak-Lee pseudo-Newtonian potential. The tori in equilibrium were perturbed with uniform sub-sonic velocity fields: vertical, radial and diagonal respectively, and allowed to evolve in time. We identify the eigenmodes corresponding to those of slender tori. The results of our simulations are relevant to the investigation of high-frequency quasi-periodic oscillations observed in stellar-mass black hole binaries.

Strong winds from massive stars are a topic of interest to a wide range of astrophysical fields. In High-Mass X-ray Binaries the presence of an accreting compact object on the one side allows to infer wind parameters from studies of the varying properties of the emitted X-rays; but on the other side the accretor's gravity and ionizing radiation can strongly influence the wind flow. Based on a collaborative effort of astronomers both from the stellar wind and the X-ray community, this presentation attempts to review our current state of knowledge and indicate avenues for future progress.

In this paper we report on a long multi-wavelength observational campaign of the supergiant fast X-ray transient prototype IGR J17544-2619. A 150 ks-long observation was carried out simultaneously with XMM-Newton and NuSTAR, catching the source in an initial faint X-ray state and then undergoing a bright X-ray outburst lasting about 7 ks. We studied the spectral variability during outburst and quiescence by using a thermal and bulk Comptonization model that is typically adopted to describe the X-ray spectral energy distribution of young pulsars in high mass X-ray binaries. Although the statistics of the collected X-ray data were relatively high we could neither confirm the presence of a cyclotron line in the broad-band spectrum of the source (0.5-40 keV), nor detect any of the previously reported tentative detection of the source spin period. The monitoring carried out with Swift/XRT during the same orbit of the system observed by XMM-Newton and NuSTAR revealed that the source remained in a low emission state for most of the time, in agreement with the known property of all supergiant fast X-ray transients being significantly sub-luminous compared to other supergiant X-ray binaries. Optical and infrared observations were carried out for a total of a few thousands of seconds during the quiescence state of the source detected by XMM-Newton and NuSTAR. The measured optical and infrared magnitudes were slightly lower than previous values reported in the literature, but compatible with the known micro-variability of supergiant stars. UV observations obtained with the UVOT telescope on-board Swift did not reveal significant changes in the magnitude of the source in this energy domain compared to previously reported values.

T CrB is a symbiotic recurrent nova known to exhibit active phases, characterised by apparent increases in the hot component temperature and the appearance of flickering, i.e. changes in the observed flux on the time-scale of minutes. Historical UV observations have ruled out orbital variability as an explanation for flickering and instead suggest flickering is caused by variable mass transfer. We have analysed optical and X-ray observations to investigate the nature of the flickering as well as the active phases in T CrB. The spectroscopic and photometric observations confirm that the active phases follow two periods of ~1000d and ~5000d. Flickering in the X-rays is detected and follows an amplitude-flux relationship similar to that observed in the optical. The flickering is most prominent at harder X-ray energies, suggesting that it originates in the boundary layer between the accretion disc and the white dwarf. The X-ray radiation from the boundary layer is then reprocessed by a thick accretion disc or a nebula into UV radiation. A more detailed understanding of flickering would benefit from long-term simultaneous X-ray and optical monitoring of the phenomena in symbiotic recurrent novae and related systems such as Z And type symbiotic stars.

We performed axisymmetric hydrodynamical simulations of oscillating tori orbiting a non-rotating black hole. The tori in equilibrium were constructed with a constant distribution of angular momentum in a pseudo-Newtonian potential (Klu{\'z}niak-Lee). Motions of the torus were triggered by adding sub-sonic velocity fields: radial, vertical and diagonal to the tori in equilibrium. As the perturbed tori evolved in time, we measured $L_{2}$ norm of density and obtained the power spectrum of $L_{2}$ norm which manifested eigenfrequencies of tori modes. The most prominent modes of oscillation excited in the torus by a quasi-random perturbation are the breathing mode and the radial and vertical epicyclic modes. The radial and the plus modes, as well as the vertical and the breathing modes will have frequencies in an approximate 3:2 ratio if the torus is several Schwarzschild radii away from the innermost stable circular orbit. Results of our simulations may be of interest in the context of high-frequency quasi-periodic oscillations (HF QPOs) observed in stellar-mass black hole binaries, as well as in supermassive black holes.

We simulate an oscillating purely hydrodynamical torus with constant specific angular mo- mentum around a Schwarzschild black hole. The goal is to search for quasi-periodic oscil- lations (QPOs) in the light curve of the torus. The initial torus setup is subjected to radial, vertical and diagonal (combination of radial and vertical) velocity perturbations. The hydro- dynamical simulations are performed using the general relativistic magnetohydrodynamics code Cosmos++ and ray-traced using the GYOTO code. We found that a horizontal velocity perturbation triggers the radial and plus modes, while a vertical velocity perturbation trig- gers the vertical and X modes. The diagonal perturbation gives a combination of the modes triggered in the radial and vertical perturbations.

The early acceleration of stellar winds in massive stars is poorly constrained. The scattering of hard X-ray photons emitted by the pulsar in the high-mass X-ray binary Vela X-1 can be used to probe the stellar wind velocity and density profile close to the surface of its supergiant companion HD 77581. We built a high signal-to-noise and high resolution hard X-ray lightcurve of Vela X-1 measured by Swift/BAT over 300 orbital periods of the system and compared it with the predictions of a grid of hydrodynamic simulations. We obtain a very good agreement between observations and simulations for a narrow set of parameters, implying that the wind velocity close to the stellar surface is twice larger than usually assumed with the standard beta law. Locally a velocity gradient of $\beta\sim0.5$ is favoured. Even if still incomplete, hydrodynamic simulations are successfully reproducing several observational properties of Vela X-1.

This is a White Paper in support of the mission concept of the Large Observatory for X-ray Timing (LOFT), proposed as a medium-sized ESA mission. We discuss the potential of LOFT for the study of high-mass X-ray binaries and ultra-luminous X-ray sources. For a summary, we refer to the paper.

This is a White Paper in support of the mission concept of the Large Observatory for X-ray Timing (LOFT), proposed as a medium-sized ESA mission. We discuss the potential of LOFT for the study of active galactic nuclei. For a summary, we refer to the paper.

$Context$: Vela X-1 is the prototype of the classical sgHMXB systems. Recent continuous and long monitoring campaigns revealed a large hard X-rays variability amplitude with strong flares and off-states. This activity has been interpreted invoking clumpy stellar winds and/or magnetic gating mechanisms. $Aims$: We are probing if the observed behaviour could be explained by unstable hydrodynamic flows close to the neutron star rather than the more exotic phenomena. $Methods$: We have used the hydrodynamic code VH-1 to simulate the flow of the stellar wind with high temporal resolution and to compare the predicted accretion rate with the observed light-curves. $Results$: The simulation results are similar to the observed variability. Off-states are predicted with a duration of 5 to 120 minutes corresponding to transient low density bubbles forming around the neutron star. Oscillations of the accretion rate with a typical period of $\sim$ 6800 sec are generated in our simulations and observed. They correspond to the complex motion of a bow shock, moving either towards or away from the neutron star. Flares are also produced by the simulations up to a level of $10^{37}$ erg/s. $Conclusions$: We have qualitatively reproduced the hard X-ray variations observed in Vela X-1 with hydrodynamic instabilities predicted by a simple model. More sophisticated phenomena, such as clumpy winds or the magnetic gating mechanism, are not excluded but not required to explain the basic phenomenology.

The Large Observatory For x-ray Timing (LOFT) was studied within ESA M3 Cosmic Vision framework and participated in the final down-selection for a launch slot in 2022-2024. Thanks to the unprecedented combination of effective area and spectral resolution of its main instrument, LOFT will study the behaviour of matter under extreme conditions, such as the strong gravitational field in the innermost regions of accretion flows close to black holes and neutron stars, and the supra-nuclear densities in the interior of neutron stars. The science payload is based on a Large Area Detector (LAD, 10 m 2 effective area, 2-30 keV, 240 eV spectral resolution, 1 deg collimated field of view) and a WideField Monitor (WFM, 2-50 keV, 4 steradian field of view, 1 arcmin source location accuracy, 300 eV spectral resolution). The WFM is equipped with an on-board system for bright events (e.g. GRB) localization. The trigger time and position of these events are broadcast to the ground within 30 s from discovery. In this paper we present the status of the mission at the end of its Phase A study.

Supergiant High Mass X-ray Binary systems (sgHMXBs) consist of a massive, late type, star and a neutron star. The massive stars exhibits strong, radiatively driven, stellar winds. Wind accretion onto compact object triggers X-ray emission, which alters the stellar wind significantly. Hydrodynamic simulation has been used to study the neutron star - stellar wind interaction it two sgHMXBs: i) A heavily obscured sgHMXB (IGR J17252-3616) discovered by INTEGRAL. To account for observable quantities (i.e., absorbing column density) we have to assume a very slow wind terminal velocity of about 500 km/s and a rather massive neutron star. If confirmed in other obscured systems, this could provide a completely new stellar wind diagnostics. ii) A classical sgHMXB (Vela X-1) has been studied in depth to understand the origin of the off-states observed in this system. Among many models used to account for this observed behavior (clumpy wind, gating mechanism) we propose that self-organized criticality of the accretion stream is the likely reason for the observed behavior. In conclusion, the neutron star, in these two examples, acts very efficiently as a probe to study stellar winds.

A population of obscured supergiant High Mass X-ray Binaries (sgHMXBs) has been discovered by INTEGRAL. X-ray wind tomography of IGR J17252-3616 inferred a slow wind velocity to account for the enhanced obscuration. The main goal of this study is to understand under which conditions high obscuration could occur. We have used an hydrodynamical code to simulate the flow of the stellar wind around the neutron star. A grid of simulations was used to study the dependency of the absorbing column density and of the X-ray light-curves on the model parameters. A comparison between the simulation results and the observations of IGR J17252-3616 provides an estimate on these parameters. We have constrained the wind terminal velocity to 500-600 km/s and the neutron star mass to 1.75-2.15 solar masses. We have confirmed that the initial hypothesis of a slow wind velocity with a moderate mass loss rate is valid. The mass of the neutron star can be constrained by studying its impact on the accretion flow.

The LOFT mission concept is one of four candidates selected by ESA for the M3 launch opportunity as Medium Size missions of the Cosmic Vision programme. The launch window is currently planned for between 2022 and 2024. LOFT is designed to exploit the diagnostics of rapid X-ray flux and spectral variability that directly probe the motion of matter down to distances very close to black holes and neutron stars, as well as the physical state of ultra-dense matter. These primary science goals will be addressed by a payload composed of a Large Area Detector (LAD) and a Wide Field Monitor (WFM). The LAD is a collimated (<1 degree field of view) experiment operating in the energy range 2-50 keV, with a 10 m^2 peak effective area and an energy resolution of 260 eV at 6 keV. The WFM will operate in the same energy range as the LAD, enabling simultaneous monitoring of a few-steradian wide field of view, with an angular resolution of <5 arcmin. The LAD and WFM experiments will allow us to investigate variability from submillisecond QPO's to year-long transient outbursts. In this paper we report the current status of the project.

We report on all the INTEGRAL and Swift data collected during the first outburst observed from IGRJ18179-1621. The broad-band spectral analysis showed that the X-ray emission from the source is heavily absorbed (N_H~10^23 cm^-2), and well described by a flat power-law with a high energy rollover (cutoff energy 9-12 keV, e-folding energy 4-7 keV). We found some evidence of a cyclotron absorption feature at 22\pm1 keV. Together with the pulsations at 11.8s discovered in the XRT data, this evidence would suggest that IGRJ18179-1621 is an obscured magnetized accreting neutron star, possibly part of a supergiant high mass X-ray binary or a Be X-ray binary system.

INTEGRAL played a key role in discovering obscured sgHMXB in the Galaxy. We used XMM-Newton to perform X-ray wind tomography of a specific of these systems, IGR J17252-3616, featuring eclipses of the accreting pulsar. The X-ray band (0.2-10 keV) reveals vital information on the geometry of the surrounding gas probing simultaneously the absorption and the fluorescence emission. The XMM observations were scheduled to cover as many orbital phases as possible. Timing analysis allows the derivation of an accurate orbital solution and of the system parameters. Spectral analysis revealed remarkable variations of the absorbing column density along the orbit and of the Fe K$\alpha$ fluorescence line around the eclipse. The combination of these observables revealed a highly asymmetric and unprecedentedly extended structure in the stellar wind extending up to 2-3 stellar radii. The observations can be modeled in terms of three independent components: i) the unperturbed stellar wind ii) the contribution of a highly asymmetric hydrodynamic wind tail-like structure and iii) a cusp of material close to the neutron star. These dynamical structures are imaged for the first time in a sgHMXB and explain the source of the high obscuration.

AX J1910.7+0917 is a still unidentified source discovered with ASCA and observed more recently with IBIS/ISGRI, mainly noticeable for its rather hard spectrum. We analyzed all the public available data on this source, and we took advantage of the recent improvements performed in the INTEGRAL data analysis software to fully exploit the IBIS/ISGRI data. In the data collected from INTEGRAL, XMM-Newton, Chandra and ASCA the source is clearly variable. The spectrum can be modelled as an absorbed powerlaw (NH~6x10^22 cm^-2, Gamma=1.5) with an iron line at 6.4 keV. The present data still do not allow for a unique classification of the source. In the IBIS/ISGRI field of view around AX J1910.7+0917, we discovered three new sources: IGR J19173+0747, IGR J19294+1327 and IGR J19149+1036, where the latter is positionally coincident with the Einstein source 2E 1912.5+1031. For the first two sources we report results obtained from follow-up observations carried out with Swift/XRT.

We take advantage of the high sensitivity of the IBIS/ISGRI telescope and the improvements in the data analysis software to investigate the nature of the still poorly known X-ray source AX J1910.7+0917, and search for close-by previously undetected objects. We analyze all publicly available INTEGRAL data of AX J1910.7+0917, together with a number of archival observations that were carried out in the direction of the source with Chandra, XMM-Newton, and ASCA. In the IBIS/ISGRI field-of-view around AX J1910.7+0917, we discovered three new sources: IGR J19173+0747, IGR J19294+1327 and IGR J19149+1036; the latter is positionally coincident with the Einstein source 2E 1912.5+1031. For the first two sources, we also report the results of follow-up observations carried out with Swift/XRT. AX J1910.7+0917 features a clear variability in the X-rays. Its spectrum can be well described with an absorbed (N_H~6x10^(22) cm^(-2)) power-law ({\Gamma}~1.5) model plus an iron line at ~6.4 keV. We also obtained a refined position and report on possible infrared counterparts. The present data do not allow for a unique classification of the sources. Based on the property of its X-ray emission and the analysis of a likely infrared counterpart, we investigate different possibilities for the nature of AX J1910.7+0917.

Our goal is to understand the specificities of highly absorbed sgHMXB and in particular of the companion stellar wind, thought to be responsible for the strong absorption. We have monitored IGR J17252-3616, a highly absorbed system featuring eclipses, with XMM-Newton to study the vari- ability of the column density and of the Fe K{\alpha} emission line along the orbit and during the eclipses. We also built a 3D model of the structure of the stellar wind to reproduce the observed variability. We first derived a refined orbital solution built from INTEGRAL, RXTE and XMM data. The XMM monitoring campaign revealed significant variation of intrinsic absorbing column density along the orbit and of the Fe K{\alpha} line equivalent width around the eclipses. The origin of the soft X-ray absorption is modeled with an dense and extended hydrodynamical tail, trailing the neutron star. This structure extends along most of the orbit, indicating that the stellar wind is strongly disrupted by the neutron star. The variability of the absorbing column density suggests that the terminal velocity of the wind is smaller (~400 km/s) than observed in classical systems. This can also explain the much stronger density perturbation inferred from the observations. Most of the Fe K{\alpha} emission is generated in the most inner region of the hydrodynamical tail. This region, that extends over a few accretion radii, is ionized and does not contribute to the soft X-ray absorption. We have built a qualitative model of the stellar wind of IGR J17252-3616 that can represent the observations and suggest that highly absorbed systems have a lower wind velocity than classical sgHMXB. This proposal could be tested with de- tailed numerical simulations and high-resolution infrared/optical observations. If confirmed, it may turn out that half of the persistent sgHMXB have low stellar wind speeds.

Aims. Soft X-ray excesses have been detected in several Be/X-ray binaries and interpreted as the signature of hard X-ray reprocessing in the inner accretion disk. The system XMMU J054134.7-682550, located in the LMC, featured a giant Type II outburst in August 2007. The geometry of this system can be understood by studying the response of the soft excess emission to the hard X-ray pulses. Methods. We have analyzed series of simultaneous observations obtained with XMM-Newton/EPIC-MOS and RXTE/PCA in order to derive spectral and temporal characteristics of the system, before, during and after the giant outburst. Spectral fits were performed and a timing analysis has been carried out. Spectral variability, spin period evolution and energy dependent pulse shapes are analysed. Results. The outburst (L_X = 3* 10^38 erg/s \approx L_EDD) spectrum could be modeled successfully using a cutoff powerlaw, a cold disk emission, a hot blackbody, and a cyclotron absorption line. The magnetic field and magnetospheric radius could be constrained. The thickness of the inner accretion disk is broadened to a width of 75 km. The hot blackbody component features sinusoidal modulations indicating that the bulk of the hard X-ray emission is emitted preferentially along the magnetic equator. The spin period of the pulsar decreased very significantly during the outburst. This is consistent with a variety of neutron star equations of state and indicates a very high accretion rate.

Classical novae are quite frequent in M~31. However, very few spectra of M31 novae have been studied to date, especially during the early decline phase. Our aim is to study the photometric and spectral evolution of a M31 nova event close to outburst. We present photometric and spectroscopic observations of M31N 2005-09c, a classical nova in the disk of M31, using the 1.3m telescope of the Skinakas Observatory in Crete (Greece), starting on the 28th September, i.e. about 5 days after outburst, and ending on the 5th October 2005, i.e. about 12 days after outburst. We also have supplementary photometric observations from the La Sagra Observatory in Northern Andalucia, Spain, on September 29 and 30, October 3, 6 and 9 and November 1, 2005. The wavelength range covered by the spectra is from 3565 A to 8365 A. The spectra are of high S/N allowing the study of the evolution of the equivalent widths of the Balmer lines, as well as the identification of non-Balmer lines. The nova displays a typical early decline spectrum that is characterized by many weak FeII multiplet emissions. It is classified as a P$_{fe}$ nova. From the nova light curve, we have also derived its speed class, t=14+-2.5 days. As the nova evolved the Balmer lines became stronger and narrower. The early decline of the expansion velocity of the nova follows a power law in time with an exponent of \~-0.2.

We present the results of our search for optical counterparts to high-mass X-ray transient sources discovered by various X-ray missions. We obtained CCD images of the X-ray fields through BVR and Halpha filters in order to identify early-type stars in the R-Halpha versus B-V colour-colour diagram. We also obtained medium-resolution spectroscopy of the candidates in order to confirm the presence of H\alpha emission and perform spectral classification. We report on the discovery of the optical counterparts to two X-ray sources: XTE J1858+034 and IGR J01363+6610, and the follow-up observations of another two, newly identified by our group: SAX J2103.5+4545 and GRO J2058+42. For another source, IGR J00370+6122, we present the first detailed optical spectral analysis. The optical photometry and spectroscopy reveal B-type companions in all five sources; GRO J2058+42, SAX J2103.5+4545 and IGR J01363+6610 are positively identified with Be/X-ray binaries, IGR J00370+6122 with a supergiant X-ray binary, while the nature of XTE J1858+034 is uncertain. We also study the relationship between the optical and X-ray emission during quiescent states.