Showing 1–30 of 30 results for author: Karjalainen, M

The PLATO Mission

Authors: Heike Rauer, Conny Aerts, Juan Cabrera, Magali Deleuil, Anders Erikson, Laurent Gizon, Mariejo Goupil, Ana Heras, Jose Lorenzo-Alvarez, Filippo Marliani, Cesar Martin-Garcia, J. Miguel Mas-Hesse, Laurence O'Rourke, Hugh Osborn, Isabella Pagano, Giampaolo Piotto, Don Pollacco, Roberto Ragazzoni, Gavin Ramsay, Stéphane Udry, Thierry Appourchaux, Willy Benz, Alexis Brandeker, Manuel Güdel, Eduardo Janot-Pacheco , et al. (801 additional authors not shown)

Abstract: PLATO (PLAnetary Transits and Oscillations of stars) is ESA's M3 mission designed to detect and characterise extrasolar planets and perform asteroseismic monitoring of a large number of stars. PLATO will detect small planets (down to <2 R_(Earth)) around bright stars (<11 mag), including terrestrial planets in the habitable zone of solar-like stars. With the complement of radial velocity observati… ▽ More PLATO (PLAnetary Transits and Oscillations of stars) is ESA's M3 mission designed to detect and characterise extrasolar planets and perform asteroseismic monitoring of a large number of stars. PLATO will detect small planets (down to <2 R_(Earth)) around bright stars (<11 mag), including terrestrial planets in the habitable zone of solar-like stars. With the complement of radial velocity observations from the ground, planets will be characterised for their radius, mass, and age with high accuracy (5 %, 10 %, 10 % for an Earth-Sun combination respectively). PLATO will provide us with a large-scale catalogue of well-characterised small planets up to intermediate orbital periods, relevant for a meaningful comparison to planet formation theories and to better understand planet evolution. It will make possible comparative exoplanetology to place our Solar System planets in a broader context. In parallel, PLATO will study (host) stars using asteroseismology, allowing us to determine the stellar properties with high accuracy, substantially enhancing our knowledge of stellar structure and evolution. The payload instrument consists of 26 cameras with 12cm aperture each. For at least four years, the mission will perform high-precision photometric measurements. Here we review the science objectives, present PLATO's target samples and fields, provide an overview of expected core science performance as well as a description of the instrument and the mission profile at the beginning of the serial production of the flight cameras. PLATO is scheduled for a launch date end 2026. This overview therefore provides a summary of the mission to the community in preparation of the upcoming operational phases. △ Less

Submitted 8 June, 2024; originally announced June 2024.

MELCHIORS: The Mercator Library of High Resolution Stellar Spectroscopy

Authors: P. Royer, T. Merle, K. Dsilva, S. Sekaran, H. Van Winckel, Y. Frémat, M. Van der Swaelmen, S. Gebruers, A. Tkachenko, M. Laverick, M. Dirickx, G. Raskin, H. Hensberge, M. Abdul-Masih, B. Acke, M. L. Alonso, S. Bandhu Mahato, P. G. Beck, N. Behara, S. Bloemen, B. Buysschaert, N. Cox, J. Debosscher, P. De Cat, P. Degroote , et al. (49 additional authors not shown)

Abstract: Over the past decades, libraries of stellar spectra have been used in a large variety of science cases, including as sources of reference spectra for a given object or a given spectral type. Despite the existence of large libraries and the increasing number of projects of large-scale spectral surveys, there is to date only one very high-resolution spectral library offering spectra from a few hundr… ▽ More Over the past decades, libraries of stellar spectra have been used in a large variety of science cases, including as sources of reference spectra for a given object or a given spectral type. Despite the existence of large libraries and the increasing number of projects of large-scale spectral surveys, there is to date only one very high-resolution spectral library offering spectra from a few hundred objects from the southern hemisphere (UVES-POP) . We aim to extend the sample, offering a finer coverage of effective temperatures and surface gravity with a uniform collection of spectra obtained in the northern hemisphere. Between 2010 and 2020, we acquired several thousand echelle spectra of bright stars with the Mercator-HERMES spectrograph located in the Roque de Los Muchachos Observatory in La Palma, whose pipeline offers high-quality data reduction products. We have also developed methods to correct for the instrumental response in order to approach the true shape of the spectral continuum. Additionally, we have devised a normalisation process to provide a homogeneous normalisation of the full spectral range for most of the objects. We present a new spectral library consisting of 3256 spectra covering 2043 stars. It combines high signal-to-noise and high spectral resolution over the entire range of effective temperatures and luminosity classes. The spectra are presented in four versions: raw, corrected from the instrumental response, with and without correction from the atmospheric molecular absorption, and normalised (including the telluric correction). △ Less

