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

We present a systematic analysis of the reaction of dark matter distribution to galaxy formation across more than eight orders of magnitude in stellar mass. We extend the previous work presented in the NIHAO-IV paper (Tollet et al.) by adding 46 new high resolution simulations of massive galaxies performed with the inclusion of Black Hole feedback. We show that outflows generated by the AGN are able to partially counteract the dark matter contraction due to the large central stellar component in massive haloes. The net effect is to relax the central dark matter distribution that moves to a less cuspy density profiles at halo masses larger than $\approx 3 \times 10^{12}$ M$_{\odot}$. The scatter around the mean value of the density profile slope ($\alpha$) is fairly constant ($\Delta\alpha \approx 0.3$), with the exception of galaxies with halo masses around $10^{12}$ M$_{\odot}$, at the transition from stellar to AGN feedback dominated systems, where the scatter increases by almost a factor three. We provide useful fitting formulas for the slope of the dark matter density profiles at few percent of the virial radius for the whole stellar mass range: $10^5-10^{12}$ M$_{\odot}$ ($2 \times 10^9-5 \times 10^{13}$ M$_{\odot}$ in halo mass).

Mrk501 is a very high-energy (VHE) gamma-ray blazar located at z=0.034. During a period of two weeks in July 2014, the highest X-ray activity of Mrk501 was observed in ~14 years of operation of the Neil Gehrels Swift Gamma-ray Burst Observatory. We characterize the broadband variability of Mrk501 from radio to VHE gamma rays, and evaluate whether it can be interpreted within theoretical scenarios widely used to explain the broadband emission from blazars. The temporal evolution of the most prominent and variable segments of the SED is described with a one-zone synchrotron self-Compton model with variations in the break energy of the electron energy distribution (EED), and with some adjustments in the magnetic field strength and spectral shape of the EED. A narrow feature at ~3 TeV was observed in the VHE spectrum measured on 2014 July 19 (MJD 56857.98), which is the day with the highest X-ray flux ($>0.3$ keV) measured during the entire Swift mission. This feature is inconsistent with the classical analytic functions to describe the measured VHE spectra (power law, log-parabola, and log-parabola with exponential cutoff) at more than 3$\sigma$. A fit with a log-parabola plus a narrow component is preferred over the fit with a single log-parabola at more than 4$\sigma$, and a dedicated Monte Carlo simulation estimated the significance of this extra component to be larger than 3$\sigma$. Under the assumption that this VHE spectral feature is real, we show that it can be reproduced with three distinct theoretical scenarios: a) a pileup in the EED due to stochastic acceleration; b) a structured jet with two-SSC emitting regions, with one region dominated by an extremely narrow EED; and c) an emission from an IC pair cascade.

Mrk 421 is a high-synchrotron-peaked blazar featuring bright and persistent GeV and TeV emission. We use the longest and densest ongoing unbiased observing campaign obtained at TeV and GeV energies during 5.5 years with the FACT telescope and the Fermi-LAT detector. The contemporaneous multi-wavelength observations were used to characterize the variability of the source and to constrain the underlying physical mechanisms. We study and correlate light curves obtained by nine different instruments from radio to gamma rays and found two significant results. The TeV and X-ray light curves are very well correlated with lag, if any, shorter than a day. The GeV light curve varies independently and accurately leads the variations observed at long wavelengths, in particular in the radio band. We find that the observations match the predictions of leptonic models and suggest that the physical conditions vary along the jet, when the emitting region moves outwards.

We introduce algorithms for black hole physics, i.e., black hole formation, accretion and feedback, into the NIHAO (Numerical Investigation of a Hundred Astrophysical Objects) project of galaxy simulations. This enables us to study high mass, elliptical galaxies, where feedback from the central black hole is generally thought to have a significant effect on their evolution. We furthermore extend the NIHAO suite by 45 simulations that encompass $z=0$ halo masses from $1 \times 10^{12}$ to $4 \times 10^{13}\,\mathrm{M}_{\odot}$, and resimulate five galaxies from the original NIHAO sample with black hole physics, which have $z=0$ halo masses from $8 \times 10^{11}$ to $3 \times 10^{12}\,\mathrm{M}_{\odot}$. Now NIHAO contains 144 different galaxies and thus has the largest sample of zoom-in simulations of galaxies, spanning $z=0$ halo masses from $9 \times 10^{8}$ to $4 \times 10^{13}\,\mathrm{M}_{\odot}$. In this paper we focus on testing the algorithms and calibrating their free parameters against the stellar mass versus halo mass relation and the black hole mass versus stellar mass relation. We also investigate the scatter of these relations, which we find is a decreasing function with time and thus in agreement with observations. For our fiducial choice of parameters we successfully quench star formation in objects above a $z=0$ halo mass of $10^{12}\,\mathrm{M}_{\odot}$, thus transforming them into red and dead galaxies.

