Showing 1–50 of 282 results for author: López, G

A detailed study of the very-high-energy Crab pulsar emission with the LST-1

Authors: CTA-LST Project, :, K. Abe, S. Abe, A. Abhishek, F. Acero, A. Aguasca-Cabot, I. Agudo, N. Alvarez Crespo, L. A. Antonelli, C. Aramo, A. Arbet-Engels, C. Arcaro, M. Artero, K. Asano, P. Aubert, A. Baktash, A. Bamba, A. Baquero Larriva, L. Baroncelli, U. Barres de Almeida, J. A. Barrio, I. Batkovic, J. Baxter, J. Becerra González , et al. (272 additional authors not shown)

Abstract: Context: There are currently three pulsars firmly detected by imaging atmospheric Cherenkov telescopes (IACTs), two of them reaching TeV energies, challenging models of very-high-energy (VHE) emission in pulsars. More precise observations are needed to better characterize pulsar emission at these energies. The LST-1 is the prototype of the Large-Sized Telescope, that will be part of the Cherenkov… ▽ More Context: There are currently three pulsars firmly detected by imaging atmospheric Cherenkov telescopes (IACTs), two of them reaching TeV energies, challenging models of very-high-energy (VHE) emission in pulsars. More precise observations are needed to better characterize pulsar emission at these energies. The LST-1 is the prototype of the Large-Sized Telescope, that will be part of the Cherenkov Telescope Array Observatory (CTAO). Its improved performance over previous IACTs makes it well suited for studying pulsars. Aims: To study the Crab pulsar emission with the LST-1, improving and complementing the results from other telescopes. These observations can also be used to characterize the potential of the LST-1 to study other pulsars and detect new ones. Methods: We analyzed a total of $\sim$103 hours of gamma-ray observations of the Crab pulsar conducted with the LST-1 in the period from September 2020 to January 2023. The observations were carried out at zenith angles less than 50 degrees. A new analysis of the Fermi-LAT data was also performed, including $\sim$14 years of observations. Results: The Crab pulsar phaseogram, long-term light-curve, and phase-resolved spectra are reconstructed with the LST-1 from 20 GeV to 450 GeV for P1 and up to 700 GeV for P2. The pulsed emission is detected with a significance of 15.2$σ$. The two characteristic emission peaks of the Crab pulsar are clearly detected (>10$σ$), as well as the so-called bridge emission (5.7$σ$). We find that both peaks are well described by power laws, with spectral indices of $\sim$3.44 and $\sim$3.03 respectively. The joint analysis of Fermi-LAT and LST-1 data shows a good agreement between both instruments in the overlapping energy range. The detailed results obtained in the first observations of the Crab pulsar with LST-1 show the potential that CTAO will have to study this type of sources. △ Less

Submitted 2 July, 2024; originally announced July 2024.

Comments: Accepted by A&A

MINDS. A multi-instrument investigation into the molecule-rich JWST-MIRI spectrum of the DF Tau binary system

Authors: Sierra L. Grant, Nicolas T. Kurtovic, Ewine F. van Dishoeck, Thomas Henning, Inga Kamp, Hugo Nowacki, Karine Perraut, Andrea Banzatti, Milou Temmink, Valentin Christiaens, Matthias Samland, Danny Gasman, Benoît Tabone, Manuel Güdel, Pierre-Olivier Lagage, Aditya M. Arabhavi, David Barrado, Alessio Caratti o Garatti, Adrian M. Glauser, Hyerin Jang, Jayatee Kanwar, Fred Lahuis, Maria Morales-Calderón, Göran Olofsson, Giulia Perotti , et al. (4 additional authors not shown)

Abstract: Most stars form in multiple systems whose properties can significantly impact circumstellar disk evolution. We investigate the physical and chemical properties of the equal-mass, small separation (~66 mas, ~9 au) DF Tau binary system. Previous observations indicated that only DF Tau A has a circumstellar disk. We present JWST-MIRI MRS observations of DF Tau. The MIRI spectrum shows a forest of H2O… ▽ More Most stars form in multiple systems whose properties can significantly impact circumstellar disk evolution. We investigate the physical and chemical properties of the equal-mass, small separation (~66 mas, ~9 au) DF Tau binary system. Previous observations indicated that only DF Tau A has a circumstellar disk. We present JWST-MIRI MRS observations of DF Tau. The MIRI spectrum shows a forest of H2O lines and emission from CO, C2H2, HCN, CO2, and OH. LTE slab models are used to determine the properties of the gas, and we analyze high angular spatial and spectral resolution data from ALMA, VLTI-GRAVITY, and IRTF-iSHELL to aid in the interpretation of the JWST data. The 1.3 mm ALMA continuum data show two equal-brightness sources of compact (R<3 au) emission, with separations and movement consistent with astrometry from VLTI-GRAVITY and the known orbit. This is interpreted as a robust detection of a disk around DF Tau B, which we suggest may host a small (~1 au) cavity to reconcile all observations. The disk around DF Tau A is expected to be a full disk, and spatially and spectrally resolved dust and gas emission points to hot, close-in (<0.2 au) material. Hot (~500-1000 K) H2O, HCN, and C2H2 emission in the MIRI data likely originate in the DF Tau A disk, while a cold (<200 K) H2O component with an extended emitting area is consistent with an origin from both disks. Despite the very compact outer disks, the inner disk composition and conditions are similar to isolated systems, suggesting that the close binary nature is not a driving factor in setting the inner disk chemistry. However, constraining the geometry of the disks, for instance, via higher resolution ALMA observations, would provide additional insight into the mid-infrared gas emission. JWST observations of spatially resolved binaries will be important for understanding the impact of binarity on inner disk chemistry more generally. △ Less

Submitted 11 July, 2024; v1 submitted 14 June, 2024; originally announced June 2024.

Comments: Accepted to A&A on June 15th, 2024. Calibrated ALMA data are available at https://zenodo.org/records/11215003. Reduced JWST-MIRI MRS spectra will be available at Spexodisks.com and available by request to the corresponding author

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.

High contrast at short separation with VLTI/GRAVITY: Bringing Gaia companions to light

Authors: N. Pourré, T. O. Winterhalder, J. -B. Le Bouquin, S. Lacour, A. Bidot, M. Nowak, A. -L. Maire, D. Mouillet, C. Babusiaux, J. Woillez, R. Abuter, A. Amorim, R. Asensio-Torres, W. O. Balmer, M. Benisty, J. -P. Berger, H. Beust, S. Blunt, A. Boccaletti, M. Bonnefoy, H. Bonnet, M. S. Bordoni, G. Bourdarot, W. Brandner, F. Cantalloube , et al. (151 additional authors not shown)

Abstract: Since 2019, GRAVITY has provided direct observations of giant planets and brown dwarfs at separations of down to 95 mas from the host star. Some of these observations have provided the first direct confirmation of companions previously detected by indirect techniques (astrometry and radial velocities). We want to improve the observing strategy and data reduction in order to lower the inner working… ▽ More Since 2019, GRAVITY has provided direct observations of giant planets and brown dwarfs at separations of down to 95 mas from the host star. Some of these observations have provided the first direct confirmation of companions previously detected by indirect techniques (astrometry and radial velocities). We want to improve the observing strategy and data reduction in order to lower the inner working angle of GRAVITY in dual-field on-axis mode. We also want to determine the current limitations of the instrument when observing faint companions with separations in the 30-150 mas range. To improve the inner working angle, we propose a fiber off-pointing strategy during the observations to maximize the ratio of companion-light-to-star-light coupling in the science fiber. We also tested a lower-order model for speckles to decouple the companion light from the star light. We then evaluated the detection limits of GRAVITY using planet injection and retrieval in representative archival data. We compare our results to theoretical expectations. We validate our observing and data-reduction strategy with on-sky observations; first in the context of brown dwarf follow-up on the auxiliary telescopes with HD 984 B, and second with the first confirmation of a substellar candidate around the star Gaia DR3 2728129004119806464. With synthetic companion injection, we demonstrate that the instrument can detect companions down to a contrast of $8\times 10^{-4}$ ($Δ\mathrm{K}= 7.7$ mag) at a separation of 35 mas, and a contrast of $3\times 10^{-5}$ ($Δ\mathrm{K}= 11$ mag) at 100 mas from a bright primary (K<6.5), for 30 min exposure time. With its inner working angle and astrometric precision, GRAVITY has a unique reach in direct observation parameter space. This study demonstrates the promising synergies between GRAVITY and Gaia for the confirmation and characterization of substellar companions. △ Less

Submitted 6 June, 2024; originally announced June 2024.

Comments: 16 pages, 14 figures. Submitted to A&A

PENELLOPE\,VI. -- Searching the PENELLOPE/UVES sample with spectro-astrometry: Two new microjets of Sz 103 and XX Cha

Authors: T. Sperling, J. Eislöffel, C. F. Manara, J. Campbell-White, C. Schneider, A. Frasca, K. Maucó, M. Siwak, B. Fuhrmeister, R. Garcia Lopez