Submitted 5 November, 2023; originally announced November 2023.

Comments: 17 pages, 18 figures Preview and access to the library: https://www.royer.se/melchiors.html

Journal ref: A&A 681, A107 (2024)

The GW Vir instability strip in the light of new observations of PG 1159 stars. Discovery of pulsations in the central star of Abell 72 and variability of RX J0122.9-7521

Authors: Paulina Sowicka, Gerald Handler, David Jones, John A. R. Caldwell, Francois van Wyk, Ernst Paunzen, Karolina Bąkowska, Luis Peralta de Arriba, Lucía Suárez-Andrés, Klaus Werner, Marie Karjalainen, Daniel L. Holdsworth

Abstract: We present the results of new time series photometric observations of 29 pre-white dwarf stars of PG 1159 spectral type, carried out in the years 2014-2022. For the majority of stars, a median noise level in Fourier amplitude spectra of 0.5-1.0 mmag was achieved. This allowed the detection of pulsations in the central star of planetary nebula Abell 72, consistent with g-modes excited in GW Vir sta… ▽ More We present the results of new time series photometric observations of 29 pre-white dwarf stars of PG 1159 spectral type, carried out in the years 2014-2022. For the majority of stars, a median noise level in Fourier amplitude spectra of 0.5-1.0 mmag was achieved. This allowed the detection of pulsations in the central star of planetary nebula Abell 72, consistent with g-modes excited in GW Vir stars, and variability in RX J0122.9-7521 that could be due to pulsations, binarity or rotation. For the remaining stars from the sample that were not observed to vary, we placed upper limits for variability. After combination with literature data, our results place the fraction of pulsating PG 1159 stars within the GW Vir instability strip at 36%. An updated list of all known PG 1159 stars is provided, containing astrometric measurements from the recent Gaia DR3 data, as well as information on physical parameters, variability, and nitrogen content. Those data are used to calculate luminosities for all PG 1159 stars to place the whole sample on the theoretical Hertzsprung-Russell diagram for the first time in that way. The pulsating stars are discussed as a group, and arguments are given that the traditional separation of GW Vir pulsators in "DOV" and "PNNV" stars is misleading and should not be used. △ Less

Submitted 28 September, 2023; originally announced September 2023.

Comments: Accepted for publication in ApJS

High-resolution Emission Spectroscopy of the Ultrahot Jupiter KELT-9b: Little Variation in Day- and Nightside Emission Line Contrasts

Authors: Andrew Ridden-Harper, Ernst de Mooij, Ray Jayawardhana, Neale Gibson, Raine Karjalainen, Marie Karjalainen

Abstract: The transmission spectrum of the ultrahot Jupiter KELT-9b ($T_{eq}$ $\sim$ 4000 K) exhibits absorption by several metal species. We searched for atomic and molecular lines in its emission spectrum by observing partial phase curves with the CARMENES spectrograph ($R$ $\sim$ 80,000 $-$ 95,000). We find evidence for emission by Si I in the atmosphere of KELT-9b for the first time. Additionally we fin… ▽ More The transmission spectrum of the ultrahot Jupiter KELT-9b ($T_{eq}$ $\sim$ 4000 K) exhibits absorption by several metal species. We searched for atomic and molecular lines in its emission spectrum by observing partial phase curves with the CARMENES spectrograph ($R$ $\sim$ 80,000 $-$ 95,000). We find evidence for emission by Si I in the atmosphere of KELT-9b for the first time. Additionally we find evidence for emission by Mg I and Ca II, which were previously detected in transmission, and confirmed earlier detections of Fe I emission. Conversely, we find no evidence for dayside emission from Al I, Ca I, Cr I, FeH, Fe II, K I, Li I, Mg II, Na I, OH, Ti I, TiO, V I, V II, VO, and Y I. By employing likelihood mapping, we find indications of there being little variation in emission line contrast between the day- and nightsides $-$suggesting that KELT-9b may harbor iron emission on its nightside. Our results demonstrate that high-resolution ground-based emission spectroscopy can provide valuable insights into exoplanet atmospheres. △ Less

Submitted 6 April, 2023; originally announced April 2023.

Comments: Accepted for publication in AJ

Companions to Kepler giant stars: A long-period eccentric substellar companion to KIC 3526061 and a stellar companion to HD 187878

