Search

We present the analysis of publicly available NuSTAR, Suzaku and XMM-Newton observations of the symbiotic recurrent nova T CrB covering the 2006.77-2022.66 yr period. The X-ray spectra are analysed by adopting a model that includes a reflection component produced by the presence of a disk that mimics the accretion disk and the immediate surrounding medium. Our best-fit model requires this disk to have a radius of 1 AU, effective thickness of 0.1 AU, averaged column density 10$^{25}$ cm$^{-2}$ and orientation of 50$^{\circ}$ with respect to the line of sight. This disk is about a factor of two larger than recent estimations for the accretion disk and its presence contributes significantly via reflection to the total X-ray flux detected from T CrB, which naturally produces the emission of the 6.4 keV Fe line. Our analysis suggests that the temperature of the boundary layer evolved from 14.8 keV in the steady-state phase (before 2016), to 2.8 keV in the 2017.24 epoch, to finally stabilise to about $\sim$8 keV in the subsequent epochs. These variations in the plasma temperature of the boundary layer are attributed to the evolution of the mass accretion rate ($\dot{M}_\mathrm{acc}$), which is estimated to have an averaged value of $\dot{M}_\mathrm{acc}$=2.6$\times10^{-8}$ M$_\odot$ yr$^{-1}$ for the current active phase. The presence of emission lines in the XMM-Newton RGS spectrum of 2017.24 prevents from adopting a black body emission model to fit the soft X-ray range. Instead, we use plasma emission models that suggest the presence of adiabatically-shocked gas produced by gas velocities of 110-200 km s$^{-1}$, very likely tracing jet-like ejections similar to what is found in other symbiotic systems. The analysis of X-ray and optical data together show that T CrB has a similar evolution as black hole binaries, accreting neutron stars and AGN in the hardness-intensity diagram., Comment: 12 pages, 7 figures, 2 tables; accepted to MNRAS

This work presents a methodological approach to generate realistic $\gamma$-ray light curves of pulsars, resembling reasonably well the observational ones observed by the Fermi-Large Area Telescope instrument, fitting at the same time their high-energy spectra. The theoretical light curves are obtained from a spectral and geometrical model of the synchro-curvature emission. Despite our model relies on a few effective physical parameters, the synthetic light curves present the same main features observed in the observational $\gamma$-ray light curve zoo, such as the different shapes, variety in the number of peaks, and a diversity of peak widths. The morphological features of the light curves allows us to statistically compare the observed properties. In particular, we find that the proportion on the number of peaks found in our synthetic light curves is in agreement with the observational one provided by the third Fermi-LAT pulsar catalog. We also found that the detection probability due to beaming is much higher for orthogonal rotators (approaching 100%) than for small inclination angles (less than 20%).The small variation on the synthetic skymaps generated for different pulsars indicates that the geometry dominates over timing and spectral properties in shaping the gamma-ray light curves. This means that geometrical parameters like the inclination angle can be in principle constrained by gamma-ray data alone independently on the specific properties of a pulsar. At the same time, we find that $\gamma$-ray spectra seen by different observers can slightly differ, opening the door to constraining the viewing angle of a particular pulsar., Comment: 14 pages, 8 figures. Accepted for publication in Monthly Notices of the Royal Astronomical Society

The detection of gravitational waves from a binary neutron star merger by Advanced LIGO and Advanced Virgo (GW170817), along with the discovery of the electromagnetic counterparts of this gravitational wave event, ushered in a new era of multimessenger astronomy, providing the first direct evidence that BNS mergers are progenitors of short gamma-ray bursts (GRBs). Such events may also produce very-high-energy (VHE, > 100GeV) photons which have yet to be detected in coincidence with a gravitational wave signal. The Cherenkov Telescope Array (CTA) is a next-generation VHE observatory which aims to be indispensable in this search, with an unparalleled sensitivity and ability to slew anywhere on the sky within a few tens of seconds. New observing modes and follow-up strategies are being developed for CTA to rapidly cover localization areas of gravitational wave events that are typically larger than the CTA field of view. This work will evaluate and provide estimations on the expected number of of gravitational wave events that will be observable with CTA, considering both on- and off-axis emission. In addition, we will present and discuss the prospects of potential follow-up strategies with CTA., Comment: Presented at the 38th International Cosmic Ray Conference (ICRC 2023), 2023 (arXiv:2309.08219)

