Search

OT 081 is a well-known, luminous blazar that is remarkably variable in many energy bands. We present the first broadband study of the source which includes very-high-energy (VHE, $E>$100\,GeV) $\gamma$-ray data taken by the MAGIC and H.E.S.S. imaging Cherenkov telescopes. The discovery of VHE $\gamma$-ray emission happened during a high state of $\gamma$-ray activity in July 2016, observed by many instruments from radio to VHE $\gamma$-rays. We identify four states of activity of the source, one of which includes VHE $\gamma$-ray emission. Variability in the VHE domain is found on daily timescales. The intrinsic VHE spectrum can be described by a power-law with index $3.27\pm0.44_{\rm stat}\pm0.15_{\rm sys}$ (MAGIC) and $3.39\pm0.58_{\rm stat}\pm0.64_{\rm sys}$ (H.E.S.S.) in the energy range of 55--300\,GeV and 120--500\,GeV, respectively. The broadband emission cannot be sucessfully reproduced by a simple one-zone synchrotron self-Compton model. Instead, an additional external Compton component is required. We test a lepto-hadronic model that reproduces the dataset well and a proton-synchrotron dominated model that requires an extreme proton luminosity. Emission models that are able to successfully represent the data place the emitting region well outside of the Broad Line Region (BLR) to a location at which the radiative environment is dominated by the infrared thermal radiation field of the dusty torus. In the scenario described by this flaring activity, the source appears to be an FSRQ, in contrast with past categorizations. This suggests that the source can be considered to be a transitional blazar, intermediate between BL~Lac and FSRQ objects., Comment: Accepted on MNRAS Corresponding authors: M. Manganaro, J. Becerra Gonz\'alez, M. Seglar-Arroyo, D. A. Sanchez

We present optical and IR observations from maximum light until around 600 d of SN 2022jli, a peculiar SE SN showing two maxima, each one with a peak luminosity of about 3 x 10^{42} erg/s and separated by 50 d. The second maximum is followed by periodic undulations with a period of P ~ 12.5 days. The spectra and the photometric evolution of the first maximum are consistent with the behaviour of a standard SE SN with an ejecta mass of 1.5 +/- 0.4 Msun, and a nickel mass of 0.12 +/- 0.01 Msun. The optical spectra after 400 d correspond to a standard SN Ic event, and at late times SN 2022jli exhibits a significant drop in the optical luminosity implying that the physical phenomena that produced the secondary maximum has ceased to power the SN light curve. One possibility is that the second maximum is powered by a magnetar with an initial spin period of P=48.5 ms and a magnetic field of B = 8.5x10^{14} G, while the light curve periodic undulations could be produced by accretion of material from a companion star onto the neutron star in a binary system. The near-IR spectra shows clear 1st CO overtone emission from about 190 d after the first maximum, and it becomes undetected at 400 d. A significant near-IR excess from hot dust emission is detected at 238 d produced by either newly formed dust in the SN ejecta or due to a strong near-IR dust echo. Depending on the assumptions of the dust composition, the estimated dust mass is 2-16 x 10^{-4} Msun. The magnetar power of the second maximum can fit in a more general picture where magnetars are the power source of super-luminous SNe, could produce their frequent bumps and undulations, and where pulsars could produce the late time excess observed in some SE SNe. The detection of CO and the potential detection of dust formed in the ejecta of SN2022jli are important to understand the formation molecules and dust in the ejecta of SE SNe., Comment: Submitted to A&A (27 pages and 22 figures)

Magnetic fields are considered to be key components of massive stars, with a far-reaching impact on their evolution and ultimate fate. A magnetic mechanism was suggested for the collimated explosion of massive stars, relevant for long-duration gamma-ray bursts, X-ray flashes, and asymmetric core collapse supernovae. However, the origin of the observed stable, globally organized magnetic fields in massive stars is still a matter of debate: it has been argued that they can be fossil, dynamo generated, or generated by strong binary interactions or merging events. Taking into account that multiplicity is a fundamental characteristic of massive stars, observational evidence is accumulating that the magnetism originates through interaction between the system components, both during the initial mass transfer or when the stellar cores merge., Comment: 6 pages, 4 figures, to be published in the International Astronomical Union Proceedings Series for the IAU GA 2024 (Astronomy in Focus, Focus Meeting 8: Advances and Challenges in Understanding the Solar and Stellar Dynamos)

Trans-Neptunian objects (TNOs) are considered to be among the most primitive objects in our Solar System. Knowledge of their primary physical properties is essential for understanding their origin and the evolution of the outer Solar System. We predicted a stellar occultation by this TNO for 2020 October 23 UT and ran a specific campaign to investigate this event. We derived the projected profile shape and size from the occultation observations by means of an elliptical fit to the occultation chords. We also performed photometric observations of (143707) 2003 UY117 to obtain the absolute magnitude and the rotational period from the observed rotational light curve. Finally, we combined these results to derive the three-dimensional shape, volume-equivalent diameter, and geometric albedo for this TNO. From the stellar occultation, we obtained a projected ellipse with axes of $(282 \pm 18) \times (184 \pm 32)$ km. The area-equivalent diameter for this ellipse is $D_\textrm{eq,A} = 228 \pm 21$ km. From our photometric $R$ band observations, we derived an absolute magnitude of $H_V = 5.97 \pm 0.07$ mag using $V-R = 0.46 \pm 0.07$ mag, which was derived from a $V$ band subset of these data. The rotational light curve has a peak-to-valley amplitude of $\Delta m = 0.36 \pm 0.13$ mag. We find the most likely rotation period to be $P = 12.376 \pm 0.0033$ hours. By combining the occultation with the rotational light curve results and assuming a triaxial ellipsoid, we derived axes of $a \times b \times c = (332 \pm 24)$ km $\times$ $(216 \pm 24)$ km $\times$ $(180\substack{+28\\-24})$ km for this ellipsoid, and therefore a volume-equivalent diameter of $D_\textrm{eq,V} = 235 \pm 25$ km. Finally, the values for the absolute magnitude and for the area-equivalent diameter yield a geometric albedo of $p_V = 0.139 \pm 0.027$., Comment: 9 pages, 7 figures, accepted for publication in Astronomy and Astrophysics on Sept 13, 2024

