Showing 1–50 of 462 results for author: Stella, L

Discovery of Polarized X-Ray Emission from the Accreting Millisecond Pulsar SRGA J144459.2-604207

Authors: Alessandro Papitto, Alessandro Di Marco, Juri Poutanen, Tuomo Salmi, Giulia Illiano, Fabio La Monaca, Filippo Ambrosino, Anna Bobrikova, Maria Cristina Baglio, Caterina Ballocco, Luciano Burderi, Sergio Campana, Francesco Coti Zelati, Tiziana Di Salvo, Riccardo La Placa, Vladislav Loktev, Sinan Long, Christian Malacaria, Arianna Miraval Zanon, Mason Ng, Maura Pilia, Andrea Sanna, Luigi Stella, Tod Strohmayer, Silvia Zane

Abstract: We report on the discovery of polarized X-ray emission from an accreting millisecond pulsar. During a 10-day-long coverage of the February 2024 outburst of SRGA J144459.2-604207, the Imaging X-ray Polarimetry Explorer (IXPE) detected an average polarization degree of the 2-8 keV emission of 2.3% +/- 0.4% at an angle of 59° +/- 6° (East of North; uncertainties quoted at the 1$σ$ confidence level).… ▽ More We report on the discovery of polarized X-ray emission from an accreting millisecond pulsar. During a 10-day-long coverage of the February 2024 outburst of SRGA J144459.2-604207, the Imaging X-ray Polarimetry Explorer (IXPE) detected an average polarization degree of the 2-8 keV emission of 2.3% +/- 0.4% at an angle of 59° +/- 6° (East of North; uncertainties quoted at the 1$σ$ confidence level). The polarized signal shows a significant energy dependence with a degree of 4.0% +/- 0.5% between 3 and 6 keV and < 2% (90% c.l.) in the 2-3 keV range. We used NICER, XMM-Newton, and NuSTAR observations to obtain an accurate pulse timing solution and perform a phase-resolved polarimetric analysis of IXPE data. We did not detect any significant variability of the Stokes parameters Q and U with the spin and the orbital phases. We used the relativistic rotating vector model to show that a moderately fan-beam emission from two point-like spots at a small magnetic obliquity ($\simeq$ 10°) is compatible with the observed pulse profile and polarization properties. IXPE also detected 52 type-I X-ray bursts, with a recurrence time $Δt_{rec}$ increasing from 2 to 8 h as a function of the observed count rate $C$ as as $Δt_{rec} \simeq C^{-0.8}$ We stacked the emission observed during all the bursts and obtained an upper limit on the polarization degree of 8.5% (90% c.l.). △ Less

Submitted 1 August, 2024; originally announced August 2024.

Comments: 13 pages, 9 figures, submitted to ApJL

Skipping a beat: discovery of persistent quasi-periodic oscillations associated with pulsed fraction drop of the spin signal in M51 ULX-7

Authors: Matteo Imbrogno, Sara Elisa Motta, Roberta Amato, Gian Luca Israel, Guillermo Andres Rodríguez Castillo, Murray Brightman, Piergiorgio Casella, Matteo Bachetti, Felix Fürst, Luigi Stella, Ciro Pinto, Fabio Pintore, Francesco Tombesi, Andrés Gúrpide, Matthew J. Middleton, Chiara Salvaggio, Andrea Tiengo, Andrea Belfiore, Andrea De Luca, Paolo Esposito, Anna Wolter, Hannah P. Earnshaw, Dominic J. Walton, Timothy P. Roberts, Luca Zampieri , et al. (2 additional authors not shown)

Abstract: The discovery of pulsations in (at least) six ultraluminous X-ray sources (ULXs) has shown that neutron stars can accrete at (highly) super-Eddington rates, challenging the standard accretion theories. M51 ULX-7, with a spin signal of $P\simeq2.8$ s, is the pulsating ULX (PULX) with the shortest known orbital period ($P_\mathrm{orb}\simeq2$ d) and has been observed multiple times by XMM-Newton, Ch… ▽ More The discovery of pulsations in (at least) six ultraluminous X-ray sources (ULXs) has shown that neutron stars can accrete at (highly) super-Eddington rates, challenging the standard accretion theories. M51 ULX-7, with a spin signal of $P\simeq2.8$ s, is the pulsating ULX (PULX) with the shortest known orbital period ($P_\mathrm{orb}\simeq2$ d) and has been observed multiple times by XMM-Newton, Chandra, and NuSTAR. We report on the timing and spectral analyses of three XMM-Newton observations of M51 ULX-7 performed between the end of 2021 and the beginning of 2022, together with a timing re-analysis of XMM-Newton, Chandra, and NuSTAR archival observations. We investigated the spin signal by applying accelerated search techniques and studied the power spectrum through the fast Fourier transform, looking for (a)periodic variability in the source flux. We analysed the energy spectra of the 2021-2022 observations and compared them to the older ones. We report the discovery of a recurrent, significant ($>$3$σ$) broad complex at mHz frequencies in the power spectra of M51 ULX-7. We did not detect the spin signal, setting a 3$σ$ upper limit on the pulsed fraction of $\lesssim10\%$ for the single observation. The complex is significantly detected also in five Chandra observations performed in 2012. M51 ULX-7 represents the second PULX for which we have a significant detection of mHz-QPOs at super-Eddington luminosities. These findings suggest that one should avoid using the observed QPO frequency to infer the mass of the accretor in a ULX. The absence of spin pulsations when the broad complex is detected suggests that the mechanism responsible for the aperiodic modulation also dampens the spin signal's pulsed fraction. If true, this represents an additional obstacle in the detection of new PULXs, suggesting an even larger occurrence of PULXs among ULXs. △ Less

Submitted 12 July, 2024; originally announced July 2024.

Comments: 14 pages (12 main text + Appendix), 5 figures. Accepted for publication on A&A

Universal and non-universal signatures in the scaling functions of critical variables

Authors: Gianluca Teza, Attilio L. Stella

Abstract: The view that the probability density function (PDF) of a key statistical variable, anomalously scaled by size or time, could furnish a hallmark of universal behavior contrasts with the circumstance that such density sensibly depends on non-universal features. We solve this apparent contradiction by demonstrating that both non-universal amplitudes and universal exponents of leading critical singul… ▽ More The view that the probability density function (PDF) of a key statistical variable, anomalously scaled by size or time, could furnish a hallmark of universal behavior contrasts with the circumstance that such density sensibly depends on non-universal features. We solve this apparent contradiction by demonstrating that both non-universal amplitudes and universal exponents of leading critical singularities in large deviation functions are determined by the PDF tails, whose form is argued on extensivity. This unexplored scenario implies a universal form of central limit theorem at criticality and is confirmed by exact calculations for mean field Ising models in equilibrium and for anomalous diffusion models. △ Less

Submitted 10 July, 2024; originally announced July 2024.

Investigating the interplay between the coronal properties and the hard X-ray variability of active galactic nuclei with NuSTAR

Authors: Roberto Serafinelli, Alessandra De Rosa, Alessia Tortosa, Luigi Stella, Fausto Vagnetti, Stefano Bianchi, Claudio Ricci, Elias Kammoun, Pierre-Olivier Petrucci, Riccardo Middei, Giorgio Lanzuisi, Andrea Marinucci, Francesco Ursini, Giorgio Matt

Abstract: Active galactic nuclei (AGN) are extremely variable in the X-ray band down to very short timescales. However, the driver behind the X-ray variability is still poorly understood. Previous results suggest that the hot corona responsible for the primary Comptonized emission observed in AGN is expected to play an important role in driving the X-ray variability. In this work, we investigate the connect… ▽ More Active galactic nuclei (AGN) are extremely variable in the X-ray band down to very short timescales. However, the driver behind the X-ray variability is still poorly understood. Previous results suggest that the hot corona responsible for the primary Comptonized emission observed in AGN is expected to play an important role in driving the X-ray variability. In this work, we investigate the connection between the X-ray amplitude variability and the coronal physical parameters; namely, the temperature ($kT$) and optical depth ($τ$). We present the spectral and timing analysis of 46 {\it NuSTAR} observations corresponding to a sample of 20 AGN. For each source, we derived the coronal temperature and optical depth through X-ray spectroscopy and computed the normalized excess variance for different energy bands on a timescale of $10$ ks. We find a strong inverse correlation between $kT$ and $τ$, with correlation coefficient of $r<-0.9$ and negligible null probability. No clear dependence was found among the temperature and physical properties, such as the black hole mass or the Eddington ratio. We also see that the observed X-ray variability is not correlated with either the coronal temperature or optical depth under the thermal equilibrium assumption, whereas it is anticorrelated with the black hole mass. These results can be interpreted through a scenario where the observed X-ray variability could primarily be driven by variations in the coronal physical properties on a timescale of less than $10$~ks; whereas we assume thermal equilibrium on such timescales in this work, given the capability of the currently available hard X-ray telescopes. Alternatively, it is also possible that the X-ray variability is mostly driven by the absolute size of the corona, which depends on the supermassive black hole mass, rather than resulting from any of its physical properties. △ Less

Submitted 24 July, 2024; v1 submitted 9 July, 2024; originally announced July 2024.

Comments: Accepted for publication in Astronomy & Astrophysics, 21 pages, 9 figures, 4 tables

Gravitational self-lensing of Fast Radio Bursts in neutron star magnetospheres: I. The model

Authors: Simone Dall'Osso, Riccardo La Placa, Luigi Stella, Pavel Bakala, Andrea Possenti

Abstract: Fast Radio Bursts (FRBs) are cosmological sub-second bursts of coherent radio emission, whose source is still unknown. To date, the galactic magnetar SGR 1935+2154 is the only astrophysical object known to emit radio bursts akin to FRBs, albeit less powerful, supporting suggestions that FRBs originate from magnetars. Many remarkable properties of FRBs, e.g. the dichotomy between repeaters and one-… ▽ More Fast Radio Bursts (FRBs) are cosmological sub-second bursts of coherent radio emission, whose source is still unknown. To date, the galactic magnetar SGR 1935+2154 is the only astrophysical object known to emit radio bursts akin to FRBs, albeit less powerful, supporting suggestions that FRBs originate from magnetars. Many remarkable properties of FRBs, e.g. the dichotomy between repeaters and one-off sources, and their power-law energy distributions (with typical index $\sim 2-3$), are not well understood yet. Moreover, the huge radio power released by the most active repeaters is challenging even for the magnetic energy reservoir of magnetars. Here we assume that FRBs originate from co-rotating hot-spots anchored in neutron star magnetospheres and get occasionally amplified by large factors via gravitational self-lensing in the strong NS field. We evaluate the probability of amplification and show that (i) a power-law energy distribution of events $\propto E^{-(2- 3)}$ is generally expected, (ii) all FRB sources may be regarded as repeating, their appearance as one-off sources or repeaters being determined by the critical dependence of the amplification probability on the emission geometry and source orientation relative to Earth and (iii) the most active repeaters, in particular, correspond to extremely rare and finely-tuned orientations ($\sim$ one in $10^6$), leading to large probabilities of amplification which make their bursts frequently detectable. At the same time, their power release appears enhanced, typically by factors $\gtrsim 10$, easing their energy budget problem. △ Less

Submitted 4 July, 2024; originally announced July 2024.

Comments: 23 pages (17 main text + Appendices), 13 Figures, accepted for publication on The Astrophysical Journal

An Evolutionary Framework for Connect-4 as Test-Bed for Comparison of Advanced Minimax, Q-Learning and MCTS

Authors: Henry Taylor, Leonardo Stella