Magnetars are neutron stars with extremely strong magnetic fields, frequently powering high-energy activity in X-rays. Pulsed radio emission following some X-ray outbursts have been detected, albeit its physical origin is unclear. It has long been speculated that the origin of magnetars' radio signals is different from those from canonical pulsars, although convincing evidence is still lacking. Five months after magnetar SGR 1935+2154's X-ray outburst and its associated Fast Radio Burst (FRB) 20200428, a radio pulsar phase was discovered. Here we report the discovery of X-ray spectral hardening associated with the emergence of periodic radio pulsations from SGR 1935+2154 and a detailed analysis of the properties of the radio pulses. The complex radio pulse morphology, which contains both narrow-band emission and frequency drifts, has not been seen before in other magnetars, but is similar to those of repeating FRBs - even though the luminosities are many orders of magnitude different. The observations suggest that radio emission originates from the outer magnetosphere of the magnetar, and the surface heating due to the bombardment of inward-going particles from the radio emission region is responsible for the observed X-ray spectral hardening., Comment: 47 pages, 11 figures

We present 294 pulsars found in GeV data from the Large Area Telescope (LAT) on the Fermi Gamma-ray Space Telescope. Another 33 millisecond pulsars (MSPs) discovered in deep radio searches of LAT sources will likely reveal pulsations once phase-connected rotation ephemerides are achieved. A further dozen optical and/or X-ray binary systems co-located with LAT sources also likely harbor gamma-ray MSPs. This catalog thus reports roughly 340 gamma-ray pulsars and candidates, 10% of all known pulsars, compared to $\leq 11$ known before Fermi. Half of the gamma-ray pulsars are young. Of these, the half that are undetected in radio have a broader Galactic latitude distribution than the young radio-loud pulsars. The others are MSPs, with 6 undetected in radio. Overall, >235 are bright enough above 50 MeV to fit the pulse profile, the energy spectrum, or both. For the common two-peaked profiles, the gamma-ray peak closest to the magnetic pole crossing generally has a softer spectrum. The spectral energy distributions tend to narrow as the spindown power $\dot E$ decreases to its observed minimum near $10^{33}$ erg s$^{-1}$, approaching the shape for synchrotron radiation from monoenergetic electrons. We calculate gamma-ray luminosities when distances are available. Our all-sky gamma-ray sensitivity map is useful for population syntheses. The electronic catalog version provides gamma-ray pulsar ephemerides, properties and fit results to guide and be compared with modeling results., Comment: 142 pages. Accepted by the Astrophysical Journal Supplement

Observational evidence points to a red supergiant (RSG) progenitor for SN 2023ixf. The progenitor candidate has been detected in archival images at wavelengths (>0.6 micron) where RSGs typically emit profusely. This object is distinctly variable in the infrared (IR). We characterize the variability using pre-explosion mid-IR (3.6 and 4.5 micron) Spitzer and ground-based near-IR (JHKs) archival data jointly covering 19 yr. The IR light curves exhibit significant variability with RMS amplitudes in the range of 0.2-0.4 mag, increasing with decreasing wavelength. From a robust period analysis of the more densely sampled Spitzer data, we measure a period of 1091+/-71 days. We demonstrate using Gaussian Process modeling that this periodicity is also present in the near-IR light curves, thus indicating a common physical origin, which is likely pulsational instability. We use a period-luminosity relation for RSGs to derive a value of M_K=-11.58+/-0.31 mag. Assuming a late M spectral type, this corresponds to log(L/L_sun)=5.27+/-0.12 at T_eff=3200 K and to log(L/L_sun)=5.37+/-0.12 at T_eff=3500 K. This gives an independent estimate of the progenitor's luminosity, unaffected by uncertainties in extinction and distance. Assuming the progenitor candidate underwent enhanced dust-driven mass-loss during the time of these archival observations, and using an empirical period-luminosity-based mass-loss prescription, we obtain a mass-loss rate of around (2-4)x10^-4 M_sun/yr. Comparing the above luminosity with stellar evolution models, we infer an initial mass for the progenitor candidate of 20+/-4 M_sun, making this one of the most massive progenitors for a Type II SN detected to-date., Comment: 20 pages, 2 tables, accepted to ApJ

We propose an experimental scheme for performing sensitive, high-precision laser spectroscopy studies on fast exotic isotopes. By inducing a step-wise resonant ionization of the atoms travelling inside an electric field and subsequently detecting the ion and the corresponding electron, time- and position-sensitive measurements of the resulting particles can be performed. Using a Mixture Density Network (MDN), we can leverage this information to predict the initial energy of individual atoms and thus apply a Doppler correction of the observed transition frequencies on an event-by-event basis. We conduct numerical simulations of the proposed experimental scheme and show that kHz-level uncertainties can be achieved for ion beams produced at extreme temperatures ($> 10^8$ K), with energy spreads as large as $10$ keV and non-uniform velocity distributions. The ability to perform in-flight spectroscopy, directly on highly energetic beams, offers unique opportunities to studying short-lived isotopes with lifetimes in the millisecond range and below, produced in low quantities, in hot and highly contaminated environments, without the need for cooling techniques. Such species are of marked interest for nuclear structure, astrophysics, and new physics searches.