Over the past few years alone, the lensing community has discovered thousands of strong lens candidates, and spectroscopically confirmed hundreds of them. In this time of abundance, it becomes pragmatic to focus our time and resources on the few extraordinary systems, in order to most efficiently study the universe. In this paper, we present such a system: DESI-090.9854-35.9683, a cluster-scale lens at $z_{\rm l} = 0.49$, with seven observed lensed sources around the core, and additional lensed sources further out in the cluster. From the number and the textbook configuration of the lensed images, a tight constraint on the mass potential of the lens is possible. This would allow for detailed analysis on the dark and luminous matter content within galaxy clusters, as well as a probe into dark energy and high-redshift galaxies. We present our spatially resolved kinematic measurements of this system from the Very Large Telescope Multi Unit Spectroscopic Explorer, which confirm five of these source galaxies (in ascending order, at $z_{\rm s} = 0.962, 0.962, 1.166, 1.432,$ and $1.432$). With previous Hubble Space Telescope imaging in the F140W and F200LP bands, we also present a simple two power-law profile flux-based lens model that, for a cluster lens, well models the five lensed arc families with redshifts. We determine the mass to be $M(< \theta_{\rm E}) = 4.78\times10^{13} M_{\odot}$ for the primary mass potential. From the model, we extrapolate the redshift of one of the two source galaxies not yet spectroscopically confirmed to be at $z_{\rm s}=4.52^{+1.03}_{-0.71}$., Comment: 14 pages, 6 figures

The efficiency of science observation Short-Term Scheduling (STS) can be defined as being a function of how many highly ranked observations are completed per unit time. Current STS at ESO's Paranal observatory is achieved through filtering and ranking observations via well-defined algorithms, leading to a proposed observation at time t. This Paranal STS model has been successfully employed for more than a decade. Here, we summarise the current VLT(I) STS model, and outline ongoing efforts of optimising the scientific return of both the VLT(I) and future ELT. We describe the STS simulator we have built that enables us to evaluate how changes in model assumptions affect STS effectiveness. Such changes include: using short-term predictions of atmospheric parameters instead of assuming their constant time evolution; assessing how the ranking weights on different observation parameters can be changed to optimise the scheduling; changing STS to be more `dynamic' to consider medium-term scheduling constraints. We present specific results comparing how machine learning predictions of the seeing can improve STS efficiency when compared to the current model of using the last 10\,min median of the measured seeing for observation selection., Comment: 14 pages, 8 figures, SPIE conference proceedings 13096

A peculiar subtype of Type Ia supernovae (SNe), 03fg-like (super-Chandrasekhar) SNe, show different observational properties from prototypical Type Ia SNe, typically having high luminosity at the light-curve peak, low expansion velocities, and strong carbon features. The origin of this class of Type Ia SNe has been actively debated. Recent nebular-phase infrared observations of the 03fg-like Type Ia SN 2022pul using the James Webb Space Telescope revealed large-scale asymmetries in the ejecta and the presence of the strong [Ne II] line at 12.81 $\mu$m, suggesting a violent merger of two white dwarfs as its origin. Polarimetry is another powerful tool to study overall ejecta asymmetries of spatially-unresolved SNe. Here, we aim to check the universality of the violent merger scenario as the origin of the 03fg-like Type Ia SNe, by studying their explosion geometries using polarimetry. In this letter, we present imaging-polarimetric observations of the two 03fg-like Type Ia SNe 2021zny and 2022ilv. SNe 2021zny and 2022ilv show high intrinsic polarization ($\sim1$ % -$\sim2$ %), which might be composed of multiple components with different polarization angles. This indicates that they have complex aspherical structures in their ejecta, supporting the violent merger scenario for their origin. Our observations provide the first clear evidence from polarimetry for such aspherical structures., Comment: 6 pages, 4 figures, submitted to A&A Letters

Recent observations of core-collapse supernovae (SNe) suggest aspherical explosions. Globally aspherical structures in SN explosions are regarded as the key for understanding their explosion mechanism. However, the exact explosion geometries from the inner cores to the outer envelopes are poorly understood. Here, we present photometric, spectroscopic and polarimetric observations of the Type IIP SN 2021yja and discuss its explosion geometry, in comparison to those of other Type IIP SNe that show large-scale aspherical structures in their hydrogen envelopes (SNe 2012aw, 2013ej and 2017gmr). During the plateau phase, SNe 2012aw and 2021yja exhibit high continuum polarization characterized by two components with perpendicular polarization angles. This behavior can be interpreted to be due to a bipolar explosion, composed of a polar (energetic) and an equatorial (bulk) components of the SN ejecta. In such a bipolar explosion, an aspherical axis created by the polar ejecta would be dominating at early phases, while the perpendicular axis along the equatorial ejecta would emerge at late phases after the receding of the photosphere in the polar ejecta. The interpretation of the bipolar explosions in SNe 2012aw and 2021yja is also supported by other observational properties, including the time evolution of the line velocities and the line shapes in the nebular spectra. The polarization of other Type IIP SNe that show large-scale aspherical structures in the hydrogen envelope (SNe 2013ej and 2017gmr) is also consistent with the bipolar-explosion scenario, although this is not conclusive., Comment: 16 pages, 10 figures, accepted for publication in A&A letter