Authors: Marie Karjalainen, Raine Karjalainen, Artie P. Hatzes, Holger Lehmann, Pierre Kervella, Saskia Hekker, Hans Van Winckel, Jakub Überlauer, Michaela Vítková, Marek Skarka, Petr Kabáth, Saskia Prins, Andrew Tkachenko, William D. Cochran, Alain Jorissen

Abstract: Context. Our knowledge of populations and occurrence of planets orbiting evolved intermediate-mass stars is still incomplete. In 2010 we started a planet-search program among 95 giant stars observed by the Kepler mission to increase the sample of giant stars with planets and with reliable estimates of stellar masses and radii. Aims. We present the two systems KIC 3526061 and HD 187878 from our pla… ▽ More Context. Our knowledge of populations and occurrence of planets orbiting evolved intermediate-mass stars is still incomplete. In 2010 we started a planet-search program among 95 giant stars observed by the Kepler mission to increase the sample of giant stars with planets and with reliable estimates of stellar masses and radii. Aims. We present the two systems KIC 3526061 and HD 187878 from our planet-search program for which we could characterise their companions. Methods. We used precise stellar radial velocity measurements taken with four different echelle spectrographs to derive an orbital solution. We used Gaia astrometric measurements to obtain the inclination of the HD 187878 system and Kepler photometric observations to estimate the stellar mass and radius. Results. We report the discovery of a sub-stellar and a stellar companion around two intermediate-mass red giant branch stars. KIC 3526061 b is most likely a brown dwarf with a minimum mass of 18.15 Jupiter masses in a long-period eccentric orbit, with the orbital period 3552 d and orbital eccentricity 0.85. It is the most evolved system found having a sub-stellar companion with such a large eccentricity and wide separation. HD 187878 B has a minimum mass of 78.4 Jupiter masses. Combining the spectroscopic orbital parameters with the astrometric proper motion anomaly we derived an orbital inclination 9.8 deg, which corresponds to the companion's mass in the stellar regime of 0.51 Sun mass. Conclusions. A sub-stellar companion of KIC 3526061 extends the sample of known red giant branch stars with sub-stellar companions on very eccentric wide orbits and might provide a probe of the dynamical evolution of such systems over time. △ Less

Submitted 22 September, 2022; originally announced September 2022.

Comments: 14 pages, 11 figures. Accepted to A&A

Journal ref: A&A 668, A26 (2022)

Periodic variable A-F spectral type stars in the northern TESS continuous viewing zone

Authors: M. Skarka, J. Žák, M. Fedurco, E. Paunzen, Z. Henzl, M. Mašek, R. Karjalainen, J. P. Sanchez Arias, Á. Sódor, R. F. Auer, P. Kabáth, M. Karjalainen, J. Liška, D. Štegner

Abstract: The goal of our study is to provide a reliable classification of variability of A-F stars brighter than 11 mag located in the northern TESS continuous viewing zone. We also aim at thorough discussion about issues in the classification related to the data characteristics and the issues arising from the similar light curve shape generated by different physical mechanisms. We used TESS long- and shor… ▽ More The goal of our study is to provide a reliable classification of variability of A-F stars brighter than 11 mag located in the northern TESS continuous viewing zone. We also aim at thorough discussion about issues in the classification related to the data characteristics and the issues arising from the similar light curve shape generated by different physical mechanisms. We used TESS long- and short-cadence photometric data and corresponding Fourier transform to classify the variability type of the stars. We present a clear and concise classification system that is demonstrated on many examples. We found clear signs of variability in 3025 of 5923 studied stars (51 %). For 1813 of these 3025 stars, we provide a classification. From the classified stars, 64.5 % are pulsating stars of GDOR and DSCT types and their hybrids. We realised that the long- and short-cadence PDCSAP data can differ significantly not only in amplitudes but also in the content of instrumental/data reduction artifacts making the long-cadence data less reliable. We identified a new group of stars showing stable light curves and characteristic frequency spectra pattern (8.5 % of the classified stars). According to the position in the Hertzsprung-Russell diagram, these stars are likely GDOR stars but are about 200 K cooler than GDORs on average and have smaller amplitudes and longer periods in average. We show that without spectroscopic observations, it can be impossible to unambiguously distinguish between ellipsoidal variability and rotational variability. We also apply our methodology to three previous studies and found significant discrepancies in the classification. We demonstrate how difficult the classification of variable A-F stars can be when using only photometric data. △ Less

Submitted 5 August, 2022; v1 submitted 26 July, 2022; originally announced July 2022.

Comments: Accepted for publication in A&A, 23 pages, 17 figures, 5 tables

Journal ref: A&A 666, A142 (2022)