In this work, we introduce a quantitative methodology to define what is the main trunk and what are the significant branches of a minimum spanning tree (MST). We apply it to the pulsar tree, i.e. the MST of the pulsar population constructed upon a Euclidean distance over the pulsar's intrinsic properties. Our method makes use of the betweenness centrality estimator, as well as of non-parametric tests to establish the distinct character of the defined branches. Armed with these concepts, we study how the pulsar population has evolved throughout history, and analyze how to judge whether a new class of pulsars appears in new data, future surveys, or new incarnations of pulsar catalogs., Comment: An updated app of the web application we develop for visualizing the pulsar tree can be found at http://www.pulsartree.ice.csic.es

The Cherenkov Telescope Array will provide the deepest survey of the Galactic Plane performed at very-high-energy gamma-rays. Consequently, this survey will unavoidably face the challenge of source confusion, i.e., the non-unique attribution of signal to a source due to multiple overlapping sources. Among the known populations of Galactic gamma-ray sources and given their extension and number, pulsar wind nebulae (PWNe, and PWN TeV halos) will be the most affected. We aim to probe source confusion of TeV PWNe in forthcoming CTA data. For this purpose, we performed and analyzed simulations of artificially confused PWNe with CTA. As a basis for our simulations, we applied our study to TeV data collected from the H.E.S.S. Galactic Plane Survey for ten extended and two point-like firmly identified PWNe, probing various configurations of source confusion involving different projected separations, relative orientations, flux levels, and extensions among sources. Source confusion, defined here to appear when the sum of the Gaussian width of two sources is larger than the separation between their centroids, occurred in $\sim$30% of the simulations. For this sample and 0.5$\deg$ of average separation between sources, we found that CTA can likely resolve up to 60% of those confused sources above 500 GeV. Finally, we also considered simulations of isolated extended sources to see how well they could be matched to a library of morphological templates. The outcome of the simulations indicates a remarkable capability (more than 95% of the cases studied) to match a simulation with the correct input template in its proper orientation., Comment: In press in MNRAS

The Large High Altitude Air Shower Observatory has reported the detection of cosmic-ray sources in Milky Way that can accelerate particles up to PeV (= 10$^{15}$ eV) energies. These sources, so called ``PeVatrons'', are mostly unidentified. Several classes of sources, such as supernova remnants, pulsar wind nebula, or young stellar clusters can potentially be the counterparts of these PeVatrons. The aim of this work is to study a pulsar wind nebula interpretation of one of these PeVatrons, LHAASO J2226+6057, which has a relatively well covered multi-frequency spectrum. We have performed a leptonic, time-dependent modeling of the pulsar wind nebula (PWN) associated with PSR J2229+6114 considering a time-energy-dependent diffusion-loss equation. Injection, energy losses, as well as escape of particles were considered to balance the time-dependent lepton population. We have also included the dynamics of the PWN and the associated supernova remnant (SNR) and their interaction via the reverse shock to study the reverberation phase of the system. We have considered different values of braking index ($n$) and true age ($t_{age}$) for the fitting of the multi-wavelength (MWL) spectral energy distribution (SED) of LHAASO J2226+6057. The best-fit PWN model parameters and their 1$\sigma$ confidence intervals were evaluated. We have also demonstrated the impact of reverberation on the MWL SED with increasing time. Additionally, we have discussed the resultant large radius and low magnetic field associated with the PWN in question, as caveats for the possible physical connection of the pulsar as the origin of this high energy source., Comment: 10 pages, 5 figures, 1 table; accepted for publication in Astronomy & Astrophysics (A&A)

An incremental version of the fourth catalog of active galactic nuclei (AGNs) detected by the Fermi-Large Area Telescope is presented. This version (4LAC-DR3) derives from the third data release of the 4FGL catalog based on 12 years of E>50 MeV gamma-ray data, where the spectral parameters, spectral energy distributions (SEDs), yearly light curves, and associations have been updated for all sources. The new reported AGNs include 587 blazar candidates and four radio galaxies. We describe the properties of the new sample and outline changes affecting the previously published one. We also introduce two new parameters in this release, namely the peak energy of the SED high-energy component and the corresponding flux. These parameters allow an assessment of the Compton dominance, the ratio of the Inverse-Compton to the synchrotron peak luminosities, without relying on X-ray data., Comment: 13 pages, 7 figures, accepted for publication in The Astrophysical Journal Supplement Series. Fits files are available at https://fermi.gsfc.nasa.gov/ssc/data/access/lat/4LACDR3/ or alternatively at https://www.ssdc.asi.it/fermi4lac-DR3/table-4LAC-DR3-h.fits and https://www.ssdc.asi.it/fermi4lac-DR3/table-4LAC-DR3-l.fits

