Showing 1–50 of 176 results for author: Schmidt, K

Constraints on Lorentz invariance violation from the extraordinary Mrk 421 flare of 2014 using a novel analysis method

Authors: MAGIC Collaboration, S. Abe, J. Abhir, A. Abhishek, V. A. Acciari, A. Aguasca-Cabot, I. Agudo, T. Aniello, S. Ansoldi, L. A. Antonelli, A. Arbet Engels, C. Arcaro, M. Artero, K. Asano, A. Babić, A. Baquero, U. Barres de Almeida, J. A. Barrio, I. Batković, A. Bautista, J. Baxter, J. Becerra González, W. Bednarek, E. Bernardini, J. Bernete , et al. (192 additional authors not shown)

Abstract: The Lorentz Invariance Violation (LIV), a proposed consequence of certain quantum gravity (QG) scenarios, could instigate an energy-dependent group velocity for ultra-relativistic particles. This energy dependence, although suppressed by the massive QG energy scale $E_\mathrm{QG}$, expected to be on the level of the Planck energy $1.22 \times 10^{19}$ GeV, is potentially detectable in astrophysica… ▽ More The Lorentz Invariance Violation (LIV), a proposed consequence of certain quantum gravity (QG) scenarios, could instigate an energy-dependent group velocity for ultra-relativistic particles. This energy dependence, although suppressed by the massive QG energy scale $E_\mathrm{QG}$, expected to be on the level of the Planck energy $1.22 \times 10^{19}$ GeV, is potentially detectable in astrophysical observations. In this scenario, the cosmological distances traversed by photons act as an amplifier for this effect. By leveraging the observation of a remarkable flare from the blazar Mrk\,421, recorded at energies above 100 GeV by the MAGIC telescopes on the night of April 25 to 26, 2014, we look for time delays scaling linearly and quadratically with the photon energies. Using for the first time in LIV studies a binned-likelihood approach we set constraints on the QG energy scale. For the linear scenario, we set $95\%$ lower limits $E_\mathrm{QG}>2.7\times10^{17}$ GeV for the subluminal case and $E_\mathrm{QG}> 3.6 \times10^{17}$ GeV for the superluminal case. For the quadratic scenario, the $95\%$ lower limits for the subluminal and superluminal cases are $E_\mathrm{QG}>2.6 \times10^{10}$ GeV and $E_\mathrm{QG}>2.5\times10^{10}$ GeV, respectively. △ Less

Submitted 11 June, 2024; originally announced June 2024.

Apparent non-variable stars from the Kepler mission

Authors: E. Paunzen, F. Binder, A. Cyniburk, M. N. Duffek, F. Haberhauer, C. Heinreichsberger, H. Kohlhofer, L. Kueß, H. M. Maitzen, T. Saalmann, A. M. Schanz, S. Schauer, K. Schmidt, A. Tokareva, I. Wizani

Abstract: The analysis of non-variable stars is generally neglected in the literature. However, such objects are needed for many calibration processes and for testing pulsational models. The photometric time series of the Kepler satellite mission still stand as the most accurate data available today and are excellently suited to the search for non-variable stars. We analysed all long-cadence light curves fo… ▽ More The analysis of non-variable stars is generally neglected in the literature. However, such objects are needed for many calibration processes and for testing pulsational models. The photometric time series of the Kepler satellite mission still stand as the most accurate data available today and are excellently suited to the search for non-variable stars. We analysed all long-cadence light curves for stars not reported as a variable so far from the Kepler satellite mission. Using the known characteristics and flaws of these data sets, we defined three different frequency ranges where we searched for non-variability. We used the Lomb-Scargle periodogram and the false-alarm probability (FAP) to analyse the cleaned data sets of 138 451 light curves. We then used log FAP > -2 to define a star as "non-variable" in the ranges below 0.1 c/d, 0.1 to 2.0 c/d, and 2.0 to 25.0 c/d, respectively. Furthermore, we also calculated the standard deviation of the mean light curve to obtain another parameter. In total, we found 14 154 stars that fulfil the set △ Less

Submitted 10 June, 2024; originally announced June 2024.

Comments: 7 pages, 5 figures, accepted for Astronomy & Astrophysics

The variability patterns of the TeV blazar PG 1553+113 from a decade of MAGIC and multi-band observations

Authors: MAGIC Collaboration, H. Abe, S. Abe, J. Abhir, 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, I. Batković, J. Baxter, J. Becerra González, E. Bernardini, J. Bernete, A. Berti, J. Besenrieder, C. Bigongiari , et al. (242 additional authors not shown)

Abstract: PG 1553+113 is one of the few blazars with a convincing quasi-periodic emission in the gamma-ray band. The source is also a very high-energy (VHE; >100 GeV) gamma-ray emitter. To better understand its properties and identify the underlying physical processes driving its variability, the MAGIC Collaboration initiated a multiyear, multiwavelength monitoring campaign in 2015 involving the OVRO 40-m a… ▽ More PG 1553+113 is one of the few blazars with a convincing quasi-periodic emission in the gamma-ray band. The source is also a very high-energy (VHE; >100 GeV) gamma-ray emitter. To better understand its properties and identify the underlying physical processes driving its variability, the MAGIC Collaboration initiated a multiyear, multiwavelength monitoring campaign in 2015 involving the OVRO 40-m and Medicina radio telescopes, REM, KVA, and the MAGIC telescopes, Swift and Fermi satellites, and the WEBT network. The analysis presented in this paper uses data until 2017 and focuses on the characterization of the variability. The gamma-ray data show a (hint of a) periodic signal compatible with literature, but the X-ray and VHE gamma-ray data do not show statistical evidence for a periodic signal. In other bands, the data are compatible with the gamma-ray period, but with a relatively high p-value. The complex connection between the low and high-energy emission and the non-monochromatic modulation and changes in flux suggests that a simple one-zone model is unable to explain all the variability. Instead, a model including a periodic component along with multiple emission zones is required. △ Less

Submitted 4 March, 2024; originally announced March 2024.

Comments: Accepted for publication in Monthly Notices of the Royal Astronomical Society. 19 pages, 9 figures. Corresponding authors: Elisa Prandini, Antonio Stamerra, Talvikki Hovatta

Performance and first measurements of the MAGIC Stellar Intensity Interferometer

Authors: MAGIC Collaboration, S. Abe, J. Abhir, V. A. Acciari, A. Aguasca-Cabot, I. Agudo, T. Aniello, S. Ansoldi, L. A. Antonelli, A. Arbet Engels, C. Arcaro, M. Artero, K. Asano, A. Babić, A. Baquero, U. Barres de Almeida, J. A. Barrio, I. Batković, A. Bautista, J. Baxter, J. Becerra González, E. Bernardini, M. Bernardos, J. Bernete, A. Berti , et al. (195 additional authors not shown)

Abstract: In recent years, a new generation of optical intensity interferometers has emerged, leveraging the existing infrastructure of Imaging Atmospheric Cherenkov Telescopes (IACTs). The MAGIC telescopes host the MAGIC-SII system (Stellar Intensity Interferometer), implemented to investigate the feasibility and potential of this technique on IACTs. After the first successful measurements in 2019, the sys… ▽ More In recent years, a new generation of optical intensity interferometers has emerged, leveraging the existing infrastructure of Imaging Atmospheric Cherenkov Telescopes (IACTs). The MAGIC telescopes host the MAGIC-SII system (Stellar Intensity Interferometer), implemented to investigate the feasibility and potential of this technique on IACTs. After the first successful measurements in 2019, the system was upgraded and now features a real-time, dead-time-free, 4-channel, GPU-based correlator. These hardware modifications allow seamless transitions between MAGIC's standard very-high-energy gamma-ray observations and optical interferometry measurements within seconds. We establish the feasibility and potential of employing IACTs as competitive optical Intensity Interferometers with minimal hardware adjustments. The measurement of a total of 22 stellar diameters are reported, 9 corresponding to reference stars with previous comparable measurements, and 13 with no prior measurements. A prospective implementation involving telescopes from the forthcoming Cherenkov Telescope Array Observatory's northern hemisphere array, such as the first prototype of its Large-Sized Telescopes, LST-1, is technically viable. This integration would significantly enhance the sensitivity of the current system and broaden the UV-plane coverage. This advancement would enable the system to achieve competitive sensitivity with the current generation of long-baseline optical interferometers over blue wavelengths. △ Less

Submitted 7 February, 2024; originally announced February 2024.

Comments: 18 pages, 13 figures, submitted to MNRAS

Insights into the broad-band emission of the TeV blazar Mrk 501 during the first X-ray polarization measurements

Authors: S. Abe, J. Abhir, V. A. Acciari, A. Aguasca-Cabot, I. Agudo, T. Aniello, S. Ansoldi, L. A. Antonelli, A. Arbet Engels, C. Arcaro, K. Asano, A. Babić, A. Baquero, U. Barres de Almeida, J. A. Barrio, I. Batković, A. Bautista, J. Baxter, J. Becerra González, W. Bednarek, E. Bernardini, M. Bernardos, J. Bernete, A. Berti, J. Besenrieder , et al. (239 additional authors not shown)

Abstract: We present the first multi-wavelength study of Mrk 501 including very-high-energy (VHE) gamma-ray observations simultaneous to X-ray polarization measurements from the Imaging X-ray Polarimetry Explorer (IXPE). We use radio-to-VHE data from a multi-wavelength campaign organized between 2022-03-01 and 2022-07-19. The observations were performed by MAGIC, Fermi-LAT, NuSTAR, Swift (XRT and UVOT), and… ▽ More We present the first multi-wavelength study of Mrk 501 including very-high-energy (VHE) gamma-ray observations simultaneous to X-ray polarization measurements from the Imaging X-ray Polarimetry Explorer (IXPE). We use radio-to-VHE data from a multi-wavelength campaign organized between 2022-03-01 and 2022-07-19. The observations were performed by MAGIC, Fermi-LAT, NuSTAR, Swift (XRT and UVOT), and several instruments covering the optical and radio bands. During the IXPE pointings, the VHE state is close to the average behavior with a 0.2-1 TeV flux of 20%-50% the emission of the Crab Nebula. Despite the average VHE activity, an extreme X-ray behavior is measured for the first two IXPE pointings in March 2022 with a synchrotron peak frequency >1 keV. For the third IXPE pointing in July 2022, the synchrotron peak shifts towards lower energies and the optical/X-ray polarization degrees drop. The X-ray polarization is systematically higher than at lower energies, suggesting an energy-stratification of the jet. While during the IXPE epochs the polarization angle in the X-ray, optical and radio bands align well, we find a clear discrepancy in the optical and radio polarization angles in the middle of the campaign. We model the broad-band spectra simultaneous to the IXPE pointings assuming a compact zone dominating in the X-rays and VHE, and an extended zone stretching further downstream the jet dominating the emission at lower energies. NuSTAR data allow us to precisely constrain the synchrotron peak and therefore the underlying electron distribution. The change between the different states observed in the three IXPE pointings can be explained by a change of magnetization and/or emission region size, which directly connects the shift of the synchrotron peak to lower energies with the drop in polarization degree. △ Less

Submitted 16 January, 2024; originally announced January 2024.

Comments: Accepted for publication in Astronomy & Astrophysics. 19 pages, 9 figures. Corresponding authors: Lea Heckmann, Axel Arbet Engels, David Paneque

Constraints on axion-like particles with the Perseus Galaxy Cluster with MAGIC

Authors: MAGIC Collaboration, H. Abe, S. Abe, J. Abhir, 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, J. Bernete, A. Berti , et al. (189 additional authors not shown)

Abstract: Axion-like particles (ALPs) are pseudo-Nambu-Goldstone bosons that emerge in various theories beyond the standard model. These particles can interact with high-energy photons in external magnetic fields, influencing the observed gamma-ray spectrum. This study analyzes 41.3 hrs of observational data from the Perseus Galaxy Cluster collected with the MAGIC telescopes. We focused on the spectra the r… ▽ More Axion-like particles (ALPs) are pseudo-Nambu-Goldstone bosons that emerge in various theories beyond the standard model. These particles can interact with high-energy photons in external magnetic fields, influencing the observed gamma-ray spectrum. This study analyzes 41.3 hrs of observational data from the Perseus Galaxy Cluster collected with the MAGIC telescopes. We focused on the spectra the radio galaxy in the center of the cluster: NGC 1275. By modeling the magnetic field surrounding this target, we searched for spectral indications of ALP presence. Despite finding no statistical evidence of ALP signatures, we were able to exclude ALP models in the sub-micro electronvolt range. Our analysis improved upon previous work by calculating the full likelihood and statistical coverage for all considered models across the parameter space. Consequently, we achieved the most stringent limits to date for ALP masses around 50 neV, with cross sections down to $g_{aγ} = 3 \times 10^{-12}$ GeV$^{-1}$. △ Less