TOI-2046b, TOI-1181b and TOI-1516b, three new hot Jupiters from \textit{TESS}: planets orbiting a young star, a subgiant and a normal star

Authors: Petr Kabáth, Priyanka Chaturvedi, Phillip J. MacQueen, Marek Skarka, Ján Šubjak, Massimilliano Esposito, William D. Cochran, Salvatore E. Bellomo, Raine Karjalainen, Eike W. Guenther, Michael Endl, Szilárd Csizmadia, Marie Karjalainen, Artie Hatzes, Jiří Žák, Davide Gandolfi, Henri M. J. Boffin, Jose I. Vines, John H. Livingston, Rafael A. García, Savita Mathur, Lucía González-Cuesta, Martin Blažek, Douglas A. Caldwell, Knicole D. Colón , et al. (32 additional authors not shown)

Abstract: We present the confirmation and characterization of three hot Jupiters, TOI-1181b, TOI-1516b, and TOI-2046b, discovered by the TESS space mission. The reported hot Jupiters have orbital periods between 1.4 and 2.05 days. The masses of the three planets are $1.18\pm0.14$ M$_{\mathrm{J}}$, $3.16\pm0.12$\, M$_{\mathrm{J}}$, and 2.30 $\pm 0.28$ M$_{\mathrm{J}}$, for TOI-1181b, TOI-1516b, and TOI-2046b… ▽ More We present the confirmation and characterization of three hot Jupiters, TOI-1181b, TOI-1516b, and TOI-2046b, discovered by the TESS space mission. The reported hot Jupiters have orbital periods between 1.4 and 2.05 days. The masses of the three planets are $1.18\pm0.14$ M$_{\mathrm{J}}$, $3.16\pm0.12$\, M$_{\mathrm{J}}$, and 2.30 $\pm 0.28$ M$_{\mathrm{J}}$, for TOI-1181b, TOI-1516b, and TOI-2046b, respectively. The stellar host of TOI-1181b is a F9IV star, whereas TOI-1516b and TOI-2046b orbit F main sequence host stars. The ages of the first two systems are in the range of 2-5 Gyrs. However, TOI-2046 is among the few youngest known planetary systems hosting a hot Jupiter, with an age estimate of 100-400 Myrs. The main instruments used for the radial velocity follow-up of these three planets are located at Ondřejov, Tautenburg and McDonald Observatory, and all three are mounted on 2-3 meter aperture telescopes, demonstrating that mid-aperture telescope networks can play a substantial role in the follow-up of gas giants discovered by \textit{TESS} and in the future by \textit{PLATO}. △ Less

Submitted 3 May, 2022; originally announced May 2022.

Comments: MNRAS accepted

TOI-1670 b and c: An Inner Sub-Neptune with an Outer Warm Jupiter Unlikely to have Originated from High-Eccentricity Migration

Authors: Quang H. Tran, Brendan P. Bowler, Michael Endl, William D. Cochran, Phillip J. MacQueen, Davide Gandolfi, Carina M. Persson, Malcolm Fridlund, Enric Palle, Grzegorz Nowak, Hans J. Deeg, Rafael Luque, John H. Livingston, Petr Kabáth, Marek Skarka, Ján Šubjak, Steve B. Howell, Simon H. Albrecht, Karen A. Collins, Massimiliano Esposito, Vincent Van Eylen, Sascha Grziwa, Elisa Goffo, Chelsea X. Huang, Jon M. Jenkins , et al. (16 additional authors not shown)

