Showing 1–50 of 254 results for author: Lopez, M

Testing particle acceleration in blazar jets with continuous high-cadence optical polarization observations

Authors: Ioannis Liodakis, Sebastian Kiehlmann, Alan P. Marscher, Haocheng Zhang, Dmitry Blinov, Svetlana G. Jorstad, Iván Agudo, Erika Benítez, Andrei Berdyugin, Giacomo Bonnoli, Carolina Casadio, Chien-Ting Chen, Wen-Ping Chen, Steven R. Ehlert, Juan Escudero, Tatiana S. Grishina, David Hiriart, Angela Hsu, Ryo Imazawa, Helen E. Jermak, Jincen Jose, Philip Kaaret, Evgenia N. Kopatskaya, Bhavana Lalchand, Elena G. Larionova , et al. (22 additional authors not shown)

Abstract: Variability can be the pathway to understanding the physical processes in astrophysical jets, however, the high-cadence observations required to test particle acceleration models are still missing. Here we report on the first attempt to produce continuous, >24 hour polarization light curves of blazars using telescopes distributed across the globe and the rotation of the Earth to avoid the rising S… ▽ More Variability can be the pathway to understanding the physical processes in astrophysical jets, however, the high-cadence observations required to test particle acceleration models are still missing. Here we report on the first attempt to produce continuous, >24 hour polarization light curves of blazars using telescopes distributed across the globe and the rotation of the Earth to avoid the rising Sun. Our campaign involved 16 telescopes in Asia, Europe, and North America. We observed BL Lacertae and CGRaBS J0211+1051 for a combined 685 telescope hours. We find large variations in the polarization degree and angle for both sources in sub-hour timescales as well as a ~180 degree rotation of the polarization angle in CGRaBS J0211+1051 in less than two days. We compared our high-cadence observations to Particle-In-Cell magnetic reconnection and turbulent plasma simulations. We find that although the state of the art simulation frameworks can produce a large fraction of the polarization properties, they do not account for the entirety of the observed polarization behavior in blazar jets. △ Less

Submitted 21 June, 2024; originally announced June 2024.

Comments: 20 pages, 15 figures, 2 tables, accepted for publication in A&A. The data used in the paper are available here: https://doi.org/10.7910/DVN/IETSXS

Confirming the Explosive Dispersal Outflow in DR21 with ALMA

Authors: E. Guzmán Ccolque, M. Fernández López, L. A. Zapata, J. Bally, P. R. Rivera-Ortiz

Abstract: We present Atacama Large Millimeter/submillimeter Array (ALMA) 1.3 mm continuum and CO(2-1) line emission observations toward the high-mass star formation region DR21. Five new continuum sources are found. We identify eighteen outflow streamers detected in CO emission radially arising from a common origin. The velocity spread of the outflow streamers range between $-$100 to $+$70 km s$^{-1}$. The… ▽ More We present Atacama Large Millimeter/submillimeter Array (ALMA) 1.3 mm continuum and CO(2-1) line emission observations toward the high-mass star formation region DR21. Five new continuum sources are found. We identify eighteen outflow streamers detected in CO emission radially arising from a common origin. The velocity spread of the outflow streamers range between $-$100 to $+$70 km s$^{-1}$. The radial velocities of each outflow roughly follow linear gradients (Hubble-Lemaitre-like expansion motions). Using the CO emission of the whole ensemble of streamers we estimate a total outflow mass of 120-210M$_{\odot}$. Additionally, we derived the dynamical age (8600 yr), momentum ($\sim10^{3}$ M$_{\odot}$ km s$^{-1}$), and kinetic energy ($\sim10^{48}$ erg) of the outflow. The morphology and kinematics presented by the CO outflow streamers confirm the presence of an explosive dispersal outflow at the heart of DR21. Five dispersal explosive outflows associated with massive star-forming regions have been confirmed in our Galaxy (Orion BN/KL, G5.89-0.39, S106-IR, IRAS16076-5134 and IRAS 12326-6245). However, their frequency of occurrence in the Galaxy and the originating nature are still uncertain. △ Less

Submitted 5 June, 2024; originally announced June 2024.

Comments: 19 pages, 15 figures including 4 in the appendix, accepted for publication to A&A

Stellar angular momentum of intermediate redshift galaxies in MUSE surveys

Authors: Constanza Muñoz López, Davor Krajnović, B. Epinat, Y. Herrero-Alonso, T. Urrutia, W. Mercier, N. F. Bouché, L. A. Boogaard, T. Contini, L. Michel-Dansac, I. Pessa

Abstract: We quantify the stellar rotation of galaxies by computing the $λ_{R}$ parameter, a proxy for the stellar angular momentum in a sample of 106 galaxies with redshift 0.1 $<$ z $<$ 0.8 and stellar masses from $\sim$10$^{7.5}$ to 10$^{11.8}$ M$_{\odot}$. The sample is located in the CANDELS/GOODS-S and COSMOS fields, and it was observed by various MUSE surveys. We create stellar velocity and velocity… ▽ More We quantify the stellar rotation of galaxies by computing the $λ_{R}$ parameter, a proxy for the stellar angular momentum in a sample of 106 galaxies with redshift 0.1 $<$ z $<$ 0.8 and stellar masses from $\sim$10$^{7.5}$ to 10$^{11.8}$ M$_{\odot}$. The sample is located in the CANDELS/GOODS-S and COSMOS fields, and it was observed by various MUSE surveys. We create stellar velocity and velocity dispersion maps using a full-spectrum fitting technique, covering spatially $\sim$2$R_{e}$ for the galaxies. We study the impact of the atmospheric seeing on the spin parameter and apply corrections when pertinent. Through the analysis of the $λ_{R}-ε$ diagram, we notice that the fraction of round and massive galaxies increases with redshift. We lack galaxies with $λ_{R}$ < 0.1 in the sample and we find only one potential, but uncertain, low-mass slow rotator at z $\sim0.3$. Moreover, we do not see an evident evolution or trend in the stellar angular momentum with redshift. We characterize the sample environment using two indicators: a local estimator based on the Voronoi tesselation method, and a global estimator derived by the use of the Friends-of-Friends algorithm. We find no correlation between the environment and $λ_{R}$ given that we are not probing dense regions or massive galaxy structures. We also analyze the kinematic maps of the sample finding that about 40$\%$ of galaxies are consistent with being regular rotators, having rotating stellar discs with flat velocity dispersion maps, while $\sim20\%$ of galaxies have complex velocity maps and can be identified as non-regular rotators in spite of their $λ_{R}$ values. For the remaining galaxies the classification is uncertain. As we lack galaxies with $λ_{R}$< 0.1, we are not able to identify when galaxies become slow rotators within the surveyed environments, area and redshift range. △ Less

Submitted 16 April, 2024; originally announced April 2024.

Comments: 26 pages, 18 figures, full abstract on the paper

Episodicity in accretion-ejection processes associated with IRAS 15398-3359

Authors: E. Guzmán Ccolque, M. Fernández López, M. M. Vazzano, I. de Gregorio, A. Plunkett, A. Santamaría-Miranda

Abstract: The protostar IRAS 15398-3359 is associated with a bipolar molecular outflow ejected in an nearly northeast-southwest (NE-SW) direction which has been extensively studied. It has been suggested previous episodic accretion events by this source. Furthermore, the analysis of the morphology and kinematics of the molecular outflow revealed the presence of four $^{12}$CO(2-1) bipolar elliptical shock-l… ▽ More The protostar IRAS 15398-3359 is associated with a bipolar molecular outflow ejected in an nearly northeast-southwest (NE-SW) direction which has been extensively studied. It has been suggested previous episodic accretion events by this source. Furthermore, the analysis of the morphology and kinematics of the molecular outflow revealed the presence of four $^{12}$CO(2-1) bipolar elliptical shock-like structures identified in both lobes. These structures seem to trace different ejections inclined $\sim$10°on the plane of the sky from each other. This led to the hypothesis that the outflow axis likely precesses and launches material episodically. We analyze ALMA archive observations in Band 6, revealing the presence of low-velocity ($<3.5$km s$^{-1}$) emission from the line $^{12}$CO(2-1) to the south and north of the protostar. We study the morphology and kinematics of the gas, which seems to support the hypothesis of a precessing episodic outflow. The ALMA observations reveal a north-south (N-S) outflow most likely associated with the IRAS 15398-3359 protostellar system. This outflow could be older than the well-studied NE-SW outflow. The orientation of the N-S outflow is 50°- 60°on the plane of the sky away from that of the NE-SW outflow. We also analyze the Spectral Energy Distribution of a far away young star and preliminary discard it as the driver of the SE outflow remnants. The new observations support the hypothesis of strong episodic accretion-ejection events in IRAS 15398-3359, accompanied by dramatic changes in the orientation of its ejection axis, implying that all the outflows in the region may have been driven by the same protostar. △ Less

Submitted 19 March, 2024; originally announced March 2024.

Comments: 7 pages, 5 figures. Accepted by A&A

Journal ref: A&A 686, A143 (2024)

A Big Ring on the Sky

Authors: Alexia M. Lopez, Roger G. Clowes, Gerard M. Williger

Abstract: We present the discovery of `A Big Ring on the Sky' (BR), the second ultra-large large-scale structure (uLSS) found in MgII-absorber catalogues, following the previously reported Giant Arc (GA). In cosmological terms the BR is close to the GA - at the same redshift $z \sim 0.8$ and with a separation on the sky of only $\sim 12^\circ$. Two extraordinary uLSSs in such close configuration raises the… ▽ More We present the discovery of `A Big Ring on the Sky' (BR), the second ultra-large large-scale structure (uLSS) found in MgII-absorber catalogues, following the previously reported Giant Arc (GA). In cosmological terms the BR is close to the GA - at the same redshift $z \sim 0.8$ and with a separation on the sky of only $\sim 12^\circ$. Two extraordinary uLSSs in such close configuration raises the possibility that together they form an even more extraordinary cosmological system. The BR is a striking circular, annulus-like, structure of diameter $\sim 400$ Mpc (proper size, present epoch). The method of discovery is as described in the GA paper, but here using the new MgII-absorber catalogues restricted to DR16Q quasars. Using the Convex Hull of Member Spheres (CHMS) algorithm, we estimate that the annulus and inner absorbers of the BR have departures from random expectations, at the density of the control field, of up to $5.2σ$. We present the discovery of the BR, assess its significance using the CHMS, Minimal Spanning Tree (MST), FilFinder and Cuzick & Edwards (CE) methods, show it in the context of the GA+BR system, and suggest some implications for the origins of uLSS and for our understanding of cosmology. For example, it may be that unusual geometric patterns, such as these uLSSs, have an origin in cosmic strings. △ Less

Submitted 10 July, 2024; v1 submitted 12 February, 2024; originally announced February 2024.

Comments: Accepted for publication in JCAP

Kilometric type II radio emissions in Wind/WAVES TNR data and association with interplanetary structures near Earth

Authors: Franco Manini, Hebe Cremades, Fernando M. López, Teresa Nieves-Chinchilla

Abstract: Type II radio bursts arise as a consequence of shocks typically driven by coronal mass ejections (CMEs). When these shocks propagate outward from the Sun, their associated radio emissions drift down in frequency as excited particles emit at the local plasma frequency, creating the usual Type II patterns. In this work, we use dynamic spectra from the Wind/WAVES Thermal Noise Receiver (TNR) to ident… ▽ More Type II radio bursts arise as a consequence of shocks typically driven by coronal mass ejections (CMEs). When these shocks propagate outward from the Sun, their associated radio emissions drift down in frequency as excited particles emit at the local plasma frequency, creating the usual Type II patterns. In this work, we use dynamic spectra from the Wind/WAVES Thermal Noise Receiver (TNR) to identify Type II radio emissions in the kilometric wavelength range (kmTIIs, f < 300 kHz) between 1 January 2000 and 31 December 2012, i.e. over a solar cycle. We identified 134 kmTII events and compiled various characteristics for each of them. Of particular importance is the finding of 45 kmTII events not reported by the official Wind/WAVES catalog (based on RAD1 and RAD2 data). We search for associations with interplanetary structures and analyze their main characteristics, in order to reveal distinctive attributes that may correlate with the occurrence of kmTII emission. We find that the fraction of interplanetary coronal mass ejections (ICMEs) classified as magnetic clouds (MCs) that are associated with kmTIIs is roughly similar to that of MCs not associated with kmTIIs. Conversely, the fraction of ICMEs with bidirectional electrons is significantly larger for those ICMEs associated with kmTIIs (74% vs. 48%). Likewise, ICMEs associated with kmTIIs are on average 23% faster. The disturbance storm time (DsT) mean value is almost twice as large for kmTII-associated ICMEs, indicating that they tend to produce intense geomagnetic storms. In addition, the proportion of ICMEs producing moderate to intense geomagnetic storms is twice as large for the kmTII-associated ICMEs. After this investigation, TNR data prove to be valuable not only as complementary data for the analysis of kmTII events but also for forecasting the arrival of shocks at Earth. △ Less

Submitted 23 November, 2023; v1 submitted 14 November, 2023; originally announced November 2023.

Comments: Revised version accepted for publication in Solar Physics on 2023-Nov-23 (21 pages, 5 figures)

Detection of anomalies amongst LIGO's glitch populations with autoencoders

Authors: Paloma Laguarta, Robin van der Laag, Melissa Lopez, Tom Dooney, Andrew L. Miller, Stefano Schmidt, Marco Cavaglia, Sarah Caudill, Kurt Driessens, Jöel Karel, Roy Lenders, Chris Van Den Broeck

Abstract: Gravitational-wave (GW) interferometers are able to detect a change in distance of $\sim$ 1/10,000th the size of a proton. Such sensitivity leads to large appearance rates of non-Gaussian transient noise bursts in the main detector strain, also known as glitches. These glitches come in a wide range of frequency-amplitude-time morphologies and are caused by environmental or instrumental processes,… ▽ More Gravitational-wave (GW) interferometers are able to detect a change in distance of $\sim$ 1/10,000th the size of a proton. Such sensitivity leads to large appearance rates of non-Gaussian transient noise bursts in the main detector strain, also known as glitches. These glitches come in a wide range of frequency-amplitude-time morphologies and are caused by environmental or instrumental processes, hindering searches for all sources of gravitational waves. Current approaches for their identification use supervised models to learn their morphology in the main strain, but do not consider relevant information provided by auxiliary channels that monitor the state of the interferometers nor provide a flexible framework for novel glitch morphologies. In this work, we present an unsupervised algorithm to find anomalous glitches. We encode a subset of auxiliary channels from LIGO Livingston in the fractal dimension, a measure for the complexity of the data, and learn the underlying distribution of the data using an auto-encoder with periodic convolutions. In this way, we uncover unknown glitch morphologies, and overlaps in time between different glitches and misclassifications. This led to the discovery of anomalies in $6.6 \%$ of the input data. The results of this investigation stress the learnable structure of auxiliary channels encoded in fractal dimension and provide a flexible framework to improve the state-of-the-art of glitch identification algorithms. △ Less