Abstract: The main goal of this study is to screen the PENELLOPE/UVES targets for outflow activity and find microjets via spectro-astrometry in, e.g., the [OI]$λ$6300 line. In total, 34 T\,Tauri stars of the PENELLOPE survey have been observed with the high resolution slit spectrograph UVES in three different slit positions rotated by $120^\text{o}$. Our spectro-astrometric analysis in the [OI]$λ$6300 wind… ▽ More The main goal of this study is to screen the PENELLOPE/UVES targets for outflow activity and find microjets via spectro-astrometry in, e.g., the [OI]$λ$6300 line. In total, 34 T\,Tauri stars of the PENELLOPE survey have been observed with the high resolution slit spectrograph UVES in three different slit positions rotated by $120^\text{o}$. Our spectro-astrometric analysis in the [OI]$λ$6300 wind line reveals two newly discovered microjets associated with Sz\,103 and XX\,Cha. Both microjets have an extent of about $0.04$ arcseconds, that is, $<10\,\text{au}$, and we confined their orientation by the three slit observations. Furthermore, we confirm the binary nature of VW\,Cha and CVSO\,109. We present (further) evidence that DK\,Tau\,B and CVSO\,104\,A are spectroscopic binaries. Sz\,115 is tentatively a spectroscopic binary. The origin of the LVC, that is, MHD winds versus photoevaporative winds, of the Sz\,103 and XX\,Cha microjets remains unclear. △ Less

Submitted 18 May, 2024; originally announced May 2024.

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

Four-of-a-kind? Comprehensive atmospheric characterisation of the HR 8799 planets with VLTI/GRAVITY

Authors: E. Nasedkin, P. Mollière, S. Lacour, M. Nowak, L. Kreidberg, T. Stolker, J. J. Wang, W. O. Balmer, J. Kammerer, J. Shangguan, R. Abuter, A. Amorim, R. Asensio-Torres, M. Benisty, J. -P. Berger, H. Beust, S. Blunt, A. Boccaletti, M. Bonnefoy, H. Bonnet, M. S. Bordoni, G. Bourdarot, W. Brandner, F. Cantalloube, P. Caselli , et al. (73 additional authors not shown)

Abstract: With four companions at separations from 16 to 71 au, HR 8799 is a unique target for direct imaging, presenting an opportunity for the comparative study of exoplanets with a shared formation history. Combining new VLTI/GRAVITY observations obtained within the ExoGRAVITY program with archival data, we perform a systematic atmospheric characterisation of all four planets. We explore different levels… ▽ More With four companions at separations from 16 to 71 au, HR 8799 is a unique target for direct imaging, presenting an opportunity for the comparative study of exoplanets with a shared formation history. Combining new VLTI/GRAVITY observations obtained within the ExoGRAVITY program with archival data, we perform a systematic atmospheric characterisation of all four planets. We explore different levels of model flexibility to understand the temperature structure, chemistry and clouds of each planet using both petitRADTRANS atmospheric retrievals and fits to self-consistent radiative-convective equilibrium models. Using Bayesian Model Averaging to combine multiple retrievals, we find that the HR 8799 planets are highly enriched in metals, with [M/H] $\gtrsim$1, and have stellar to super-stellar C/O ratios. The C/O ratio increases with increasing separation from $0.55^{+0.12}_{-0.10}$ for d to $0.78^{+0.03}_{-0.04}$ for b, with the exception of the innermost planet which has a C/O ratio of $0.87\pm0.03$. By retrieving a quench pressure and using a disequilibrium chemistry model we derive vertical mixing strengths compatible with predictions for high-metallicity, self-luminous atmospheres. Bayesian evidence comparisons strongly favour the presence of HCN in HR 8799 c and e, as well as CH$_{4}$ in HR 8799 c, with detections at $>5σ$ confidence. All of the planets are cloudy, with no evidence for patchiness. The clouds of c, d and e are best fit by silicate clouds lying above a deep iron cloud layer, while the clouds of the cooler HR 8799 b are more likely composed of Na$_{2}$S. With well defined atmospheric properties, future exploration of this system is well positioned to unveil further detail in these planets, extending our understanding of the composition, structure, and formation history of these siblings. △ Less

Submitted 4 April, 2024; originally announced April 2024.

Comments: 45 pages, 25 figures, submitted to Astronomy & Astrophysics

Combining Gaia and GRAVITY: Characterising five new Directly Detected Substellar Companions

Authors: T. O. Winterhalder, S. Lacour, A. Mérand, A. -L. Maire, J. Kammerer, T. Stolker, N. Pourré, C. Babusiaux, R. Abuter, A. Amorim, R. Asensio-Torres, W. O. Balmer, M. Benisty, J. -P. Berger, H. Beust, S. Blunt, A. Boccaletti, M. Bonnefoy, H. Bonnet, M. S. Bordoni, G. Bourdarot, W. Brandner, F. Cantalloube, P. Caselli, B. Charnay , et al. (74 additional authors not shown)

Abstract: Precise mass constraints are vital for the characterisation of brown dwarfs and exoplanets. Here we present how the combination of data obtained by Gaia and GRAVITY can help enlarge the sample of substellar companions with measured dynamical masses. We show how the Non-Single-Star (NSS) two-body orbit catalogue contained in Gaia DR3 can be used to inform high-angular-resolution follow-up observati… ▽ More Precise mass constraints are vital for the characterisation of brown dwarfs and exoplanets. Here we present how the combination of data obtained by Gaia and GRAVITY can help enlarge the sample of substellar companions with measured dynamical masses. We show how the Non-Single-Star (NSS) two-body orbit catalogue contained in Gaia DR3 can be used to inform high-angular-resolution follow-up observations with GRAVITY. Applying the method presented in this work to eight Gaia candidate systems, we detect all eight predicted companions, seven of which were previously unknown and five are of a substellar nature. Among the sample is Gaia DR3 2728129004119806464 B, which - detected at an angular separation of (34.01 $\pm$ 0.15) mas from the host - is the closest substellar companion ever imaged. This translates to a semi-major axis of (0.938 $\pm$ 0.023) AU. WT 766 B, detected at a greater angular separation, was confirmed to be on an orbit exhibiting an even smaller semi-major axis of (0.676 $\pm$ 0.008) AU. The GRAVITY data were then used to break the host-companion mass degeneracy inherent to the Gaia NSS orbit solutions as well as to constrain the orbital solutions of the respective target systems. Knowledge of the companion masses enabled us to further characterise them in terms of their ages, effective temperatures, and radii via the application of evolutionary models. The inferred ages exhibit a distinct bias towards values younger than what is to be expected based on the literature. The results serve as an independent validation of the orbital solutions published in the NSS two-body orbit catalogue and show that the combination of astrometric survey missions and high-angular-resolution direct imaging holds great promise for efficiently increasing the sample of directly imaged companions in the future, especially in the light of Gaia's upcoming DR4 and the advent of GRAVITY+. △ Less

Submitted 24 June, 2024; v1 submitted 19 March, 2024; originally announced March 2024.

Comments: Accepted for publication in Astronomy and Astrophysics

Ultralight vector dark matter search using data from the KAGRA O3GK run

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, A. G. Abac, R. Abbott, H. Abe, I. Abouelfettouh, F. Acernese, K. Ackley, C. Adamcewicz, S. Adhicary, N. Adhikari, R. X. Adhikari, V. K. Adkins, V. B. Adya, C. Affeldt, D. Agarwal, M. Agathos, O. D. Aguiar, I. Aguilar, L. Aiello, A. Ain, P. Ajith, T. Akutsu, S. Albanesi , et al. (1778 additional authors not shown)

Abstract: Among the various candidates for dark matter (DM), ultralight vector DM can be probed by laser interferometric gravitational wave detectors through the measurement of oscillating length changes in the arm cavities. In this context, KAGRA has a unique feature due to differing compositions of its mirrors, enhancing the signal of vector DM in the length change in the auxiliary channels. Here we prese… ▽ More Among the various candidates for dark matter (DM), ultralight vector DM can be probed by laser interferometric gravitational wave detectors through the measurement of oscillating length changes in the arm cavities. In this context, KAGRA has a unique feature due to differing compositions of its mirrors, enhancing the signal of vector DM in the length change in the auxiliary channels. Here we present the result of a search for $U(1)_{B-L}$ gauge boson DM using the KAGRA data from auxiliary length channels during the first joint observation run together with GEO600. By applying our search pipeline, which takes into account the stochastic nature of ultralight DM, upper bounds on the coupling strength between the $U(1)_{B-L}$ gauge boson and ordinary matter are obtained for a range of DM masses. While our constraints are less stringent than those derived from previous experiments, this study demonstrates the applicability of our method to the lower-mass vector DM search, which is made difficult in this measurement by the short observation time compared to the auto-correlation time scale of DM. △ Less

Submitted 5 March, 2024; originally announced March 2024.

Comments: 20 pages, 5 figures

Report number: LIGO-P2300250

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

A catalogue of dual-field interferometric binary calibrators

Authors: M. Nowak, S. Lacour, R. Abuter, A. Amorim, R. Asensio-Torres, W. O. Balmer, M. Benisty, J. -P. Berger, H. Beust, S. Blunt, A. Boccaletti, M. Bonnefoy, H. Bonnet, M. S. Bordoni, G. Bourdarot, W. Brandner, F. Cantalloube, B. Charnay, G. Chauvin, A. Chavez, E. Choquet, V. Christiaens, Y. Clénet, V. Coudé du Foresto, A. Cridland , et al. (75 additional authors not shown)

