Search

Using the First Light And Reionisation Epoch Simulations (\textsc{Flares}) we explore the dust driven properties of massive high-redshift galaxies at $z\in[5,10]$. By post-processing the galaxy sample using the radiative transfer code \textsc{skirt} we obtain the full spectral energy distribution. We explore the resultant luminosity functions, IRX-$\beta$ relations as well as the luminosity-weighted dust temperatures in the Epoch of Reionisation (EoR). We find that most of our results are in agreement with the current set of observations, but under-predict the number densities of bright IR galaxies, which are extremely biased towards the most overdense regions. We see that the \textsc{Flares} IRX-$\beta$ relation (for $5\le z\le8$) predominantly follows the local starburst relation. The IRX shows an increase with stellar mass, plateauing at the high-mass end ($\sim10^{10}$M$_{\odot}$) and shows no evolution in the median normalisation with redshift. We also look at the dependence of the peak dust temperature ($T_{\mathrm{peak}}$) on various galaxy properties including the stellar mass, IR luminosity and sSFR, finding the correlation to be strongest with sSFR. The luminosity-weighted dust temperatures increase towards higher redshifts, with the slope of the $T_{\mathrm{peak}}$ - redshift relation showing a higher slope than the lower redshift relations obtained from previous observational and theoretical works. The results from \textsc{Flares}, which is able to provide a better statistical sample of high-redshift galaxies compared to other simulations, provides a distinct vantage point for the high-redshift Universe., Comment: 19 pages, 16 figures, Published in MNRAS. Data link: https://flaresimulations.github.io/data.html

We focus our analysis on 55 BL Lac objects with a hard Fermi gamma-ray spectrum, and for which a redshift or a lower limit to it has been determined by a previous study of ours. We extrapolate the spectral fits given by the 4FGL catalogue to the VHE band (>0.1 TeV), which can be explored by imaging atmospheric Cherenkov telescopes. Furthermore, we take into account the absorption due to the extragalactic background light, strongly depending on the redshift. Finally, we compare our results with publicly available sensitivity curves for a selection of imaging atmospheric Cherenkov telescope arrays currently operating or under construction. From our extrapolations and simulations we find a large number of promising candidates for observation with the forthcoming Cherenkov Telescope Array observatory., 16 pages, 3 figures, accepted for publication in MNRAS

The diffuse extragalactic background light (EBL) is formed by ultraviolet (UV), optical, and infrared (IR) photons mainly produced by star formation processes over the history of the Universe, and contains essential information about galaxy evolution and cosmology. Here, we present a new determination of the evolving EBL spectral energy distribution using a novel approach purely based on galaxy data aiming to reduce current uncertainties on the higher redshifts and IR intensities. Our calculations use multiwavelength observations from the UV to the far-IR of a sample of approximately 150,000 galaxies detected up to $z\sim 6$ in the five fields of the Cosmic Assembly Near-Infrared Deep Extragalactic Legacy Survey (CANDELS) from the Hubble Space Telescope. This is one of the most comprehensive and deepest multi-wavelength galaxy datasets ever obtained. These unprecedented resources allow us to derive the overall EBL evolution up to $z\sim 6$ and its uncertainties. Our results agree with cosmic observables estimated from galaxy surveys and $\gamma$-ray attenuation such as monochromatic luminosity densities, including those in the far-IR, and star formation rate densities, also at the highest redshits. Optical depths from our EBL approximation, which will be robust at high redshifts and for $\gamma$ rays up to tens of TeV, will be reported in a companion paper., Comment: 17 pages, 13 figures, 4 tables; accepted for publication in MNRAS on August 6, 2021. Online material available at https://www.ucm.es/blazars/ebl

Synchrotron radiation from a decelerating blastwave is a widely accepted model of radio to X-ray afterglow emission from gamma-ray bursts (GRBs). GeV gamma-ray emission detected by the Fermi Large Area Telescope (LAT) and the duration of which extends beyond the prompt gamma-ray emission phase, is also compatible with broad features of afterglow emission. We revisit the synchrotron self-Compton (SSC) emission model from a decelerating blastwave to fit multiwavelength data from three bright GRBs, namely GRB~190114C, GRB~130427A and GRB~090510. We constrain the afterglow model parameters using the simultaneous fit of the spectral energy distributions at different times and light curves at different frequencies for these bursts. We find that a constant density interstellar medium is favored for the short GRB~090510, while a wind-type environment is favored for the long GRB~130427A and GRB~190114C. The sub-TeV component in GRB~190114C detected by MAGIC is the SSC emission in our modelling. Furthermore we find that the SSC emission in the Thomson regime is adequate to fit the spectra and light curves of GRB~190114C. For the other two GRBs, lacking sub-TeV detection, the SSC emissions are also modeled in the Thomson regime. For the model parameters we have used, the $\gamma\gamma$ attenuation in the blastwave is negligible in the sub-TeV range compared to the redshift-dependent $\gamma\gamma$ attenuation in the extragalactic background light., Comment: 14 pages, 4 figures, Accepted for publication in MNRAS, Comments welcome

The study of gamma-ray blazars is usually hindered due to the lack of information on their redshifts and on their low energy photon fields. This information is key to understand the effect on the gamma-ray absorption due to either extragalactic background light and/or intrinsic absorption and emission processes. All this information has also an impact on the determination of the location of the emitting region within the relativistic jets. In this work a new optical spectroscopic characterization is presented for three gamma-ray blazars: S4 0954+65, TXS 1515-273 and RX J0812.0+0237. For all the three targets the redshift determination is successful, and for the first time in the case of TXS 1515-273 and RX J0812.0+0237. Their classification as BL~Lac type is confirmed based on these new optical spectra. For S4 0954+65 (z=$0.3694\pm0.0011$) an estimation on the disk, broad line region and torus luminosities is performed based on the observed optical emission lines. The results from this study are compatible with the nature of S4 0954+65 as a transitional blazar. In the case of TXS 1515-273 ($z=0.1281\pm 0.0004$), although its optical spectrum is dominated by the continuum emission from the jet, applying the pPXF technique, the stellar population can be unveiled and is compatible with an old and metallic population. It is also the case of RX J0812.0+0237 ($z=0.1721\pm 0.0002$). Moreover, this work confirms that the optical spectrum from RX J0812.0+0237 is compatible with an extreme blazar classification., Submitted to MNRAS