Submitted 5 October, 2023; originally announced October 2023.

Comments: 14 pages, 13 figues, to be published in Classical and Quantum Gravity

First light of VLT/HiRISE: High-resolution spectroscopy of young giant exoplanets

Authors: A. Vigan, M. El Morsy, M. Lopez, G. P. P. L. Otten, J. Garcia, J. Costes, E. Muslimov, A. Viret, Y. Charles, G. Zins, G. Murray, A. Costille, J. Paufique, U. Seemann, M. Houllé, H. Anwand-Heerwart, M. Phillips, A. Abinanti, P. Balard, I. Baraffe, J. -A. Benedetti, P. Blanchard, L. Blanco, J. -L. Beuzit, E. Choquet , et al. (24 additional authors not shown)

Abstract: A major endeavor of this decade is the direct characterization of young giant exoplanets at high spectral resolution to determine the composition of their atmosphere and infer their formation processes and evolution. Such a goal represents a major challenge owing to their small angular separation and luminosity contrast with respect to their parent stars. Instead of designing and implementing comp… ▽ More A major endeavor of this decade is the direct characterization of young giant exoplanets at high spectral resolution to determine the composition of their atmosphere and infer their formation processes and evolution. Such a goal represents a major challenge owing to their small angular separation and luminosity contrast with respect to their parent stars. Instead of designing and implementing completely new facilities, it has been proposed to leverage the capabilities of existing instruments that offer either high contrast imaging or high dispersion spectroscopy, by coupling them using optical fibers. In this work we present the implementation and first on-sky results of the HiRISE instrument at the very large telescope (VLT), which combines the exoplanet imager SPHERE with the recently upgraded high resolution spectrograph CRIRES using single-mode fibers. The goal of HiRISE is to enable the characterization of known companions in the $H$ band, at a spectral resolution of the order of $R = λ/Δλ= 100\,000$, in a few hours of observing time. We present the main design choices and the technical implementation of the system, which is constituted of three major parts: the fiber injection module inside of SPHERE, the fiber bundle around the telescope, and the fiber extraction module at the entrance of CRIRES. We also detail the specific calibrations required for HiRISE and the operations of the instrument for science observations. Finally, we detail the performance of the system in terms of astrometry, temporal stability, optical aberrations, and transmission, for which we report a peak value of $\sim$3.9% based on sky measurements in median observing conditions. Finally, we report on the first astrophysical detection of HiRISE to illustrate its potential. △ Less

Submitted 22 November, 2023; v1 submitted 21 September, 2023; originally announced September 2023.

Comments: 17 pages, 15 figures, 3 tables. Accepted in A&A on 24 October 2023

Gravitational-Wave Searches for Cosmic String Cusps in Einstein Telescope Data using Deep Learning

Authors: Quirijn Meijer, Melissa Lopez, Daichi Tsuna, Sarah Caudill

Abstract: Gravitational-wave searches for cosmic strings are currently hindered by the presence of detector glitches, some classes of which strongly resemble cosmic string signals. This confusion greatly reduces the efficiency of searches. A deep-learning model is proposed for the task of distinguishing between gravitational wave signals from cosmic string cusps and simulated blip glitches in design sensiti… ▽ More Gravitational-wave searches for cosmic strings are currently hindered by the presence of detector glitches, some classes of which strongly resemble cosmic string signals. This confusion greatly reduces the efficiency of searches. A deep-learning model is proposed for the task of distinguishing between gravitational wave signals from cosmic string cusps and simulated blip glitches in design sensitivity data from the future Einstein Telescope. The model is an ensemble consisting of three convolutional neural networks, achieving an accuracy of 79%, a true positive rate of 76%, and a false positive rate of 18%. This marks the first time convolutional neural networks have been trained on a realistic population of Einstein Telescope glitches. On a dataset consisting of signals and glitches, the model is shown to outperform matched filtering, specifically being better at rejecting glitches. The behaviour of the model is interpreted through the application of several methods, including a novel technique called waveform surgery, used to quantify the importance of waveform sections to a classification model. In addition, a method to visualise convolutional neural network activations for one-dimensional time series is proposed and used. These analyses help further the understanding of the morphological differences between cosmic string cusp signals and blip glitches. Because of its classification speed in the order of magnitude of milliseconds, the deep-learning model is suitable for future use as part of a real-time detection pipeline. The deep-learning model is transverse and can therefore potentially be applied to other transient searches. △ Less

Submitted 22 January, 2024; v1 submitted 23 August, 2023; originally announced August 2023.

Comments: 15 pages, 20 figures

Multi-messenger characterization of Mrk 501 during historically low X-ray and $γ$-ray activity

Authors: MAGIC collaboration, H. Abe, S. Abe, V. A. Acciari, I. Agudo, T. Aniello, S. Ansoldi, L. A. Antonelli, A. Arbet Engels, C. Arcaro, M. Artero, K. Asano, D. Baack, A. Babić, A. Baquero, U. Barres de Almeida, J. A. Barrio, I. Batković, J. Baxter, J. Becerra González, W. Bednarek, E. Bernardini, M. Bernardos, A. Berti, J. Besenrieder , et al. (300 additional authors not shown)

Abstract: We study the broadband emission of Mrk 501 using multi-wavelength observations from 2017 to 2020 performed with a multitude of instruments, involving, among others, MAGIC, Fermi-LAT, NuSTAR, Swift, GASP-WEBT, and OVRO. Mrk 501 showed an extremely low broadband activity, which may help to unravel its baseline emission. Nonetheless, significant flux variations are detected at all wavebands, with the… ▽ More We study the broadband emission of Mrk 501 using multi-wavelength observations from 2017 to 2020 performed with a multitude of instruments, involving, among others, MAGIC, Fermi-LAT, NuSTAR, Swift, GASP-WEBT, and OVRO. Mrk 501 showed an extremely low broadband activity, which may help to unravel its baseline emission. Nonetheless, significant flux variations are detected at all wavebands, with the highest occurring at X-rays and very-high-energy (VHE) $γ$-rays. A significant correlation ($>$3$σ$) between X-rays and VHE $γ$-rays is measured, supporting leptonic scenarios to explain the variable parts of the emission, also during low activity. This is further supported when we extend our data from 2008 to 2020, and identify, for the first time, significant correlations between Swift-XRT and Fermi-LAT. We additionally find correlations between high-energy $γ$-rays and radio, with the radio lagging by more than 100 days, placing the $γ$-ray emission zone upstream of the radio-bright regions in the jet. Furthermore, Mrk 501 showed a historically low activity in X-rays and VHE $γ$-rays from mid-2017 to mid-2019 with a stable VHE flux ($>$0.2 TeV) of 5% the emission of the Crab Nebula. The broadband spectral energy distribution (SED) of this 2-year-long low-state, the potential baseline emission of Mrk 501, can be characterized with one-zone leptonic models, and with (lepto)-hadronic models fulfilling neutrino flux constraints from IceCube. We explore the time evolution of the SED towards the low-state, revealing that the stable baseline emission may be ascribed to a standing shock, and the variable emission to an additional expanding or traveling shock. △ Less

Submitted 5 March, 2023; v1 submitted 5 October, 2022; originally announced October 2022.

Comments: 55 pages, 30 figures, 14 tables, accepted by APJS. Corresponding authors are L. Heckmann, D. Paneque, S. Gasparyan, M. Cerruti, and N. Sahakyan

Journal ref: ApJS 266 37 (2023)

Connecting SPHERE and CRIRES+ for the characterisation of young exoplanets at high spectral resolution: status update of VLT/HiRISE

Authors: A. Vigan, M. Lopez, M. El Morsy, E. Muslimov, A. Viret, G. Zins, G. Murray, A. Costille, G. P. P. L. Otten, U. Seemann, H. Anwand-Heerwart, K. Dohlen, P. Blanchard, J. Garcia, Y. Charles, N. Tchoubaklian, T. Ely, M. Phillips, J. Paufique, J. -L. Beuzit, M. Houllé, J. Costes, R. Pourcelot, I. Baraffe, R. Dorn , et al. (10 additional authors not shown)

Abstract: New generation exoplanet imagers on large ground-based telescopes are highly optimised for the detection of young giant exoplanets in the near-infrared, but they are intrinsically limited for their characterisation by the low spectral resolution of their integral field spectrographs ($R<100$). High-dispersion spectroscopy at $R \gg 10^4$ would be a powerful tool for the characterisation of these p… ▽ More New generation exoplanet imagers on large ground-based telescopes are highly optimised for the detection of young giant exoplanets in the near-infrared, but they are intrinsically limited for their characterisation by the low spectral resolution of their integral field spectrographs ($R<100$). High-dispersion spectroscopy at $R \gg 10^4$ would be a powerful tool for the characterisation of these planets, but there is currently no high-resolution spectrograph with extreme adaptive optics and coronagraphy that would enable such characterisation. With project HiRISE we propose to use fiber coupling to combine the capabilities of two flagship instruments at the Very Large Telescope in Chile: the exoplanet imager SPHERE and the high-resolution spectrograph CRIRES+. The coupling will be implemented at the telescope in early 2023. We provide a general overview of the implementation of HiRISE, of its assembly, integration and testing (AIT) phase in Europe, and a brief assessment of its expected performance based on the final hardware. △ Less

Submitted 13 July, 2022; originally announced July 2022.

Comments: Submitted to SPIE Astronomical Telescopes + Instrumentation 2022, Adaptive Optics Systems VIII, Paper 12185-27

Is the Observable Universe Consistent with the Cosmological Principle?

Authors: Pavan Kumar Aluri, Paolo Cea, Pravabati Chingangbam, Ming-Chung Chu, Roger G. Clowes, Damien Hutsemékers, Joby P. Kochappan, Alexia M. Lopez, Lang Liu, Niels C. M. Martens, C. J. A. P. Martins, Konstantinos Migkas, Eoin Ó Colgáin, Pratyush Pranav, Lior Shamir, Ashok K. Singal, M. M. Sheikh-Jabbari, Jenny Wagner, Shao-Jiang Wang, David L. Wiltshire, Shek Yeung, Lu Yin, Wen Zhao

Abstract: The Cosmological Principle (CP) -- the notion that the Universe is spatially isotropic and homogeneous on large scales -- underlies a century of progress in cosmology. It is conventionally formulated through the Friedmann-Lemaître-Robertson-Walker (FLRW) cosmologies as the spacetime metric, and culminates in the successful and highly predictive $Λ$-Cold-Dark-Matter ($Λ$CDM) model. Yet, tensions ha… ▽ More The Cosmological Principle (CP) -- the notion that the Universe is spatially isotropic and homogeneous on large scales -- underlies a century of progress in cosmology. It is conventionally formulated through the Friedmann-Lemaître-Robertson-Walker (FLRW) cosmologies as the spacetime metric, and culminates in the successful and highly predictive $Λ$-Cold-Dark-Matter ($Λ$CDM) model. Yet, tensions have emerged within the $Λ$CDM model, most notably a statistically significant discrepancy in the value of the Hubble constant, $H_0$. Since the notion of cosmic expansion determined by a single parameter is intimately tied to the CP, implications of the $H_0$ tension may extend beyond $Λ$CDM to the CP itself. This review surveys current observational hints for deviations from the expectations of the CP, highlighting synergies and disagreements that warrant further study. Setting aside the debate about individual large structures, potential deviations from the CP include variations of cosmological parameters on the sky, discrepancies in the cosmic dipoles, and mysterious alignments in quasar polarizations and galaxy spins. While it is possible that a host of observational systematics are impacting results, it is equally plausible that precision cosmology may have outgrown the FLRW paradigm, an extremely pragmatic but non-fundamental symmetry assumption. △ Less

Submitted 27 February, 2023; v1 submitted 12 July, 2022; originally announced July 2022.

Comments: extended contents and references, 73 pages (excluding references), 30 figures, version accepted for publication in Class. Quant. Grav. "Focus issue on the Hubble constant tension"

Journal ref: Classical and Quantum Gravity, Vol. 40, Issue No. 9, Page No. 094001 (2023)

Validation of strategies for coupling exoplanet PSFs into single-mode fibres for high-dispersion coronagraphy

Authors: M. El Morsy, A. Vigan, M. Lopez, G. P. P. L. Otten, E. Choquet, F. Madec, A. Costille, J. -F. Sauvage, K. Dohlen, E. Muslimov, R. Pourcelot, J. Floriot, J. -A. Benedetti, P. Blanchard, P. Balard, G. Murray

Abstract: On large ground-based telescopes, the combination of extreme adaptive optics (ExAO) and coronagraphy with high-dispersion spectroscopy (HDS), sometimes referred to as high-dispersion coronagraphy (HDC), is starting to emerge as a powerful technique for the direct characterisation of giant exoplanets. The high spectral resolution not only brings a major gain in terms of accessible spectral features… ▽ More On large ground-based telescopes, the combination of extreme adaptive optics (ExAO) and coronagraphy with high-dispersion spectroscopy (HDS), sometimes referred to as high-dispersion coronagraphy (HDC), is starting to emerge as a powerful technique for the direct characterisation of giant exoplanets. The high spectral resolution not only brings a major gain in terms of accessible spectral features but also enables a better separation of the stellar and planetary signals. Ongoing projects such as Keck/KPIC, Subaru/REACH, and VLT/HiRISE base their observing strategy on the use of a few science fibres, one of which is dedicated to sampling the planet's signal, while the others sample the residual starlight in the speckle field. The main challenge in this approach is to blindly centre the planet's point spread function (PSF) accurately on the science fibre, with an accuracy of less than 0.1 $λ/D$ to maximise the coupling efficiency. In the context of the HiRISE project, three possible centring strategies are foreseen, either based on retro-injecting calibration fibres to localise the position of the science fibre or based on a dedicated centring fibre. We implemented these three approaches, and we compared their centring accuracy using an upgraded setup of the MITHiC high-contrast imaging testbed, which is similar to the setup that will be adopted in HiRISE. Our results demonstrate that reaching a specification accuracy of 0.1 $λ/D$ is extremely challenging regardless of the chosen centring strategy. It requires a high level of accuracy at every step of the centring procedure, which can be reached with very stable instruments. We studied the contributors to the centring error in the case of MITHiC and we propose a quantification for some of the most impacting terms. △ Less