Studies of the magnetic characteristics of massive stars have recently received significant attention because they are progenitors of highly magnetised compact objects. Stars initially more massive than about 8M_sun leave behind neutron stars and black holes by the end of their evolution. The merging of binary compact remnant systems produces astrophysical transients detectable by gravitational wave observatories. Studies of magnetic fields in massive stars with low metallicities are of particular interest because they provide important information on the role of magnetic fields in the star formation of the early Universe. While several detections of massive Galactic magnetic stars have been reported in the last few decades, the impact of a low-metallicity environment on the occurrence and strength of stellar magnetic fields has not yet been explored. Because of the similarity between Of?p stars in the Magellanic Clouds (MCs) and Galactic magnetic Of?p stars, which possess globally organised magnetic fields, we searched for magnetic fields in Of?p stars in the MCs. Additionally, we observed the massive contact binary Cl* NGC 346 SSN7 in the Small Magellanic Cloud to test the theoretical scenario that the origin of magnetic fields involves a merger event or a common envelope evolution. We obtained and analysed measurements of the magnetic field in four massive Of?p stars in the MCs and the binary Cl* NGC 346 SSN7 using the ESO/VLT FORS2 spectrograph in spectropolarimetric mode. We detected kilogauss-scale magnetic fields in two Of?p-type stars and in the contact binary Cl* NGC 346 SSN7. These results suggest that the impact of low metallicity on the occurrence and strength of magnetic fields in massive stars is low. However, because the explored stellar sample is very small, additional observations of massive stars in the MCs are necessary., Comment: 7 pages, 2 tables, 8 figures, published as a letter in A&A

We have learnt in the last decades that the majority of galaxies belong to high density regions interconnected in a sponge-like fashion. This large-scale structure is characterised by clusters, filaments, walls, where most galaxies concentrate, but also under-dense regions, called voids. The void regions and the galaxies within represent an ideal place for the study of galaxy formation and evolution as they are largely unaffected by the complex physical processes that transform galaxies in high-density environments. These void galaxies can hold the key as well to answer current challenges to the $\Lambda$CDM paradigm. The Calar Alto Void Integral-field Treasury surveY (CAVITY) is a Legacy project approved by the Calar Alto Observatory to obtain spatially resolved spectroscopic information of $\sim300$ void galaxies in the Local Universe (0.005 < z < 0.050) covering from -17.0 to -21.5 in $\rm r$ band absolute magnitude. It officially started in January 2021 and has been awarded 110 useful dark observing nights at the 3.5 m telescope using the PMAS spectrograph. Complementary follow-up projects including deep optical imaging, integrated, as well as resolved CO data, and integrated HI spectra, have joint the PMAS observations and naturally complete the scientific aim of characterising galaxies in cosmic voids. The extension data has been denominated CAVITY+. The data will be available to the whole community in different data releases, the first of which is planned for July 2024, and it will provide the community with PMAS data cubes for around 100 void galaxies through a user friendly, and well documented, database platform. We present here the survey, sample selection, data reduction, quality control schemes, science goals, and some examples of the scientific power of the CAVITY and CAVITY+ data., Comment: 19 pages, 11 figures. Accepted for publication in A&A - Replacement after A&A minor language edition

