-
Aims: To investigate the extension of the very-high-energy spectral tail of
the Crab pulsar at energies above 400 GeV. Methods: We analyzed $\sim$320 hours
of good quality data of Crab with the MAGIC telescope, obtained from February
2007 until April 2014. Results: We report the most energetic pulsed emission
ever detected from the Crab pulsar reaching up to 1.5 TeV. The pulse profile
shows two narrow peaks synchronized with the ones measured in the GeV energy
range. The spectra of the two peaks follow two different power-law functions
from 70 GeV up to 1.5 TeV and connect smoothly with the spectra measured above
10 GeV by the Large Area Telescope (LAT) on board of the Fermi satellite. When
making a joint fit of the LAT and MAGIC data, above 10 GeV, the photon indices
of the spectra differ by 0.5$\pm$0.1. Conclusions: We measured with the MAGIC
telescopes the most energetic pulsed photons from a pulsar to date. Such TeV
pulsed photons require a parent population of electrons with a Lorentz factor
of at least $5\times 10^6$. These results strongly suggest IC scattering off
low energy photons as the emission mechanism and a gamma-ray production region
in the vicinity of the light cylinder.
-
We present the first results from a new, high resolution, $^{12}$CO(1-0),
$^{13}$CO(1-0), and C$^{18}$O(1-0) molecular line survey of the Orion A cloud,
hereafter referred to as the CARMA-NRO Orion Survey. CARMA observations have
been combined with single-dish data from the Nobeyama 45m telescope to provide
extended images at about 0.01 pc resolution, with a dynamic range of
approximately 1200 in spatial scale. Here we describe the practical details of
the data combination in uv space, including flux scale matching, the conversion
of single dish data to visibilities, and joint deconvolution of single dish and
interferometric data. A $\Delta$-variance analysis indicates that no artifacts
are caused by combining data from the two instruments. Initial analysis of the
data cubes, including moment maps, average spectra, channel maps,
position-velocity diagrams, excitation temperature, column density, and line
ratio maps provides evidence of complex and interesting structures such as
filaments, bipolar outflows, shells, bubbles, and photo-eroded pillars. The
implications for star formation processes are profound and follow-up scientific
studies by the CARMA-NRO Orion team are now underway. We plan to make all the
data products described here generally accessible; some are already available
at https://dataverse.harvard.edu/dataverse/CARMA-NRO-Orion
-
The spherical Jeans equation (SJE) is widely used in dynamical modelling of
the Milky Way (MW) halo potential. We use haloes and galaxies from the
cosmological Millennium-II simulation and hydrodynamical APOSTLE simulations to
investigate the performance of the SJE in recovering the underlying mass
profiles of MW mass haloes. The best-fitting halo mass and concentration
parameters scatter by 25% and 40% around their input values, respectively, when
dark matter particles are used as tracers. This scatter becomes as large as a
factor of 3 when using star particles instead. This is significantly larger
than the estimated statistical uncertainty associated with the use of the SJE.
The existence of correlated phase-space structures that violate the steady
state assumption of the SJE as well as non-spherical geometries are the
principal sources of the scatter. Binary haloes show larger scatter because
they are more aspherical in shape and have a more perturbed dynamical state.
Our results confirm the previous study of Wang et al. (2017) that the number of
independent phase-space structures sets an intrinsic limiting precision on
dynamical inferences based on the steady state assumption. Modelling with a
radius-independent velocity anisotropy, or using tracers within a limited outer
radius, result in significantly larger scatter, but the ensemble-averaged
measurement over the whole halo sample is approximately unbiased.
-
We present a new complete Near-Infrared (NIR, $JHK_s$) census of RR Lyrae
stars (RRLs) in the globular $\omega$ Cen (NGC 5139). We collected 15,472
$JHK_s$ images with 4-8m class telescopes over 15 years (2000-2015) covering a
sky area around the cluster center of 60x34 arcmin$^2$. These images provided
calibrated photometry for 182 out of the 198 cluster RRL candidates with ten to
sixty measurements per band. We also provide new homogeneous estimates of the
photometric amplitude for 180 ($J$), 176 ($H$) and 174 ($K_s$) RRLs. These data
were supplemented with single-epoch $JK_s$ magnitudes from VHS and with
single-epoch $H$ magnitudes from 2MASS. Using proprietary optical and NIR data
together with new optical light curves (ASAS-SN) we also updated pulsation
periods for 59 candidate RRLs. As a whole, we provide $JHK_s$ magnitudes for 90
RRab (fundamentals), 103 RRc (first overtones) and one RRd (mixed--mode
pulsator). We found that NIR/optical photometric amplitude ratios increase when
moving from first overtone to fundamental and to long-period (P>0.7 days)
fundamental RRLs. Using predicted Period-Luminosity-Metallicity relations, we
derive a true distance modulus of 13.674$\pm$0.008$\pm$0.038 mag (statistical
error and standard deviation of the median)---based on spectroscopic iron
abundances---and of 13.698$\pm$0.004$\pm$0.048 mag---based on photometric iron
abundances. We also found evidence of possible systematics at the 5-10% level
in the zero-point of the PLs based on the five calibrating RRLs whose
parallaxes had been determined with HST
-
We present the latest results of a semi-analytic model of galaxy formation,
"New Numerical Galaxy Catalogue", which is combined with large cosmological
N-body simulations. This model can reproduce statistical properties of galaxies
at z < 6.0. We focus on the properties of active galactic nuclei (AGNs) and
supermassive black holes, especially on the accretion timescale onto black
holes. We find that the number density of AGNs at z < 1.5 and at hard X-ray
luminosity 10^{ 44 }< erg/s is underestimated compared with recent
observational estimates when we assume the exponentially decreasing accretion
rate and the accretion timescale which is proportional to the dynamical time of
the host halo or the bulge, as is often assumed in semi-analytic models. We
show that to solve this discrepancy, the accretion timescale of such less
luminous AGNs instead should be a function of the black hole mass and the
accreted gas mass. This timescale can be obtained from a phenomenological
modelling of the gas angular momentum loss in the circumnuclear torus and/or
the accretion disc. Such models predict a longer accretion timescale for less
luminous AGNs at z < 1.0 than bright QSOs whose accretion timescale would be
10^{ 7-8 } yr. With this newly introduced accretion timescale, our model can
explain the observed luminosity functions of AGNs at z < 6.0.
-
The energy of the ultrahigh energy spectral cutoff was measured, integrating
over the northern hemisphere sky, by the Telescope Array (TA) collaboration, to
be $10^{19.78 \pm 0.06}$ eV, in agreement with the High Resolution Fly's Eye
(HiRes) experiment, whereas the Pierre Auger experiment, integrating over the
southern hemisphere sky, measured the cutoff to be at 10$^{19.62 \pm 0.02}$ eV.
An 11% energy scale difference between the TA and Auger does not account for
this difference. However, in comparing the spectra of the Telescope Array and
Pierre Auger experiments in the band of declination common to both experiments
($-15.7^{\circ} < \delta < 24.8^{\circ}$) we have found agreement in the energy
of the spectral cutoff. While the Auger result is essentially unchanged, the TA
cutoff energy has changed to $10^{19.59 \pm 0.06}$ eV. In this paper we argue
that this is an astrophysical effect.
-
The flat-spectrum radio quasar PKS 1441+25 at a redshift of z = 0.940 is
detected between 40 and 250 GeV with a significance of 25.5 {\sigma} using the
MAGIC telescopes. Together with the gravitationally lensed blazar QSO B0218+357
(z = 0.944), PKS 1441+25 is the most distant very high energy (VHE) blazar
detected to date. The observations were triggered by an outburst in 2015 April
seen at GeV energies with the Large Area Telescope on board Fermi.
Multi-wavelength observations suggest a subdivision of the high state into two
distinct flux states. In the band covered by MAGIC, the variability time scale
is estimated to be 6.4 +/- 1.9 days. Modeling the broadband spectral energy
distribution with an external Compton model, the location of the emitting
region is understood as originating in the jet outside the broad line region
(BLR) during the period of high activity, while being partially within the BLR
during the period of low (typical) activity. The observed VHE spectrum during
the highest activity is used to probe the extragalactic background light at an
unprecedented distance scale for ground-based gamma-ray astronomy.
-
We developed control and visualization programs, YUI and HANA, for High-
Resolution Chopper spectrometer (HRC) installed at BL12 in MLF, J-PARC. YUI is
a comprehensive program to control DAQ-middleware, the accessories, and sample
environment devices. HANA is a program for the data transformation and
visualization of inelastic neutron scattering spectra. In this paper, we
describe the basic system structures and unique functions of these programs
from the viewpoint of users.
-
We present the cross-correlation between 151 luminous quasars ($M_{
\mathrm{UV}} < -26$) and 179 protocluster candidates at $z \sim 3.8$, extracted
from the Wide imaging survey ($ \sim 121~ $deg$^2$) performed with a part of
the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP). We find that only two
out of 151 quasars reside in regions that are more overdense compared to the
average field at $ > 4 \sigma $. The distributions of the distance between
quasars and the nearest protoclusters and the significance of the overdensity
at the position of quasars are statistically identical to those found for
$g$-dropout galaxies, suggesting that quasars tend to reside in almost the same
environment as star-forming galaxies at this redshift. Using stacking analysis,
we find that the average density of $g$-dropout galaxies around quasars is
slightly higher than that around $g$-dropout galaxies on $1.0 - 2.5$ pMpc
scales, while at $ < 0.5$ pMpc that around quasars tends to be lower. We also
find that quasars with higher UV-luminosity or with more massive black holes
tend to avoid the most overdense regions, and that the quasar near zone sizes
are anti-correlated with overdensity. These findings are consistent with a
scenario in which the luminous quasar at $z \sim4 $ resides in structures that
are less massive than those expected for the progenitors of today's rich
clusters of galaxies, and possibly that luminous quasars may be suppressing
star formation in their close vicinity.
-
We present the cross-correlation between 151 luminous quasars ($M_{
\mathrm{UV}} < -26$) and 179 protocluster candidates at $z \sim 3.8$, extracted
from the Wide imaging survey ($ \sim 121~ $deg$^2$) performed with a part of
the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP). We find that only two
out of 151 quasars reside in regions that are more overdense compared to the
average field at $ > 4 \sigma $. The distributions of the distance between
quasars and the nearest protoclusters and the significance of the overdensity
at the position of quasars are statistically identical to those found for
$g$-dropout galaxies, suggesting that quasars tend to reside in almost the same
environment as star-forming galaxies at this redshift. Using stacking analysis,
we find that the average density of $g$-dropout galaxies around quasars is
slightly higher than that around $g$-dropout galaxies on $1.0 - 2.5$ pMpc
scales, while at $ < 0.5$ pMpc that around quasars tends to be lower. We also
find that quasars with higher UV-luminosity or with more massive black holes
tend to avoid the most overdense regions, and that the quasar near zone sizes
are anti-correlated with overdensity. These findings are consistent with a
scenario in which the luminous quasar at $z \sim4 $ resides in structures that
are less massive than those expected for the progenitors of today's rich
clusters of galaxies, and possibly that luminous quasars may be suppressing
star formation in their close vicinity.
-
Spontaneous breaking of Lorentz symmetry at energies on the order of the
Planck energy or lower is predicted by many quantum gravity theories, implying
non-trivial dispersion relations for the photon in vacuum. Consequently,
gamma-rays of different energies, emitted simultaneously from astrophysical
sources, could accumulate measurable differences in their time of flight until
they reach the Earth. Such tests have been carried out in the past using fast
variations of gamma-ray flux from pulsars, and more recently from active
galactic nuclei and gamma-ray bursts. We present new constraints studying the
gamma-ray emission of the galactic Crab Pulsar, recently observed up to TeV
energies by the MAGIC collaboration. A profile likelihood analysis of pulsar
events reconstructed for energies above 400GeV finds no significant variation
in arrival time as their energy increases. Ninety-five percent~CL limits are
obtained on the effective Lorentz invariance violating energy scale at the
level of $E_{\mathrm{QG}_1} > 5.5\cdot 10^{17}$GeV ($4.5\cdot 10^{17}$GeV) for
a linear, and $E_{\mathrm{QG}_2} > 5.9\cdot 10^{10}$GeV ($5.3\cdot 10^{10}$GeV)
for a quadratic scenario, for the subluminal and the superluminal cases,
respectively. A substantial part of this study is dedicated to calibration of
the test statistic, with respect to bias and coverage properties. Moreover, the
limits take into account systematic uncertainties, found to worsen the
statistical limits by about 36--42\%. Our constraints would have resulted much
more competitive if the intrinsic pulse shape of the pulsar between 200GeV and
400GeV was understood in sufficient detail and allowed inclusion of events well
below 400GeV.
-
The microquasar Cygnus X-1 displays the two typical soft and hard X-ray
states of a black-hole transient. During the latter, Cygnus X-1 shows a
one-sided relativistic radio-jet. Recent detection of the system in the high
energy (HE; $E\gtrsim60$ MeV) gamma-ray range with \textit{Fermi}-LAT
associates this emission with the outflow. Former MAGIC observations revealed a
hint of flaring activity in the very high-energy (VHE; $E\gtrsim100$ GeV)
regime during this X-ray state. We analyze $\sim97$ hr of Cygnus X-1 data taken
with the MAGIC telescopes between July 2007 and October 2014. To shed light on
the correlation between hard X-ray and VHE gamma rays as previously suggested,
we study each main X-ray state separately. We perform an orbital phase-folded
analysis to look for variability in the VHE band. Additionally, to place this
variability behavior in a multiwavelength context, we compare our results with
\textit{Fermi}-LAT, \textit{AGILE}, \textit{Swift}-BAT, \textit{MAXI},
\textit{RXTE}-ASM, AMI and RATAN-600 data. We do not detect Cygnus X-1 in the
VHE regime. We establish upper limits for each X-ray state, assuming a
power-law distribution with photon index $\Gamma=3.2$. For steady emission in
the hard and soft X-ray states, we set integral upper limits at 95\% confidence
level for energies above 200 GeV at $2.6\times10^{-12}$~photons
cm$^{-2}$s$^{-1}$ and $1.0\times10^{-11}$~photons cm$^{-2}$s$^{-1}$,
respectively. We rule out steady VHE gamma-ray emission above this energy
range, at the level of the MAGIC sensitivity, originating in the interaction
between the relativistic jet and the surrounding medium, while the emission
above this flux level produced inside the binary still remains a valid
possibility.
-
The large jet kinetic power and non-thermal processes occurring in the
microquasar SS 433 make this source a good candidate for a very high-energy
(VHE) gamma-ray emitter. Gamma-ray fluxes have been predicted for both the
central binary and the interaction regions between jets and surrounding nebula.
Also, non-thermal emission at lower energies has been previously reported. We
explore the capability of SS 433 to emit VHE gamma rays during periods in which
the expected flux attenuation due to periodic eclipses and precession of the
circumstellar disk periodically covering the central binary system is expected
to be at its minimum. The eastern and western SS433/W50 interaction regions are
also examined. We aim to constrain some theoretical models previously developed
for this system. We made use of dedicated observations from MAGIC and H.E.S.S.
from 2006 to 2011 which were combined for the first time and accounted for a
total effective observation time of 16.5 h. Gamma-ray attenuation does not
affect the jet/medium interaction regions. The analysis of a larger data set
amounting to 40-80 h, depending on the region, was employed. No evidence of VHE
gamma-ray emission was found. Upper limits were computed for the combined data
set. We place constraints on the particle acceleration fraction at the inner
jet regions and on the physics of the jet/medium interactions. Our findings
suggest that the fraction of the jet kinetic power transferred to relativistic
protons must be relatively small to explain the lack of TeV and neutrino
emission from the central system. At the SS433/W50 interface, the presence of
magnetic fields greater 10$\mu$G is derived assuming a synchrotron origin for
the observed X-ray emission. This also implies the presence of high-energy
electrons with energies up to 50 TeV, preventing an efficient production of
gamma-ray fluxes in these interaction regions.
-
Context. PKS 1510-089 is one of only a few flat spectrum radio quasars
detected in the VHE (very-high-energy, > 100 GeV) gamma-ray band. Aims. We
study the broadband spectral and temporal properties of the PKS 1510-089
emission during a high gamma-ray state. Methods. We performed VHE gamma-ray
observations of PKS 1510-089 with the MAGIC telescopes during a long high
gamma-ray state in May 2015. In order to perform broadband modelling of the
source, we have also gathered contemporaneous multiwavelength data in radio,
IR, optical photometry and polarization, UV, X-ray and GeV gamma-ray ranges. We
construct a broadband spectral energy distribution (SED) in two periods,
selected according to VHE gamma-ray state. Results. PKS 1510-089 has been
detected by MAGIC during a few day-long observations performed in the middle of
a long, high optical and gamma-ray state, showing for the first time a
significant VHE gamma-ray variability. Similarly to the optical and gamma-ray
high state of the source detected in 2012, it was accompanied by a rotation of
the optical polarization angle and the emission of a new jet component observed
in radio. However, due to large uncertainty on the knot separation time, the
association with the VHE gamma-ray emission cannot be firmly established. The
spectral shape in the VHE band during the flare is similar to the ones obtained
during previous measurements of the source. The observed flux variability sets
for the first time constraints on the size of the region from which VHE gamma
rays are emitted. We model the broadband SED in the framework of the external
Compton scenario and discuss the possible emission site in view of
multiwavelength data and alternative emission models.
-
In this work we present data from observations with the MAGIC telescopes of
SN 2014J detected in January 21 2014, the closest Type Ia supernova since
Imaging Air Cherenkov Telescopes started to operate. We probe the possibility
of very-high-energy (VHE; $E\geq100$ GeV) gamma rays produced in the early
stages of Type Ia supernova explosions. We performed follow-up observations
after this supernova explosion for 5 days, between January 27 and February 2 in
2014. We search for gamma-ray signal in the energy range between 100 GeV and
several TeV from the location of SN 2014J using data from a total of $\sim5.5$
hours of observations. Prospects for observing gamma-rays of hadronic origin
from SN 2014J in the near future are also being addressed. No significant
excess was detected from the direction of SN 2014J. Upper limits at 95$\%$
confidence level on the integral flux, assuming a power-law spectrum,
d$F/$d$E\propto E^{-\Gamma}$, with a spectral index of $\Gamma=2.6$, for
energies higher than 300 GeV and 700 GeV, are established at
$1.3\times10^{-12}$ and $4.1\times10^{-13}$ photons~cm$^{-2}$s$^{-1}$,
respectively. For the first time, upper limits on the VHE emission of a Type Ia
supernova are established. The energy fraction isotropically emitted into TeV
gamma rays during the first $\sim10$ days after the supernova explosion for
energies greater than 300 GeV is limited to $10^{-6}$ of the total available
energy budget ($\sim 10^{51}$ erg). Within the assumed theoretical scenario,
the MAGIC upper limits on the VHE emission suggest that SN 2014J will not be
detectable in the future by any current or planned generation of Imaging
Atmospheric Cherenkov Telescopes.
-
We present the first detection of the nearby (z=0.084) low-luminosity BL Lac
object 1ES 1741+196 in the very high energy (VHE: E$>$100 GeV) band. This
object lies in a triplet of interacting galaxies. Early predictions had
suggested 1ES 1741+196 to be, along with several other high-frequency BL Lac
sources, within the reach of MAGIC detectability. Its detection by MAGIC, later
confirmed by VERITAS, helps to expand the small population of known TeV BL
Lacs. The source was observed with the MAGIC telescopes between 2010 April and
2011 May, collecting 46 h of good quality data. These observations led to the
detection of the source at 6.0 $\sigma$ confidence level, with a steady flux
$\mathrm{F}(> 100 {\rm GeV}) = (6.4 \pm 1.7_{\mathrm{stat}}\pm
2.6_{\mathrm{syst}}) \cdot 10^{-12}$ ph cm$^{-2}$ s$^{-1}$ and a differential
spectral photon index $\Gamma = 2.4 \pm 0.2_{\mathrm{stat}} \pm
0.2_{\mathrm{syst}}$ in the range of $\sim$80 GeV - 3 TeV. To study the
broad-band spectral energy distribution (SED) simultaneous with MAGIC
observations, we use KVA, Swift/UVOT and XRT, and Fermi/LAT data. One-zone
synchrotron-self-Compton (SSC) modeling of the SED of 1ES 1741+196 suggests
values for the SSC parameters that are quite common among known TeV BL Lacs
except for a relatively low Doppler factor and slope of electron energy
distribution. A thermal feature seen in the SED is well matched by a giant
elliptical's template. This appears to be the signature of thermal emission
from the host galaxy, which is clearly resolved in optical observations.
-
We present an extensive study of the BL Lac object Mrk 501 based on a data
set collected during the multi-instrument campaign spanning from 2009 March 15
to 2009 August 1 which includes, among other instruments, MAGIC, VERITAS,
Whipple 10-m, Fermi-LAT, RXTE, Swift, GASP-WEBT and VLBA. We find an increase
in the fractional variability with energy, while no significant interband
correlations of flux changes are found in the acquired data set. The higher
variability in the very high energy (>100 GeV, VHE) gamma-ray emission and the
lack of correlation with the X-ray emission indicate that the highest-energy
electrons that are responsible for the VHE gamma-rays do not make a dominant
contribution to the ~1 keV emission. Alternatively, there could be a very
variable component contributing to the VHE gamma-ray emission in addition to
that coming from the synchrotron self-Compton (SSC) scenarios. The space of SSC
model parameters is probed following a dedicated grid-scan strategy, allowing
for a wide range of models to be tested and offering a study of the degeneracy
of model-to-data agreement in the individual model parameters. We find that
there is some degeneracy in both the one-zone and the two-zone SSC scenarios
that were probed, with several combinations of model parameters yielding a
similar model-to-data agreement, and some parameters better constrained than
others. The SSC model grid-scan shows that the flaring activity around 2009 May
22 cannot be modeled adequately with a one-zone SSC scenario, while it can be
suitably described within a two-independent-zone SSC scenario. The observation
of an electric vector polarization angle rotation coincident with the gamma-ray
flare from 2009 May 1 resembles those reported previously for low frequency
peaked blazars, hence suggesting that there are many similarities in the
flaring mechanisms of blazars with different jet properties.
-
We present optical observations of SN 2013dx, related to the Fermi burst GRB
130702A occurred at a redshift z = 0.145. It is the second-best sampled GRB-SN
after SN~1998bw: the observational light curves contain more than 280 data
points in uBgrRiz filters until 88 day after the burst, and the data were
collected from our observational collaboration (Maidanak Observatory,
Abastumani Observatory, Crimean Astrophysical Observatory, Mondy Observatory,
National Observatory of Turkey, Observatorio del Roque de los Muchachos) and
from the literature. We model numerically the multicolour light curves using
the one-dimensional radiation hydrodynamical code STELLA, previously widely
implemented for the modelling of typical non-GRB SNe. The best-fitted model has
the following parameters: pre-supernova star mass M = 25 M_Sun, mass of a
compact remnant M_CR = 6 M_Sun, total energy of the outburst E_oburst = 3.5 x
10^(52) erg, pre-supernova star radius R = 100 R_Sun, M_56Ni = 0.2 M_Sun which
is totally mixed through the ejecta; M_O = 16.6 M_Sun, M_Si = 1.2 M_Sun, and
M_Fe = 1.2 M_Sun, and the radiative efficiency of the SN is 0.1 per cent.
-
The MAGIC (Major Atmospheric Gamma-ray Imaging Cherenkov) telescopes observed
the BL Lac object H1722+119 (redshift unknown) for six consecutive nights
between 2013 May 17 and 22, for a total of 12.5 h. The observations were
triggered by high activity in the optical band measured by the KVA (Kungliga
Vetenskapsakademien) telescope. The source was for the first time detected in
the very high energy (VHE, $E > 100$ GeV) $\gamma$-ray band with a statistical
significance of 5.9 $\sigma$. The integral flux above 150 GeV is estimated to
be $(2.0\pm 0.5)$ per cent of the Crab Nebula flux. We used contemporaneous
high energy (HE, 100 MeV $ < E < 100$ GeV) $\gamma$-ray observations from
Fermi-LAT (Large Area Telescope) to estimate the redshift of the source. Within
the framework of the current extragalactic background light models, we estimate
the redshift to be $z = 0.34 \pm 0.15$. Additionally, we used contemporaneous
X-ray to radio data collected by the instruments on board the Swift satellite,
the KVA, and the OVRO (Owens Valley Radio Observatory) telescope to study
multifrequency characteristics of the source. We found no significant temporal
variability of the flux in the HE and VHE bands. The flux in the optical and
radio wavebands, on the other hand, did vary with different patterns. The
spectral energy distribution (SED) of H1722+119 shows surprising behaviour in
the $\sim 3\times10^{14} - 10^{18}$ Hz frequency range. It can be modelled
using an inhomogeneous helical jet synchrotron self-Compton model.
-
Context. QSO B0218+357 is a gravitationally lensed blazar located at a
redshift of 0.944. The gravitational lensing splits the emitted radiation into
two components, spatially indistinguishable by gamma-ray instruments, but
separated by a 10-12 day delay. In July 2014, QSO B0218+357 experienced a
violent flare observed by the Fermi-LAT and followed by the MAGIC telescopes.
Aims. The spectral energy distribution of QSO B0218+357 can give information on
the energetics of z ~ 1 very high energy gamma- ray sources. Moreover the
gamma-ray emission can also be used as a probe of the extragalactic background
light at z ~ 1. Methods. MAGIC performed observations of QSO B0218+357 during
the expected arrival time of the delayed component of the emission. The MAGIC
and Fermi-LAT observations were accompanied by quasi-simultaneous optical data
from the KVA telescope and X-ray observations by Swift-XRT. We construct a
multiwavelength spectral energy distribution of QSO B0218+357 and use it to
model the source. The GeV and sub-TeV data, obtained by Fermi-LAT and MAGIC,
are used to set constraints on the extragalactic background light. Results.
Very high energy gamma-ray emission was detected from the direction of QSO
B0218+357 by the MAGIC telescopes during the expected time of arrival of the
trailing component of the flare, making it the farthest very high energy
gamma-ray sources detected to date. The observed emission spans the energy
range from 65 to 175 GeV. The combined MAGIC and Fermi-LAT spectral energy
distribution of QSO B0218+357 is consistent with current extragalactic
background light models. The broad band emission can be modeled in the
framework of a two zone external Compton scenario, where the GeV emission comes
from an emission region in the jet, located outside the broad line region.
-
The stellar initial mass function (IMF) of early-type galaxies is the
combination of the IMF of the stellar population formed in-situ and that of
accreted stellar populations. Using as an observable the effective IMF
$\alpha_{IMF}$, defined as the ratio between the true stellar mass of a galaxy
and the stellar mass inferred assuming a Salpeter IMF, we present a theoretical
model for its evolution as a result of dry mergers. We use a simple dry merger
evolution model, based on cosmological $N$-body simulations, together with
empirically motivated prescriptions for the IMF to make predictions for how the
effective IMF of massive early-type galaxies changes from $z=2$ to $z=0$. We
find that the IMF normalization of individual galaxies becomes lighter with
time. At fixed velocity dispersion, $\alpha_{IMF}$ is predicted to be constant
with redshift. Current constraints on the evolution of the IMF are in slight
tension with this prediction, even though systematic uncertainties prevent a
conclusive statement. The correlation of $\alpha_{IMF}$ with stellar mass
becomes shallower with time, while the correlation between $\alpha_{IMF}$ and
velocity dispersion is mostly preserved by dry mergers. We also find that dry
mergers can mix the dependence of the IMF on stellar mass and velocity
dispersion, making it challenging to infer, from $z=0$ observations of global
galactic properties, what is the quantity that is originally coupled with the
IMF.
-
Upgraded electronics, improved water system dynamics, better calibration and
analysis techniques allowed Super-Kamiokande-IV to clearly observe very
low-energy 8B solar neutrino interactions, with recoil electron kinetic
energies as low as 3.49 MeV. Super-Kamiokande-IV data-taking began in September
of 2008; this paper includes data until February 2014, a total livetime of 1664
days. The measured solar neutrino flux is (2.308+-0.020(stat.) +
0.039-0.040(syst.)) x 106/(cm2sec) assuming no oscillations. The observed
recoil electron energy spectrum is consistent with no distortions due to
neutrino oscillations. An extended maximum likelihood fit to the amplitude of
the expected solar zenith angle variation of the neutrino-electron elastic
scattering rate in SK-IV results in a day/night asymmetry of
(-3.6+-1.6(stat.)+-0.6(syst.))%. The SK-IV solar neutrino data determine the
solar mixing angle as sin2 theta_12 = 0.327+0.026-0.031, all SK solar data
(SK-I, SK-II, SK III and SKIV) measures this angle to be sin2 theta_12 =
0.334+0.027-0.023, the determined mass-squared splitting is Delta m2_21 =
4.8+1.5-0.8 x10-5 eV2.
-
We explore the effect of varying the mass of the seed black hole on the
resulting black hole mass - bulge mass relation at z ~ 0, using a semi-analytic
model of galaxy formation combined with large cosmological N-body simulations.
We constrain our model by requiring the observed properties of galaxies at z ~
0 are reproduced. In keeping with previous semi-analytic models, we place a
seed black hole immediately after a galaxy forms. When the mass of the seed is
set at 10^5 M_sun, we find that the model results become inconsistent with
recent observational results of the black hole mass - bulge mass relation for
dwarf galaxies. In particular, the model predicts that bulges with ~ 10^9 M_sun
harbour larger black holes than observed. On the other hand, when we employ
seed black holes with 10^3 M_sun, or randomly select their mass within a
10^(3-5) M_sun range, the resulting relation is consistent with observation
estimates, including the observed dispersion. We find that to obtain stronger
constraints on the mass of seed black holes, observations of less massive
bulges at z ~ 0 are a more powerful comparison than the relations at higher
redshifts.
-
We study the multi-band variability and correlations of the TeV blazar Mrk
421 on year time scales, which can bring additional insight on the processes
responsible for its broadband emission. We observed Mrk 421 in the very high
energy (VHE) gamma-ray range with the Cherenkov telescope MAGIC-I from March
2007 to June 2009 for a total of 96 hours of effective time after quality cuts.
The VHE flux variability is quantified with several methods, including the
Bayesian Block algorithm, which is applied to data from Cherenkov telescopes
for the first time. The 2.3 year long MAGIC light curve is complemented with
data from the Swift/BAT and RXTE/ASM satellites and the KVA, GASP-WEBT, OVRO,
and Mets\"ahovi telescopes from February 2007 to July 2009, allowing for an
excellent characterisation of the multi-band variability and correlations over
year time scales. Mrk 421 was found in different gamma-ray emission states
during the 2.3 year long observation period. Flares and different levels of
variability in the gamma-ray light curve could be identified with the Bayesian
Block algorithm. The same behaviour of a quiet and active emission was found in
the X-ray light curves measured by Swift/BAT and the RXTE/ASM, with a direct
correlation in time. The behaviour of the optical light curve of GASP-WEBT and
the radio light curves by OVRO and Mets\"ahovi are different as they show no
coincident features with the higher energetic light curves and a less variable
emission. The fractional variability is overall increasing with energy. The
comparable variability in the X-ray and VHE bands and their direct correlation
during both high- and low-activity periods spanning many months show that the
electron populations radiating the X-ray and gamma-ray photons are either the
same, as expected in the Synchrotron-Self-Compton mechanism, or at least
strongly correlated, as expected in electromagnetic cascades.
-
The gamma-ray binary LS I +61$^{\circ}$303 is a well established source from
centimeter radio up to very high energy (VHE; E$>$100 GeV). Its broadband
emission shows a periodicity of $\sim$26.5 days, coincident with the orbital
period. A longer (super-orbital) period of 1667 $\pm$ 8 days was discovered in
radio and confirmed in optical and high energy (HE; E>100 MeV) gamma-ray
observations. We present a four-year campaign performed by MAGIC together with
archival data concentrating on a search for a long timescale signature in the
VHE emission. We focus on the search for super-orbital modulation of the VHE
peak and on the search for correlations between TeV emission and optical
determination of the extension of the circumstellar disk. A four-year campaign
has been carried out by MAGIC. The source was observed during the orbital
phases when the periodic VHE outbursts have occurred ($\phi$=0.55-0.75).
Additionally, we included archival MAGIC observations and data published by the
VERITAS collaboration in these studies. For the correlation studies, LS I
+61$^{\circ}$303 has also been observed during the orbital phases where
sporadic VHE emission had been detected in the past ($\phi$=0.75-1.0). These
MAGIC observations were simultaneous with optical spectroscopy from the
LIVERPOOL telescope. The TeV flux of the periodical outburst in orbital phases
$\phi$=0.5--0.75 was found to show yearly variability consistent with the
$\sim$4.5 years long-term modulation found in the radio band. This modulation
of the TeV flux can be well described by a sine function with the best fit
period of $1610\pm 58$ days. The complete dataset span two super-orbital
periods. There is no evidence for a correlation between the TeV emission and
the mass-loss rate of the Be star but this may be affected by the strong, short
timescale (as short as intra-day) variation displayed by the H$\alpha$ fluxes.