The detected high-energy pulsars' population is growing in number, and thus, having agile and physically relevant codes to analyze it consistently is important. Here, we update our existing synchro-curvature radiation model by including a better treatment of the particle injection, particularly where the large pitch angle particles dominate the spectra, and by implementing a fast and accurate minimization technique. The latter allows a large improvement in computational cost, needed to test model enhancements and to apply the model to a larger pulsar population. We successfully fit the sample of pulsars with X-ray and $\gamma$-ray data. Our results indicate that, for every emitting particle, the spatial extent of their trajectory where the pitch angle is large and most of the detected X-ray radiation is produced is a small fraction of the light cylinder. We also confirm with this new approach that synchrotron radiation is not negligible for most of the gamma-ray pulsars detected. In addition, with the results obtained, we argue that J0357+3205 and J2055+2539 are MeV-pulsar candidates and are suggested for exhaustive observations in this energy band., Comment: 12 pages, 7 figures. Accepted for publication in Monthly Notices of the Royal Astronomical Society

The $P\dot P$ diagram is a cornerstone of pulsar research. It is used in multiple ways for classifying the population, understanding evolutionary tracks, identifying issues in our theoretical reach, and more. However, we have been looking at the same plot for more than five decades. A fresh appraisal may be healthy. Is the $P\dot P$-diagram the most useful or complete way to visualize the pulsars we know? Here we pose a fresh look at the information we have on the pulsar population. First, we use principal components analysis over magnitudes depending on the intrinsic pulsar's timing properties (proxies to relevant physical pulsar features), to analyze whether the information contained by the pulsar's period and period derivative is enough to describe the variety of the pulsar population. Even when the variables of interest depend on $P$ and $\dot P$, we show that $P\dot P$ are not principal components. Thus, any distance ranking or visualization based only on $P$ and $\dot P$ is potentially misleading. Next, we define and compute a properly normalized distance to measure pulsar nearness, calculate the minimum spanning tree of the population, and discuss possible applications. The pulsar tree hosts information about pulsar similarities that go beyond $P$ and $\dot P$, and are thus naturally difficult to read from the $P\dot P$-diagram. We use this work to introduce the pulsar tree website http://www.pulsartree.ice.csic.es containing visualization tools and data to allow users to gather information in terms of MST and distance ranking., Comment: In press in MNRAS. The pulsar tree website is at http://www.pulsartree.ice.csic.es

The period and the period derivative of a pulsar are critical magnitudes for defining the properties of the magnetospheric size and plasma dynamics. The pulsar light cylinder, the magnetic field intensity nearby it, and the curvature radius all depend on these timing properties, and shape the observed high-energy synchro-curvature emission. Therefore, the radiative properties of pulsars are inextricably linked to them. This fact poses the question of how well does a given pulsar's spectral energy distribution embeds information of the timing parameters, and if so, whether we can deduce them if they have not been measured directly. This is relevant to possibly constrain the timing properties of potential pulsar candidates among unidentified $\gamma$-ray sources. We consider well-measured pulsar spectra blinding us from the knowledge of their timing properties, and address this question by using our radiative synchro-curvature model that was proven able to fit the observed spectra of the pulsar population. We find that in the majority of the cases studied (8/13), the spin period is constrained within a range of about one order of magnitude, within which the real period lies. In the other cases, there is degeneracy and no period range can be constrained. This can be used to facilitate the blind search of pulsed signals., Comment: 8 pages, 4 figures. Accepted for publication in The Astrophysical Journal

We present an incremental version (4FGL-DR3, for Data Release 3) of the fourth Fermi-LAT catalog of gamma-ray sources. Based on the first twelve years of science data in the energy range from 50 MeV to 1 TeV, it contains 6658 sources. The analysis improves on that used for the 4FGL catalog over eight years of data: more sources are fit with curved spectra, we introduce a more robust spectral parameterization for pulsars, and we extend the spectral points to 1 TeV. The spectral parameters, spectral energy distributions, and associations are updated for all sources. Light curves are rebuilt for all sources with 1 yr intervals (not 2 month intervals). Among the 5064 original 4FGL sources, 16 were deleted, 112 are formally below the detection threshold over 12 yr (but are kept in the list), while 74 are newly associated, 10 have an improved association, and seven associations were withdrawn. Pulsars are split explicitly between young and millisecond pulsars. Pulsars and binaries newly detected in LAT sources, as well as more than 100 newly classified blazars, are reported. We add three extended sources and 1607 new point sources, mostly just above the detection threshold, among which eight are considered identified, and 699 have a plausible counterpart at other wavelengths. We discuss degree-scale residuals to the global sky model and clusters of soft unassociated point sources close to the Galactic plane, which are possibly related to limitations of the interstellar emission model and missing extended sources., Comment: accepted in ApJS; follow-up paper to 1902.10045