We implement an optically thin approximation for the effects of the local radiation field from stars and hot gas on the gas heating and cooling in the N-body SPH code GASOLINE2. We resimulate three galaxies from the NIHAO project: one dwarf, one Milky Way-like and one massive spiral, and study what are the local radiation field effects on various galaxy properties. We also study the effects of varying the Ultra Violet Background (UVB) model, by running the same galaxies with two different UVBs. Galaxy properties at $z=0$ like stellar mass, stellar effective mass radius, HI mass, and radial extent of the HI disc, show significant changes between the models with and without the local radiation field, and smaller differences between the two UVB models. The intrinsic effect of the local radiation field through cosmic time is to increase the equilibrium temperature at the interface between the galaxies and their circumgalactic media (CGM), moving this boundary inwards, while leaving relatively unchanged the gas inflow rate. Consequently, the temperature of the inflow increases when considering the local radiation sources. This temperature increase is a function of total galaxy mass, with a median CGM temperature difference of one order of magnitude for the massive spiral. The local radiation field suppresses the stellar mass growth by $\sim$20 per cent by $z=0$ for all three galaxies, while the HI mass is roughly halfed. The differences in the gas phase diagrams, significantly impact the HI column densities, shifting their peaks in the distributions towards lower $N_{\rm HI}$.

Mrk 421 is a high-synchrotron-peaked blazar featuring bright and persistent GeV and TeV emission. We use multi-wavelength light curves of Mrk 421 spanning 5.5 years with FACT (TeV) and Fermi LAT (GeV) in the gamma rays, Swift BAT, Swift XRT and MAXI in the X-rays, together with optical and radio data and investigate the physical processes driving the emission and variability. Observations by FACT are continuous and not triggered, so the source was found in a wide range of flux states and more than 30 flares were identified from X-rays to TeV. The light curves in TeV and X-rays feature very similar flares with rise and decay times of a few days and zero lag, characteristic for electron processes. At least two parameters per flare, the amplitude and the cut-off energy, are required to explain the observed variability. In addition, the GeV light curve leads and is strongly correlated with the optical and radio light curves as expected from SSC emitting shock propagating in a conical jet.

We use thirty-eight high-resolution simulations of galaxy formation between redshift 10 and 5 to study the impact of a 3 keV warm dark matter (WDM) candidate on the high-redshift Universe. We focus our attention on the stellar mass function and the global star formation rate and consider the consequences for reionization, namely the neutral hydrogen fraction evolution and the electron scattering optical depth. We find that three different effects contribute to differentiate warm and cold dark matter (CDM) predictions: WDM suppresses the number of haloes with mass less than few $10^9$ M$_{\odot}$; at a fixed halo mass, WDM produces fewer stars than CDM; and finally at halo masses below $10^9$ M$_{\odot}$, WDM has a larger fraction of dark haloes than CDM post-reionization. These three effects combine to produce a lower stellar mass function in WDM for galaxies with stellar masses at and below $\sim 10^7$ M$_{\odot}$. For $z > 7$, the global star formation density is lower by a factor of two in the WDM scenario, and for a fixed escape fraction, the fraction of neutral hydrogen is higher by 0.3 at $z \sim 6$. This latter quantity can be partially reconciled with CDM and observations only by increasing the escape fraction from 23 per cent to 34 per cent. Overall, our study shows that galaxy formation simulations at high redshift are a key tool to differentiate between dark matter candidates given a model for baryonic physics.