Submitted 15 January, 2024; originally announced January 2024.

Comments: 25 pages, 10 figures, accepted for publication in Physics of the Dark Universe

First characterization of the emission behavior of Mrk421 from radio to VHE gamma rays with simultaneous X-ray polarization measurements

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

Abstract: We perform the first broadband study of Mrk421 from radio to TeV gamma rays with simultaneous measurements of the X-ray polarization from IXPE. The data were collected within an extensive multiwavelength campaign organized between May and June 2022 using MAGIC, Fermi-LAT, NuSTAR, XMM-Newton, Swift, and several optical and radio telescopes to complement IXPE. During the IXPE exposures, the measured… ▽ More We perform the first broadband study of Mrk421 from radio to TeV gamma rays with simultaneous measurements of the X-ray polarization from IXPE. The data were collected within an extensive multiwavelength campaign organized between May and June 2022 using MAGIC, Fermi-LAT, NuSTAR, XMM-Newton, Swift, and several optical and radio telescopes to complement IXPE. During the IXPE exposures, the measured 0.2-1 TeV flux is close to the quiescent state and ranges from 25% to 50% of the Crab Nebula without intra-night variability. Throughout the campaign, the VHE and X-ray emission are positively correlated at a $4σ$ significance level. The IXPE measurements unveil a X-ray polarization degree that is a factor of 2-5 higher than in the optical/radio bands; that implies an energy-stratified jet in which the VHE photons are emitted co-spatially with the X-rays, in the vicinity of a shock front. The June 2022 observations exhibit a rotation of the X-ray polarization angle. Despite no simultaneous VHE coverage being available during a large fraction of the swing, the Swift-XRT monitoring unveils an X-ray flux increase with a clear spectral hardening. It suggests that flares in high synchrotron peaked blazars can be accompanied by a polarization angle rotation, as observed in some flat spectrum radio quasars. Finally, during the polarization angle rotation, NuSTAR data reveal two contiguous spectral hysteresis loops in opposite directions (clockwise and counter-clockwise), implying important changes in the particle acceleration efficiency on $\sim$hour timescales. △ Less

Submitted 17 December, 2023; originally announced December 2023.

Comments: Accepted for publication in Astronomy & Astrophysics. 29 pages, 22 figures. Corresponding authors: Axel Arbet Engels, Felix Schmuckermaier, David Paneque

A Collection of German Science Interests in the Next Generation Very Large Array

Authors: M. Kadler, D. A. Riechers, J. Agarwal, A. -K. Baczko, H. Beuther, F. Bigiel, T. Birnstiel, B. Boccardi, D. J. Bomans, L. Boogaard, T. T. Braun, S. Britzen, M. Brüggen, A. Brunthaler, P. Caselli, D. Elsässer, S. von Fellenberg, M. Flock, C. M. Fromm, L. Fuhrmann, P. Hartogh, M. Hoeft, R. P. Keenan, Y. Kovalev, K. Kreckel , et al. (66 additional authors not shown)

Abstract: The Next Generation Very Large Array (ngVLA) is a planned radio interferometer providing unprecedented sensitivity at wavelengths between 21 cm and 3 mm. Its 263 antenna element array will be spatially distributed across North America to enable both superb low surface brightness recovery and sub-milliarcsecond angular resolution imaging. The project was developed by the international astronomy com… ▽ More The Next Generation Very Large Array (ngVLA) is a planned radio interferometer providing unprecedented sensitivity at wavelengths between 21 cm and 3 mm. Its 263 antenna element array will be spatially distributed across North America to enable both superb low surface brightness recovery and sub-milliarcsecond angular resolution imaging. The project was developed by the international astronomy community under the lead of the National Radio Astronomy Observatory (NRAO), and is anticipated to be built between 2027 and 2037. Two workshops have been held in 2022 and 2023 with the goal to discuss and consolidate the scientific interests in the ngVLA within the German astronomical community. This community paper constitutes a collection of 48 science ideas which the German community aims to pursue with the ngVLA in the 2030s. This is not a complete list and the ideas are not developed at the level of a "Science Book", such that the present document is mainly meant provide a basis for further discussion within the community. As such, additional contributions are welcome, and will be considered for inclusion in future revisions. △ Less

Submitted 18 June, 2024; v1 submitted 16 November, 2023; originally announced November 2023.

Comments: Version 2.0 (status June 18, 2024): 169 pages, comments and future contributions welcome [v2.0: 7 new science cases added, some minor revisions to other chapters]

MAGIC detection of GRB 201216C at $z=1.1$

Authors: 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, J. Bernete, A. Berti, J. Besenrieder, C. Bigongiari , et al. (195 additional authors not shown)

Abstract: Gamma-ray bursts (GRBs) are explosive transient events occurring at cosmological distances, releasing a large amount of energy as electromagnetic radiation over several energy bands. We report the detection of the long GRB~201216C by the MAGIC telescopes. The source is located at $z=1.1$ and thus it is the farthest one detected at very high energies. The emission above \SI{70}{\GeV} of GRB~201216C… ▽ More Gamma-ray bursts (GRBs) are explosive transient events occurring at cosmological distances, releasing a large amount of energy as electromagnetic radiation over several energy bands. We report the detection of the long GRB~201216C by the MAGIC telescopes. The source is located at $z=1.1$ and thus it is the farthest one detected at very high energies. The emission above \SI{70}{\GeV} of GRB~201216C is modelled together with multi-wavelength data within a synchrotron and synchrotron-self Compton (SSC) scenario. We find that SSC can explain the broadband data well from the optical to the very-high-energy band. For the late-time radio data, a different component is needed to account for the observed emission. Differently from previous GRBs detected in the very-high-energy range, the model for GRB~201216C strongly favors a wind-like medium. The model parameters have values similar to those found in past studies of the afterglows of GRBs detected up to GeV energies. △ Less

Submitted 10 October, 2023; originally announced October 2023.

Comments: 13 pages, 6 figures, 2 tables. Accepted for publication in Monthly Notices of the Royal Astronomical Society

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 527, Issue 3, January 2024, Pages 5856-5867

Multi-year characterisation of the broad-band emission from the intermittent extreme BL Lac 1ES~2344+514

Authors: 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, I. Batković, J. Baxter, J. Becerra González, E. Bernardini, J. Bernete, A. Berti, J. Besenrieder, C. Bigongiari, A. Biland, O. Blanch , et al. (210 additional authors not shown)

Abstract: The BL Lac 1ES 2344+514 is known for temporary extreme properties (e.g., a shift of the synchrotron SED peak energy $ν_{synch,p}$ above 1keV). While those extreme states were so far observed only during high flux levels, additional multi-year observing campaigns are required to achieve a coherent picture. Here, we report the longest investigation of the source from radio to VHE performed so far, f… ▽ More The BL Lac 1ES 2344+514 is known for temporary extreme properties (e.g., a shift of the synchrotron SED peak energy $ν_{synch,p}$ above 1keV). While those extreme states were so far observed only during high flux levels, additional multi-year observing campaigns are required to achieve a coherent picture. Here, we report the longest investigation of the source from radio to VHE performed so far, focusing on a systematic characterisation of the intermittent extreme states. While our results confirm that 1ES 2344+514 typically exhibits $ν_{synch,p}>$1keV during elevated flux periods, we also find periods where the extreme state coincides with low flux activity. A strong spectral variability thus happens in the quiescent state, and is likely caused by an increase of the electron acceleration efficiency without a change in the electron injection luminosity. We also report a strong X-ray flare (among the brightest for 1ES 2344+514) without a significant shift of $ν_{synch,p}$. During this particular flare, the X-ray spectrum is among the softest of the campaign. It unveils complexity in the spectral evolution, where the common harder-when-brighter trend observed in BL Lacs is violated. During a low and hard X-ray state, we find an excess of the UV flux with respect to an extrapolation of the X-ray spectrum to lower energies. This UV excess implies that at least two regions contribute significantly to the infrared/optical/ultraviolet/X-ray emission. Using the simultaneous MAGIC, XMM-Newton, NuSTAR, and AstroSat observations, we argue that a region possibly associated with the 10 GHz radio core may explain such an excess. Finally, we investigate a VHE flare, showing an absence of simultaneous variability in the 0.3-2keV band. Using a time-dependent leptonic modelling, we show that this behaviour, in contradiction to single-zone scenarios, can instead be explained by a two-component model. △ Less

Submitted 5 October, 2023; originally announced October 2023.

Comments: Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 682, A114 (2024)

Performance of the joint LST-1 and MAGIC observations evaluated with Crab Nebula data

Authors: H. Abe, K. Abe, S. Abe, V. A. Acciari, A. Aguasca-Cabot, I. Agudo, N. Alvarez Crespo, T. Aniello, S. Ansoldi, L. A. Antonelli, C. Aramo, A. Arbet-Engels, C. Arcaro, M. Artero, K. Asano, P. Aubert, D. Baack, A. Babić, A. Baktash, A. Bamba, A. Baquero Larriva, L. Baroncelli, U. Barres de Almeida, J. A. Barrio, I. Batković , et al. (344 additional authors not shown)

Abstract: Aims. LST-1, the prototype of the Large-Sized Telescope for the upcoming Cherenkov Telescope Array Observatory, is concluding its commissioning in Observatorio del Roque de los Muchachos on the island of La Palma. The proximity of LST-1 (Large-Sized Telescope 1) to the two MAGIC (Major Atmospheric Gamma Imaging Cherenkov) telescopes permits observations of the same gamma-ray events with both syste… ▽ More Aims. LST-1, the prototype of the Large-Sized Telescope for the upcoming Cherenkov Telescope Array Observatory, is concluding its commissioning in Observatorio del Roque de los Muchachos on the island of La Palma. The proximity of LST-1 (Large-Sized Telescope 1) to the two MAGIC (Major Atmospheric Gamma Imaging Cherenkov) telescopes permits observations of the same gamma-ray events with both systems. Methods. We describe the joint LST-1+MAGIC analysis pipeline and use simultaneous Crab Nebula observations and Monte Carlo simulations to assess the performance of the three-telescope system. The addition of the LST-1 telescope allows the recovery of events in which one of the MAGIC images is too dim to survive analysis quality cuts. Results. Thanks to the resulting increase in the collection area and stronger background rejection, we find a significant improvement in sensitivity, allowing the detection of 30% weaker fluxes in the energy range between 200 GeV and 3 TeV. The spectrum of the Crab Nebula, reconstructed in the energy range ~60 GeV to ~10 TeV, is in agreement with previous measurements. △ Less

Submitted 3 October, 2023; originally announced October 2023.

Comments: Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 680, A66 (2023)

Deep learning-based radiointerferometric imaging with GAN-aided training

Authors: F. Geyer, K. Schmidt, J. Kummer, M. Brüggen, H. W. Edler, D. Elsässer, F. Griese, A. Poggenpohl, L. Rustige, W. Rhode

Abstract: Radio interferometry invariably suffers from an incomplete coverage of the spatial Fourier space, which leads to imaging artifacts. The current state-of-the-art technique is to create an image by Fourier-transforming the incomplete visibility data and to clean the systematic effects originating from incomplete data in Fourier space. Previously, we have shown how super-resolution methods based on c… ▽ More Radio interferometry invariably suffers from an incomplete coverage of the spatial Fourier space, which leads to imaging artifacts. The current state-of-the-art technique is to create an image by Fourier-transforming the incomplete visibility data and to clean the systematic effects originating from incomplete data in Fourier space. Previously, we have shown how super-resolution methods based on convolutional neural networks can reconstruct sparse visibility data. Our previous work has suffered from a low realism of the training data. The aim of this work is to build a whole simulation chain for realistic radio sources that then leads to a vastly improved neural net for the reconstruction of missing visibilities. This method offers considerable improvements in terms of speed, automatization and reproducibility over the standard techniques. Here we generate large amounts of training data by creating images of radio galaxies with a generative adversarial network (GAN) that has been trained on radio survey data. Then, we applied the Radio Interferometer Measurement Equation (RIME) in order to simulate the measurement process of a radio interferometer. We show that our neural network can reconstruct faithfully images of realistic radio galaxies. The reconstructed images agree well with the original images in terms of the source area, integrated flux density, peak flux density, and the multi-scale structural similarity index. Finally, we show how the neural net can be adapted to estimate the uncertainties in the imaging process. △ Less

Submitted 26 July, 2023; originally announced July 2023.

Comments: Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 677, A167 (2023)

Improved $S$-factor of the $^{12}$C(p,$γ$)$^{13}$N reaction at $E\,=\,$320-620~keV and the 422~keV resonance