Star-forming regions have been proposed as potential Galactic cosmic-ray accelerators for decades. Cosmic ray acceleration can be probed through observations of gamma-rays produced in inelastic proton-proton collisions, at GeV and TeV energies. We analyze more than 11 years of Fermi-LAT data from the direction of Westerlund 2, one of the most massive and best-studied star-forming regions in our Galaxy. The spectral and morphological characteristics of the LAT source agree with the ones in the TeV regime (HESS J1023-575), allowing the description of the gamma-ray source from a few hundreds of MeV to a few tens of TeVs. We will present the results and discuss the implications of the identification with the stellar cluster and the radiation mechanism involved., Comment: Proceeding for contribution to the 37th International Cosmic Ray Conference (ICRC2021)

We report on the analysis of a deep Chandra observation of the high-magnetic field pulsar (PSR) J1119-6127 and its compact pulsar wind nebula (PWN) taken in October 2019, three years after the source went into outburst. The 0.5-7 keV post-outburst (2019) spectrum of the pulsar is best described by a two-component blackbody plus powerlaw model with a temperature of 0.2\pm0.1 keV, photon index of 1.8\pm0.4 and X-ray luminosity of ~1.9e33 erg s^{-1}, consistent with its pre-burst quiescent phase. We find that the pulsar has gone back to quiescence. The compact nebula shows a jet-like morphology elongated in the north-south direction, similar to the pre-burst phase. The post-outburst PWN spectrum is best fit by an absorbed powerlaw with a photon index of 2.3\pm0.5 and flux of ~3.2e-14 erg cm^{-2} s^{-1} (0.5-7 keV). The PWN spectrum shows evidence of spectral softening in the post-outburst phase, with the pre-burst photon index of 1.2\pm0.4 changing to 2.3\pm0.5, and pre-burst luminosity of ~1.5e32 erg s^{-1} changing to 2.7e32 erg s^{-1} in the 0.5-7 keV band, suggesting magnetar outbursts can impact PWNe. The observed timescale for returning to quiescence, of just a few years, implies a rather fast cooling process and favors a scenario where J1119 is temporarily powered by magnetic energy following the magnetar outburst, in addition to its spin-down energy., Comment: 10 pages, 4 figures, 2 tables, Accepted for publication in the Astrophysical Journal

Star-forming regions have been proposed as potential Galactic cosmic-ray accelerators for decades. Cosmic-ray acceleration can be probed through observations of gamma-rays produced in inelastic proton-proton collisions, at GeV and TeV energies. In this paper, we analyze more than 11 years of Fermi-LAT data from the direction of Westerlund 2, one of the most massive and best-studied star-forming regions in our Galaxy. In particular, we investigate the characteristics of the bright pulsar PSR J1023-5746 that dominates the gamma-ray emission below a few GeV at the position of Westerlund 2, and the underlying extended source FGES J1023.3-5747. The analysis results in a clear identification of FGES J1023.3-5747 as the GeV counterpart of the TeV source HESS J1023-575, through its morphological and spectral properties. This identification provides new clues about the origin of the HESS J1023-575 gamma-ray emission, favouring a hadronic origin of the emission, powered by Westerlund 2, rather than a leptonic origin related to either the pulsar wind nebula associated with PSR J1023-5746 or the cluster itself. This result indirectly supports the hypothesis that star-forming regions can contribute to the cosmic-ray sea observed in our Galaxy, Comment: In press in MNRAS

The unidentified TeV source MGRO J1908+06, with emission extending from hundreds of GeV to beyond 100TeV, is one of the most intriguing sources in the Galactic plane. MGRO J1908+06 spatially associates with an IceCube hotspot of neutrino emission. Although the hotspot is not significant yet, this suggests a possible hadronic origin of the observed gamma-ray radiation. Here we describe a multiwavelength analysis on MGRO J1908+06 to determine its nature. We identify, for the first time, an extended GeV source as the counterpart of MGRO J1908+06, discovering possibly associated molecular clouds (MCs). The GeV spectrum shows two well-differentiated components: a soft spectral component below $\sim10$ GeV, and a hard one ($\Gamma\sim1.6$) above these energies. The lower-energy part is likely associated with the dense MCs surrounding the supernova remnant SNR G40.5$-$0.5, whereas the higher-energy component, which connects smoothly with the spectrum observed in TeV range, resembles the inverse Compton emission observed in relic pulsar wind nebulae. This simple scenario seems to describe the data satisfactorily, but raises questions about the interpretation of the emission at hundreds of TeV. In this scenario, no detectable neutrino flux would be expected., Comment: ApJL in press