We use cosmological hydrodynamical galaxy formation simulations from the NIHAO project to investigate the impact of the threshold for star formation on the response of the dark matter (DM) halo to baryonic processes. The fiducial NIHAO threshold, $n=10\, {\rm cm}^{-3}$, results in strong expansion of the DM halo in galaxies with stellar masses in the range $10^{7.5} < M_{star} < 10^{9.5} M_{\odot}$. We find that lower thresholds such as $n=0.1$ (as employed by the EAGLE/APOSTLE and Illustris/AURIGA projects) do not result in significant halo expansion at any mass scale. Halo expansion driven by supernova feedback requires significant fluctuations in the local gas fraction on sub-dynamical times (i.e., < 50 Myr at galaxy half-light radii), which are themselves caused by variability in the star formation rate. At one per cent of the virial radius, simulations with $n=10$ have gas fractions of $\simeq 0.2$ and variations of $\simeq 0.1$, while $n=0.1$ simulations have order of magnitude lower gas fractions and hence do not expand the halo. The observed DM circular velocities of nearby dwarf galaxies are inconsistent with CDM simulations with $n=0.1$ and $n=1$, but in reasonable agreement with $n=10$. Star formation rates are more variable for higher $n$, lower galaxy masses, and when star formation is measured on shorter time scales. For example, simulations with $n=10$ have up to 0.4 dex higher scatter in specific star formation rates than simulations with $n=0.1$. Thus observationally constraining the sub-grid model for star formation, and hence the nature of DM, should be possible in the near future.

In this third paper of the series, we investigate the effects of warm dark matter with a particle mass of $m_\mathrm{WDM}=3\,\mathrm{keV}$ on the smallest galaxies in our Universe. We present a sample of 21 hydrodynamical cosmological simulations of dwarf galaxies and 20 simulations of satellite-host galaxy interaction that we performed both in a Cold Dark Matter (CDM) and Warm Dark Matter (WDM) scenario. In the WDM simulations, we observe a higher critical mass for the onset of star formation. Structure growth is delayed in WDM, as a result WDM haloes have a stellar population on average two Gyrs younger than their CDM counterparts. Nevertheless, despite this delayed star formation, CDM and WDM galaxies are both able to reproduce the observed scaling relations for velocity dispersion, stellar mass, size, and metallicity at $z=0$. WDM satellite haloes in a Milky Way mass host are more susceptible to tidal stripping due to their lower concentrations, but their galaxies can even survive longer than the CDM counterparts if they live in a dark matter halo with a steeper central slope. In agreement with our previous CDM satellite study we observe a steepening of the WDM satellites' central dark matter density slope due to stripping. The difference in the average stellar age for satellite galaxies, between CDM and WDM, could be used in the future for disentangling these two models.

Flux distribution is an important tool to understand the variability processes in active galactic nuclei. We now have available a great deal of observational evidences pointing towards the presence of log-normal components in the high energy light curves, and different models have been proposed to explain these data. Here, we collect some of the recent developments on this topic using the well-known blazar Mrk 501 as example of complex and interesting aspects coming from its flux distribution in different energy ranges and at different timescales. The observational data we refer to are those collected in a complementary manner by Fermi-LAT over multiple years, and by the First G-APD Cherenkov Telescope (FACT) telescope and the H.E.S.S. array in correspondence of the bright flare of June 2014.

A multiwavelength campaign was organized to take place between March and July of 2012. Excellent temporal coverage was obtained with more than 25 instruments, including the MAGIC, FACT and VERITAS Cherenkov telescopes, the instruments on board the Swift and Fermi spacecraft, and the telescopes operated by the GASP-WEBT collaboration. Mrk 501 showed a very high energy (VHE) gamma-ray flux above 0.2 TeV of $\sim$0.5 times the Crab Nebula flux (CU) for most of the campaign. The highest activity occurred on 2012 June 9, when the VHE flux was $\sim$3 CU, and the peak of the high-energy spectral component was found to be at $\sim$2 TeV. This study reports very hard X-ray spectra, and the hardest VHE spectra measured to date for Mrk 501. The fractional variability was found to increase with energy, with the highest variability occurring at VHE, and a significant correlation between the X-ray and VHE bands. The unprecedentedly hard X-ray and VHE spectra measured imply that their low- and high-energy components peaked above 5 keV and 0.5 TeV, respectively, during a large fraction of the observing campaign, and hence that Mrk 501 behaved like an extreme high-frequency- peaked blazar (EHBL) throughout the 2012 observing season. This suggests that being an EHBL may not be a permanent characteristic of a blazar, but rather a state which may change over time. The one-zone synchrotron self-Compton (SSC) scenario can successfully describe the segments of the SED where most energy is emitted, with a significant correlation between the electron energy density and the VHE gamma-ray activity, suggesting that most of the variability may be explained by the injection of high-energy electrons. The one-zone SSC scenario used reproduces the behaviour seen between the measured X-ray and VHE gamma-ray fluxes, and predicts that the correlation becomes stronger with increasing energy of the X-rays.