Abstract: We report the discovery of two transiting planets around the bright ($V=9.9$ mag) main sequence F7 star TOI-1670 by the Transiting Exoplanet Survey Satellite. TOI-1670 b is a sub-Neptune ($R_\mathrm{b} = 2.06_{-0.15}^{+0.19}$ $R_\oplus$) on a 10.9-day orbit and TOI-1670 c is a warm Jupiter ($R_\mathrm{c} = 0.987_{-0.025}^{+0.025}$ $R_\mathrm{Jup}$) on a 40.7-day orbit. Using radial velocity observ… ▽ More We report the discovery of two transiting planets around the bright ($V=9.9$ mag) main sequence F7 star TOI-1670 by the Transiting Exoplanet Survey Satellite. TOI-1670 b is a sub-Neptune ($R_\mathrm{b} = 2.06_{-0.15}^{+0.19}$ $R_\oplus$) on a 10.9-day orbit and TOI-1670 c is a warm Jupiter ($R_\mathrm{c} = 0.987_{-0.025}^{+0.025}$ $R_\mathrm{Jup}$) on a 40.7-day orbit. Using radial velocity observations gathered with the Tull coudé Spectrograph on the Harlan J. Smith telescope and HARPS-N on the Telescopio Nazionale Galileo, we find a planet mass of $M_\mathrm{c} = 0.63_{-0.08}^{+0.09}$ $M_\mathrm{Jup}$ for the outer warm Jupiter, implying a mean density of $ρ_c = 0.81_{-0.11}^{+0.13}$ g cm$^{-3}$. The inner sub-Neptune is undetected in our radial velocity data ($M_\mathrm{b} < 0.13$ $M_\mathrm{Jup}$ at the 99% confidence level). Multi-planet systems like TOI-1670 hosting an outer warm Jupiter on a nearly circular orbit ($e_\mathrm{c} = 0.09_{-0.04}^{+0.05}$) and one or more inner coplanar planets are more consistent with "gentle" formation mechanisms such as disk migration or $in$ $situ$ formation rather than high-eccentricity migration. Of the 11 known systems with a warm Jupiter and a smaller inner companion, 8 (73%) are near a low-order mean-motion resonance, which can be a signature of migration. TOI-1670 joins two other systems (27% of this subsample) with period commensurabilities greater than 3, a common feature of $in$ $situ$ formation or halted inward migration. TOI-1670 and the handful of similar systems support a diversity of formation pathways for warm Jupiters. △ Less

Submitted 8 March, 2022; originally announced March 2022.

Comments: 23 pages, 9 figures, accepted for publication in AJ

TOI-1268b: the youngest, hot, Saturn-mass transiting exoplanet

Authors: J. Šubjak, M. Endl, P. Chaturvedi, R. Karjalainen, W. D. Cochran, M. Esposito, D. Gandolfi, K. W. F. Lam, K. Stassun, J. Žák, N. Lodieu, H. M. J. Boffin, P. J. MacQueen, A. Hatzes, E. W. Guenther, I. Georgieva, S. Grziwa, H. Schmerling, M. Skarka, M. Blažek, M. Karjalainen, M. Špoková, H. Isaacson, A. W. Howard, C. J. Burke , et al. (19 additional authors not shown)

Abstract: We report the discovery of TOI-1268b, a transiting Saturn-mass planet from the TESS space mission. With an age of less than one Gyr, derived from various age indicators, TOI-1268b is the youngest Saturn-mass planet known to date and contributes to the small sample of well characterised young planets. It has an orbital period of $P\,=\,8.1577080\pm0.0000044$ days, and transits an early K dwarf star… ▽ More We report the discovery of TOI-1268b, a transiting Saturn-mass planet from the TESS space mission. With an age of less than one Gyr, derived from various age indicators, TOI-1268b is the youngest Saturn-mass planet known to date and contributes to the small sample of well characterised young planets. It has an orbital period of $P\,=\,8.1577080\pm0.0000044$ days, and transits an early K dwarf star with a mass of $M_\star$ = $ 0.96 \pm 0.04$ $M_{\odot}$, a radius of $R_\star$ = $ 0.92 \pm 0.06$ $R_{\odot}$, an effective temperature of $T_\mathrm{eff}\,=\,5300\pm100$ K, and a metallicity of $0.36\pm0.06$ dex. By combining TESS photometry with high-resolution spectra acquired with the Tull spectrograph at McDonald observatory, and the high-resolution spectrographs at Tautenburg and Ondrejov observatories, we measured a planetary mass of $M_\mathrm{p}\,=\,96.4 \pm 8.3\,M_{\oplus}$ and a radius of $R_\mathrm{p}\,=\,9.1 \pm 0.6\,R_{\oplus}$. TOI-1268 is an ideal system to study the role of star-planet tidal interactions for non-inflated Saturn-mass planets. We used system parameters derived in this paper to constrain the planet tidal quality factor to the range of $10^{4.5-5.3}$. When compared with the sample of other non-inflated Saturn-mass planets, TOI-1268b is one of the best candidates for transmission spectroscopy studies. △ Less

Submitted 23 February, 2022; v1 submitted 31 January, 2022; originally announced January 2022.

Comments: 21 pages, 17 figures, accepted for publication in Astronomy & Astrophysics; see independent work by Dong et al. for Rossiter-McLaughlin measurement of TOI-1268b

Journal ref: A&A 662, A107 (2022)

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

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

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

Submitted 29 November, 2021; originally announced November 2021.

Comments: 17 pages, 5 figures

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

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

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

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

Submitted 25 November, 2021; originally announced November 2021.

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

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

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

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