Since 2009, several rapid and bright flares have been observed at high energies (>100 MeV) from the direction of the Crab Nebula. Several hypotheses have been put forward to explain this phenomenon, but the origin is still unclear. The detection of counterparts at higher energies with the next generation of Cherenkov telescopes will be determinant to constrain the underlying emission mechanisms. We aim at studying the capability of the Cherenkov Telescope Array (CTA) to explore the physics behind the flares, by performing simulations of the Crab Nebula spectral energy distribution, both in flaring and steady state, for different parameters related to the physical conditions in the nebula. In particular, we explore the data recorded by Fermi during two particular flares that occurred in 2011 and 2013. The expected GeV and TeV gamma-ray emission is derived using different radiation models. The resulting emission is convoluted with the CTA response and tested for detection, obtaining an exclusion region for the space of parameters that rule the different flare emission models. Our simulations show different scenarios that may be favourable for achieving the detection of the flares in Crab with CTA, in different regimes of energy. In particular, we find that observations with low sub-100 GeV energy threshold telescopes could provide the most model-constraining results., Comment: 11 pages, 9 figures; typos corrected

High mass X-ray binaries are among the brightest X-ray sources in the Milky Way, as well as in nearby Galaxies. Thanks to their highly variable emissions and complex phenomenology, they have attracted the interest of the high energy astrophysical community since the dawn of X-ray Astronomy. In more recent years, they have challenged our comprehension of physical processes in many more energy bands, ranging from the infrared to very high energies. In this review, we provide a broad but concise summary of the physical processes dominating the emission from high mass X-ray binaries across virtually the whole electromagnetic spectrum. These comprise the interaction of stellar winds with the high gravitational and magnetic fields of compact objects, the behaviour of matter under extreme magnetic and gravity conditions, and the perturbation of the massive star evolutionary processes by presence in a binary system. We highlight the role of the INTEGRAL mission in the discovery of many of the most interesting objects in the high mass X-ray binary class and its contribution in reviving the interest for these sources over the past two decades. We show how the INTEGRAL discoveries have not only contributed to significantly increase the number of high mass X-ray binaries known, thus advancing our understanding of the population as a whole, but also have opened new windows of investigation that stimulated the multi-wavelength approach nowadays common in most astrophysical research fields. We conclude the review by providing an overview of future facilities being planned from the X-ray to the very high energy domain that will hopefully help us in finding an answer to the many questions left open after more than 18 years of INTEGRAL scientific observations., Comment: 61 pages, 22 figures, 1 table, will be published in New Astronomy Reviews

Microquasars, the local siblings of extragalactic quasars, are binary systems comprising a compact object and a companion star. By accreting matter from their companions, microquasars launch powerful winds and jets, influencing the interstellar environment around them. Steady gamma-ray emission is expected to rise from their central objects, or from interactions between their outflows and the surrounding medium. The latter prediction was recently confirmed with the detection of SS 433 at high (TeV) energies. In this report, we analyze more than ten years of GeV gamma-ray data from the Fermi Gamma-ray Space Telescope on this source. Detailed scrutiny of the data reveal emission in the SS 433 vicinity, co-spatial with a gas enhancement, and hints for emission possibly associated with a terminal lobe of one of the jets. Both gamma-ray excesses are relatively far from the central binary, and the former shows evidence for a periodic variation at the precessional period of SS 433, linking it with the microquasar. This result challenges obvious interpretations and is unexpected from any previously published theoretical models. It provides us with a chance to unveil the particle transport from SS 433 and to probe the structure of the local magnetic field in its vicinity., Comment: published in Nature Astronomy

Insight-HXMT is the first Chinese X-ray astronomical mission, launched successfully on June 15, 2017, from China's Jiuquan Satellite Launch Center. Insight-HXMT was designed to have a broad energy coverage in X-rays, from 1-250 keV, with excellent timing and adequate energy resolution at soft X-rays, and the largest effective area at hard X-rays. Here, we present a collection of papers of the Journal of High Energy Astrophysics on the Early Results of China's 1st X-ray Astronomy Satellite Insight-HXMT. These papers cover the in-orbit performance, the background model, and all calibration results, together with several first results on observations of binaries and details on the Galactic plane survey., Comment: Published to introduce a set of 20+ papers from the HXMT collaboration. Most of the referred papers can be accessed in astro-ph already, and all will be avaialble from JHEAp's webpage too, at https://www.sciencedirect.com/journal/journal-of-high-energy-astrophysics/special-issue/10B0FLKPH6C