The cosmic spectral energy distribution (CSED) is the total emissivity as a function of wavelength of galaxies in a given cosmic volume. We compare the observed CSED from the UV to the submm to that computed from the EAGLE cosmological hydrodynamical simulation, post-processed with stellar population synthesis models and including dust radiative transfer using the SKIRT code. The agreement with the data is better than 0.15 dex over the entire wavelength range at redshift $z=0$, except at UV wavelengths where the EAGLE model overestimates the observed CSED by up to a factor 2. Global properties of the CSED as inferred from CIGALE fits, such as the stellar mass density, mean star formation density, and mean dust-to-stellar-mass ratio, agree to within better than 20 per cent. At higher redshift, EAGLE increasingly underestimates the CSED at optical-NIR wavelengths with the FIR/submm emissivity underestimated by more than a factor of 5 by redshift $z=1$. We believe that these differences are due to a combination of incompleteness of the EAGLE-SKIRT database, the small simulation volume and the consequent lack of luminous galaxies, and our lack of knowledge on the evolution of the characteristics of the interstellar dust in galaxies. The impressive agreement between the simulated and observed CSED at lower $z$ confirms that the combination of EAGLE and SKIRT dust processing yields a fairly realistic representation of the local Universe., 14 pages, 5 figures, accepted for publication in MNRAS

Interactions between the extragalactic background light (EBL) and very high energy {\gamma} rays (VHE; E > 10 GeV) from cosmological sources alter their {\gamma}-ray spectrum. The stronger absorption of harder {\gamma} rays causes a steepening of the observed {\gamma}-ray spectrum at the high energy tail which can be expressed by an increase of the power index. The effect provides a link between high energy astrophysics and the evolution of galaxies. In this work we develop a new hybrid EBL model by augmenting our previous analytic model with information from semi-analytic galaxy catalogues. The model allows us to study the {\gamma}-ray opacities of individual {\gamma}-ray paths through the (simulated) universe and evaluate the effect of local fluctuations in the EBL intensity along each path. We confirm an order of magnitude fluctuations in the local EBL (based on the cumulative light from model galaxy distributions within spheres of 50 $h^{-1}$Mpc). However, the effect of these fluctuations on the VHE {\gamma}-ray opacity is insignificant due to the overall very small contribution of the local EBL to the total EBL. We also investigate the effect which a galaxy crossing the line of sight of the {\gamma}-ray source may have on the VHE spectrum. We find that galaxies with stellar masses of M$ > 10^{11}$ M$_{\odot}$ could have significant effect on the {\gamma}-ray absorption. It is unlikely that the observed variation in the spectral index is caused by the proximity of single galaxies to the {\gamma}-ray paths as these are extremely rare., Comment: 10 pages, 11 figures

We present the 1-100 GeV spectral energy distribution for a population of 148 high-synchrotron-peaked blazars (HSPs) recently detected with Fermi-LAT as part of the First Brazil-ICRANet Gamma-ray Blazar catalogue (1BIGB). Most of the 1BIGB sources do not appear in previous Fermi-LAT catalogues and their gamma-ray spectral properties are presented here for the first time, representing a significant new extension of the gamma-ray blazar population. Since our sample was originally selected from an excess signal in the 0.3-500 GeV band, the sources stand out as promising TeV blazar candidates, potentially in reach of the forthcoming very-high-energy (VHE) gamma-ray observatory, CTA. The flux estimates presented here are derived considering PASS8 data, integrating over more than 9 years of Fermi-LAT observations. We also review the full broadband fit between 0.3-500 GeV presented in the original 1BIGB paper for all sources, updating the power-law parameters with currently available Fermi-LAT dataset. The importance of these sources in the context of VHE population studies with both current instruments and the future CTA is evaluated. To do so, we select a subsample of 1BIGB sources and extrapolate their gamma-ray SEDs to the highest energies, properly accounting for absorption due to the extragalactic background light. We compare those extrapolations to the published CTA sensitivity curves and estimate their detectability by CTA. Two notable sources from our sample, namely 1BIGB J224910.6-130002 and 1BIGB J194356.2+211821, are discussed in greater detail. All gamma-ray SEDs, which are shown here for the first time, are made publicly available via the Brazilian Science Data Center (BSDC) service, maintained at CBPF, in Rio de Janeiro., Comment: 13 pages, 6 figures

This study investigates the impact of the fluctuations in the extra galactic background light (EBL) on the attenuation of the hard $\gamma$-ray spectra of distant blazars. EBL fluctuations occur on the scales up to 100 Mpc and are caused by clustering of galaxies. The EBL photons interact with high energy $\gamma$-rays via the electron-positron pair production mechanism: $\gamma + \gamma' \rightarrow e^+ + e^-$. The attenuation of $\gamma$-rays depends on their energy and the density of the intervening EBL photon field. Using a simple model for the evolution of the mean EBL photon density, we implement an analytical description of the EBL fluctuations. We find that the amplitudes of the EBL energy density can vary by $\pm 1\%$ as a function of environment. The EBL fluctuations lead to mild alterations of the optical depth or equivalently the transmissivity for $\gamma$-rays from distant blazars. Our model predicts maximum changes of $\pm 10\%$ in the $\gamma$-ray transmissivity. However, this translates into marginal differences in the power law slopes of currently observed $\gamma$-ray spectra. The slopes of deabsorbed $\gamma$-ray spectra differ by not more than $\pm 1\%$ if EBL fluctuations are included.