Submitted 24 September, 2021; originally announced September 2021.

Comments: Accepted for publication by The Astrophysical Journal

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

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

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

Submitted 25 August, 2021; originally announced August 2021.

Comments: 19 pages, 9 figures, Accepted to ApJ

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

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

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

Submitted 20 July, 2021; originally announced July 2021.

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

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

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

Submitted 10 June, 2021; originally announced June 2021.

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

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

A Near-Infrared Chemical Inventory of the Atmosphere of 55 Cancri e

Authors: Emily K. Deibert, Ernst J. W. de Mooij, Ray Jayawardhana, Andrew Ridden-Harper, Suresh Sivanandam, Raine Karjalainen, Marie Karjalainen

Abstract: We present high-resolution near-infrared spectra taken during eight transits of 55 Cancri e, a nearby low-density super-Earth with a short orbital period (< 18 hours). While this exoplanet's bulk density indicates a possible atmosphere, one has not been detected definitively. Our analysis relies on the Doppler cross-correlation technique, which takes advantage of the high spectral resolution and b… ▽ More We present high-resolution near-infrared spectra taken during eight transits of 55 Cancri e, a nearby low-density super-Earth with a short orbital period (< 18 hours). While this exoplanet's bulk density indicates a possible atmosphere, one has not been detected definitively. Our analysis relies on the Doppler cross-correlation technique, which takes advantage of the high spectral resolution and broad wavelength coverage of our data, to search for the thousands of absorption features from hydrogen-, carbon-, and nitrogen-rich molecular species in the planetary atmosphere. Although we are unable to detect an atmosphere around 55 Cancri e, we do place strong constraints on the levels of HCN, NH${}_3$, and C${}_2$H${}_2$ that may be present. In particular, at a mean molecular weight of 5 amu we can rule out the presence of HCN in the atmosphere down to a volume mixing ratio (VMR) of 0.02%, NH${}_3$ down to a VMR of 0.08%, and C${}_2$H${}_2$ down to a VMR of 1.0%. If the mean molecular weight is relaxed to 2 amu, we can rule out the presence of HCN, NH${}_3$, and C${}_2$H${}_2$ down to VMRs of 0.001%, 0.0025%, and 0.08% respectively. Our results reduce the parameter space of possible atmospheres consistent with the analysis of HST/WFC3 observations by Tsiaras et al. (2016), and indicate that if 55 Cancri e harbors an atmosphere, it must have a high mean molecular weight and/or clouds. △ Less

Submitted 17 February, 2021; originally announced February 2021.

Comments: 32 pages, 34 figures. Accepted for publication in AJ

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

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

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

Submitted 21 December, 2020; originally announced December 2020.

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

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

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

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

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

Submitted 18 December, 2020; originally announced December 2020.

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

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

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

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

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

Submitted 13 December, 2020; originally announced December 2020.

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

Journal ref: ApJ 908, 90 (2021)

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

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

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

Submitted 2 December, 2020; originally announced December 2020.

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

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

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

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

Submitted 20 November, 2020; originally announced November 2020.

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

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

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

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

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

Submitted 29 October, 2020; originally announced October 2020.

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

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

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

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

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

Submitted 17 July, 2020; originally announced July 2020.

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

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

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

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

Submitted 8 June, 2020; originally announced June 2020.

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

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

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

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

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

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

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

Detection of ionized calcium in the atmosphere of the ultra-hot Jupiter KELT-9b

Authors: Jake D. Turner, Ernst J. W. de Mooij, Ray Jayawardhana, Mitchell E. Young, Luca Fossati, Tommi Koskinen, Joshua D. Lothringer, Raine Karjalainen, Marie Karjalainen

Abstract: With a day-side temperature in excess of 4500K, comparable to a mid-K-type star, KELT-9b is the hottest planet known. Its extreme temperature makes KELT-9b a particularly interesting test bed for investigating the nature and diversity of gas giant planets. We observed the transit of KELT-9b at high spectral resolution (R$\sim$94,600) with the CARMENES instrument on the Calar Alto 3.5-m telescope.… ▽ More With a day-side temperature in excess of 4500K, comparable to a mid-K-type star, KELT-9b is the hottest planet known. Its extreme temperature makes KELT-9b a particularly interesting test bed for investigating the nature and diversity of gas giant planets. We observed the transit of KELT-9b at high spectral resolution (R$\sim$94,600) with the CARMENES instrument on the Calar Alto 3.5-m telescope. Using these data, we detect for the first time ionized calcium (CaII triplet) absorption in the atmosphere of KELT-9b; this is the second time that CaII has been observed in a hot Jupiter. Our observations also reveal prominent H$α$ absorption, confirming the presence of an extended hydrogen envelope around KELT-9b. We compare our detections with an atmospheric model and find that all four lines form between atmospheric temperatures of 6100 K and 8000 K and that the CaII lines form at pressures between 10 and 50 nbar while the H$α$ line forms at a lower pressure ($\sim$6 nbar), higher up in the atmosphere. The altitude that the core of H$α$ line forms is found to be $\sim$1.4 R$_{p}$, well within the planetary Roche lobe ($\sim$1.9 R$_{p}$). Therefore, rather than probing the escaping upper atmosphere directly, the H$α$ line and the other observed Balmer and metal lines serve as atmospheric thermometers enabling us to probe the planet's temperature profile, thus energy budget. △ Less