Abstract: A major challenge in decision making domains with large state spaces is to effectively select actions which maximize utility. In recent years, approaches such as reinforcement learning (RL) and search algorithms have been successful to tackle this issue, despite their differences. RL defines a learning framework that an agent explores and interacts with. Search algorithms provide a formalism to se… ▽ More A major challenge in decision making domains with large state spaces is to effectively select actions which maximize utility. In recent years, approaches such as reinforcement learning (RL) and search algorithms have been successful to tackle this issue, despite their differences. RL defines a learning framework that an agent explores and interacts with. Search algorithms provide a formalism to search for a solution. However, it is often difficult to evaluate the performances of such approaches in a practical way. Motivated by this problem, we focus on one game domain, i.e., Connect-4, and develop a novel evolutionary framework to evaluate three classes of algorithms: RL, Minimax and Monte Carlo tree search (MCTS). The contribution of this paper is threefold: i) we implement advanced versions of these algorithms and provide a systematic comparison with their standard counterpart, ii) we develop a novel evaluation framework, which we call the Evolutionary Tournament, and iii) we conduct an extensive evaluation of the relative performance of each algorithm to compare our findings. We evaluate different metrics and show that MCTS achieves the best results in terms of win percentage, whereas Minimax and Q-Learning are ranked in second and third place, respectively, although the latter is shown to be the fastest to make a decision. △ Less

Submitted 26 May, 2024; originally announced May 2024.

Comments: 8 pages, 4 figures, 1 table

The Orbit of NGC 5907 ULX-1

Authors: Andrea Belfiore, Ruben Salvaterra, Lara Sidoli, Gian Luca Israel, Luigi Stella, Andrea De Luca, Sandro Mereghetti, Paolo Esposito, Fabio Pintore, Antonino D'Aì, Guillermo Rodrìguez Castillo, Dominic J. Walton, Felix Fürst, Danilo Magistrali, Anna Wolter, Matteo Imbrogno

Abstract: We report on the orbit of the binary system powering the most extreme ultraluminous X-ray pulsar known to date: NGC 5907 ULX-1 (hereafter ULX1). ULX1 has been the target of a substantial multi-instrument campaign, mainly in the X-ray band, but no clear counterparts are known in other bands. Although ULX1 is highly variable and pulsations can be transient (regardless of the source flux), the timing… ▽ More We report on the orbit of the binary system powering the most extreme ultraluminous X-ray pulsar known to date: NGC 5907 ULX-1 (hereafter ULX1). ULX1 has been the target of a substantial multi-instrument campaign, mainly in the X-ray band, but no clear counterparts are known in other bands. Although ULX1 is highly variable and pulsations can be transient (regardless of the source flux), the timing data collected so far allow us to investigate the orbit of this system. We find an orbital period $P_{orb}=5.7^{+0.1}_{-0.6}\text{ d}$ and a projected semi-axis $A_1 =3.1^{+0.8}_{-0.9}\text{ lts}$. The most likely ephemeris is: $P_{orb}=5.6585(6)\text{ d}$, $A_1 = 3.1(4)\text{ lts}$, and the epoch of ascending nodes passage is: $T_{asc} = 57751.37(5)\text{ MJD}$. However, there are 6 similar solutions, acceptable within $3\,σ$. We find further indications that ULX1 is a high-mass X-ray binary. This implies that we are observing its orbit face-on, with an inclination $<5\text{ deg}$. △ Less

Submitted 7 May, 2024; originally announced May 2024.

Comments: 11 pages, 2 tables, 6 figures

Journal ref: ApJ 965, 78B (2024)

Chronos: Learning the Language of Time Series

Authors: Abdul Fatir Ansari, Lorenzo Stella, Caner Turkmen, Xiyuan Zhang, Pedro Mercado, Huibin Shen, Oleksandr Shchur, Syama Sundar Rangapuram, Sebastian Pineda Arango, Shubham Kapoor, Jasper Zschiegner, Danielle C. Maddix, Hao Wang, Michael W. Mahoney, Kari Torkkola, Andrew Gordon Wilson, Michael Bohlke-Schneider, Yuyang Wang

Abstract: We introduce Chronos, a simple yet effective framework for pretrained probabilistic time series models. Chronos tokenizes time series values using scaling and quantization into a fixed vocabulary and trains existing transformer-based language model architectures on these tokenized time series via the cross-entropy loss. We pretrained Chronos models based on the T5 family (ranging from 20M to 710M… ▽ More We introduce Chronos, a simple yet effective framework for pretrained probabilistic time series models. Chronos tokenizes time series values using scaling and quantization into a fixed vocabulary and trains existing transformer-based language model architectures on these tokenized time series via the cross-entropy loss. We pretrained Chronos models based on the T5 family (ranging from 20M to 710M parameters) on a large collection of publicly available datasets, complemented by a synthetic dataset that we generated via Gaussian processes to improve generalization. In a comprehensive benchmark consisting of 42 datasets, and comprising both classical local models and deep learning methods, we show that Chronos models: (a) significantly outperform other methods on datasets that were part of the training corpus; and (b) have comparable and occasionally superior zero-shot performance on new datasets, relative to methods that were trained specifically on them. Our results demonstrate that Chronos models can leverage time series data from diverse domains to improve zero-shot accuracy on unseen forecasting tasks, positioning pretrained models as a viable tool to greatly simplify forecasting pipelines. △ Less

Submitted 2 May, 2024; v1 submitted 12 March, 2024; originally announced March 2024.

Comments: Code and model checkpoints available at https://github.com/amazon-science/chronos-forecasting

Deep Non-Parametric Time Series Forecaster

Authors: Syama Sundar Rangapuram, Jan Gasthaus, Lorenzo Stella, Valentin Flunkert, David Salinas, Yuyang Wang, Tim Januschowski

Abstract: This paper presents non-parametric baseline models for time series forecasting. Unlike classical forecasting models, the proposed approach does not assume any parametric form for the predictive distribution and instead generates predictions by sampling from the empirical distribution according to a tunable strategy. By virtue of this, the model is always able to produce reasonable forecasts (i.e.,… ▽ More This paper presents non-parametric baseline models for time series forecasting. Unlike classical forecasting models, the proposed approach does not assume any parametric form for the predictive distribution and instead generates predictions by sampling from the empirical distribution according to a tunable strategy. By virtue of this, the model is always able to produce reasonable forecasts (i.e., predictions within the observed data range) without fail unlike classical models that suffer from numerical stability on some data distributions. Moreover, we develop a global version of the proposed method that automatically learns the sampling strategy by exploiting the information across multiple related time series. The empirical evaluation shows that the proposed methods have reasonable and consistent performance across all datasets, proving them to be strong baselines to be considered in one's forecasting toolbox. △ Less

Submitted 22 December, 2023; originally announced December 2023.

Sub-second infrared variability from the archetypal accreting neutron star 4U~1728-34

Authors: F. M. Vincentelli, P. Casella, A. Borghese, Y. Cavecchi, G. Mastroserio, L. Stella, D. Altamirano, M. Armas Padilla, M. C. Baglio, T. M. Belloni, J. Casares, V. A. Cúneo, N. Degenaar, M. Díaz Trigo, R. Fender, T. Maccarone, J. Malzac, D. Mata Sánchez, M. Middleton, S. Migliari, T. Muñoz-Darias, K. O'Brien, G. Panizo-Espinar, J. Sánchez-Sierras, D. M. Russell , et al. (1 additional authors not shown)

Abstract: We report on the first simultaneous high-time resolution X-ray and infrared (IR) observations of a neutron star low mass X-ray binary in its hard state. We performed $\approx 2\,$h of simultaneous observations of 4U 1728-34 using HAWK-I@VLT, XMM-Newton and NuSTAR. The source displayed significant X-ray and IR variability down to sub-second timescales. By measuring the cross-correlation function be… ▽ More We report on the first simultaneous high-time resolution X-ray and infrared (IR) observations of a neutron star low mass X-ray binary in its hard state. We performed $\approx 2\,$h of simultaneous observations of 4U 1728-34 using HAWK-I@VLT, XMM-Newton and NuSTAR. The source displayed significant X-ray and IR variability down to sub-second timescales. By measuring the cross-correlation function between the infrared and X-ray lightcurves, we discovered a significant correlation with an infrared lead of $\approx 30-40\,$ms with respect to the X-rays. We analysed the X-ray energy dependence of the lag, finding a marginal increase towards higher energies. Given the sign of the lag, we interpret this as possible evidence of Comptonization from external seed photons. We discuss the origin of the IR seed photons in terms of cyclo-synchrotron radiation from an extended hot flow. Finally, we also observed the IR counterpart of a type-I X-ray burst, with a delay of $\approx7.2\,$s. Although some additional effects may be at play, by assuming that this lag is due to light travel time between the central object and the companion star, we find that 4U 1728-34 must have an orbital period longer than $3\,$h and an inclination higher than 8$^\circ$. △ Less

Submitted 29 August, 2023; originally announced August 2023.

Comments: Accepted for publication in MNRAS, 11 pages 7 figures

Discovery of a magnetar candidate X-ray pulsar in the Large Magellanic Cloud

Authors: M. Imbrogno, G. L. Israel, G. A. Rodríguez Castillo, D. A. H. Buckley, F. Coti Zelati, N. Rea, I. M. Monageng, P. Casella, L. Stella, F. Haberl, P. Esposito, F. Tombesi, A. De Luca, A. Tiengo

Abstract: During a systematic search for new X-ray pulsators in the XMM-Newton archive, we discovered a high amplitude ($PF\simeq86\%$) periodic ($P\simeq7.25\,\mathrm{s}$) modulation in the X-ray flux of 4XMM J045626.3-694723 (J0456 hereafter), a previously unclassified source in the Large Magellanic Cloud (LMC). The period of the modulation is strongly suggestive of a spinning neutron star (NS). The sourc… ▽ More During a systematic search for new X-ray pulsators in the XMM-Newton archive, we discovered a high amplitude ($PF\simeq86\%$) periodic ($P\simeq7.25\,\mathrm{s}$) modulation in the X-ray flux of 4XMM J045626.3-694723 (J0456 hereafter), a previously unclassified source in the Large Magellanic Cloud (LMC). The period of the modulation is strongly suggestive of a spinning neutron star (NS). The source was detected only during one out of six observations in 2018-2022. Based on an absorbed power-law spectral model with photon slope of $Γ\simeq 1.9$, we derive a 0.3-10 keV luminosity of $L_\mathrm{X}\simeq2.7\times10^{34}$ erg cm$^{-2}$ s$^{-1}$ for a distance of 50 kpc. The X-ray properties of J0456 are at variance with those of variable LMC X-ray pulsars hosted in high-mass X-ray binary systems with a Be-star companion. Based on SALT spectroscopic observations of the only optical object that matches the X-ray uncertainty region, we cannot completely rule out that J0456 is a NS accreting from a late-type (G8-K3) star, an as-yet-unobserved binary evolutionary outcome in the MCs. We show that the source properties are in better agreement with those of magnetars. J0456 may thus be second known magnetar in the LMC after SGR 0526-66. △ Less

Submitted 18 July, 2023; originally announced July 2023.

Comments: 9 pages, 6 figures, 1 table. Accepted for publication in MNRAS

Cascading Failures in the Global Financial System: A Dynamical Model

Authors: Leonardo Stella, Dario Bauso, Franco Blanchini, Patrizio Colaneri

Abstract: In this paper, we propose a dynamical model to capture cascading failures among interconnected organizations in the global financial system. Failures can take the form of bankruptcies, defaults, and other insolvencies. The network that underpins the financial interdependencies between different organizations constitutes the backbone of the financial system. A failure in one or more of these organi… ▽ More In this paper, we propose a dynamical model to capture cascading failures among interconnected organizations in the global financial system. Failures can take the form of bankruptcies, defaults, and other insolvencies. The network that underpins the financial interdependencies between different organizations constitutes the backbone of the financial system. A failure in one or more of these organizations can lead the propagation of the financial collapse onto other organizations in a domino effect. Paramount importance is therefore given to the mitigation of these failures. Motivated by the relevance of this problem and recent prominent events connected to it, we develop a framework that allows us to investigate under what conditions organizations remain healthy or are involved in the propagation of the failures in the network. The contribution of this paper is the following: i) we develop a dynamical model that describes the equity values of financial organizations and their evolution over time given an initial condition; ii) we characterize the equilibria for this model by proving the existence and uniqueness of these equilibria, and by providing an explicit expression for them; and iii) we provide a computational method via sign-space iteration to analyze the propagation of failures and the attractive equilibrium point. △ Less