The Pierre Auger Collaboration has provided a compelling indication for a possible correlation between the arrival directions of ultrahigh-energy cosmic rays and nearby starburst galaxies. Herein we show how the latest large-scale modeling of starburst galaxies is compatible with the cosmic rays producing the anisotropy signal being accelerated at the terminal shock of superwinds., Comment: Matching version to be published in Physical Review D; 9 pages revtex, 5 .pdf figures

This chapter provides a phenomenological appraisal of the high-energy emission of millisecond pulsars. We comment on some of their properties as a population, as well as consider the especial cases of transitional pulsars, other redbacks, and black widow systems., Comment: Draft for a chapter of the book 'Millisecond Pulsars', of the Astrophysics and Space Science Library (ASSL) series edited by Sudip Bhattacharyya, Alessandro Papitto and Dipankar Bhattacharya

We report on a detailed spectral characterization of the non-thermal X-ray emission for a large sample of gamma-ray pulsars in the second Fermi-LAT catalogue. We outline the criteria adopted for the selection of our sample, its completeness, and critically describe different approaches to estimate the spectral shape and flux of pulsars. We perform a systematic modelling of the pulsars' X-ray spectra using archival observations with XMM-Newton, Chandra, and NuSTAR and extract the corresponding non-thermal X-ray spectral distributions. This set of data is made available online and is useful to confront with predictions of theoretical models., Comment: MNRAS in press

The Cherenkov Telescope Array (CTA) will be the next generation instrument for the very high energy gamma-ray astrophysics domain. With its enhanced sensitivity in comparison with the current facilities, CTA is expected to shed light on a varied population of sources. In particular, we will achieve a deeper knowledge of the Crab nebula and pulsar, which are the best characterized pulsar wind nebula and rotation powered pulsar, respectively. We aim at studying the capabilities of CTA regarding these objects through simulations, using the main tools currently in development for the CTA future data analysis: Gammapy and ctools. We conclude that, even using conservative Instrument Response Functions, CTA will be able to resolve many uncertainties regarding the spectrum and morphology of the pulsar and its nebula. The large energy range covered by CTA will allow us to disentangle the nebula spectral shape among different hypotheses, corresponding to different underlying emitting mechanisms. In addition, resolving internal structures (smaller than ~ 0.02 degrees in size) in the nebula and unveiling their location, would provide crucial information about the propagation of particles in the magnetized medium. We used a theoretical asymmetric model to characterise the morphology of the nebula and we showed that if predictions of such morphology exist, for instance as a result of hydrodynamical or magneto-hydrodynamical simulations, it can be directly compared with CTA results. We also tested the capability of CTA to detect periodic radiation from the Crab pulsar obtaining a precise measurement of different light curves shapes., Comment: in press in MNRAS

In this work, we introduce the use of the differential geometry Frenet-Serret equations to describe a magnetic line in a pulsar magnetosphere. These equations, which need to be solved numerically, fix the magnetic line in terms of their tangent, normal, and binormal vectors at each position, given assumptions on the radius of curvature and torsion. Once the representation of the magnetic line is defined, we provide the relevant set of transformations between reference frames; the ultimate aim is to express the map of the emission directions in the star co-rotating frame. In this frame, an emission map can be directly read as a light curve seen by observers located at a certain fixed angle with respect to the rotational axis. We provide a detailed step-by-step numerical recipe to obtain the emission map for a given emission process, and give a set of simplified benchmark tests. Key to our approach is that it offers a setting to achieve an effective description of the system's geometry {\it together} with the radiation spectrum. This allows to compute multi-frequency light curves produced by a specific radiation process (and not just geometry) in the pulsar magnetosphere, and intimately relates with averaged observables such as the spectral energy distribution., Comment: Accepted for publication in MNRAS