Submitted 18 December, 2019; v1 submitted 15 October, 2019; originally announced October 2019.

Comments: Accepted to ApJL (December 11 2019), 11 pages, 3 figures, 2 tables, 2 Appendices (with 2 figures)

Time resolved spectroscopy of dust and gas from extrasolar planetesimals orbiting WD 1145+017

Authors: Marie Karjalainen, Ernst J. W. de Mooij, Raine Karjalainen, Neale P. Gibson

Abstract: Multiple long and variable transits caused by dust from possibly disintegrating asteroids were detected in light curves of WD 1145+017. We present time-resolved spectroscopic observations of this target with QUCAM CCDs mounted in the Intermediate dispersion Spectrograph and Imaging System at the 4.2-m William Herschel Telescope in two different spectral arms: the blue arm covering 3800-4025 Å and… ▽ More Multiple long and variable transits caused by dust from possibly disintegrating asteroids were detected in light curves of WD 1145+017. We present time-resolved spectroscopic observations of this target with QUCAM CCDs mounted in the Intermediate dispersion Spectrograph and Imaging System at the 4.2-m William Herschel Telescope in two different spectral arms: the blue arm covering 3800-4025 Å and the red arm covering 7000-7430 Å. When comparing individual transits in both arms, our observations show with 20 σ significance an evident colour difference between the in- and out-of-transit data of the order of 0.05-0.1 mag, where transits are deeper in the red arm. We also show with > 6 σ significance that spectral lines in the blue arm are shallower during transits than out-of-transit. For the circumstellar lines it also appears that during transits the reduction in absorption is larger on the red side of the spectral profiles. Our results confirm previous findings showing the u'-band excess and a decrease in line absorption during transits. Both can be explained by an opaque body blocking a fraction of the gas disc causing the absorption, implying that the absorbing gas is between the white dwarf and the transiting objects. Our results also demonstrate the capability of EMCCDs to perform high-quality time resolved spectroscopy of relatively faint targets. △ Less

Submitted 10 October, 2018; originally announced October 2018.

Comments: 9 pages, 5 figures. Accepted to MNRAS

Occultations from an active accretion disk in a 72 day detached post-Algol system detected by K2

Authors: G. Zhou, S. Rappaport, L. Nelson, C. X. Huang, A. Senhadji, J. E. Rodriguez, A. Vanderburg, S. Quinn, C. I. Johnson, D. W. Latham, G. Torres, B. L. Gary, T. G. Tan, M. C. Johnson, J. Burt, M. H. Kristiansen, T. L. Jacobs, D. LaCourse, H. M. Schwengeler, I. Terentev, A. Bieryla, G. A. Esquerdo, P. Berlind, M. L. Calkins, J. Bento , et al. (6 additional authors not shown)