Submitted 24 May, 2022; originally announced May 2022.

Journal ref: A&A 667, A171 (2022)

Simulating Transient Noise Bursts in LIGO with gengli

Authors: Melissa Lopez, Vincent Boudart, Stefano Schmidt, Sarah Caudill

Abstract: In the field of gravitational-wave (GW) interferometers, the most severe limitation to the detection of transient signals from astrophysical sources comes from transient noise artefacts, known as glitches, that happens at a rate around $1$ per minute. Because glitches reduce the amount of scientific data available, there is a need for better modelling and inclusion of glitches in large-scale studi… ▽ More In the field of gravitational-wave (GW) interferometers, the most severe limitation to the detection of transient signals from astrophysical sources comes from transient noise artefacts, known as glitches, that happens at a rate around $1$ per minute. Because glitches reduce the amount of scientific data available, there is a need for better modelling and inclusion of glitches in large-scale studies, such as stress testing the search pipelines and increasing the confidence of detection. In this work, we employ a Generative Adversarial Network (GAN) to produce a particular class of glitches ({\it blip}) in the time domain. We share the trained network through a user-friendly open-source software package called \texttt{gengli} and provide practical examples of its usage. △ Less

Submitted 10 November, 2022; v1 submitted 18 May, 2022; originally announced May 2022.

Comments: 6 pages, 4 figures

Simulating Transient Noise Bursts in LIGO with Generative Adversarial Networks

Authors: Melissa Lopez, Vincent Boudart, Kerwin Buijsman, Amit Reza, Sarah Caudill

Abstract: The noise of gravitational-wave (GW) interferometers limits their sensitivity and impacts the data quality, hindering the detection of GW signals from astrophysical sources. For transient searches, the most problematic are transient noise artifacts, known as glitches, that happen at a rate around $ 1 \text{ min}^{-1}$, and can mimic GW signals. Because of this, there is a need for better modeling… ▽ More The noise of gravitational-wave (GW) interferometers limits their sensitivity and impacts the data quality, hindering the detection of GW signals from astrophysical sources. For transient searches, the most problematic are transient noise artifacts, known as glitches, that happen at a rate around $ 1 \text{ min}^{-1}$, and can mimic GW signals. Because of this, there is a need for better modeling and inclusion of glitches in large-scale studies, such as stress testing the pipelines. In this proof-of concept work we employ Generative Adversarial Networks (GAN), a state-of-the-art Deep Learning algorithm inspired by Game Theory, to learn the underlying distribution of blip glitches and to generate artificial populations. We reconstruct the glitch in the time-domain, providing a smooth input that the GAN can learn. With this methodology, we can create distributions of $\sim 10^{3}$ glitches from Hanford and Livingston detectors in less than one second. Furthermore, we employ several metrics to measure the performance of our methodology and the quality of its generations. This investigation will be extended in the future to different glitch classes with the final goal of creating an open-source interface for mock data generation. △ Less

Submitted 21 March, 2022; v1 submitted 12 March, 2022; originally announced March 2022.

Comments: 17 pages, 21 figures to be published in Physical Review D

A Giant Arc on the Sky

Authors: Alexia M. Lopez, Roger G. Clowes, Gerard M. Williger

Abstract: We present the serendipitous discovery of a `Giant Arc on the Sky' at $z \sim 0.8$. The Giant Arc (GA) spans $\sim 1$ Gpc (proper size, present epoch), and appears to be almost symmetrical on the sky. It was discovered via intervening MgII absorbers in the spectra of background quasars, using the catalogues of Zhu \& Ménard. The use of MgII absorbers represents a new approach to the investigation… ▽ More We present the serendipitous discovery of a `Giant Arc on the Sky' at $z \sim 0.8$. The Giant Arc (GA) spans $\sim 1$ Gpc (proper size, present epoch), and appears to be almost symmetrical on the sky. It was discovered via intervening MgII absorbers in the spectra of background quasars, using the catalogues of Zhu \& Ménard. The use of MgII absorbers represents a new approach to the investigation of large-scale structures (LSSs) at redshifts $0.45 \lesssim z \lesssim 2.25$. We present the observational properties of the GA, and we assess it statistically using methods based on: (i) single-linkage hierarchical clustering ($\sim 4.5σ$); (ii) the Cuzick-Edwards test ($\sim 3.0σ$); and (iii) power spectrum analysis ($\sim 4.8σ$). Each of these methods has distinctive attributes and powers, and we advise considering the evidence from the ensemble. We discuss our approaches to mitigating any {\it post-hoc} aspects of analysing significance after discovery. The overdensity of the GA is $δρ/ ρ\sim 1.3 \pm 0.3$. The GA is the newest and one of the largest of a steadily accumulating set of very large LSSs that may (cautiously) challenge the Cosmological Principle, upon which the `standard model' of cosmology is founded. Conceivably, the GA is the precursor of a structure like the Sloan Great Wall (but the GA is about twice the size), seen when the Universe was about half its present age. △ Less

Submitted 2 August, 2022; v1 submitted 18 January, 2022; originally announced January 2022.

Comments: Accepted for MNRAS

A solar flare driven by thermal conduction observed in mid-infrared

Authors: Fernando M. López, C. Guillermo Giménez de Castro, Cristina H. Mandrini, Paulo J. A. Simões, Germán D. Cristiani, Dale E. Gary, Carlos Francile, Pascal Démoulin

Abstract: The mid-infrared (mid-IR) range has been mostly unexplored for the investigation of solar flares. It is only recently that new mid-IR flare observations have begun opening a new window into the response and evolution of the solar chromosphere. These new observations have been mostly performed by the AR30T and BR30T telescopes that are operating in Argentina and Brazil, respectively. We present the… ▽ More The mid-infrared (mid-IR) range has been mostly unexplored for the investigation of solar flares. It is only recently that new mid-IR flare observations have begun opening a new window into the response and evolution of the solar chromosphere. These new observations have been mostly performed by the AR30T and BR30T telescopes that are operating in Argentina and Brazil, respectively. We present the analysis of SOL2019-05-15T19:24, a GOES class C2.0 solar flare observed at 30~THz (10$\ μ$m) by the ground-based telescope AR30T. Our aim is to characterize the evolution of the flaring atmosphere and the energy transport mechanism in the context of mid-IR emission. We performed a multi-wavelength analysis of the event by complementing the mid-IR data with diverse ground- and space-based data from the Solar Dynamics Observatory (SDO), the H--$α$ Solar Telescope for Argentina (HASTA), and the Expanded Owens Valley Solar Array (EOVSA). Our study includes the analysis of the magnetic field evolution of the flaring region and of the development of the flare. The mid-IR images from AR30T show two bright and compact flare sources that are spatially associated with the flare kernels observed in ultraviolet (UV) by SDO. We confirm that the temporal association between mid-IR and UV fluxes previously reported for strong flares is also observed for this small flare. The EOVSA microwave data revealed flare spectra consistent with thermal free-free emission, which lead us to dismiss the existence of a significant number of non-thermal electrons. We thus consider thermal conduction as the primary mechanism responsible for energy transport. Our estimates for the thermal conduction energy and total radiated energy fall within the same order of magnitude, reinforcing our conclusions. △ Less

Submitted 29 October, 2021; originally announced October 2021.

Comments: 10 pages, 8 figures. Accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 657, A51 (2022)

MATISSE, the VLTI mid-infrared imaging spectro-interferometer

Authors: B. Lopez, S. Lagarde, R. G. Petrov, W. Jaffe, P. Antonelli, F. Allouche, P. Berio, A. Matter, A. Meilland, F. Millour, S. Robbe-Dubois, Th. Henning, G. Weigelt, A. Glindemann, T. Agocs, Ch. Bailet, U. Beckmann, F. Bettonvil, R. van Boekel, P. Bourget, Y. Bresson, P. Bristow, P. Cruzalèbes, E. Eldswijk, Y. Fanteï Caujolle , et al. (128 additional authors not shown)

Abstract: Context:Optical interferometry is at a key development stage. ESO's VLTI has established a stable, robust infrastructure for long-baseline interferometry for general astronomical observers. The present second-generation instruments offer a wide wavelength coverage and improved performance. Their sensitivity and measurement accuracy lead to data and images of high reliability. Aims:We have develope… ▽ More Context:Optical interferometry is at a key development stage. ESO's VLTI has established a stable, robust infrastructure for long-baseline interferometry for general astronomical observers. The present second-generation instruments offer a wide wavelength coverage and improved performance. Their sensitivity and measurement accuracy lead to data and images of high reliability. Aims:We have developed MATISSE, the Multi AperTure mid-Infrared SpectroScopic Experiment, to access high resolution imaging in a wide spectral domain and explore topics such: stellar activity and mass loss; planet formation and evolution in the gas and dust disks around young stars; accretion processes around super massive black holes in AGN. Methods:The instrument is a spectro-interferometric imager covering three atmospheric bands (L,M,N) from 2.8 to 13.0 mu, combining four optical beams from the VLTI's telscopes. Its concept, related observing procedure, data reduction and calibration approach are the product of 30 years of instrumental research. The instrument utilizes a multi-axial beam combination that delivers spectrally dispersed fringes. The signal provides the following quantities at several spectral resolutions: photometric flux, coherent fluxes, visibilities, closure phases, wavelength differential visibilities and phases, and aperture-synthesis imaging. Results:We provide an overview of the physical principle of the instrument and its functionalities, the characteristics of the delivered signal, a description of the observing modes and of their performance limits. An ensemble of data and reconstructed images are illustrating the first acquired key observations. Conclusion:The instrument has been in operation at Cerro Paranal, ESO, Chile since 2018, and has been open for science use by the international community since April 2019. The first scientific results are being published now. △ Less

Submitted 2 March, 2022; v1 submitted 29 October, 2021; originally announced October 2021.

Comments: 24 pages, 26 figures, submitted to Astronomy & Astrophysics

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

Chaotic Inflation and Reheating in Generalized Scalar-Tensor Gravity

Authors: Matías López, Giovanni Otalora, Nelson Videla

Abstract: In the present work, we study slow-roll inflation in scalar-tensor gravity theories in the presence of both the non-minimal coupling between the scalar field and curvature, and the Galileon self-interaction of the scalar field. Furthermore, we give predictions for the duration of reheating as well as for the reheating temperature after inflation. After working out the expressions for the power spe… ▽ More In the present work, we study slow-roll inflation in scalar-tensor gravity theories in the presence of both the non-minimal coupling between the scalar field and curvature, and the Galileon self-interaction of the scalar field. Furthermore, we give predictions for the duration of reheating as well as for the reheating temperature after inflation. After working out the expressions for the power spectra of scalar and tensor perturbations in the case of a general non-minimal coupling function that depends solely on the scalar field and a general scalar potential, we focus on the special cases of the power-law coupling function and chaotic quadratic inflation. Thus, under the slow-roll approximation we confront the predictions of the model with the current PLANCK constraints on the spectral index $n_s$ and the tensor-to-scalar ratio $r$ using the $n_{s}-r$ plane. We found that the combination of the non-minimal coupling and Galileon self-interaction effects allows us to obtain better results for $r$ than in the case in which each effect is considered separately. Particularly, we obtained that the predictions of the model are in agreement with the current observational bounds on $n_{s}$ and $r$ within the $95 \%$ C.L region and also slightly inside the $68 \%$ C.L region. Also, we investigate the oscillatory regime after the end of inflation by solving the full background equations, and then we determine the upper bound for the Galileon and non-minimal coupling parameters under the condition that the scalar field oscillates coherently during reheating. Finally, after approximating reheating by a constant equation of state, we derive the relations between the reheating duration, the temperature at the end of reheating, its equation of state, and the number of $e$-folds of inflation and then we relate them all to the inflationary observables. △ Less

Submitted 13 September, 2021; v1 submitted 15 July, 2021; originally announced July 2021.

Comments: 26 pages, 8 figures, 2 tables. Accepted version for publication in JCAP

Journal ref: JCAP10(2021)021

Searching for ring-like structures in the Cosmic Microwave Background

Authors: Melissa Lopez, Pietro Bonizzi, Kurt Driessens, Gideon Koekoek, Jacco de Vries, Ronald Westra

Abstract: In this research, we present an alternative methodology to search for ring-like structures in the sky with unusually large temperature gradients, namely Hawking points (HP), in the Cosmic Microwave Background (CMB), which are possible observational effects associated with Conformal Cyclic Cosmology (CCC). To assess the performance of our method, we constructed an artificial data set of HP, accordi… ▽ More In this research, we present an alternative methodology to search for ring-like structures in the sky with unusually large temperature gradients, namely Hawking points (HP), in the Cosmic Microwave Background (CMB), which are possible observational effects associated with Conformal Cyclic Cosmology (CCC). To assess the performance of our method, we constructed an artificial data set of HP, according to CCC, and we were able to retrieve $95 \%$ of ring-like anomalies from it. Furthermore, we scanned the \textit{Planck} CMB sky map and compared it to simulations according to $ΛCDM$, where we applied robust statistical tests to assess the existence of HP. Even though no significant ring-like structures were observed, we report the largest excess of HP candidates found at $α= $1\% significance level for the analyzed sky maps (CMB at 70GHz, SEVEM, SMICA, and Commander-Ruler), and we stress the need to continue the theoretical and experimental research in this direction. △ Less

Submitted 6 January, 2023; v1 submitted 9 May, 2021; originally announced May 2021.

Comments: 9 pages, 5 figures, 1 table, 1 algorithm

Mid-infrared circumstellar emission of the long-period Cepheid l Carinae resolved with VLTI/MATISSE

Authors: V. Hocdé, N. Nardetto, A. Matter, E. Lagadec, A. Mérand, P. Cruzalèbes, A. Meilland, F. Millour, B. Lopez, P. Berio, G. Weigelt, R. Petrov, J. W. Isbell, W. Jaffe, P. Kervella, A. Glindemann, M. Schöller, F. Allouche, A. Gallenne, A. Domiciano de Souza, G. Niccolini, E. Kokoulina, J. Varga, S. Lagarde, J. -C. Augereau , et al. (129 additional authors not shown)