Submitted 10 November, 2023; v1 submitted 7 July, 2023; originally announced July 2023.

First detection of X-ray polarization from the accreting neutron star 4U 1820-303

Authors: Alessandro Di Marco, Fabio La Monaca, Juri Poutanen, Thomas D. Russell, Alessio Anitra, Ruben Farinelli, Guglielmo Mastroserio, Fabio Muleri, Fei Xie, Matteo Bachetti, Luciano Burderi, Francesco Carotenuto, Melania Del Santo, Tiziana Di Salvo, Michal Dovciak, Andrea Gnarini, Rosario Iaria, Jari J. E. Kajava, Kuan Liu, Riccardo Middei, Stephen L. O'Dell, Maura Pilia, John Rankin, Andrea Sanna, Jakob van den Eijnden , et al. (94 additional authors not shown)

Abstract: This paper reports the first detection of polarization in the X-rays for atoll-source 4U 1820-303, obtained with the Imaging X-ray Polarimetry Explorer (IXPE) at 99.999% confidence level (CL). Simultaneous polarimetric measurements were also performed in the radio with the Australia Telescope Compact Array (ATCA). The IXPE observations of 4U 1820-303 were coordinated with Swift-XRT, NICER, and NuS… ▽ More This paper reports the first detection of polarization in the X-rays for atoll-source 4U 1820-303, obtained with the Imaging X-ray Polarimetry Explorer (IXPE) at 99.999% confidence level (CL). Simultaneous polarimetric measurements were also performed in the radio with the Australia Telescope Compact Array (ATCA). The IXPE observations of 4U 1820-303 were coordinated with Swift-XRT, NICER, and NuSTAR aiming to obtain an accurate X-ray spectral model covering a broad energy interval. The source shows a significant polarization above 4 keV, with a polarization degree of 2.0(0.5)% and a polarization angle of -55(7) deg in the 4-7 keV energy range, and a polarization degree of 10(2)% and a polarization angle of -67(7) deg in the 7-8 keV energy bin. This polarization also shows a clear energy trend with polarization degree increasing with energy and a hint for a position-angle change of about 90 deg at 96% CL around 4 keV. The spectro-polarimetric fit indicates that the accretion disk is polarized orthogonally to the hard spectral component, which is presumably produced in the boundary/spreading layer. We do not detect linear polarization from the radio counterpart, with a 99.97% upper limit of 50% at 7.25 GHz. △ Less

Submitted 21 August, 2023; v1 submitted 14 June, 2023; originally announced June 2023.

Journal ref: ApJL 953 L22 (2023)

Magnetar central engines in gamma-ray busts follow the universal relation of accreting magnetic stars

Authors: Simone Dall'Osso, Giulia Stratta, Rosalba Perna, Giovanni de Cesare, Luigi Stella

Abstract: Gamma-ray bursts (GRBs), both long and short, are explosive events whose inner engine is generally expected to be a black hole or a highly magnetic neutron star (magnetar) accreting high density matter. Recognizing the nature of GRB central engines, and in particular the formation of neutron stars (NSs), is of high astrophysical significance. A possible signature of NSs in GRBs is the presence of… ▽ More Gamma-ray bursts (GRBs), both long and short, are explosive events whose inner engine is generally expected to be a black hole or a highly magnetic neutron star (magnetar) accreting high density matter. Recognizing the nature of GRB central engines, and in particular the formation of neutron stars (NSs), is of high astrophysical significance. A possible signature of NSs in GRBs is the presence of a plateau in the early X-ray afterglow. Here we carefully select a subset of long and short GRBs with a clear plateau, and look for an additional NS signature in their prompt emission, namely a transition between accretion and propeller in analogy with accreting, magnetic compact objects in other astrophysical sources. We estimate from the prompt emission the minimum accretion luminosity below which the propeller mechanism sets in, and the NS magnetic field and spin period from the plateau. We demonstrate that these three quantities obey the same universal relation in GRBs as in other accreting compact objects switching from accretion to propeller. This relation provides also an estimate of the radiative efficiency of GRBs, which we find to be several times lower than radiatively efficient accretion in X-ray binaries and in agreement with theoretical expectations. These results provide additional support to the idea that at least some GRBs are powered by magnetars surrounded by an accretion disc. △ Less

Submitted 28 April, 2023; originally announced May 2023.

Comments: 15 pages, 5 figures, accepted for publication in The Astrophysical Journal Letters

Adaptive proximal algorithms for convex optimization under local Lipschitz continuity of the gradient

Authors: Puya Latafat, Andreas Themelis, Lorenzo Stella, Panagiotis Patrinos

Abstract: Backtracking linesearch is the de facto approach for minimizing continuously differentiable functions with locally Lipschitz gradient. In recent years, it has been shown that in the convex setting it is possible to avoid linesearch altogether, and to allow the stepsize to adapt based on a local smoothness estimate without any backtracks or evaluations of the function value. In this work we propose… ▽ More Backtracking linesearch is the de facto approach for minimizing continuously differentiable functions with locally Lipschitz gradient. In recent years, it has been shown that in the convex setting it is possible to avoid linesearch altogether, and to allow the stepsize to adapt based on a local smoothness estimate without any backtracks or evaluations of the function value. In this work we propose an adaptive proximal gradient method, adaPG, that uses novel estimates of the local smoothness modulus which leads to less conservative stepsize updates and that can additionally cope with nonsmooth terms. This idea is extended to the primal-dual setting where an adaptive three-term primal-dual algorithm, adaPD, is proposed which can be viewed as an extension of the PDHG method. Moreover, in this setting the "essentially" fully adaptive variant adaPD$^+$ is proposed that avoids evaluating the linear operator norm by invoking a backtracking procedure, that, remarkably, does not require extra gradient evaluations. Numerical simulations demonstrate the effectiveness of the proposed algorithms compared to the state of the art. △ Less

Submitted 13 March, 2024; v1 submitted 11 January, 2023; originally announced January 2023.

MSC Class: 65K05; 90C06; 90C25; 90C30; 90C47

Timing analysis of the 2022 outburst of the accreting millisecond X-ray pulsar SAX J1808.4$-$3658: hints of an orbital shrinking

Authors: Giulia Illiano, Alessandro Papitto, Andrea Sanna, Peter Bult, Filippo Ambrosino, Arianna Miraval Zanon, Francesco Coti Zelati, Luigi Stella, Diego Altamirano, Maria Cristina Baglio, Enrico Bozzo, Luciano Burderi, Domitilla de Martino, Alessandro Di Marco, Tiziana di Salvo, Carlo Ferrigno, Vladislav Loktev, Alessio Marino, Mason Ng, Maura Pilia, Juri Poutanen, Tuomo Salmi

Abstract: We present a pulse timing analysis of NICER observations of the accreting millisecond X-ray pulsar SAX J1808.4$-$3658 during the outburst that started on 2022 August 19. Similar to previous outbursts, after decaying from a peak luminosity of $\simeq 1\times10^{36} \, \mathrm{erg \, s^{-1}}$ in about a week, the pulsar entered in a $\sim 1$ month-long reflaring stage. Comparison of the average puls… ▽ More We present a pulse timing analysis of NICER observations of the accreting millisecond X-ray pulsar SAX J1808.4$-$3658 during the outburst that started on 2022 August 19. Similar to previous outbursts, after decaying from a peak luminosity of $\simeq 1\times10^{36} \, \mathrm{erg \, s^{-1}}$ in about a week, the pulsar entered in a $\sim 1$ month-long reflaring stage. Comparison of the average pulsar spin frequency during the outburst with those previously measured confirmed the long-term spin derivative of $\dotν_{\textrm{SD}}=-(1.15\pm0.06)\times 10^{-15} \, \mathrm{Hz\,s^{-1}}$, compatible with the spin-down torque of a $\approx 10^{26} \, \mathrm{G \, cm^3}$ rotating magnetic dipole. For the first time in the last twenty years, the orbital phase evolution shows evidence for a decrease of the orbital period. The long-term behaviour of the orbit is dominated by a $\sim 11 \, \mathrm{s}$ modulation of the orbital phase epoch consistent with a $\sim 21 \, \mathrm{yr}$ period. We discuss the observed evolution in terms of a coupling between the orbit and variations in the mass quadrupole of the companion star. △ Less

Submitted 22 December, 2022; v1 submitted 19 December, 2022; originally announced December 2022.

Comments: 11 pages, 3 figures, 1 table. Accepted for publication in ApJ Letters

Thermoelectric properties of cement composite analogues from first principles calculations

Authors: Esther Orisakwe, Conrad Johnston, Ruchita Jani, Xiaoli Liu, Lorenzo Stella, Jorge Kohanoff, Niall Holmes, Brian Norton, Ming Qu, Hongxi Yin, Kazuaki Yazawa

Abstract: Buildings are responsible for a considerable fraction of the energy wasted globally every year, and as a result, excess carbon emissions. While heat is lost directly in colder months and climates, resulting in increased heating loads, in hot climates cooling and ventilation is required. One avenue towards improving the energy efficiency of buildings is to integrate thermoelectric devices and mater… ▽ More Buildings are responsible for a considerable fraction of the energy wasted globally every year, and as a result, excess carbon emissions. While heat is lost directly in colder months and climates, resulting in increased heating loads, in hot climates cooling and ventilation is required. One avenue towards improving the energy efficiency of buildings is to integrate thermoelectric devices and materials within the fabric of the building to exploit the temperature gradient between the inside and outside to do useful work. Cement-based materials are ubiquitous in modern buildings and present an interesting opportunity to be functionalised. We present a systematic investigation of the electronic transport coefficients relevant to the thermoelectric materials of the calcium silicate hydrate (C-S-H) gel analogue, tobermorite, using Density Functional Theory calculations with the Boltzmann transport method. The calculated values of the Seebeck coefficient are within the typical magnitude (200 - 600 $μV/K$) indicative of a good thermoelectric material. The tobermorite models are predicted to be intrinsically $p$-type thermoelectric material because of the presence of large concentration of the Si-O tetrahedra sites. The calculated electronic $ZT$ for the tobermorite models have their optimal values of 0.983 at (400 $\mathrm{K}$ and $10^{17}$ $\mathrm{cm^{-3}}$) for tobermorite 9 Å, 0.985 at (400 $\mathrm{K}$ and $10^{17}$ $\mathrm{cm^{-3}}$) for tobermorite 11 Å and 1.20 at (225 $\mathrm{K}$ and $10^{19}$ $\mathrm{cm^{-3}}$) for tobermorite 14 Å, respectively. △ Less

Submitted 30 November, 2022; originally announced November 2022.

Comments: 34 pages, 10 figures

Universal singularities of anomalous diffusion in the Richardson class

Authors: Attilio L. Stella, Aleksei Chechkin, Gianluca Teza

Abstract: Inhomogeneous environments are rather ubiquitous in nature, often implying anomalies resulting in deviation from Gaussianity of diffusion processes. While sub- and superdiffusion are usually due to conversing environmental features (hindering or favoring the motion, respectively), they are both observed in systems ranging from the micro- to the cosmological scale. Here we show how a model encompas… ▽ More Inhomogeneous environments are rather ubiquitous in nature, often implying anomalies resulting in deviation from Gaussianity of diffusion processes. While sub- and superdiffusion are usually due to conversing environmental features (hindering or favoring the motion, respectively), they are both observed in systems ranging from the micro- to the cosmological scale. Here we show how a model encompassing sub- and superdiffusion in an inhomogeneous environment exhibits a critical singularity in the normalized generator of the cumulants. The singularity originates directly from the asymptotics of the non-Gaussian scaling function of displacement, which we prove to be independent of other details and hence to retain a universal character. Our analysis, based on the method first applied in [A. L. Stella et al., arXiv:2209.02042 (2022)], further allows to establish a relation between the asympototics and diffusion exponents characteristic of processes in the Richardson class. Extensive numerical tests fully confirm the results. △ Less

Submitted 27 November, 2022; originally announced November 2022.