Authors: J. Skowronski, E. Masha, D. Piatti, M. Aliotta, H. Babu, D. Bemmerer, A. Boeltzig, R. Depalo, A. Caciolli, F. Cavanna, L. Csedreki, Z. Fülöp, G. Imbriani, D. Rapagnani, S. Rümmler, K. Schmidt, R. S. Sidhu, T. Szücs, S. Turkat, A. Yadav

Abstract: The 12C(p,γ)13N reaction is the onset process of both the CNO and Hot CNO cycles that drive massive star, Red and Asymptotic Giant Branch star and novae nucleosynthesis. The 12C(p,γ)13N rate affects the final abundances of the stable 12,13C nuclides, with ramifications for meteoritic carbon isotopic abundances and the s-process neutron source strength. Here, a new underground measurement of the 12… ▽ More The 12C(p,γ)13N reaction is the onset process of both the CNO and Hot CNO cycles that drive massive star, Red and Asymptotic Giant Branch star and novae nucleosynthesis. The 12C(p,γ)13N rate affects the final abundances of the stable 12,13C nuclides, with ramifications for meteoritic carbon isotopic abundances and the s-process neutron source strength. Here, a new underground measurement of the 12C(p,γ)13N cross-section is reported. The present data, obtained at the Felsenkeller shallow-underground laboratory in Dresden (Germany), encompass the 320-620 keV center of mass energy range to include the wide and poorly constrained E = 422 keV resonance that dominates the rate at high temperatures. This work S-factor results, lower than literature by 25%, are included in a new comprehensive R-matrix fit, and the energy of the 1+ first excited state of 13N is found to be 2369.6(4) keV, with radiative and proton width of 0.49(3) eV and 34.9(2) keV respectively. A new reaction rate, based on present R-matrix fit and extrapolation, is suggested. △ Less

Submitted 15 June, 2023; originally announced June 2023.

Comments: Submitted to Phys. Rev. C (Letter)

Journal ref: Phys. Rev. C 107, L062801 (2023)

The Magnificent Five Images of Supernova Refsdal: Time Delay and Magnification Measurements

Authors: Patrick L. Kelly, Steven Rodney, Tommaso Treu, Simon Birrer, Vivien Bonvin, Luc Dessart, Ryan J. Foley, Alexei V. Filippenko, Daniel Gilman, Saurabh Jha, Jens Hjorth, Kaisey Mandel, Martin Millon, Justin Pierel, Stephen Thorp, Adi Zitrin, Tom Broadhurst, Wenlei Chen, Jose M. Diego, Alan Dressler, Or Graur, Mathilde Jauzac, Matthew A. Malkan, Curtis McCully, Masamune Oguri , et al. (6 additional authors not shown)

Abstract: In late 2014, four images of Supernova (SN) "Refsdal," the first known example of a strongly lensed SN with multiple resolved images, were detected in the MACS J1149 galaxy-cluster field. Following the images' discovery, the SN was predicted to reappear within hundreds of days at a new position ~8 arcseconds away in the field. The observed reappearance in late 2015 makes it possible to carry out R… ▽ More In late 2014, four images of Supernova (SN) "Refsdal," the first known example of a strongly lensed SN with multiple resolved images, were detected in the MACS J1149 galaxy-cluster field. Following the images' discovery, the SN was predicted to reappear within hundreds of days at a new position ~8 arcseconds away in the field. The observed reappearance in late 2015 makes it possible to carry out Refsdal's (1964) original proposal to use a multiply imaged SN to measure the Hubble constant H0, since the time delay between appearances should vary inversely with H0. Moreover, the position, brightness, and timing of the reappearance enable a novel test of the blind predictions of galaxy-cluster models, which are typically constrained only by the positions of multiply imaged galaxies. We have developed a new photometry pipeline that uses DOLPHOT to measure the fluxes of the five images of SN Refsdal from difference images. We apply four separate techniques to perform a blind measurement of the relative time delays and magnification ratios (mu_i/mu_1) between the last image SX and the earlier images S1-S4. We measure the relative time delay of SX-S1 to be 376.0+5.6-5.5 days and the relative magnification to be 0.30+0.05-0.03. This corresponds to a 1.5% precision on the time delay and 17% precision for the magnification ratios, and includes uncertainties due to millilensing and microlensing. In an accompanying paper, we place initial and blind constraints on the value of the Hubble constant. △ Less

Submitted 10 May, 2023; originally announced May 2023.

Comments: Published in ApJ. Companion paper presenting H0 constraints published in Science (DOI: 10.1126/science.abh1322)

Constraints on the Hubble constant from Supernova Refsdal's reappearance

Authors: Patrick L. Kelly, Steven Rodney, Tommaso Treu, Masamune Oguri, Wenlei Chen, Adi Zitrin, Simon Birrer, Vivien Bonvin, Luc Dessart, Jose M. Diego, Alexei V. Filippenko, Ryan J. Foley, Daniel Gilman, Jens Hjorth, Mathilde Jauzac, Kaisey Mandel, Martin Millon, Justin Pierel, Keren Sharon, Stephen Thorp, Liliya Williams, Tom Broadhurst, Alan Dressler, Or Graur, Saurabh Jha , et al. (5 additional authors not shown)

Abstract: The gravitationally lensed Supernova Refsdal appeared in multiple images, produced through gravitational lensing by a massive foreground galaxy cluster. After the supernova appeared in 2014, lens models of the galaxy cluster predicted an additional image of the supernova would appear in 2015, which was subsequently observed. We use the time delays between the images to perform a blinded measuremen… ▽ More The gravitationally lensed Supernova Refsdal appeared in multiple images, produced through gravitational lensing by a massive foreground galaxy cluster. After the supernova appeared in 2014, lens models of the galaxy cluster predicted an additional image of the supernova would appear in 2015, which was subsequently observed. We use the time delays between the images to perform a blinded measurement of the expansion rate of the Universe, quantified by the Hubble constant (H0). Using eight cluster lens models, we infer H0 = 64.8 +4.4-4.3 km / s / Mpc, where Mpc is the megaparsec. Using the two models most consistent with the observations, we find H0 = 66.6 +4.1-3.3 km / s / Mpc. The observations are best reproduced by models that assign dark-matter halos to individual galaxies and the overall cluster. △ Less

Submitted 2 September, 2023; v1 submitted 10 May, 2023; originally announced May 2023.

Comments: Published in Science on May 11, 2023; this version updated to reflect minor edits to galley proofs. Companion paper presenting time-delay and relative magnification measurements published in ApJ (DOI: 10.3847/1538-4357/ac4ccb)

Journal ref: Science, Volume 380, Issue 6649, article id. abh1322, May 11, 2023

Search for Gamma-ray Spectral Lines from Dark Matter Annihilation up to 100 TeV towards the Galactic Center with MAGIC

Authors: MAGIC Collaboration, H. Abe, S. Abe, V. A. Acciari, 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, W. Bhattacharyya , et al. (188 additional authors not shown)

Abstract: Line-like features in TeV $γ$-rays constitute a ''smoking gun'' for TeV-scale particle dark matter and new physics. Probing the Galactic Center region with ground-based Cherenkov telescopes enables the search for TeV spectral features in immediate association with a dense dark matter reservoir at a sensitivity out of reach for satellite $γ$-ray detectors, and direct detection and collider experime… ▽ More Line-like features in TeV $γ$-rays constitute a ''smoking gun'' for TeV-scale particle dark matter and new physics. Probing the Galactic Center region with ground-based Cherenkov telescopes enables the search for TeV spectral features in immediate association with a dense dark matter reservoir at a sensitivity out of reach for satellite $γ$-ray detectors, and direct detection and collider experiments. We report on 223 hours of observations of the Galactic Center region with the MAGIC stereoscopic telescope system reaching $γ$-ray energies up to 100 TeV. We improved the sensitivity to spectral lines at high energies using large-zenith-angle observations and a novel background modeling method within a maximum-likelihood analysis in the energy domain. No line-like spectral feature is found in our analysis. Therefore, we constrain the cross section for dark matter annihilation into two photons to $\langle σv \rangle \lesssim 5 \times 10^{-28}\,\mathrm{cm^3\,s^{-1}}$ at 1 TeV and $\langle σv \rangle \lesssim 1 \times 10^{-25}\,\mathrm{cm^3\,s^{-1}}$ at 100 TeV, achieving the best limits to date for a dark matter mass above 20 TeV and a cuspy dark matter profile at the Galactic Center. Finally, we use the derived limits for both cuspy and cored dark matter profiles to constrain supersymmetric wino models. △ Less

Submitted 20 December, 2022; originally announced December 2022.

Comments: Accepted in Phys. Rev. Lett., includes Supplemental materials. 9+4 pages, 3+3 figures, 3 tables. Corresponding authors: T. Inada, D. Kerszberg, M. Hütten

Report number: KEK-TH-2487, KEK-Cosmo-0307

MAGIC observations provide compelling evidence of the hadronic multi-TeV emission from the putative PeVatron SNR G106.3+2.7

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. (192 additional authors not shown)

Abstract: The SNR G106.3+2.7, detected at 1--100 TeV energies by different $γ$-ray facilities, is one of the most promising PeVatron candidates. This SNR has a cometary shape which can be divided into a head and a tail region with different physical conditions. However, it is not identified in which region the 100 TeV emission is produced due to the limited position accuracy and/or angular resolution of exi… ▽ More The SNR G106.3+2.7, detected at 1--100 TeV energies by different $γ$-ray facilities, is one of the most promising PeVatron candidates. This SNR has a cometary shape which can be divided into a head and a tail region with different physical conditions. However, it is not identified in which region the 100 TeV emission is produced due to the limited position accuracy and/or angular resolution of existing observational data. Additionally, it remains unclear whether the origin of the $γ$-ray emission is leptonic or hadronic. With the better angular resolution provided by these new MAGIC data compared to earlier $γ$-ray datasets, we aim to reveal the acceleration site of PeV particles and the emission mechanism by resolving the SNR G106.3+2.7 with 0.1$^\circ$ resolution at TeV energies. We detected extended $γ$-ray emission spatially coincident with the radio continuum emission at the head and tail of SNR G106.3+2.7. The fact that we detected a significant $γ$-ray emission with energies above 6.0 TeV from the tail region only suggests that the emissions above 10 TeV, detected with air shower experiments (Milagro, HAWC, Tibet AS$γ$ and LHAASO), are emitted only from the SNR tail. Under this assumption, the multi-wavelength spectrum of the head region can be explained with either hadronic or leptonic models, while the leptonic model for the tail region is in contradiction with the emission above 10 TeV and X-rays. In contrast, the hadronic model could reproduce the observed spectrum at the tail by assuming a proton spectrum with a cutoff energy of $\sim 1$ PeV for the tail region. Such a high energy emission in this middle-aged SNR (4--10 kyr) can be explained by considering the scenario that protons escaping from the SNR in the past interact with surrounding dense gases at present. △ Less

Submitted 28 November, 2022; originally announced November 2022.

Comments: 13 pages, 7 figures, Accepted for publication in A&A

Journal ref: A&A 671, A12 (2023)

Long-term multi-wavelength study of 1ES 0647+250

Authors: MAGIC Collaboration, V. A. Acciari, 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. Becerra González, W. Bednarek, E. Bernardini, M. Bernardos, A. Berti, J. Besenrieder, W. Bhattacharyya, C. Bigongiari, A. Biland, O. Blanch , et al. (195 additional authors not shown)

Abstract: The BL Lac object 1ES 0647+250 is one of the few distant $γ$-ray emitting blazars detected at very high energies (VHE, $\gtrsim$100 GeV) during a non-flaring state. It was detected with the MAGIC telescopes during its low activity in the years 2009-2011, as well as during three flaring activities in the years 2014, 2019 and 2020, with the highest VHE flux in the latter epoch. An extensive multi-in… ▽ More The BL Lac object 1ES 0647+250 is one of the few distant $γ$-ray emitting blazars detected at very high energies (VHE, $\gtrsim$100 GeV) during a non-flaring state. It was detected with the MAGIC telescopes during its low activity in the years 2009-2011, as well as during three flaring activities in the years 2014, 2019 and 2020, with the highest VHE flux in the latter epoch. An extensive multi-instrument data set was collected within several coordinated observing campaigns throughout these years. We aim to characterise the long-term multi-band flux variability of 1ES 0647+250, as well as its broadband spectral energy distribution (SED) during four distinct activity states selected in four different epochs, in order to constrain the physical parameters of the blazar emission region under certain assumptions. We evaluate the variability and correlation of the emission in the different energy bands with the fractional variability and the Z-transformed Discrete Correlation Function, as well as its spectral evolution in X-rays and $γ$ rays. Owing to the controversy in the redshift measurements of 1ES 0647+250 reported in the literature, we also estimate its distance in an indirect manner through the comparison of the GeV and TeV spectra from simultaneous observations with Fermi-LAT and MAGIC during the strongest flaring activity detected to date. Moreover, we interpret the SEDs from the four distinct activity states within the framework of one-component and two-component leptonic models, proposing specific scenarios that are able to reproduce the available multi-instrument data. △ Less