Abstract: Dual-field interferometric observations with VLTI/GRAVITY sometimes require the use of a "binary calibrator", a binary star whose individual components remain unresolved by the interferometer, with a separation between 400 and 2000 mas for observations with the Units Telescopes (UTs), or 1200 to 3000 mas for the Auxiliary Telescopes (ATs). The separation vector also needs to be predictable to with… ▽ More Dual-field interferometric observations with VLTI/GRAVITY sometimes require the use of a "binary calibrator", a binary star whose individual components remain unresolved by the interferometer, with a separation between 400 and 2000 mas for observations with the Units Telescopes (UTs), or 1200 to 3000 mas for the Auxiliary Telescopes (ATs). The separation vector also needs to be predictable to within 10 mas for proper pointing of the instrument. Up until now, no list of properly vetted calibrators was available for dual-field observations with VLTI/GRAVITY on the UTs. Our objective is to compile such a list, and make it available to the community. We identify a list of candidates from the Washington Double Star (WDS) catalogue, all with appropriate separations and brightness, scattered over the Southern sky. We observe them as part of a dedicated calibration programme, and determine whether these objects are true binaries (excluding higher multiplicities resolved interferometrically but unseen by imaging), and extract measurements of the separation vectors. We combine these new measurements with those available in the WDS to determine updated orbital parameters for all our vetted calibrators. We compile a list of 13 vetted binary calibrators for observations with VLTI/GRAVITY on the UTs, and provide orbital estimates and astrometric predictions for each of them. We show that our list guarantees that there are always at least two binary calibrators at airmass < 2 in the sky over the Paranal observatory, at any point in time. Any Principal Investigator wishing to use the dual-field mode of VLTI/GRAVITY with the UTs can now refer to this list to select an appropriate calibrator. We encourage the use of "whereistheplanet" to predict the astrometry of these calibrators, which seamlessly integrates with "p2Gravity" for VLTI/GRAVITY dual-field observing material preparation. △ Less

Submitted 7 February, 2024; originally announced February 2024.

Comments: 16 pages, 6 figures. Submitted to A&A

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

The GRAVITY young stellar object survey XIII. Tracing the time-variable asymmetric disk structure in the inner AU of the Herbig star HD98922

Authors: GRAVITY Collaboration, V. Ganci, L. Labadie, K. Perraut, A. Wojtczak, J. Kaufhold, M. Benisty, E. Alecian, G. Bourdarot, W. Brandner, A. Caratti o Garatti, C. Dougados, R. Garcia Lopez, J. Sanchez-Bermudez, A. Soulain, A. Amorim, J. -P. Berger, P. Caselli, Y. Clénet, A. Drescher, A. Eckart, F. Eisenhauer, M. Fabricius, H. Feuchtgruber, P. Garcia , et al. (30 additional authors not shown)

Abstract: Temporal variability in the photometric and spectroscopic properties of protoplanetary disks is common in YSO. However, evidence pointing toward changes in their morphology over short timescales has only been found for a few sources, mainly due to a lack of high cadence observations at mas resolution. We combine GRAVITY multi-epoch observations of HD98922 at mas resolution with PIONIER archival da… ▽ More Temporal variability in the photometric and spectroscopic properties of protoplanetary disks is common in YSO. However, evidence pointing toward changes in their morphology over short timescales has only been found for a few sources, mainly due to a lack of high cadence observations at mas resolution. We combine GRAVITY multi-epoch observations of HD98922 at mas resolution with PIONIER archival data covering a total time span of 11 years. We interpret the interferometric visibilities and spectral energy distribution with geometrical models and through radiative transfer techniques. We investigated high-spectral-resolution quantities to obtain information on the properties of the HI BrG-line-emitting region. The observations are best fitted by a model of a crescent-like asymmetric dust feature located at 1 au and accounting for 70% of the NIR emission. The feature has an almost constant magnitude and orbits the central star with a possible sub-Keplerian period of 12 months, although a 9 month period is another, albeit less probable, solution. The radiative transfer models show that the emission originates from a small amount of carbon-rich (25%) silicates, or quantum-heated particles located in a low-density region. Among different possible scenarios, we favor hydrodynamical instabilities in the inner disk that can create a large vortex. The high spectral resolution differential phases in the BrG-line show that the hot-gas component is offset from the star and in some cases is located between the star and the crescent feature. The scale of the emission does not favor magnetospheric accretion as a driving mechanism. The scenario of an asymmetric disk wind or a massive accreting substellar or planetary companion is discussed. With this unique observational data set for HD98922, we reveal morphological variability in the innermost 2 au of its disk region. △ Less

Submitted 31 January, 2024; originally announced January 2024.

Comments: 45 pages, 20 figures, accepted by and to be published in Astronomy & Astrophysics (A&A)

A dynamical measure of the black hole mass in a quasar 11 billion years ago

Authors: R. Abuter, F. Allouche, A. Amorim, C. Bailet, A. Berdeu, J. -P. Berger, P. Berio, A. Bigioli, O. Boebion, M. -L. Bolzer, H. Bonnet, G. Bourdarot, P. Bourget, W. Brandner, Y. Cao, R. Conzelmann, M. Comin, Y. Clénet, B. Courtney-Barrer, R. Davies, D. Defrère, A. Delboulbé, F. Delplancke-Ströbele, R. Dembet, J. Dexter , et al. (102 additional authors not shown)

Abstract: Tight relationships exist in the local universe between the central stellar properties of galaxies and the mass of their supermassive black hole. These suggest galaxies and black holes co-evolve, with the main regulation mechanism being energetic feedback from accretion onto the black hole during its quasar phase. A crucial question is how the relationship between black holes and galaxies evolves… ▽ More Tight relationships exist in the local universe between the central stellar properties of galaxies and the mass of their supermassive black hole. These suggest galaxies and black holes co-evolve, with the main regulation mechanism being energetic feedback from accretion onto the black hole during its quasar phase. A crucial question is how the relationship between black holes and galaxies evolves with time; a key epoch to probe this relationship is at the peaks of star formation and black hole growth 8-12 billion years ago (redshifts 1-3). Here we report a dynamical measurement of the mass of the black hole in a luminous quasar at a redshift of 2, with a look back time of 11 billion years, by spatially resolving the broad line region. We detect a 40 micro-arcsecond (0.31 pc) spatial offset between the red and blue photocenters of the H$α$ line that traces the velocity gradient of a rotating broad line region. The flux and differential phase spectra are well reproduced by a thick, moderately inclined disk of gas clouds within the sphere of influence of a central black hole with a mass of 3.2x10$^{8}$ solar masses. Molecular gas data reveal a dynamical mass for the host galaxy of 6x10$^{11}$ solar masses, which indicates an under-massive black hole accreting at a super-Eddington rate. This suggests a host galaxy that grew faster than the supermassive black hole, indicating a delay between galaxy and black hole formation for some systems. △ Less

Submitted 25 January, 2024; originally announced January 2024.

Comments: 5 pages Main text, 8 figures, 2 tables, to be published in Nature, under embargo until 29 January 2024 16:00 (London)

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

The GRAVITY young stellar object survey XII. The hot gas disk component in Herbig Ae/Be stars

Authors: GRAVITY Collaboration, R. Garcia Lopez, A. Natta, R. Fedriani, A. Caratti o Garatti, J. Sanchez-Bermudez, K. Perraut, C. Dougados, Y. -I. Bouarour, J. Bouvier, W. Brandner, P. Garcia, M. Koutoulaki, L. Labadie, H. Linz, E. Al'ecian, M. Benisty, J. -P. Berger, G. Bourdarot, P. Caselli, Y. Clenet, P. T. de Zeeuw, R. Davies, A. Eckart, F. Eisenhauer , et al. (24 additional authors not shown)

Abstract: The region of protoplanetary disks closest to a star (within 1-2\,au) is shaped by a number of different processes, from accretion of the disk material onto the central star to ejection in the form of winds and jets. Optical and near-IR emission lines are potentially good tracers of inner disk processes if very high spatial and/or spectral resolution are achieved. In this paper, we exploit the cap… ▽ More The region of protoplanetary disks closest to a star (within 1-2\,au) is shaped by a number of different processes, from accretion of the disk material onto the central star to ejection in the form of winds and jets. Optical and near-IR emission lines are potentially good tracers of inner disk processes if very high spatial and/or spectral resolution are achieved. In this paper, we exploit the capabilities of the VLTI-GRAVITY near-IR interferometer to determine the location and kinematics of the hydrogen emission line Bracket gamma. We present VLTI-GRAVITY observations of the Bracket gamma line for a sample of 26 stars of intermediate mass (HAEBE), the largest sample so far analysed with near-IR interferometry. The Bracket gamma line was detected in 17 objects. The emission is very compact (in most cases only marginally resolved), with a size of 10-30R* (1-5 mas). About half of the total flux comes from even smaller regions, which are unresolved in our data. For eight objects, it was possible to determine the position angle (PA) of the line-emitting region, which is generally in agreement with that of the inner-dusty disk emitting the K-band continuum. The position-velocity pattern of the Bracket gamma line-emitting region of the sampled objects is roughly consistent with Keplerian rotation. The exception is HD~45677, which shows more extended emission and more complex kinematics. The most likely scenario for the Bracket gamma origin is that the emission comes from an MHD wind launched very close to the central star, in a region well within the dust sublimation radius. An origin in the bound gas layer at the disk surface cannot be ruled out, while accreting matter provides only a minor fraction of the total flux. These results show the potential of near-IR spectro-interferometry to study line emission in young stellar objects. △ Less

Submitted 15 January, 2024; originally announced January 2024.

Comments: A&A accepted

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