Abstract: The nature of circumstellar envelopes (CSE) around Cepheids is still a matter of debate. The physical origin of their infrared (IR) excess could be either a shell of ionized gas, or a dust envelope, or both. This study aims at constraining the geometry and the IR excess of the environment of the long-period Cepheid $\ell$ Car (P=35.5 days) at mid-IR wavelengths to understand its physical nature. W… ▽ More The nature of circumstellar envelopes (CSE) around Cepheids is still a matter of debate. The physical origin of their infrared (IR) excess could be either a shell of ionized gas, or a dust envelope, or both. This study aims at constraining the geometry and the IR excess of the environment of the long-period Cepheid $\ell$ Car (P=35.5 days) at mid-IR wavelengths to understand its physical nature. We first use photometric observations in various bands and Spitzer Space Telescope spectroscopy to constrain the IR excess of $\ell$ Car. Then, we analyze the VLTI/MATISSE measurements at a specific phase of observation, in order to determine the flux contribution, the size and shape of the environment of the star in the L band. We finally test the hypothesis of a shell of ionized gas in order to model the IR excess. We report the first detection in the L band of a centro-symmetric extended emission around l Car, of about 1.7$R_\star$ in FWHM, producing an excess of about 7.0\% in this band. In the N band, there is no clear evidence for dust emission from VLTI/MATISSE correlated flux and Spitzer data. On the other side, the modeled shell of ionized gas implies a more compact CSE ($1.13\pm0.02\,R_\star$) and fainter (IR excess of 1\% in the L band). We provide new evidences for a compact CSE of $\ell$ Car and we demonstrate the capabilities of VLTI/MATISSE for determining common properties of CSEs. While the compact CSE of $\ell$ Car is probably of gaseous nature, the tested model of a shell of ionized gas is not able to simultaneously reproduce the IR excess and the interferometric observations. Further Galactic Cepheids observations with VLTI/MATISSE are necessary for determining the properties of CSEs, which may also depend on both the pulsation period and the evolutionary state of the stars. △ Less

Submitted 31 March, 2021; originally announced March 2021.

Comments: 13 pages, 8 figures, accepted in Astronomy and Astrophysics

Calibration of residual aberrations in exoplanet imagers with large numbers of degrees of freedom

Authors: Raphaël Pourcelot, Arthur Vigan, Kjetil Dohlen, Bastien Rouzé, Jean-François Sauvage, Mona El Morsy, Maxime Lopez, Mamadou N'Diaye, Amandine Caillat, Élodie Choquet, Gilles P. P. L. Otten, Alain Abbinanti, Philippe Balard, Marcel Carbillet, Philippe Blanchard, Jérémy Hulin, Émilie Robert

Abstract: Imaging faint objects, such as exoplanets or disks, around nearby stars is extremely challenging because host star images are dominated by the telescope diffraction pattern. Using a coronagraph is an efficient solution for removing diffraction but requires an incoming wavefront with good quality to maximize starlight rejection. On the ground, the most advanced exoplanet imagers use extreme adaptiv… ▽ More Imaging faint objects, such as exoplanets or disks, around nearby stars is extremely challenging because host star images are dominated by the telescope diffraction pattern. Using a coronagraph is an efficient solution for removing diffraction but requires an incoming wavefront with good quality to maximize starlight rejection. On the ground, the most advanced exoplanet imagers use extreme adaptive optics (ExAO) systems that are based on a deformable mirror (DM) with a large number of actuators to efficiently compensate for high-order aberrations and provide diffraction-limited images. While several exoplanet imagers with DMs using around 1500 actuators are now routinely operating on large telescopes to observe gas giant planets, future systems may require a tenfold increase in the number of degrees of freedom to look for rocky planets. In this paper, we explore wavefront correction with a secondary adaptive optics system that controls a very large number of degrees of freedom that are not corrected by the primary ExAO system. Using Marseille Imaging Testbed for High Contrast (MITHiC), we implement a second stage of adaptive optics with ZELDA, a Zernike wavefront sensor, and a spatial light modulator (SLM) to compensate for the phase aberrations of the bench downstream residual aberrations from adaptive optics. We demonstrate that their correction up to 137 cycles per pupil with nanometric accuracy is possible, provided there is a simple distortion calibration of the pupil and a moderate wavefront low-pass filtering. We also use ZELDA for a fast compensation of ExAO residuals, showing its promising implementation as a second-stage correction for the observation of rocky planets around nearby stars. Finally, we present images with a classical Lyot coronagraph on MITHiC and validate our ability to reach its theoretical performance with our calibration. △ Less

Submitted 29 March, 2021; originally announced March 2021.

Journal ref: A&A 649, A170 (2021)

Observation of a sudden cessation of a very-high-energy gamma-ray flare in PKS 1510-089 with H.E.S.S. and MAGIC in May 2016

Authors: H. E. S. S. Collaboration, H. Abdalla, R. Adam, F. Aharonian, F. Ait Benkhali, E. O. Angüner, C. Arcaro, C. Arm, T. Armstrong, H. Ashkar, M. Backes, V. Baghmanyan, V. Barbosa Martins, A. Barnacka, M. Barnard, Y. Becherini, D. Berge, K. Bernlöhr, B. Bi, M. Böttcher, C. Boisson, J. Bolmont, S. Bonnefoy, M. de Bony de Lavergne, J. Bregeon , et al. (409 additional authors not shown)

Abstract: The flat spectrum radio quasar (FSRQ) PKS 1510-089 is known for its complex multiwavelength behavior, and is one of only a few FSRQs detected at very high energy (VHE, $E>100\,$GeV) $γ$-rays. VHE $γ$-ray observations with H.E.S.S. and MAGIC during late May and early June 2016 resulted in the detection of an unprecedented flare, which reveals for the first time VHE $γ$-ray intranight variability in… ▽ More The flat spectrum radio quasar (FSRQ) PKS 1510-089 is known for its complex multiwavelength behavior, and is one of only a few FSRQs detected at very high energy (VHE, $E>100\,$GeV) $γ$-rays. VHE $γ$-ray observations with H.E.S.S. and MAGIC during late May and early June 2016 resulted in the detection of an unprecedented flare, which reveals for the first time VHE $γ$-ray intranight variability in this source. While a common variability timescale of $1.5\,$hr is found, there is a significant deviation near the end of the flare with a timescale of $\sim 20\,$min marking the cessation of the event. The peak flux is nearly two orders of magnitude above the low-level emission. For the first time, curvature is detected in the VHE $γ$-ray spectrum of PKS 1510-089, which is fully explained through absorption by the extragalactic background light. Optical R-band observations with ATOM reveal a counterpart of the $γ$-ray flare, even though the detailed flux evolution differs from the VHE ightcurve. Interestingly, a steep flux decrease is observed at the same time as the cessation of the VHE flare. In the high energy (HE, $E>100\,$MeV) $γ$-ray band only a moderate flux increase is observed with Fermi-LAT, while the HE $γ$-ray spectrum significantly hardens up to a photon index of 1.6. A search for broad-line region (BLR) absorption features in the $γ$-ray spectrum indicates that the emission region is located outside of the BLR. Radio VLBI observations reveal a fast moving knot interacting with a standing jet feature around the time of the flare. As the standing feature is located $\sim 50\,$pc from the black hole, the emission region of the flare may have been located at a significant distance from the black hole. If this correlation is indeed true, VHE $γ$ rays have been produced far down the jet where turbulent plasma crosses a standing shock. △ Less

Submitted 18 December, 2020; originally announced December 2020.

Comments: 25 pages, 17 figures. Accepted for publication by Astronomy & Astrophysics. Corresponding authors: M. Zacharias, J. Sitarek, D. Sanchez, T. Terzic

Journal ref: A&A 648, A23 (2021)

The asymmetric inner disk of the Herbig Ae star HD 163296 in the eyes of VLTI/MATISSE: evidence for a vortex?

Authors: J. Varga, M. Hogerheijde, R. van Boekel, L. Klarmann, R. Petrov, L. B. F. M. Waters, S. Lagarde, E. Pantin, Ph. Berio, G. Weigelt, S. Robbe-Dubois, B. Lopez, F. Millour, J. -C. Augereau, H. Meheut, A. Meilland, Th. Henning, W. Jaffe, F. Bettonvil, P. Bristow, K. -H. Hofmann, A. Matter, G. Zins, S. Wolf, F. Allouche , et al. (111 additional authors not shown)

Abstract: Context. The inner few au region of planet-forming disks is a complex environment. High angular resolution observations have a key role in understanding the disk structure and the dynamical processes at work. Aims. In this study we aim to characterize the mid-infrared brightness distribution of the inner disk of the young intermediate-mass star HD 163296, from VLTI/MATISSE observations. Methods. W… ▽ More Context. The inner few au region of planet-forming disks is a complex environment. High angular resolution observations have a key role in understanding the disk structure and the dynamical processes at work. Aims. In this study we aim to characterize the mid-infrared brightness distribution of the inner disk of the young intermediate-mass star HD 163296, from VLTI/MATISSE observations. Methods. We use geometric models to fit the data. Our models include a smoothed ring, a flat disk with inner cavity, and a 2D Gaussian. The models can account for disk inclination and for azimuthal asymmetries as well. We also perform numerical hydro-dynamical simulations of the inner edge of the disk. Results. Our modeling reveals a significant brightness asymmetry in the L-band disk emission. The brightness maximum of the asymmetry is located at the NW part of the disk image, nearly at the position angle of the semimajor axis. The surface brightness ratio in the azimuthal variation is $3.5 \pm 0.2$. Comparing our result on the location of the asymmetry with other interferometric measurements, we confirm that the morphology of the $r<0.3$ au disk region is time-variable. We propose that this asymmetric structure, located in or near the inner rim of the dusty disk, orbits the star. For the physical origin of the asymmetry, we tested a hypothesis where a vortex is created by Rossby wave instability, and we find that a unique large scale vortex may be compatible with our data. The half-light radius of the L-band emitting region is $0.33\pm 0.01$ au, the inclination is ${52^\circ}^{+5^\circ}_{-7^\circ}$, and the position angle is $143^\circ \pm 3^\circ$. Our models predict that a non-negligible fraction of the L-band disk emission originates inside the dust sublimation radius for $μ$m-sized grains. Refractory grains or large ($\gtrsim 10\ μ$m-sized) grains could be the origin for this emission. △ Less

Submitted 10 December, 2020; originally announced December 2020.

Comments: accepted for publication in A&A

Journal ref: A&A 647, A56 (2021)

Multiwavelength variability and correlation studies of Mrk 421 during historically low X-ray and $γ$-ray activity in 2015$-$2016

Authors: MAGIC Collaboration, V. A. Acciari, S. Ansoldi, L. A. Antonelli, K. Asano, A. Babić, B. Banerjee, A. Baquero, U. Barres de Almeida, J. A. Barrio, J. Becerra González, W. Bednarek, L. Bellizzi, E. Bernardini, M. Bernardos, A. Berti, J. Besenrieder, W. Bhattacharyya, C. Bigongiari, O. Blanch, G. Bonnoli, Ž. Bošnjak, G. Busetto, R. Carosi, G. Ceribella , et al. (205 additional authors not shown)