Submitted 23 November, 2022; originally announced November 2022.

Comments: 20 pages, 7 figures. Accepted in A&A. Corresponding authors: Jorge Otero-Santos; Daniel Morcuende; Vandad Fallah Ramazani; Daniela Dorner; David Paneque (mailto: contact.magic@mpp.mpg.de)

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

A lower bound on intergalactic magnetic fields from time variability of 1ES 0229+200 from MAGIC and Fermi/LAT observations

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

Abstract: Extended and delayed emission around distant TeV sources induced by the effects of propagation of gamma rays through the intergalactic medium can be used for the measurement of the intergalactic magnetic field (IGMF). We search for delayed GeV emission from the hard-spectrum TeV blazar 1ES 0229+200 with the goal to detect or constrain the IGMF-dependent secondary flux generated during the propagat… ▽ More Extended and delayed emission around distant TeV sources induced by the effects of propagation of gamma rays through the intergalactic medium can be used for the measurement of the intergalactic magnetic field (IGMF). We search for delayed GeV emission from the hard-spectrum TeV blazar 1ES 0229+200 with the goal to detect or constrain the IGMF-dependent secondary flux generated during the propagation of TeV gamma rays through the intergalactic medium. We analyze the most recent MAGIC observations over a 5 year time span and complement them with historic data of the H.E.S.S. and VERITAS telescopes along with a 12-year long exposure of the Fermi/LAT telescope. We use them to trace source evolution in the GeV-TeV band over one-and-a-half decade in time. We use Monte Carlo simulations to predict the delayed secondary gamma-ray flux, modulated by the source variability, as revealed by TeV-band observations. We then compare these predictions for various assumed IGMF strengths to all available measurements of the gamma-ray flux evolution. We find that the source flux in the energy range above 200 GeV experiences variations around its average on the 14 years time span of observations. No evidence for the flux variability is found in 1-100 GeV energy range accessible to Fermi/LAT. Non-detection of variability due to delayed emission from electromagnetic cascade developing in the intergalactic medium imposes a lower bound of B>1.8e-17 G for long correlation length IGMF and B>1e-14 G for an IGMF of the cosmological origin. Though weaker than the one previously derived from the analysis of Fermi/LAT data, this bound is more robust, being based on a conservative intrinsic source spectrum estimate and accounting for the details of source variability in the TeV energy band. We discuss implications of this bound for cosmological magnetic fields which might explain the baryon asymmetry of the Universe. △ Less

Submitted 7 October, 2022; originally announced October 2022.

Comments: 10 pages, 5 figures, accepted to A&A. Corresponding authors: Ievgen Vovk, Paolo Da Vela (mailto:contact.magic@mpp.mpg.de) and Andrii Neronov (mailto:Andrii.Neronov@unige.ch)

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

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)

Measurement of Hadron Production in $π^-$-C Interactions at 158 and 350 GeV/c with NA61/SHINE at the CERN SPS

Authors: NA61/SHINE Collaboration, :, H. Adhikary, K. K. Allison, N. Amin, E. V. Andronov, T. Antićić, I. -C. Arsene, Y. Balkova, M. Baszczyk, D. Battaglia, S. Bhosale, A. Blondel, M. Bogomilov, Y. Bondar, N. Bostan, A. Brandin, A. Bravar, W. Bryliński, J. Brzychczyk, M. Buryakov, M. Ćirković, M. Csanad, J. Cybowska, T. Czopowicz , et al. (134 additional authors not shown)

Abstract: We present a measurement of the momentum spectra of $π^\pm$, K$^\pm$, p$^\pm$, $Λ$, $\barΛ$ and K$^{0}_{S}$ produced in interactions of negatively charged pions with carbon nuclei at beam momenta of 158 and 350 GeV/c. The total production cross sections are measured as well. The data were collected with the large-acceptance spectrometer of the fixed target experiment NA61/SHINE at the CERN SPS. Th… ▽ More We present a measurement of the momentum spectra of $π^\pm$, K$^\pm$, p$^\pm$, $Λ$, $\barΛ$ and K$^{0}_{S}$ produced in interactions of negatively charged pions with carbon nuclei at beam momenta of 158 and 350 GeV/c. The total production cross sections are measured as well. The data were collected with the large-acceptance spectrometer of the fixed target experiment NA61/SHINE at the CERN SPS. The obtained double-differential $p$-$p_T$ spectra provide a unique reference data set with unprecedented precision and large phase-space coverage to tune models used for the simulation of particle production in extensive air showers in which pions are the most numerous projectiles. △ Less

Submitted 21 September, 2022; originally announced September 2022.

Gamma-ray observations of MAXI J1820+070 during the 2018 outburst

Authors: H. Abe, S. Abe, V. A. Acciari, 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, W. Bhattacharyya, C. Bigongiari , et al. (418 additional authors not shown)

Abstract: MAXI J1820+070 is a low-mass X-ray binary with a black hole as a compact object. This binary underwent an exceptionally bright X-ray outburst from March to October 2018, showing evidence of a non-thermal particle population through its radio emission during this whole period. The combined results of 59.5 hours of observations of the MAXI J1820+070 outburst with the H.E.S.S., MAGIC and VERITAS expe… ▽ More MAXI J1820+070 is a low-mass X-ray binary with a black hole as a compact object. This binary underwent an exceptionally bright X-ray outburst from March to October 2018, showing evidence of a non-thermal particle population through its radio emission during this whole period. The combined results of 59.5 hours of observations of the MAXI J1820+070 outburst with the H.E.S.S., MAGIC and VERITAS experiments at energies above 200 GeV are presented, together with Fermi-LAT data between 0.1 and 500 GeV, and multiwavelength observations from radio to X-rays. Gamma-ray emission is not detected from MAXI J1820+070, but the obtained upper limits and the multiwavelength data allow us to put meaningful constraints on the source properties under reasonable assumptions regarding the non-thermal particle population and the jet synchrotron spectrum. In particular, it is possible to show that, if a high-energy gamma-ray emitting region is present during the hard state of the source, its predicted flux should be at most a factor of 20 below the obtained Fermi-LAT upper limits, and closer to them for magnetic fields significantly below equipartition. During the state transitions, under the plausible assumption that electrons are accelerated up to ~ 500 GeV, the multiwavelength data and the gamma-ray upper limits lead consistently to the conclusion that a potential high-energy and very-high-energy gamma-ray emitting region should be located at a distance from the black hole ranging between 10^11 and 10^13 cm. Similar outbursts from low-mass X-ray binaries might be detectable in the near future with upcoming instruments such as CTA. △ Less

Submitted 6 October, 2022; v1 submitted 20 September, 2022; originally announced September 2022.

Comments: 15 pages, 5 figures. Accepted for publication in MNRAS

Deep Learning-based Imaging in Radio Interferometry

Authors: Kevin Schmidt, Felix Geyer, Stefan Fröse, Paul-Simon Blomenkamp, Marcus Brüggen, Francesco de Gasperin, Dominik Elsässer, Wolfgang Rhode

Abstract: The sparse layouts of radio interferometers result in an incomplete sampling of the sky in Fourier space which leads to artifacts in the reconstructed images. Cleaning these systematic effects is essential for the scientific use of radiointerferometric images. Established reconstruction methods are often time-consuming, require expert-knowledge, and suffer from a lack of reproducibility. We have d… ▽ More The sparse layouts of radio interferometers result in an incomplete sampling of the sky in Fourier space which leads to artifacts in the reconstructed images. Cleaning these systematic effects is essential for the scientific use of radiointerferometric images. Established reconstruction methods are often time-consuming, require expert-knowledge, and suffer from a lack of reproducibility. We have developed a prototype Deep Learning-based method that generates reproducible images in an expedient fashion. To this end, we take advantage of the efficiency of Convolutional Neural Networks to reconstruct image data from incomplete information in Fourier space. The Neural Network architecture is inspired by super-resolution models that utilize residual blocks. Using simulated data of radio galaxies that are composed of Gaussian components we train Deep Learning models whose reconstruction capability is quantified using various measures. The reconstruction performance is evaluated on clean and noisy input data by comparing the resulting predictions with the true source images. We find that source angles and sizes are well reproduced, while the recovered fluxes show substantial scatter, albeit not worse than existing methods without fine-tuning. Finally, we propose more advanced approaches using Deep Learning that include uncertainty estimates and a concept to analyze larger images. △ Less

Submitted 22 March, 2022; originally announced March 2022.

Comments: Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 664, A134 (2022)

Proton acceleration in thermonuclear nova explosions revealed by gamma rays

Authors: MAGIC Collaboration, V. A. Acciari, S. Ansoldi, L. A. Antonelli, A. Arbet Engels, M. Artero, K. Asano, D. Baack, A. Babić, A. Baquero, U. Barres de Almeida, J. A. Barrio, I. Batković, J. Becerra González, W. Bednarek, L. Bellizzi, E. Bernardini, M. Bernardos, A. Berti, J. Besenrieder, W. Bhattacharyya, C. Bigongiari, A. Biland, O. Blanch, H. Bökenkamp , et al. (186 additional authors not shown)

Abstract: Classical novae are cataclysmic binary star systems in which the matter of a companion star is accreted on a white dwarf (WD). Accumulation of hydrogen in a layer eventually causes a thermonuclear explosion on the surface of the WD, brightening the WD to ~10^5 solar luminosities and triggering ejection of the accumulated matter.They provide extreme conditions required to accelerate particles, elec… ▽ More Classical novae are cataclysmic binary star systems in which the matter of a companion star is accreted on a white dwarf (WD). Accumulation of hydrogen in a layer eventually causes a thermonuclear explosion on the surface of the WD, brightening the WD to ~10^5 solar luminosities and triggering ejection of the accumulated matter.They provide extreme conditions required to accelerate particles, electrons or protons, to high energies. Here we present the detection of gamma rays by the MAGIC telescopes from the 2021 outburst of RS Ophiuchi (RS Oph), a recurrent nova with a red giant (RG) companion, that allowed us, for the first time, to accurately characterize the emission from a nova in the 60 GeV to 250 GeV energy range. The theoretical interpretation of the combined Fermi-LAT and MAGIC data suggests that protons are accelerated to hundreds of GeV in the nova shock. Such protons should create bubbles of enhanced Cosmic Ray density, on the order of 10 pc, from the recurrent novae. △ Less

Submitted 10 November, 2022; v1 submitted 15 February, 2022; originally announced February 2022.

Comments: THIS VERSION IS AFTER THE JOURNAL REVIEW

Journal ref: Nature Astronomy, Volume 6, p. 689-697 (2022)

Equivalent widths of Lyman $α$ emitters in MUSE-Wide and MUSE-Deep

Authors: J. Kerutt, L. Wisotzki, A. Verhamme, K. B. Schmidt, F. Leclercq, E. C. Herenz, T. Urrutia, T. Garel, T. Hashimoto, M. Maseda, J. Matthee, H. Kusakabe, J. Schaye, J. Richard, B. Guiderdoni, V. Mauerhofer, T. Nanayakkara, E. Vitte

Abstract: The aim of this study is to better understand the connection between the Lyman $α$ rest-frame equivalent width (EW$_0$) and spectral properties as well as ultraviolet (UV) continuum morphology by obtaining reliable EW$_0$ histograms for a statistical sample of galaxies and by assessing the fraction of objects with large equivalent widths. We used integral field spectroscopy from MUSE combined with… ▽ More The aim of this study is to better understand the connection between the Lyman $α$ rest-frame equivalent width (EW$_0$) and spectral properties as well as ultraviolet (UV) continuum morphology by obtaining reliable EW$_0$ histograms for a statistical sample of galaxies and by assessing the fraction of objects with large equivalent widths. We used integral field spectroscopy from MUSE combined with broad-band data from the Hubble Space Telescope (HST) to measure EW$_0$. We analysed the emission lines of $1920$ Lyman $α$ emitters (LAEs) detected in the full MUSE-Wide (one hour exposure time) and MUSE-Deep (ten hour exposure time) surveys and found UV continuum counterparts in archival HST data. We fitted the UV continuum photometric images using the Galfit software to gain morphological information on the rest-UV emission and fitted the spectra obtained from MUSE to determine the double peak fraction, asymmetry, full-width at half maximum, and flux of the Lyman $α$ line. The two surveys show different histograms of Lyman $α$ EW$_0$. In MUSE-Wide, $20\%$ of objects have EW$_0 > 240$ Å, while this fraction is only $11\%$ in MUSE-Deep and $\approx 16\%$ for the full sample. This includes objects without HST continuum counterparts (one-third of our sample), for which we give lower limits for EW$_0$. The object with the highest securely measured EW$_0$ has EW$_0=589 \pm 193$ Å (the highest lower limit being EW$_0=4464$ Å). We investigate the connection between EW$_0$ and Lyman $α$ spectral or UV continuum morphological properties. The survey depth has to be taken into account when studying EW$_0$ distributions. We find that in general, high EW$_0$ objects can have a wide range of spectral and UV morphological properties, which might reflect that the underlying causes for high EW$_0$ values are equally varied. (abridged) △ Less