In the first paper we presented 27 hydrodynamical cosmological simulations of galaxies with total masses between $5 \times 10^8$ and $10^{10}\,\mathrm{M}_\odot$. In this second paper we use a subset of these cosmological simulations as initial conditions (ICs) for more than forty hydrodynamical simulations of satellite and host galaxy interaction. Our cosmological ICs seem to suggest that galaxies on these mass scales have very little rotational support and are velocity dispersion ($\sigma$) dominated. Accretion and environmental effects increase the scatter in the galaxy scaling relations (e.g. size - velocity dispersion) in very good agreement with observations. Star formation is substantially quenched after accretion. Mass removal due to tidal forces has several effects: it creates a very flat stellar velocity dispersion profiles, and it reduces the dark matter content at all scales (even in the centre), which in turn lowers the stellar velocity on scales around 0.5 kpc even when the galaxy does not lose stellar mass. Satellites that start with a cored dark matter profile are more prone to either be destroyed or to end up in a very dark matter poor galaxy. Finally, we found that tidal effects always increase the "cuspyness" of the dark matter profile, even for haloes that infall with a core.

The satellites of the Milky Way and Andromeda represent the smallest galaxies we can observe in our Universe. In this series of papers we aim to shed light on their formation and evolution using cosmological hydrodynamical simulations. In this first paper we focus on the galaxy properties before accretion, by simulating twenty seven haloes with masses between $5\times 10^8$ and $10^{10} M_\odot$. Out of this set nineteen haloes successfully form stars, while eight remain dark. The simulated galaxies match quite well present day observed scaling relations between stellar mass, size and metallicity, showing that such relations are in place before accretion. Our galaxies show a large variety of star formation histories, from extended star formation periods to single bursts. As in more massive galaxies, large star formation bursts are connected with major mergers events, which greatly contribute to the overall stellar mass build up. The intrinsic stochasticity of mergers induces a large scatter in the stellar mass halo mass relation, up to two orders of magnitude. Despite the bursty star formation history, on these mass scales baryons are very ineffective in modifying the dark matter profiles, and galaxies with a stellar mass below $\approx 10^6 M_\odot$ retain their cuspy central dark matter distribution, very similar to results from pure N-body simulations.

Imaging Atmospheric Cherenkov Telescopes (IACTs) need imaging optics with large apertures to map the faint Cherenkov light emitted in extensive air showers onto their image sensors. Segmented reflectors fulfill these needs using mass produced and light weight mirror facets. However, as the overall image is the sum of the individual mirror facet images, alignment is important. Here we present a method to determine the mirror facet positions on a segmented reflector in a very direct way. Our method reconstructs the mirror facet positions from photographs and a laser distance meter measurement which goes from the center of the image sensor plane to the center of each mirror facet. We use our method to both align the mirror facet positions and to feed the measured positions into our IACT simulation. We demonstrate our implementation on the 4 m First Geiger-mode Avalanche Cherenkov Telescope (FACT).

For studying variable sources like blazars, it is crucial to achieve unbiased monitoring, either with dedicated telescopes in pointing mode or survey instruments. At TeV energies, the High Altitude Water Cherenkov (HAWC) observatory monitors approximately two thirds of the sky every day. It uses the water Cherenkov technique, which provides an excellent duty cycle independent of weather and season. The First G-APD Cherenkov Telescope (FACT) monitors a small sample of sources with better sensitivity, using the imaging air Cherenkov technique. Thanks to its camera with silicon-based photosensors, FACT features an excellent detector performance and stability and extends its observations to times with strong moonlight, increasing the duty cycle compared to other imaging air Cherenkov telescopes. As FACT and HAWC have overlapping energy ranges, a joint study can exploit the longer daily coverage given that the observatories' locations are offset by 5.3 hours. Furthermore, the better sensitivity of FACT adds a finer resolution of features on hour-long time scales, while the continuous duty cycle of HAWC ensures evenly sampled long-term coverage. Thus, the two instruments complement each other to provide a more complete picture of blazar variability. In this presentation, the first joint study of light curves from the two instruments will be shown, correlating long-term measurements with daily sampling between air and water Cherenkov telescopes. The presented results focus on the study of the variability of the bright blazars Mrk 421 and Mrk 501 during the last two years featuring various flaring activities.