TeV photons from extragalactic sources are absorbed in the intergalactic medium and initiate electromagnetic cascades. These cascades offer a unique tool to probe the properties of the universe at cosmological scales. We present a new Monte Carlo code dedicated to the physics of such cascades. This code has been tested against both published results and analytical approximations, and is made publicly available. Using this numerical tool, we investigate the main cascade properties (spectrum, halo extension, time delays), and study in detail their dependence on the physical parameters (extra-galactic magnetic field, extra-galactic background light, source redshift, source spectrum and beaming emission). The limitations of analytical solutions are emphasised. In particular, analytical approximations account only for the first generation of photons and higher branches of the cascade tree are neglected., Comment: 18 pages, 21 figures, Accepted by MNRAS; this preprint matches accepted version

Monthly notices of the Royal Astronomical Society 465(3), 3767 - 3774 (2017). doi:10.1093/mnras/stw3032, The absorption of TeV $\gamma$-ray photons produced in relativistic jets by surrounding soft photon fields is a long-standing problem of jet physics. In some cases, the most likely emission site close to the central black hole is ruled out because of the high opacity caused by strong optical and infrared photon sources, such as the broad-line region. Mostly neglected for jet modelling is the absorption of $\gamma$-rays in the starlight photon field of the host galaxy. Analysing the absorption for arbitrary locations and observation angles of the $\gamma$-ray emission site within the host galaxy, we find that the distance to the galaxy centre, the observation angle, and the distribution of starlight in the galaxy are crucial for the amount of absorption. We derive the absorption value for a sample of 20 TeV-detected blazars with a redshift $z_r$ < 0.2. The absorption value of the $\gamma$-ray emission located in the galaxy centre may be as high as 20 per cent, with an average value of 6 per cent. This is important in order to determine the intrinsic blazar parameters. We see no significant trends in our sample between the degree of absorption and host properties, such as starlight emissivity, galactic size, half-light radius, and redshift. While the uncertainty of the spectral properties of the extragalactic background light exceeds the effect of absorption by stellar light from the host galaxy in distant objects, the latter is a dominant effect in nearby sources. It may also be revealed in a differential comparison of sources with similar redshifts., Published by Oxford Univ. Press, Oxford

We investigate the formation of TeV spectrum of a distant blazar PKS 1424+240 residing at a redshift $z\geq0.6$ in two scenarios in the frame of a lepto-hadronic jet model taking both the uncertainties of the extragalactic background light (EBL) and its redshift into account. In the first scenario, TeV emission is attributed to the synchrotron emission of pair cascades resulting from $p\gamma$ interaction; in the second scenario, TeV emission is attributed to the proton-synchrotron emission, and an internal absorption due to interaction with the photons around the jet is included. The results show that in the first scenario the 68\% upper limit of its redshift within which this scenario can explain the VERITAS TeV spectrum in 2009 well is $\sim0.75$; in the second scenario, this upper limit of the redshift becomes $\sim1.03$. However, the second scenario can be excluded because it requires an unreasonable photon field around the jet with a luminosity of $\sim10^{43}\rm \ erg\ s^{-1}$. In conclusion, the jet model can explain its TeV spectrum with a low EBL density if $0.6, Comment: accepted by MNRAS

We present very-high-energy $\gamma$-ray observations of the BL Lac object 1ES 2344+514 taken by the Very Energetic Radiation Imaging Telescope Array System (VERITAS) between 2007 and 2015. 1ES 2344+514 is detected with a statistical significance above background of $20.8\sigma$ in $47.2$ hours (livetime) of observations, making this the most comprehensive very-high-energy study of 1ES 2344+514 to date. Using these observations the temporal properties of 1ES 2344+514 are studied on short and long times scales. We fit a constant flux model to nightly- and seasonally-binned light curves and apply a fractional variability test, to determine the stability of the source on different timescales. We reject the constant-flux model for the 2007-2008 and 2014-2015 nightly-binned light curves and for the long-term seasonally-binned light curve at the $> 3\sigma$ level. The spectra of the time-averaged emission before and after correction for attenuation by the extragalactic background light are obtained. The observed time-averaged spectrum above 200 GeV is satisfactorily fitted (${\chi^2/NDF = 7.89/6}$) by a power-law function with index $\Gamma = 2.46 \pm 0.06_{stat} \pm 0.20_{sys} $ and extends to at least 8 TeV. The extragalactic-background-light-deabsorbed spectrum is adequately fit (${\chi^2/NDF = 6.73/6}$) by a power-law function with index $\Gamma = 2.15 \pm 0.06_{stat} \pm 0.20_{sys} $ while an F-test indicates that the power-law with exponential cutoff function provides a marginally-better fit ($\chi^2/NDF $ = $2.56 / 5 $) at the 2.1$\sigma$ level. The source location is found to be consistent with the published radio location and its spatial extent is consistent with a point source., Comment: 7 pages, 2 figures. Published in Monthly Notices of the Royal Astronomical Society

We present an analysis of the blank sky spectra observed with the Faint Object Spectrograph on board the Hubble Space Telescope. We study the diffuse sky emission from ultraviolet to optical wavelengths, which is composed of the zodiacal light (ZL), diffuse Galactic light (DGL), and residual emission. The observations were performed toward 54 fields distributed widely over the sky, with the spectral coverage from 0.2 to 0.7 um. In order to avoid contaminating light from the earthshine, we use the data collected only in orbital nighttime. The observed intensity is decomposed into the ZL, DGL, and residual emission, in eight photometric bands spanning our spectral coverage. We found that the derived ZL reflectance spectrum is flat in the optical, which indicates major contribution of C-type asteroids to the interplanetary dust (IPD). In addition, the ZL reflectance spectrum has an absorption feature at ~0.3 um. The shape of the DGL spectrum is consistent with those found in earlier measurements and model predictions. While the residual emission contains a contribution from the extragalactic background light, we found that the spectral shape of the residual looks similar to the ZL spectrum. Moreover, its optical intensity is much higher than that measured from beyond the IPD cloud by Pioneer10/11, and also than that of the integrated galaxy light. These findings may indicate the presence of an isotropic ZL component, which is missed in the conventional ZL models., This paper is dedicated to the memory of Dr. Kimiaki Kawara; accepted for publication in PASJ