We present a JWST MIRI/MRS spectrum (5-27 $\mathrm{\mu}$m) of the Type Ia supernova (SN Ia), SN 2021aefx at $+415$ days past $B$-band maximum. The spectrum, which was obtained during the iron-dominated nebular phase, has been analyzed in combination with previous JWST observations of SN 2021aefx, to provide the first JWST time series analysis of an SN Ia. We find the temporal evolution of the [Co III] 11.888 $\mathrm{\mu}$m feature directly traces the decay of $^{56}$Co. The spectra, line profiles, and their evolution are analyzed with off-center delayed-detonation models. Best fits were obtained with White Dwarf (WD) central densities of $\rho_c=0.9-1.1\times 10^9$g cm$^{-3}$, a WD mass of M$_{\mathrm{WD}}$=1.33-1.35M$_\odot$, a WD magnetic field of $\approx10^6$G, and an off-center deflagration-to-detonation transition at $\approx$ 0.5 $M_\odot$ seen opposite to the line of sight of the observer (-30). The inner electron capture core is dominated by energy deposition from $\gamma$-rays whereas a broader region is dominated by positron deposition, placing SN 2021aefx at +415 d in the transitional phase of the evolution to the positron-dominated regime. The formerly `flat-tilted' profile at 9 $\mathrm{\mu}$m now has significant contribution from [Ni IV], [Fe II], and [Fe III] and less from [Ar III], which alters the shape of the feature as positrons excite mostly the low-velocity Ar. Overall, the strength of the stable Ni features in the spectrum is dominated by positron transport rather than the Ni mass. Based on multi-dimensional models, our analysis is consistent with a single-spot, close-to-central ignition with an indication for a pre-existing turbulent velocity field, and excludes a multiple-spot, off-center ignition., Comment: Accepted for publication in ApJ

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., Comment: Accepted for publication in Monthly Notices of the Royal Astronomical Society. 19 pages, 9 figures. Corresponding authors: Elisa Prandini, Antonio Stamerra, Talvikki Hovatta

We present JWST spectral and photometric observations of the Type IIP supernova (SN) 2022acko at ~50 days past explosion. These data are the first JWST spectral observations of a core-collapse SN. We identify ~30 different H I features, other features associated with products produced from the CNO cycle, and s-process elements such as Sc II and Ba II. By combining the JWST spectra with ground-based optical and NIR spectra, we construct a full Spectral Energy Distribution from 0.4 to 25 microns and find that the JWST spectra are fully consistent with the simultaneous JWST photometry. The data lack signatures of CO formation and we estimate a limit on the CO mass of < 10^{-8} solar mass. We demonstrate how the CO fundamental band limits can be used to probe underlying physics during stellar evolution, explosion, and the environment. The observations indicate little mixing between the H envelope and C/O core in the ejecta and show no evidence of dust. The data presented here set a critical baseline for future JWST observations, where possible molecular and dust formation may be seen.

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\gamma} = 3 \times 10^{-12}$ GeV$^{-1}$., Comment: 25 pages, 10 figures, accepted for publication in Physics of the Dark Universe

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\sigma$ 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., Comment: Accepted for publication in Astronomy & Astrophysics. 29 pages, 22 figures. Corresponding authors: Axel Arbet Engels, Felix Schmuckermaier, David Paneque

This work aims at constraining the size, shape, and geometric albedo of the dwarf planet candidate 2002 MS4 through the analysis of nine stellar occultation events. Using multichord detection, we also studied the object's topography by analyzing the obtained limb and the residuals between observed chords and the best-fitted ellipse. We predicted and organized the observational campaigns of nine stellar occultations by 2002 MS4 between 2019 and 2022, resulting in two single-chord events, four double-chord detections, and three events with three to up to sixty-one positive chords. Using 13 selected chords from the 8 August 2020 event, we determined the global elliptical limb of 2002 MS4. The best-fitted ellipse, combined with the object's rotational information from the literature, constrains the object's size, shape, and albedo. Additionally, we developed a new method to characterize topography features on the object's limb. The global limb has a semi-major axis of 412 $\pm$ 10 km, a semi-minor axis of 385 $\pm$ 17 km, and the position angle of the minor axis is 121 $^\circ$ $\pm$ 16$^\circ$. From this instantaneous limb, we obtained 2002 MS4's geometric albedo and the projected area-equivalent diameter. Significant deviations from the fitted ellipse in the northernmost limb are detected from multiple sites highlighting three distinct topographic features: one 11 km depth depression followed by a 25$^{+4}_{-5}$ km height elevation next to a crater-like depression with an extension of 322 $\pm$ 39 km and 45.1 $\pm$ 1.5 km deep. Our results present an object that is $\approx$138 km smaller in diameter than derived from thermal data, possibly indicating the presence of a so-far unknown satellite. However, within the error bars, the geometric albedo in the V-band agrees with the results published in the literature, even with the radiometric-derived albedo.

Gravitational lensing provides unique insights into astrophysics and cosmology, including the determination of galaxy mass profiles and constraining cosmological parameters. We present spectroscopic confirmation and lens modeling of the strong lensing system DESI-253.2534+26.8843, discovered in the Dark Energy Spectroscopic Instrument (DESI) Legacy Imaging Surveys data. This system consists of a massive elliptical galaxy surrounded by four blue images forming an Einstein Cross pattern. We obtained spectroscopic observations of this system using the Multi Unit Spectroscopic Explorer (MUSE) on ESO's Very Large Telescope (VLT) and confirmed its lensing nature. The main lens, which is the elliptical galaxy, has a redshift of $z_{L1} = 0.636\pm 0.001$, while the spectra of the background source images are typical of a starburst galaxy and have a redshift of $z_s = 2.597 \pm 0.001$. Additionally, we identified a faint galaxy foreground of one of the lensed images, with a redshift of $z_{L2} = 0.386$. We employed the GIGA-Lens modeling code to characterize this system and determined the Einstein radius of the main lens to be $\theta_{E} =2.520{''}_{-0.031}^{+0.032}$, which corresponds to a velocity dispersion of $\sigma$ = 379 $\pm$ 2 km s$^{-1}$. Our study contributes to a growing catalog of this rare kind of strong lensing systems and demonstrates the effectiveness of spectroscopic integral field unit observations and advanced modeling techniques in understanding the properties of these systems., Comment: Accepted for publication in ApJL

We present a JWST mid-infrared spectrum of the under-luminous Type Ia Supernova (SN Ia) 2022xkq, obtained with the medium-resolution spectrometer on the Mid-Infrared Instrument (MIRI) $\sim130$ days post-explosion. We identify the first MIR lines beyond 14 $\mu$m in SN Ia observations. We find features unique to under-luminous SNe Ia, including: isolated emission of stable Ni, strong blends of [Ti II], and large ratios of singly ionized to doubly ionized species in both [Ar] and [Co]. Comparisons to normal-luminosity SNe Ia spectra at similar phases show a tentative trend between the width of the [Co III] 11.888 $\mu$m feature and the SN light curve shape. Using non-LTE-multi-dimensional radiation hydro simulations and the observed electron capture elements we constrain the mass of the exploding white dwarf. The best-fitting model shows that SN 2022xkq is consistent with an off-center delayed-detonation explosion of a near-Chandrasekhar mass WD (M$_{\rm ej}$ $\approx 1.37$ M$_{\odot}$) of high-central density ($\rho_c \geq 2.0\times10^{9}$ g cm$^{-3}$) seen equator on, which produced M($^{56}$Ni) $= 0.324$ M$_{\odot}$ and M($^{58}$Ni) $\geq 0.06$ M$_{\odot}$. The observed line widths are consistent with the overall abundance distribution; and the narrow stable Ni lines indicate little to no mixing in the central regions, favoring central ignition of sub-sonic carbon burning followed by an off-center DDT beginning at a single point. Additional observations may further constrain the physics revealing the presence of additional species including Cr and Mn. Our work demonstrates the power of using the full coverage of MIRI in combination with detailed modeling to elucidate the physics of SNe Ia at a level not previously possible., Comment: 31 pages, 18 figures, accepted to ApJ; updated to accepted version

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., Comment: 13 pages, 6 figures, 2 tables. Accepted for publication in Monthly Notices of the Royal Astronomical Society

The BL Lac 1ES 2344+514 is known for temporary extreme properties (e.g., a shift of the synchrotron SED peak energy $\nu_{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 $\nu_{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 $\nu_{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., Comment: Accepted for publication in Astronomy & Astrophysics

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., Comment: Accepted for publication in Astronomy & Astrophysics

Accurate distance determination to astrophysical objects is essential for the understanding of their intrinsic brightness and size. The distance to SN 1987A has been previously measured by the expanding photosphere method, and by using the angular size of the circumstellar rings with absolute sizes derived from light curves of narrow UV emission lines, with reported distances ranging from 46.77 kpc to 55 kpc. In this study, we independently determined the distance to SN 1987A using photometry and imaging polarimetry observations of AT 2019xis, a light echo of SN 1987A, by adopting a radiative transfer model of the light echo developed in Ding et al. (2021). We obtained distances to SN 1987A in the range from 49.09 $\pm$ 2.16 kpc to 59.39 $\pm$ 3.27 kpc, depending on the interstellar polarization and extinction corrections, which are consistent with the literature values. This study demonstrates the potential of using light echoes as a tool for distance determination to astrophysical objects in the Milky Way, up to kiloparsec level scales., Comment: Accepted for publication in ApJL

The explosion processes of supernovae (SNe) are imprinted in their explosion geometries. Here, we study the intrinsic polarization of 15 hydrogen-rich core-collapse SNe and explore the relation with the photometric and spectroscopic properties. Our sample shows diverse properties of the continuum polarization. The polarization of most SNe has a low degree at early phases but shows a sudden rise to $\sim 1$ \% degree at certain points during the photospheric phase as well as a slow decline during the tail phase, with a constant polarization angle. The variation in the timing of peak polarisation values implies diversity in the explosion geometry: some SNe have aspherical structures only in their helium cores, while in other SNe these reach out to a significant part of the outer hydrogen envelope with a common axis from the helium core to the hydrogen envelope. Other SNe show high polarization from early phases and a change of the polarization angle around the middle of the photospheric phase. This implies that the ejecta are significantly aspherical to the outermost layer and have multi-directional aspherical structures. Exceptionally, the Type~IIL SN~2017ahn shows low polarization at both the photospheric and tail phases. Our results show that the timing of the polarization rise in Type~IIP SNe is likely correlated with their brightness, velocity and the amount of radioactive Ni produced: brighter SNe with faster ejecta velocity and a larger $^{56}$Ni mass have more extended-aspherical explosion geometries. In particular, there is a clear correlation between the timing of the polarization rise and the explosion energy, that is, the explosion asphericity is proportional to the explosion energy. This implies that the development of a global aspherical structure, e.g., a jet, might be the key to realising an energetic SN in the mechanism of SN explosions., Comment: 27 pages, 30 figures

Aims. A transiting planet candidate with a sub-Neptune radius orbiting the nearby ($d$ = 51.9$\pm$0.07 pc) M1.5 V star TOI-1470 with a period of $\sim$2.5 d was announced by the NASA Transiting Exoplanet Survey Satellite (TESS), which observed the field of TOI-1470 in four different sectors. We aim to validate its planetary nature using precise radial velocities (RVs) taken with the CARMENES spectrograph. Methods. We obtained 44 RV measurements with CARMENES spanning eight months between 3 June 2020 and 17 January 2021. For a better characterization of the parent star activity, we also collected contemporaneous optical photometric observations at the Joan Or\'o and Sierra Nevada Observatories, and we retrieved archival photometry from the literature. We used ground-based photometric observations from MuSCAT and also from MuSCAT2 and MuSCAT3 to confirm the planetary transit signals. We performed a combined photometric and spectroscopic analysis by including Gaussian processes and Keplerian orbits to simultaneously account for the stellar activity and planetary signals. Results. We estimate that TOI-1470 has a rotation period of 29$\pm$3 d based on photometric and spectroscopic data. The combined analysis confirms the discovery of the announced transiting planet, TOI-1470 b, with an orbital period of 2.527093$\pm$0.000003 d, a mass of $7.32^{+1.21}_{-1.24}$ M$_{\oplus}$, and a radius of $2.18^{+0.04}_{-0.04}$ R$_{\oplus}$. We also discover a second transiting planet that was not announced previously by TESS, TOI-1470 c, with an orbital period of 18.08816$\pm$0.00006 d, a mass of $7.24^{+2.87}_{-2.77}$ M$_{\oplus}$, and a radius of $2.47^{+0.02}_{-0.02}$ R$_{\oplus}$. The two planets are placed on the same side of the radius valley of M dwarfs and lie between TOI-1470 and the inner border of its habitable zone., Comment: arXiv admin note: substantial text overlap with arXiv:2111.14602

Studies of magnetic fields in the most evolved massive stars, the Wolf-Rayet stars, are of special importance because they are progenitors of certain types of supernovae. The first detection of a magnetic field of the order of a few hundred Gauss in the WN7 star WR55, based on a few FORS2 low-resolution spectropolarimetric observations, was reported in 2020. In this work we present new FORS2 observations allowing us to detect magnetic and spectroscopic variability with a period of 11.90 h. No significant frequencies were detected in TESS and ASAS-SN photometric observations. Importantly, magnetic field detections are achieved currently only in two Wolf-Rayet stars, WR6 and WR55, both showing the presence of corotating interacting regions., Comment: Accepted to be published in MNRAS Letters

Abridged - Stars with ZAMS masses between 140 and $260 M_\odot$ are thought to explode as pair-instability supernovae (PISNe). During their thermonuclear runaway, PISNe can produce up to several tens of solar masses of radioactive nickel, resulting in luminous transients similar to some superluminous supernovae (SLSNe). Yet, no unambiguous PISN has been discovered so far. SN2018ibb is a H-poor SLSN at $z=0.166$ that evolves extremely slowly compared to the hundreds of known SLSNe. Between mid 2018 and early 2022, we monitored its photometric and spectroscopic evolution from the UV to the NIR with 2-10m class telescopes. SN2018ibb radiated $>3\times10^{51} \rm erg$ during its evolution, and its bolometric light curve reached $>2\times10^{44} \rm erg\,s^{-1}$ at peak. The long-lasting rise of $>93$ rest-frame days implies a long diffusion time, which requires a very high total ejected mass. The PISN mechanism naturally provides both the energy source ($^{56}$Ni) and the long diffusion time. Theoretical models of PISNe make clear predictions for their photometric and spectroscopic properties. SN2018ibb complies with most tests on the light curves, nebular spectra and host galaxy, potentially all tests with the interpretation we propose. Both the light curve and the spectra require 25-44 $M_\odot$ of freshly nucleosynthesised $^{56}$Ni, pointing to the explosion of a metal-poor star with a He-core mass of 120-130 $M_\odot$ at the time of death. This interpretation is also supported by the tentative detection of [Co II]$\lambda$1.025$\mu$m, which has never been observed in any other PISN candidate or SLSN before. Powering by a central engine, such as a magnetar or a black hole, can be excluded with high confidence. This makes SN2018ibb by far the best candidate for being a PISN, to date., Comment: Accepted in A&A, the revised version includes a PISN rate estimate and an additional test with PISN models. 47 pages, main text 41 pages, 38 figures, 16 Tables

We present linear polarimetry for seven hydrogen-poor superluminous supernovae (SLSNe-I). For SN 2017gci, for which we present two epochs of spectropolarimetry at +3 d and +29 d post-peak in rest frame, accompanied by four epochs of imaging polarimetry up to +108 d. The spectropolarimetry at +3 d shows increasing polarisation degree P towards the redder wavelengths and exhibits signs of axial symmetry, but at +29 d P=0 throughout the spectrum implying that the photosphere of SN 2017gci evolved from a slightly aspherical configuration to a more spherical one in the first month post-peak. However, an increase of P to 0.5% at +55 d accompanied by a different orientation of the axial symmetry compared to +3 d implies the presence of additional sources of polarisation at this phase. The increase in polarisation is possibly caused by interaction with circumstellar matter as already suggested by a knee in the light curve and a possible detection of broad Ha emission. We also analysed the sample of all 16 SLSNe-I with polarimetry to date. The data taken during the early spectroscopic phase show consistently low P indicating spherical photospheres. No clear relation between the polarimetry and spectral phase was seen when the spectra resemble Type Ic SNe during the photospheric and nebular phases. The light curve decline rate also shows no clear relation with the polarisation properties. While only slow-evolving SLSNe-I have shown non-zero P, the fast-evolving ones have not been observed at sufficiently late times to conclude that none of them exhibit changing P. However, the four SLSNe-I with increasing polarisation degree also have irregular light curve declines. For up to half of them, the photometric, spectroscopic and polarimetric properties are affected by CSM interaction. As such CSM interaction clearly plays an important role in understanding the polarimetric evolution of SLSNe-I., Comment: Accepted to A&A on 18/03/2023

Optical spectropolarimetry of the normal thermonuclear supernova SN2019np from -14.5 to +14.5 days relative to B-band maximum detected an intrinsic continuum polarization, p(cont), of 0.21+-0.09% at the first epoch. Between days -11.5 to +05, p(cont) remained about 0 and by day +14.5 was again significant at 0.19+-0.10%. Not considering the first epoch, the dominant axis of SiII(6355A) was roughly constant, staying close to the continuum until both rotated in opposite directions on day +14.5. Detailed radiation-hydrodynamical simulations produce a very steep density slope in the outermost ejecta so that the low first-epoch p(const) of about 0.2% nevertheless suggests a separate structure with an axis ratio of about 2 in the outer carbon-rich 3...5E-3 Mo. Large-amplitude fluctuations in the polarization profiles and a flocculent appearance of the polar diagram for the CaII near-infrared triplet (NIR3) may be related by a common origin. The temporal evolution of the polarization spectra agrees with an off-center delayed detonation. The late-time increase in polarization and the possible change in position angle are also consistent with an aspherical 56Ni core. The p(cont) and the absorptions due to Si II(6355A) and the CaII NIR3 form in the same region of the extended photosphere, with an interplay between line occultation and thermalization producing p. Small-scale polarization features may be due to small-scale structures, but many could be related to atomic patterns of the quasi-continuum; the hardly have an equivalent in the total-flux spectrum. We compare SN2019np to other SNe and develop future objectives and strategies for SNIa spectropolarimetry., Comment: 23 pages, 18 figures, 2 tables, MNRAS, submitted 11/9/22, revised 12/9/22

We present a JWST/MIRI low-resolution mid-infrared (MIR) spectroscopic observation of the normal Type Ia supernova (SN Ia) SN 2021aefx at +323 days past rest-frame B-band maximum light. The spectrum ranges from 4-14 um, and shows many unique qualities including a flat-topped [Ar III] 8.991 um profile, a strongly tilted [Co III] 11.888 um feature, and multiple stable Ni lines. These features provide critical information about the physics of the explosion. The observations are compared to synthetic spectra from detailed NLTE multi-dimensional models. The results of the best-fitting model are used to identify the components of the spectral blends and provide a quantitative comparison to the explosion physics. Emission line profiles and the presence of electron capture (EC) elements are used to constrain the mass of the exploding white dwarf (WD) and the chemical asymmetries in the ejecta. We show that the observations of SN 2021aefx are consistent with an off-center delayed-detonation explosion of a near-Chandrasekhar mass (Mch) WD at a viewing angle of -30 degrees relative to the point of the deflagration-to-detonation transition. From the strength of the stable Ni lines we determine that there is little to no mixing in the central regions of the ejecta. Based on both the presence of stable Ni and the Ar velocity distributions, we obtain a strict lower limit of 1.2 Msun of the initial WD, implying that most sub-Mch explosions models are not viable models for SN 2021aefx. The analysis here shows the crucial importance of MIR spectra for distinguishing between explosion scenarios for SNe Ia., Comment: 21 pages, 9 figures, 4 tables, accepted to ApJL; updated to accepted version

The introduction of deep wide-field surveys in recent years and the adoption of machine learning techniques have led to the discoveries of $\mathcal{O}(10^4)$ strong gravitational lensing systems and candidates. However, the discovery of multiply lensed transients remains a rarity. Lensed transients and especially lensed supernovae are invaluable tools to cosmology as they allow us to constrain cosmological parameters via lens modeling and the measurements of their time delays. In this paper, we develop a pipeline to perform a targeted lensed transient search. We apply this pipeline to 5807 strong lenses and candidates, identified in the literature, in the DESI Legacy Imaging Surveys Data Release 9 (DR9) footprint. For each system, we analyze every exposure in all observed bands (DECam $g$, $r$, and $z$). Our pipeline finds, groups, and ranks detections that are in sufficient proximity temporally and spatially. After the first round of inspection, for promising candidate systems, we further examine the newly available DR10 data (with additional $i$ and $\textrm{Y}$ bands). Here we present our targeted lensed supernova search pipeline and seven new lensed supernova candidates, including a very likely lensed supernova $-$ probably a Type Ia $-$ in a system with an Einstein radius of $\sim 1.5''$., Comment: 53 pages, 50 figures, 3 tables

Tidal disruption events (TDEs) occur when a star orbiting a massive black hole is sufficiently close to be tidally ripped apart by the black hole. AT 2022cmc is the first relativistic TDE that was observed (and discovered) as an optically bright and fast transient, showing signatures of non-thermal radiation induced by a jet which is oriented towards the Earth. In this work, we present optical linear and circular polarization measurements, observed with VLT/FORS2 in the $R$-band (which corresponds to the blue/UV part of the spectrum in rest frame), $\sim$ 7.2 and $\sim$ 12.2 rest-frame days after the first detection, respectively, when the light curve of the transient had settled in a bright blue plateau. Both linear and circular polarization are consistent with zero, $p_{lin}$ = 0.14 $\pm$ 0.73 % and $p_{cir}$ = $-$0.30 $\pm$ 0.53 %. This is the highest S/N linear polarization measurement obtained for a relativistic TDE and the first circular polarimetry for such a transient. The non detection of the linear and circular polarization is consistent with the scenario of AT 2022cmc being a TDE where the thermal component (disk+outflows) is viewed pole-on, assuming an axially symmetric geometry. The presence and effect of a jet and/or external shocks are, however, difficult to disentangle., Comment: Accepted for publication in ApJL

Line-like features in TeV $\gamma$-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 $\gamma$-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 $\gamma$-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 \sigma v \rangle \lesssim 5 \times 10^{-28}\,\mathrm{cm^3\,s^{-1}}$ at 1 TeV and $\langle \sigma 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., Comment: Accepted in Phys. Rev. Lett., includes Supplemental materials. 9+4 pages, 3+3 figures, 3 tables. Corresponding authors: T. Inada, D. Kerszberg, M. H\"utten

The SNR G106.3+2.7, detected at 1--100 TeV energies by different $\gamma$-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 $\gamma$-ray emission is leptonic or hadronic. With the better angular resolution provided by these new MAGIC data compared to earlier $\gamma$-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 $\gamma$-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 $\gamma$-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$\gamma$ 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., Comment: 13 pages, 7 figures, Accepted for publication in A&A

The BL Lac object 1ES 0647+250 is one of the few distant $\gamma$-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 $\gamma$ 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., Comment: 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)

Investigating interstellar (IS) dust properties in external galaxies is important not only to infer the intrinsic properties of astronomical objects but also to understand the star/planet formation in the galaxies. From the non-Milky-Way-like extinction and interstellar polarization (ISP) observed in reddened Type Ia supernovae (SNe), it has been suggested that their host galaxies contain dust grains whose properties are substantially different from the Milky-Way (MW) dust. It is important to investigate the universality of such non-MW-like dust in the universe. Here we report spectropolarimetry of two highly-extinguished Type II SNe (SN 2022aau and SN 2022ame). SN 2022aau shows a polarization maximum at a shorter wavelength than MW stars, which is also observed in some Type Ia SNe. This is clear evidence for the existence of non-MW-like dust in its host galaxy (i.e., NGC 1672). This fact implies that such non-MW-like dust might be more common in some environments than expected, and thus it might affect the picture of the star/planet formation. On the other hand, SN 2022ame shows MW-like ISP, implying the presence of MW-like dust in its host galaxy (i.e., NGC 1255). Our findings confirm that dust properties of galaxies are diverse, either locally or globally. The present work demonstrates that further investigation of IS dust properties in external galaxies using polarimetry of highly-reddened SNe is promising, providing a great opportunity to study the universality of such non-MW-like dust grains in the universe., Comment: 11 pages, 5 figures, accepted for publication in ApJ Letters

A rare class of supernovae (SNe) is characterized by strong interaction between the ejecta and several solar masses of circumstellar matter (CSM) as evidenced by strong Balmer-line emission. Within the first few weeks after the explosion, they may display spectral features similar to overluminous Type Ia SNe, while at later phase their observation properties exhibit remarkable similarities with some extreme case of Type IIn SNe that show strong Balmer lines years after the explosion. We present polarimetric observations of SN2018evt obtained by the ESO Very Large Telescope from 172 to 219 days after the estimated time of peak luminosity to study the geometry of the CSM. The nonzero continuum polarization decreases over time, suggesting that the mass loss of the progenitor star is aspherical. The prominent H$\alpha$ emission can be decomposed into a broad, time-evolving component and an intermediate-width, static component. The former shows polarized signals, and it is likely to arise from a cold dense shell (CDS) within the region between the forward and reverse shocks. The latter is significantly unpolarized, and it is likely to arise from shocked, fragmented gas clouds in the H-rich CSM. We infer that SN2018evt exploded inside a massive and aspherical circumstellar cloud. The symmetry axes of the CSM and the SN appear to be similar. SN\,2018evt shows observational properties common to events that display strong interaction between the ejecta and CSM, implying that they share similar circumstellar configurations. Our preliminary estimate also suggests that the circumstellar environment of SN2018evt has been significantly enriched at a rate of $\sim0.1$ M$_\odot$ yr$^{-1}$ over a period of $>100$ yr., Comment: 29 pages, 20 figures, 10 tables, accepted for publication in MNRAS

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., Comment: 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)

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) $\gamma$-rays. A significant correlation ($>$3$\sigma$) between X-rays and VHE $\gamma$-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 $\gamma$-rays and radio, with the radio lagging by more than 100 days, placing the $\gamma$-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 $\gamma$-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., Comment: 55 pages, 30 figures, 14 tables, accepted by APJS. Corresponding authors are L. Heckmann, D. Paneque, S. Gasparyan, M. Cerruti, and N. Sahakyan

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., Comment: 15 pages, 5 figures. Accepted for publication in MNRAS

We report spectropolarimetric observations of the Type Ia supernova (SN) 2021rhu at four epochs: $-$7, +0, +36, and +79 days relative to its $B$-band maximum luminosity. A wavelength-dependent continuum polarization peaking at $3890 \pm 93$ Angstroms and reaching a level of $p_{\rm max}=1.78% \pm 0.02$% was found. The peak of the polarization curve is bluer than is typical in the Milky Way, indicating a larger proportion of small dust grains along the sightline to the SN. After removing the interstellar polarization, we found a pronounced increase of the polarization in the CaII near-infrared triplet, from $\sim$0.3% at day $-$7 to $\sim$2.5% at day +79. No temporal evolution in high-resolution flux spectra across the NaID and CaIIH&K features was seen from days +39 to +74, indicating that the late-time increase in polarization is intrinsic to the SN as opposed to being caused by scattering of SN photons in circumstellar or interstellar matter. We suggest that an explanation for the late-time rise of the CaII near-infrared triplet polarization may be the alignment of calcium atoms in a weak magnetic field through optical excitation/pumping by anisotropic radiation from the SN., Comment: 14 pages, 9 figures, accepted for publication in ApJ

The time delay between multiple images of strongly lensed quasars is a powerful tool for measuring the Hubble constant (H0). To achieve H0 measurements with higher precision and accuracy using the time delay, it is crucial to expand the sample of lensed quasars. We conduct a search for strongly lensed quasars in the Dark Energy Spectroscopic Instrument (DESI) Legacy Imaging Surveys. The DESI Legacy Surveys comprise 19,000 deg2 of the extragalactic sky observed in three optical bands (g, r, and z), making it well suited for the discovery of new strongly lensed quasars. We apply an autocorrelation algorithm to ~5 million objects classified as quasars in the DESI Quasar Sample. These systems are visually inspected and ranked. Here, we present 436 new multiply lensed and binary quasar candidates, 65 of which have redshifts from Sloan Digital Sky Survey Data Release 16. We provide redshifts for an additional 18 candidates from the SuperNova Integral Field Spectrograph., Comment: 16 pages, 11 figures, 2 tables

A star crossing the tidal radius of a supermassive black hole will be spectacularly ripped apart with an accompanying burst of radiation. A few tens of such tidal disruption events (TDEs) have now been identified in the optical wavelengths, but the exact origin of the strong optical emission remains inconclusive. Here we report polarimetric observations of three TDEs. The continuum polarization is independent of wavelength, while emission lines are partially depolarized. These signatures are consistent with optical photons being scattered and polarized in an envelope of free electrons. An almost axisymmetric photosphere viewed from different angles is in broad agreement with the data, but there is also evidence for deviations from axial symmetry before the peak of the flare and significant time evolution at early times, compatible with the rapid formation of an accretion disk. By combining a super-Eddington accretion model with a radiative transfer code we generate predictions for the degree of polarization as a function of disk mass and viewing angle, and we show that the predicted levels are compatible with the observations, for extended reprocessing envelopes of $\sim$1000 gravitational radii. Spectropolarimetry therefore constitutes a new observational test for TDE models, and opens an important new line of exploration in the study of TDEs., Comment: Author's version of paper to appear in Nature Astronomy. In the journal version the detailed discussion on the ISP determination will be moved from the Methods section to a Supplementary Information section. 58 pages in double spacing format, including 5 Figures, 10 Extended Data Figures and 2 Tables

MARCOT Pathfinder is a precursor for MARCOT (Multi Array of Combined Telescopes) at Calar Alto Observatory (CAHA) in Spain. MARCOT is intended to provide CARMENES, currently fiber-fed from the CAHA 3.5m Telescope, with a 5-15m light collecting area from a battery of several tens of small telescopes that are incoherently fed into the final joint single fiber feed of the spectrograph. The modular concept, based on commercially available telescopes, results in cost estimates that are a fraction of the ones for extremely large telescopes (ELT). As a novel approach, MARCOT will employ Multi-Mode Photonic Lanterns (MM-PL) that are being developed as a variant of classical photonic lanterns, to combine the light from the individual telescopes to a single fiber feed to the instrument. This progress report presents the overall concept of MARCOT, the pathfinder telescope and enclosure that is being commissioned at CAHA, the concept of MM-PL, and the next step of installing the Potsdam Multiplex Raman Spectrograph (MRS). MARCOT Pathfinder will be used to validate the conceptual design and predicted performance of MM-PL on sky with a 7 unit telescope prototype., Comment: 10 pages, 7 figures, submitted to SPIE

This work presents a set of studies addressing the use of the low-dispersion slitless prism on Roman for SN spectroscopy as part of the Roman High Latitude Time Domain Survey (HLTDS). We find SN spectral energy distributions including prism data carry more information than imaging alone at fixed total observing time, improving redshift measurements and sub-typing of SNe. The Roman field of view will typically include ~ 10 SNe Ia at observable redshifts at a range of phases (the multiplexing of host galaxies is much greater as they are always present), building up SN spectral time series without targeted observations. We show that fitting these time series extracts more information than stacking the data over all the phases, resulting in a large improvement in precision for SN Ia subclassification measurements. A prism on Roman thus significantly enhances scientific opportunities for the mission, and is particularly important for the Roman SN cosmology program to provide the systematics-controlled measurement that is a focus of the Roman dark energy mission. Optimizing the prism parameters, we conclude that the blue cutoff should be set as blue as the prism image quality allows (~ 7500A), the red cutoff should be set to ~ 18000A to minimize thermal background, and the two-pixel dispersion should be >~ 70., Comment: Submitted to ApJ

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., Comment: THIS VERSION IS AFTER THE JOURNAL REVIEW

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 $\sigma$ level to a $\gamma$-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\"ahovi, OVRO), the optical/UV (ASAS-SN, KVA, REM, Swift/UVOT), the X-rays (Swift/XRT, NuSTAR), the high-energy $\gamma$ rays (Fermi/LAT) and the very-high-energy (VHE) $\gamma$ rays (MAGIC). In $\gamma$ 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$\sigma$ 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 $\gamma$ rays. According to the model presented here, the December 2018 $\gamma$-ray flare was connected to a neutrino emission that was too brief and not bright enough to be detected by current neutrino instruments., Comment: 18 pages, 6 figures; V2 to match published version, typo fixed