Watching the Universe's acceleration era with the SKAO

Authors: C. M. J. Marques, C. J. A. P. Martins, B. Gilabert López

Abstract: The possibility of watching the Universe expand in real time and in a model-independent way, first envisaged by Allan Sandage more than 60 years ago and known as the redshift drift, is within reach of forthcoming astrophysical facilities, particularly the Extremely Large Telescope (ELT) and the Square Kilometre Array Observatory (SKAO). The latter, probing lower redshifts, enables us to watch the… ▽ More The possibility of watching the Universe expand in real time and in a model-independent way, first envisaged by Allan Sandage more than 60 years ago and known as the redshift drift, is within reach of forthcoming astrophysical facilities, particularly the Extremely Large Telescope (ELT) and the Square Kilometre Array Observatory (SKAO). The latter, probing lower redshifts, enables us to watch the Universe's acceleration era in real time, while the former does the same for the matter era. We use Fisher Matrix Analysis techniques, which we show to give comparable results to those of a Markov Chain Monte Carlo approach, to discuss forecasts for SKAO measurements of the redshift drift and their cosmological impact. We consider specific fiducial cosmological models but mainly rely on a more agnostic cosmographic series (which includes the deceleration and jerk parameters), and we also discuss prospects for measurements of the drift of the drift. Overall, our analysis shows that SKAO measurements, with a reasonable amount of observing time, can provide a competitive probe of the low-redshift accelerating Universe. △ Less

Submitted 15 December, 2023; originally announced December 2023.

Comments: 13 pages, 5 figures, 7 tables; MNRAS (in press)

The GRAVITY young stellar object survey: XI. Imaging the hot gas emission around the Herbig Ae star HD 58647

Authors: Y. -I. Bouarour, R. Garcia Lopez, J. Sanchez-Bermudez, A. Caratti o Garatti, K. Perraut, N. Aimar, A. Amorim, J. -P. Berger, G. Bourdarot, W. Brandner, Y. Clénet, P. T. de Zeeuw, C. Dougados, A. Drescher, A. Eckart, F. Eisenhauer, M. Flock, P. Garcia, E. Gendron, R. Genzel, S. Gillessen, S. Grant, G. Heißel, Th. Henning, L. Jocou , et al. (23 additional authors not shown)

Abstract: We aim to investigate the origin of the HI Br$γ$ emission in young stars by using GRAVITY to image the innermost region of circumstellar disks, where important physical processes such as accretion and winds occur. With high spectral and angular resolution, we focus on studying the continuum and the HI Br$γ$-emitting area of the Herbig star HD58647. Using VLTI-GRAVITY, we conducted observations of… ▽ More We aim to investigate the origin of the HI Br$γ$ emission in young stars by using GRAVITY to image the innermost region of circumstellar disks, where important physical processes such as accretion and winds occur. With high spectral and angular resolution, we focus on studying the continuum and the HI Br$γ$-emitting area of the Herbig star HD58647. Using VLTI-GRAVITY, we conducted observations of HD58647 with both high spectral and high angular resolution. Thanks to the extensive $uv$ coverage, we were able to obtain detailed images of the circumstellar environment at a sub-au scale, specifically capturing the continuum and the Br$γ$-emitting region. Through the analysis of velocity-dispersed images and photocentre shifts, we were able to investigate the kinematics of the HI Br$γ$-emitting region. The recovered continuum images show extended emission where the disk major axis is oriented along a position angle of 14\degr. The size of the continuum emission at 5-sigma levels is $\sim$ 1.5 times more extended than the sizes reported from geometrical fitting (3.69 mas $\pm$ 0.02 mas). This result supports the existence of dust particles close to the stellar surface, screened from the stellar radiation by an optically thick gaseous disk. Moreover, for the first time with GRAVITY, the hot gas component of HD58647 traced by the Br$γ$ ,has been imaged. This allowed us to constrain the size of the Br$γ$-emitting region and study the kinematics of the hot gas; we find its velocity field to be roughly consistent with gas that obeys Keplerian motion. The velocity-dispersed images show that the size of the hot gas emission is from a more compact region than the continuum (2.3 mas $\pm$ 0.2 mas). Finally, the line phases show that the emission is not entirely consistent with Keplerian rotation, hinting at a more complex structure in the hot gaseous disk. △ Less

Submitted 14 December, 2023; originally announced December 2023.

Comments: accepted in A&A on 20/11/2023

VLTI/GRAVITY Provides Evidence the Young, Substellar Companion HD 136164 Ab formed like a "Failed Star"

Authors: William O. Balmer, L. Pueyo, S. Lacour, J. J. Wang, T. Stolker, J. Kammerer, N. Pourré, M. Nowak, E. Rickman, S. Blunt, A. Sivaramakrishnan, D. Sing, K. Wagner, G. -D. Marleau, A. -M. Lagrange, R. Abuter, A. Amorim, R. Asensio-Torres, J. -P. Berger, H. Beust, A. Boccaletti, A. Bohn, M. Bonnefoy, H. Bonnet, M. S. Bordoni , et al. (71 additional authors not shown)

Abstract: Young, low-mass Brown Dwarfs orbiting early-type stars, with low mass ratios ($q\lesssim0.01$), appear intrinsically rare and present a formation dilemma: could a handful of these objects be the highest mass outcomes of ``planetary" formation channels (bottom up within a protoplanetary disk), or are they more representative of the lowest mass ``failed binaries" (formed via disk fragmentation, or c… ▽ More Young, low-mass Brown Dwarfs orbiting early-type stars, with low mass ratios ($q\lesssim0.01$), appear intrinsically rare and present a formation dilemma: could a handful of these objects be the highest mass outcomes of ``planetary" formation channels (bottom up within a protoplanetary disk), or are they more representative of the lowest mass ``failed binaries" (formed via disk fragmentation, or core fragmentation)? Additionally, their orbits can yield model-independent dynamical masses, and when paired with wide wavelength coverage and accurate system age estimates, can constrain evolutionary models in a regime where the models have a wide dispersion depending on initial conditions. We present new interferometric observations of the $16\,\mathrm{Myr}$ substellar companion HD~136164~Ab (HIP~75056~Ab) with VLTI/GRAVITY and an updated orbit fit including proper motion measurements from the Hipparcos-Gaia Catalogue of Accelerations. We estimate a dynamical mass of $35\pm10\,\mathrm{M_J}$ ($q\sim0.02$), making HD~136164~Ab the youngest substellar companion with a dynamical mass estimate. The new mass and newly constrained orbital eccentricity ($e=0.44\pm0.03$) and separation ($22.5\pm1\,\mathrm{au}$) could indicate that the companion formed via the low-mass tail of the Initial Mass Function. Our atmospheric fit to the \texttt{SPHINX} M-dwarf model grid suggests a sub-solar C/O ratio of $0.45$, and $3\times$ solar metallicity, which could indicate formation in the circumstellar disk via disk fragmentation. Either way, the revised mass estimate likely excludes ``bottom-up" formation via core accretion in the circumstellar disk. HD~136164~Ab joins a select group of young substellar objects with dynamical mass estimates; epoch astrometry from future \textit{Gaia} data releases will constrain the dynamical mass of this crucial object further. △ Less

Submitted 13 December, 2023; originally announced December 2023.

Comments: Accepted for publication in The Astronomical Journal. 9 figures, 3 tables

On the rise and fall of galactic ionizing output at the end of reionization

Authors: Christopher Cain, Anson D'Aloisio, Garett Lopez, Nakul Gangolli, Joshua T. Roth

Abstract: Quasar absorption spectra measurements suggest that reionization proceeded rapidly, ended late at $z \sim 5.5$, and was followed by a flat evolution of the ionizing background. Simulations that can reproduce this behavior often rely on a fine-tuned galaxy ionizing emissivity, which peaks at $z \sim 6 - 7$ and drops by a factor of $1.5-2.5$ by $z \sim 5$. This is puzzling since the abundance of gal… ▽ More Quasar absorption spectra measurements suggest that reionization proceeded rapidly, ended late at $z \sim 5.5$, and was followed by a flat evolution of the ionizing background. Simulations that can reproduce this behavior often rely on a fine-tuned galaxy ionizing emissivity, which peaks at $z \sim 6 - 7$ and drops by a factor of $1.5-2.5$ by $z \sim 5$. This is puzzling since the abundance of galaxies has been observed to grow monotonically during this period. Explanations for this include effects such as dust obscuration of ionizing photon escape and feedback due to photo-heating of the IGM. We explore the possibility that this drop in emissivity is instead an artifact of one or more modeling deficiencies in reionization simulations. These include possibly incorrect assumptions about the ionizing spectrum and/or inaccurate modeling of the clumpiness of the IGM. Our results suggest that the need for a drop could be alleviated if simulations are underestimating the IGM opacity from massive, star-forming halos. Other potential modeling issues either have a small effect or require a steeper drop when remedied. We construct an illustrative model in which the emissivity is nearly flat at the end of reionization, evolving only $\sim 0.05$ dex at $5 < z < 7$. More realistic scenarios, however, require a $\sim 0.1-0.3$ dex drop. We also study the evolution of the Ly$α$ effective optical depth distribution in these scenarios and compare them to recent measurements. We find models that feature a hard ionizing spectrum and/or are driven by faint, low-bias sources can most easily reproduce the mean transmission and optical depth distribution of the forest simultaneously. Lastly, we show that the reduced speed of light approximation and low spatial resolution in the forest can lead to erroneous conclusions about the end of reionization. △ Less

Submitted 22 November, 2023; originally announced November 2023.

Comments: 20 pages, 15 figures, to be submitted to MNRAS. Comments welcome

The CARMENES search for exoplanets around M dwarfs. Telluric absorption corrected high S/N optical and near-infrared template spectra of 382 M dwarf stars

Authors: E. Nagel, S. Czesla, A. Kaminski, M. Zechmeister, L. Tal-Or, J. H. M. M. Schmitt, A. Reiners, A. Quirrenbach, A. García López, J. A. Caballero, I. Ribas, P. J. Amado, V. J. S. Béjar, M. Cortés-Contreras, S. Dreizler, A. P. Hatzes, Th. Henning, S. V. Jeffers, M. Kürster, M. Lafarga, M. López-Puertas, D. Montes, J. C. Morales, S. Pedraz, A. Schweitzer

Abstract: Light from celestial objects interacts with the molecules of the Earth's atmosphere, resulting in the production of telluric absorption lines in ground-based spectral data. Correcting for these lines, which strongly affect red and infrared wavelengths, is often needed in a wide variety of scientific applications. Here, we present the template division telluric modeling (TDTM) technique, a method f… ▽ More Light from celestial objects interacts with the molecules of the Earth's atmosphere, resulting in the production of telluric absorption lines in ground-based spectral data. Correcting for these lines, which strongly affect red and infrared wavelengths, is often needed in a wide variety of scientific applications. Here, we present the template division telluric modeling (TDTM) technique, a method for accurately removing telluric absorption lines in stars that exhibit numerous intrinsic features. Based on the Earth's barycentric motion throughout the year, our approach is suited for disentangling telluric and stellar spectral components. By fitting a synthetic transmission model, telluric-free spectra are derived. We demonstrate the performance of the TDTM technique in correcting telluric contamination using a high-resolution optical spectral time series of the feature-rich M3.0 dwarf star Wolf 294 that was obtained with the CARMENES spectrograph. We apply the TDTM approach to the CARMENES survey sample, which consists of 382 targets encompassing 22357 optical and 20314 near-infrared spectra, to correct for telluric absorption. The corrected spectra are coadded to construct template spectra for each of our targets. This library of telluric-free, high signal-to-noise ratio, high-resolution (R>80000) templates comprises the most comprehensive collection of spectral M-dwarf data available to date, both in terms of quantity and quality, and is available at the project website (http://carmenes.cab.inta-csic.es). △ Less

Submitted 23 October, 2023; originally announced October 2023.

Comments: 31 pages, 24 figures, 3 tables, accepted for publication by A&A

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

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)