There is strong observational evidence indicating a time lag of order of some 100 Myr between the onset of starburst and AGN activity in galaxies. Dynamical time lags have been invoked to explain this. We extend this approach by introducing a viscous time lag the gas additionally needs to flow through the AGN's accretion disc before it reaches the central black hole. Our calculations reproduce the observed time lags and are in accordance with the observed correlation between black hole mass and stellar velocity dispersion.

Imaging Atmospheric Cherenkov Telescopes (IACTs) need imaging optics with large apertures and high image intensities to map the faint Cherenkov light emitted from cosmic ray air showers onto their image sensors. Segmented reflectors fulfill these needs, and as they are composed from mass production mirror facets they are inexpensive and lightweight. However, as the overall image is a superposition of the individual facet images, alignment is a challenge. Here we present a computer vision based star tracking alignment method, which also works for limited or changing star light visibility. Our method normalizes the mirror facet reflection intensities to become independent of the reference star's intensity or the cloud coverage. Using two CCD cameras, our method records the mirror facet orientations asynchronously of the telescope drive system, and thus makes the method easy to integrate into existing telescopes. It can be combined with remote facet actuation, but does not require one to work. Furthermore, it can reconstruct all individual mirror facet point spread functions without moving any mirror. We present alignment results on the 4 meter First Geiger-mode Avalanche Cherenkov Telescope (FACT).

Imaging Atmospheric Cherenkov Telescopes (IACTs) need imaging optics with large apertures and high image intensities to map the faint Cherenkov light emitted from cosmic ray air showers onto their image sensors. Segmented reflectors fulfill these needs, and composed from mass production mirror facets they are inexpensive and lightweight. However, as the overall image is a superposition of the individual facet images, alignment remains a challenge. Here we present a simple, yet extendable method, to align a segmented reflector using its Bokeh. Bokeh alignment does not need a star or good weather nights but can be done even during daytime. Bokeh alignment optimizes the facet orientations by comparing the segmented reflectors Bokeh to a predefined template. The optimal Bokeh template is highly constricted by the reflector's aperture and is easy accessible. The Bokeh is observed using the out of focus image of a near by point like light source in a distance of about 10 focal lengths. We introduce Bokeh alignment on segmented reflectors and demonstrate it on the First Geiger-mode Avalanche Cherenkov Telescope (FACT) on La Palma, Spain.

The circumnuclear disc (CND) orbiting the Galaxy's central black hole is a reservoir of material that can ultimately provide energy through accretion, or form stars in the presence of the black hole, as evidenced by the stellar cluster that is presently located at the CND's centre. In this paper, we report the results of a computational study of the dynamics of the CND. The results lead us to question two paradigms that are prevalent in previous research on the Galactic Centre. The first is that the disc's inner cavity is maintained by the interaction of the central stellar cluster's strong winds with the disc's inner rim, and second, that the presence of unstable clumps in the disc implies that the CND is a transient feature. Our simulations show that, in the absence of a magnetic field, the interaction of the wind with the inner disc rim actually leads to a filling of the inner cavity within a few orbital time-scales, contrary to previous expectations. However, including the effects of magnetic fields stabilizes the inner disc rim against rapid inward migration. Furthermore, this interaction causes instabilities that continuously create clumps that are individually unstable against tidal shearing. Thus the occurrence of such unstable clumps does not necessarily mean that the disc is itself a transient phenomenon. The next steps in this investigation are to explore the effect of the magnetorotational instability on the disc evolution and to test whether the results presented here persist for longer time-scales than those considered here.

We derive analytical expressions for the mass flow rates in and from accretion discs, taking into account the Eddington limit. This allows us to connect the basic properties of outflows with those of the accretion flow. As an example, we derive the radial surface density distribution in the accretion disc of Mrk 231.