Investigating the origin of optical and X-ray pulsations of the transitional millisecond pulsar PSR J1023+0038

Authors: G. Illiano, A. Papitto, F. Ambrosino, A. Miraval Zanon, F. Coti Zelati, L. Stella, L. Zampieri, A. Burtovoi, S. Campana, P. Casella, M. Cecconi, D. de Martino, M. Fiori, A. Ghedina, M. Gonzales, M. Hernandez Diaz, G. L. Israel, F. Leone, G. Naletto, H. Perez Ventura, C. Riverol, L. Riverol, D. F. Torres, M. Turchetta

Abstract: PSR J1023+0038 is the first millisecond pulsar that was ever observed as an optical and UV pulsar. So far, it is the only optical transitional millisecond pulsar. The rotation- and accretion-powered emission mechanisms hardly individually explain the observed characteristics of optical pulsations. A synergistic model, combining these standard emission processes, was proposed to explain the origin… ▽ More PSR J1023+0038 is the first millisecond pulsar that was ever observed as an optical and UV pulsar. So far, it is the only optical transitional millisecond pulsar. The rotation- and accretion-powered emission mechanisms hardly individually explain the observed characteristics of optical pulsations. A synergistic model, combining these standard emission processes, was proposed to explain the origin of the X-ray/UV/optical pulsations. We study the phase lag between the pulses in the optical and X-ray bands to gain insight into the physical mechanisms that cause it. We performed a detailed timing analysis of simultaneous or quasi-simultaneous observations in the X-ray band, acquired with the XMM-Newton and NICER satellites, and in the optical band, with the fast photometers SiFAP2 (mounted at the 3.6 m Telescopio Nazionale Galileo) and Aqueye+ (mounted at the 1.8 m Copernicus Telescope). We estimated the time lag of the optical pulsation with respect to that in the X-rays by modeling the folded pulse profiles with two harmonic components. Optical pulses lag the X-ray pulses by $\sim$ 150 $μ$s in observations acquired with instruments (NICER and Aqueye+) whose absolute timing uncertainty is much smaller than the measured lag. We also show that the phase lag between optical and X-ray pulsations lies in a limited range of values, $δφ\in$ (0 $-$ 0.15), which is maintained over timescales of about five years. This indicates that both pulsations originate from the same region, and it supports the hypothesis of a common emission mechanism. Our results are interpreted in the shock-driven mini pulsar nebula scenario. This scenario suggests that optical and X-ray pulses are produced by synchrotron emission from the shock that formed within a few light cylinder radii away ($\sim$ 100 km) from the pulsar, where its striped wind encounters the accretion disk inflow. △ Less

Submitted 23 November, 2022; originally announced November 2022.

Comments: 13 pages, 6 figures, 5 tables, accepted for publication in A&A

Journal ref: A&A 669, A26 (2023)

Black hole mass and spin measurements through the Relativistic Precession Model: XTE J1859+226

Authors: S. E Motta, T. Belloni, L. Stella, G. Pappas, J. A. Casares, T. Muñoz-Darias, M. A. P. Torres, I. V. Yanes-Rizo

Abstract: The X-ray light curves of accreting black holes and neutron stars in binary systems show various types of quasi-periodic oscillations (QPOs), the origin of which is still debated. The Relativistic Precession Model identifies the QPO frequencies with fundamental time scales from General Relativity, and has been proposed as a possible explanation of certain types of such oscillations. Under specific… ▽ More The X-ray light curves of accreting black holes and neutron stars in binary systems show various types of quasi-periodic oscillations (QPOs), the origin of which is still debated. The Relativistic Precession Model identifies the QPO frequencies with fundamental time scales from General Relativity, and has been proposed as a possible explanation of certain types of such oscillations. Under specific conditions (i.e., the detection of a particular QPOs triplet) such a model can be used to obtain self-consistent measurements of the mass and spin of the compact object. So far this has been possible only in the black hole binary GRO J1655-40. In the RXTE/PCA data from the 1999-2000 outburst of the black hole transient XTE J1859+226 we found a QPO triplet, and used the the Relativistic Precession Model to obtain high-precision measurements of the black hole mass and spin - M = (7.85+/-0.46) Msun, a* = 0.149+/-0.005 - the former being consistent with the most recent dynamical mass determination from optical measurements. Similarly to what has been already observed in other black hole systems, the frequencies of the QPOs and broad-band noise components match the general relativistic frequencies of particle motion close to the compact object predicted by the model. Our findings confirm previous results and further support the validity of the Relativistic Precession Model, which is the only electromagnetic-measurement-based method that so far has consistently yielded spins close to those from the gravitational waves produced by merging binary black holes. △ Less

Submitted 21 September, 2022; originally announced September 2022.

Comments: 8 pages, 5 figures, accepted for publication in MNRAs

Anomalous dynamical scaling determines universal critical singularities

Authors: Attilio L. Stella, Aleksei Chechkin, Gianluca Teza

Abstract: Anomalous diffusion phenomena occur on length scales spanning from intracellular to astrophysical ranges. A specific form of decay at large argument of the probability density function of rescaled displacement (scaling function) is derived and shown to imply universal singularities in the normalized cumulant generator. Exact calculations for continuous time random walks provide paradigmatic exampl… ▽ More Anomalous diffusion phenomena occur on length scales spanning from intracellular to astrophysical ranges. A specific form of decay at large argument of the probability density function of rescaled displacement (scaling function) is derived and shown to imply universal singularities in the normalized cumulant generator. Exact calculations for continuous time random walks provide paradigmatic examples connected with singularities of second order phase transitions. In the biased case scaling is restricted to displacements in the drift direction and singularities have no equilibrium analogue. △ Less

Submitted 24 October, 2022; v1 submitted 5 September, 2022; originally announced September 2022.

UV and X-ray pulse amplitude variability in the transitional millisecond pulsar PSR J1023+0038

Authors: A. Miraval Zanon, F. Ambrosino, F. Coti Zelati, S. Campana, A. Papitto, G. Illiano, G. L. Israel, L. Stella, P. D'Avanzo, M. C. Baglio

Abstract: The transitional millisecond pulsar PSR\,J1023+0038 is the first millisecond pulsar discovered to emit UV and optical pulses. Here we present the results of the UV and X-ray phase-resolved timing analysis of observations performed with the Hubble Space Telescope, \textit{XMM-Newton} and NuSTAR satellites between 2014 and 2021. Ultraviolet pulsations are detected in the high luminosity mode and dis… ▽ More The transitional millisecond pulsar PSR\,J1023+0038 is the first millisecond pulsar discovered to emit UV and optical pulses. Here we present the results of the UV and X-ray phase-resolved timing analysis of observations performed with the Hubble Space Telescope, \textit{XMM-Newton} and NuSTAR satellites between 2014 and 2021. Ultraviolet pulsations are detected in the high luminosity mode and disappear during low and flaring modes, similar to what is observed in the X-ray band. In the high mode, we find variability in both the UV and X-ray pulse amplitudes. The root mean square pulsed amplitude in the UV band ranges from $\sim$2.1\% down to $\sim$0.7\%, while it oscillates in the interval $5.5-12\%$ in the X-ray band. This variability is not correlated with the orbital phase, like what has been observed in the optical band. Notwithstanding the rather low statistics, we have marginal evidence that variations in the pulse amplitude do not occur simultaneously in the UV and X-ray bands. When the UV pulsed amplitude decreases below the detection threshold, no significant variation in the X-ray pulsed amplitude is observed. These oscillations in the pulse amplitude could be caused by small random variations in the mass accretion rate leading to a variation in the size of the intra-binary shock region. Finally, we find that the pulsed flux spectral distribution from the X-ray to the UV band is well fitted using a power-law relation of the form $νF_ν^{pulsed} \sim ν^{0.4}$. This supports the hypothesis of a common physical mechanism underlying the X-ray, UV, and optical pulsed emissions in PSR\,J1023+0038. △ Less

Submitted 4 March, 2022; originally announced March 2022.

Comments: 8 pages, 4 figures, accepted for publication in A&A

Journal ref: A&A 660, A63 (2022)

The supernova of the MAGIC GRB190114C

Authors: A. Melandri, L. Izzo, E. Pian, D. B. Malesani, M. Della Valle, A. Rossi, P. D'Avanzo, D. Guetta, P. A. Mazzali, S. Benetti, N. Masetti, E. Palazzi, S. Savaglio, L. Amati, L. A. Antonelli, C. Ashall, M. G. Bernardini, S. Campana, R. Carini, S. Covino, V. D'Elia, A. de Ugarte Postigo, M. De Pasquale, A. V. Filippenko, A. S. Fruchter , et al. (20 additional authors not shown)

Abstract: We observed GRB190114C (redshift z = 0.4245), the first GRB ever detected at TeV energies, at optical and near-infrared wavelengths with several ground-based telescopes and the Hubble Space Telescope, with the primary goal of studying its underlying supernova, SN2019jrj. The monitoring spanned the time interval between 1.3 and 370 days after the burst, in the observer frame. We find that the after… ▽ More We observed GRB190114C (redshift z = 0.4245), the first GRB ever detected at TeV energies, at optical and near-infrared wavelengths with several ground-based telescopes and the Hubble Space Telescope, with the primary goal of studying its underlying supernova, SN2019jrj. The monitoring spanned the time interval between 1.3 and 370 days after the burst, in the observer frame. We find that the afterglow emission can be modelled with a forward shock propagating in a uniform medium modified by time-variable extinction along the line of sight. A jet break could be present after 7 rest-frame days, and accordingly the maximum luminosity of the underlying SN ranges between that of stripped-envelope corecollapse supernovae (SNe) of intermediate luminosity, and that of the luminous GRB-associated SN2013dx. The observed spectral absorption lines of SN2019jrj are not as broad as in classical GRB-SNe, and are rather more similar to those of less-luminous core-collapse SNe. Taking the broad-lined stripped-envelope core-collapse SN2004aw as an analogue, we tentatively derive the basic physical properties of SN2019jrj. We discuss the possibility that a fraction of the TeV emission of this source might have had a hadronic origin and estimate the expected high-energy neutrino detection level with IceCube. △ Less

Submitted 9 December, 2021; originally announced December 2021.

Comments: 9 pages, 4 figures, accepted for publication in A&A

Journal ref: A&A 659, A39 (2022)

Digital Twinning Remote Laboratories for Online Practical Learning

Authors: Claire Palmer, Ben Roullier, Muhammad Aamir, Frank McQuade, Leonardo Stella, Ashiq Anjum

Abstract: The COVID19 pandemic has demonstrated a need for remote learning and virtual learning applications such as virtual reality (VR) and tablet-based solutions. Creating complex learning scenarios by developers is highly time-consuming and can take over a year. It is also costly to employ teams of system analysts, developers and 3D artists. There is a requirement to provide a simple method to enable le… ▽ More The COVID19 pandemic has demonstrated a need for remote learning and virtual learning applications such as virtual reality (VR) and tablet-based solutions. Creating complex learning scenarios by developers is highly time-consuming and can take over a year. It is also costly to employ teams of system analysts, developers and 3D artists. There is a requirement to provide a simple method to enable lecturers to create their own content for their laboratory tutorials. Research has been undertaken into developing generic models to enable the semi-automatic creation of a virtual learning tools for subjects that require practical interactions with the lab resources. In addition to the system for creating digital twins, a case study describing the creation of a virtual learning application for an electrical laboratory tutorial has been presented. △ Less

Submitted 21 July, 2022; v1 submitted 1 December, 2021; originally announced December 2021.

Comments: 56 pages, 14 figures. arXiv admin note: text overlap with arXiv:2106.09344. Added version accepted for publication

Journal ref: Production and Manufacturing Research, 20(1), 519-545 (2022)

Virtual Reality based Digital Twin System for remote laboratories and online practical learning

Authors: Claire Palmer, Ben Roullier, Muhammad Aamir, Leonardo Stella, Uchenna Diala, Ashiq Anjum, Frank Mcquade, Keith Cox, Alex Calvert