First VLTI/GRAVITY Observations of HIP 65426 b: Evidence for a Low or Moderate Orbital Eccentricity

Authors: S. Blunt, W. O. Balmer, J. J. Wang, S. Lacour, S. Petrus, G. Bourdarot, J. Kammerer, N. Pourré, E. Rickman, J. Shangguan, T. Winterhalder, R. Abuter, A. Amorim, R. Asensio-Torres, M. Benisty, J. -P. Berger, H. Beust, A. Boccaletti, A. Bohn, M. Bonnefoy, H. Bonnet, W. Brandner, F. Cantalloube, P. Caselli, B. Charnay , et al. (73 additional authors not shown)

Abstract: Giant exoplanets have been directly imaged over orders of magnitude of orbital separations, prompting theoretical and observational investigations of their formation pathways. In this paper, we present new VLTI/GRAVITY astrometric data of HIP 65426 b, a cold, giant exoplanet which is a particular challenge for most formation theories at a projected separation of 92 au from its primary. Leveraging… ▽ More Giant exoplanets have been directly imaged over orders of magnitude of orbital separations, prompting theoretical and observational investigations of their formation pathways. In this paper, we present new VLTI/GRAVITY astrometric data of HIP 65426 b, a cold, giant exoplanet which is a particular challenge for most formation theories at a projected separation of 92 au from its primary. Leveraging GRAVITY's astrometric precision, we present an updated eccentricity posterior that disfavors large eccentricities. The eccentricity posterior is still prior-dependent, and we extensively interpret and discuss the limits of the posterior constraints presented here. We also perform updated spectral comparisons with self-consistent forward-modeled spectra, finding a best fit ExoREM model with solar metallicity and C/O=0.6. An important caveat is that it is difficult to estimate robust errors on these values, which are subject to interpolation errors as well as potentially missing model physics. Taken together, the orbital and atmospheric constraints paint a preliminary picture of formation inconsistent with scattering after disk dispersal. Further work is needed to validate this interpretation. Analysis code used to perform this work is available at https://github.com/sblunt/hip65426. △ Less

Submitted 6 October, 2023; v1 submitted 29 September, 2023; originally announced October 2023.

Comments: 28 pages, 14 figures. Revised and resubmitted to AJ

VLTI/GRAVITY Observations and Characterization of the Brown Dwarf Companion HD 72946 B

Authors: W. O. Balmer, L. Pueyo, T. Stolker, H. Reggiani, S. Lacour, A. -L. Maire, P. Mollière, M. Nowak, D. Sing, N. Pourré, S. Blunt, J. J. Wang, E. Rickman, Th. Henning, K. Ward-Duong, R. Abuter, A. Amorim, R. Asensio-Torres, M. Benisty, J. -P. Berger, H. Beust, A. Boccaletti, A. Bohn, M. Bonnefoy, H. Bonnet , et al. (74 additional authors not shown)

Abstract: Tension remains between the observed and modeled properties of substellar objects, but objects in binary orbits, with known dynamical masses can provide a way forward. HD 72946 B is a recently imaged brown dwarf companion to the nearby, solar type star. We achieve $\sim100~μ\mathrm{as}$ relative astrometry of HD 72946 B in the K-band using VLTI/GRAVITY, unprecedented for a benchmark brown dwarf. W… ▽ More Tension remains between the observed and modeled properties of substellar objects, but objects in binary orbits, with known dynamical masses can provide a way forward. HD 72946 B is a recently imaged brown dwarf companion to the nearby, solar type star. We achieve $\sim100~μ\mathrm{as}$ relative astrometry of HD 72946 B in the K-band using VLTI/GRAVITY, unprecedented for a benchmark brown dwarf. We fit an ensemble of measurements of the orbit using orbitize! and derive a strong dynamical mass constraint $\mathrm{M_B}=69.5\pm0.5~\mathrm{M_{Jup}}$ assuming a strong prior on the host star mass $\mathrm{M_A}=0.97\pm0.01~\mathrm{M_\odot}$ from an updated stellar analysis. We fit the spectrum of the companion to a grid of self-consistent BT-Settl-CIFIST model atmospheres, and perform atmospheric retrievals using petitRADTRANS. A dynamical mass prior only marginally influences the sampled distribution on effective temperature, but has a large influence on the surface gravity and radius, as expected. The dynamical mass alone does not strongly influence retrieved pressure-temperature or cloud parameters within our current retrieval setup. Independent of cloud prescription and prior assumptions, we find agreement within $\pm2\,σ$ between the C/O ratio of the host ($0.52\pm0.05)$ and brown dwarf ($0.43$ to $0.63$), as expected from a molecular cloud collapse formation scenario, but our retrieved metallicities are implausibly high ($0.6-0.8$) in light of an excellent agreement of the data with the solar abundance model grid. Future work on our retrieval framework will seek to resolve this tension. Additional study of low surface-gravity objects is necessary to assess the influence of a dynamical mass prior on atmospheric analysis. △ Less

Submitted 15 September, 2023; v1 submitted 8 September, 2023; originally announced September 2023.

Comments: Accepted for publication in the Astrophysical Journal. 32 pages, 14 figures, 7 tables. v2&3 correct errors in co-author's affiliations, figure rendering, and some grant acknowledgements

Search for Eccentric Black Hole Coalescences during the Third Observing Run of LIGO and Virgo

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, A. G. Abac, R. Abbott, H. Abe, F. Acernese, K. Ackley, C. Adamcewicz, S. Adhicary, N. Adhikari, R. X. Adhikari, V. K. Adkins, V. B. Adya, C. Affeldt, D. Agarwal, M. Agathos, O. D. Aguiar, I. Aguilar, L. Aiello, A. Ain, P. Ajith, T. Akutsu, S. Albanesi, R. A. Alfaidi , et al. (1750 additional authors not shown)

Abstract: Despite the growing number of confident binary black hole coalescences observed through gravitational waves so far, the astrophysical origin of these binaries remains uncertain. Orbital eccentricity is one of the clearest tracers of binary formation channels. Identifying binary eccentricity, however, remains challenging due to the limited availability of gravitational waveforms that include effect… ▽ More Despite the growing number of confident binary black hole coalescences observed through gravitational waves so far, the astrophysical origin of these binaries remains uncertain. Orbital eccentricity is one of the clearest tracers of binary formation channels. Identifying binary eccentricity, however, remains challenging due to the limited availability of gravitational waveforms that include effects of eccentricity. Here, we present observational results for a waveform-independent search sensitive to eccentric black hole coalescences, covering the third observing run (O3) of the LIGO and Virgo detectors. We identified no new high-significance candidates beyond those that were already identified with searches focusing on quasi-circular binaries. We determine the sensitivity of our search to high-mass (total mass $M>70$ $M_\odot$) binaries covering eccentricities up to 0.3 at 15 Hz orbital frequency, and use this to compare model predictions to search results. Assuming all detections are indeed quasi-circular, for our fiducial population model, we place an upper limit for the merger rate density of high-mass binaries with eccentricities $0 < e \leq 0.3$ at $0.33$ Gpc$^{-3}$ yr$^{-1}$ at 90\% confidence level. △ Less