Submitted 14 February, 2022; originally announced February 2022.

Comments: 28 pages, 21 + 1 figures, 7 + 1 tables, accepted for publication in A&A

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

Investigating the blazar TXS 0506+056 through sharp multi-wavelength eyes during 2017-2019

Authors: MAGIC Collaboration, V. A. Acciari, T. Aniello, S. Ansoldi, L. A. Antonelli, A. Arbet Engels, M. Artero, K. Asano, D. Baack, A. Babić, A. Baquero, U. Barres de Almeida, J. A. Barrio, I. Batković, J. Becerra González, W. Bednarek, E. Bernardini, M. Bernardos, A. Berti, J. Besenrieder, W. Bhattacharyya, C. Bigongiari, A. Biland, O. Blanch, H. Bökenkamp , et al. (192 additional authors not shown)

Abstract: The blazar TXS 0506+056 got into the spotlight of the astrophysical community in September 2017, when a high-energy neutrino detected by IceCube (IceCube-170922A) was associated at the 3 $σ$ level to a $γ$-ray flare from this source. This multi-messenger photon-neutrino association remains, as per today, the most significant one ever observed. TXS 0506+056 was a poorly studied object before the Ic… ▽ More The blazar TXS 0506+056 got into the spotlight of the astrophysical community in September 2017, when a high-energy neutrino detected by IceCube (IceCube-170922A) was associated at the 3 $σ$ level to a $γ$-ray flare from this source. This multi-messenger photon-neutrino association remains, as per today, the most significant one ever observed. TXS 0506+056 was a poorly studied object before the IceCube-170922A event. To better characterize its broad-band emission, we organized a multi-wavelength campaign lasting 16 months (November 2017 to February 2019), covering the radio-band (Metsähovi, OVRO), the optical/UV (ASAS-SN, KVA, REM, Swift/UVOT), the X-rays (Swift/XRT, NuSTAR), the high-energy $γ$ rays (Fermi/LAT) and the very-high-energy (VHE) $γ$ rays (MAGIC). In $γ$ rays, the behaviour of the source was significantly different from the 2017 one: MAGIC observations show the presence of flaring activity during December 2018, while the source only shows an excess at the 4$σ$ level during the rest of the campaign (74 hours of accumulated exposure); Fermi/LAT observations show several short (days-to-week timescale) flares, different from the long-term brightening of 2017. No significant flares are detected at lower energies. The radio light curve shows an increasing flux trend, not seen in other wavelengths. We model the multi-wavelength spectral energy distributions in a lepto-hadronic scenario, in which the hadronic emission emerges as Bethe-Heitler and pion-decay cascade in the X-rays and VHE $γ$ rays. According to the model presented here, the December 2018 $γ$-ray flare was connected to a neutrino emission that was too brief and not bright enough to be detected by current neutrino instruments. △ Less

Submitted 1 May, 2022; v1 submitted 5 February, 2022; originally announced February 2022.

Comments: 18 pages, 6 figures; V2 to match published version, typo fixed

Combined searches for dark matter in dwarf spheroidal galaxies observed with the MAGIC telescopes, including new data from Coma Berenices and Draco

Authors: MAGIC Collaboration, V. A. Acciari, S. Ansoldi, L. A. Antonelli, A. Arbet Engels, M. Artero, K. Asano, D. Baack, A. Babić, A. Baquero, U. Barres de Almeida, J. A. Barrio, I. Batković, J. Becerra González, W. Bednarek, L. Bellizzi, E. Bernardini, M. Bernardos, A. Berti, J. Besenrieder, W. Bhattacharyya, C. Bigongiari, A. Biland, O. Blanch, H. Bökenkamp , et al. (169 additional authors not shown)

Abstract: Milky Way dwarf spheroidal galaxies (dSphs) are among the best candidates to search for signals of dark matter annihilation with Imaging Atmospheric Cherenkov Telescopes, given their high mass-to-light ratios and the fact that they are free of astrophysical gamma-ray emitting sources. Since 2011, MAGIC has performed a multi-year observation program in search for Weakly Interacting Massive Particle… ▽ More Milky Way dwarf spheroidal galaxies (dSphs) are among the best candidates to search for signals of dark matter annihilation with Imaging Atmospheric Cherenkov Telescopes, given their high mass-to-light ratios and the fact that they are free of astrophysical gamma-ray emitting sources. Since 2011, MAGIC has performed a multi-year observation program in search for Weakly Interacting Massive Particles (WIMPs) in dSphs. Results on the observations of Segue 1 and Ursa Major II dSphs have already been published and include some of the most stringent upper limits (ULs) on the velocity-averaged cross-section $\langle σ_{\mathrm{ann}} v \rangle$ of WIMP annihilation from observations of dSphs. In this work, we report on the analyses of 52.1 h of data of Draco dSph and 49.5 h of Coma Berenices dSph observed with the MAGIC telescopes in 2018 and in 2019 respectively. No hint of a signal has been detected from either of these targets and new constraints on the $\langle σ_{\mathrm{ann}} v \rangle$ of WIMP candidates have been derived. In order to improve the sensitivity of the search and reduce the effect of the systematic uncertainties due to the $J$-factor estimates, we have combined the data of all dSphs observed with the MAGIC telescopes. Using 354.3 h of dSphs good quality data, 95 % CL ULs on $\langle σ_{\mathrm{ann}} v \rangle$ have been obtained for 9 annihilation channels. For most of the channels, these results reach values of the order of $10^{-24} $cm$^3$/s at ${\sim}1$ TeV and are the most stringent limits obtained with the MAGIC telescopes so far. △ Less

Submitted 29 November, 2021; originally announced November 2021.

Comments: 17 pages, 5 figures

Journal ref: Physics of the Dark Universe, Volume 35, March 2022, 100912

Multiwavelength study of the gravitationally lensed blazar QSO B0218+357 between 2016 and 2020

Authors: MAGIC Collaboration, V. A. Acciari, S. Ansoldi, L. A. Antonelli, A. Arbet Engels, M. Artero, K. Asano, D. Baack, A. Babić, A. Baquero, U. Barres de Almeida, J. A. Barrio, I. Batković, 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 , et al. (186 additional authors not shown)

Abstract: We report multiwavelength observations of the gravitationally lensed blazar QSO B0218+357 in 2016-2020. Optical, X-ray and GeV flares were detected. The contemporaneous MAGIC observations do not show significant very-high-energy (VHE, >= 100 GeV) gamma-ray emission. The lack of enhancement in radio emission measured by OVRO indicates the multi-zone nature of the emission from this object. We const… ▽ More We report multiwavelength observations of the gravitationally lensed blazar QSO B0218+357 in 2016-2020. Optical, X-ray and GeV flares were detected. The contemporaneous MAGIC observations do not show significant very-high-energy (VHE, >= 100 GeV) gamma-ray emission. The lack of enhancement in radio emission measured by OVRO indicates the multi-zone nature of the emission from this object. We constrain the VHE duty cycle of the source to be < 16 2014-like flares per year (95% confidence). For the first time for this source, a broadband low-state SED is constructed with a deep exposure up to the VHE range. A flux upper limit on the low-state VHE gamma-ray emission of an order of magnitude below that of the 2014 flare is determined. The X-ray data are used to fit the column density of (8.10 +- 0.93 stat ) x 10^21 cm^-2 of the dust in the lensing galaxy. VLBI observations show a clear radio core and jet components in both lensed images, yet no significant movement of the components is seen. The radio measurements are used to model the source-lens-observer geometry and determine the magnifications and time delays for both components. The quiescent emission is modeled with the high-energy bump explained as a combination of synchrotron-self-Compton and external Compton emission from a region located outside of the broad line region. The bulk of the low-energy emission is explained as originating from a tens-of-parsecs scale jet. △ Less

Submitted 25 November, 2021; originally announced November 2021.

Comments: 21 pages, 12 figures, 5 tables, accepted for publication in MNRAS

Observation of the gamma-ray binary HESS J0632+057 with the H.E.S.S., MAGIC, and VERITAS telescopes

Authors: C. B. Adams, W. Benbow, A. Brill, J. H. Buckley, M. Capasso, A. J. Chromey, M. Errando, A. Falcone, K. A. Farrell, Q. Feng, J P. Finley, G. Foote, L. Fortson, A. Furniss, A. Gent, G. H. Gillanders, C. Giuri, O. Gueta, D. Hanna, T. Hassan, O. Hervet, J. Holder, B. Hona, T. B. Humensky, W. Jin , et al. (387 additional authors not shown)

Abstract: The results of gamma-ray observations of the binary system HESS J0632+057 collected during 450 hours over 15 years, between 2004 and 2019, are presented. Data taken with the atmospheric Cherenkov telescopes H.E.S.S., MAGIC, and VERITAS at energies above 350 GeV were used together with observations at X-ray energies obtained with Swift-XRT, Chandra, XMM-Newton, NuSTAR, and Suzaku. Some of these obs… ▽ More The results of gamma-ray observations of the binary system HESS J0632+057 collected during 450 hours over 15 years, between 2004 and 2019, are presented. Data taken with the atmospheric Cherenkov telescopes H.E.S.S., MAGIC, and VERITAS at energies above 350 GeV were used together with observations at X-ray energies obtained with Swift-XRT, Chandra, XMM-Newton, NuSTAR, and Suzaku. Some of these observations were accompanied by measurements of the Hα emission line. A significant detection of the modulation of the VHE gamma-ray fluxes with a period of 316.7+-4.4 days is reported, consistent with the period of 317.3+-0.7 days obtained with a refined analysis of X-ray data. The analysis of data of four orbital cycles with dense observational coverage reveals short timescale variability, with flux-decay timescales of less than 20 days at very high energies. Flux variations observed over the time scale of several years indicate orbit-to-orbit variability. The analysis confirms the previously reported correlation of X-ray and gamma-ray emission from the system at very high significance, but can not find any correlation of optical Hα parameters with X-ray or gamma-ray energy fluxes in simultaneous observations. The key finding is that the emission of HESS J0632+057 in the X-ray and gamma-ray energy bands is highly variable on different time scales. The ratio of gamma-ray to X-ray flux shows the equality or even dominance of the gamma-ray energy range. This wealth of new data is interpreted taking into account the insufficient knowledge of the ephemeris of the system, and discussed in the context of results reported on other gamma-ray binary systems. △ Less

Submitted 24 September, 2021; originally announced September 2021.

Comments: Accepted for publication by The Astrophysical Journal

Search for Very High-Energy Emission from the millisecond pulsar PSR J0218+4232

Authors: V. A. Acciari, S. Ansoldi, L. A. Antonelli, A. Arbet Engels, M. Artero, K. Asano, D. Baack, A. Babić, A. Baquero, U. Barres de Almeida, J. A. Barrio, I. Batković, 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 , et al. (176 additional authors not shown)

Abstract: PSR J0218+4232 is one of the most energetic millisecond pulsars known and has long been considered as one of the best candidates for very high-energy (VHE; >100 GeV) gamma-ray emission. Using 11.5 years of Fermi Large Area Telescope (LAT) data between 100 MeV and 870 GeV, and ~90 hours of MAGIC observations in the 20 GeV to 20 TeV range, we have searched for the highest energy gamma-ray emission f… ▽ More PSR J0218+4232 is one of the most energetic millisecond pulsars known and has long been considered as one of the best candidates for very high-energy (VHE; >100 GeV) gamma-ray emission. Using 11.5 years of Fermi Large Area Telescope (LAT) data between 100 MeV and 870 GeV, and ~90 hours of MAGIC observations in the 20 GeV to 20 TeV range, we have searched for the highest energy gamma-ray emission from PSR J0218+4232. Based on the analysis of the LAT data, we find evidence for pulsed emission above 25 GeV, but see no evidence for emission above 100 GeV (VHE) with MAGIC. We present the results of searches for gamma-ray emission, along with theoretical modeling, to interpret the lack of VHE emission. We conclude that, based on the experimental observations and theoretical modeling, it will remain extremely challenging to detect VHE emission from PSR J0218+4232 with the current generation of Imaging Atmospheric Cherenkov Telescopes (IACTs), and maybe even with future ones, such as the Cherenkov Telescope Array (CTA). △ Less