The BL Lacertae object 1ES 1440+122 was observed in the energy range from 85 GeV to 30 TeV by the VERITAS array of imaging atmospheric Cherenkov telescopes. The observations, taken between 2008 May and 2010 June and totalling 53 hours, resulted in the discovery of $\gamma$-ray emission from the blazar, which has a redshift $z$=0.163. 1ES 1440+122 is detected at a statistical significance of 5.5 standard deviations above the background with an integral flux of (2.8$\pm0.7_{\mathrm{stat}}\pm0.8_{\mathrm{sys}}$) $\times$ 10$^{-12}$ cm$^{-2}$ s$^{-1}$ (1.2\% of the Crab Nebula's flux) above 200 GeV. The measured spectrum is described well by a power law from 0.2 TeV to 1.3 TeV with a photon index of 3.1 $\pm$ 0.4$_{\mathrm{stat}}$ $\pm$ 0.2$_{\mathrm{sys}}$. Quasi-simultaneous multi-wavelength data from the Fermi Large Area Telescope (0.3--300 GeV) and the Swift X-ray Telescope (0.2--10 keV) are additionally used to model the properties of the emission region. A synchrotron self-Compton model produces a good representation of the multi-wavelength data. Adding an external-Compton or a hadronic component also adequately describes the data., Comment: 8 pages, 4 figures. Accepted for publication in MNRAS

The extragalactic background light (EBL), i.e., the diffuse meta-galactic photon field in the ultraviolet to infrared, is dominated by the emission from stars in galaxies. It is, therefore, intimately connected with the integrated star formation rate density (SFRD). In this paper, the SFRD is constrained using recent limits on the EBL density derived from observations of distant sources of high and very-high energy gamma-rays. The stellar EBL contribution is modeled utilizing simple stellar population spectra including dust attenuation and emission. For modeling the SFRD up to z=4 a broken power law function in z+1 is assumed. A wide range of values for the different model parameters (SFRD(z), metallicity, dust absorption) is investigated and their impact on the resulting EBL is studied. The calculated EBL densities are compared with the specific EBL density limits and constraints on the SFRD are derived. For the fiducial model, adopting a Chabrier initial mass function (IMF) and a second power law index for the SFRD of beta=0.3, the SFRD is constrained to

The detections of X-ray emission from the kiloparsec-scale jets of blazars and radio galaxies may imply the existence of high energy electrons in these extended jets, and these electrons could produce high energy emission through inverse Compton (IC) process. In this paper we study the non-variable hard TeV emission from a blazar. The multi-band emission consists of two components: one is the traditional synchrotron self-Compton (SSC) emission from the inner jet, and the other is the emission produced via SSC and IC scattering of cosmic microwave background (CMB) photons (IC/CMB) and extragalactic background light (EBL) photons (IC/EBL) by relativistic electrons in the extended jet under the stochastic acceleration scenario. Such a model is applied to 1ES 1101-232. The results indicate that (1) the non-variable hard TeV emission of 1ES 1101-232 can be reproduced well, which is dominated by IC/CMB emission from the extended jet, using three characteristic values of Doppler factor ($\delta_{\rm D}=5,10,15$) for the TeV emitting region in the extended jet; and (2) in the cases of $\delta_{\rm D}=15$ and 10, the physical parameters can achieve the equipartition (or quasi-equipartition) between the relativistic electrons and the magnetic field; In contrast, the physical parameters largely deviate from the equipartition for the case of $\delta_{\rm D}=5$. It is therefore concluded that the TeV emission region of 1ES 1101-232 in the extended jet should be moderately or highly beamed.

BL Lacertae (Lac) objects that are detected at very-high energies (VHE) are of fundamental importance to study multiple astrophysical processes, including the physics of jets, the properties of the extragalactic background light and the strength of the intergalactic magnetic field. Unfortunately, since most blazars have featureless optical spectra that preclude a redshift determination, a substantial fraction of these VHE extragalactic sources cannot be used for cosmological studies. To assess whether molecular lines are a viable way to establish distances, we have undertaken a pilot program at the IRAM 30m telescope to search for CO lines in three BL Lac objects with known redshifts. We report a positive detection of M_H2 ~ 3x10^8 Msun toward 1ES 1959+650, but due to the poor quality of the baseline, this value is affected by a large systematic uncertainty. For the remaining two sources, W Comae and RGB J0710+591, we derive 3sigma upper limits at, respectively, M_H2 < 8.0x10^8 Msun and M_H2 < 1.6x10^9 Msun, assuming a line width of 150 km/s and a standard conversion factor alpha=4 M_sun/(K km/s pc^2). If these low molecular gas masses are typical for blazars, blind redshift searches in molecular lines are currently unfeasible. However, deep observations are still a promising way to obtain precise redshifts for sources whose approximate distances are known via indirect methods. Our observations further reveal a deficiency of molecular gas in BL Lac objects compared to quasars, suggesting that the host galaxies of these two types of active galactic nuclei (AGN) are not drawn from the same parent population. Future observations are needed to assess whether this discrepancy is statistically significant, but our pilot program shows how studies of the interstellar medium in AGN can provide key information to explore the connection between the active nuclei and the host galaxies.

Attenuation of high-energy gamma rays by pair-production with UV, op- tical and IR extragalactic background light (EBL) photons provides a link between the history of galaxy formation and high-energy astrophysics. We present results from our latest semi-analytic models (SAMs), which employ the main ingredients thought to be important to galaxy formation and evo- lution, as well as an improved model for reprocessing of starlight by dust to mid- and far-IR wavelengths. These SAMs are based upon aCDM hi- erarchical structural formation scenario, and are successful in reproducing a large variety of observational constraints such as number counts, luminosity and mass functions, and color bimodality. Our fiducial model is based upon a WMAP5 cosmology, and treats dust emission using empirical templates. This model predicts a background flux considerably lower than optical and near -IR measurements that rely on subtraction of zodiacal and galactic foregrounds, and near the lower bounds set by number counts of resolvable sources at a large number of wavelengths. We also show the results of varying several ingredi- ents of our SAM, including the cosmological parameters, the dust attenuation model, and the dust emission templates. For each of these we show how the optical depth from pair-production is affected by redshift and gamma-ray en- ergy, and the effect of gamma-ray absorption on the spectra of a variety of extragalactic sources. We conclude with a discussion of the implications of our work, comparisons to other models and key measurements of the EBL and a discussion of how the burgeoning science of gamma-ray astronomy will con- tinue to help constrain cosmology. The low EBL flux predicted by our fiducial model suggests an optimistic future for further studies of distant gamma-ray sources.

We combine a semi-analytic model of galaxy formation with simple analytic recipes describing the absorption and re-emission of starlight by dust in the interstellar medium of galaxies. We use the resulting models to predict galaxy counts and luminosity functions from the far-ultraviolet to the sub-mm, from redshift five to the present, and compare with an extensive compilation of observations. We find that in order to reproduce the rest-UV and optical luminosity functions at high redshift, we must assume an evolving normalization in the dust-to-metal ratio, implying that galaxies of a given bolometric luminosity (or metal column density) must be less extinguished than their local counterparts. In our best-fit model, we find remarkably good agreement with observations from rest-frame 1500 Angstroms to 250 microns. At longer wavelengths, most dramatically in the sub-mm, our models underpredict the number of bright galaxies by a large factor. The models reproduce the observed total IR luminosity function fairly well. We show the results of varying several ingredients of the models, including various aspects of the dust attenuation recipe, the dust emission templates, and the cosmology. We use our models to predict the integrated Extragalactic Background Light (EBL), and compare with an observationally-motivated EBL model and with other available observational constraints.

We study the quasi-simultaneous near-IR, optical, UV, and X-ray photometry of eleven gamma-ray selected blazars for which redshift estimates larger than 1.2 have been recently provided. Four of these objects turn out to be high-power blazars with the peak of their synchrotron emission between ~ 3 x 10^15 and ~ 10^16 Hz, and therefore of a kind predicted to exist but never seen before. This discovery has important implications for our understanding of physical processes in blazars, including the so-called "blazar sequence", and might also help constraining the extragalactic background light through gamma-ray absorption since two sources are strongly detected even in the 10 - 100 GeV Fermi-LAT band. Based on our previous work and their high powers, these sources are very likely high-redshift flat-spectrum radio quasars with their emission lines swamped by the non-thermal continuum.

The Fermi satellite has detected GeV emission from a number of gamma-ray bursts and active galactic nuclei at high redshift, z > 1.5. We examine the constraints that the detections of gamma rays from several of these sources place on the contribution of population-III stars to the extragalactic background light. Emission from these primordial stars, particularly redshifted Lyman-alpha emission, can interact with gamma rays to produce electron-positron pairs and create an optical depth to the propagation of gamma-ray emission, and the detection of emission above 10 GeV can therefore constrain the production of this background. We consider two initial mass functions for the early stars, and use derived SEDs for each to put upper limits on the star-formation rate density of massive early stars from redshifts 6 to 10. Our limits are complementary to those set on a high near-IR background flux by ground-based TeV-scale observations, and show that current data can limit star-formation in the late stages of reionization to less than 0.5 M_solar yr^-1 Mpc^-3. Our results also show that the total background flux from population-III stars must be considerably less than that from resolved galaxies at wavelengths below 1.5 microns.

Recent measurements by the Pierre Auger Observatory suggest that the composition of ultra-high energy cosmic rays (UHECRs) becomes dominated by heavy nuclei at high energies. However, until now there has been no astrophysical motivation for considering a source highly enriched in heavy elements. Here we demonstrate that the outflows from Gamma-Ray Bursts (GRBs) may indeed be composed primarily of nuclei with masses A ~ 40-200, which are synthesized as hot material expands away from the central engine. In particular, if the jet is magnetically-dominated (rather than a thermally-driven fireball) its low entropy enables heavy elements to form efficiently. Adopting the millisecond proto-magnetar model for the GRB central engine, we show that heavy nuclei are both synthesized in proto-magnetar winds and can in principle be accelerated to energies >1e20 eV in the shocks or regions of magnetic reconnection that are responsible for powering the GRB. Similar results may apply to accretion-powered GRB models if the jet originates from a magnetized disk wind. Depending on the precise distribution of nuclei synthesized, we predict that the average primary mass may continue to increase beyond Fe group elements at the highest energies, possibly reaching the A ~ 90 (Zirconium), A ~ 130 (Tellurium), or even A ~ 195 (Platinum) peaks. Future measurements of the UHECR composition at energies >~ 1e20 eV can thus confirm or constrain our model and, potentially, probe the nature of GRB outflows. The longer attenuation length of ultra-heavy nuclei through the extragalactic background light greatly expands the volume of accesible sources and alleviates the energetic constraints on GRBs as the source of UHECRs.

The prospects for future blazar surveys by next-generation very high energy (VHE) gamma-ray telescopes, such as Advanced Gamma-ray Imaging System (AGIS) and Cherenkov Telescope Array (CTA), are investigated using the latest model of blazar luminosity function and its evolution which is in good agreement with the flux and redshift distribution of observed blazars as well as the extragalactic gamma-ray background. We extend and improve the template of spectral energy distributions (SEDs) based on the blazar SED sequence paradigm, to make it reliable also in the VHE bands (above 100 GeV) by comparing with the existing VHE blazar data. Assuming the planned CTA sensitivities, a blind survey using a total survey time of ~100 hrs could detect ~3 VHE blazars, with larger expected numbers for wider/shallower surveys. We also discuss following-up of Fermi blazars. Detectability of VHE blazars in the plane of Fermi flux and redshift is presented, which would be useful for future survey planning. Prospects and strategies are discussed to constrain the extragalactic background light (EBL) by using the absorption feature of brightest blazar spectra, as well as cut-offs in the redshift distribution. We will be able to get useful constraints on EBL by VHE blazars at different redshifts ranging 0.3-1 TeV corresponding to z=0.10-0.36., 11 pages, 11 figures, accepted for publication in PASJ

There are a total of 1451 gamma-ray emitting objects in the Fermi Large Area Telescope First Source Catalogue. The point source location accuracy of typically a few arcminutes has allowed the counterparts for many of these sources to be found at other wavelengths, but even so there are 630 which are described as having no plausible counterpart at 80% confidence. In order to help identify the unknown objects, we have cross-correlated the positions of these sources with the Rosat All Sky Survey Bright Source Catalogue. In this way, for Fermi sources which have a possible counterpart in soft X-rays, we can use the, much smaller, Rosat error box to search for identifications. We find a strong correlation between the two samples and calculate that there are about 60 sources with a Rosat counterpart. Using the Rosat error boxes we provide tentative associations for half of them, demonstrate that the majority of these are either blazars or blazar candidates and give evidence that most belong to the BL Lac class. Given that they are X-ray selected and most are high synchrotron peaked objects, which indicates the presence of high energy electrons, these sources are also good candidates for TeV emission, and therefore good probes of the extragalactic background light.

We discuss a method to constrain the distance of blazars with unknown redshift using combined observations in the GeV and TeV regimes. We assume that the VHE spectrum corrected for the absorption through the interaction with the Extragalactic Background Light can not be harder than the spectrum in the Fermi/LAT band. Starting from the observed VHE spectral data we derive the EBL-corrected spectra as a function of the redshift z and fit them with power laws to be compared with power law fits to the LAT data. We apply the method to all TeV blazars detected by LAT with known distance and derive an empirical law describing the relation between the upper limits and the true redshifts that can be used to estimate the distance of unknown redshift blazars. Using different EBL models leads to systematic changes in the derived upper limits. Finally, we use this relation to infer the distance of the unknown redshift blazar PKS 1424+240.

Gamma rays from extragalactic sources are attenuated by pair-production interactions with diffuse photons of the extragalactic background light (EBL). Gamma-ray bursts (GRBs) are a source of high-redshift photons above 10 GeV, and could be therefore useful as a probe of the evolving UV background radiation. In this paper, we develop a simple phenomenological model for the number and redshift distribution of gamma-ray bursts that can be seen at GeV energies with the Fermi satellite and MAGIC atmospheric Cherenkov telescope. We estimate the observed number of gamma rays per year, and show how this result is modified by considering interactions with different realizations of the evolving EBL. We also discuss the bright Fermi GRB 080916C in the context of this model. We find that the LAT on Fermi can be expected to see a small number of photons above 10 GeV each year from distant GRBs. Annual results for ground-based instruments like MAGIC are highly variable due to the low duty cycle and sky coverage of the telescope. However, successfully viewing a bright or intermediate GRB from the ground could provide hundreds of photons from high redshift, which would almost certainly be extremely useful in constraining both GRB physics and the high-redshift EBL.

The MAGIC collaboration reported the detection of a new VHE source, MAGIC J0223+430, located close to the position of the blazar 3C66A, considered a candidate TeV blazar since a long time. A careful analysis showed that the events with energies above 150 GeV are centered on the position of the FRI radiogalaxy 3C66B (at 6 arcmin from 3C66A), with a probability of 95.4% (85.4% including systematic uncertainties) that the source is not related to 3C66A. We present a model for the possible emission of 3C66B based on the structured jet model already used to interpret the TeV emission of the radiogalaxy M87. The model requires parameters similar to those used for M87 but a larger luminosity for the layer, to account for the more luminous TeV emission. We also show that the spectrum obtained by MAGIC can be interpreted as the combined emission of 3C66B, dominating above ~200 GeV, and of 3C66A. The high-energy emission from the latter source, being strongly attenuated by the interaction with the extragalactic background light, can only contribute at low energies. If we were to see the jet emission of 3C 66B at small viewing angles we would see a spectral energy distribution closely resembling the one of S5 0716+714, a typical blazar., Comment: 5 pages, 2 figures, accepted for publication by MNRAS

We revisit the limits of the level of the extragalactic background light (EBL) recently reported by the MAGIC collaboration based on the observed gamma-ray spectrum of the quasar 3C279, considering the impact of absorption of high-energy gamma-ray photons inside the broad line region (BLR) of the quasar. We use the photoionization code CLOUDY to calculate the expected optical-UV radiation field inside the BLR and the optical depth to gamma-rays for a relatively extended set of the parameters. We found that the absorption of gamma-ray photons, though important for the estimate of the true radiative output of the source, does not produce an important hardening of the spectrum of 3C279 in the energy band accessible by MAGIC, supporting the method used to infer the upper limits to the level of the EBL.

Recent observations with the low threshold Cherenkov telescopes proved that sub-TeV gamma-rays are able to arrive from active galaxies at relatively large distances in spite of expected severe absorption in the extragalactic background light (EBL). We calculate the gamma-ray spectra at TeV energies from two example OVV quasars, 3C 273 and 3C 279, assuming that gamma-rays are injected in the inner parts of the jets launched by the accretion disks. It is assumed that gamma-rays in the broad energy range (from MeV up to TeV) are produced in these blazars with the power law spectrum with the spectral index as observed from these objects by the EGRET telescope at GeV energies. We take into account the internal absorption of these gamma-rays by considering a few different models for the radiation field surrounding the jet. The classical picture of the relativistic blob in jet model for the injection of primary gamma-rays is considered with the injection rate of gamma-rays as observed by the EGRET telescope in the GeV energy range. The results of calculations are compared with positive detection and the upper limits on the sub-TeV gamma-ray fluxes from these two sources. It is concluded that even with the Stecker EBL model, the level of the gamma-ray emission from 3C~279 is close to the recent measurements in the sub-TeV gamma-ray energies provided that the injected gamma-ray spectrum extends from the GeV energies over the next two decades with this same spectral index. We also suggest that a few day time scale flare from 3C 273 can be detected by the MAGIC II stereo telescopes.

The energy spectra of TeV gamma-rays from blazars, after being corrected for intergalatic absorption in the Extragalactic Background Light (EBL), appear unusually hard, a fact that poses challenges to the conventional models of particle acceleration in TeV blazars and/or to the EBL models. In this paper we show that the internal absorption of gamma-rays caused by interactions with dense narrow-band radiation fields in the vicinity of compact gamma-ray production regions can lead to the formation of gamma-ray spectra of an almost arbitrary hardness. This allows significant relaxation of the current tight constraints on particle acceleration and radiation models, although at the expense of enhanced requirements to the available nonthermal energy budget. The latter, however, is not a critical issue, as long as it can be largely compensated by the Doppler boosting, assuming very large ($\geq 30$) Doppler factors of the relativistically moving gamma-ray production regions. The suggested scenario of formation of hard gamma-ray spectra predicts detectable synchrotron radiation of secondary electron-positron pairs which might require a revision of the current ``standard paradigm'' of spectral energy distributions of gamma-ray blazars. If the primary gamma-rays are of hadronic origin related to $pp$ or $p \gamma$ interactions, the ``internal gamma-ray absorption'' model predicts neutrino fluxes close to the detection threshold of the next generation high energy neutrino detectors.

We present a comparison between the SCUBA Half Degree Extragalactic Survey (SHADES) at 450 and 850 microns in the Lockman Hole East with a deep Spitzer Space Telescope survey at 3.6-24 microns conducted in Guaranteed Time. Using stacking analyses we demonstrate a striking correspondence between the galaxies contributing the submm extragalactic background light, with those likely to dominate the backgrounds at Spitzer wavelengths. Using a combination BRIzK plus Spitzer photometric redshifts, we show that at least a third of the Spitzer-identified submm galaxies at 1 and specific star formation rates are in striking disagreement with some semi-analytic predictions from the Millenium simulation. The discrepancy could either be resolved with a top-heavy initial mass function, or a significant component of the submm flux heated by the interstellar radiation field., Comment: Accepted by MNRAS. SMG environments analysis extended. 17 pages, 18 figures. Includes BoxedEPS

Since 2002, the number of detected blazars at gamma-ray energies above 100 GeV has more than doubled. I study 17 blazars currently known to emit E>100 GeV gamma rays. Their intrinsic energy spectra are reconstructed by removing extragalactic background light attenuation effects. Luminosity and spectral slope in the E>100 GeV region are then compared and correlated among each other, with X-ray, optical and radio data, and with the estimated black hole (BH) masses of the respective host galaxies. According to expectations from synchrotron self-Compton emission models, a correlation on the 3.6-sigma significance level between gamma-ray and X-ray fluxes is found, while correlations between gamma-ray and optical/radio fluxes are less pronounced. Further, a general hardening of the E>100 GeV spectra with increasing gamma-ray luminosity is observed. This goes in line with a correlation of the gamma-ray luminosity and the synchrotron peak frequency, which is also seen. Tests for possible selection effects reveal a hardening of the spectra with increasing redshift. The gamma-ray emission might depend on the mass of the central BH. The studied blazars show no correlation of the BH masses with the spectral index and the luminosity in the E>100 GeV region. Also temporal properties of the X-ray and E>100 GeV gamma-ray flux are considered. No general trends are found, except that the blazars with the most massive BHs do not show particularly high duty cycles. In general, VHE flare time-scales are not found to scale with the BH mass. As a specific application of the luminosity study, a constraint for the still undetermined redshift of the blazar PG 1553+113 is discussed.

We present new determinations of the local sub-mm luminosity functions, solving the ``sub-mm Olbers' Paradox.'' We also present predictions of source counts and luminosity functions in current and future far-infrared to sub-mm surveys. Using the sub-mm colour temperature relations from the SCUBA Local Universe Galaxy Survey, and the discovery of excess 450 micron excess emission in these galaxies, we interpolate and extrapolate the IRAS detections to make predictions of the SEDs of all 15411 PSC-z galaxies from 50-1300 microns. Despite the long extrapolations we find excellent agreement with (a) the 90 micron luminosity function of Serjeant et al. (2001), (b) the 850 micron luminosity function of Dunne et al. (2000), (c) the mm-wave photometry of Andreani & Franceschini (1996); (d) the asymptotic differential and integral source count predictions at 50-1300 microns by Rowan-Robinson (2001). We find the local 850 micron sub-mm luminosity density converges to (7.3+/-0.2)x10^{19} h_{65} W/Hz/Mpc^3. Remarkably, the local spectral luminosity density and the extragalactic background light together strongly constrain the cosmic star formation history for a wide class of evolutionary assumptions. We find that the extragalactic background light, the 850 micron 8mJy source counts, and the Omega_* constraints all independently point to a decline in the comoving star formation rate at z>1. In order to reconcile this with direct determinations, we suggest either there is a top-heavy initial mass function at high redshifts, and/or there is stronger evolution in the more luminous far-infrared galaxies than seen in the population as a whole.

The large majority of sources detected by the Fermi Gamma-ray Space Telescope are blazars, belonging in particular to the blazar subclass of BL Lacertae objects (BL Lacs). BL Lacs often have featureless optical spectra, which make it difficult and sometimes impossible to determine their redshifts. This presents a severe impediment for the use of BL Lacs to measure the spectrum of the extragalactic background light through its interaction with high-energy gamma-ray photons. I present here high-quality optical spectroscopy of two of the brightest gamma-ray blazars, namely, PKS 0447-439 and PMN J0630-24. The medium-resolution and high signal-to-noise ratio optical spectra show clear absorption lines, which place these BL Lacs at relatively high redshifts of z>=1.246 for PKS 0447-439 and z>=1.238 for PMN J0630-24.

We present a multi-wavelength analysis of five submillimeter sources (S_1.1mm = 0.54-2.02 mJy) that were detected during our 1.1-mm-deep continuum survey in the SXDF-UDS-CANDELS field (2 arcmin^2, 1sigma = 0.055 mJy beam^-1) using the Atacama Large Millimeter/submillimeter Array (ALMA). The two brightest sources correspond to a known single-dish (AzTEC) selected bright submillimeter galaxy (SMG), whereas the remaining three are faint SMGs newly uncovered by ALMA. If we exclude the two brightest sources, the contribution of the ALMA-detected faint SMGs to the infrared extragalactic background light is estimated to be ~ 4.1^{+5.4}_{-3.0} Jy deg^{-2}, which corresponds to ~ 16^{+22}_{-12}% of the infrared extragalactic background light. This suggests that their contribution to the infrared extragalactic background light is as large as that of bright SMGs. We identified multi-wavelength counterparts of the five ALMA sources. One of the sources (SXDF-ALMA3) is extremely faint in the optical to near-infrared region despite its infrared luminosity (L_IR ~ 1e12 L_sun or SFR ~ 100 M_sun yr^{-1}). By fitting the spectral energy distributions (SEDs) at the optical-to-near-infrared wavelengths of the remaining four ALMA sources, we obtained the photometric redshifts (z_photo) and stellar masses (M_*): z_photo ~ 1.3-2.5, M_* ~ (3.5-9.5)e10 M_sun. We also derived their star formation rates (SFRs) and specific SFRs (sSFRs) as ~ 30-200 M_sun yr^{-1} and ~ 0.8-2 Gyr^{-1}, respectively. These values imply that they are main-sequence star-forming galaxies., PASJ accepted, 15 pages, 6 figures, 2 tables

The stellar haloes of galaxies can currently be studied either through observations of resolved halo stars or through surface photometry. Curiously, the two methods appear to give conflicting results, as a number of surface photometry measurements have revealed integrated colours that are too red to be reconciled with the halo properties inferred from the study of resolved stars. Several explanations for this anomaly have been proposed - including dust photoluminescence, extinction of extragalactic background light and a bottom-heavy stellar initial mass function. A decisive test is, however, still lacking. Here, we explain how observations of the halo of a nearby galaxy, involving a combination of both surface photometry and bright star counts, can be used to distinguish between the proposed explanations. We derive the observational requirements for this endeavour and find that star counts in filters VI and surface photometry in filters VIJ appears to be the optimal strategy. Since the required halo star counts are already available for many nearby galaxies, the most challenging part of this test is likely to be the optical surface photometry, which requires several nights of exposure time on a 4-8 m telescope, and the near-IR surface photometry, which is most readily carried out using the upcoming James Webb Space Telescope.

In this paper, we describe a search for very low-surface-brightness features in deep CCD frames. We show how cross-correlation with templates can be used to detect low-surface-brightness features, which would not be detected using more typical connected-pixel algorithms. The search technique is specifically designed to search for large, low-surface- brightness galaxies like Malin 1. Although we identify 19 potentially interesting, extended, low-surface-brightness objects, we conclude that Malin 1 type objects per se are less prolific than `normal' L^*^ galaxies, and that galaxies with -22 < M_v_ < 19 and μ_0_ ~26 Vmu are at least an order of magnitude less common than their normal-surface- brightness counterparts.

TeV photons from blazars at relatively large distances, interacting with the optical-IR cosmic background, are efficiently converted into electron-positron pairs. The produced pairs are extremely relativistic (Lorentz factors of the order of 1e6 1e7 and promptly loose their energy through inverse Compton scatterings with the photons of the microwave cosmic background, producing emission in the GeV band. The spectrum and the flux level of this reprocessed emission is critically dependent on the intensity of the intergalactic magnetic field, B, that can deflect the pairs diluting the intrinsic emission over a large solid angle. We derive a simple relation for the reprocessed spectrum expected from a steady source. We apply this treatment to the blazar 1ES 0229+200, whose intrinsic very hard TeV spectrum is expected to be approximately steady. Comparing the predicted reprocessed emission with the upper limits measured by the Fermi/Large Area Telescope, we constrain the value of the intergalactic magnetic field to be larger than $B \simeq 5\times 10^{-15}$ Gauss, depending on the model of extragalactic background light.

The problem concerning the origin of the Intergalactic Magnetic Field (IGMF) is one of the long-standing problems of astrophysics and cosmology, and direct measurements are difficult. TeV photons emitted by TeV blazars produce electron-positron pairs because of interactions with the extragalactic background light (EBL). These pairs emit secondary cascade gamma-rays via Inverse Compton scattering of the cosmic microwave background (CMB) photons. In this process, the trajectories of the pairs are deviated by the IGMF, and the cascade gamma-ray emission appears as extended emission around TeV source. We used the EBL, synchrotron self-Compton (SSC) model and the observed limits on Lorentz factor of electron-positron pairs to calculate the cascade-radiation spectrum, and then to fit the observed GeV to TeV and multi-waveband spectra of TeV blazars to constrain the IGMF. We obtained the GeV energy spectra of three TeV blazars by analyzing the Fermi-LAT data of the past similar to 3 yr. The flux upper limits of Fermi-LAT in the 90% significant level of 1ES 0229+200 suggests that the IGMF is stronger than 2 x 10(-18) G for an engine time of TeV activity with three years. The relationships between the deduced lower limits of IGMF and various engine times for 1ES 0229+200 and 1ES 0347-121 are presented by us.