Submitted 7 August, 2023; originally announced August 2023.

Comments: 24 pages, 5 figures

Report number: LIGO-P2300080

Observations of the Crab Nebula and Pulsar with the Large-Sized Telescope Prototype of the Cherenkov Telescope Array

Authors: CTA-LST Project, :, H. Abe, K. Abe, S. Abe, A. Aguasca-Cabot, I. Agudo, N. Alvarez Crespo, L. A. Antonelli, C. Aramo, A. Arbet-Engels, C. Arcaro, M. Artero, K. Asano, P. Aubert, A. Baktash, A. Bamba, A. Baquero Larriva, L. Baroncelli, U. Barres de Almeida, J. A. Barrio, I. Batkovic, J. Baxter, J. Becerra González, E. Bernardini , et al. (467 additional authors not shown)

Abstract: CTA (Cherenkov Telescope Array) is the next generation ground-based observatory for gamma-ray astronomy at very-high energies. The Large-Sized Telescope prototype (LST-1) is located at the Northern site of CTA, on the Canary Island of La Palma. LSTs are designed to provide optimal performance in the lowest part of the energy range covered by CTA, down to $\simeq 20$ GeV. LST-1 started performing a… ▽ More CTA (Cherenkov Telescope Array) is the next generation ground-based observatory for gamma-ray astronomy at very-high energies. The Large-Sized Telescope prototype (LST-1) is located at the Northern site of CTA, on the Canary Island of La Palma. LSTs are designed to provide optimal performance in the lowest part of the energy range covered by CTA, down to $\simeq 20$ GeV. LST-1 started performing astronomical observations in November 2019, during its commissioning phase, and it has been taking data since then. We present the first LST-1 observations of the Crab Nebula, the standard candle of very-high energy gamma-ray astronomy, and use them, together with simulations, to assess the basic performance parameters of the telescope. The data sample consists of around 36 hours of observations at low zenith angles collected between November 2020 and March 2022. LST-1 has reached the expected performance during its commissioning period - only a minor adjustment of the preexisting simulations was needed to match the telescope behavior. The energy threshold at trigger level is estimated to be around 20 GeV, rising to $\simeq 30$ GeV after data analysis. Performance parameters depend strongly on energy, and on the strength of the gamma-ray selection cuts in the analysis: angular resolution ranges from 0.12 to 0.40 degrees, and energy resolution from 15 to 50%. Flux sensitivity is around 1.1% of the Crab Nebula flux above 250 GeV for a 50-h observation (12% for 30 minutes). The spectral energy distribution (in the 0.03 - 30 TeV range) and the light curve obtained for the Crab Nebula agree with previous measurements, considering statistical and systematic uncertainties. A clear periodic signal is also detected from the pulsar at the center of the Nebula. △ Less

Submitted 19 July, 2023; v1 submitted 22 June, 2023; originally announced June 2023.

Comments: Accepted in ApJ. v3: updated author list and acknowledgements, fixed typos and other minor issues

The GRAVITY young stellar object survey -- XI. Probing the inner disk and magnetospheric accretion region of CI Tau

Authors: GRAVITY Collaboration, A. Soulain, K. Perraut, J. Bouvier, G. Pantolmos, A. Caratti o Garatti, P. Caselli, P. Garcia, R. Garcia Lopez

Abstract: Aims: We aim at spatially and spectrally resolving the innermost scale of the young stellar object CI Tau to constrain the inner disk properties and better understand the magnetospheric accretion phenomenon. Methods: The high sensitivity offered by the combination of the four 8-m telescopes of the VLTI allied with the spectral resolution of the K-band beam combiner GRAVITY offers a unique capabili… ▽ More Aims: We aim at spatially and spectrally resolving the innermost scale of the young stellar object CI Tau to constrain the inner disk properties and better understand the magnetospheric accretion phenomenon. Methods: The high sensitivity offered by the combination of the four 8-m telescopes of the VLTI allied with the spectral resolution of the K-band beam combiner GRAVITY offers a unique capability to probe the sub-au scale of the CI Tau system, tracing both dust and gas emission regions. We develop a geometrical model to fit the interferometric observables and constrain the physical properties of the inner dusty disk. The continuum-corrected pure line visibilities have been used to estimate the size of the Br$γ$ emitting region. Results: From the K-band continuum study, we report an highly inclined resolved inner dusty disk, with an inner edge located at a distance of $21\pm2\,R_\star$ from the central star, which is significantly larger than the dust sublimation radius (R$_{sub}= 4.3$ to $8.6\,R_\star$). The inner disk appears misaligned compared to the outer disk observed by ALMA and the non-zero closure phase indicates the presence of a bright asymmetry on the south-west side. From the differential visibilities across the Br$γ$ line, we resolve the line emitting region, and measure a size of $4.8^{+0.8}_{-1.0}$ $R_\star$. Conclusions: The extended inner disk edge compared to the dust sublimation radius is consistent with the claim of an inner planet, CI Tau b, orbiting close-in. The inner-outer disk misalignment may be induced by gravitational torques or magnetic warping. The size of the Br$γ$ emitting region is consistent with the magnetospheric accretion process. Assuming it corresponds to the magnetospheric radius, it is significantly smaller than the co-rotation radius, which suggests an unstable accretion regime that is consistent with CI Tau being a burster. △ Less

Submitted 14 May, 2023; originally announced May 2023.

Comments: 12 pages, 13 figures

Journal ref: A&A 674, A203 (2023)

Study of the bipolar jet of the YSO Th 28 with VLT/SINFONI: Jet morphology and H$_2$ emission

Authors: S. Yu. Melnikov, P. A. Boley, N. S. Nikonova, A. Caratti o Garatti, R. Garcia Lopez, B. Stecklum, J. Eislöffel, G. Weigelt

Abstract: $Context.$ The YSO Th 28 possesses a highly collimated jet, which clearly exhibits an asymmetric brightness of its jet lobes at optical and NIR wavelengths. There may be asymmetry in the jet plasma parameters in opposite jet lobes (e.g. electron density, temperature, and outflow velocity). $Aims.… ▽ More $Context.$ The YSO Th 28 possesses a highly collimated jet, which clearly exhibits an asymmetric brightness of its jet lobes at optical and NIR wavelengths. There may be asymmetry in the jet plasma parameters in opposite jet lobes (e.g. electron density, temperature, and outflow velocity). $Aims.$ We examined the Th 28 jet in a 3"x3" where the jet material is collimated and accelerated. Our goal is to map the morphology and determine its physical parameters to determine the physical origin of such asymmetries. $Methods.$ We present $JHK$-spectra of Th 28 obtained with the SINFONI on the (VLT, ESO) in June-July 2015. $Results.$ The [Fe II] emission originates in collimated jet lobes. Two new axial knots are detected at 1" in the blue lobe and 1".2 in the red lobe. The H$_2$ radiation is emitted from an extended region with a radius of $\gtrsim270$ au, which is perpendicular to the jet. The PV diagrams of the bright H$_2$ lines reveal faint H$_2$ emission along both jet lobes as well. The compact and faint H I emission (Pa$β$ and Br$γ$) comes from two regions, namely from a spherical region around the star and from the jet lobes. The size of the jet launching region is derived as 0".015 ($\sim$3 au at 185 pc), and the initial opening angle of the Th 28 jet is $\sim28^0$, which makes this jet substantially less collimated than most jets from other CTTs. $Conclusions.$ The emission in [Fe II], H$_2$, and H I lines suggests a morphology in which the ionised gas in the disc appears to be disrupted by the jet. The resolved disc-like H$_2$ emission most likely arises in the disc atmosphere from shocks caused by a radial uncollimated wind. The asymmetry of the [Fe II] photocentre shifts with respect to the jet source arises in the immediate vicinity of the driving source of Th28 and suggests that the observed brightness asymmetry is intrinsic as well. △ Less

Submitted 26 April, 2023; v1 submitted 25 April, 2023; originally announced April 2023.

Comments: 18 pages, 14 figures

Journal ref: A&A 673, A156 (2023)

Forbidden emission lines in protostellar outflows and jets with MUSE

Authors: Lizxandra Flores-Rivera, Mario Flock, Nicolás Kurtovic, Bernd Husemann, Andrea Banzatti, Simon C. Ringqvist, Sebastian Kamann, André Müller, Christian Fendt, Rebeca Garcia Lopez, Gabriel-Dominique Marleau, Thomas Henning, Carlos Carrasco-Gonzalez, Roy van Boekel, Miriam Keppler, Ralf Launhardt, Yuhiko Aoyama