Abstract: Disks in binary systems can cause exotic eclipsing events. MWC 882 (BD-22 4376, EPIC 225300403) is such a disk-eclipsing system identified from observations during Campaign 11 of the K2 mission. We propose that MWC 882 is a post-Algol system with a B7 donor star of mass $0.542\pm0.053\,M_\odot$ in a 72 day period orbit around an A0 accreting star of mass $3.24\pm0.29\,M_\odot$. The… ▽ More Disks in binary systems can cause exotic eclipsing events. MWC 882 (BD-22 4376, EPIC 225300403) is such a disk-eclipsing system identified from observations during Campaign 11 of the K2 mission. We propose that MWC 882 is a post-Algol system with a B7 donor star of mass $0.542\pm0.053\,M_\odot$ in a 72 day period orbit around an A0 accreting star of mass $3.24\pm0.29\,M_\odot$. The $59.9\pm6.2\,R_\odot$ disk around the accreting star occults the donor star once every orbit, inducing 19 day long, 7% deep eclipses identified by K2, and subsequently found in pre-discovery ASAS and ASAS-SN observations. We coordinated a campaign of photometric and spectroscopic observations for MWC 882 to measure the dynamical masses of the components and to monitor the system during eclipse. We found the photometric eclipse to be gray to $\approx 1$%. We found the primary star exhibits spectroscopic signatures of active accretion, and observed gas absorption features from the disk during eclipse. We suggest MWC 882 initially consisted of a $\approx 3.6\,M_\odot$ donor star transferring mass via Roche lobe overflow to a $\approx 2.1\,M_\odot$ accretor in a $\approx 7$ day initial orbit. Through angular momentum conservation, the donor star is pushed outward during mass transfer to its current orbit of 72 days. The observed state of the system corresponds with the donor star having left the Red Giant Branch ~0.3 Myr ago, terminating active mass transfer. The present disk is expected to be short-lived ($10^2$ years) without an active feeding mechanism, presenting a challenge to this model. △ Less

Submitted 18 January, 2018; originally announced January 2018.

Comments: 19 pages, accepted for publication in ApJ

The Radial Velocity Variability of the K-giant Gamma Draconis: Stellar Variability Masquerading as a Planet

Authors: A. P. Hatzes, M. Endl, W. D. Cochran, P. J. MacQueen, I. Han, B. -C. Lee, K. -M. Kim, D. Mkrtichian, M. Doellinger, M. Hartmann, M. Karjalainen, S. Dreizler

Abstract: We present precise stellar radial velocity measurements of Gamma Dra taken from 2003 to 2017. The data from 2003 to 2011 show coherent, long-lived variations with a period of 702 d. These variations are consistent with the presence of a planetary companion having m sin i = 10.7 M_Jup whose orbital properties are typical for giant planets found around evolved stars. An analysis of the Hipparcos pho… ▽ More We present precise stellar radial velocity measurements of Gamma Dra taken from 2003 to 2017. The data from 2003 to 2011 show coherent, long-lived variations with a period of 702 d. These variations are consistent with the presence of a planetary companion having m sin i = 10.7 M_Jup whose orbital properties are typical for giant planets found around evolved stars. An analysis of the Hipparcos photometry, Ca II S-index measurements, and measurements of the spectral line shapes during this time show no variations with the radial velocity of the planet which seems to "confirm"' the presence of the planet. However, radial velocity measurements taken 2011 -- 2017 seem to refute this. From 2011 to 2013 the radial velocity variations virtually disappear only to return in 2014, but with a noticeable phase shift. The total radial velocity variations are consistent either with amplitude variations on timescales of ~ 10.6 yr, or the beating effect between two periods of 666 d and 801 d. It seems unlikely that both these signals stem from a two-planet system. A simple dynamical analysis indicates that there is only a 1-2 % chance that the two-planet is stable. Rather, we suggest that this multi-periodic behavior may represent a new form of stellar variability, possibly related to oscillatory convective modes. If such intrinsic stellar variability is common around K giant stars and is attributed to planetary companions, then the planet occurrence rate among these stars may be significantly lower than thought. △ Less

Submitted 16 January, 2018; originally announced January 2018.

Comments: 37 pages, 9 figures. Accepted for publication in The Astronomical Journal

Vortices in equilibrium scalar electrodynamics

Authors: A. Rajantie, K. Kajantie, M. Karjalainen, M. Laine, J. Peisa

Abstract: Scalar electrodynamics can be used to investigate the formation of cosmic strings in the early universe. We present the results of lattice Monte Carlo simulations of an effective three-dimensional U(1)+Higgs theory that describes the equilibrium properties of finite-temperature scalar electrodynamics near the transition. A gauge-invariant criterion for the existence of a vortex is used in measur… ▽ More Scalar electrodynamics can be used to investigate the formation of cosmic strings in the early universe. We present the results of lattice Monte Carlo simulations of an effective three-dimensional U(1)+Higgs theory that describes the equilibrium properties of finite-temperature scalar electrodynamics near the transition. A gauge-invariant criterion for the existence of a vortex is used in measuring the properties of the vortex network in the equilibrium state both in the Coulomb and in the Higgs phase of the system. The naive definition of the vortex density becomes meaningless in the continuum limit and special care is needed in extracting physical quantities. Numerical evidence for a physical discontinuity in the vortex density is given. △ Less

Submitted 30 July, 1998; originally announced July 1998.

Comments: 4 pages. Talk given by A. Rajantie at PASCOS-98, March 1998

Report number: HIP-1998-47/TH