Abstract: There is a need for remote learning and virtual learning applications such as virtual reality (VR) and tablet-based solutions which the current pandemic has demonstrated. Creating complex learning scenarios by developers is highly time-consuming and can take over a year. There is a need to provide a simple method to enable lecturers to create their own content for their laboratory tutorials. Resea… ▽ More There is a need for remote learning and virtual learning applications such as virtual reality (VR) and tablet-based solutions which the current pandemic has demonstrated. Creating complex learning scenarios by developers is highly time-consuming and can take over a year. There is a need to provide a simple method to enable lecturers to create their own content for their laboratory tutorials. Research is currently being undertaken into developing generic models to enable the semi-automatic creation of a virtual learning application. A case study describing the creation of a virtual learning application for an electrical laboratory tutorial is presented. △ Less

Submitted 17 June, 2021; originally announced June 2021.

Comments: 6 pages, 4 figures, accepted for publication ICMR2021 18th International Conference in Manufacturing Research Virtual Conference hosted by the University of Derby, UK 7 - 10 September 2021

The Impact of Irrational Behaviours in the Optional Prisoner's Dilemma with Game-Environment Feedback

Authors: Leonardo Stella, Dario Bauso

Abstract: In the optional prisoner's dilemma (OPD), players can choose to cooperate and defect as usual, but can also abstain as a third possible strategy. This strategy models the players' participation in the game and is a relevant aspect in many settings, e.g. social networks or opinion dynamics where abstention is an option during an election. In this paper, we provide a formulation of the OPD where we… ▽ More In the optional prisoner's dilemma (OPD), players can choose to cooperate and defect as usual, but can also abstain as a third possible strategy. This strategy models the players' participation in the game and is a relevant aspect in many settings, e.g. social networks or opinion dynamics where abstention is an option during an election. In this paper, we provide a formulation of the OPD where we consider irrational behaviours in the population inspired by prospect theory. Prospect theory has gained increasing popularity in recent times thanks to its ability to capture aspects such as reference dependence or loss aversion which are common in human behaviour. This element is original in our formulation of the game and is incorporated in our framework through pairwise comparison dynamics. Recently, the impact of the environment has been studied in the form of feedback on the population dynamics. Another element of novelty in our work is the extension of the game-environment feedback to the OPD in two forms of dynamics, the replicator and the pairwise comparison. The contribution of this paper is threefold. First, we propose a modelling framework where prospect theory is used to capture irrational behaviours in an evolutionary game with game-environment feedback. Second, we carry out the stability analysis of the system equilibria and discuss the oscillating behaviours arising from the game-environment feedback. Finally, we extend our previous results to the OPD and we discuss the main differences between the model resulting from the replicator dynamics and the one resulting from the pairwise comparison dynamics. △ Less

Submitted 13 June, 2021; originally announced June 2021.

Comments: 12 pages, 5 figures

Entropic measure unveils country competitiveness and product specialization in the World trade web

Authors: Gianluca Teza, Michele Caraglio, Attilio L. Stella

Abstract: We show how the Shannon entropy function can be used as a basis to set up complexity measures weighting the economic efficiency of countries and the specialization of products beyond bare diversification. This entropy function guarantees the existence of a fixed point which is rapidly reached by an iterative scheme converging to our self-consistent measures. Our approach naturally allows to decomp… ▽ More We show how the Shannon entropy function can be used as a basis to set up complexity measures weighting the economic efficiency of countries and the specialization of products beyond bare diversification. This entropy function guarantees the existence of a fixed point which is rapidly reached by an iterative scheme converging to our self-consistent measures. Our approach naturally allows to decompose into inter-sectorial and intra-sectorial contributions the country competitivity measure if products are partitioned into larger categories. Besides outlining the technical features and advantages of the method, we describe a wide range of results arising from the analysis of the obtained rankings and we benchmark these observations against those established with other economical parameters. These comparisons allow to partition countries and products into various main typologies, with well-revealed characterizing features. Our methods have wide applicability to general problems of ranking in bipartite networks. △ Less

Submitted 3 June, 2021; originally announced June 2021.

Journal ref: Sci Rep 11, 10189 (2021)

Exploring higher order images with Fe K_α-lines from relativistic disks: black hole spin determination and bias

Authors: M. Falanga, P. Bakala, R. La Placa, V. De Falco, A. De Rosa, L. Stella

Abstract: We study the contributions to the relativistic Fe $K_α$ line profile from higher order images (HOIs) produced by strongly deflected rays from the disk which cross the plunging region, located between the innermost stable circular orbit (ISCO) radius and the event horizon of a Kerr black hole. We investigate the characteristics features imprinted by the HOIs in the line profile for different black… ▽ More We study the contributions to the relativistic Fe $K_α$ line profile from higher order images (HOIs) produced by strongly deflected rays from the disk which cross the plunging region, located between the innermost stable circular orbit (ISCO) radius and the event horizon of a Kerr black hole. We investigate the characteristics features imprinted by the HOIs in the line profile for different black hole spins, disk emissivity laws and inclinations. We find that they extend from the red wing of the profile up to energies slightly lower than those of the blue peak, adding $\sim 0.4 - 1.3$\% to the total line flux. The contribution to the specific flux is often in the $\sim 1$\% to 7\% range, with the highest values attained for low and negative spin ($a\lesssim 0.3$) black holes surrounded by intermediate inclination angle ($i\sim40^{\circ}$) disks. We simulate future observations of a black hole X-ray binary system with the Large Area Detector of the planned X-ray astronomy \emph{enhanced X-ray Timing and Polarimetry Mission} (eXTP) and find that the \fekal\ of systems accreting at $\lesssim 1 $\% the Eddington rate are affected by the HOI features for a range of parameters. This would provide evidence of the extreme gravitational lensing of HOI rays. Our simulations show also that not accounting for HOI contributions to the Fe $K_α$ line profile may systematically bias measurements of the black hole spin parameter towards values higher by up to $\sim 0.3$ than the inputted ones. △ Less

Submitted 15 April, 2021; originally announced April 2021.

Comments: accepted for publication in MNRAS, 11 pages, 7 figures, 2 tables

Millisecond Magnetars

Authors: Simone Dall'Osso, Luigi Stella

Abstract: Two classes of X-ray/$γ$-ray sources, the Soft Gamma Repeaters and the Anomalous X-ray Pulsars have been identified with isolated, slowly spinning magnetars, neutron stars whose emission draws energy from their extremely strong magnetic field ($\sim 10^{15}-10^{16}$ G). Magnetars are believed to form with millisecond spin period and to represent an important fraction of the whole population of you… ▽ More Two classes of X-ray/$γ$-ray sources, the Soft Gamma Repeaters and the Anomalous X-ray Pulsars have been identified with isolated, slowly spinning magnetars, neutron stars whose emission draws energy from their extremely strong magnetic field ($\sim 10^{15}-10^{16}$ G). Magnetars are believed to form with millisecond spin period and to represent an important fraction of the whole population of young neutron stars. Newborn magnetars can convert very quickly their rotational energy into electromagnetic and/or gravitational waves, by virtue of their strong magnetic fields and fast spins. This chapter provides a brief summary of astrophysical problems and scenarios in which millisecond magnetars are believed to play a key role: these include Gamma Ray Bursts, Supernovae, Gravitational Wave events and Fast Radio Bursts. △ Less

Submitted 19 March, 2021; originally announced March 2021.

Comments: 31 pages, 3 figures, Preprint of 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

Optical and ultraviolet pulsed emission from an accreting millisecond pulsar

Authors: F. Ambrosino, A. Miraval Zanon, A. Papitto, F. Coti Zelati, S. Campana, P. D'Avanzo, L. Stella, T. Di Salvo, L. Burderi, P. Casella, A. Sanna, D. de Martino, M. Cadelano, A. Ghedina, F. Leone, F. Meddi, P. Cretaro, M. C. Baglio, E. Poretti, R. P. Mignani, D. F. Torres, G. L. Israel, M. Cecconi, D. M. Russell, M. D. Gonzalez Gomez , et al. (6 additional authors not shown)

Abstract: Millisecond spinning, low magnetic field neutron stars are believed to attain their fast rotation in a 0.1-1 Gyr-long phase during which they accrete matter endowed with angular momentum from a low-mass companion star. Despite extensive searches, coherent periodicities originating from accreting neutron star magnetospheres have been detected only at X-ray energies and in ~10% of the presently know… ▽ More Millisecond spinning, low magnetic field neutron stars are believed to attain their fast rotation in a 0.1-1 Gyr-long phase during which they accrete matter endowed with angular momentum from a low-mass companion star. Despite extensive searches, coherent periodicities originating from accreting neutron star magnetospheres have been detected only at X-ray energies and in ~10% of the presently known systems. Here we report the detection of optical and ultraviolet coherent pulsations at the X-ray period of the transient low mass X-ray binary system SAX J1808.4-3658, during an accretion outburst that occurred in August 2019. At the time of the observations, the pulsar was surrounded by an accretion disc, displayed X-ray pulsations and its luminosity was consistent with magnetically funneled accretion onto the neutron star. Current accretion models fail to account for the luminosity of both optical and ultraviolet pulsations; these are instead more likely driven by synchro-curvature radiation in the pulsar magnetosphere or just outside of it. This interpretation would imply that particle acceleration can take place even when mass accretion is going on, and opens up new perspectives in the study of coherent optical/UV pulsations from fast spinning accreting neutron stars in low-mass X-ray binary systems. △ Less

Submitted 23 February, 2021; originally announced February 2021.

Comments: 47 pages, 9 figures. The first two authors contributed equally to this work; Nature Astronomy (2021), published on-line on February 22, 2021; doi:10.1038/s41550-021-01308-0

Journal ref: Nature Astronomy (2021)

X-ray and Radio Bursts from the Magnetar 1E1547.0-5408

Authors: G. L. Israel, M. Burgay, N. Rea, P. Esposito, A. Possenti, S. Dall'Osso, L. Stella, M. Pilia, A. Tiengo, A. Ridnaia, A. Y. Lien, D. D. Frederiks, F. Bernardini

Abstract: We report on simultaneous radio and X-ray observations of the radio-emitting magnetar 1E1547.0-5408 on 2009 January 25 and February 3, with the 64-m Parkes radio telescope and the Chandra and XMM-Newton X-ray observatories. The magnetar was observed in a period of intense X-ray bursting activity and enhanced X-ray emission. We report here on the detection of two radio bursts from 1E1547.0-5408, re… ▽ More We report on simultaneous radio and X-ray observations of the radio-emitting magnetar 1E1547.0-5408 on 2009 January 25 and February 3, with the 64-m Parkes radio telescope and the Chandra and XMM-Newton X-ray observatories. The magnetar was observed in a period of intense X-ray bursting activity and enhanced X-ray emission. We report here on the detection of two radio bursts from 1E1547.0-5408, reminiscent of Fast Radio Bursts (FRBs). One of the radio bursts was anticipated by ~1s (about half a rotation period of the pulsar) by a bright SGR-like X-ray burst, resulting in a F_radio/F_X ~ 10^-9. Radio pulsations were not detected during the observation showing the FRB-like radio bursts, while they were detected in the previous radio observation. We also found that the two radio bursts are neither aligned with the latter radio pulsations nor with the peak of the X-ray pulse profile (phase shift of ~0.2). Comparing the luminosity of these FRB-like bursts and those reported from SGR1935+2154, we find that the wide range in radio efficiency and/or luminosity of magnetar bursts in the Galaxy may bridge the gap between "ordinary" pulsar radio bursts and the extragalactic FRB phenomenon. △ Less

Submitted 12 November, 2020; originally announced November 2020.

Comments: Accepted for publication in ApJ. 9 pages, 3 figures

Topological disentanglement of linear polymers under tension

Authors: Michele Caraglio, Boris Marcone, Fulvio Baldovin, Enzo Orlandini, Attilio L. Stella