Submitted 25 August, 2021; originally announced August 2021.

Comments: 19 pages, 9 figures, Accepted to ApJ

Recovery and analysis of rest-frame UV emission lines in 2052 galaxies observed with MUSE at $1.5 < z < 6.4$

Authors: K. B. Schmidt, J. Kerutt, L. Wisotzki, T. Urrutia, A. Feltre, M. V. Maseda, T. Nanayakkara, R. Bacon, L. A. Boogaard, S. Conseil, T. Contini, E. C. Herenz, W. Kollatschny, M. Krumpe, F. Leclercq, G. Mahler, J. Matthee, V. Mauerhofer, J. Richard, J. Schaye

Abstract: [Abbreviated] Rest-frame UV emission lines probe physical parameters of the emitting star-forming galaxies and their environments. The strongest main UV line, Ly$α$, has been instrumental in advancing the general knowledge of galaxy formation in the early universe. However, observing Ly$α$ emission becomes increasingly challenging at $z \gtrsim 6$ when the neutral hydrogen fraction of the CGM and… ▽ More [Abbreviated] Rest-frame UV emission lines probe physical parameters of the emitting star-forming galaxies and their environments. The strongest main UV line, Ly$α$, has been instrumental in advancing the general knowledge of galaxy formation in the early universe. However, observing Ly$α$ emission becomes increasingly challenging at $z \gtrsim 6$ when the neutral hydrogen fraction of the CGM and IGM increases. Secondary weaker UV emission lines provide important alternative methods for studying galaxy properties at high redshift. We present a large sample of rest-frame UV emission line sources at intermediate redshift for calibrating and exploring the connection between secondary UV lines and the emitting galaxies' physical properties and their Ly$α$ emission. The sample of 2052 emission line sources with $1.5 < z < 6.4$ was selected through untargeted source detection in three-dimensional MUSE data cubes. We searched optimally extracted 1D spectra of the full sample for UV emission features via emission line template matching, resulting in a sample of more than 100 rest-frame UV emission line detections. We show that the detection efficiency of (non-Ly$α$) UV emission lines increases with survey depth, and that the UV emission line strength often correlate with the strength of Ciii]. We measured the velocity offsets of resonant emission lines with respect to systemic tracers as well as the electron density and the gas-phase abundance. Lastly, using "PhotoIonization Model Probability Density Functions" we find that the UV line emitters generally have ionization parameter log10(U) $\approx$ -2.5 and metal mass fractions that scatter around Z $\approx$ 10$^{-2}$, that is Z $\approx$ 0.66Z$\odot$. Value-added catalogs of the full sample of MUSE objects studied in this work and a collection of UV line emitters from the literature are provided with this paper. △ Less

Submitted 3 August, 2021; originally announced August 2021.

Comments: Accepted by A&A. Contains 54 pages, 21 figures, and appendices. Catalogs described in Tables B.1 and C.1 will be made available with published paper

Journal ref: A&A 654, A80 (2021)

First detection of VHE gamma-ray emission from TXS~1515--273, study of its X-ray variability and spectral energy distribution

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

Abstract: We report here on the first multi-wavelength (MWL) campaign on the blazar TXS 1515-273, undertaken in 2019 and extending from radio to very-high-energy gamma rays (VHE). Up until now, this blazar had not been the subject of any detailed MWL observations. It has a rather hard photon index at GeV energies and was considered a candidate extreme high-synchrotronpeaked source. MAGIC observations result… ▽ More We report here on the first multi-wavelength (MWL) campaign on the blazar TXS 1515-273, undertaken in 2019 and extending from radio to very-high-energy gamma rays (VHE). Up until now, this blazar had not been the subject of any detailed MWL observations. It has a rather hard photon index at GeV energies and was considered a candidate extreme high-synchrotronpeaked source. MAGIC observations resulted in the first-time detection of the source in VHE with a statistical significance of 7.6$σ$. The average integral VHE flux of the source is 6 $\pm$ 1% of the Crab nebula flux above 400 GeV. X-ray coverage was provided by Swift-XRT, XMMNewton, and NuSTAR. The long continuous X-ray observations were separated by $\sim$ 9 h, both showing clear hour scale flares. In the XMM-Newton data, both the rise and decay timescales are longer in the soft X-ray than in the hard X-ray band, indicating the presence of a particle cooling regime. The X-ray variability timescales were used to constrain the size of the emission region and the strength of the magnetic field. The data allowed us to determine the synchrotron peak frequency and classify the source as a flaring high, but not extreme, synchrotron peaked object. Considering the constraints and variability patterns from the X-ray data, we model the broad-band spectral energy distribution. We applied a simple one-zone model, which could not reproduce the radio emission and the shape of the optical emission, and a two-component leptonic model with two interacting components, enabling us to reproduce the emission from radio to VHE band. △ Less

Submitted 20 July, 2021; originally announced July 2021.

The MUSE-Wide survey: Three-dimensional clustering analysis of Lyman-$α$ emitters at $3.3<z<6$

Authors: Y. Herrero Alonso, M. Krumpe, L. Wisotzki, T. Miyaji, T. Garel, K. B. Schmidt, C. Diener, T. Urrutia, J. Kerutt, E. C. Herenz, J. Schaye, G. Pezzulli, M. V. Maseda, L. Boogaard, J. Richard

Abstract: We present an analysis of the spatial clustering of 695 Ly$α$-emitting galaxies (LAE) in the MUSE-Wide survey. All objects have spectroscopically confirmed redshifts in the range $3.3<z<6$. We employ the K-estimator of Adelberger et al. (2005), adapted and optimized for our sample. We also explore the standard two-point correlation function approach, which is however less suited for a pencil-beam… ▽ More We present an analysis of the spatial clustering of 695 Ly$α$-emitting galaxies (LAE) in the MUSE-Wide survey. All objects have spectroscopically confirmed redshifts in the range $3.3<z<6$. We employ the K-estimator of Adelberger et al. (2005), adapted and optimized for our sample. We also explore the standard two-point correlation function approach, which is however less suited for a pencil-beam survey such as ours. The results from both approaches are consistent. We parametrize the clustering properties by, (i) modelling the clustering signal with a power law (PL), and (ii) adopting a Halo Occupation Distribution (HOD) model. Applying HOD modeling, we infer a large-scale bias of $b_{\rm{HOD}}=2.80^{+0.38}_{-0.38}$ at a median redshift of the number of galaxy pairs $\langle z_{\rm pair}\rangle\simeq3.82$, while the PL analysis results in $b_{\rm{PL}}=3.03^{+1.51}_{-0.52}$ ($r_0=3.60^{+3.10}_{-0.90}\;h^{-1}$Mpc and $γ=1.30^{+0.36}_{-0.45}$). The implied typical dark matter halo (DMH) mass is $\log(M_{\rm{DMH}}/[h^{-1}\rm{M}_\odot])=11.34^{+0.23}_{-0.27}$. We study possible dependencies of the clustering signal on object properties by bisecting the sample into disjoint subsets, considering Ly$α$ luminosity, UV absolute magnitude, Ly$α$ equivalent width, and redshift as variables. We find a suggestive trend of more luminous Ly$α$ emitters residing in more massive DMHs than their lower Ly$α$ luminosity counterparts. We also compare our results to mock LAE catalogs based on a semi-analytic model of galaxy formation and find a stronger clustering signal than in our observed sample. By adopting a galaxy-conserving model we estimate that the LAEs in the MUSE-Wide survey will typically evolve into galaxies hosted by halos of $\log(M_{\rm{DMH}}/[h^{-1}\rm{M}_\odot])\approx13.5$ at redshift zero, suggesting that we observe the ancestors of present-day galaxy groups. △ Less

Submitted 8 July, 2021; originally announced July 2021.

Comments: Accepted for publication in A&A. 22 pages, 20 figures, 4 tables

Journal ref: A&A 653, A136 (2021)

Disentangling the physical origin of emission line ratio offsets at high redshift with spatially resolved spectroscopy

Authors: Jessie Hirtenstein, Tucker Jones, Ryan L. Sanders, Crystal L. Martin, M. C. Cooper, Gabriel Brammer, Tommaso Treu, Kasper Schmidt, Alice Shapley