Abstract: Forbidden emission lines in protoplanetary disks are a key diagnostic in studies of the evolution of the disk and the host star. We report spatially resolved emission lines, [OI] 6300, 6363, [NII] 6548, 6583, H$α$, and [SII] 6716, 6730 Angstrom that are believed to be associated with jets and magnetically driven winds in the inner disks. Observations were carried out with the optical integral fiel… ▽ More Forbidden emission lines in protoplanetary disks are a key diagnostic in studies of the evolution of the disk and the host star. We report spatially resolved emission lines, [OI] 6300, 6363, [NII] 6548, 6583, H$α$, and [SII] 6716, 6730 Angstrom that are believed to be associated with jets and magnetically driven winds in the inner disks. Observations were carried out with the optical integral field spectrograph of the Multi Unit Spectroscopic Explorer (MUSE), at the Very Large Telescope (VLT). With a resolution of 0.025 X 0.025 arcsec$^{2}$, we aim to derive the position angle of the outflow/jet (PA$_{outflow/jet}$) that is connected with the inner disk. The forbidden emission lines analyzed here have their origin at the inner parts of the protoplanetary disk. From the maximum intensity emission along the outflow/jet in DL Tau, CI Tau, DS Tau, IP Tau, and IM Lup, we were able to reliably measure the PA$_{outflow/jet}$ for most of the identified lines. We found that our estimates agree with PA$_{dust}$ for most of the disks. These estimates depend on the signal-to-noise level and the collimation of the outflow (jet). The outflows/jets in CIDA 9, GO Tau, and GW Lup are too compact for a PA$_{outflow/jet}$ to be estimated. Based on our kinematics analysis, we confirm that DL Tau and CI Tau host a strong outflow/jet with line-of-sight velocities much greater than 100 km s$^{-1}$, whereas DS Tau, IP Tau, and IM Lup velocities are lower and their structures encompass low-velocity components to be more associated with winds. Our estimates for the mass-loss rate, $\dot{M}_{loss}$, range between (1.1-6.5)X10$^{-7}$-10$^{-8}$ $M_{\odot}$ yr$^{-1}$ for the disk-outflow/jet systems analyzed here. The outflow/jet systems analyzed here are aligned within around 1 degree between the inner and outer disk. Further observations are needed to confirm a potential misalignment in IM Lup. △ Less

Submitted 6 January, 2023; originally announced January 2023.

Comments: Accepted on December 19 2022

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)

The GRAVITY Young Stellar Object survey -- IX. Spatially resolved kinematics of hot hydrogen gas in the star/disk interaction region of T Tauri stars

Authors: GRAVITY Collaboration, J. A. Wojtczak, L. Labadie, K. Perraut, B. Tessore, A. Soulain, V. Ganci, J. Bouvier, C. Dougados, E. Alécian, H. Nowacki, G. Cozzo, W. Brandner, A. Caratti o Garatti, P. Garcia, R. Garcia Lopez, J. Sanchez-Bermudez, A. Amorim, M. Benisty, J. -P. Berger, G. Bourdarot, P. Caselli, Y. Clénet, P. T. de Zeeuw, R. Davies , et al. (36 additional authors not shown)

Abstract: Aims: We aim to spatially and spectrally resolve the Br-gamma hydrogen emission line with the methods of interferometry in order to examine the kinematics of the hydrogen gas emission region in the inner accretion disk of a sample of solar-like young stellar objects. The goal is to identify trends and categories among the sources of our sample and to discuss whether or not they can be tied to diff… ▽ More Aims: We aim to spatially and spectrally resolve the Br-gamma hydrogen emission line with the methods of interferometry in order to examine the kinematics of the hydrogen gas emission region in the inner accretion disk of a sample of solar-like young stellar objects. The goal is to identify trends and categories among the sources of our sample and to discuss whether or not they can be tied to different origin mechanisms associated with Br-gamma emission in T Tauri stars, chiefly and most prominently magnetospheric accretion. Methods: We observed a sample of seven T Tauri stars for the first time with VLTI GRAVITY, recording spectra and spectrally dispersed interferometric quantities across the Br-gamma line in the NIR K-band. We use them to extract the size of the Br-gamma emission region and the photocenter shifts. To assist in the interpretation, we also make use of radiative transfer models of magnetospheric accretion to establish a baseline of expected interferometric signatures if accretion is the primary driver of Br-gamma emission. Results: From among our sample, we find that five of the seven T~Tauri stars show an emission region with a half-flux radius in the range broadly expected for magnetospheric truncation. Two of the five objects also show Br-gamma emission primarily originating from within the corotation radius, while two other objects exhibit extended emission on a scale beyond 10 R$_*$, one of them even beyond the K~band continuum half-flux radius of 11.3 R$_*$. Conclusions: We find strong evidence to suggest that for the two weakest accretors in the sample, magnetospheric accretion is the primary driver of Br-gamma radiation. The results for the remaining sources imply either partial or strong contributions coming from spatially extended emission components in the form of outflows, such as stellar or disk winds. △ Less

Submitted 23 November, 2022; v1 submitted 24 October, 2022; originally announced October 2022.

Comments: Version corresponds to the one in A&A production. The author list has been amended with "GRAVITY Collaboration" as first entry in keeping with the convention established for papers published by the GRAVITY consortium in recent years

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

Multi-wavelength study of the galactic PeVatron candidate LHAASO J2108+5157

Authors: S. Abe, A. Aguasca-Cabot, I. Agudo, N. Alvarez Crespo, L. A. Antonelli, C. Aramo, A. Arbet-Engels, M. Artero, K. Asano, P. Aubert, A. Baktash, A. Bamba, A. Baquero Larriva, L. Baroncelli, U. Barres de Almeida, J. A. Barrio, I. Batkovic, J. Baxter, J. Becerra González, E. Bernardini, M. I. Bernardos, J. Bernete Medrano, A. Berti, P. Bhattacharjee, N. Biederbeck , et al. (245 additional authors not shown)

Abstract: LHAASO J2108+5157 is one of the few known unidentified Ultra-High-Energy (UHE) gamma-ray sources with no Very-High-Energy (VHE) counterpart, recently discovered by the LHAASO collaboration. We observed LHAASO J2108+5157 in the X-ray band with XMM-Newton in 2021 for a total of 3.8 hours and at TeV energies with the Large-Sized Telescope prototype (LST-1), yielding 49 hours of good quality data. In… ▽ More LHAASO J2108+5157 is one of the few known unidentified Ultra-High-Energy (UHE) gamma-ray sources with no Very-High-Energy (VHE) counterpart, recently discovered by the LHAASO collaboration. We observed LHAASO J2108+5157 in the X-ray band with XMM-Newton in 2021 for a total of 3.8 hours and at TeV energies with the Large-Sized Telescope prototype (LST-1), yielding 49 hours of good quality data. In addition, we analyzed 12 years of Fermi-LAT data, to better constrain emission of its High-Energy (HE) counterpart 4FGL J2108.0+5155. We found an excess (3.7 sigma) in the LST-1 data at energies E > 3 TeV. Further analysis in the whole LST-1 energy range assuming a point-like source, resulted in a hint (2.2 sigma) of hard emission which can be described with a single power law with photon index Gamma = 1.6 +- 0.2 between 0.3 - 100 TeV. We did not find any significant extended emission which could be related to a Supernova Remnant (SNR) or Pulsar Wind Nebula (PWN) in the XMM-Newton data, which puts strong constraints on possible synchrotron emission of relativistic electrons. The LST-1 and LHAASO observations can be explained as inverse Compton-dominated leptonic emission of relativistic electrons with a cutoff energy of $100^{+70}_{-30}$ TeV. The low magnetic field in the source imposed by the X-ray upper limits on synchrotron emission is compatible with a hypothesis of a PWN or a TeV halo. The lack of a pulsar in the neighborhood of the UHE source is a challenge to the PWN/TeV-halo scenario. The UHE gamma rays can also be explained as $π^0$ decay-dominated hadronic emission due to interaction of relativistic protons with one of the two known molecular clouds in the direction of the source. The hard spectrum in the LST-1 band is compatible with protons escaping a shock around a middle-aged SNR because of their high low-energy cut-off. △ Less

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

Comments: 26 pages, 9 figures, 5 tables; Astronomy and Astrophysics (accepted on 6 March 2023)

Journal ref: A&A 673, A75 (2023)

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

Direct discovery of the inner exoplanet in the HD206893 system. Evidence for deuterium burning in a planetary-mass companion

Authors: S. Hinkley, S. Lacour, G. -D. Marleau, A. M. Lagrange, J. J. Wang, J. Kammerer, A. Cumming, M. Nowak, L. Rodet, T. Stolker, W. -O. Balmer, S. Ray, M. Bonnefoy, P. Mollière, C. Lazzoni, G. Kennedy, C. Mordasini, R. Abuter, S. Aigrain, A. Amorim, R. Asensio-Torres, C. Babusiaux, M. Benisty, J. -P. Berger, H. Beust , et al. (89 additional authors not shown)

