Georgios Vasilopoulos
Astrophysicist at National and Kapodistrian University of Athens
Latest Updates
Broadband study of the Be X-ray binary RX J0520.5-6932 during its outburst in 2024
Dec 2024
A new giant outburst of the Be X-ray binary RX J0520.5-6932 was detected and subsequently observed with several space-borne and ground-based instruments. This study presents a comprehensive analysis of the optical and X-ray data, focusing on the spectral and timing characteristics of selected X-ray observations. A joint fit of spectra from simultaneous observations performed by the X-ray telescope (XRT) on the Neil Gehrels Swift Observatory (Swift) and Nuclear Spectroscopic Telescope ARray (NuSTAR) provides broadband parameter constraints, including a cyclotron resonant scattering feature (CRSF) at 32.2(+0.8/-0.7) keV with no significant energy change since 2014, and a weaker Fe line. Independent spectral analyses of observations by the Lobster Eye Imager for Astronomy (LEIA), Einstein Probe (EP), Swift-XRT, and NuSTAR demonstrate the consistency of parameters across different bands. Luminosity variations during the current outburst were tracked. The light curve of the Optical Gravitational Lensing Experiment (OGLE) aligns with the X-ray data in both 2014 and 2024. Spin evolution over 10 years is studied after adding Fermi Gamma-ray Burst Monitor (GBM) data, improving the orbital parameters, with an estimated orbital period of 24.39 days, slightly differing from OGLE data. Despite intrinsic spin-up during outbursts, a spin-down of ~0.04s over 10.3 years is suggested. For the new outburst, the pulse profiles indicate a complicated energy-dependent shape, with decreases around 15 keV and 25 keV in the pulsed fraction, a first for an extragalactic source. Phase-resolved NuSTAR data indicate variations in parameters such as flux, photon index, and CRSF energy with rotation phase.
[read the article; Yang et al.]
Spectral Studies of Super-Eddington Accreting Neutron Stars in the Magellanic Clouds
November 2024
Major outbursts of BeXRBs offer a unique laboratory for studying accretion onto magnetized neutron stars (NSs) over a large dynamic range. The accreting material is entrained from the accretion disk by the strong magnetic field, and then channeled onto the NS, forming a so-called accretion column (AC). Physical simulation of the AC requires consideration of various physical processes occurring in strongly magnetized plasma, including complex multi-dimensional radiative transfer and the presence of a radiation-dominated shock. Analytical models based on the Becker and Wolff (2007) model have proven successful at reproducing the observed AC spectra in super-critical sources in which radiation pressure plays the dominant role in controlling the dynamics of the accreting material. In this study, we will apply the model to obtain spectra observed during super-Eddington outbursts of BeXRBs in the Magellanic Clouds.
[read the article; Vasilopoulos et al.]
Searching for Hadronic Signatures in the Time Domain of Blazar Emission: The Case of Mrk 501
September 2024
Blazars—a subclass of active galaxies—are intrinsically time-variable broadband sources of electromagnetic radiation. In this contribution, we explored relativistic proton (hadronic) signatures in the time domain blazar emission and searched for those parameter combinations that unveil their presence during flaring epochs. We generated time series for key model parameters, like magnetic field strength and the power-law index of radiating particles, which were motivated from a simulated time series with statistical properties describing the observed GeV gamma-ray flux. We chose the TeV blazar Mrk 501 as our test case, as it had been the study ground for extensive investigations during individual flaring events. Using the code LeHaMoC, we computed the electromagnetic and neutrino emissions for a period of several years that contained several flares of interest. We show that for both of those particle distributions the power-law index variations that were tied to moderate changes in the magnetic field strength of the emitting region might naturally lead to hard X-ray flares with very-high-energy γ-ray counterparts. We found spectral differences measurable by the Cherenkov Telescope Array Observatory at sub-TeV energies, and we computed the neutrino fluence over 14.5 years. The latter predicted ∼0.2 muon and anti-muon neutrinos, consistent with the non-detection of high-energy neutrinos from Mrk 501.
[read the article; Chatzis et al.]
TeV pion bumps in the gamma-ray spectra of flaring blazars
April 2024
Very high-energy (VHE, $E>100$ GeV) observations of the blazar Mrk 501 with MAGIC in 2014 provided evidence for an unusual narrow spectral feature at about 3 TeV during an extreme X-ray flaring activity. The one-zone synchrotron-self Compton scenario, widely used in blazar broadband spectral modeling, fails to explain the narrow TeV component. Motivated by this rare observation, we propose an alternative model for the production of narrow features in the VHE spectra of flaring blazars. These spectral features may result from the decay of neutral pions ($π^0$ bumps) that are in turn produced via interactions of protons (of tens of TeV energy) with energetic photons, whose density increases during hard X-ray flares. We explored the conditions needed for the emergence of narrow $π^0$ bumps in VHE blazar spectra during X-ray flares reaching synchrotron energies $\sim 100$ keV using time-dependent radiative transfer calculations. We focused on high-synchrotron peaked (HSP) blazars, which comprise the majority of VHE-detected extragalactic sources. We find that synchrotron-dominated flares with peak energies $\gtrsim100$ keV can be ideal periods for the search of $π^0$ bumps in the VHE spectra of HSP blazars. The flaring region is optically thin to photopion production, its energy content is dominated by the relativistic proton population, and the inferred jet power is highly super-Eddington. Application of the model to the spectral energy distribution of Mrk 501 on MJD 56857.98 shows that the VHE spectrum of the flare is described well by the sum of a synchrotron self-Compton (SSC) component and a distinct $π^0$ bump centered at 3 TeV. Spectral fitting of simulated SSC+$π^0$ spectra for the Cherenkov Telescope Array (CTA) show that a $π^0$ bump could be detected at a 5$σ$ significance level with a 30-min exposure.
[read the article; Petropoulou et al.]
eRASSU J060839.5-704014: A double degenerate ultra-compact binary in the direction of the LMC
March 2024
Context: During four all-sky surveys (eRASS1-4), eROSITA, the soft X-ray instrument aboard Spektrum-Roentgen-Gamma (SRG) detected a new supersoft X-ray source, eRASSU J060839.5−704014, in the direction of the Large Magellanic Cloud (LMC).
Aims: We arranged follow-up observations in the X-ray and optical wavelengths and further searched in archival observations to reveal the nature of the object.
Methods: Using X-ray observations with XMM-Newton we investigated the temporal and spectral behaviour of the source.
Results: We discover pulsations at 374 s with a pulse profile consistent with 100% modulation. We identify two other periodicities in the eROSITA data, which we establish as aliases due to the sampling of the eROSITA light curve. We identify a multi-wavelength counterpart to the X-ray source in UVW1 and g, r, i, and z images obtained by the optical/UV monitor on XMM-Newton and the Dark Energy Camera at the Cerro Tololo Inter-American Observatory. The timing and spectral characteristics of the source are consistent with a double degenerate ultra-compact binary system in the foreground of the LMC. eRASSU J060839.5−704014 belongs to a rare class of AM CVns, which are important to study in the context of progenitors of SN Ia and for persistent gravitational wave detection.
Conclusions: We identify eRASSU J060839.5−704014 as a new double degenerate ultra-compact binary located in the foreground of the LMC.
[read the article; Maitra et al.]
Application of neural networks to synchro-Compton blazar emission models
March 2024
Context: Jets from supermassive black holes at the centers of active galaxies are the most powerful and persistent sources of electromagnetic radiation in the Universe. To infer the physical conditions in the otherwise out-of-reach regions of extragalactic jets, we usually rely on fitting their spectral energy distributions (SEDs). The calculation of radiative models for the jet's non-thermal emission usually relies on numerical solvers of coupled partial differential equations.
Aims: In this work, we use machine learning to tackle the problem of high computational complexity to significantly reduce the SED model evaluation time, which is necessary for SED fittings carried out with Bayesian inference methods.
Methods: We computed the SEDs based on the synchrotron self-Compton model for blazar emission using the radiation code ATHEvA. We used them to train neural networks (NNs) to explore whether they can replace the original code, which is computationally expensive.
Results: We find that a NN with gated recurrent unit neurons (GRUN) can effectively replace the ATHEvA leptonic code for this application, while it can be efficiently coupled with Markov chain Monte Carlo (MCMC) and nested sampling algorithms for fitting purposes. We demonstrate this approach through an application to simulated data sets, as well as a subsequent application to observational data.
Conclusions: We present a proof-of-concept application of NNs to blazar science as the first step in a list of future applications involving hadronic processes and even larger parameter spaces. We offer this tool to the community through a public repository.
The results of our work are available in GitHub; Blazar_ML. This includes: (a) the NN and accompanied code produced to train them, (b) code for visualization of results in python and jupyter notebooks with instructions, and (c) part of the ATHEvA datasets that can be used for evaluation and plotting examples.
[read the article; Tzavelas et al.]
LeHaMoC: A versatile time-dependent lepto-hadronic modeling code for high-energy astrophysical sources
March 2024
Context. Recent associations of high-energy neutrinos with active galactic nuclei (AGN) have revived the interest in leptohadronic models of radiation from astrophysical sources. The rapid increase in the amount of acquired multi-messenger data will require fast numerical models that may be applied to large source samples.
Aims: We develop a time-dependent leptohadronic code, LeHaMoC, that offers several notable benefits compared to other existing codes, such as versatility and speed.
Methods: LeHaMoC solves the Fokker-Planck equations of photons and relativistic particles (i.e. electrons, positrons, protons, and neutrinos) produced in a homogeneous magnetized source that may also be expanding. The code utilizes a fully implicit difference scheme that allows fast computation of steady-state and dynamically evolving physical problems.
Results: We first present test cases where we compare the numerical results obtained with LeHaMoC against exact analytical solutions and numerical results computed with ATHEvA, a well-tested code of similar philosophy but a different numerical implementation. We find a good agreement (within 10-30%) with the numerical results obtained with ATHEvA without evidence of systematic differences. We then demonstrate the capabilities of the code through illustrative examples. First, we fit the spectral energy distribution from a jetted AGN in the context of a synchrotron-self Compton model and a proton-synchrotron model using Bayesian inference. Second, we compute the high-energy neutrino signal and the electromagnetic cascade induced by hadronic interactions in the corona of NGC 1068.
Conclusions: LeHaMoC is easily customized to model a variety of high-energy astrophysical sources and has the potential to become a widely utilized tool in multi-messenger astrophysics.
Instructions for downloading the code, accessing online documentation, and reproducing applications presented in this paper can be found at LeHaMoC Github repository.
[read the article; Stathopoulos et al.]
Theoretical Analysis of the RX J0209.6-7427 X-ray Spectrum During Giant Outburst
Feb 2024
We model the spectral formation occurring in the binary X-ray pulsar RX~J0209.6-7427 during the 2019 super-Eddington outburst. Using a theoretical model previously developed by the authors, we are able to produce spectra that closely resemble the phase-averaged X-ray spectra observed using NuSTAR and Insight-HXMT during low and high luminosity states of the outburst, respectively. The theoretical model simulates the accretion of fully ionized gas in a dipole magnetic field, and includes a complete description of the radiation hydrodynamics, matter distribution, and spectral formation. Type II X-ray outbursts provide an opportunity to study accretion over a large range of luminosities for the same neutron star. The analysis performed here represents the first time both the outburst low and high states of an accretion-powered X-ray pulsar are modeled using a physics-based model rather than standard phenomenological fitting with arbitrary mathematical functions. We find the outer polar cap radius remains constant and the column is more fully-filled with increasing luminosity, Comptonized bremsstrahlung dominates the formation of the phase-averaged X-ray spectrum, and a negative correlation exists between cyclotron centroid energy and luminosity, as expected. The super-Eddington nature of the outburst is rendered possible due to the low scattering cross section for photons propagating parallel to the magnetic field. We also find emission through the column top dominates in both the low and high states, implying the pulse profiles should have a roughly sinusoidal shape, which agrees with observed properties of ultra-luminous X-ray pulsars.
[read the article; West et al.]
The Luminosity Phase Space of Galactic and Extragalactic X-Ray Transients Out to Intermediate Redshifts
Dec 2023
We present a detailed compilation and analysis of the X-ray phase space of low- to intermediate-redshift (0 ≤ z ≤ 1) transients that consolidates observed light curves (and theory where necessary) for a large variety of classes of transient/variable phenomena in the 0.3-10 keV energy band. We include gamma-ray burst afterglows, supernovae, supernova shock breakouts and shocks interacting with the environment, tidal disruption events and active galactic nuclei, fast blue optical transients, cataclysmic variables, magnetar flares/outbursts and fast radio bursts, cool stellar flares, X-ray binary outbursts, and ultraluminous X-ray sources. Our overarching goal is to offer a comprehensive resource for the examination of these ephemeral events, extending the X-ray duration-luminosity phase space (DLPS) to show luminosity evolution. We use existing observations (both targeted and serendipitous) to characterize the behavior of various transient/variable populations. Contextualizing transient signals in the larger DLPS serves two primary purposes: to identify areas of interest (i.e., regions in the parameter space where one would expect detections, but in which observations have historically been lacking), and to provide initial qualitative guidance in classifying newly discovered transient signals. We find that while the most luminous (largely extragalactic) and least luminous (largely Galactic) part of the phase space is well populated at t > 0.1 days, intermediate-luminosity phenomena (L X = 1034-1042 erg s-1) represent a gap in the phase space. We thus identify L X = 1034-1042 erg s-1 and t = 10-4 to 0.1 days as a key discovery phase space in transient X-ray astronomy.
[read the article; Polzin et al.]
The high energy X-ray probe (HEX-P): a new window into neutron star accretion
Dec 2023
Accreting neutron stars (NSs) represent a unique laboratory for probing the physics of accretion in the presence of strong magnetic fields (B ≳ 108 G). Additionally, the matter inside the NS itself exists in an ultra-dense, cold state that cannot be reproduced in Earth-based laboratories. Hence, observational studies of these objects are a way to probe the most extreme physical regimes. Here we present an overview of the field and discuss the most important outstanding problems related to NS accretion. We show how these open questions regarding accreting NSs in both low-mass and high-mass X-ray binary systems can be addressed with the High-Energy X-ray Probe (HEX-P) via simulated data. In particular, with the broad X-ray passband and improved sensitivity afforded by a low X-ray background, HEX-P will be able to 1) distinguish between competing continuum emission models; 2) provide tighter upper limits on NS radii via reflection modeling techniques that are independent and complementary to other existing methods; 3) constrain magnetic field geometry, plasma parameters, and accretion column emission patterns by characterizing fundamental and harmonic cyclotron lines and exploring their behavior with pulse phase; 4) directly measure the surface magnetic field strength of highly magnetized NSs at the lowest accretion luminosities; as well as 5) detect cyclotron line features in extragalactic sources and probe their dependence on luminosity in the super-Eddington regime in order to distinguish between geometrical evolution and accretion-induced decay of the magnetic field. In these ways HEX-P will provide an essential new tool for exploring the physics of NSs, their magnetic fields, and the physics of extreme accretion.
[read the article; Ludlam et al.]
From Stellar Death to Cosmic Revelations: Zooming in on Compact Objects, Relativistic Outflows and Supernova Remnants with AXIS
Nov 2023
Compact objects and supernova remnants provide nearby laboratories to probe the fate of stars after they die, and the way they impact, and are impacted by, their surrounding medium. The past five decades have significantly advanced our understanding of these objects, and showed that they are most relevant to our understanding of some of the most mysterious energetic events in the distant Universe, including Fast Radio Bursts and Gravitational Wave sources. However, many questions remain to be answered. These include: What powers the diversity of explosive phenomena across the electromagnetic spectrum? What are the mass and spin distributions of neutron stars and stellar mass black holes? How do interacting compact binaries with white dwarfs - the electromagnetic counterparts to gravitational wave LISA sources - form and behave? Which objects inhabit the faint end of the X-ray luminosity function? How do relativistic winds impact their surroundings? What do neutron star kicks reveal about fundamental physics and supernova explosions? How do supernova remnant shocks impact cosmic magnetism? This plethora of questions will be addressed with AXIS - the Advanced X-ray Imaging Satellite - a NASA Probe Mission Concept designed to be the premier high-angular resolution X-ray mission for the next decade. AXIS, thanks to its combined (a) unprecedented imaging resolution over its full field of view, (b) unprecedented sensitivity to faint objects due to its large effective area and low background, and (c) rapid response capability, will provide a giant leap in discovering and identifying populations of compact objects (isolated and binaries), particularly in crowded regions such as globular clusters and the Galactic Center, while addressing science questions and priorities of the US Decadal Survey for Astronomy and Astrophysics (Astro2020).
[read the article; Safi-Harb et al.]
The high energy X-ray probe (HEX-P): studying extreme accretion with ultraluminous X-ray sources
Nov 2023
Ultraluminous X-ray sources (ULXs) represent an extreme class of accreting compact objects: from the identification of some of the accretors as neutron stars to the detection of powerful winds travelling at 0.1–0.2 c, the increasing evidence points towards ULXs harbouring stellar-mass compact objects undergoing highly super-Eddington accretion. Measuring their intrinsic properties, such as the accretion rate onto the compact object, the outflow rate, the masses of accretor/companion-hence their progenitors, lifetimes, and future evolution-is challenging due to ULXs being mostly extragalactic and in crowded fields. Yet ULXs represent our best opportunity to understand super-Eddington accretion physics and the paths through binary evolution to eventual double compact object binaries and gravitational-wave sources. Methods: Through a combination of end-to-end and single-source simulations, we investigate the ability of HEX-P to study ULXs in the context of their host galaxies and compare it to XMM-Newton and NuSTAR, the current instruments with the most similar capabilities.Results: HEX-P's higher sensitivity, which is driven by its narrow point-spread function and low background, allows it to detect pulsations and broad spectral features from ULXs better than XMM-Newton and NuSTAR.Discussion: We describe the value of HEX-P in understanding ULXs and their associated key physics, through a combination of broadband sensitivity, timing resolution, and angular resolution, which make the mission ideal for pulsation detection and low-background, broadband spectral studies.
[read the article; Bachetti et al.]
Constraining the Evolution of the Unstable Accretion Disk in SMC X-1 with NICER
Aug 2023
Neutron star high-mass X-ray binaries with superorbital modulations in luminosity host warped inner accretion disks that occult the neutron star during precession. In SMC X-1, the instability in the warped disk geometry causes superorbital period "excursions": times of instability when the superorbital period decreases from its typical value of 55 to ~40 days. Disk instability makes SMC X-1 an ideal system in which to investigate the effects of variable disk geometry on the inner accretion flow. Using the high-resolution spectral and timing capabilities of the Neutron Star Interior Composition Explorer, we examined the high state of four different superorbital cycles of SMC X-1 to search for changes in spectral shape and connections to the unstable disk geometry. We performed pulse phase-averaged and phase-resolved spectroscopy to closely compare the changes in spectral shape and any cycle-to-cycle variations. While some parameters, including the photon index and absorbing column density, show slight variations with superorbital phase, these changes are most evident during the intermediate state of the superorbital cycle. Few spectral changes are observed within the high state of the superorbital cycle, possibly indicating the disk instability does not significantly change SMC X-1's accretion process.
[read the article; Brumback et al.]
On the cyclotron absorption line and evidence of the spectral transition in SMC X-2 during 2022 giant outburst
May 2023
We report comprehensive spectral and temporal properties of the Be/X-ray binary pulsar SMC X-2 using X-ray observations during the 2015 and 2022 outbursts. The pulse profile of the pulsar is unique and strongly luminosity dependent. It evolves from a broad-humped into a double-peaked profile above luminosity 3 × 1038 erg s-1. The pulse fraction of the pulsar is found to be a linear function of luminosity as well as energy. We also studied the spectral evolution of the source during the latest 2022 outburst with NICER. The observed photon index shows a negative and positive correlation below and above the critical luminosity, respectively, suggesting evidence of spectral transition from the sub-critical to supercritical regime. The broad-band spectroscopy of four sets of NuSTAR and XRT/NICER data from both outbursts can be described using a cut-off power-law model with a blackbody component. In addition to the 6.4 keV iron fluorescence line, an absorption-like feature is clearly detected in the spectra. The cyclotron line energy observed during the 2015 outburst is below 29.5 keV, however latest estimates in the 2022 outburst suggest a value of 31.5 keV. Moreover, an increase of 3.4 keV is detected in the cyclotron line energy at equal levels of luminosity observed in 2022 with respect to 2015. The observed cyclotron line energy variation is explored in terms of accretion induced screening mechanism or geometrical variation in line forming region.
[read the article; Jaisawal et al.]
A shared accretion instability for black holes and neutron stars
March 2023
Accretion disks around compact objects are expected to enter an unstable phase at high luminosity1. One instability may occur when the radiation pressure generated by accretion modifies the disk viscosity, resulting in the cyclic depletion and refilling of the inner disk on short timescales2. Such a scenario, however, has only been quantitatively verified for a single stellar-mass black hole3-5. Although there are hints of these cycles in a few isolated cases6-10, their apparent absence in the variable emission of most bright accreting neutron stars and black holes has been a continuing puzzle11. Here we report the presence of the same multiwavelength instability around an accreting neutron star. Moreover, we show that the variability across the electromagnetic spectrum—from radio to X-ray—of both black holes and neutron stars at high accretion rates can be explained consistently if the accretion disks are unstable, producing relativistic ejections during transitions that deplete or refill the inner disk. Such a new association allows us to identify the main physical components responsible for the fast multiwavelength variability of highly accreting compact objects.
[read the article; Vincentelli et al.]
A Bayesian approach for torque modelling of BeXRB pulsars with application to super-Eddington accretors
March 2023
In this study, we present a method to estimate posterior distributions for standard accretion torque model parameters and binary orbital parameters for X-ray binaries using a nested sampling algorithm for Bayesian parameter estimation. We study the spin evolution of two Be X-ray binary systems in the Magellanic Clouds, RX J0520.5-6932 and RX J0209-7427, during major outbursts, in which they surpassed the Eddington limit. Moreover, we apply our method to the recently discovered Swift J0243.6+6124, the only known Galactic pulsating ultra-luminous X-ray source. This is an excellent candidate for studying the disc evolution at super-Eddington accretion rates, because its luminosity spans several orders of magnitude during its outburst, with a maximum LX that exceeded the Eddington limit by a factor of ~10. Our method, when applied to RX J0520.5-6932 and RX J0209-7427, is able to identify the more favourable torque model for each system, while yielding meaningful ranges for the NS and orbital parameters. Our analysis for Swift J0243.6+6124 illustrates that, contrary to the standard torque model predictions, the magnetospheric radius (Rm) and the Alfvén radius (RA) are not proportional to each other when surpassing the Eddington limit. Reported distance estimates of this source range between 5 and 7 kpc. Smaller distances require non-typical neutron star properties (i.e. mass and radius) and possibly lower radiative efficiency of the accretion column.
[read the article; Karaferias et al.]