Abstract: We present spatially resolved Hubble Space Telescope grism spectroscopy of 15 galaxies at $z\sim0.8$ drawn from the DEEP2 survey. We analyze H$α$+[N II], [S II] and [S III] emission on kpc scales to explore which mechanisms are powering emission lines at high redshifts, testing which processes may be responsible for the well-known offset of high redshift galaxies from the $z\sim0$ locus in the [O… ▽ More We present spatially resolved Hubble Space Telescope grism spectroscopy of 15 galaxies at $z\sim0.8$ drawn from the DEEP2 survey. We analyze H$α$+[N II], [S II] and [S III] emission on kpc scales to explore which mechanisms are powering emission lines at high redshifts, testing which processes may be responsible for the well-known offset of high redshift galaxies from the $z\sim0$ locus in the [O III]/H$β$ versus [N II]/H$α$ BPT (Baldwin-Phillips-Terlevich) excitation diagram. We study spatially resolved emission line maps to examine evidence for active galactic nuclei (AGN), shocks, diffuse ionized gas (DIG), or escaping ionizing radiation, all of which may contribute to the BPT offsets observed in our sample. We do not find significant evidence of AGN in our sample and quantify that, on average, AGN would need to contribute $\sim$25% of the H$α$ flux in the central resolution element in order to cause the observed BPT offsets. We find weak ($2σ$) evidence of DIG emission at low surface brightnesses, yielding an implied total DIG emission fraction of $\sim$20%, which is not significant enough to be the dominant emission line driver in our sample. In general we find that the observed emission is dominated by star forming H II regions. We discuss trends with demographic properties and the possible role of $α$-enhanced abundance patterns in the emission spectra of high redshift galaxies. Our results indicate that photo-ionization modeling with stellar population synthesis inputs is a valid tool to explore the specific star formation properties which may cause BPT offsets, to be explored in future work. △ Less

Submitted 17 August, 2021; v1 submitted 25 June, 2021; originally announced June 2021.

Comments: 27 pages, 26 figures. Accepted to ApJ

Investigation of the correlation patterns and the Compton dominance variability of Mrk 421 in 2017

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

Abstract: We present a detailed characterisation and theoretical interpretation of the broadband emission of the paradigmatic TeV blazar Mrk 421, with special focus on the multi-band flux correlations. The dataset has been collected through an extensive multiwavelength campaign organised between 2016 December and 2017 June. The instruments involved are MAGIC, FACT, Fermi-LAT, Swift, GASP-WEBT, OVRO, Medicin… ▽ More We present a detailed characterisation and theoretical interpretation of the broadband emission of the paradigmatic TeV blazar Mrk 421, with special focus on the multi-band flux correlations. The dataset has been collected through an extensive multiwavelength campaign organised between 2016 December and 2017 June. The instruments involved are MAGIC, FACT, Fermi-LAT, Swift, GASP-WEBT, OVRO, Medicina and Metsähovi. Additionally, four deep exposures (several hours long) with simultaneous MAGIC and NuSTAR observations allowed a precise measurement of the falling segments of the two spectral components. The very-high-energy (VHE; E > 100 GeV) gamma rays and X-rays are positively correlated at zero time lag, but the strength and characteristics of the correlation change substantially across the various energy bands probed. The VHE versus X-ray fluxes follow different patterns, partly due to substantial changes in the Compton dominance during a few days without a simultaneous increase in the X-ray flux (i.e. orphan gamma-ray activity). Studying the broadband spectral energy distribution (SED) during the days including NuSTAR observations, we show that these changes can be explained within a one-zone leptonic model with a blob that increases its size over time. Our multi-band correlation study also hints at an anti-correlation between UV/optical and X-ray at a significance higher than 3 sigmas. A VHE flare observed on 2017 February 4 shows gamma-ray variability on multi-hour timescales, with a factor 10 increase in the TeV flux but only a moderate increase in the keV flux. The related broadband SED is better described by a two-zone leptonic scenario rather than by a one-zone scenario. We find that the flare can be produced by the appearance of a compact second blob populated by high energetic electrons spanning a narrow range of Lorentz factors. △ Less

Submitted 10 June, 2021; originally announced June 2021.

Comments: 38 pages, 23 figures, accepted for publication in Astronomy and Astrophysics. Corresponding authors: Axel Arbet-Engels, David Paneque

Journal ref: A&A 655, A89 (2021)

The MUSE Extremely Deep Field: the Cosmic Web in Emission at High Redshift

Authors: Roland Bacon, David Mary, Thibault Garel, Jeremy Blaizot, Michael Maseda, Joop Schaye, Lutz Wisotzki, Simon Conseil, Jarle Brinchmann, Floriane Leclercq, Valentina Abril-Melgarejo, Leindert Boogaard, Nicolas Bouché, Thierry Contini, Anna Feltre, Bruno Guiderdoni, Christian Herenz, Wolfram Kollatschny, Haruka Kusakabe, Jorryt Matthee, Léo Michel-Dansac, Themiya Nanayakkara, Johan Richard, Martin Roth, Kasper B. Schmidt , et al. (7 additional authors not shown)

Abstract: We report the discovery of diffuse extended Ly-alpha emission from redshift 3.1 to 4.5, tracing cosmic web filaments on scales of 2.5-4 comoving Mpc. These structures have been observed in overdensities of Ly-alpha emitters in the MUSE Extremely Deep Field, a 140 hour deep MUSE observation located in the Hubble Ultra Deep Field. Among the 22 overdense regions identified, 5 are likely to harbor ver… ▽ More We report the discovery of diffuse extended Ly-alpha emission from redshift 3.1 to 4.5, tracing cosmic web filaments on scales of 2.5-4 comoving Mpc. These structures have been observed in overdensities of Ly-alpha emitters in the MUSE Extremely Deep Field, a 140 hour deep MUSE observation located in the Hubble Ultra Deep Field. Among the 22 overdense regions identified, 5 are likely to harbor very extended Ly-alpha emission at high significance with an average surface brightness of $\mathrm{5 \times 10^{-20} erg s^{-1} cm^{-2} arcsec^{-2}}$. Remarkably, 70% of the total Ly-alpha luminosity from these filaments comes from beyond the circumgalactic medium of any identified Ly-alpha emitters. Fluorescent Ly-alpha emission powered by the cosmic UV background can only account for less than 34% of this emission at z$\approx$3 and for not more than 10% at higher redshift. We find that the bulk of this diffuse emission can be reproduced by the unresolved Ly-alpha emission of a large population of ultra low luminosity Ly-alpha emitters ($\mathrm{<10^{40} erg s^{-1}}$), provided that the faint end of the Ly-alpha luminosity function is steep ($α\lessapprox -1.8$), it extends down to luminosities lower than $\mathrm{10^{38} - 10^{37} erg s^{-1}}$ and the clustering of these Ly-alpha emitters is significant (filling factor $< 1/6$). If these Ly-alpha emitters are powered by star formation, then this implies their luminosity function needs to extend down to star formation rates $\mathrm{< 10^{-4} M_\odot yr^{-1}}$. These observations provide the first detection of the cosmic web in Ly-alpha emission in typical filamentary environments and the first observational clue for the existence of a large population of ultra low luminosity Ly-alpha emitters at high redshift. △ Less

Submitted 10 February, 2021; originally announced February 2021.

Comments: 28 pages, 19 figures

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

Improving $z\sim7-11$ Galaxy Property Estimates with JWST/NIRCam Medium-Band Photometry

Authors: Guido Roberts-Borsani, Tommaso Treu, Charlotte Mason, Kasper B. Schmidt, Tucker Jones, Adriano Fontana

Abstract: The past decade has seen impressive progress in the detection of $z>7$ galaxies with the Hubble Space Telescope, however little is known about their properties. The James Webb Space Telescope will revolutionise the high-$z$ field by providing NIR (i.e., rest-frame optical) data of unprecedented depth and spatial resolution. Measuring galaxy quantities such as resolved stellar ages or gas metallici… ▽ More The past decade has seen impressive progress in the detection of $z>7$ galaxies with the Hubble Space Telescope, however little is known about their properties. The James Webb Space Telescope will revolutionise the high-$z$ field by providing NIR (i.e., rest-frame optical) data of unprecedented depth and spatial resolution. Measuring galaxy quantities such as resolved stellar ages or gas metallicity gradients traditionally requires spectroscopy, as broad-band imaging filters are generally too coarse to fully isolate diagnostics such as the 4000 Å (rest-frame) break, continuum emission from aged stars, and key emission lines (e.g., [OII], [OIII], H$β$). However, in this paper, we show that adding NIRCam images through a strategically chosen medium-band filter to common wide-band filters sets adopted by ERS and GTO programs delivers tighter constraints on these galactic properties. To constrain the choice of filter, we perform a systematic investigation of which combinations of wide-band filters from ERS and GTO programs and single medium-band filters offer the tightest constraints on several galaxy properties at redshifts $z\sim7-11$. We employ the JAGUAR extragalactic catalogs to construct statistical samples of physically-motivated mock photometry and conduct SED-fitting procedures to evaluate the accuracy of galaxy property (and photo-$z$) recovery with a simple star-formation history model. We find that adding $>4.1 μ$m medium filters at comparable depth to the broad-band filters can significantly improve photo-$z$s and yield close to order-of-magnitude improvements in the determination of quantities such as stellar ages, metallicities, SF-related quantities and emission line fluxes at $z\sim8$. For resolved sources, the proposed approach enables spatially-resolved determination of these quantities that would be prohibitive with slit spectroscopy. △ Less

Submitted 8 February, 2021; originally announced February 2021.

Comments: 13 pages, 7 figures. Accepted for publication in ApJ

Ultra-faint [CII] emission in a redshift = 2 gravitationally-lensed metal-poor dwarf galaxy

Authors: M. Rybak, E. da Cunha, B. Groves, J. A. Hodge, M. Aravena, M. Maseda, L. Boogaard, D. Berg, S. Charlot, R. Decarli, D. K. Erb, E. Nelson, C. Pacifici, K. B. Schmidt, F. Walter, A. van der Wel

Abstract: Extreme emission-line galaxies (EELGs) at redshift z=1-2 provide a unique view of metal-poor, starburst sources that are the likely drivers of the cosmic reionization at z$\geq6$. However, the molecular gas reservoirs of EELGs - the fuel for their intense star-formation - remain beyond the reach of current facilities. We present ALMA [CII] and PdBI CO(2-1) observations of a z=1.8, strongly lensed… ▽ More Extreme emission-line galaxies (EELGs) at redshift z=1-2 provide a unique view of metal-poor, starburst sources that are the likely drivers of the cosmic reionization at z$\geq6$. However, the molecular gas reservoirs of EELGs - the fuel for their intense star-formation - remain beyond the reach of current facilities. We present ALMA [CII] and PdBI CO(2-1) observations of a z=1.8, strongly lensed EELG SL2S 0217, a bright Lyman-$α$ emitter with a metallicity 0.05 $Z_\odot$. We obtain a tentative (3-4$σ$) detection of the [CII] line and set an upper limit on the [CII]/SFR ratio of $\leq1\times10^6$ $L_\odot$/($M_\odot$ yr$^{-1}$), based on the synthesized images and visibility-plane analysis. The CO(2-1) emission is not detected. Photoionization modelling indicates that up to 80% of the [CII] emission originates from neutral or molecular gas, although we can not rule out that the gas is fully ionized. The very faint [CII] emission is in line with both nearby metal-poor dwarfs and high-redshift Lyman-$α$ emitters, and predictions from hydrodynamical simulations. However, the [CII] line is 30$\times$ fainter than predicted by the De Looze et al. [CII]-SFR relation for local dwarfs, illustrating the danger of extrapolating locally-calibrated relations to high-redshift, metal-poor galaxies. △ Less

Submitted 28 March, 2021; v1 submitted 4 January, 2021; originally announced January 2021.

Comments: Corrected to match the published version. We thank Y. Harikane and I. Smail for their comments on the previous arXiv version

VHE gamma-ray detection of FSRQ QSO B1420+326 and modeling of its enhanced broadband state in 2020

Authors: V. A. Acciari, S. Ansoldi, L. A. Antonelli, A. Arbet Engels, M. Artero, 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. (209 additional authors not shown)

Abstract: Context. QSO B1420+326 is a blazar classified as a Flat Spectrum Radio Quasar (FSRQ). In the beginning of 2020 it underwent an enhanced flux state. An extensive multiwavelength campaign allowed us to trace the evolution of the flare. Aims. We search for VHE gamma-ray emission from QSO B1420+326 during this flaring state. We aim to characterize and model the broadband emission of the source over di… ▽ More Context. QSO B1420+326 is a blazar classified as a Flat Spectrum Radio Quasar (FSRQ). In the beginning of 2020 it underwent an enhanced flux state. An extensive multiwavelength campaign allowed us to trace the evolution of the flare. Aims. We search for VHE gamma-ray emission from QSO B1420+326 during this flaring state. We aim to characterize and model the broadband emission of the source over different phases of the flare. Methods. The source was observed with a number of instruments in radio, near infrared, optical (including polarimetry and spectroscopy), ultra-violet, X-ray and gamma-ray bands. We use dedicated optical spectroscopy results to estimate the accretion disk and the dust torus luminosity. We perform spectral energy distribution modeling in the framework of combined Synchrotron-Self-Compton and External Compton scenario in which the electron energy distribution is partially determined from acceleration and cooling processes. Results. During the enhanced state the flux of both SED components drastically increased and the peaks were shifted to higher energies. Follow up observations with the MAGIC telescopes led to the detection of very-high-energy gamma-ray emission from this source, making it one of only a handful of FSRQs known in this energy range. Modeling allows us to constrain the evolution of the magnetic field and electron energy distribution in the emission region. The gamma-ray flare was accompanied by a rotation of the optical polarization vector during a low polarization state. Also, a new, superluminal radio knot contemporaneously appeared in the radio image of the jet. The optical spectroscopy shows a prominent FeII bump with flux evolving together with the continuum emission and a MgII line with varying equivalent width. △ Less

Submitted 21 December, 2020; originally announced December 2020.

Comments: 20 pages, 15 figures, 7 tables, accepted for publication in A&A

Journal ref: A&A 647, A163 (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)

MAGIC observations of the nearby short gamma-ray burst GRB 160821B

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. (166 additional authors not shown)

Abstract: The coincident detection of GW170817 in gravitational waves and electromagnetic radiation spanning the radio to MeV gamma-ray bands provided the first direct evidence that short gamma-ray bursts (GRBs) can originate from binary neutron star (BNS) mergers. On the other hand, the properties of short GRBs in high-energy gamma rays are still poorly constrained, with only $\sim$20 events detected in th… ▽ More The coincident detection of GW170817 in gravitational waves and electromagnetic radiation spanning the radio to MeV gamma-ray bands provided the first direct evidence that short gamma-ray bursts (GRBs) can originate from binary neutron star (BNS) mergers. On the other hand, the properties of short GRBs in high-energy gamma rays are still poorly constrained, with only $\sim$20 events detected in the GeV band, and none in the TeV band. GRB~160821B is one of the nearest short GRBs known at $z=0.162$. Recent analyses of the multiwavelength observational data of its afterglow emission revealed an optical-infrared kilonova component, characteristic of heavy-element nucleosynthesis in a BNS merger. Aiming to better clarify the nature of short GRBs, this burst was automatically followed up with the MAGIC telescopes, starting from 24 seconds after the burst trigger. Evidence of a gamma-ray signal is found above $\sim$0.5 TeV at a significance of $\sim3\,σ$ during observations that lasted until 4 hours after the burst. Assuming that the observed excess events correspond to gamma-ray emission from GRB 160821B, in conjunction with data at other wavelengths, we investigate its origin in the framework of GRB afterglow models. The simplest interpretation with one-zone models of synchrotron-self-Compton emission from the external forward shock has difficulty accounting for the putative TeV flux. Alternative scenarios are discussed where the TeV emission can be relatively enhanced. The role of future GeV-TeV observations of short GRBs in advancing our understanding of BNS mergers and related topics is briefly addressed. △ Less

Submitted 13 December, 2020; originally announced December 2020.

Comments: 13 pages, 4 figures, Accepted for publication in ApJ. Corresponding authors: K. Noda, L. Nava, and S. Inoue

Journal ref: ApJ 908, 90 (2021)

MusE GAs Flow and Wind (MEGAFLOW) VI. A study of CIV and MgII absorbing gas surrounding [OII] emitting galaxies

Authors: Ilane Schroetter, Nicolas F. Bouché, Johannes Zabl, Hadi Rahmani, Martin Wendt, Sowgat Muzahid, Thierry Contini, Joop Schaye, Kasper B. Schmidt, Lutz Wisotzki

Abstract: Using the MEGAFLOW survey, which consists of a combination of MUSE and UVES observations of 22 quasar fields selected to contain strong MgII absorbers, we measure covering fractions of CIV and MgII as a function of impact parameter $b$ using a novel Bayesian logistic regression method on unbinned data, appropriate for small samples. We also analyse how the CIV and MgII covering fractions evolve wi… ▽ More Using the MEGAFLOW survey, which consists of a combination of MUSE and UVES observations of 22 quasar fields selected to contain strong MgII absorbers, we measure covering fractions of CIV and MgII as a function of impact parameter $b$ using a novel Bayesian logistic regression method on unbinned data, appropriate for small samples. We also analyse how the CIV and MgII covering fractions evolve with redshift. In the MUSE data, we found 215 $z=1-1.5$ [OII] emitters with fluxes $>10^{-17}$ erg\,s$^{-1}$\,cm$^{-2}$ and within 250 kpc of quasar sight-lines. Over this redshift path $z=1-1.5$, we have 19 (32) CIV (MgII) absorption systems with rest-frame equivalent width (REW) $W_r>$0.05Å associated with at least one [OII] emitter. The covering fractions of $z\approx1.2$ CIV (MgII) absorbers with mean $W_r\approx$0.7Å (1.0Å), exceeds 50\% within 23$^{+62}_{-16}$ (46$^{+18}_{-13}$) kpc. Together with published studies, our results suggest that the covering fraction of CIV (MgII) becomes larger (smaller) with time, respectively. For absorption systems that have CIV but not MgII, we find in 73\% of the cases no [OII] counterpart. This may indicate that the CIV comes from the intergalactic medium (IGM), i.e. beyond 250 kpc, or that it is associated with lower-mass or quiescent galaxies. △ Less

Submitted 4 June, 2021; v1 submitted 9 December, 2020; originally announced December 2020.

Comments: 9 pages, 9 figures, accepted for publication in MNRAS

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

Detection of the Geminga pulsar with MAGIC hints at a power-law tail emission beyond 15 GeV

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. (168 additional authors not shown)

Abstract: We report the detection of pulsed gamma-ray emission from the Geminga pulsar (PSR J0633+1746) between $15\,$GeV and $75\,$GeV. This is the first time a middle-aged pulsar has been detected up to these energies. Observations were carried out with the MAGIC telescopes between 2017 and 2019 using the low-energy threshold Sum-Trigger-II system. After quality selection cuts, $\sim 80\,$hours of observa… ▽ More We report the detection of pulsed gamma-ray emission from the Geminga pulsar (PSR J0633+1746) between $15\,$GeV and $75\,$GeV. This is the first time a middle-aged pulsar has been detected up to these energies. Observations were carried out with the MAGIC telescopes between 2017 and 2019 using the low-energy threshold Sum-Trigger-II system. After quality selection cuts, $\sim 80\,$hours of observational data were used for this analysis. To compare with the emission at lower energies below the sensitivity range of MAGIC, $11$ years of Fermi-LAT data above $100\,$MeV were also analysed. From the two pulses per rotation seen by Fermi-LAT, only the second one, P2, is detected in the MAGIC energy range, with a significance of $6.3\,σ$. The spectrum measured by MAGIC is well-represented by a simple power law of spectral index $Γ= 5.62\pm0.54$, which smoothly extends the Fermi-LAT spectrum. A joint fit to MAGIC and Fermi-LAT data rules out the existence of a sub-exponential cut-off in the combined energy range at the $3.6\,σ$ significance level. The power-law tail emission detected by MAGIC is interpreted as the transition from curvature radiation to Inverse Compton Scattering of particles accelerated in the northern outer gap. △ Less

Submitted 20 November, 2020; originally announced November 2020.

Comments: Corresponding authors: Marcos López-Moya, Giovanni Ceribella, Kouichi Hirotani and Thomas Schweizer

Journal ref: A&A 643, L14 (2020)

Velocity distribution of larger meteoroids and small asteroids impacting Earth

Authors: Esther Drolshagen, Theresa Ott, Detlef Koschny, Gerhard Drolshagen, Anna Kristiane Schmidt, Björn Poppe

Abstract: Various meteor and fireball networks exist worldwide. Most data sets which include ground-based observational data of meteors are affected by biases. The larger and faster the entering meteoroid, the brighter is the produced meteor. Hence, small and slow objects often stay undetected. This bias of meteor observations towards faster meteoroids is a challenge if quantitative population and flux mode… ▽ More Various meteor and fireball networks exist worldwide. Most data sets which include ground-based observational data of meteors are affected by biases. The larger and faster the entering meteoroid, the brighter is the produced meteor. Hence, small and slow objects often stay undetected. This bias of meteor observations towards faster meteoroids is a challenge if quantitative population and flux models are derived. In this work the velocity distribution of objects in space is analysed by using different data sets that are not affected by this velocity bias since they include only large objects, like the near-Earth object (NEO) risk list of ESA's (European Space Agency) SSA (Space Situational Awareness) near-Earth object Coordination Centre (NEOCC), and the fireballs in NASA's (National Aeronautics and Space Administration) CNEOS (Center for near-Earth object Studies) JPL (Jet Propulsion Laboratory) fireball database. Additionally, when only the largest objects recorded with the CILBO (Canary Island Long-Baseline Observatory) camera setup were analysed, a very similar distribution was shown. These velocity distributions are in good agreement with a widely used velocity distribution for smaller sporadic meteoroids in free space which was adopted as reference by the ECSS (European Cooperation for Space Standardisation) Space Environment Standard. △ Less

Submitted 16 November, 2020; originally announced November 2020.

Journal ref: Planetary and Space Science Volume 184, May 2020, 104869

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)