Abstract: Long term precise radial velocity (RV) monitoring of the nearby star HD206893, as well as anomalies in the system proper motion, have suggested the presence of an additional, inner companion in the system. Here we describe the results of a multi-epoch search for the companion responsible for this RV drift and proper motion anomaly using the VLTI/GRAVITY instrument. Utilizing information from ongoi… ▽ More Long term precise radial velocity (RV) monitoring of the nearby star HD206893, as well as anomalies in the system proper motion, have suggested the presence of an additional, inner companion in the system. Here we describe the results of a multi-epoch search for the companion responsible for this RV drift and proper motion anomaly using the VLTI/GRAVITY instrument. Utilizing information from ongoing precision RV measurements with the HARPS spectrograph, as well as Gaia host star astrometry, we report a high significance detection of the companion HD206893c over three epochs, with clear evidence for Keplerian orbital motion. Our astrometry with $\sim$50-100 $μ$arcsec precision afforded by GRAVITY allows us to derive a dynamical mass of 12.7$^{+1.2}_{-1.0}$ M$_{\rm Jup}$ and an orbital separation of 3.53$^{+0.08}_{-0.06}$ au for HD206893c. Our fits to the orbits of both companions in the system utilize both Gaia astrometry and RVs to also provide a precise dynamical estimate of the previously uncertain mass of the B component, and therefore derive an age of $155\pm15$ Myr. We find that theoretical atmospheric/evolutionary models incorporating deuterium burning for HD206893c, parameterized by cloudy atmospheres provide a good simultaneous fit to the luminosity of both HD206893B and c. In addition to utilizing long-term RV information, this effort is an early example of a direct imaging discovery of a bona fide exoplanet that was guided in part with Gaia astrometry. Utilizing Gaia astrometry is expected to be one of the primary techniques going forward to identify and characterize additional directly imaged planets. Lastly, this discovery is another example of the power of optical interferometry to directly detect and characterize extrasolar planets where they form at ice-line orbital separations of 2-4\,au. △ Less

Submitted 3 April, 2023; v1 submitted 9 August, 2022; originally announced August 2022.

Comments: Accepted to A&A

A detailed analysis of the Gl 486 planetary system

Authors: J. A. Caballero, E. Gonzalez-Alvarez, M. Brady, T. Trifonov, T. G. Ellis, C. Dorn, C. Cifuentes, K. Molaverdikhani, J. L. Bean, T. Boyajian, E. Rodriguez, J. Sanz-Forcada, M. R. Zapatero Osorio, C. Abia, P. J. Amado, N. Anugu, V. J. S. Bejar, C. L. Davies, S. Dreizler, F. Dubois, J. Ennis, N. Espinoza, C. D. Farrington, A. Garcia Lopez, T. Gardner , et al. (42 additional authors not shown)

Abstract: The Gl 486 system consists of a very nearby, relatively bright, weakly active M3.5 V star at just 8 pc with a warm transiting rocky planet of about 1.3 R_Terra and 3.0 M_Terra that is ideal for both transmission and emission spectroscopy and for testing interior models of telluric planets. To prepare for future studies, we collected light curves of seven new transits observed with the CHEOPS space… ▽ More The Gl 486 system consists of a very nearby, relatively bright, weakly active M3.5 V star at just 8 pc with a warm transiting rocky planet of about 1.3 R_Terra and 3.0 M_Terra that is ideal for both transmission and emission spectroscopy and for testing interior models of telluric planets. To prepare for future studies, we collected light curves of seven new transits observed with the CHEOPS space mission and new radial velocities obtained with MAROON-X/Gemini North and CARMENES/Calar Alto telescopes, together with previously published spectroscopic and photometric data from the two spectrographs and TESS. We also performed interferometric observations with the CHARA Array and new photometric monitoring with a suite of smaller telescopes. From interferometry, we measure a limb-darkened disc angular size of the star Gl 486. Together with a corrected Gaia EDR3 parallax, we obtain a stellar radius. We also measure a stellar rotation period at P_rot ~ 49.9 d, an upper limit to its XUV (5-920 AA) flux with new Hubble/STIS data, and, for the first time, a variety of element abundances (Fe, Mg, Si, V, Sr, Zr, Rb) and C/O ratio. Besides, we impose restrictive constraints on the presence of additional components, either stellar or substellar, in the system. With the input stellar parameters and the radial-velocity and transit data, we determine the radius and mass of the planet Gl 486 b at R_p = 1.343+/0.063 R_Terra and M_p = 3.00+/-0.13 M_Terra. From the planet parameters and the stellar element abundances, we infer the most probable models of planet internal structure and composition, which are consistent with a relatively small metallic core with respect to the Earth, a deep silicate mantle, and a thin volatile upper layer. With all these ingredients, we outline prospects for Gl 486 b atmospheric studies, especially with forthcoming James Webb Space Telescope observations (abridged). △ Less

Submitted 20 June, 2022; originally announced June 2022.

Comments: A&A, in press. See https://carmenes.caha.es/ext/pressreleases/GJ486/Exoearth.Video.mp4

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

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

The mass of Beta Pictoris c from Beta Pictoris b orbital motion

Authors: S. Lacour, J. J. Wang, L. Rodet, M. Nowak, J. Shangguan, H. Beust, A. -M. Lagrange, R. Abuter, A. Amorim, R. Asensio-Torres, M. Benisty, J. -P. Berger, S. Blunt, A. Boccaletti, A. Bohn, M. -L. Bolzer, M. Bonnefoy, H. Bonnet, G. Bourdarot, W. Brandner, F. Cantalloube, P. Caselli, B. Charnay, G. Chauvin, E. Choquet , et al. (74 additional authors not shown)

Abstract: We aim to demonstrate that the presence and mass of an exoplanet can now be effectively derived from the astrometry of another exoplanet. We combined previous astrometry of $β$ Pictoris b with a new set of observations from the GRAVITY interferometer. The orbital motion of $β$ Pictoris b is fit using Markov chain Monte Carlo simulations in Jacobi coordinates. The inner planet, $β$ Pictoris c, was… ▽ More We aim to demonstrate that the presence and mass of an exoplanet can now be effectively derived from the astrometry of another exoplanet. We combined previous astrometry of $β$ Pictoris b with a new set of observations from the GRAVITY interferometer. The orbital motion of $β$ Pictoris b is fit using Markov chain Monte Carlo simulations in Jacobi coordinates. The inner planet, $β$ Pictoris c, was also reobserved at a separation of 96\,mas, confirming the previous orbital estimations. From the astrometry of planet b only, we can (i) detect the presence of $β$ Pictoris c and (ii) constrain its mass to $10.04^{+4.53}_{-3.10}\,M_{\rm Jup}$. If one adds the astrometry of $β$ Pictoris c, the mass is narrowed down to $9.15^{+1.08}_{-1.06}\,M_{\rm Jup}$. The inclusion of radial velocity measurements does not affect the orbital parameters significantly, but it does slightly decrease the mass estimate to $8.89^{+0.75}_{-0.75}\,M_{\rm Jup}$. With a semimajor axis of $2.68\pm0.02$\,au, a period of $1221\pm15$ days, and an eccentricity of $0.32\pm0.02$, the orbital parameters of $β$ Pictoris c are now constrained as precisely as those of $β$ Pictoris b. The orbital configuration is compatible with a high-order mean-motion resonance (7:1). The impact of the resonance on the planets' dynamics would then be negligible with respect to the secular perturbations, which might have played an important role in the eccentricity excitation of the outer planet. △ Less

Submitted 22 September, 2021; originally announced September 2021.

Comments: Accepted in A&A Letter

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

The GRAVITY Young Stellar Object survey VIII. Gas and dust faint inner rings in the hybrid disk of HD141569

Authors: GRAVITY Collaboration, V. Ganci, L. Labadie, L. Klarmann, A. de Valon, K. Perraut, M. Benisty, W. Brandner, A. Caratti o Garatti, C. Dougados, F. Eupen, R. Garcia Lopez, R. Grellmann, J. Sanchez-Bermudez, A. Wojtczak, P. Garcia, A. Amorim, M. Bauböck, J. -P. Berger, P. Caselli, Y. Clénet, V. Coudé du Foresto, P. T. de Zeeuw, A. Drescher, G. Duvert , et al. (38 additional authors not shown)

Abstract: The formation and evolution of planetary systems impact the primordial accretion disk. HD141569 is the only known pre-main sequence star characterized by a hybrid disk. Observations probed the outer-disk structure showing a complex system of rings and interferometric observations attempted to characterize its inner 5 au region, but derived limited constraints. The goal of this work was to explore… ▽ More The formation and evolution of planetary systems impact the primordial accretion disk. HD141569 is the only known pre-main sequence star characterized by a hybrid disk. Observations probed the outer-disk structure showing a complex system of rings and interferometric observations attempted to characterize its inner 5 au region, but derived limited constraints. The goal of this work was to explore with new high-resolution interferometric observations the properties of the dust and gas in the internal regions of HD141569. We observed HD141569 on mas scales with GRAVITY/VLTI in the near-infrared at low and high spectral resolution. We interpreted the visibilities and spectral energy distribution with geometrical models and radiative transfer techniques to constrain the dust emission. We analyzed the high spectral resolution quantities to investigate the properties of the Br-Gamma line emitting region. Thanks to the combination of three different epochs, GRAVITY resolves the inner dusty disk in the K band. Data modeling shows that an IR excess of about 6% is spatially resolved and that the origin of this emission is confined in a ring of material located at a radius of 1 au from the star with a width smaller than 0.3 au. The MCMax modeling suggests that this emission could originate from a small amount of QHPs, while large silicate grain models cannot reproduce at the same time the observational constraints on the properties of near-IR and mid-IR fluxes. The differential phases in the Br-Gamma line clearly show an S-shape that can be best reproduced witha gas disk in Keplerian rotation, confined within 0.09 au. This is also hinted at by the double-peaked Br-Gamma emission line shape. The modeling of the continuum and gas emission shows that the inclination and position angle of these two components are consistent with a system showing relatively coplanar rings on all scales. △ Less

Submitted 22 September, 2021; v1 submitted 21 September, 2021; originally announced September 2021.

Comments: Accepted for publication in A&A; 25 pages, 15 figures, 5 tables

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

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