Abstract: We develop a theoretical description of the topological disentanglement occurring when torus knots reach the ends of a semi-flexible polymer under tension. These include decays into simpler knots and total unknotting. The minimal number of crossings and the minimal knot contour length are the topological invariants playing a key role in the model. The crossings behave as particles diffusing along… ▽ More We develop a theoretical description of the topological disentanglement occurring when torus knots reach the ends of a semi-flexible polymer under tension. These include decays into simpler knots and total unknotting. The minimal number of crossings and the minimal knot contour length are the topological invariants playing a key role in the model. The crossings behave as particles diffusing along the chain and the application of appropriate boundary conditions at the ends of the chain accounts for the knot disentanglement. Starting from the number of particles and their positions, suitable rules allow reconstructing the type and location of the knot moving on the chain. Our theory is extensively benchmarked with corresponding Molecular Dynamics simulations and the results show a remarkable agreement between the simulations and the theoretical predictions of the model. △ Less

Submitted 3 November, 2020; v1 submitted 20 October, 2020; originally announced October 2020.

Comments: 12 pages, 11 figures

Journal ref: Polymers (2020) 12 (11), 2580

The Role of Asymptomatic Individuals in the COVID-19 Pandemic via Complex Networks

Authors: Leonardo Stella, Alejandro Pinel Martínez, Dario Bauso, Patrizio Colaneri

Abstract: Recent seroprevalence studies have tried to estimate the real number of asymptomatic cases affected by COVID-19. It is of paramount importance to understand the impact of these infections in order to prevent a second wave. This study aims to model the interactions in the population by means of a complex network and to shed some light on the effectiveness of localised control measures in Italy in r… ▽ More Recent seroprevalence studies have tried to estimate the real number of asymptomatic cases affected by COVID-19. It is of paramount importance to understand the impact of these infections in order to prevent a second wave. This study aims to model the interactions in the population by means of a complex network and to shed some light on the effectiveness of localised control measures in Italy in relation to the school opening in mid-September. The formulation of an epidemiological predictive model is given: the advantage of using this model lies in that it discriminates between asymptomatic and symptomatic cases of COVID-19 as the interactions with these two categories of infected individuals are captured separately, allowing for a study on the impact of asymptomatic cases. This model is then extended to a structured nonhomogeneous version by means of the Watts-Strogatz complex network, which is adopted widely to model societal interactions as it holds the small world property. Finally, a case study on the situation in Italy is given: first the homogeneous model is used to compare the official data with the data of the recent seroprevalence study from Istat; second, in view of the return to school in mid-September, a study at regional level is conducted. The results of this study highlight the importance of coordinating the deployment of appropriate control measures that take into account the role of asymptomatic infections, especially in younger individuals, and inter-regional connectivity in Italy. △ Less

Submitted 8 September, 2020; originally announced September 2020.

Comments: 38 pages, journal

Anomaly Detection at Scale: The Case for Deep Distributional Time Series Models

Authors: Fadhel Ayed, Lorenzo Stella, Tim Januschowski, Jan Gasthaus

Abstract: This paper introduces a new methodology for detecting anomalies in time series data, with a primary application to monitoring the health of (micro-) services and cloud resources. The main novelty in our approach is that instead of modeling time series consisting of real values or vectors of real values, we model time series of probability distributions over real values (or vectors). This extension… ▽ More This paper introduces a new methodology for detecting anomalies in time series data, with a primary application to monitoring the health of (micro-) services and cloud resources. The main novelty in our approach is that instead of modeling time series consisting of real values or vectors of real values, we model time series of probability distributions over real values (or vectors). This extension to time series of probability distributions allows the technique to be applied to the common scenario where the data is generated by requests coming in to a service, which is then aggregated at a fixed temporal frequency. Our method is amenable to streaming anomaly detection and scales to monitoring for anomalies on millions of time series. We show the superior accuracy of our method on synthetic and public real-world data. On the Yahoo Webscope data set, we outperform the state of the art in 3 out of 4 data sets and we show that we outperform popular open-source anomaly detection tools by up to 17% average improvement for a real-world data set. △ Less

Submitted 30 July, 2020; originally announced July 2020.

Douglas-Rachford splitting and ADMM for nonconvex optimization: Accelerated and Newton-type linesearch algorithms

Authors: Andreas Themelis, Lorenzo Stella, Panagiotis Patrinos

Abstract: Although the performance of popular optimization algorithms such as Douglas-Rachford splitting (DRS) and the ADMM is satisfactory in small and well-scaled problems, ill conditioning and problem size pose a severe obstacle to their reliable employment. Expanding on recent convergence results for DRS and ADMM applied to nonconvex problems, we propose two linesearch algorithms to enhance and robustif… ▽ More Although the performance of popular optimization algorithms such as Douglas-Rachford splitting (DRS) and the ADMM is satisfactory in small and well-scaled problems, ill conditioning and problem size pose a severe obstacle to their reliable employment. Expanding on recent convergence results for DRS and ADMM applied to nonconvex problems, we propose two linesearch algorithms to enhance and robustify these methods by means of quasi-Newton directions. The proposed algorithms are suited for nonconvex problems, require the same black-box oracle of DRS and ADMM, and maintain their (subsequential) convergence properties. Numerical evidence shows that the employment of L-BFGS in the proposed framework greatly improves convergence of DRS and ADMM, making them robust to ill conditioning. Under regularity and nondegeneracy assumptions at the limit point, superlinear convergence is shown when quasi-Newton Broyden directions are adopted. △ Less

Submitted 3 November, 2021; v1 submitted 20 May, 2020; originally announced May 2020.

MSC Class: 90C06; 90C25; 90C26; 49J52; 49J53

Journal ref: Comput Optim Appl 82, 395-440 (2022)

Deep Learning for Time Series Forecasting: Tutorial and Literature Survey

Authors: Konstantinos Benidis, Syama Sundar Rangapuram, Valentin Flunkert, Yuyang Wang, Danielle Maddix, Caner Turkmen, Jan Gasthaus, Michael Bohlke-Schneider, David Salinas, Lorenzo Stella, Francois-Xavier Aubet, Laurent Callot, Tim Januschowski

Abstract: Deep learning based forecasting methods have become the methods of choice in many applications of time series prediction or forecasting often outperforming other approaches. Consequently, over the last years, these methods are now ubiquitous in large-scale industrial forecasting applications and have consistently ranked among the best entries in forecasting competitions (e.g., M4 and M5). This pra… ▽ More Deep learning based forecasting methods have become the methods of choice in many applications of time series prediction or forecasting often outperforming other approaches. Consequently, over the last years, these methods are now ubiquitous in large-scale industrial forecasting applications and have consistently ranked among the best entries in forecasting competitions (e.g., M4 and M5). This practical success has further increased the academic interest to understand and improve deep forecasting methods. In this article we provide an introduction and overview of the field: We present important building blocks for deep forecasting in some depth; using these building blocks, we then survey the breadth of the recent deep forecasting literature. △ Less

Submitted 15 June, 2022; v1 submitted 21 April, 2020; originally announced April 2020.

Comments: 33 pages, 6 figures

ACM Class: A.1

Journal ref: ACM Computing Surveys (2022)

A very young radio-loud magnetar

Authors: P. Esposito, N. Rea, A. Borghese, F. Coti Zelati, D. Viganò, G. L. Israel, A. Tiengo, A. Ridolfi, A. Possenti, M. Burgay, D. Götz, F. Pintore, L. Stella, C. Dehman, M. Ronchi, S. Campana, A. Garcia-Garcia, V. Graber, S. Mereghetti, R. Perna, G. A. Rodríguez Castillo, R. Turolla, S. Zane

Abstract: The magnetar Swift ,J1818.0-1607 was discovered in March 2020 when Swift detected a 9 ms hard X-ray burst and a long-lived outburst. Prompt X-ray observations revealed a spin period of 1.36 s, soon confirmed by the discovery of radio pulsations. We report here on the analysis of the Swift burst and follow-up X-ray and radio observations. The burst average luminosity was… ▽ More The magnetar Swift ,J1818.0-1607 was discovered in March 2020 when Swift detected a 9 ms hard X-ray burst and a long-lived outburst. Prompt X-ray observations revealed a spin period of 1.36 s, soon confirmed by the discovery of radio pulsations. We report here on the analysis of the Swift burst and follow-up X-ray and radio observations. The burst average luminosity was $L_{\rm burst} \sim2\times 10^{39}$ erg/s (at 4.8 kpc). Simultaneous observations with XMM-Newton and NuSTAR three days after the burst provided a source spectrum well fit by an absorbed blackbody ($N_{\rm H} = (1.13\pm0.03) \times 10^{23}$ cm$^{-2}$ and $kT = 1.16\pm0.03$ keV) plus a power-law ($Γ=0.0\pm1.3$) in the 1-20 keV band, with a luminosity of $\sim$$8\times10^{34}$ erg/s, dominated by the blackbody emission. From our timing analysis, we derive a dipolar magnetic field $B \sim 7\times10^{14}$ G, spin-down luminosity $\dot{E}_{\rm rot} \sim 1.4\times10^{36}$ erg/s and characteristic age of 240 yr, the shortest currently known. Archival observations led to an upper limit on the quiescent luminosity $<$$5.5\times10^{33}$ erg/s, lower than the value expected from magnetar cooling models at the source characteristic age. A 1 hr radio observation with the Sardinia Radio Telescope taken about 1 week after the X-ray burst detected a number of strong and short radio pulses at 1.5 GHz, in addition to regular pulsed emission; they were emitted at an average rate 0.9 min$^{-1}$ and accounted for $\sim$50% of the total pulsed radio fluence. We conclude that Swift ,J1818.0-1607 is a peculiar magnetar belonging to the small, diverse group of young neutron stars with properties straddling those of rotationally and magnetically powered pulsars. Future observations will make a better estimation of the age possible by measuring the spin-down rate in quiescence. △ Less

Submitted 25 May, 2020; v1 submitted 8 April, 2020; originally announced April 2020.

Comments: 10 pages, 4 figures, 1 table; revised version accepted for publication in The Astrophysical Journal Letters

Journal ref: The Astrophysical Journal Letters, Volume 896, L30 (2020)

Exact Coarse Graining Preserves Entropy Production out of Equilibrium

Authors: Gianluca Teza, Attilio L. Stella

Abstract: The entropy production rate associated with broken time-reversal symmetry provides an essential characterization of nanosystems out of equilibrium, from driven colloidal particles to molecular motors. Limited access to the dynamical states is generally expected to hinder the correct estimation of this observable. Here we show how memoryless jump processes can be coarse grained exactly preserving i… ▽ More The entropy production rate associated with broken time-reversal symmetry provides an essential characterization of nanosystems out of equilibrium, from driven colloidal particles to molecular motors. Limited access to the dynamical states is generally expected to hinder the correct estimation of this observable. Here we show how memoryless jump processes can be coarse grained exactly preserving its average and fluctuations at stationarity. This supports univocal applicability of fluctuation theorems for entropy and allows inference of the genuine thermodynamics together with inaccessible process details. △ Less

Submitted 24 November, 2020; v1 submitted 19 March, 2020; originally announced March 2020.

Journal ref: Phys. Rev. Lett. 125, 110601 (2020)

Neutron Star Radius-to-mass Ratio from Partial Accretion Disc Occultation as Measured through Fe K$α$ Line Profiles

Authors: Riccardo La Placa, Luigi Stella, Alessando Papitto, Pavel Bakala, Tiziana Di Salvo, Maurizio Falanga, Vittorio De Falco, Alessandra De Rosa

Abstract: We present a new method to measure the radius-to-mass ratio (R/M) of weakly magnetic, disc-accreting neutron stars by exploiting the occultation of parts of the inner disc by the star itself. This occultation imprints characteristic features on the X-ray line profile that are unique and are expected to be present in low mass X-ray binary systems seen under inclinations higher than ~65 degrees. We… ▽ More We present a new method to measure the radius-to-mass ratio (R/M) of weakly magnetic, disc-accreting neutron stars by exploiting the occultation of parts of the inner disc by the star itself. This occultation imprints characteristic features on the X-ray line profile that are unique and are expected to be present in low mass X-ray binary systems seen under inclinations higher than ~65 degrees. We analyse a NuSTAR observation of a good candidate system, 4U 1636-53, and find that X-ray spectra from current instrumentation are unlikely to single out the occultation features owing to insufficient signal-to-noise. Based on an extensive set of simulations we show that large-area X-ray detectors of the future generation could measure R/M to ~2÷3% precision over a range of inclinations. Such is the precision in radius determination required to derive tight constraints on the equation of state of ultradense matter and it represents the goal that other methods too aim to achieve in the future. △ Less