Abstract: We report a characterization of the multi-band flux variability and correlations of the nearby (z=0.031) blazar Markarian 421 (Mrk 421) using data from Metsähovi, Swift, Fermi-LAT, MAGIC, FACT and other collaborations and instruments from November 2014 till June 2016. Mrk 421 did not show any prominent flaring activity, but exhibited periods of historically low activity above 1 TeV (F… ▽ More We report a characterization of the multi-band flux variability and correlations of the nearby (z=0.031) blazar Markarian 421 (Mrk 421) using data from Metsähovi, Swift, Fermi-LAT, MAGIC, FACT and other collaborations and instruments from November 2014 till June 2016. Mrk 421 did not show any prominent flaring activity, but exhibited periods of historically low activity above 1 TeV (F$_{>1\mathrm{TeV}}<$ 1.7$\times$10$^{-12}$ ph cm$^{-2}$ s$^{-1}$) and in the 2-10 keV (X-ray) band (F$_{2-10 \mathrm{keV}}<$3.6$\times$10$^{-11}$ erg cm$^{-2}$ s$^{-1}$), during which the Swift-BAT data suggests an additional spectral component beyond the regular synchrotron emission. The highest flux variability occurs in X-rays and very-high-energy (E$>$0.1 TeV) $γ$-rays, which, despite the low activity, show a significant positive correlation with no time lag. The HR$_\mathrm{keV}$ and HR$_\mathrm{TeV}$ show the harder-when-brighter trend observed in many blazars, but the trend flattens at the highest fluxes, which suggests a change in the processes dominating the blazar variability. Enlarging our data set with data from years 2007 to 2014, we measured a positive correlation between the optical and the GeV emission over a range of about 60 days centered at time lag zero, and a positive correlation between the optical/GeV and the radio emission over a range of about 60 days centered at a time lag of $43^{+9}_{-6}$ days.This observation is consistent with the radio-bright zone being located about 0.2 parsec downstream from the optical/GeV emission regions of the jet. The flux distributions are better described with a LogNormal function in most of the energy bands probed, indicating that the variability in Mrk 421 is likely produced by a multiplicative process. △ Less

Submitted 2 December, 2020; originally announced December 2020.

Comments: 45 pages, 30 figures, 6 tables, accepted for publication in MNRAS

On the collisional disalignment of dust grains in illuminated and shaded regions of IC 63

Authors: Archana Soam, B-G Andersson, Jose Acosta-Pulido, Manuel Fernández López, John E. Vaillancourt, Susanna L. Widicus Weaver, Vilppu Piirola, Michael S. Gordon

Abstract: Interstellar dust grain alignment causes polarization from UV to mm wavelengths, allowing the study of the geometry and strength of the magnetic field. Over last couple of decades observations and theory have led to the establishment of the Radiative Alignment Torque (RAT) mechanism as leading candidate to explain the effect. With a quantitatively well constrained theory, polarization can be used… ▽ More Interstellar dust grain alignment causes polarization from UV to mm wavelengths, allowing the study of the geometry and strength of the magnetic field. Over last couple of decades observations and theory have led to the establishment of the Radiative Alignment Torque (RAT) mechanism as leading candidate to explain the effect. With a quantitatively well constrained theory, polarization can be used not only to study the interstellar magnetic field, but also the dust and other environmental parameters. Photo-dissociation Regions (PDRs), with their intense, anisotropic radiation fields, consequent rapid $\rm H_{2}$ formation, and high spatial density-contrast provide a rich environment for such studies. Here we discuss an expanded optical, NIR, and mm-wave study of the IC\,63 nebula, showing strong $\rm H_{2}$ formation-enhanced alignment and the first direct empirical evidence for disalignment due to gas-grain collisions using high-resolution $\rm HCO^{+}$(J=1-0) observations. We find that relative amount of polarization is marginally anti-correlated with column density of $\rm HCO^{+}$. However, separating the lines of sight of optical polarimetry into those behind, or in front of, a dense clump as seen from $γ$ Cas, the distribution separates into two well defined sets, with data corresponding to \enquote{shaded} gas having a shallower slope. This is expected if the decrease in polarization is caused by collisions since collisional disalignment rate is proportional to R$_C\propto n\sqrt{T}$. Ratios of the best-fit slopes for the \enquote{illuminated} and \enquote{shaded} samples of lines of sight agrees, within the uncertainties, with the square-root of the two-temperature H$_2$ excitation in the nebula seen by Thi et al. (2009). △ Less

Submitted 18 November, 2020; originally announced November 2020.

Comments: 14 pages, 6 figures, Accepted for publication in ApJ

The dual nature of blazar fast variability. Space and ground observations of S5 0716+714

Authors: C. M. Raiteri, M. Villata, D. Carosati, E. Benítez, S. O. Kurtanidze, A. C. Gupta, D. O. Mirzaqulov, F. D'Ammando, V. M. Larionov, T. Pursimo, J. A. Acosta-Pulido, G. V. Baida, B. Balmaverde, G. Bonnoli, G. A. Borman, M. I. Carnerero, W. -P. Chen, V. Dhiman, A. Di Maggio, S. A. Ehgamberdiev, D. Hiriart, G. N. Kimeridze, O. M. Kurtanidze, C. S. Lin, J. M. Lopez , et al. (17 additional authors not shown)

Abstract: Blazar S5 0716+714 is well-known for its short-term variability, down to intra-day time-scales. We here present the 2-min cadence optical light curve obtained by the TESS space telescope in 2019 December - 2020 January and analyse the object fast variability with unprecedented sampling. Supporting observations by the Whole Earth Blazar Telescope Collaboration in B, V , R, and I bands allow us to i… ▽ More Blazar S5 0716+714 is well-known for its short-term variability, down to intra-day time-scales. We here present the 2-min cadence optical light curve obtained by the TESS space telescope in 2019 December - 2020 January and analyse the object fast variability with unprecedented sampling. Supporting observations by the Whole Earth Blazar Telescope Collaboration in B, V , R, and I bands allow us to investigate the spectral variability during the TESS pointing. The spectral analysis is further extended in frequency to the UV and X-ray bands with data from the Neil Gehrels Swift Observatory. We develop a new method to unveil the shortest optical variability time-scales. This is based on progressive de-trending of the TESS light curve by means of cubic spline interpolations through the binned fluxes, with decreasing time bins. The de-trended light curves are then analysed with classical tools for time-series analysis (periodogram, auto-correlation and structure functions). The results show that below 3 d there are significant characteristic variability time-scales of about 1.7, 0.5, and 0.2 d. Variability on time-scales < 0.2 d is strongly chromatic and must be ascribed to intrinsic energetic processes involving emitting regions, likely jet sub-structures, with dimension less than about milliparsec. In contrast, flux changes on time-scales > 0.5 d are quasi-achromatic and are probably due to Doppler factor changes of geometric origin. △ Less

Submitted 12 November, 2020; originally announced November 2020.

Comments: 17 pages, 18 figures, accepted by MNRAS

Study of the GeV to TeV morphology of the $γ$-Cygni SNR (G78.2+2.1) with MAGIC and Fermi-LAT

Authors: MAGIC Collaboration, V. A. Acciari, S. Ansoldi, L. A. Antonelli, A. Arbet Engels, D. Baack, A. Babić, B. Banerjee, U. Barres de Almeida, J. A. Barrio, J. Becerra González, W. Bednarek, L. Bellizzi, E. Bernardini, A. Berti, J. Besenrieder, W. Bhattacharyya, C. Bigongiari, A. Biland, O. Blanch, G. Bonnoli, Ž. Bošnjak, G. Busetto, R. Carosi, G. Ceribella , et al. (161 additional authors not shown)

Abstract: Context. Diffusive shock acceleration (DSA) is the most promising mechanism to accelerate Galactic cosmic rays (CRs) in the shocks of supernova remnants (SNRs). The turbulence upstream is supposedly generated by the CRs, but this process is not well understood. The dominant mechanism may depend on the evolutionary state of the shock and can be studied via the CRs escaping upstream into the interst… ▽ More Context. Diffusive shock acceleration (DSA) is the most promising mechanism to accelerate Galactic cosmic rays (CRs) in the shocks of supernova remnants (SNRs). The turbulence upstream is supposedly generated by the CRs, but this process is not well understood. The dominant mechanism may depend on the evolutionary state of the shock and can be studied via the CRs escaping upstream into the interstellar medium (ISM). Aims. Previous observations of the $γ$-Cygni SNR showed a difference in morphology between GeV and TeV energies. Since this SNR has the right age and is at the evolutionary stage for a significant fraction of CRs to escape, we aim to understand $γ$-ray emission in the vicinity of the $γ$-Cygni SNR. Methods. We observed the region of the $γ$-Cygni SNR with the MAGIC Imaging Atmospheric Cherenkov telescopes between May 2015 and September 2017 recording 87 h of good-quality data. Additionally we analysed Fermi-LAT data to study the energy dependence of the morphology as well as the energy spectrum in the GeV to TeV range. The energy spectra and morphology were compared against theoretical predictions, which include a detailed derivation of the CR escape process and their $γ$-ray generation. Results. The MAGIC and Fermi-LAT data allowed us to identify three emission regions, which can be associated with the SNR and dominate at different energies. Our hadronic emission model accounts well for the morphology and energy spectrum of all source components. It constrains the time-dependence of the maximum energy of the CRs at the shock, the time-dependence of the level of turbulence, and the diffusion coefficient immediately outside the SNR shock. While in agreement with the standard picture of DSA, the time-dependence of the maximum energy was found to be steeper than predicted and the level of turbulence was found to change over the lifetime of the SNR. △ Less

Submitted 29 October, 2020; originally announced October 2020.

Comments: 20 pages, 10 figures. Accepted for publication in A&A. Corresponding authors: M. Strzys, G. Morlino, S. Masuda, and I. Vovk

Journal ref: A&A 670, A8 (2023)

Sensitivity of the Cherenkov Telescope Array for probing cosmology and fundamental physics with gamma-ray propagation

Authors: The Cherenkov Telescope Array Consortium, :, H. Abdalla, H. Abe, F. Acero, A. Acharyya, R. Adam, I. Agudo, A. Aguirre-Santaella, R. Alfaro, J. Alfaro, C. Alispach, R. Aloisio, R. Alves B, L. Amati, E. Amato, G. Ambrosi, E. O. Angüner, A. Araudo, T. Armstrong, F. Arqueros, L. Arrabito, K. Asano, Y. Ascasíbar, M. Ashley , et al. (474 additional authors not shown)

Abstract: The Cherenkov Telescope Array (CTA), the new-generation ground-based observatory for $γ$-ray astronomy, provides unique capabilities to address significant open questions in astrophysics, cosmology, and fundamental physics. We study some of the salient areas of $γ$-ray cosmology that can be explored as part of the Key Science Projects of CTA, through simulated observations of active galactic nucle… ▽ More The Cherenkov Telescope Array (CTA), the new-generation ground-based observatory for $γ$-ray astronomy, provides unique capabilities to address significant open questions in astrophysics, cosmology, and fundamental physics. We study some of the salient areas of $γ$-ray cosmology that can be explored as part of the Key Science Projects of CTA, through simulated observations of active galactic nuclei (AGN) and of their relativistic jets. Observations of AGN with CTA will enable a measurement of $γ$-ray absorption on the extragalactic background light with a statistical uncertainty below 15% up to a redshift $z=2$ and to constrain or detect $γ$-ray halos up to intergalactic-magnetic-field strengths of at least 0.3pG. Extragalactic observations with CTA also show promising potential to probe physics beyond the Standard Model. The best limits on Lorentz invariance violation from $γ$-ray astronomy will be improved by a factor of at least two to three. CTA will also probe the parameter space in which axion-like particles could constitute a significant fraction, if not all, of dark matter. We conclude on the synergies between CTA and other upcoming facilities that will foster the growth of $γ$-ray cosmology. △ Less

Submitted 26 February, 2021; v1 submitted 3 October, 2020; originally announced October 2020.

Comments: 71 pages (including affiliations and references), 13 figures, 6 tables. Accepted in JCAP; matches published version. Corresponding authors: Jonathan Biteau, Julien Lefaucheur, Humberto Martinez-Huerta, Manuel Meyer, Santiago Pita, Ievgen Vovk

Journal ref: JCAP 02 (2021) 048

Direct characterization of young giant exoplanets at high spectral resolution by coupling SPHERE and CRIRES+

Authors: G. P. P. L. Otten, A. Vigan, E. Muslimov, M. N'Diaye, E. Choquet, U. Seemann, K. Dohlen, M. Houllé, P. Cristofari, M. W. Phillips, Y. Charles, I. Baraffe, J. -L. Beuzit, A. Costille, R. Dorn, M. El Morsy, M. Kasper, M. Lopez, C. Mordasini, R. Pourcelot, A. Reiners, J. -F. Sauvage

Abstract: Studies of atmospheres of directly imaged exoplanets with high-resolution spectrographs have shown that their characterization is predominantly limited by noise on the stellar halo at the location of the studied exoplanet. An instrumental combination of high-contrast imaging and high spectral resolution that suppresses this noise and resolves the spectral lines can therefore yield higher quality s… ▽ More Studies of atmospheres of directly imaged exoplanets with high-resolution spectrographs have shown that their characterization is predominantly limited by noise on the stellar halo at the location of the studied exoplanet. An instrumental combination of high-contrast imaging and high spectral resolution that suppresses this noise and resolves the spectral lines can therefore yield higher quality spectra. We study the performance of the proposed HiRISE fiber coupling between the SPHERE and CRIRES+ at the VLT for spectral characterization of directly imaged planets. Using end-to-end simulations of HiRISE we determine the S/N of the detection of molecular species for known exoplanets in $H$ and $K$ bands, and compare them to CRIRES+. We investigate the ultimate detection limits of HiRISE as a function of stellar magnitude, and we quantify the impact of different coronagraphs and of the system transmission. We find that HiRISE largely outperforms CRIRES+ for companions around bright hosts like $β$ Pic or 51 Eri. For an $H=3.5$ host, we observe a gain of a factor of up to 16 in observing time with HiRISE to reach the same S/N on a companion at 200 mas. More generally, HiRISE provides better performance than CRIRES+ in two-hour integration times between 50-350 mas for hosts with $H<8.5$ and between 50-700 mas for $H<7$. For fainter hosts like PDS 70 and HIP 65426, no significant improvements are observed. We find that using no coronagraph yields the best S/N when characterizing known exoplanets due to higher transmission and fiber-based starlight suppression. We demonstrate that the overall transmission of the system is in fact the main driver of performance. Finally, we show that HiRISE outperforms the best detection limits of SPHERE for bright stars, opening major possibilities for the characterization of future planetary companions detected by other techniques. △ Less

Submitted 10 December, 2020; v1 submitted 3 September, 2020; originally announced September 2020.

Comments: 21 pages, 17 figures. 2 tables. Originally accepted for publication in A&A. Revised version after error in sky background was discovered (too high) and corrected in proofing stage. Revisions made in Abstract, Sections 3-6 and Figures 6-12,C1,C2 that address the changes in performance. Resubmitted to A&A

Sensitivity of the Cherenkov Telescope Array to a dark matter signal from the Galactic centre

Authors: The Cherenkov Telescope Array Consortium, :, A. Acharyya, R. Adam, C. Adams, I. Agudo, A. Aguirre-Santaella, R. Alfaro, J. Alfaro, C. Alispach, R. Aloisio, R. Alves Batista, L. Amati, G. Ambrosi, E. O. Angüner, L. A. Antonelli, C. Aramo, A. Araudo, T. Armstrong, F. Arqueros, K. Asano, Y. Ascasíbar, M. Ashley, C. Balazs, O. Ballester , et al. (427 additional authors not shown)

Abstract: We provide an updated assessment of the power of the Cherenkov Telescope Array (CTA) to search for thermally produced dark matter at the TeV scale, via the associated gamma-ray signal from pair-annihilating dark matter particles in the region around the Galactic centre. We find that CTA will open a new window of discovery potential, significantly extending the range of robustly testable models giv… ▽ More We provide an updated assessment of the power of the Cherenkov Telescope Array (CTA) to search for thermally produced dark matter at the TeV scale, via the associated gamma-ray signal from pair-annihilating dark matter particles in the region around the Galactic centre. We find that CTA will open a new window of discovery potential, significantly extending the range of robustly testable models given a standard cuspy profile of the dark matter density distribution. Importantly, even for a cored profile, the projected sensitivity of CTA will be sufficient to probe various well-motivated models of thermally produced dark matter at the TeV scale. This is due to CTA's unprecedented sensitivity, angular and energy resolutions, and the planned observational strategy. The survey of the inner Galaxy will cover a much larger region than corresponding previous observational campaigns with imaging atmospheric Cherenkov telescopes. CTA will map with unprecedented precision the large-scale diffuse emission in high-energy gamma rays, constituting a background for dark matter searches for which we adopt state-of-the-art models based on current data. Throughout our analysis, we use up-to-date event reconstruction Monte Carlo tools developed by the CTA consortium, and pay special attention to quantifying the level of instrumental systematic uncertainties, as well as background template systematic errors, required to probe thermally produced dark matter at these energies. "Full likelihood tables complementing our analysis are provided here [ https://doi.org/10.5281/zenodo.4057987 ]" △ Less

Submitted 30 January, 2021; v1 submitted 31 July, 2020; originally announced July 2020.

Comments: 68 pages (including references) and 26 figures; text identical to the version published in JCAP

Journal ref: JCAP01(2021)057

Studying the nature of the unidentified gamma-ray source HESS J1841-055 with the MAGIC telescopes

Authors: MAGIC Collaboration, V. A. Acciari, S. Ansoldi, L. A. Antonelli, A. Arbet Engels, K. Asano, D. Baack, A. Babić, A. Baquero, U. Barres de Almeida, J. A. Barrio, J. Becerra González, W. Bednarek, L. Bellizzi, E. Bernardini, M. Bernardos, A. Berti, J. Besenrieder, W. Bhattacharyya, C. Bigongiari, A. Biland, O. Blanch, G. Bonnoli, Ž. Bošnjak, G. Busetto , et al. (162 additional authors not shown)

Abstract: We investigate the physical nature and origin of the gamma-ray emission from the extended source HESS J1841-055 observed at TeV and GeV energies. We observed HESS J1841-055 at TeV energies for a total effective time of 43 hours with the MAGIC telescopes, in 2012 and 2013. Additionally, we analysed the GeV counterpart making use of about 10 years of Fermi-LAT data. Using both Fermi-LAT and MAGIC, w… ▽ More We investigate the physical nature and origin of the gamma-ray emission from the extended source HESS J1841-055 observed at TeV and GeV energies. We observed HESS J1841-055 at TeV energies for a total effective time of 43 hours with the MAGIC telescopes, in 2012 and 2013. Additionally, we analysed the GeV counterpart making use of about 10 years of Fermi-LAT data. Using both Fermi-LAT and MAGIC, we study both the spectral and energy-dependent morphology of the source for almost four decades of energy. The origin of the gamma-ray emission from this region is investigated using multi-waveband information on sources present in this region, suggested to be associated with this unidentified gamma-ray source. We find that the extended emission at GeV-TeV energies is best described by more than one source model. We also perform the first energy-dependent analysis of the HESS J1841-055 region at GeV-TeV. We find that the emission at lower energies comes from a diffuse or extended component, while the major contribution of gamma rays above 1 TeV arises from the southern part of the source. Moreover, we find that a significant curvature is present in the combined observed spectrum of MAGIC and Fermi-LAT. The first multi-wavelength spectral energy distribution of this unidentified source shows that the emission at GeV-TeV energies can be well explained with both leptonic and hadronic models. For the leptonic scenario, bremsstrahlung is the dominant emission compared to inverse Compton. On the other hand, for the hadronic model, gamma-ray resulting from the decay of neutral pions ($π^0$) can explain the observed spectrum. The presence of dense molecular clouds overlapping with HESS J1841-055 makes both bremsstrahlung and $π^0$-decay processes the dominant emission mechanisms for the source. △ Less

Submitted 17 July, 2020; originally announced July 2020.

Comments: 11 pages, 8 figures, 5 tables. Accepted in MNRAS

Observation of inverse Compton emission from a long $γ$-ray burst

Authors: V. A. Acciari, S. Ansoldi, L. A. Antonelli, A. Arbet Engels, D. Baack, A. Babić, B. Banerjee, U. Barres de Almeida, J. A. Barrio, J. Becerra González, W. Bednarek, L. Bellizzi, E. Bernardini, A. Berti, J. Besenrieder, W. Bhattacharyya, C. Bigongiari, A. Biland, O. Blanch, G. Bonnoli, Ž. Bošnjak, G. Busetto, R. Carosi, G. Ceribella, Y. Chai , et al. (279 additional authors not shown)

Abstract: Long-duration gamma-ray bursts (GRBs) originate from ultra-relativistic jets launched from the collapsing cores of dying massive stars. They are characterised by an initial phase of bright and highly variable radiation in the keV-MeV band that is likely produced within the jet and lasts from milliseconds to minutes, known as the prompt emission. Subsequently, the interaction of the jet with the ex… ▽ More Long-duration gamma-ray bursts (GRBs) originate from ultra-relativistic jets launched from the collapsing cores of dying massive stars. They are characterised by an initial phase of bright and highly variable radiation in the keV-MeV band that is likely produced within the jet and lasts from milliseconds to minutes, known as the prompt emission. Subsequently, the interaction of the jet with the external medium generates external shock waves, responsible for the afterglow emission, which lasts from days to months, and occurs over a broad energy range, from the radio to the GeV bands. The afterglow emission is generally well explained as synchrotron radiation by electrons accelerated at the external shock. Recently, an intense, long-lasting emission between 0.2 and 1 TeV was observed from the GRB 190114C. Here we present the results of our multi-frequency observational campaign of GRB~190114C, and study the evolution in time of the GRB emission across 17 orders of magnitude in energy, from $5\times10^{-6}$ up to $10^{12}$\,eV. We find that the broadband spectral energy distribution is double-peaked, with the TeV emission constituting a distinct spectral component that has power comparable to the synchrotron component. This component is associated with the afterglow, and is satisfactorily explained by inverse Compton upscattering of synchrotron photons by high-energy electrons. We find that the conditions required to account for the observed TeV component are not atypical, supporting the possibility that inverse Compton emission is commonly produced in GRBs. △ Less

Submitted 12 June, 2020; originally announced June 2020.

Journal ref: Nature 575 (2019) 459-463

Teraelectronvolt emission from the $γ$-ray burst GRB 190114C

Authors: MAGIC Collaboration, V. A. Acciari, S. Ansoldi, L. A. Antonelli, A. Arbet Engels, D. Baack, A. Babić, B. Banerjee, U. Barres de Almeida, J. A. Barrio, J. Becerra González, W. Bednarek, L. Bellizzi, E. Bernardini, A. Berti, J. Besenrieder, W. Bhattacharyya, C. Bigongiari, A. Biland, O. Blanch, G. Bonnoli, Ž. Bošnjak, G. Busetto, A. Carosi, R. Carosi , et al. (155 additional authors not shown)

Abstract: Gamma-ray bursts (GRBs) of the long-duration class are the most luminous sources of electromagnetic radiation known in the Universe. They are generated by outflows of plasma ejected at near the speed of light by newly formed neutron stars or black holes of stellar mass at cosmological distances. Prompt flashes of MeV gamma rays are followed by longer-lasting afterglow emission from radio waves to… ▽ More Gamma-ray bursts (GRBs) of the long-duration class are the most luminous sources of electromagnetic radiation known in the Universe. They are generated by outflows of plasma ejected at near the speed of light by newly formed neutron stars or black holes of stellar mass at cosmological distances. Prompt flashes of MeV gamma rays are followed by longer-lasting afterglow emission from radio waves to GeV gamma rays, due to synchrotron radiation by energetic electrons in accompanying shock waves. Although emission of gamma rays at even higher, TeV energies by other radiation mechanisms had been theoretically predicted, it had never been detected previously. Here we report the clear detection of GRB 190114C in the TeV band, achieved after many years of dedicated searches for TeV emission from GRBs. Gamma rays in the energy range 0.2--1 TeV are observed from about 1 minute after the burst (at more than 50 standard deviations in the first 20 minutes). This unambiguously reveals a new emission component in the afterglow of a GRB, whose power is comparable to that of the synchrotron component. The observed similarity in the radiated power and temporal behaviour of the TeV and X-ray bands points to processes such as inverse Compton radiation as the mechanism of the TeV emission, while processes such as synchrotron emission by ultrahigh-energy protons are disfavoured due to their low radiative efficiency. △ Less

Submitted 12 June, 2020; originally announced June 2020.

Journal ref: Nature 575 (2019) 455-458

An intermittent extreme BL Lac: MWL study of 1ES 2344+514 in an enhanced state

Authors: MAGIC Collaboration, V. A. Acciari, S. Ansoldi, L. A. Antonelli, A. Arbet Engels, D. Baack, A. Babić, B. Banerjee, U. Barres de Almeida, J. A. Barrio, J. Becerra González, W. Bednarek, L. Bellizzi, E. Bernardini, A. Berti, J. Besenrieder, W. Bhattacharyya, C. Bigongiari, O. Blanch, G. Bonnoli, Ž. Bošnjak, G. Busetto, R. Carosi, G. Ceribella, M. Cerruti , et al. (196 additional authors not shown)

Abstract: Extreme High-frequency BL~Lacs (EHBL) feature their synchrotron peak of the broadband spectral energy distribution (SED) at $ν_{\rm s} \geq $10$^{17}$\,Hz. The BL~Lac object 1ES~2344+514 was included in the EHBL family because of its impressive shift of the synchrotron peak in 1996. During the following years, the source appeared to be in a low state without showing any extreme behaviours. In Augu… ▽ More Extreme High-frequency BL~Lacs (EHBL) feature their synchrotron peak of the broadband spectral energy distribution (SED) at $ν_{\rm s} \geq $10$^{17}$\,Hz. The BL~Lac object 1ES~2344+514 was included in the EHBL family because of its impressive shift of the synchrotron peak in 1996. During the following years, the source appeared to be in a low state without showing any extreme behaviours. In August 2016, 1ES~2344+514 was detected with the ground-based $γ$-ray telescope FACT during a high $γ$-ray state, triggering multi-wavelength (MWL) observations. We studied the MWL light curves of 1ES~2344+514 during the 2016 flaring state, using data from radio to VHE $γ$ rays taken with OVRO, KAIT, KVA, NOT, some telescopes of the GASP-WEBT collaboration at the Teide, Crimean, and St. Petersburg observatories, \textit{Swift}-UVOT, \textit{Swift}-XRT, \textit{Fermi}-LAT, FACT and MAGIC. With simultaneous observations of the flare, we built the broadband SED and studied it in the framework of a leptonic and an hadronic model. The VHE $γ$-ray observations show a flux level of 55\% of the Crab Nebula flux above 300\,GeV, similar to the historical maximum of 1995. The combination of MAGIC and \textit{Fermi}-LAT spectra provides an unprecedented characterization of the inverse-Compton peak for this object during a flaring episode. The $Γ$ index of the intrinsic spectrum in the VHE $γ$-ray band is $2.04\pm0.12_{\rm stat}\pm0.15_{\rm sys}$. We find the source in an extreme state with a shift of the position of the synchrotron peak to frequencies above or equal to $10^{18}$\,Hz △ Less

Submitted 1 August, 2020; v1 submitted 11 June, 2020; originally announced June 2020.

Comments: Accepted in MNRAS. Corresponding authors: M. Manganaro, A. Arbet Engels, D. Dorner

Journal ref: MNRAS 496, 3912 (2020)

Testing two-component models on very-high-energy gamma-ray emitting BL Lac objects

Authors: MAGIC Collaboration, V. A. Acciari, S. Ansoldi, L. A. Antonelli, A. Arbet Engels, D. Baack, A. Babić, B. Banerjee, U. Barres de Almeida, J. A. Barrio, J. Becerra González, W. Bednarek, L. Bellizzi, E. Bernardini, A. Berti, J. Besenrieder, W. Bhattacharyya, C. Bigongiari, A. Biland, O. Blanch, G. Bonnoli, Ž. Bošnjak, G. Busetto, R. Carosi, G. Ceribella , et al. (179 additional authors not shown)

Abstract: Context. It has become evident that one-zone synchrotron self-Compton models are not always adequate for very-high-energy (VHE) gamma-ray emitting blazars. While two-component models are performing better, they are difficult to constrain due to the large number of free parameters. Aims. In this work, we make a first attempt to take into account the observational constraints from Very Long Baseline… ▽ More Context. It has become evident that one-zone synchrotron self-Compton models are not always adequate for very-high-energy (VHE) gamma-ray emitting blazars. While two-component models are performing better, they are difficult to constrain due to the large number of free parameters. Aims. In this work, we make a first attempt to take into account the observational constraints from Very Long Baseline Interferometry (VLBI) data, long-term light curves (radio, optical, and X-rays) and optical polarisation to limit the parameter space for a two-component model and test if it can still reproduce the observed spectral energy distribution (SED) of the blazars. Methods. We selected five TeV BL Lac objects based on the availability of VHE gamma-ray and optical polarisation data. We collected constraints for the jet parameters from VLBI observations. We evaluated the contributions of the two components to the optical flux by means of decomposition of long-term radio and optical light curves as well as modeling of the optical polarisation variability of the objects. We selected eight epochs for these five objects, based on the variability observed at VHE gamma rays, for which we constructed the SEDs that we then modeled with a two-component model. Results. We found parameter sets which can reproduce the broadband SED of the sources in the framework of two-component models considering all available observational constraints from VLBI observations. Moreover, the constraints obtained from the long-term behavior of the sources in the lower energy bands could be used to determine the region where the emission in each band originates. Finally, we attempted to use optical polarisation data to shed new light on the behavior of the two components in the optical band. Our observationally constrained two zone model allows explanation of the entire SED from radio to VHE with two co-located emission regions. △ Less

Submitted 8 June, 2020; originally announced June 2020.

Journal ref: A&A 640, A132 (2020)

NIR Counterparts to ULXs (III): Completing the photometric survey and selected spectroscopic results

Authors: K. M. López, M. Heida, P. G. Jonker, M. A. P. Torres, T. P. Roberts, D. J. Walton, D. -S. Moon, F. A. Harrison

Abstract: We present results from the remaining sources in our search for near-infrared (NIR) candidate counterparts to ultraluminous X-ray sources (ULXs) within $\simeq$ 10 Mpc. We observed 23 ULXs in 15 galaxies and detected NIR candidate counterparts to six of them. Two of these have an absolute magnitude consistent with a single red supergiant (RSG). Three counterparts are too bright for a RSG and spati… ▽ More We present results from the remaining sources in our search for near-infrared (NIR) candidate counterparts to ultraluminous X-ray sources (ULXs) within $\simeq$ 10 Mpc. We observed 23 ULXs in 15 galaxies and detected NIR candidate counterparts to six of them. Two of these have an absolute magnitude consistent with a single red supergiant (RSG). Three counterparts are too bright for a RSG and spatially extended, and thus we classify them as stellar clusters. The other candidate is too faint for a RSG. Additionally, we present the results of our NIR spectroscopic follow-up of five sources: four originally classified as RSG and one as a stellar cluster on the basis of previous photometry. The stellar cluster candidate is actually a nebula. Of the four RSGs candidates, one source has a broad H$α$ emission line redshifted by $\sim z = 1$, making it a background AGN. Two other sources show stellar spectra consistent with them being RSGs. The final RSG candidate is too faint to classify, but does not show strong (nebular) emission lines in its spectrum. After our search for NIR counterparts to 113 ULXs, where we detected a candidate counterpart for 38 ULXs, we have spectroscopically confirmed the nature of 12: five sources are nebulae, one source is not classified, one source is an AGN and five are RSGs. These possible five ULX-RSG binary systems would constitute $\simeq$ $(4 \pm 2)\%$ of the observed ULXs, a fraction almost four times larger than what was predicted by binary evolution simulations. △ Less

Submitted 4 June, 2020; originally announced June 2020.

Comments: 20 pages, 10 figures, 8 tables, accepted for publication in MNRAS

MAGIC observations of the diffuse $γ$-ray emission in the vicinity of the Galactic Centre

Authors: MAGIC Collaboration, V. A. Acciari, S. Ansoldi, L. A. Antonelli, A. Arbet Engels, D. Baack, A. Babić, B. Banerjee, U. Barres de Almeida, J. A. Barrio, J. Becerra González, W. Bednarek, L. Bellizzi, E. Bernardini, A. Berti, J. Besenrieder, W. Bhattacharyya, C. Bigongiari, A. Biland, O. Blanch, G. Bonnoli, Ž. Bošnjak, G. Busetto, R. Carosi, G. Ceribella , et al. (151 additional authors not shown)

Abstract: Aims: $γ$ rays can be used as a tracer in the search of sources of Galactic cosmic rays (CRs). We present deep observations of the Galactic Centre (GC) region with the MAGIC telescopes, which we use for inferring the underlying CR distribution. Methods: We observed the GC region for ${\approx}100$ hours with the MAGIC telescopes from 2012 to 2017, at high zenith angles (58-70~deg). This implies a… ▽ More Aims: $γ$ rays can be used as a tracer in the search of sources of Galactic cosmic rays (CRs). We present deep observations of the Galactic Centre (GC) region with the MAGIC telescopes, which we use for inferring the underlying CR distribution. Methods: We observed the GC region for ${\approx}100$ hours with the MAGIC telescopes from 2012 to 2017, at high zenith angles (58-70~deg). This implies a larger energy threshold, but also an increased effective collection area compared to low zenith observations. Using new software, we derive instrument response and background models, enabling us to study the diffuse emission in the region. We use pre-existing data of the gas distribution in the GC region to derive the underlying distribution of CRs. Results: We obtain a significant detection for all four model components used to fit our data (Sgr~A*, ``Arc'', G0.9+0.1, and an extended component for the Galactic Ridge). We find that the diffuse component is best described as a power-law with index 2 and an exponential cut-off at around 20~TeV with the significance of the cut-off being only 2~$σ$. The derived cosmic-ray profile hints to a peak at the GC position, with a measured profile index of $1.2 \pm 0.3$, supporting the hypothesis of a CR accelerator at the GC. We argue that the measurements of this profile are presently limited by our knowledge of the gas distribution in the GC vicinity. △ Less

Submitted 31 May, 2020; originally announced June 2020.

Comments: 10 pages, 10 figures. Accepted for publication by A&A. Corresponding authors: Christian Fruck (fruck@mpp.mpg.de), Ievgen Vovk (vovk@icrr.u-tokyo.ac.jp), Yuki Iwamura (iwamura@icrr.u-tokyo.ac.jp) and Marcel Strzys (strzys@icrr.u-tokyo.ac.jp)

Journal ref: A&A 642, A190 (2020)

Morphological and spectral study of 4FGL J1115.1-6118 in the region of the young massive stellar cluster NGC 3603

Authors: L. Saha, A. Domínguez, L. Tibaldo, S. Marchesi, M. Ajello, M. Lemoine-Goumard, M. López

Abstract: We report a detailed study of an unidentified gamma-ray source located in the region of the compact stellar cluster NGC 3603. This is a star-forming region (SFR) powered by a massive cluster of OB stars. A dedicated analysis of about 10 years of data from 10 GeV - 1 TeV, provided by the Large Area Telescope (LAT) onboard the Fermi Gamma-ray Space Telescope, yields the detection of a pointlike sour… ▽ More We report a detailed study of an unidentified gamma-ray source located in the region of the compact stellar cluster NGC 3603. This is a star-forming region (SFR) powered by a massive cluster of OB stars. A dedicated analysis of about 10 years of data from 10 GeV - 1 TeV, provided by the Large Area Telescope (LAT) onboard the Fermi Gamma-ray Space Telescope, yields the detection of a pointlike source at a significance of 9$σ$. The source photon spectrum can be described by a power-law model with best fit spectral index of $2.35 \pm 0.03$. In addition, the analysis of a deep Chandra image in the 0.5 - 7 keV band reliably rules out an extragalactic origin for the gamma rays. We also conclude that the broadband spectral energy distribution of the point source can be explained well with both leptonic and hadronic models. No firm evidence of association with any other classes of known gamma-ray emitters is found, therefore we speculate that 4FGL J1115.1-6118 is a gamma-ray emitting SFR. △ Less

Submitted 30 May, 2020; originally announced June 2020.

Comments: 16 pages, 9 figures, 7 tables. Submitted to ApJ

Search for the imprint of axion-like particles in the highest-energy photons of hard $γ$-ray blazars

Authors: Rolf Buehler, Galo Gallardo, Gernot Maier, Alberto Domínguez, Marcos López, Manuel Meyer

Abstract: Axion-like particles (ALPs), predicted in theories beyond the Standard Model, can have observational effects on the transparency of the Universe to $γ$ rays in the presence of magnetic fields. In this work, we search for effects compatible with the existence of ALPs with 80 months of data from the ${\it Fermi}$ Large Area Telescope, by comparing the distributions of observed highest energy photons… ▽ More Axion-like particles (ALPs), predicted in theories beyond the Standard Model, can have observational effects on the transparency of the Universe to $γ$ rays in the presence of magnetic fields. In this work, we search for effects compatible with the existence of ALPs with 80 months of data from the ${\it Fermi}$ Large Area Telescope, by comparing the distributions of observed highest energy photons from sources beyond redshifts of z $\geq$ 0.1 with theoretical predictions in the presence of ALPs. We find no evidence for an increased $γ$-ray transparency due to ALPs and therefore we set limits on the ALPs parameters assuming a value of the intergalactic magnetic field strength of 1 nG. Photon-ALP couplings above $10^{-11}$ GeV$^{-1}$ are excluded for ALP masses $m_{a}$ $\lesssim3.0$ neV. These limits exclude a region of the parameter space not covered by other $γ$-ray telescopes and are compatible with constraints imposed by other experiments. △ Less

Submitted 28 October, 2020; v1 submitted 20 April, 2020; originally announced April 2020.

Comments: 23 pages, 8 figures, 2 tables

HATS: A Ground-Based Telescope to Explore the THz Domain

Authors: C. G. Giménez de Castro, J. -P. Raulin, A. Valio, G. Alaia, V. Alvarenga, E. C. Bortolucci, S. H. Fernandes, C. Francile, T. Giorgetti, A. S. Kudaka, F. M. López, R. Marcon, A. Marún, M. Zaquela

Abstract: The almost unexplored frequency window from submillimeter to mid-infrared (mid-IR) may bring new clues about the particle acceleration and transport processes and the atmospheric thermal response during solar flares. Because of its technical complexity and the special atmospheric environment needed, observations at these frequencies are very sparse. The High Altitude THz Solar Photometer (HATS) is… ▽ More The almost unexplored frequency window from submillimeter to mid-infrared (mid-IR) may bring new clues about the particle acceleration and transport processes and the atmospheric thermal response during solar flares. Because of its technical complexity and the special atmospheric environment needed, observations at these frequencies are very sparse. The High Altitude THz Solar Photometer (HATS) is a full-Sun ground-based telescope designed to observe the continuum from the submillimeter to the mid-IR. It has a 457-mm spherical mirror with the sensor in its primary focus. The sensor is a Golay cell with high sensitivity in a very wide frequency range. The telescope has a polar mount, and a custom-built data acquisition system based on a 32 ksamples per second, 24 bits (72 dB dynamic range), 8 channels analog-to-digital board. Changing only the composition of the low- and band-pass filters in front of the Golay cell, the telescope can be setup to detect very different frequency bands; making the instrument very versatile. In this article we describe the telescope characteristics and its development status. Moreover, we give estimates of the expected fluxes during flares. △ Less

Submitted 28 March, 2020; originally announced March 2020.

Comments: Accepted for publication in Solar Physics. 13 pages, 3 Figures and 2 Tables

A search for dark matter in Triangulum II with the MAGIC telescopes

Authors: MAGIC Collaboration, V. A. Acciari, S. Ansoldi, L. A. Antonelli, A. Arbet Engels, D. Baack, A. Babić, B. Banerjee, U. Barres de Almeida, J. A. Barrio, J. Becerra González, W. Bednarek, L. Bellizzi, E. Bernardini, A. Berti, J. Besenrieder, W. Bhattacharyya, C. Bigongiari, A. Biland, O. Blanch, G. Bonnoli, Ž. Bošnjak, G. Busetto, R. Carosi, G. Ceribella , et al. (150 additional authors not shown)

Abstract: We present the first results from very-high-energy observations of the dwarf spheroidal satellite candidate Triangulum II with the MAGIC telescopes from 62.4 hours of good-quality data taken between August 2016 and August 2017. We find no gamma-ray excess in the direction of Triangulum II, and upper limits on both the differential and integral gamma-ray flux are presented. Currently, the kinematic… ▽ More We present the first results from very-high-energy observations of the dwarf spheroidal satellite candidate Triangulum II with the MAGIC telescopes from 62.4 hours of good-quality data taken between August 2016 and August 2017. We find no gamma-ray excess in the direction of Triangulum II, and upper limits on both the differential and integral gamma-ray flux are presented. Currently, the kinematics of Triangulum II are affected by large uncertainties leading to a bias in the determination of the properties of its dark matter halo. Using a scaling relation between the annihilation J-factor and heliocentric distance of well-known dwarf spheroidal galaxies, we estimate an annihilation J-factor for Triangulum II for WIMP dark matter of $\log[J_{\text{ann}}({0.5^{\circ}})/$ GeV$^{2}$ cm$^{-5}] = 19.35 \pm 0.37$. We also derive a dark matter density profile for the object relying on results from resolved simulations of Milky Way sized dark matter halos. We obtain 95% confidence-level limits on the thermally averaged annihilation cross section for WIMP annihilation into various Standard Model channels. The most stringent limits are obtained in the $τ^{+}τ^{-}$ final state, where a cross section for annihilation down to $\langle σ_{\text{ann}} v \rangle = 3.05 \times 10^{-24}$ cm$^{3}$ s$^{-1}$ is excluded. △ Less

Submitted 11 March, 2020; originally announced March 2020.

Comments: 14 pages, 8 figures, accepted for publication in Physics of the Dark Universe

Journal ref: Physics of the Dark Universe Volume 28, May 2020, 100529

A ringed pole-on outflow from DO Tauri revealed by ALMA

Authors: M. Fernández López, L. A. Zapata, L. F. Rodríguez, M. M. Vazzano, A. E. Guzmán, R. López

Abstract: We present new ALMA Band 6 observations including the CO(2-1) line and 1.3 mm continuum emission from the surroundings of the young stellar object DO Tauri. The ALMA CO molecular data show three different series of rings at different radial velocities. These rings have radii around 220 au and 800 au. We make individual fits to the rings and note that their centers are aligned with DO Tauri and its… ▽ More We present new ALMA Band 6 observations including the CO(2-1) line and 1.3 mm continuum emission from the surroundings of the young stellar object DO Tauri. The ALMA CO molecular data show three different series of rings at different radial velocities. These rings have radii around 220 au and 800 au. We make individual fits to the rings and note that their centers are aligned with DO Tauri and its optical high-velocity jet. In addition, we notice that the velocity of these structures increases with the separation from the young star. We discuss the data under the hypothesis that the rings represent velocity cuts through three outflowing shells that are possibly driven by a wide-angle wind, dragging the environment material along a direction close to the line of sight (i=19°). We estimate the dynamical ages, the mass, the momentum and the energy of each individual outflow shell and those of the whole outflow. The results are in agreement with those found in outflows from Class II sources. We make a rough estimate for the size of the jet/wind launching region, which needs to be of <15 au. We report the physical characteristics of DO Tauri's disk continuum emission (almost face-on and with a projected major axis in the north-south direction) and its velocity gradient orientation (north-south), indicative of disk rotation for a 1-2 Msun central star. Finally we show an HST [SII] image of the optical jet and report a measurement of its orientation in the plane of the sky. △ Less

Submitted 2 March, 2020; originally announced March 2020.

Comments: 26 pages, 11 figures, accepted by ApJ

The Great Markarian 421 Flare of February 2010: Multiwavelength variability and correlation studies

Authors: A. U. Abeysekara, W. Benbow, R. Bird, A. Brill, R. Brose, M. Buchovecky, J. H. Buckley, J. L. Christiansen, A. J. Chromey, M. K. Daniel, J. Dumm, A. Falcone, Q. Feng, J. P. Finley, L. Fortson, A. Furniss, N. Galante, A. Gent, G. H. Gillanders, C. Giuri, O. Gueta, T. Hassan, O. Hervet, J. Holder, G. Hughes , et al. (234 additional authors not shown)

Abstract: We report on variability and correlation studies using multiwavelength observations of the blazar Mrk 421 during the month of February, 2010 when an extraordinary flare reaching a level of $\sim$27~Crab Units above 1~TeV was measured in very-high-energy (VHE) $γ$-rays with the VERITAS observatory. This is the highest flux state for Mrk 421 ever observed in VHE $γ$-rays. Data are analyzed from a co… ▽ More We report on variability and correlation studies using multiwavelength observations of the blazar Mrk 421 during the month of February, 2010 when an extraordinary flare reaching a level of $\sim$27~Crab Units above 1~TeV was measured in very-high-energy (VHE) $γ$-rays with the VERITAS observatory. This is the highest flux state for Mrk 421 ever observed in VHE $γ$-rays. Data are analyzed from a coordinated campaign across multiple instruments including VHE $γ$-ray (VERITAS, MAGIC), high-energy (HE) $γ$-ray (Fermi-LAT), X-ray (Swift}, RXTE, MAXI), optical (including the GASP-WEBT collaboration and polarization data) and radio (Metsähovi, OVRO, UMRAO). Light curves are produced spanning multiple days before and after the peak of the VHE flare, including over several flare `decline' epochs. The main flare statistics allow 2-minute time bins to be constructed in both the VHE and optical bands enabling a cross-correlation analysis that shows evidence for an optical lag of $\sim$25-55 minutes, the first time-lagged correlation between these bands reported on such short timescales. Limits on the Doppler factor ($δ\gtrsim 33$) and the size of the emission region ($ δ^{-1}R_B \lesssim 3.8\times 10^{13}\,\,\mbox{cm}$) are obtained from the fast variability observed by VERITAS during the main flare. Analysis of 10-minute-binned VHE and X-ray data over the decline epochs shows an extraordinary range of behavior in the flux-flux relationship: from linear to quadratic to lack of correlation to anti-correlation. Taken together, these detailed observations of an unprecedented flare seen in Mrk 421 are difficult to explain by the classic single-zone synchrotron self-Compton model. △ Less

Submitted 10 February, 2020; originally announced February 2020.

Comments: 41 pages including 3 appendices, 13 figures; version accepted to Astrophysical Journal

Journal ref: The Astrophysical Journal, Volume 890, Number 2 (2020)

Broadband characterisation of the very intense TeV flares of the blazar 1ES 1959+650 in 2016

Authors: MAGIC Collaboration, V. A. Acciari, S. Ansoldi, L. A. Antonelli, A. Arbet Engels, D. Baack, A. Babić, B. Banerjee, U. Barres de Almeida, J. A. Barrio, J. Becerra González, W. Bednarek, L. Bellizzi, E. Bernardini, A. Berti, J. Besenrieder, W. Bhattacharyya, C. Bigongiari, A. Biland, O. Blanch, G. Bonnoli, Z. Bosnjak, G. Busetto, R. Carosi, G. Ceribella , et al. (142 additional authors not shown)

Abstract: 1ES 1959+650 is a bright TeV high-frequency-peaked BL Lac object exhibiting interesting features like "orphan" TeV flares and a broad emission in the high-energy regime, that are difficult to interpret using conventional one-zone Synchrotron Self-Compton (SSC) scenarios. We report the results from the Major Atmospheric Gamma Imaging Cherenkov (MAGIC) observations in 2016 along with the multi-wavel… ▽ More 1ES 1959+650 is a bright TeV high-frequency-peaked BL Lac object exhibiting interesting features like "orphan" TeV flares and a broad emission in the high-energy regime, that are difficult to interpret using conventional one-zone Synchrotron Self-Compton (SSC) scenarios. We report the results from the Major Atmospheric Gamma Imaging Cherenkov (MAGIC) observations in 2016 along with the multi-wavelength data from the Fermi Large Area Telescope (LAT) and Swift instruments. MAGIC observed 1ES 1959+650 with different emission levels in the very-high-energy (VHE, E >100 GeV) gamma-ray band during 2016. In the long-term data, the X-ray spectrum becomes harder with increasing flux and a hint of a similar trend is also visible in the VHE band. An exceptionally high VHE flux reaching ~ 3 times the Crab Nebula flux was measured by MAGIC on the 13th, 14th of June and 1st July 2016 (the highest flux observed since 2002). During these flares, the high-energy peak of the spectral energy distribution (SED) lies in the VHE domain and extends up to several TeV. The spectrum in the gamma-ray (both Fermi-LAT and VHE bands) and the X-ray bands are quite hard. On 13th June and 1st July 2016, the source showed rapid variations of the VHE flux within timescales of less than an hour. A simple one-zone SSC model can describe the data during the flares requiring moderate to high values of the Doppler factors (>=30-60). Alternatively, the high-energy peak of the SED can be explained by a purely hadronic model attributed to proton-synchrotron radiation with jet power L_{jet}~10^{46} erg/s and under high values of the magnetic field strength (~100 G) and maximum proton energy (~few EeV). Mixed lepto-hadronic models require super-Eddington values of the jet power. We conclude that it is difficult to get detectable neutrino emission from the source during the extreme VHE flaring period of 2016. △ Less

Submitted 31 January, 2020; originally announced February 2020.

Comments: 18 pages, 9 figures, accepted to A&A

Journal ref: A&A 638, A14 (2020)

Bounds on Lorentz invariance violation from MAGIC observation of GRB 190114C

Authors: MAGIC Collaboration, V. A. Acciari, S. Ansoldi, L. A. Antonelli, A. Arbet Engels, D. Baack, A. Babić, B. Banerjee, U. Barres de Almeida, J. A. Barrio, J. Becerra González, W. Bednarek, L. Bellizzi, E. Bernardini, A. Berti, J. Besenrieder, W. Bhattacharyya, C. Bigongiari, A. Biland, O. Blanch, G. Bonnoli, Ž. Bošnjak, G. Busetto, R. Carosi, G. Ceribella , et al. (159 additional authors not shown)

Abstract: On January 14, 2019, the Major Atmospheric Gamma Imaging Cherenkov telescopes detected GRB 190114C above 0.2 TeV, recording the most energetic photons ever observed from a gamma-ray burst. We use this unique observation to probe an energy dependence of the speed of light in vacuo for photons as predicted by several quantum gravity models. Based on a set of assumptions on the possible intrinsic spe… ▽ More On January 14, 2019, the Major Atmospheric Gamma Imaging Cherenkov telescopes detected GRB 190114C above 0.2 TeV, recording the most energetic photons ever observed from a gamma-ray burst. We use this unique observation to probe an energy dependence of the speed of light in vacuo for photons as predicted by several quantum gravity models. Based on a set of assumptions on the possible intrinsic spectral and temporal evolution, we obtain competitive lower limits on the quadratic leading order of speed of light modification. △ Less

Submitted 10 July, 2020; v1 submitted 27 January, 2020; originally announced January 2020.

Comments: 8 pages, 5 figures, Published in Physical Review Letters, Corresponding authors: Tomislav Terzić (tterzic@phy.uniri.hr), Giacomo D'Amico (damico@mppmu.mpg.de), Daniel Kerszberg (dkerszberg@ifae.es), Cedric Perennes (perennes.cedric@gmail.com)

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

MAGIC very large zenith angle observations of the Crab Nebula up to 100 TeV

Authors: MAGIC Collaboration, V. A. Acciari, S. Ansoldi, L. A. Antonelli, A. Arbet Engels, D. Baack, A. Babić, B. Banerjee, U. Barres de Almeida, J. A. Barrio, J. Becerra González, W. Bednarek, L. Bellizzi, E. Bernardini, A. Berti, J. Besenrieder, W. Bhattacharyya, C. Bigongiari, A. Biland, O. Blanch, G. Bonnoli, Ž. Bošnjak, G. Busetto, R. Carosi, G. Ceribella , et al. (147 additional authors not shown)

Abstract: Aims: We aim to measure the Crab Nebula gamma-ray spectral energy distribution in the ~100 TeV energy domain and test the validity of existing leptonic emission models at these high energies. Methods: We use the novel very large zenith angle observations with the MAGIC telescope system to increase the collection area above 10 TeV. We also develop an auxiliary procedure of monitoring atmospheric… ▽ More Aims: We aim to measure the Crab Nebula gamma-ray spectral energy distribution in the ~100 TeV energy domain and test the validity of existing leptonic emission models at these high energies. Methods: We use the novel very large zenith angle observations with the MAGIC telescope system to increase the collection area above 10 TeV. We also develop an auxiliary procedure of monitoring atmospheric transmission in order to assure proper calibration of the accumulated data. This employs recording of optical images of the stellar field next to the source position, which provides a better than 10% accuracy for the transmission measurements. Results: We demonstrate that MAGIC very large zenith angle observations yield a collection area larger than a square kilometer. In only ~56 hr of observations, we detect the gamma-ray emission from the Crab Nebula up to 100 TeV, thus providing the highest energy measurement of this source to date with Imaging Atmospheric Cherenkov Telescopes. Comparing accumulated and archival MAGIC and Fermi/LAT data with some of the existing emission models, we find that none of them provides an accurate description of the 1 GeV to 100 TeV gamma-ray signal. △ Less

Submitted 26 January, 2020; originally announced January 2020.

Comments: 10 pages, 10 figures, accepted to A&A

Journal ref: A&A 635, A158 (2020)

Unravelling the complex behavior of Mrk 421 with simultaneous X-ray and VHE observations during an extreme flaring activity in April 2013

Authors: MAGIC collaboration, V. A. Acciari, S. Ansoldi, L. A. Antonelli, A. Arbet Engels, D. Baack, A. Babic, B. Banerjee, U. Barres de Almeida, J. A. Barrio, J. Becerra Gonzalez, W. Bednarek, L. Bellizzi, E. Bernardini, A. Berti, J. Besenrieder, W. Bhattacharyya, C. Bigongiari, A. Biland, O. Blanch, G. Bonnoli, Z. Bosnjak, G. Busetto, R. Carosi, G. Ceribella , et al. (215 additional authors not shown)

Abstract: We report on a multi-band variability and correlation study of the TeV blazar Mrk 421 during an exceptional flaring activity observed from 2013 April 11 to 2013 April 19. The study uses, among others, data from GASP-WEBT, Swift, NuSTAR, Fermi-LAT, VERITAS, and MAGIC. The large blazar activity, and the 43 hours of simultaneous NuSTAR and MAGIC/VERITAS observations, permitted variability studies on… ▽ More We report on a multi-band variability and correlation study of the TeV blazar Mrk 421 during an exceptional flaring activity observed from 2013 April 11 to 2013 April 19. The study uses, among others, data from GASP-WEBT, Swift, NuSTAR, Fermi-LAT, VERITAS, and MAGIC. The large blazar activity, and the 43 hours of simultaneous NuSTAR and MAGIC/VERITAS observations, permitted variability studies on 15 minute time bins, and over three X-ray bands (3-7 keV, 7-30 keV and 30-80 keV) and three very-high-energy (>0.1 TeV, hereafter VHE) gamma-ray bands (0.2-0.4 TeV, 0.4-0.8 TeV and >0.8 TeV). We detected substantial flux variations on multi-hour and sub-hour timescales in all the X-ray and VHE gamma-ray bands. The characteristics of the sub-hour flux variations are essentially energy-independent, while the multi-hour flux variations can have a strong dependence on the energy of the X-ray and the VHE gamma rays. The three VHE bands and the three X-ray bands are positively correlated with no time-lag, but the strength and the characteristics of the correlation changes substantially over time and across energy bands. Our findings favour multi-zone scenarios for explaining the achromatic/chromatic variability of the fast/slow components of the light curves, as well as the changes in the flux-flux correlation on day-long timescales. We interpret these results within a magnetic reconnection scenario, where the multi-hour flux variations are dominated by the combined emission from various plasmoids of different sizes and velocities, while the sub-hour flux variations are dominated by the emission from a single small plasmoid moving across the magnetic reconnection layer. △ Less

Submitted 23 January, 2020; originally announced January 2020.

Comments: 46 pages, 17 figures, 9 tables, accepted for publication in ApJ Supplements

Stellar Properties of the Host Galaxy of an Ultraluminous X-ray Source in NGC 5252

Authors: Minjin Kim, Kristhell M. López, Peter G. Jonker, Luis C. Ho, Myungshin Im

Abstract: An ultraluminous X-ray source (ULX) in NGC 5252 has been known as a strong candidate for an off-nuclear intermediate-mass black hole. We present near-infrared imaging data of the ULX obtained with the William Herschel Telescope. Using this data we estimate a stellar mass associated with the ULX of $\approx 10^{7.9\pm0.1}M_{\rm \odot}$ , suggesting that it could be (the remnant of) a dwarf galaxy t… ▽ More An ultraluminous X-ray source (ULX) in NGC 5252 has been known as a strong candidate for an off-nuclear intermediate-mass black hole. We present near-infrared imaging data of the ULX obtained with the William Herschel Telescope. Using this data we estimate a stellar mass associated with the ULX of $\approx 10^{7.9\pm0.1}M_{\rm \odot}$ , suggesting that it could be (the remnant of) a dwarf galaxy that is in the process of merging with NGC 5252. Based on a correlation between the mass of the central black hole (BH) and host galaxy, the ULX is powered by a $10^5M_{\rm \odot}$ black hole. Alternatively, if the BH mass is $\approx 10^6M_{\rm \odot}$ or larger, the host galaxy of the ULX must have been heavily stripped during the merger. The ULX $K_s$-band luminosity is two orders of magnitude smaller than that expected from an ordinary active galactic nucleus with the observed [O III] luminosity, which also suggests the ULX lacks a dusty torus. We discuss how these findings provide suggestive evidence that the ULX is hosting an intermediate-mass black hole. △ Less

Submitted 22 January, 2020; originally announced January 2020.

Comments: To appear in MNRAS

Study of the variable broadband emission of Markarian 501 during the most extreme Swift X-ray activity

Authors: MAGIC Collaboration, V. A. Acciari, S. Ansoldi, L. A. Antonelli, A. Babić, B. Banerjee, U. Barres de Almeida, J. A. Barrio, J. Becerra González, W. Bednarek, E. Bernardini, A. Berti, J. Besenrieder, W. Bhattacharyya, C. Bigongiari, O. Blanch, G. Bonnoli, G. Busetto, R. Carosi, G. Ceribella, S. Cikota, S. M. Colak, P. Colin, E. Colombo, J. L. Contreras , et al. (163 additional authors not shown)

Abstract: Mrk501 is a very high-energy (VHE) gamma-ray blazar located at z=0.034. During a period of two weeks in July 2014, the highest X-ray activity of Mrk501 was observed in ~14 years of operation of the Neil Gehrels Swift Gamma-ray Burst Observatory. We characterize the broadband variability of Mrk501 from radio to VHE gamma rays, and evaluate whether it can be interpreted within theoretical scenarios… ▽ More Mrk501 is a very high-energy (VHE) gamma-ray blazar located at z=0.034. During a period of two weeks in July 2014, the highest X-ray activity of Mrk501 was observed in ~14 years of operation of the Neil Gehrels Swift Gamma-ray Burst Observatory. We characterize the broadband variability of Mrk501 from radio to VHE gamma rays, and evaluate whether it can be interpreted within theoretical scenarios widely used to explain the broadband emission from blazars. The temporal evolution of the most prominent and variable segments of the SED is described with a one-zone synchrotron self-Compton model with variations in the break energy of the electron energy distribution (EED), and with some adjustments in the magnetic field strength and spectral shape of the EED. A narrow feature at ~3 TeV was observed in the VHE spectrum measured on 2014 July 19 (MJD 56857.98), which is the day with the highest X-ray flux ($>0.3$ keV) measured during the entire Swift mission. This feature is inconsistent with the classical analytic functions to describe the measured VHE spectra (power law, log-parabola, and log-parabola with exponential cutoff) at more than 3$σ$. A fit with a log-parabola plus a narrow component is preferred over the fit with a single log-parabola at more than 4$σ$, and a dedicated Monte Carlo simulation estimated the significance of this extra component to be larger than 3$σ$. Under the assumption that this VHE spectral feature is real, we show that it can be reproduced with three distinct theoretical scenarios: a) a pileup in the EED due to stochastic acceleration; b) a structured jet with two-SSC emitting regions, with one region dominated by an extremely narrow EED; and c) an emission from an IC pair cascade. △ Less

Submitted 16 March, 2021; v1 submitted 21 January, 2020; originally announced January 2020.

Comments: Accepted for publication in A&A. Corresponding authors: J. Becerra González (jbecerra@iac.es), D. Paneque (dpaneque@mppmu.mpg.de), C. Wendel (cwendel@astro.uni-wuerzburg.de)

Journal ref: A&A 637, A86 (2020)