We apply a synchro-curvature spectral emission model based on characterizing the dynamics of magnetospheric particles to fit the phase-average spectra of the most extended database for the non-thermal spectra of pulsars. We consider 36 pulsars with well-determined non-thermal spectra from X-rays to gamma-rays. The sample includes Crab and Crab twin, for which the spectra extends even to the optical/ultraviolet and infrared energies. We find that the model --with just three physical parameters and a global scaling-- can fit the observations well across eight orders of magnitude for 18 of the 36 pulsars studied. Additionally, we find a set of 8 pulsars for which the model still provides arguably good fits and another set of 10 pulsars for which the model fails in reproducing the spectra. We discuss why, propose and provide physical interpretations for a simple model extension (related to the geometry of the accelerating system with regards to the observer) that allows dealing with all such cases, ultimately providing very good fits for all pulsars. The extended model is still austere, adding only two additional parameters to the former set, of the same kind of the ones previously used. We use these fits to discuss issues going from the observed spectral origin, to the extent of the dominance of synchrotron or curvature regimes, the use of a model as predictor for searching new non-thermal pulsars starting from gamma-ray surveys, and how the model offers a setting where phase shifts between X-ray and gamma-ray light curves would naturally arise., Comment: 20 pages, 15 figures, 2 tables. Accepted for publication in MNRAS on August 14

Auger data show evidence for a correlation between ultrahigh-energy cosmic rays (UHECRs) and nearby starburst galaxies. This intriguing correlation is consistent with data collected by the Telescope Array, which have revealed a much more pronounced directional ``hot spot'' in arrival directions not far from the starburst galaxy M82. In this work, we assume starbursts are sources of UHECRs and investigate the prospects to use the observed distribution of UHECR arrival directions to constrain Galactic magnetic field models. We show that if the Telescope Array hot spot indeed originates from M82, UHECR data would place a strong constraint on the coherent and turbulent components of the Galactic magnetic field., Comment: To appear in Proceedings of the 36th International Cosmic Ray Conference

In a recent work, we numerically studied the radiative properties of the reverberation phase of pulsar wind nebulae (PWNe), i.e., when the reverse shock created by the supernova explosion travels back towards the pulsar, compressing the wind bubble. We focused on several well-characterized PWNe and used them as examples for introducing the concept of superefficiency. The latter is a period of the PWN evolution, happening within reverberation, where the luminosity in a given band exceeds the spin-down power at the time. Here, we explore a broad range of PWN models to study their reverberation and superefficiency phases in a systematic way. Armed with these models we consider two aspects: On the one hand, we analyze via Monte Carlo simulations how many Galactic PWNe are expected to be reverberating or in a superefficiency stage at any given time, providing the first such estimations. On the other hand, we focus on searching for observational signatures of such periods. We analyze archival observations and check for the existence of possible candidates for superefficient PWNe. We also provide predictions for the future evolution of the magnetar nebula J1834.9-0846 (which we consider to be starting its reverberation period) along the next 50 years. Using our simulations as input we study how sensitive current and future X-ray satellites (like eXTP or Athena) will be to observe such evolution, concluding that they will be able to track it in detail., Comment: Accepted for publication in MNRAS

We report on a multi-wavelength study of the unclassified X-ray source CXOU J110926.4-650224 (J1109). We identified the optical counterpart as a blue star with a magnitude of $\sim$20.1 (3300-10500 $\require{mediawiki-texvc} \AA$). The optical emission was variable on timescales from hundreds to thousands of seconds. The spectrum showed prominent emission lines with variable profiles at different epochs. Simultaneous XMM-Newton and NuSTAR observations revealed a bimodal distribution of the X-ray count rates on timescales as short as tens of seconds, as well as sporadic flaring activity. The average broad-band (0.3-79 keV) spectrum was adequately described by an absorbed power law model with photon index of $\Gamma$=1.63$\pm$0.01 (at 1$\sigma$ c.l.), and the X-ray luminosity was (2.16$\pm$0.04)$\times$10$^{34}$ erg s$^{-1}$ for a distance of 4 kpc. Based on observations with different instruments, the X-ray luminosity has remained relatively steady over the past $\sim$15 years. J1109 is spatially associated with the gamma-ray source FL8Y J1109.8-6500, which was detected with Fermi at an average luminosity of (1.5$\pm$0.2)$\times$10$^{34}$ erg s$^{-1}$ (assuming the distance of J1109) over the 0.1-300 GeV energy band between 2008 and 2016. The source was undetected during ATCA radio observations that were simultaneous with NuSTAR, down to a 3$\sigma$ flux upper limit of 18 $\mu$Jy/beam (at 7.25 GHz). We show that the phenomenological properties of J1109 point to a binary transitional pulsar candidate currently in a sub-luminous accretion disk state, and that the upper limits derived for the radio emission are consistent with the expected radio luminosity for accreting neutron stars at similar X-ray luminosities., Comment: 21 pages, 15 figures, accepted for publication in A&A

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., Comment: Accepted for publication on Sci. China Phys. Mech. Astron. (2019)