MUSE observations towards the lensing cluster A2744: Intersection between the LBG and LAE populations at z $\sim$ 3-7

Authors: G. de La Vieuville, R. Pelló, J. Richard, G. Mahler, L. Lévêque, F. E. Bauer, D. J. Lagattuta, J. Blaizot, T. Contini, L. Guaita, H. Kusakabe, N. Laporte, J. Martinez, M. V. Maseda, D. Schaerer, K. B. Schmidt, A. Verhamme

Abstract: We present a study of the intersection between the populations of star forming galaxies (SFGs) selected as either Lyman break galaxies (LBGs) or Lyman-alpha emitters (LAEs) in the redshift range 2.9 - 6.7, within the same volume of universe sampled by the Multi-Unit Spectroscopic Explorer (MUSE) behind the Hubble Frontier Fields lensing cluster A2744. We define three samples of star-forming galaxi… ▽ More We present a study of the intersection between the populations of star forming galaxies (SFGs) selected as either Lyman break galaxies (LBGs) or Lyman-alpha emitters (LAEs) in the redshift range 2.9 - 6.7, within the same volume of universe sampled by the Multi-Unit Spectroscopic Explorer (MUSE) behind the Hubble Frontier Fields lensing cluster A2744. We define three samples of star-forming galaxies: LBG galaxies with an LAE counterpart (92 galaxies), LBG galaxies without LAE counterpart (408 galaxies) and LAE galaxies without an LBG counterpart (46 galaxies). All these galaxies are intrinsically faint due to the lensing nature of the sample (Muv $\ge$ -20.5). The fraction of LAEs among all selected SFGs increases with redshift up to z $\sim$ 6 and decreases for higher redshifts. The evolution of LAE/LBG populations with UV magnitude and Lya luminosity shows that the LAE selection is able to identify intrinsically UV faint galaxies with Muv $\ge$ -15 that are typically missed in the deepest lensing photometric surveys. The LBG population seems to fairly represent the total population of SFGs down to Muv$\sim$-15. Galaxies with Muv$<-17$ tend to have SFRLya$<$SFRuv, whereas the opposite trend is observed within our sample for faint galaxies with Muv$>-17$, including galaxies only detected by their Lya emission, with a large scatter. These trends, previously observed in other samples of SFGs at high-$z$, are seen here for very faint Muv$\sim -15$ galaxies, much fainter than in previous studies. There is no clear evidence, based on the present results, for an intrinsic difference on the properties of the two populations selected as LBG and/or LAE. The observed trends could be explained by a combination of several facts, like the existence of different star-formation regimes, the dust content, the relative distribution and morphology of dust and stars, or the stellar populations △ Less

Submitted 15 November, 2020; v1 submitted 27 September, 2020; originally announced September 2020.

Journal ref: A&A 644, A39 (2020)

An Atlas of MUSE Observations towards Twelve Massive Lensing Clusters

Authors: Johan Richard, Adélaïde Claeyssens, David J. Lagattuta, Lucia Guaita, Franz E. Bauer, Roser Pello, David Carton, Roland Bacon, Geneviève Soucail, Gonzalo Prieto Lyon, Jean-Paul Kneib, Guillaume Mahler, Benjamin Clément, Wilfried Mercier, Andrei Variu, Amélie Tamone, Harald Ebeling, Kasper B. Schmidt, Themiya Nanayakkara, Michael Maseda, Peter M. Weilbacher, Nicolas Bouché, Rychard J. Bouwens, Lutz Wisotzki, Geoffroy de la Vieuville , et al. (3 additional authors not shown)

Abstract: Spectroscopic surveys of massive galaxy clusters reveal the properties of faint background galaxies, thanks to the magnification provided by strong gravitational lensing. We present a systematic analysis of integral-field-spectroscopy observations of 12 massive clusters, conducted with the Multi Unit Spectroscopic Explorer (MUSE). All data were taken under very good seeing conditions (0.6") in eff… ▽ More Spectroscopic surveys of massive galaxy clusters reveal the properties of faint background galaxies, thanks to the magnification provided by strong gravitational lensing. We present a systematic analysis of integral-field-spectroscopy observations of 12 massive clusters, conducted with the Multi Unit Spectroscopic Explorer (MUSE). All data were taken under very good seeing conditions (0.6") in effective exposure times between two and 15 hrs per pointing, for a total of 125 hrs. Our observations cover a total solid angle of ~23 arcmin$^2$ in the direction of clusters, many of which were previously studied by the MACS, Frontier Fields, GLASS and CLASH programs. The achieved emission line detection limit at 5$σ$ for a point source varies between (0.77--1.5)$\times$10$^{-18}$ erg\,s$^{-1}$\,cm$^{-2}$ at 7000Å. We present our developed strategy to reduce these observational data, detect sources and determine their redshifts. We construct robust mass models for each cluster to further confirm our redshift measurements using strong-lensing constraints, and identify a total of 312 strongly lensed sources producing 939 multiple images. The final redshift catalogs contain more than 3300 robust redshifts, of which 40\% are for cluster members and $\sim$30\% for lensed Lyman-$α$ emitters. 14\% of all sources are line emitters not seen in the available HST images, even at the depth of the FFs ($\sim29$ AB). We find that the magnification distribution of the lensed sources in the high-magnification regime ($μ{=}$ 2--25) follows the theoretical expectation of $N(z)\proptoμ^{-2}$. The quality of this dataset, number of lensed sources, and number of strong-lensing constraints enables detailed studies of the physical properties of both the lensing cluster and the background galaxies. The full data products from this work are made available to the community. [abridged] △ Less

Submitted 1 January, 2021; v1 submitted 21 September, 2020; originally announced September 2020.

Comments: 47 pages, 13 figures, 18 tables. Accepted for publication in A&A. Data Release products available on CDS and at https://cral-perso.univ-lyon1.fr/labo/perso/johan.richard/MUSE_data_release

Journal ref: A&A 646, A83 (2021)

Evaluation of Experimental Constraints on the $^{44}$Ti($α$,p)$^{47}$V Reaction Cross Section Relevant for Supernovae

Authors: Kelly A. Chipps, Philip Adsley, Manoel Couder, W. Raphael Hix, Zachary Meisel, Konrad Schmidt

Abstract: Due to its importance as an astronomical observable in core-collapse supernovae (CCSNe), the reactions producing and destroying $^{44}$Ti must be well constrained. Generally, statistical model calculations such as Hauser-Feshbach are employed when experimental cross sections are not available, but the variation in such adopted rates can be large. Here, data from the literature is compared with sta… ▽ More Due to its importance as an astronomical observable in core-collapse supernovae (CCSNe), the reactions producing and destroying $^{44}$Ti must be well constrained. Generally, statistical model calculations such as Hauser-Feshbach are employed when experimental cross sections are not available, but the variation in such adopted rates can be large. Here, data from the literature is compared with statistical model calculations of the $^{44}$Ti($α$,p)$^{47}$V reaction cross section and used to constrain the possible reaction rate variation over the temperatures relevant to CCSNe. Suggestions for targeted future measurements are given. △ Less

Submitted 30 November, 2020; v1 submitted 31 July, 2020; originally announced August 2020.

Comments: 7 pages, 4 figures. PRC, published

Journal ref: Phys. Rev. C 102, 035806 (2020)

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