Submitted 17 March, 2020; originally announced March 2020.

Comments: 17 pages, 8 figures; this is a pre-print edition of an article that has been accepted for publication in the Astrophysical Journal

Journal ref: The Astrophysical Journal 893 (2020) 129

The X-ray outburst of the Galactic Center magnetar over six years of Chandra observations

Authors: N. Rea, F. Coti Zelati, D. Vigano', A. Papitto, F. Baganoff, A. Borghese, S. Campana, P. Esposito, D. Haggard, G. L. Israel, S. Mereghetti, R. Mignani, R. Perna, J. A. Pons, G. Ponti, L. Stella, D. F. Torres, R. Turolla, S. Zane

Abstract: The magnetar SGR J1745-2900 discovered at parsecs distance from the Milky Way central black hole, Sagittarius A*, represents the closest pulsar to a supermassive black hole ever detected. Furthermore, its intriguing radio emission has been used to study the environment of the black hole, as well as to derive a precise position and proper motion for this object. The discovery of SGR J1745-2900 has… ▽ More The magnetar SGR J1745-2900 discovered at parsecs distance from the Milky Way central black hole, Sagittarius A*, represents the closest pulsar to a supermassive black hole ever detected. Furthermore, its intriguing radio emission has been used to study the environment of the black hole, as well as to derive a precise position and proper motion for this object. The discovery of SGR J1745-2900 has opened interesting debates about the number, age and nature of pulsars expected in the Galactic center region. In this work, we present extensive X-ray monitoring of the outburst of SGR J1745-2900 using the Chandra X-ray Observatory, the only instrument with the spatial resolution to distinguish the magnetar from the supermassive black hole (2.4" angular distance). It was monitored from its outburst onset in April 2013 until August 2019, collecting more than fifty Chandra observations for a total of more than 2.3 Ms of data. Soon after the outburst onset, the magnetar emission settled onto a purely thermal emission state that cooled from a temperature of about 0.9 to 0.6 keV over 6 years. The pulsar timing properties showed at least two changes in the period derivative, increasing by a factor of about 4 during the outburst decay. We find that the long-term properties of this outburst challenge current models for the magnetar outbursts. △ Less

Submitted 24 March, 2020; v1 submitted 16 March, 2020; originally announced March 2020.

Comments: 11 pages, 6 figures. Accepted by ApJ

Search for the optical counterpart of the GW170814 gravitationalwave event with the VLT Survey Telescope

Authors: A. Grado, E. Cappellaro, S. Covino, F. Getman, G. Greco, L. Limatola, S. Yang, L. Amati, S. Benetti, M. Branchesi, E. Brocato, M. Botticella, S. Campana, M. Cantiello, M. Dadina, F. D'Ammando, G. De Cesare, V. D'Elia, M. Della Valle, E. Iodice, G. Longo, M. Mapelli, N. Masetti, L. Nicastro, E. Palazzi , et al. (8 additional authors not shown)

Abstract: We report on the search for the optical counterpart of the gravitational event GW170814, which was carried out with the VLT Survey Telescope (VST) by the GRAvitational Wave Inaf TeAm (GRAWITA). Observations started 17.5 hours after the LIGO/Virgo alert and we covered an area of 99 deg$^2$ that encloses $\sim 77\%$ and $\sim 59\%$ of the initial and refined localization probability regions, respect… ▽ More We report on the search for the optical counterpart of the gravitational event GW170814, which was carried out with the VLT Survey Telescope (VST) by the GRAvitational Wave Inaf TeAm (GRAWITA). Observations started 17.5 hours after the LIGO/Virgo alert and we covered an area of 99 deg$^2$ that encloses $\sim 77\%$ and $\sim 59\%$ of the initial and refined localization probability regions, respectively. A total of six epochs were secured over nearly two months. The survey reached an average limiting magnitude of 22 AB mag in the $r-$band. After assuming the model described in Perna et al. 2019 that derives as possible optical counterpart of a BBH event a transient source declining in about one day, we have computed a survey efficiency of about $5\%$. This paper describes the VST observational strategy and the results obtained by our analysis pipelines developed to search for optical transients in multi-epoch images. We report the catalogue of the candidates with possible identifications based on light-curve fitting. We have identified two dozens of SNe, nine AGNs, one QSO. Nineteen transients characterized by a single detection were not classified. We have restricted our analysis only to the candidates that fall into the refined localization map. None out of 39 left candidates could be positively associated with GW170814. This result implies that the possible emission of optical radiation from a BBH merger had to be fainter than r $\sim$ 22 ($L_{optical}$ $\sim$ $1.4 \times 10^{42}$ erg/s) on a time interval ranging from a few hours up to two months after the GW event. △ Less

Submitted 27 January, 2020; originally announced January 2020.

Comments: 24 pages, 7 figures, MNRAS, 492, 1731 (2020)

The Ultraluminous X-ray sources population of the galaxy NGC 7456

Authors: F. Pintore, M. Marelli, R. Salvaterra, G. L. Israel, G. A. Rodríguez Castillo, P. Esposito, A. Belfiore, A. De Luca, A. Wolter, S. Mereghetti, L. Stella, M. Rigoselli, H. P. Earnshaw, C. Pinto, T. P. Roberts, D. J. Walton, F. Bernardini, F. Haberl, C. Salvaggio, A. Tiengo, L. Zampieri, M. Bachetti, M. Brightman, P. Casella, D. D'Agostino , et al. (6 additional authors not shown)

Abstract: Ultraluminous X-ray sources (ULXs) are a class of accreting compact objects with X-ray luminosities above 1e39 erg/s. The ULX population counts several hundreds objects but only a minor fraction is well studied. Here we present a detailed analysis of all ULXs hosted in the galaxy NGC 7456. It was observed in X-rays only once in the past (in 2005) by XMM-Newton, but the observation was short and st… ▽ More Ultraluminous X-ray sources (ULXs) are a class of accreting compact objects with X-ray luminosities above 1e39 erg/s. The ULX population counts several hundreds objects but only a minor fraction is well studied. Here we present a detailed analysis of all ULXs hosted in the galaxy NGC 7456. It was observed in X-rays only once in the past (in 2005) by XMM-Newton, but the observation was short and strongly affected by high background. In 2018, we obtained a new, deeper (~90 ks) XMM-Newton observation that allowed us to perform a detailed characterization of the ULXs hosted in the galaxy. ULX-1 and ULX-2, the two brightest objects (Lx~(6-10)e39 erg/s), have spectra that can be described by a two-thermal component model as often found in ULXs. ULX-1 shows also one order of magnitude in flux variability on short-term timescales (hundreds to thousand ks). The other sources (ULX-3 and ULX-4) show flux changes of at least an order of magnitude, and these objects may be candidate transient ULXs although longer X-ray monitoring or further studies are required to ascribe them to the ULX population. In addition, we found a previously undetected source that might be a new candidate ULX (labelled as ULX-5) with a luminosity of ~1e39 erg/s and hard power-law spectral shape, whose nature is still unclear and for which a background Active Galactic Nucleus cannot be excluded. We discuss the properties of all the ULXs in NGC 7456 within the framework of super-Eddington accretion onto stellar mass compact objects. Although no pulsations were detected, we cannot exclude that the sources host neutron stars. △ Less

Submitted 27 January, 2020; v1 submitted 23 January, 2020; originally announced January 2020.

Comments: Accepted on ApJ; 10 pages, 5 figures, 3 tables

Three-dimensional general relativistic Poynting-Robertson effect II: Radiation field from a rigidly rotating spherical source

Authors: Pavel Bakala, Vittorio De Falco, Emmanuele Battista, Kateřina Goluchová, Debora Lančová, Maurizio Falanga, Luigi Stella

Abstract: We investigate the three-dimensional, general relativistic Poynting-Robertson effect in the case of rigidly rotating spherical source which emits radiation radially in the local comoving frame. Such radiation field is meant to approximate the field produced by the surface of a rotating neutron star, or by the central radiating hot corona of accreting black holes; it extends the purely radial radia… ▽ More We investigate the three-dimensional, general relativistic Poynting-Robertson effect in the case of rigidly rotating spherical source which emits radiation radially in the local comoving frame. Such radiation field is meant to approximate the field produced by the surface of a rotating neutron star, or by the central radiating hot corona of accreting black holes; it extends the purely radial radiation field that we considered in a previous study. Its angular momentum is expressed in terms of the rotation frequency and radius of the emitting source. For the background we adopt a Kerr spacetime geometry. We derive the equations of motion for test particles influenced by such radiation field, recovering the classical and weak-field approximation for slow rotation. We concentrate on solutions consisting of particles orbiting along circular orbits off and parallel to the equatorial plane, which are stabilized by the balance between gravitational attraction, radiation force and PR drag. Such solutions are found to lie on a critical hypersurface, whose shape may morph from prolate to oblate depending on the Kerr spin parameter and the luminosity, rotation and radius of the radiating sphere. For selected parameter ranges, the critical hypersurface intersects the radiating sphere giving rise to a bulging equatorial region or, alternatively, two lobes above the poles. We calculate the trajectories of test particles in the close vicinity of the critical hypersurface for a selected set of initial parameters and analyze the spatial and angular velocity of test particles captured on the critical hypersurface. △ Less

Submitted 1 November, 2019; originally announced November 2019.

Comments: 20 pages, 13 Figures

Diffuse X-ray emission around an ultraluminous X-ray pulsar

Authors: Andrea Belfiore, Paolo Esposito, Fabio Pintore, Giovanni Novara, Ruben Salvaterra, Andrea De Luca, Andrea Tiengo, Patrizia Caraveo, Felix Fuerst, Gian Luca Israel, Danilo Magistrali, Martino Marelli, Sandro Mereghetti, Alessandro Papitto, Guillermo Rodriguez Castillo, Chiara Salvaggio, Luigi Stella, Dominic Walton, Anna Wolter, Luca Zampieri

Abstract: Ultraluminous X-ray sources (ULXs) are extragalactic X-ray emitters located off-center of their host galaxy and with a luminosity in excess of a few ${10^{39}\text{ erg s}^{-1}}$, if emitted isotropically. The discovery of periodic modulation revealed that in some ULXs the accreting compact object is a neutron star, indicating luminosities substantially above their Eddington limit. The most extrem… ▽ More Ultraluminous X-ray sources (ULXs) are extragalactic X-ray emitters located off-center of their host galaxy and with a luminosity in excess of a few ${10^{39}\text{ erg s}^{-1}}$, if emitted isotropically. The discovery of periodic modulation revealed that in some ULXs the accreting compact object is a neutron star, indicating luminosities substantially above their Eddington limit. The most extreme object in this respect is ${NGC 5907~ULX-1}$ (ULX1), with a peak luminosity that is 500 times its Eddington limit. During a Chandra observation to probe a low state of ULX1, we detected diffuse X-ray emission at the position of ULX1. Its diameter is $2.7 \pm 1.0$ arcsec and contains 25 photons, none below 0.8 keV. We interpret this extended structure as an expanding nebula powered by the wind of ULX1. Its diameter of about ${200\text{ pc}}$, characteristic energy of ${\sim 1.9\text{ keV}}$, and luminosity of ${\sim 2\times10^{38}\text{ erg s}^{-1}}$ imply a mechanical power of ${1.3\times10^{41}\text{ erg s}^{-1}}$ and an age ${\sim 7 \times 10^{4}\text{ yr}}$. This interpretation suggests that a genuinely super-Eddington regime can be sustained for time scales much longer than the spin-up time of the neutron star powering the system. As the mechanical power from a single ULX nebula can rival the injection rate of cosmic rays of an entire galaxy, ULX nebulae could be important cosmic ray accelerators. △ Less

Submitted 25 October, 2019; originally announced October 2019.

Comments: 28 pages, 6 figures, 1 table - accepted for publication

A new approximation of photon geodesics in Schwarzschild spacetime

Authors: Riccardo La Placa, Pavel Bakala, Luigi Stella, Maurizio Falanga

Abstract: In this research note we introduce a new approximation of photon geodesics in Schwarzschild spacetime which is especially useful to describe highly bent trajectories, for which the angle between the initial emission position and the line of sight to the observer approaches $π$: this corresponds to the points behind the central mass of the Schwarzschild metric with respect to the observer. The appr… ▽ More In this research note we introduce a new approximation of photon geodesics in Schwarzschild spacetime which is especially useful to describe highly bent trajectories, for which the angle between the initial emission position and the line of sight to the observer approaches $π$: this corresponds to the points behind the central mass of the Schwarzschild metric with respect to the observer. The approximation maintains very good accuracy overall, with deviations from the exact numerical results below $1\%$ up to the innermost stable circular orbit (ISCO) located at $6~GM/c^2$. △ Less

Submitted 25 July, 2019; originally announced July 2019.

Comments: 3 pages, 1 figure, 3 subfigures

Journal ref: Research Notes of the AAS, Volume 3 (2019), 99

GluonTS: Probabilistic Time Series Models in Python

Authors: Alexander Alexandrov, Konstantinos Benidis, Michael Bohlke-Schneider, Valentin Flunkert, Jan Gasthaus, Tim Januschowski, Danielle C. Maddix, Syama Rangapuram, David Salinas, Jasper Schulz, Lorenzo Stella, Ali Caner Türkmen, Yuyang Wang

Abstract: We introduce Gluon Time Series (GluonTS, available at https://gluon-ts.mxnet.io), a library for deep-learning-based time series modeling. GluonTS simplifies the development of and experimentation with time series models for common tasks such as forecasting or anomaly detection. It provides all necessary components and tools that scientists need for quickly building new models, for efficiently runn… ▽ More We introduce Gluon Time Series (GluonTS, available at https://gluon-ts.mxnet.io), a library for deep-learning-based time series modeling. GluonTS simplifies the development of and experimentation with time series models for common tasks such as forecasting or anomaly detection. It provides all necessary components and tools that scientists need for quickly building new models, for efficiently running and analyzing experiments and for evaluating model accuracy. △ Less

Submitted 14 June, 2019; v1 submitted 12 June, 2019; originally announced June 2019.

Comments: ICML Time Series Workshop 2019

Discovery of a 2.8 s pulsar in a 2 d orbit High-Mass X-ray Binary powering the Ultraluminous X-ray source ULX-7 in M51

Authors: G. A. Rodríguez Castillo, G. L. Israel, A. Belfiore, F. Bernardini, P. Esposito, F. Pintore, A. De Luca, A. Papitto, L. Stella, A. Tiengo, L. Zampieri, M. Bachetti, M. Brightman, P. Casella, D. D'Agostino, S. Dall'Osso, H. P. Earnshaw, F. Fürst, F. Haberl, F. A. Harrison, M. Mapelli, M. Marelli, M. Middleton, C. Pinto, T. P. Roberts , et al. (4 additional authors not shown)

Abstract: We discovered 2.8 s pulsations in the X-ray emission of the ultraluminous X-ray source (ULX) M51 ULX-7 within the UNSEeN project, which was designed to hunt for new pulsating ULXs (PULXs) with XMM-Newton. The pulse shape is sinusoidal and large variations of its amplitude were observed even within single exposures (pulsed fraction from less than 5% to 20%). M51 ULX-7 is a variable source, generall… ▽ More We discovered 2.8 s pulsations in the X-ray emission of the ultraluminous X-ray source (ULX) M51 ULX-7 within the UNSEeN project, which was designed to hunt for new pulsating ULXs (PULXs) with XMM-Newton. The pulse shape is sinusoidal and large variations of its amplitude were observed even within single exposures (pulsed fraction from less than 5% to 20%). M51 ULX-7 is a variable source, generally observed at an X-ray luminosity between $10^{39}$ and $10^{40}$ erg s$^{-1}$, located in the outskirts of the spiral galaxy M51a at a distance of 8.6 Mpc. According to our analysis, the X-ray pulsar orbits in a 2-d binary with a projected semi-major axis $a_\mathrm{X} \sin i \simeq$ 28 lt-s. For a neutron star (NS) of 1.4 $M_{\odot}$, this implies a lower limit on the companion mass of 8 $M_{\odot}$, placing the system hosting M51 ULX-7 in the high-mass X-ray binary class. The barycentric pulse period decreased by $\simeq$0.4 ms in the 31 d spanned by our May -- June 2018 observations, corresponding to a spin-up rate $\dot{P} \simeq -1.5\times10^{-10}\text{s s}^{-1}$. In an archival 2005 XMM-Newton exposure, we measured a spin period of $\sim$3.3 s, indicating a secular spin-up of $\dot{P}_{\mathrm{sec}}\simeq -10^{-9}\text{ s s}^{-1}$, a value in the range of other known PULXs. Our findings suggest that the system consists of an OB giant and a moderately magnetic (dipole field component in the range $10^{12}$ G $\lesssim B_{\mathrm{dip}}\lesssim 10^{13}$G) accreting NS with weakly beamed emission ($1/12\lesssim b\lesssim1/4$). △ Less

Submitted 11 May, 2020; v1 submitted 11 June, 2019; originally announced June 2019.

Comments: Accepted for publication in The Astrophysical Journal

A new method to constrain neutron star structure from quasi-periodic oscillations

Authors: Andrea Maselli, George Pappas, Paolo Pani, Leonardo Gualtieri, Sara Motta, Valeria Ferrari, Luigi Stella

Abstract: We develop a new method to measure neutron star parameters and derive constraints on the equation of state of dense matter by fitting the frequencies of simultaneous Quasi Periodic Oscillation modes observed in the X-ray flux of accreting neutron stars in low mass X-ray binaries. To this aim we calculate the fundamental frequencies of geodesic motion around rotating neutron stars based on an accur… ▽ More We develop a new method to measure neutron star parameters and derive constraints on the equation of state of dense matter by fitting the frequencies of simultaneous Quasi Periodic Oscillation modes observed in the X-ray flux of accreting neutron stars in low mass X-ray binaries. To this aim we calculate the fundamental frequencies of geodesic motion around rotating neutron stars based on an accurate general-relativistic approximation for their external spacetime. Once the fundamental frequencies are related to the observed frequencies through a QPO model, they can be fit to the data to obtain estimates of the three parameters describing the spacetime, namely the neutron star mass, angular momentum and quadrupole moment. From these parameters we derive information on the neutron star structure and equation of state. We present a proof of principle of our method applied to pairs of kHz QPO frequencies observed from three systems (4U1608-52, 4U0614+09 and 4U1728-34). We identify the kHz QPOs with the azimuthal and the periastron precession frequencies of matter orbiting the neutron star, and via our Bayesian inference technique we derive constraints on the neutrons stars' masses and radii. This method is applicable to other geodesic-frequency-based QPO models. △ Less

Submitted 30 June, 2020; v1 submitted 12 May, 2019; originally announced May 2019.

Comments: Changes to match the version in press on ApJ

Pulsating in unison at optical and X-ray energies: simultaneous high-time resolution observations of the transitional millisecond pulsar PSR J1023+0038

Authors: A. Papitto, F. Ambrosino, L. Stella, D. F. Torres, F. Coti Zelati, A. Ghedina, F. Meddi, A. Sanna, P. Casella, Y. Dallilar, S. Eikenberry, G. L. Israel, F. Onori, S. Piranomonte, E. Bozzo, L. Burderi, S. Campana, D. de Martino, T. Di Salvo, C. Ferrigno, N. Rea, A. Riggio, S. Serrano, A. Veledina, L. Zampieri

Abstract: PSR J1023+0038 is the first millisecond pulsar discovered to pulsate in the visible band; such a detection took place when the pulsar was surrounded by an accretion disk and also showed X-ray pulsations. We report on the first high time resolution observational campaign of this transitional pulsar in the disk state, using simultaneous observations in the optical (TNG, NOT, TJO), X-ray (XMM-Newton,… ▽ More PSR J1023+0038 is the first millisecond pulsar discovered to pulsate in the visible band; such a detection took place when the pulsar was surrounded by an accretion disk and also showed X-ray pulsations. We report on the first high time resolution observational campaign of this transitional pulsar in the disk state, using simultaneous observations in the optical (TNG, NOT, TJO), X-ray (XMM-Newton, NuSTAR, NICER), infrared (GTC) and UV (Swift) bands. Optical and X-ray pulsations were detected simultaneously in the X-ray high intensity mode in which the source spends $\sim$ 70% of the time, and both disappeared in the low mode, indicating a common underlying physical mechanism. In addition, optical and X-ray pulses were emitted within a few km, had similar pulse shape and distribution of the pulsed flux density compatible with a power-law relation $F_ν \propto ν^{-0.7}$ connecting the optical and the 0.3-45 keV X-ray band. Optical pulses were detected also during flares with a pulsed flux reduced by one third with respect to the high mode; the lack of a simultaneous detection of X-ray pulses is compatible with the lower photon statistics. We show that magnetically channeled accretion of plasma onto the surface of the neutron star cannot account for the optical pulsed luminosity ($\sim 10^{31}$ erg/s). On the other hand, magnetospheric rotation-powered pulsar emission would require an extremely efficient conversion of spin-down power into pulsed optical and X-ray emission. We then propose that optical and X-ray pulses are instead produced by synchrotron emission from the intrabinary shock that forms where a striped pulsar wind meets the accretion disk, within a few light cylinder radii away, $\sim$ 100 km, from the pulsar. △ Less

Submitted 18 June, 2019; v1 submitted 23 April, 2019; originally announced April 2019.

Comments: 26 pages, 14 figures, first submitted to ApJ on 2019, January 18

Detailed X-ray spectroscopy of the magnetar 1E 2259+586

Authors: D. Pizzocaro, A. Tiengo, S. Mereghetti, R. Turolla, P. Esposito, L. Stella, S. Zane, N. Rea, F. Coti Zelati, G. Israel

Abstract: Magnetic field geometry is expected to play a fundamental role in magnetar activity. The discovery of a phase-variable absorption feature in the X-ray spectrum of SGR 0418+5729, interpreted as cyclotron resonant scattering, suggests the presence of very strong non-dipolar components in the magnetic fields of magnetars. We performed a deep XMM-Newton observation of pulsar 1E 2259+586, to search for… ▽ More Magnetic field geometry is expected to play a fundamental role in magnetar activity. The discovery of a phase-variable absorption feature in the X-ray spectrum of SGR 0418+5729, interpreted as cyclotron resonant scattering, suggests the presence of very strong non-dipolar components in the magnetic fields of magnetars. We performed a deep XMM-Newton observation of pulsar 1E 2259+586, to search for spectral features due to intense local magnetic fields. In the phase-averaged X-ray spectrum, we found evidence for a broad absorption feature at very low energy (0.7 keV). If the feature is intrinsic to the source, it might be due to resonant scattering/absorption by protons close to star surface. The line energy implies a magnetic field of ~ 10^14 G, roughly similar to the spin-down measure, ~ 6x10^13 G. Examination of the X-ray phase-energy diagram shows evidence for a further absorption feature, the energy of which strongly depends on the rotational phase (E >~ 1 keV ). Unlike similar features detected in other magnetar sources, notably SGR 0418+5729, it is too shallow and limited to a small phase interval to be modeled with a narrow phase-variable cyclotron absorption line. A detailed phase-resolved spectral analysis reveals significant phase-dependent variability in the continuum, especially above 2 keV. We conclude that all the variability with phase in 1E 2259+586 can be attributed to changes in the continuum properties which appear consistent with the predictions of the Resonant Compton Scattering model. △ Less

Submitted 30 May, 2019; v1 submitted 16 April, 2019; originally announced April 2019.

Journal ref: A&A 626, A39 (2019)