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.

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

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.

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

The Cherenkov Telescope Array (CTA) is the next generation very high energy gamma-ray observatory covering the 20 GeV - 300 TeV energy range with unprecedented sensitivity, angular and energy resolution. With a site in each hemisphere, CTA will provide full-sky coverage. Four Large Size Telescopes (LSTs) in each site will be dedicated to the lowest energy range (20 GeV - 200 GeV). The first LST prototype has been installed at the CTA Northern site (Canary Island of La Palma, Spain) in October 2018 and it had been since then in commissioning phase. LST1 is located at about 100 m from MAGIC, a system of two 17m-diameter Imaging Atmospheric Cherenkov Telescopes designed to perform gamma-ray astronomy in the energy range from 50 GeV with standard trigger (30 GeV with SumTrigger) to 50 TeV and whose performance is very well established. The co-location of LST1 and MAGIC offers the great opportunity of cross-calibrating the two systems on an event-by-event basis. It will be indeed possible to compare the parameters of the same extensive air shower reconstructed by the two instruments. We investigated the performance that could be reached with combined observations.

The Cherenkov Telescope Array (CTA) is the major next-generation observatory for ground-based very-high-energy gamma-ray astronomy. It will improve the sensitivity of current ground-based instruments by a factor of five to twenty, depending on the energy, greatly improving both their angular and energy resolutions over four decades in energy (from 20 GeV to 300 TeV). This achievement will be possible by using tens of imaging Cherenkov telescopes of three successive sizes. They will be arranged into two arrays, one per hemisphere, located on the La Palma island (Spain) and in Paranal (Chile). We present here the optimised and final telescope arrays for both CTA sites, as well as their foreseen performance, resulting from the analysis of three different large-scale Monte Carlo productions.

We investigate the properties of X-ray emission from shock breakout of a supernova in a stellar wind. We consider a simple model describing aspherical explosions, in which the shock front with an ellipsoidal shape propagates into the dense circumstellar matter. For this model, both X-ray light curves and spectra are simultaneously calculated using a Monte Carlo method. We show that the shock breakout occurs simultaneously in all directions in a steady and spherically symmetric wind. As a result, even for the aspherical explosion, the rise and decay timescales of the light curve do not significantly depend on the viewing angles. This fact suggests that the light curve of the shock breakout may be used as a probe of the wind mass loss rate. We compare our results with the observed spectrum and light curve of XRO 080109/SN 2008D. The observation can be reproduced by an explosion with a shock velocity of 60% of the speed of light and a circumstellar matter with a mass loss rate of 5.e-4 Msun/yr.

We have constructed the database of stars in the local group using the extended version of the SAGA (Stellar Abundances for Galactic Archaeology) database that contains stars in 24 dwarf spheroidal galaxies and ultra faint dwarfs. The new version of the database includes more than 4500 stars in the Milky Way, by removing the previous metallicity criterion of [Fe/H] <= -2.5, and more than 6000 stars in the local group galaxies. We examined a validity of using a combined data set for elemental abundances. We also checked a consistency between the derived distances to individual stars and those to galaxies in the literature values. Using the updated database, the characteristics of stars in dwarf galaxies are discussed. Our statistical analyses of alpha-element abundances show that the change of the slope of the [alpha/Fe] relative to [Fe/H] (so-called "knee") occurs at [Fe/H] = -1.0+-0.1 for the Milky Way. The knee positions for selected galaxies are derived by applying the same method. Star formation history of individual galaxies are explored using the slope of the cumulative metallicity distribution function. Radial gradients along the four directions are inspected in six galaxies where we find no direction dependence of metallicity gradients along the major and minor axes. The compilation of all the available data shows a lack of CEMP-s population in dwarf galaxies, while there may be some CEMP-no stars at [Fe/H] <~ -3 even in the very small sample. The inspection of the relationship between Eu and Ba abundances confirms an anomalously Ba-rich population in Fornax, which indicates a pre-enrichment of interstellar gas with r-process elements. We do not find any evidence of anti-correlations in O-Na and Mg-Al abundances, which characterises the abundance trends in the Galactic globular clusters.

We investigate the mechanism to reproduce notable spectral features at the ignition phase of nova explosion observed for a super-Eddington X-ray transient source MAXI J0158$-$744 in the Small Magellanic Cloud. These are a strong Ne IX emission line at 0.92 keV with a large equivalent width of $0.32^{+0.21}_{-0.11}$ keV and the absence of Ne X line at 1.02 keV. In this paper, we calculate the radiative transfer using a Monte Carlo code, taking into account the line blanketing effect due to transitions of N, O, Ne, Mg and Al ions in an accelerating wind emanating from a white dwarf with a structure based on a spherically symmetric stationary model. We found that the strong Ne IX line can be reproduced if the mass fraction of Ne is enhanced to $10^{-3}$ or more and that of O is reduced to $\sim5\times10^{-9}$ or less and that the absence of other lines including Ne X ions at higher energies can be also reproduced by the line blanketing effect. This enhancement of the Ne mass fraction indicates that the ejecta are enriched by Ne dredged up from the surface of the white dwarf composed of O, Ne, and Mg rather than C and O, as already pointed out in the previous work. We argue that the CNO cycle driving this nova explosion converted most of C and O into N and thus reduced the O mass fraction.

We present theoretical expectations for non-thermal emission due to the bulk Comptonization at the ultra-relativistic shock breakout. We calculate the transfer of photons emitted from the shocked matter with a Monte Carlo code fully taking into account special relativity. As a hydrodynamical model, we use a self-similar solution of Nakayama & Shigeyama (2005). Our calculations reveal that the spectral shape exhibits a double peak or a single peak depending on the shock temperature at the shock breakout. If it is significantly smaller than the rest energy of an electron, the spectrum has a double peak. We also display a few example of light curves, and estimate the total radiation energy. In comparison with observations of gamma-ray bursts, a part of the higher energy component in the spectra and the total energy can be reproduced by some parameter sets. Meanwhile, the lower energy counterpart in the Band function is not reproduced by our results and the duration time seems too short to represent the entire event of a gamma-ray burst. Therefore the subsequent phase will constitute the lower energy part in the spectrum.