Source author record

F. Prada

F. Prada appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

96works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

96 published item(s)

preprint2026arXiv

Baryon Acoustic Oscillations from the C IV Forest with DESI DR2

We present a measurement of Baryon Acoustic Oscillations (BAO) in the cross-correlation of triply ionized carbon C IV absorption with the positions of quasars (QSO) and Emission Line Galaxies (ELG). We use quasars and ELGs from the second data release (DR2) of the Dark Energy Spectroscopic Instrument (DESI) survey. Our data sample consists of 2.5 million quasars, 3.1 million ELGs, and the C IV absorption is measured along the line of sight of 1.5 million high redshift quasars with $z > 1.3$. We measure the isotropic BAO signal at 4.2$σ$ for the CIV$\times$QSO cross-correlation. This translates into a 3.0% precision measurement of the ratio of the isotropic distance scale, $D_{\rm V}$, and the sound horizon at the drag epoch, $r_{\rm d}$, with $D_{\rm V}/r_{\rm d}(z_{\rm eff} = 1.92) = 30.3 \pm 0.9$. We make the first detection of the BAO feature in the CIV$\times$ELG cross-correlation at a significance of 2.5$σ$ and find $D_{\rm V}/r_{\rm d}(z_{\rm eff} = 1.47) = 24.6 \pm 1.0$.

preprint2026arXiv

DESI Data Release 1: Stellar Catalogue

In this paper we present the stellar Value-Added Catalogue (VAC) based on the DESI Data Release 1. This VAC contains stellar parameter, abundance and radial velocity measurements for more than 4 million stars. It also contains, for the first time, measurements from individual epochs for more than a million stars with at least two observations. The main contribution to the catalogue comes from the bright program of the main survey, which includes $\sim $2.5 million stars, and the backup program, which includes $\sim $ 1 million stars. The combined magnitude range for the stars in the catalogue extends from ${\it Gaia}$ G $\sim 12$ to G $\sim 21$. For the magnitude range $17.5<G<21$ this catalogue represents a factor of 10 increase in the number of stars with radial velocity and abundance measurements compared to existing surveys. Despite DESI's resolution (R $\sim 2500-5000$), the median radial velocity uncertainty for stars in the catalogue is better than 1 km s$^{-1}$. The stellar parameters and abundances of stars in DESI are measured by two independent pipelines, and after applying a temperature-dependent calibration, [Fe/H] abundances of high signal-to-noise stars are accurate to better than $\sim$ 0.1 dex when compared to high-resolution surveys. The catalogue probes different Galactic components including a particularly large number of distant stars: tens of thousands of stars further than 10 kpc, and thousands further than 50 kpc. The catalogue also contains several thousand extremely metal-poor stars with ${\rm [Fe/H]}<-3$. The released sample of stars includes measurements for thousands of stars that are members of dwarf galaxies, open and globular clusters as well as members of several dozen stellar streams. The next public DESI data release is expected in less than two years and will contain three times as many stars as DR1.

preprint2026arXiv

Electron temperature relations and the direct N, O, Ne, S and Ar abundances of 49959 star-forming galaxies in DESI Data Release 2

We present the largest direct-method abundance catalogue of galaxies to date, containing measurements of 49$\,$959 star-forming galaxies at $z < 0.96$ from DESI Data Release 2. By directly measuring electron temperatures across multiple ionization zones, we provide constraints on a number of electron temperature relations finding good consistency with previous literature relations. Using these temperature measurements, we derive reliable abundances for N, O, Ne, S and Ar and measure the evolution of abundances and abundance ratios of as a function of metallicity and other galaxy properties. Our measurements include direct oxygen abundances for 49$\,$766 galaxies, leading to the discovery of the two most metal-poor galaxies in the nearby Universe, with oxygen abundances of $\rm 12+\log(O/H) = 6.77_{-0.03}^{+0.03}~\rm dex$ (1.2\% $\rm Z_{\odot}$) and $\rm 12+\log(O/H) = 6.81_{-0.04}^{+0.04}~\rm dex$ (1.3\% $\rm Z_{\odot}$). We identify a rare outlier population of 139 galaxies with high N/O ratios at low metallicity, reminiscent of galaxy abundances observed in the early Universe. We find these high N/O galaxies are more massive than typical galaxies at the same metallicity. We find the Ne/O ratio is constant at low metallicity but increases significantly at $\rm 12+log(O/H) > 8.105\pm0.004$ dex. We show that the S/O and Ar/O abundance ratios are strongly correlated, consistent with the expected additional Type Ia enrichment channel for S and Ar. In this work we present an initial survey of the key properties of the sample, with this dataset serving as a foundation for extensive future work on galaxy abundances at low redshift.

preprint2026arXiv

The Linear Point Standard Ruler with DESI DR1 and DR2 Data

The linear point, a purely geometric feature in the monopole of the two-point correlation function, has been proposed as an alternative standard ruler. Compared to the peak in the correlation function, it is more robust to late-time nonlinear effects at the percent level. In light of improved simulations and high quality data, we revisit the robustness of the linear point and use it as an alternative to template-based fitting approaches typically used in BAO analyses. We present the linear point measurements on galaxy samples from the first and second data releases (DR1 and DR2) of the DESI survey. We convert the linear point into a dimensionless parameter $α_{iso,LP}$, defined as the ratio of the linear point in the fiducial cosmology and the observed value, analogous to the isotropic BAO scaling parameter $α_{iso}$ used in previous BAO measurements. Using the 2nd generation of AbacusSummit mock catalogs, we find that linear point measurements are more precise when calculated in the post-reconstruction regime with 15-60% smaller uncertainties than those pre-reconstruction. We find a systematic shift in the linear point measurements compared against the isotropic BAO measurements in mocks; we attribute this to the isotropic damping parameter responsible for smearing the linear point in the nonlinear regime. We propose a sample-dependent correction that mitigates the impact of late-time nonlinear effects. While this introduces a cosmology dependence in an otherwise model-independent measurement, this is necessary given the sub-percent precision dictated by current cosmological surveys. Comparing $α_{iso,LP}$ with isotropic BAO measurements made on the DESI DR1 and DR2 galaxy samples, we find excellent agreement after applying this correction, particularly post-reconstruction. We discuss future scope regarding cosmological inference with linear point measurements.

preprint2025arXiv

The Compilation and Validation of the Spectroscopic Redshift Catalogs for the DESI-COSMOS and DESI-XMMLSS Fields

Over several dedicated programs that include targets beyond the main cosmological samples, the Dark Energy Spectroscopic Instrument (DESI) collected spectra for 304,970 unique objects in two fields centered on the COSMOS and XMM-LSS fields. In this work, we develop spectroscopic redshift robustness criteria for those spectra, validate these criteria using visual inspection, and provide two custom Value-Added Catalogs with our redshift characterizations. With these criteria, we reliably classify 212,935 galaxies below z < 1.6, 9,713 quasars and 35,222 stars. As a critical element in characterizing the selection function, we provide the description of 70 different algorithms that were used to select these targets from imaging data. To facilitate joint imaging/spectroscopic analyses, we provide row-matched photometry from the Dark Energy Camera, Hyper-Suprime Cam, and public COSMOS2020 photometric catalogs. Finally, we demonstrate example applications of these large catalogs to photometric redshift estimation, cluster finding, and completeness studies.

preprint2022arXiv

Target Selection and Validation of DESI Luminous Red Galaxies

The Dark Energy Spectroscopic Instrument (DESI) is carrying out a 5-year survey that aims to measure the redshifts of tens of millions of galaxies and quasars, including 8 million luminous red galaxies (LRGs) in the redshift range of $0.4<z<{\sim}\,1.0$. Here we present the selection of the DESI LRG sample and assess its spectroscopic performance using data from Survey Validation (SV) and the first 2 months of the Main Survey. The DESI LRG sample, selected using $g$, $r$, $z$, and $W1$ photometry from the DESI Legacy Imaging Surveys, is highly robust against imaging systematics. The sample has a target density of 605 deg$^{-2}$ and a comoving number density of $5\times10^{-4}\ h^3\mathrm{Mpc}^{-3}$ in $0.4<z<0.8$; this is a significantly higher density than previous LRG surveys (such as SDSS, BOSS and eBOSS) while also extending to $z \sim 1$. After applying a bright star veto mask developed for the sample, $98.9\%$ of the observed LRG targets yield confident redshifts (with a catastrophic failure rate of $0.2\%$ in the confident redshifts), and only $0.5\%$ of the LRG targets are stellar contamination. The LRG redshift efficiency varies with source brightness and effective exposure time, and we present a simple model that accurately characterizes this dependence. In the appendices, we describe the extended LRG samples observed during SV.

preprint2022arXiv

The Uchuu-SDSS galaxy lightcones: a clustering, RSD and BAO study

We present the data release of the Uchuu-SDSS galaxies: a set of 32 high-fidelity galaxy lightcones constructed from the large Uchuu 2.1 trillion particle $N$-body simulation using Planck cosmology. We adopt subhalo abundance matching to populate the Uchuu-box halo catalogues with SDSS galaxy luminosities. These cubic box galaxy catalogues generated at several redshifts are combined to create the set of lightcones with redshift-evolving galaxy properties. The Uchuu-SDSS galaxy lightcones are built to reproduce the footprint and statistical properties of the SDSS main galaxy survey, along with stellar masses and star formation rates. This facilitates direct comparison of the observed SDSS and simulated Uchuu-SDSS data. Our lightcones reproduce a large number of observational results, such as the distribution of galaxy properties, the galaxy clustering, the stellar mass functions, and the halo occupation distributions. Using the simulated and real data we select samples of bright red galaxies at $z_\mathrm{eff}=0.15$ to explore Redshift Space Distortions and Baryon Acoustic Oscillations (BAO) utilizing a full-shape analytical model of the two-point correlation function. We create a set of 5100 galaxy lightcones using GLAM N-body simulations to compute covariance errors. We report a $\sim 30\%$ precision increase on $fσ_8$, due to our better estimate of the covariance matrix. From our BAO-inferred $α_{\parallel}$ and $α_{\perp}$ parameters, we obtain the first SDSS measurements of the Hubble and angular diameter distances $D_\mathrm{H}(z=0.15) / r_d = 27.9^{+3.1}_{-2.7}$, $D_\mathrm{M}(z=0.15) / r_d = 5.1^{+0.4}_{-0.4}$. Overall, we conclude that the Planck LCDM cosmology nicely explains the observed large-scale structure statistics of SDSS. All data sets are made publicly available.

preprint2021arXiv

[OII] emitters in MultiDark-Galaxies and DEEP2

We use three semi-analytic models (SAMs) of galaxy formation and evolution, run on the same 1$h^{-1}$Gpc MultiDark Planck2 cosmological simulation, to investigate the properties of [OII] emission line galaxies in the redshift range $0.6<z<1.2$. We compare model predictions with different observational data sets, including DEEP2--Firefly galaxies with absolute magnitudes. We estimate the [OII] luminosity, L[OII], using simple relations derived both from the models and observations and also using a public code. This code ideally uses as input instantaneous star formation rates (SFRs), which are only provided by one of the SAMs under consideration. We use this SAM to study the feasibility of inferring galaxies' L[OII] for models that only provide average SFRs. We find that the post-processing computation of L[OII] from average SFRs is accurate for model galaxies with dust attenuated L[OII]$\lesssim10^{42.2}$erg s$^{-1}$ ($<5\%$ discrepancy). We also explore how to derive the [OII] luminosity from simple relations using global properties usually output by SAMs. Besides the SFR, the model L[OII] is best correlated with the observed-frame $u$ and $g$ broad-band magnitudes. These correlations have coefficients (r-values) above 0.64 and a dispersion that varies with L[OII]. We use these correlations and an observational one based on SFR and metallicity to derive L[OII]. These relations result in [OII] luminosity functions and halo occupation distributions with shapes that vary depending on both the model and the method used. Nevertheless, for all the considered models, the amplitude of the clustering at scales above 1$h^{-1}$Mpc remains unchanged independently of the method used to derive L[OII].

preprint2021arXiv

Testing the 10 spectrograph units for DESI: approach and results

The recently commissioned Dark Energy Spectroscopic Instrument (DESI) will measure the expansion history of the Universe using the Baryon Acoustic Oscillation technique. The spectra of 35 million galaxies and quasars over 14000 sqdeg will be measured during the life of the experiment. A new prime focus corrector for the KPNO Mayall telescope delivers light to 5000 fiber optic positioners. The fibers in turn feed ten broad-band spectrographs. A consortium of Aix-Marseille University (AMU) and CNRS laboratories (LAM, OHP and CPPM) together with LPNHE (CNRS, IN2P3, Sorbonne Université and Université de Paris) and the WINLIGHT Systems company based in Pertuis (France), were in charge of integrating and validating the performance requirements of the ten full spectrographs, equipped with their cryostats, shutters and other mechanisms. We present a summary of our activity which allowed an efficient validation of the systems in a short-time schedule. We detail the main results. We emphasize the benefits of our approach and also its limitations.

preprint2020arXiv

CASTLE: performances and science cases

We present here the Calar Alto Schmidt-Lemaitre Telescope (CASTLE) concept, a technology demonstrator for curved detectors, that will be installed at the Calar Alto Observatory (Spain). This telescope has a wide field of view (2.36x1.56 deg^2) and a design, optimised to generate a Point Spread Function with very low level wings and reduced ghost features, which makes it considerably less susceptible to several systematic effects usually affecting similar systems. These characteristics are particularly suited to study the low surface brightness Universe. CASTLE will be able to reach surface brightness orders of magnitude fainter than the sky background level and observe the extremely extended and faint features around galaxies such as tidal features, stellar halos, intra-cluster light, etc. CASTLE will also be used to search and detect astrophysical transients such as gamma ray bursts (GRB), gravitational wave optical counterparts, neutrino counterparts, etc. This will increase the number of precisely localized GRBs from 20% to 60% (in the case of Fermi/GMB GRBs).

preprint2016arXiv

A $K_s$-band selected catalogue of objects in the ALHAMBRA survey

The original ALHAMBRA catalogue contained over 400,000 galaxies selected using a synthetic F814W image, to the magnitude limit AB(F814W)$\approx$24.5. Given the photometric redshift depth of the ALHAMBRA multiband data (<z>=0.86) and the approximately $I$-band selection, there is a noticeable bias against red objects at moderate redshift. We avoid this bias by creating a new catalogue selected in the $K_s$ band. This newly obtained catalogue is certainly shallower in terms of apparent magnitude, but deeper in terms of redshift, with a significant population of red objects at $z>1$. We select objects using the $K_s$ band images, which reach an approximate AB magnitude limit $K_s \approx 22$. We generate masks and derive completeness functions to characterize the sample. We have tested the quality of the photometry and photometric redshifts using both internal and external checks. Our final catalogue includes $\approx 95,000$ sources down to $K_s \approx 22$, with a significant tail towards high redshift. We have checked that there is a large sample of objects with spectral energy distributions that correspond to that of massive, passively evolving galaxies at $z > 1$, reaching as far as $z \approx 2.5$. We have tested the possibility of combining our data with deep infrared observations at longer wavelengths, particularly Spitzer IRAC data.

preprint2016arXiv

SDSS-IV eBOSS emission-line galaxy pilot survey

The Sloan Digital Sky Survey IV extended Baryonic Oscillation Spectroscopic Survey (SDSS-IV/eBOSS) will observe 195,000 emission-line galaxies (ELGs) to measure the Baryonic Acoustic Oscillation standard ruler (BAO) at redshift 0.9. To test different ELG selection algorithms, 9,000 spectra were observed with the SDSS spectrograph as a pilot survey based on data from several imaging surveys. First, using visual inspection and redshift quality flags, we show that the automated spectroscopic redshifts assigned by the pipeline meet the quality requirements for a reliable BAO measurement. We also show the correlations between sky emission, signal-to-noise ratio in the emission lines, and redshift error. Then we provide a detailed description of each target selection algorithm we tested and compare them with the requirements of the eBOSS experiment. As a result, we provide reliable redshift distributions for the different target selection schemes we tested. Finally, we determine an target selection algorithms that is best suited to be applied on DECam photometry because they fulfill the eBOSS survey efficiency requirements.

preprint2016arXiv

The ALHAMBRA survey: evolution of galaxy spectral segregation

We study the clustering of galaxies as a function of spectral type and redshift in the range $0.35 < z < 1.1$ using data from the Advanced Large Homogeneous Area Medium Band Redshift Astronomical (ALHAMBRA) survey. The data cover 2.381 deg$^2$ in 7 fields, after applying a detailed angular selection mask, with accurate photometric redshifts [$σ_z < 0.014(1+z)$] down to $I_{AB} < 24$. From this catalog we draw five fixed number density, redshift-limited bins. We estimate the clustering evolution for two different spectral populations selected using the ALHAMBRA-based photometric templates: quiescent and star-forming galaxies. For each sample, we measure the real-space clustering using the projected correlation function. Our calculations are performed over the range $[0.03,10.0] h^{-1}$ Mpc, allowing us to find a steeper trend for $r_p \lesssim 0.2 h^{-1}$ Mpc, which is especially clear for star-forming galaxies. Our analysis also shows a clear early differentiation in the clustering properties of both populations: star-forming galaxies show weaker clustering with evolution in the correlation length over the analysed redshift range, while quiescent galaxies show stronger clustering already at high redshifts, and no appreciable evolution. We also perform the bias calculation where similar segregation is found, but now it is among the quiescent galaxies where a growing evolution with redshift is clearer. These findings clearly corroborate the well known colour-density relation, confirming that quiescent galaxies are mainly located in dark matter halos that are more massive than those typically populated by star-forming galaxies.

preprint2016arXiv

The Low Redshift survey at Calar Alto (LoRCA)

The Baryon Acoustic Oscillation (BAO) feature in the power spectrum of galaxies provides a standard ruler to measure the accelerated expansion of the Universe. To extract all available information about dark energy, it is necessary to measure a standard ruler in the local, z<0.2, universe where dark energy dominates most the energy density of the Universe. Though the volume available in the local universe is limited, it is just big enough to measure accurately the long 100 Mpc/h wave-mode of the BAO. Using cosmological N-body simulations and approximate methods based on Lagrangian perturbation theory, we construct a suite of a thousand light-cones to evaluate the precision at which one can measure the BAO standard ruler in the local universe. We find that using the most massive galaxies on the full sky (34,000 sq. deg.), i.e. a K(2MASS)<14 magnitude-limited sample, one can measure the BAO scale up to a precision of 4\% and 1.2\% using reconstruction). We also find that such a survey would help to detect the dynamics of dark energy.Therefore, we propose a 3-year long observational project, named the Low Redshift survey at Calar Alto (LoRCA), to observe spectroscopically about 200,000 galaxies in the northern sky to contribute to the construction of aforementioned galaxy sample. The suite of light-cones is made available to the public.

preprint2016arXiv

The SDSS-IV eBOSS: emission line galaxy catalogues at z=0.8 and study of systematic errors in the angular clustering

We present two wide-field catalogs of photometrically-selected emission line galaxies (ELGs) at z=0.8 covering about 2800 deg^2 over the south galactic cap. The catalogs were obtained using a Fisher discriminant technique described in a companion paper. The two catalogs differ by the imaging used to define the Fisher discriminant: the first catalog includes imaging from the Sloan Digital Sky Survey and the Wide-Field Infrared Survey Explorer, the second also includes information from the South Galactic Cap U-band Sky Survey (SCUSS). Containing respectively 560,045 and 615,601 objects, they represent the largest ELG catalogs available today and were designed for the ELG programme of the extended Baryon Oscillation Spectroscopic Survey (eBOSS). We study potential sources of systematic variation in the angular distribution of the selected ELGs due to fluctuations of the observational parameters. We model the influence of the observational parameters using a multivariate regression and implement a weighting scheme that allows effective removal of all of the systematic errors induced by the observational parameters. We show that fluctuations in the imaging zero-points of the photometric bands have minor impact on the angular distribution of objects in our catalogs. We compute the angular clustering of both catalogs and show that our weighting procedure effectively removes spurious clustering on large scales. We fit a model to the small scale angular clustering, showing that the selections have similar biases of 1.35/D_a(z) and 1.28/D_a(z). Both catalogs are publicly available.

preprint2016arXiv

The SDSS-IV extended Baryon Oscillation Spectroscopic Survey: selecting emission line galaxies using the Fisher discriminant

We present a new selection technique of producing spectroscopic target catalogues for massive spectroscopic surveys for cosmology. This work was conducted in the context of the extended Baryon Oscillation Spectroscopic Survey (eBOSS), which will use ~200 000 emission line galaxies (ELGs) at 0.6<zspec<1.0 to obtain a precise baryon acoustic oscillation measurement. Our proposed selection technique is based on optical and near-infrared broad-band filter photometry. We used a training sample to define a quantity, the Fisher discriminant (linear combination of colours), which correlates best with the desired properties of the target: redshift and [OII] flux. The proposed selections are simply done by applying a cut on magnitudes and this Fisher discriminant. We used public data and dedicated SDSS spectroscopy to quantify the redshift distribution and [OII] flux of our ELG target selections. We demonstrate that two of our selections fulfil the initial eBOSS/ELG redshift requirements: for a target density of 180 deg^2, ~70% of the selected objects have 0.6<zspec<1.0 and only ~1% of those galaxies in the range 0.6<zspec<1.0 are expected to have a catastrophic zspec estimate. Additionally, the stacked spectra and stacked deep images for those two selections show characteristic features of star-forming galaxies. The proposed approach using the Fisher discriminant could, however, be used to efficiently select other galaxy populations, based on multi-band photometry, providing that spectroscopic information is available. This technique could thus be useful for other future massive spectroscopic surveys such as PFS, DESI, and 4MOST.

preprint2015arXiv

High redshift galaxies in the ALHAMBRA survey: I. selection method and number counts based on redshift PDFs

Context. Most observational results on the high redshift restframe UV-bright galaxies are based on samples pinpointed using the so called dropout technique or Ly-alpha selection. However, the availability of multifilter data allows now replacing the dropout selections by direct methods based on photometric redshifts. In this paper we present the methodology to select and study the population of high redshift galaxies in the ALHAMBRA survey data. Aims. Our aim is to develop a less biased methodology than the traditional dropout technique to study the high redshift galaxies in ALHAMBRA and other multifilter data. Thanks to the wide area ALHAMBRA covers, we especially aim at contributing in the study of the brightest, less frequent, high redshift galaxies. Methods. The methodology is based on redshift probability distribution functions (zPDFs). It is shown how a clean galaxy sample can be obtained by selecting the galaxies with high integrated probability of being within a given redshift interval. However, reaching both a complete and clean sample with this method is challenging. Hence, a method to derive statistical properties by summing the zPDFs of all the galaxies in the redshift bin of interest is introduced. Results. Using this methodology we derive the galaxy rest frame UV number counts in five redshift bins centred at z=2.5, 3.0, 3.5, 4.0, and 4.5, being complete up to the limiting magnitude at m_UV(AB)=24. With the wide field ALHAMBRA data we especially contribute in the study of the brightest ends of these counts, sampling well the surface densities down to m_UV(AB)=21-22. Conclusions. We show that using the zPDFs it is easy to select a clean sample of high redshift galaxies. We also show that statistical analysis of the properties of galaxies is better done using a probabilistic approach, which takes into account both the incompleteness and contamination in a natural way.

preprint2015arXiv

Photometric redshifts and clustering of emission line galaxies selected jointly by DES and eBOSS

We present the results of the first test plates of the extended Baryon Oscillation Spectroscopic Survey. This paper focuses on the emission line galaxies (ELG) population targetted from the Dark Energy Survey (DES) photometry. We analyse the success rate, efficiency, redshift distribution, and clustering properties of the targets. From the 9000 spectroscopic redshifts targetted, 4600 have been selected from the DES photometry. The total success rate for redshifts between 0.6 and 1.2 is 71\% and 68\% respectively for a bright and faint, on average more distant, samples including redshifts measured from a single strong emission line. We find a mean redshift of 0.8 and 0.87, with 15 and 13\% of unknown redshifts respectively for the bright and faint samples. In the redshift range 0.6<z<1.2, for the most secure spectroscopic redshifts, the mean redshift for the bright and faint sample is 0.85 and 0.9 respectively. Star contamination is lower than 2\%. We measure a galaxy bias averaged on scales of 1 and 10~Mpc/h of 1.72 \pm 0.1 for the bright sample and of 1.78 \pm 0.12 for the faint sample. The error on the galaxy bias have been obtained propagating the errors in the correlation function to the fitted parameters. This redshift evolution for the galaxy bias is in agreement with theoretical expectations for a galaxy population with MB-5\log h < -21.0. We note that biasing is derived from the galaxy clustering relative to a model for the mass fluctuations. We investigate the quality of the DES photometric redshifts and find that the outlier fraction can be reduced using a comparison between template fitting and neural network, or using a random forest algorithm.

preprint2015arXiv

Stellar populations of galaxies in the ALHAMBRA survey up to $z \sim 1$. I. MUFFIT: A Multi-Filter Fitting code for stellar population diagnostics

We present MUFFIT, a new generic code optimized to retrieve the main stellar population parameters of galaxies in photometric multi-filter surveys, and we check its reliability and feasibility with real galaxy data from the ALHAMBRA survey. Making use of an error-weighted $χ^2$-test, we compare the multi-filter fluxes of galaxies with the synthetic photometry of mixtures of two single stellar populations at different redshifts and extinctions, to provide through a Monte Carlo method the most likely range of stellar population parameters (mainly ages and metallicities), extinctions, redshifts, and stellar masses. To improve the diagnostic reliability, MUFFIT identifies and removes from the analysis those bands that are significantly affected by emission lines. We highlight that the retrieved age-metallicity locus for a sample of $z \le 0.22$ early-type galaxies in ALHAMBRA at different stellar mass bins are in very good agreement with the ones from SDSS spectroscopic diagnostics. Moreover, a one-to-one comparison between the redshifts, ages, metallicities, and stellar masses derived spectroscopically for SDSS and by MUFFIT for ALHAMBRA reveals good qualitative agreements in all the parameters. In addition, and using as input the results from photometric-redshift codes, MUFFIT improves the photometric-redshift accuracy by $\sim 10$-$20\%$, and it also detects nebular emissions in galaxies, providing physical information about their strengths. Our results show the potential of multi-filter galaxy data to conduct reliable stellar population studies with the appropiate analysis techniques, as MUFFIT.

preprint2015arXiv

The 2009 multiwavelength campaign on Mrk 421: Variability and correlation studies

We performed a 4.5-month multi-instrument campaign (from radio to VHE gamma rays) on Mrk421 between January 2009 and June 2009, which included VLBA, F-GAMMA, GASP-WEBT, Swift, RXTE, Fermi-LAT, MAGIC, and Whipple, among other instruments and collaborations. Mrk421 was found in its typical (non-flaring) activity state, with a VHE flux of about half that of the Crab Nebula, yet the light curves show significant variability at all wavelengths, the highest variability being in the X-rays. We determined the power spectral densities (PSD) at most wavelengths and found that all PSDs can be described by power-laws without a break, and with indices consistent with pink/red-noise behavior. We observed a harder-when-brighter behavior in the X-ray spectra and measured a positive correlation between VHE and X-ray fluxes with zero time lag. Such characteristics have been reported many times during flaring activity, but here they are reported for the first time in the non-flaring state. We also observed an overall anti-correlation between optical/UV and X-rays extending over the duration of the campaign. The harder-when-brighter behavior in the X-ray spectra and the measured positive X-ray/VHE correlation during the 2009 multi-wavelength campaign suggests that the physical processes dominating the emission during non-flaring states have similarities with those occurring during flaring activity. In particular, this observation supports leptonic scenarios as being responsible for the emission of Mrk421 during non-flaring activity. Such a temporally extended X-ray/VHE correlation is not driven by any single flaring event, and hence is difficult to explain within the standard hadronic scenarios. The highest variability is observed in the X-ray band, which, within the one-zone synchrotron self-Compton scenario, indicates that the electron energy distribution is most variable at the highest energies.

preprint2015arXiv

The ALHAMBRA survey : Estimation of the clustering signal encoded in the cosmic variance

The relative cosmic variance ($σ_v$) is a fundamental source of uncertainty in pencil-beam surveys and, as a particular case of count-in-cell statistics, can be used to estimate the bias between galaxies and their underlying dark-matter distribution. Our goal is to test the significance of the clustering information encoded in the $σ_v$ measured in the ALHAMBRA survey. We measure the cosmic variance of several galaxy populations selected with $B-$band luminosity at $0.35 \leq z < 1.05$ as the intrinsic dispersion in the number density distribution derived from the 48 ALHAMBRA subfields. We compare the observational $σ_v$ with the cosmic variance of the dark matter expected from the theory, $σ_{v,{\rm dm}}$. This provides an estimation of the galaxy bias $b$. The galaxy bias from the cosmic variance is in excellent agreement with the bias estimated by two-point correlation function analysis in ALHAMBRA. This holds for different redshift bins, for red and blue subsamples, and for several $B-$band luminosity selections. We find that $b$ increases with the $B-$band luminosity and the redshift, as expected from previous work. Moreover, red galaxies have a larger bias than blue galaxies, with a relative bias of $b_{\rm rel} = 1.4 \pm 0.2$. Our results demonstrate that the cosmic variance measured in ALHAMBRA is due to the clustering of galaxies and can be used to characterise the $σ_v$ affecting pencil-beam surveys. In addition, it can also be used to estimate the galaxy bias $b$ from a method independent of correlation functions.

preprint2015arXiv

The impact from survey depth and resolution on the morphological classification of galaxies

We consistently analyse for the first time the impact of survey depth and spatial resolution on the most used morphological parameters for classifying galaxies through non-parametric methods: Abraham and Conselice-Bershady concentration indices, Gini, M20 moment of light, asymmetry, and smoothness. Three different non-local datasets are used, ALHAMBRA and SXDS (examples of deep ground-based surveys), and COSMOS (deep space-based survey). We used a sample of 3000 local, visually classified galaxies, measuring their morphological parameters at their real redshifts (z ~ 0). Then we simulated them to match the redshift and magnitude distributions of galaxies in the non-local surveys. The comparisons of the two sets allow to put constraints on the use of each parameter for morphological classification and evaluate the effectiveness of the commonly used morphological diagnostic diagrams. All analysed parameters suffer from biases related to spatial resolution and depth, the impact of the former being much stronger. When including asymmetry and smoothness in classification diagrams, the noise effects must be taken into account carefully, especially for ground-based surveys. M20 is significantly affected, changing both the shape and range of its distribution at all brightness levels.We suggest that diagnostic diagrams based on 2 - 3 parameters should be avoided when classifying galaxies in ground-based surveys, independently of their brightness; for COSMOS they should be avoided for galaxies fainter than F814 = 23.0. These results can be applied directly to surveys similar to ALHAMBRA, SXDS and COSMOS, and also can serve as an upper/lower limit for shallower/deeper ones.

preprint2014arXiv

An 8-mm diameter Fiber Robot Positioner for Massive Spectroscopy Surveys

Massive spectroscopic survey are becoming trendy in astrophysics and cosmology, as they can address new fundamental knowledge such as Galactic Archaeology and probe the nature of the mysterious Dark Energy. To enable massive spectroscopic surveys, new technology are being developed to place thousands of optical fibers at a given position on a focal plane. These technology needs to be: 1) accurate, with micrometer positional accuracy; 2) fast to minimize overhead; 3) robust to minimize failure; and 4) low cost. In this paper we present the development of a new 8-mm in diameter fiber positionner robot using two 4mm DC-brushless gearmotors, developed in the context of the Dark Energy Spectroscopic Instrument. This development was conducted by a Spanish-Swiss (ES-CH) team led by the Instituto de Física Teórica (UAM-CSIC) and the Laboratoire d'Astrophysique (EPFL), in collaboration with the AVS company in Spain and the Faulhaber group (MPS & FAULHABER-MINIMOTOR) in Switzerland.

preprint2014arXiv

Contemporaneous observations of the radio galaxy NGC 1275 from radio to very high energy gamma-rays

The radio galaxy NGC 1275, recently identified as a very high energy (VHE, >100 GeV) gamma-ray emitter by MAGIC, is one of the few non-blazar AGN detected in the VHE regime. In order to better understand the origin of the gamma-ray emission and locate it within the galaxy, we studied contemporaneous multi-frequency observations of NGC 1275 and modeled the overall spectral energy distribution (SED). We analyzed unpublished MAGIC observations carried out between Oct. 2009 and Feb. 2010, and the previously published ones taken between Aug. 2010 and Feb. 2011. We studied the multi-band variability and correlations by analyzing data of Fermi-LAT (0.1-100 GeV), as well as Chandra (X-ray), KVA (optical) and MOJAVE (radio) data taken during the same period. Using customized Monte Carlo simulations corresponding to early MAGIC stereo data, we detect NGC 1275 also in the earlier campaign. The flux level and energy spectra are similar to the results of the second campaign. The monthly light curve >100 GeV shows a hint of variability at the 3.6 sigma level. In the Fermi-LAT band, both flux and spectral shape variabilities are reported. The optical light curve is variable and shows a clear correlation with the gamma-ray flux >100 MeV. In radio, 3 compact components are resolved in the innermost part of the jet. One of them shows a similar trend as the LAT and KVA light curves. The 0.1-650 GeV spectra measured simultaneously with MAGIC and Fermi-LAT can be well fitted either by a log-parabola or by a power-law with a sub-exponential cutoff for both campaigns. A single-zone synchrotron-self-Compton model, with an electron spectrum following a power-law with an exponential cutoff, can explain the broadband SED and the multi-band behavior of the source. However, this model suggests an untypical low bulk-Lorentz factor or a velocity alignment closer to the line of sight than the pc-scale radio jet.

preprint2014arXiv

Discovery of very high energy gamma-ray emission from the blazar 1ES 0033+595 by the MAGIC telescopes

The number of known very high energy (VHE) blazars is $\sim\,50$, which is very small in comparison to the number of blazars detected in other frequencies. This situation is a handicap for population studies of blazars, which emit about half of their luminosity in the $γ$-ray domain. Moreover, VHE blazars, if distant, allow for the study of the environment that the high-energy $γ$-rays traverse in their path towards the Earth, like the extragalactic background light (EBL) and the intergalactic magnetic field (IGMF), and hence they have a special interest for the astrophysics community. We present the first VHE detection of 1ES\,0033+595 with a statistical significance of 5.5\,$σ$. The VHE emission of this object is constant throughout the MAGIC observations (2009 August and October), and can be parameterized with a power law with an integral flux above 150 GeV of $(7.1\pm1.3)\times 10^{-12} {\mathrm{ph\,cm^{-2}\,s^{-1}}}$ and a photon index of ($3.8\pm0.7$). We model its spectral energy distribution (SED) as the result of inverse Compton scattering of synchrotron photons. For the study of the SED we used simultaneous optical R-band data from the KVA telescope, archival X-ray data by \textit{Swift} as well as \textit{INTEGRAL}, and simultaneous high energy (HE, $300$\,MeV~--~$10$\,GeV) $γ$-ray data from the \textit{Fermi} LAT observatory. Using the empirical approach of Prandini et al. (2010) and the \textit{Fermi}-LAT and MAGIC spectra for this object, we estimate the redshift of this source to be $0.34\pm0.08\pm0.05$. This is a relevant result because this source is possibly one of the ten most distant VHE blazars known to date, and with further (simultaneous) observations could play an important role in blazar population studies, as well as future constraints on the EBL and IGMF.

preprint2014arXiv

First broadband characterization and redshift determination of the VHE blazar MAGIC J2001+439

We aim to characterize the broadband emission from 2FGL J2001.1+4352, which has been associated with the unknown-redshift blazar MG4 J200112+4352. Based on its gamma-ray spectral properties, it was identified as a potential very high energy (VHE; E > 100 GeV) gamma-ray emitter. The source was observed with MAGIC first in 2009 and later in 2010 within a multi-instrument observation campaign. The MAGIC observations yielded 14.8 hours of good quality stereoscopic data. The object was monitored at radio, optical and gamma-ray energies during the years 2010 and 2011. The source, named MAGIC J2001+439, is detected for the first time at VHE with MAGIC at a statistical significance of 6.3 σ (E > 70 GeV) during a 1.3-hour long observation on 2010 July 16. The multi-instrument observations show variability in all energy bands with the highest amplitude of variability in the X-ray and VHE bands. We also organized deep imaging optical observations with the Nordic Optical Telescope in 2013 to determine the source redshift. We determine for the first time the redshift of this BL Lac object through the measurement of its host galaxy during low blazar activity. Using the observational evidence that the luminosities of BL Lac host galaxies are confined to a relatively narrow range, we obtain z = 0.18 +/- 0.04. Additionally, we use the Fermi-LAT and MAGIC gamma-ray spectra to provide an independent redshift estimation, z = 0.17 +/- 0.10. Using the former (more accurate) redshift value, we adequately describe the broadband emission with a one-zone SSC model for different activity states and interpret the few-day timescale variability as produced by changes in the high-energy component of the electron energy distribution.

preprint2014arXiv

MAGIC gamma-ray and multifrequency observations of flat spectrum radio quasar PKS 1510-089 in early 2012

Among more than fifty blazars detected in very high energy (VHE, E>100GeV) gamma-rays, only three belong to the subclass of Flat Spectrum Radio Quasars (FSRQs). MAGIC observed FSRQ PKS 1510-089 in February-April 2012 during a high activity state in the high energy (HE, E>100 MeV) gamma-ray band observed by AGILE and Fermi. MAGIC observations result in the detection of a source with significance of 6.0 sigma. In agreement with the previous VHE observations of the source, we find no statistically significant variability during the MAGIC observations in daily, weekly or monthly time scales. The other two known VHE FSRQs have shown daily scale to sub-hour variability. We study the multifrequency behaviour of the source at the epoch of MAGIC observation, collecting quasi-simultaneous data at radio and optical (GASP-WEBT and F-Gamma collaborations, REM, Steward, Perkins, Liverpool, OVRO and VLBA telescopes), X-ray (Swift satellite) and HE gamma-ray frequencies. The gamma-ray SED combining AGILE, Fermi and MAGIC data joins smoothly and shows no hint of a break. The multifrequency light curves suggest a common origin for the millimeter radio and HE gamma-ray emission and the HE gamma-ray flaring starts when the new component is ejected from the 43GHz VLBA core. The quasi-simultaneous multifrequency SED is modelled with a one-zone inverse Compton model. We study two different origins of the seed photons for the inverse Compton scattering, namely the infra-red torus and a slow sheath surrounding the jet around the VLBA core. Both models fit the data well. However, the fast HE gamma-ray variability requires that within the modelled large emitting region, there must exist more compact regions. We suggest that these observed signatures would be most naturally explained by a turbulent plasma flowing at a relativistic speed down the jet and crossing a standing conical shock.

preprint2014arXiv

MAGIC long-term study of the distant TeV blazar PKS 1424+240 in a multiwavelength context

We present a study of the very high energy (VHE; E > 100 GeV) gamma-ray emission of the blazar PKS 1424+240 observed with the MAGIC telescopes. The primary aim of this paper is the multiwavelength spectral characterization and modeling of this blazar, which is made particularly interesting by the recent discovery of a lower limit of its redshift of z > 0.6 and makes it a promising candidate to be the most distant VHE source. The source has been observed with the MAGIC telescopes in VHE gamma rays for a total observation time of ~33.6 h from 2009 to 2011. The source was marginally detected in VHE gamma rays during 2009 and 2010, and later, the detection was confirmed during an optical outburst in 2011. The combined significance of the stacked sample is ~7.2 sigma. The differential spectra measured during the different campaigns can be described by steep power laws with the indices ranging from 3.5 +/- 1.2 to 5.0 +/- 1.7. The MAGIC spectra corrected for the absorption due to the extragalactic background light connect smoothly, within systematic errors, with the mean spectrum in 2009-2011 observed at lower energies by the Fermi-LAT. The absorption-corrected MAGIC spectrum is flat with no apparent turn down up to 400 GeV. The multiwavelength light curve shows increasing flux in radio and optical bands that could point to a common origin from the same region of the jet. The large separation between the two peaks of the constructed non-simultaneous spectral energy distribution also requires an extremely high Doppler factor if an one zone synchrotron self-Compton model is applied. We find that a two-component synchrotron self-Compton model describes the spectral energy distribution of the source well, if the source is located at z~0.6.

preprint2014arXiv

MAGIC reveals a complex morphology within the unidentified gamma-ray source HESS J1857+026

HESS J1857+026 is an extended TeV gamma-ray source that was discovered by H.E.S.S. as part of its Galactic plane survey. Given its broadband spectral energy distribution and its spatial coincidence with the young energetic pulsar PSR J1856+0245, the source has been put forward as a pulsar wind nebula (PWN) candidate. MAGIC has performed follow-up observations aimed at mapping the source down to energies approaching 100 GeV in order to better understand its complex morphology. HESS J1857+026 was observed by MAGIC in 2010, yielding 29 hours of good quality stereoscopic data that allowed us to map the source region in two separate ranges of energy. We detected very-high-energy gamma-ray emission from HESS J1857+026 with a significance of $12 σ$ above $150$ GeV. The differential energy spectrum between $100$ GeV and $13$ TeV is well described by a power law function $dN/dE = N_0(E/1\textrm{TeV})^{-Γ}$ with $N_0 = (5.37 \pm0.44_{stat} \pm1.5_{sys}) \times 10^{-12} (\textrm{TeV}^{-1} \textrm{cm}^{-2}$ $\textrm{ s}^{-1})$ and $Γ= 2.16\pm0.07_{stat} \pm0.15_{sys}$, which bridges the gap between the GeV emission measured by Fermi-LAT and the multi-TeV emission measured by H.E.S.S.. In addition, we present a detailed analysis of the energy-dependent morphology of this region. We couple these results with archival multi-wavelength data and outline evidence in favor of a two-source scenario, whereby one source is associated with a PWN, while the other could be linked with a molecular cloud complex containing an HII region and a possible gas cavity.

preprint2014arXiv

Multifrequency Studies of the Peculiar Quasar 4C +21.35 During the 2010 Flaring Activity

The discovery of rapidly variable Very High Energy (VHE; E > 100 GeV) gamma-ray emission from 4C +21.35 (PKS 1222+216) by MAGIC on 2010 June 17, triggered by the high activity detected by the Fermi Large Area Telescope (LAT) in high energy (HE; E > 100 MeV) gamma-rays, poses intriguing questions on the location of the gamma-ray emitting region in this flat spectrum radio quasar. We present multifrequency data of 4C +21.35 collected from centimeter to VHE during 2010 to investigate the properties of this source and discuss a possible emission model. The first hint of detection at VHE was observed by MAGIC on 2010 May 3, soon after a gamma-ray flare detected by Fermi-LAT that peaked on April 29. The same emission mechanism may therefore be responsible for both the HE and VHE emission during the 2010 flaring episodes. Two optical peaks were detected on 2010 April 20 and June 30, close in time but not simultaneous with the two gamma-ray peaks, while no clear connection was observed between the X-ray an gamma-ray emission. An increasing flux density was observed in radio and mm bands from the beginning of 2009, in accordance with the increasing gamma-ray activity observed by Fermi-LAT, and peaking on 2011 January 27 in the mm regime (230 GHz). We model the spectral energy distributions (SEDs) of 4C +21.35 for the two periods of the VHE detection and a quiescent state, using a one-zone model with the emission coming from a very compact region outside the broad line region. The three SEDs can be fit with a combination of synchrotron self-Compton and external Compton emission of seed photons from a dust torus, changing only the electron distribution parameters between the epochs. The fit of the optical/UV part of the spectrum for 2010 April 29 seems to favor an inner disk radius of <6 gravitational radii, as one would expect from a prograde-rotating Kerr black hole.

preprint2014arXiv

Multiwavelength observations of Mrk 501 in 2008

Mrk 501 is one of the brightest blazars at TeV energies and has been extensively studied since its first VHE detection in 1996. Our goal is to characterize in detail the source gamma-ray emission, together with the radio-to-X-ray emission, during the non-flaring (low) activity, which is less often studied than the occasional flaring (high) activity. We organized a multiwavelength (MW) campaign on Mrk 501 between March and May 2008. This multi-instrument effort included the most sensitive VHE gamma-ray instruments in the northern hemisphere, namely the imaging atmospheric Cherenkov telescopes MAGIC and VERITAS, as well as Swift, RXTE, the F-GAMMA, GASP-WEBT, and other collaborations and instruments. Mrk 501 was found to be in a low state of activity during the campaign, with a VHE flux in the range of 10%-20% of the Crab nebula flux. Nevertheless, significant flux variations were detected with various instruments, with a trend of increasing variability with energy. The broadband spectral energy distribution during the two different emission states of the campaign can be adequately described within the homogeneous one-zone synchrotron self-Compton model, with the (slightly) higher state described by an increase in the electron number density. This agrees with previous studies of the broadband emission of this source during flaring and non-flaring states. We report for the first time a tentative X-ray-to-VHE correlation during a low VHE activity. Although marginally significant, this positive correlation between X-ray and VHE, which has been reported many times during flaring activity, suggests that the mechanisms that dominate the X-ray/VHE emission during non-flaring-activity are not substantially different from those that are responsible for the emission during flaring activity.

preprint2014arXiv

Optimized dark matter searches in deep observations of Segue 1 with MAGIC

We present the results of stereoscopic observations of the satellite galaxy Segue 1 with the MAGIC Telescopes, carried out between 2011 and 2013. With almost 160 hours of good-quality data, this is the deepest observational campaign on any dwarf galaxy performed so far in the very high energy range of the electromagnetic spectrum. We search this large data sample for signals of dark matter particles in the mass range between 100 GeV and 20 TeV. For this we use the full likelihood analysis method, which provides optimal sensitivity to characteristic gamma-ray spectral features, like those expected from dark matter annihilation or decay. In particular, we focus our search on gamma-rays produced from different final state Standard Model particles, annihilation with internal bremsstrahlung, monochromatic lines and box-shaped signals. Our results represent the most stringent constraints to the annihilation cross-section or decay lifetime obtained from observations of satellite galaxies, for masses above few hundred GeV. In particular, our strongest limit (95% confidence level) corresponds to a ~500 GeV dark matter particle annihilating into tau+tau-, and is of order <sigma v> ~ 1.2x10^{-24} cm^3 s^{-1} - a factor ~40 above the <sigma v> thermal value.

preprint2014arXiv

Rapid and multi-band variability of the TeV-bright active nucleus of the galaxy IC 310

Context. The radio galaxy IC 310 has recently been identified as a gamma-ray emitter based on observations at GeV energies with Fermi-LAT and at very high energies (VHE, E>100GeV) with the MAGIC telescopes. Originally classified as a head-tail radio galaxy, the nature of this object is subject of controversy since its nucleus shows blazar-like behavior. Aims. In order to understand the nature of IC 310 and the origin of the VHE emission we studied the spectral and flux variability of IC 310 from the X-ray band to the VHE gamma-ray regime. Methods. The light curve of IC 310 above 300GeV has been measured with the MAGIC telescopes from Oct. 2009 to Feb. 2010. Fermi-LAT data (2008-2011) in the 10-500GeV energy range were also analyzed. In X-ray, archival observations from 2003 to 2007 with XMM, Chandra, and Swift-XRT in the 0.5-10keV band were studied. Results. The VHE light curve reveals several high-amplitude and short-duration flares. Day-to-day flux variability is clearly present. The photon index between 120GeV and 8TeV remains at the value $Γ\sim2.0$ during both low and high flux states. The VHE spectral shape does not show significant variability, whereas the flux at 1TeV changes by a factor of $\sim7$. Fermi-LAT detected only eight gamma-ray events in the energy range 10GeV-500GeV in three years of observation. The measured photon index of $Γ=1.3\pm0.5$ in the Fermi-LAT range is very hard. The X-ray measurements show strong variability in flux and photon index. The latter varied from $1.76\pm0.07$ to $2.55\pm0.07$. Conclusion. The rapid variability measured confirms the blazar-like behavior of IC 310. The TeV emission seems to originate from scales of less than 80 Schwarzschild radii within the compact core of its FRI radio jet with orientation angle 10deg-38deg. The SED resembles that of an extreme blazar, albeit the luminosity is more than two orders of magnitude lower.

preprint2014arXiv

Search for Very-High-Energy Gamma Rays from the z = 0.896 Quasar 4C +55.17 with the MAGIC telescopes

The bright gamma-ray quasar 4C +55.17 is a distant source ($z = 0.896$) with a hard spectrum at GeV energies as observed by the Large Area Telescope (LAT) on board the {\it Fermi} satellite. This source is identified as a good source candidate for very-high-energy (VHE; $> 30$ GeV) gamma rays. In general VHE gamma rays from distant sources provide an unique opportunity to study the extragalactic background light (EBL) and underlying astrophysics. The flux intensity of this source in the VHE range is investigated. Then, constraints on the EBL are derived from the attenuation of gamma-ray photons coming from the distant blazar. We searched for a gamma-ray signal from this object using the 35-hour observations taken by the MAGIC telescopes between November 2010 and January 2011. No significant VHE gamma-ray signal was detected. We computed the upper limits of the integrated gamma-ray flux at $95\%$ confidence level of $9.4 \times 10^{-12}$ cm$^{-2}$ s$^{-1}$ and $2.5 \times 10^{-12}$ cm$^{-2}$ s$^{-1}$ above $100$ GeV and $200$ GeV, respectively. The differential upper limits in four energy bins in the range from $80$ GeV to $500$ GeV are also derived. The upper limits are consistent with the attenuation predicted by low-flux EBL models on the assumption of a simple power-law spectrum extrapolated from LAT data.

preprint2014arXiv

The ALHAMBRA survey: Accurate merger fractions by PDF analysis of photometric close pairs

Our goal is to develop and test a novel methodology to compute accurate close pair fractions with photometric redshifts. We improve the current methodologies to estimate the merger fraction f_m from photometric redshifts by (i) using the full probability distribution functions (PDFs) of the sources in redshift space, (ii) including the variation in the luminosity of the sources with z in both the selection of the samples and in the luminosity ratio constrain, and (iii) splitting individual PDFs into red and blue spectral templates to deal robustly with colour selections. We test the performance of our new methodology with the PDFs provided by the ALHAMBRA photometric survey. The merger fractions and rates from the ALHAMBRA survey are in excellent agreement with those from spectroscopic work, both for the general population and for red and blue galaxies. With the merger rate of bright (M_B <= -20 - 1.1z) galaxies evolving as (1+z)^n, the power-law index n is larger for blue galaxies (n = 2.7 +- 0.5) than for red galaxies (n = 1.3 +- 0.4), confirming previous results. Integrating the merger rate over cosmic time, we find that the average number of mergers per galaxy since z = 1 is N_m = 0.57 +- 0.05 for red galaxies and N_m = 0.26 +- 0.02 for blue galaxies. Our new methodology exploits statistically all the available information provided by photometric redshift codes and provides accurate measurements of the merger fraction by close pairs only using photometric redshifts. Current and future photometric surveys will benefit of this new methodology.

preprint2014arXiv

The ALHAMBRA survey: an empirical estimation of the cosmic variance for merger fraction studies based on close pairs

Our goal is to estimate empirically, for the first time, the cosmic variance that affects merger fraction studies based on close pairs. We compute the merger fraction from photometric redshift close pairs with 10h^-1 kpc <= rp <= 50h^-1 kpc and Dv <= 500 km/s, and measure it in the 48 sub-fields of the ALHAMBRA survey. We study the distribution of the measured merger fractions, that follow a log-normal function, and estimate the cosmic variance sigma_v as the intrinsic dispersion of the observed distribution. We develop a maximum likelihood estimator to measure a reliable sigma_v and avoid the dispersion due to the observational errors (including the Poisson shot noise term). The cosmic variance of the merger fraction depends mainly on (i) the number density of the populations under study, both for the principal (n_1) and the companion (n_2) galaxy in the close pair, and (ii) the probed cosmic volume V_c. We find a significant dependence on neither the search radius used to define close companions, the redshift, nor the physical selection (luminosity or stellar mass) of the samples. We provide a parametrisation of the cosmic variance with n_1, n_2, and V_c, sigma_v = 0.48 n_1^{-0.54} V_c^{-0.48} (n_2/n_1)^{-0.37}. Thanks to this prescription, future merger fraction studies based on close pairs could account properly for the cosmic variance on their results.

preprint2014arXiv

The ALHAMBRA survey: evolution of galaxy clustering since $z \sim 1$

We study the clustering of galaxies as function of luminosity and redshift in the range $0.35 < z < 1.25$ using data from the Advanced Large Homogeneous Area Medium Band Redshift Astronomical (ALHAMBRA) survey. The ALHAMBRA data used in this work cover $2.38 \mathrm{deg}^2$ in 7 independent fields, after applying a detailed angular selection mask, with accurate photometric redshifts, $σ_z \lesssim 0.014 (1+z)$, down to $I_{\rm AB} < 24$. Given the depth of the survey, we select samples in $B$-band luminosity down to $L^{\rm th} \simeq 0.16 L^{*}$ at $z = 0.9$. We measure the real-space clustering using the projected correlation function, accounting for photometric redshifts uncertainties. We infer the galaxy bias, and study its evolution with luminosity. We study the effect of sample variance, and confirm earlier results that the COSMOS and ELAIS-N1 fields are dominated by the presence of large structures. For the intermediate and bright samples, $L^{\rm med} \gtrsim 0.6L^{*}$, we obtain a strong dependence of bias on luminosity, in agreement with previous results at similar redshift. We are able to extend this study to fainter luminosities, where we obtain an almost flat relation, similar to that observed at low redshift. Regarding the evolution of bias with redshift, our results suggest that the different galaxy populations studied reside in haloes covering a range in mass between $\log_{10}[M_{\rm h}/(h^{-1}\mathrm{M}_{\odot})] \gtrsim 11.5$ for samples with $L^{\rm med} \simeq 0.3 L^{*}$ and $\log_{10}[M_{\rm h}/(h^{-1}\mathrm{M}_{\odot})] \gtrsim 13.0$ for samples with $L^{\rm med} \simeq 2 L^{*}$, with typical occupation numbers in the range of $\sim 1 - 3$ galaxies per halo.

preprint2013arXiv

Detection of the cosmic γ-ray horizon from multiwavelength observations of blazars

The first statistically significant detection of the cosmic γ-ray horizon (CGRH) that is independent of any extragalactic background light (EBL) model is presented. The CGRH is a fundamental quantity in cosmology. It gives an estimate of the opacity of the Universe to very high energy (VHE) γ-ray photons due to photon-photon pair production with the EBL. The only estimations of the CGRH to date are predictions from EBL models and lower limits from γ-ray observations of cosmological blazars and γ-ray bursts. Here, we present homogeneous synchrotron/synchrotron self-Compton (SSC) models of the spectral energy distributions of 15 blazars based on (almost) simultaneous observations from radio up to the highest energy γ-rays taken with the Fermi satellite. These synchrotron/SSC models predict the unattenuated VHE fluxes, which are compared with the observations by imaging atmospheric Cherenkov telescopes. This comparison provides an estimate of the optical depth of the EBL, which allows a derivation of the CGRH through a maximum likelihood analysis that is EBL-model independent. We find that the observed CGRH is compatible with the current knowledge of the EBL.

preprint2013arXiv

Discovery of very high energy gamma-ray emission from the blazar 1ES 1727+502 with the MAGIC Telescopes

Motivated by the Costamante & Ghisellini (2002) predictions we investigated if the blazar 1ES 1727+502 (z=0.055) is emitting very high energy (VHE, E>100 GeV) gamma rays. We observed the BL Lac object 1ES 1727+502 in stereoscopic mode with the two MAGIC telescopes during 14 nights between May 6th and June 10th 2011, for a total effective observing time of 12.6 hours. For the study of the multiwavelength spectral energy distribution (SED) we use simultaneous optical R-band data from the KVA telescope, archival UV/optical and X-ray observations by instruments UVOT and XRT on board of the Swift satellite and high energy (HE, 0.1 GeV - 100 GeV) gamma-ray data from the Fermi-LAT instrument. We detect, for the first time, VHE gamma-ray emission from 1ES 1727+502 at a statistical significance of 5.5 sigma. The integral flux above 150 GeV is estimated to be (2.1\pm0.4)% of the Crab Nebula flux and the de-absorbed VHE spectrum has a photon index of (2.7\pm0.5). No significant short-term variability was found in any of the wavebands presented here. We model the SED using a one-zone synchrotron self-Compton model obtaining parameters typical for this class of sources.

preprint2013arXiv

Lyman break and UV-selected galaxies at $z \sim 1$ I. Stellar populations from ALHAMBRA survey

We take advantage of the exceptional photometric coverage provided by the combination of GALEX data in the UV and the ALHAMBRA survey in the optical and near-IR to analyze the physical properties of a sample of 1225 GALEX-selected Lyman break galaxies (LBGs) at $0.8 \lesssim z \lesssim 1.2$ located in the COSMOS field. This is the largest sample of LBGs studied at that redshift range so far. According to a spectral energy distribution (SED) fitting with synthetic stellar population templates, we find that LBGs at $z \sim 1$ are mostly young galaxies with a median age of 341 Myr and have intermediate dust attenuation, $\ < E_s (B-V) \ > \sim 0.20$. Due to their selection criterion, LBGs at $z \sim 1$ are UV-bright galaxies and have high dust-corrected total SFR, with a median value of 16.9 $M_\odot {\rm yr}^{-1}$. Their median stellar mass is $\log{\left(M_*/M_\odot \right)} = 9.74$. We obtain that the dust-corrected total SFR of LBGs increases with stellar mass and the specific SFR is lower for more massive galaxies. Only 2% of the galaxies selected through the Lyman break criterion have an AGN nature. LBGs at $z \sim 1$ are mostly located over the blue cloud of the color-magnitude diagram of galaxies at their redshift, with only the oldest and/or the dustiest deviating towards the green valley and red sequence. Morphologically, 69% of LBGs are disk-like galaxies, with the fraction of interacting, compact, or irregular systems being much lower, below 12%. LBGs have a median effective radius of 2.5 kpc and bigger galaxies have higher total SFR and stellar mass. Comparing to their high-redshift analogues, we find evidence that LBGs at lower redshifts are bigger, redder in the UV continuum, and have a major presence of older stellar populations in their SEDs. However, we do not find significant difference in the distributions of stellar mass or dust attenuation.

preprint2013arXiv

Lyman break and UV-selected galaxies at z ~ 1: II. PACS-100um/160um FIR detections

We report the PACS-100um/160um detections of a sample of 42 GALEX-selected and FIR-detected Lyman break galaxies (LBGs) at z ~ 1 located in the COSMOS field and analyze their ultra-violet (UV) to far-infrared (FIR) properties. The detection of these LBGs in the FIR indicates that they have a dust content high enough so that its emission can be directly detected. According to a spectral energy distribution (SED) fitting with stellar population templates to their UV-to-near-IR observed photometry, PACS-detected LBGs tend to be bigger, more massive, dustier, redder in the UV continuum, and UV-brighter than PACS-undetected LBGs. PACS-detected LBGs at z ~ 1 are mostly disk-like galaxies and are located over the green-valley and red sequence of the color-magnitude diagram of galaxies at their redshift. By using their UV and IR emission, we find that PACS-detected LBGs tend to be less dusty and have slightly higher total star-formation rates (SFRs) than other PACS-detected UV-selected galaxies within their same redshift range. As a consequence of the selection effect due to the depth of the FIR observations employed, all our PACS-detected LBGs are LIRGs. However, none of them are in the ULIRG regime, where the FIR observations are complete. The finding of ULIRGs-LBGs at higher redshifts suggests an evolution of the FIR emission of LBGs with cosmic time. In an IRX-$β$ diagram, PACS-detected LBGs at z ~ 1 tend to be located around the relation for local starburst similarly to other UV-selected PACS-detected galaxies at their same redshift. Consequently, the dust-correction factors obtained with their UV continuum slope allow to determine their total SFR, unlike at higher redshifts. However, the dust attenuation derived from UV to NIR SED fitting overestimates the total SFR for most of our PACS-detected LBGs in age-dependent way: the overestimation factor is higher in younger galaxies.

preprint2013arXiv

The ALHAMBRA survey: Discovery of a faint QSO at z = 5.41

We aim to illustrate the potentiality of the Advanced Large, Homogeneous Area, Medium-Band Redshift Astronomical (ALHAMBRA) survey to investigate the high redshift universe through the detection of quasi stellar objects (QSOs) at redshifts larger than 5. The search for z>5 QSOs candidates was done by fitting an extensive library of spectral energy distributions --including active and non-active galaxy templates as well as stars-- to the photometric database of the ALHAMBRA survey (composed of 20 optical medium-band plus the 3 broad-band JHKs filters). Our selection over ~1 square degree of ALHAMBRA data (~1/4 of the total area covered by the survey), combined with GTC/OSIRIS spectroscopy, has yielded the identification of an optically faint QSO at very high redshift (z = 5.41). The QSO has an absolute magnitude of ~-24 at the 1450Å continuum, a bolometric luminosity of ~2x10^46 erg/s and an estimated black hole mass of ~10^8 Msolar. This QSO adds itself to a reduced number of known UV faint sources at these redshifts. The preliminary derived space density is compatible with the most recent determinations of the high-z QSO luminosity functions (QLF). This new detection shows how ALHAMBRA, as well as forthcoming well designed photometric surveys, can provide a wealth of information on the origin and early evolution of this kind of objects.

preprint2013arXiv

The ALHAMBRA survey: reliable morphological catalogue of 22,051 early- and late-type galaxies

ALHAMBRA is a photometric survey designed to trace the cosmic evolution and cosmic variance. It covers a large area of ~ 4 sq. deg in 8 fields, where 7 fields overlap with other surveys, allowing to have complementary data in other wavelengths. All observations were carried out in 20 continuous, medium band (30 nm width) optical and 3 near-infrared (JHK) bands, providing the precise measurements of photometric redshifts. In addition, morphological classification of galaxies is crucial for any kind of galaxy formation and cosmic evolution studies, providing the information about star formation histories, their environment and interactions, internal perturbations, etc. We present a morphological classification of > 40,000 galaxies in the ALHAMBRA survey. We associate to every galaxy a probability to be early-type using the automated Bayesian code galSVM. Despite of the spatial resolution of the ALHAMBRA images (~ 1 arcsec), for 22,051 galaxies we obtained the contamination by other type of less than 10%. Of those, 1,640 and 10,322 galaxies are classified as early- (down to redshifts ~ 0.5) and late-type (down to redshifts ~ 1.0), respectively, with magnitudes F613W < 22.0. In addition, for magnitude range 22.0 < F613W < 23.0 we classified other 10,089 late-type galaxies with redshifts < 1.3. We show that the classified objects populate the expected regions in the colour-mass and colour-magnitude planes. The presented dataset is especially attractive given the homogeneous multi-wavelength coverage available in the ALHAMBRA fields, and is intended to be used in a variety of scientific applications. The low-contamination catalogue (< 10%) is made publicly available with the present paper.

preprint2012arXiv

CLUES on Fermi-LAT prospects for the extragalactic detection of munuSSM gravitino Dark Matter

The munuSSM is a supersymmetric model that has been proposed to solve the problems generated by other supersymmetric extensions of the standard model of particle physics. Given that R-parity is broken in the munuSSM, the gravitino is a natural candidate for decaying dark matter since its lifetime becomes much longer than the age of the Universe. In this model, gravitino dark matter could be detectable through the emission of a monochromatic gamma ray in a two-body decay. We study the prospects of the Fermi-LAT telescope to detect such monochromatic lines in 5 years of observations of the most massive nearby extragalactic objects. The dark matter halo around the Virgo galaxy cluster is selected as a reference case, since it is associated to a particularly high signal-to-noise ratio and is located in a region scarcely affected by the astrophysical diffuse emission from the galactic plane. The simulation of both signal and background gamma-ray events is carried out with the Fermi Science Tools, and the dark matter distribution around Virgo is taken from a N-body simulation of the nearby extragalactic Universe, with constrained initial conditions provided by the CLUES project. We find that a gravitino with a mass range of 0.6 to 2 GeV, and with a lifetime range of about 3x10^27 to 2x10^28 s would be detectable by the Fermi-LAT with a signal-to-noise ratio larger than 3. We also obtain that gravitino masses larger than about 4 GeV are already excluded in the munuSSM by Fermi-LAT data of the galactic halo

preprint2012arXiv

Constraining Cosmic Rays and Magnetic Fields in the Perseus Galaxy Cluster with TeV observations by the MAGIC telescopes

Galaxy clusters are being assembled today in the most energetic phase of hierarchical structure formation which manifests itself in powerful shocks that contribute to a substantial energy density of cosmic rays (CRs). Hence, clusters are expected to be luminous gamma-ray emitters since they also act as energy reservoirs for additional CR sources, such as active galactic nuclei and supernova-driven galactic winds. To detect the gamma-ray emission from CR interactions with the ambient cluster gas, we conducted the deepest to date observational campaign targeting a galaxy cluster at very high-energy gamma-rays and observed the Perseus cluster with the MAGIC Cherenkov telescopes for a total of ~85 hr of effective observing time. This campaign resulted in the detection of the central radio galaxy NGC 1275 at energies E > 100 GeV with a very steep energy spectrum. Here, we restrict our analysis to energies E > 630 GeV and detect no significant gamma-ray excess. This constrains the average CR-to-thermal pressure ratio to be <= 1-2%, depending on assumptions and the model for CR emission. Comparing these gamma-ray upper limits to predictions from cosmological cluster simulations that include CRs constrains the maximum CR acceleration efficiency at structure formation shocks to be < 50%. Alternatively, this may argue for non-negligible CR transport processes such as CR streaming and diffusion into the outer cluster regions. Finally, we derive lower limits on the magnetic field distribution assuming that the Perseus radio mini-halo is generated by secondary electrons/positrons that are created in hadronic CR interactions: assuming a spectrum of E^-2.2 around TeV energies as implied by cluster simulations, we limit the central magnetic field to be > 4-9 microG, depending on the rate of decline of the magnetic field strength toward larger radii.

preprint2012arXiv

Detection of very-high energy γ-ray emission from NGC 1275 by the MAGIC telescopes

We report on the detection of very-high energy (VHE, E>100 GeV) gamma-ray emission from NGC 1275, the central radio galaxy of the Perseus cluster of galaxies. The source has been detected by the MAGIC telescopes with a statistical significance of 6.6 sigma above 100 GeV in 46 hr of stereo observations carried out between August 2010 and February 2011. The measured differential energy spectrum between 70 GeV and 500 GeV can be described by a power law with a steep spectral index of Γ=-4.1+/-0.7stat+/-0.3syst, and the average flux above 100 GeV is F_{gamma}=(1.3+/-0.2stat+/-0.3syst) x 10^-11 cm^-2 s^-1. These results, combined with the power-law spectrum measured in the first two years of observations by the Fermi-LAT above 100 MeV, with a spectral index of Gamma ~ -2.1, strongly suggest the presence of a break or cut-off around tens of GeV in the NGC 1275 spectrum. The light curve of the source above 100 GeV does not show hints of variability on a month time scale. Finally, we report on the nondetection in the present data of the radio galaxy IC 310, previously discovered by the Fermi-LAT and MAGIC. The derived flux upper limit F^{U.L.}_{gamma} (>300 GeV)=1.2 x 10^-12 cm^-2 s^-1 is a factor ~ 3 lower than the mean flux measured by MAGIC between October 2009 and February 2010, thus confirming the year time-scale variability of the source at VHE.

preprint2012arXiv

Discovery of VHE gamma-ray emission from the BL Lac object B3 2247+381 with the MAGIC telescopes

We study the non-thermal jet emission of the BL Lac object B3 2247+381 during a high optical state. The MAGIC telescopes observed the source during 13 nights between September 30th and October 30th 2010, collecting a total of 14.2 hours of good quality very high energy (VHE) $γ$-ray data. Simultaneous multiwavelength data was obtained with X-ray observations by the Swift satellite and optical R-band observations at the KVA-telescope. We also use high energy $γ$-ray (HE, 0.1 GeV-100 GeV) data from the Fermi satellite. The BL Lac object B3 2247+381 (z=0.119) was detected, for the first time, at VHE $γ$-rays at a statistical significance of 5.6 $σ$. A soft VHE spectrum with a photon index of -3.2 $\pm$ 0.6 was determined. No significant short term flux variations were found. We model the spectral energy distribution using a one-zone SSC-model, which can successfully describe our data.

preprint2012arXiv

Discovery of VHE gamma-rays from the blazar 1ES 1215+303 with the MAGIC Telescopes and simultaneous multi-wavelength observations

Context. We present the discovery of very high energy (VHE, E > 100GeV) gamma-ray emission from the BL Lac object 1ES 1215+303 by the MAGIC telescopes and simultaneous multi-wavelength data in a broad energy range from radio to gamma-rays. Aims. We study the VHE gamma-ray emission from 1ES 1215+303 and its relation to the emissions in other wavelengths. Methods. Triggered by an optical outburst, MAGIC observed the source in January-February 2011 for 20.3 hrs. The target was monitored in the optical R-band by the KVA telescope that also performed optical polarization measurements. We triggered target of opportunity observations with the Swift satellite and obtained simultaneous and quasi-simultaneous data from the Fermi Large Area Telescope and from the Metsähovi radio telescope. We also present the analysis of older MAGIC data taken in 2010. Results. The MAGIC observations of 1ES 1215+303 carried out in January-February 2011 resulted in the first detection of the source at VHE with a statistical significance of 9.4 sigma. Simultaneously, the source was observed in a high optical and X-ray state. In 2010 the source was observed in a lower state in optical, X-ray, and VHE, while the GeV gamma-ray flux and the radio flux were comparable in 2010 and 2011. The spectral energy distribution obtained with the 2011 data can be modeled with a simple one zone SSC model, but it requires extreme values for the Doppler factor or the electron energy distribution.

preprint2012arXiv

High zenith angle observations of PKS 2155-304 with the MAGIC-I telescope

The high frequency peaked BL Lac PKS 2155-304 with a redshift of z=0.116 was discovered in 1997 in the very high energy (VHE, E >100GeV) gamma-ray range by the University of Durham Mark VI gamma-ray Cherenkov telescope in Australia with a flux corresponding to 20% of the Crab Nebula flux. It was later observed and detected with high significance by the Southern Cherenkov observatory H.E.S.S. Detection from the Northern hemisphere is difficult due to challenging observation conditions under large zenith angles. In July 2006, the H.E.S.S. collaboration reported an extraordinary outburst of VHE gamma-emission. During the outburst, the VHE gamma-ray emission was found to be variable on the time scales of minutes and with a mean flux of ~7 times the flux observed from the Crab Nebula. Follow-up observations with the MAGIC-I standalone Cherenkov telescope were triggered by this extraordinary outburst and PKS 2155-304 was observed between 28 July to 2 August 2006 for 15 hours at large zenith angles. Here we present our studies on the behavior of the source after its extraordinary flare and an enhanced analysis method for data taken at high zenith angles. We developed improved methods for event selection that led to a better background suppression. The averaged energy spectrum we derived has a spectral index of -3.5 +/- 0.2 above 400GeV, which is in good agreement with the spectral shape measured by H.E.S.S. during the major flare on MJD 53944. Furthermore, we present the spectral energy distribution modeling of PKS 2155-304. With our observations we increased the duty cycle of the source extending the light curve derived by H.E.S.S. after the outburst. Finally, we find night-by-night variability with a maximal amplitude of a factor three to four and an intranight variability in one of the nights (MJD 53945) with a similar amplitude.

preprint2012arXiv

MAGIC observations of the giant radio galaxy M87 in a low-emission state between 2005 and 2007

We present the results of a long M87 monitoring campaign in very high energy $γ$-rays with the MAGIC-I Cherenkov telescope. We aim to model the persistent non-thermal jet emission by monitoring and characterizing the very high energy $γ$-ray emission of M87 during a low state. A total of 150\,h of data were taken between 2005 and 2007 with the single MAGIC-I telescope, out of which 128.6\,h survived the data quality selection. We also collected data in the X-ray and \textit{Fermi}--LAT bands from the literature (partially contemporaneous). No flaring activity was found during the campaign. The source was found to be in a persistent low-emission state, which was at a confidence level of $7σ$. We present the spectrum between 100\,GeV and 2\,TeV, which is consistent with a simple power law with a photon index $Γ=2.21\pm0.21$ and a flux normalization at 300\,GeV of $(7.7\pm1.3) \times 10^{-8}$ TeV$^{-1}$ s$^{-1}$ m$^{-2}$. The extrapolation of the MAGIC spectrum into the GeV energy range matches the previously published \textit{Fermi}--LAT spectrum well, covering a combined energy range of four orders of magnitude with the same spectral index. We model the broad band energy spectrum with a spine layer model, which can satisfactorily describe our data.

preprint2012arXiv

Morphological and spectral properties of the W51 region measured with the MAGIC telescopes

The W51 complex hosts the supernova remnant W51C which is known to interact with the molecular clouds in the star forming region W51B. In addition, a possible pulsar wind nebula CXO J192318.5+140305 was found likely associated with the supernova remnant. Gamma-ray emission from this region was discovered by Fermi/LAT (between 0.2 and 50 GeV) and H.E.S.S. (>1 TeV). The spatial distribution of the events could not be used to pinpoint the location of the emission among the pulsar wind nebula, the supernova remnant shell and/or the molecular cloud. However, the modeling of the spectral energy distribution presented by the Fermi/LAT collaboration suggests a hadronic emission mechanism. We performed observations of the W51 complex with the MAGIC telescopes for more than 50 hours. The good angular resolution in the medium (few hundred GeV) to high (above 1 TeV) energies allow us to perform morphological studies. We detect an extended emission of very-high-energy gamma rays, with a significance of 11 standard deviations. We extend the spectrum from the highest Fermi/LAT energies to \sim 5 TeV and find that it follows a single power law with an index of 2.58 \pm 0.07stat \pm 0.22syst . The main part of the emission coincides with the shocked cloud region, while we find a feature extending towards the pulsar wind nebula. The possible contribution of the pulsar wind nebula, assuming a point-like source, shows no dependence on energy and it is about 20% of the overall emission. The broad band spectral energy distribution can be explained with a hadronic model that implies proton acceleration above 100 TeV. This result, together with the morphology of the source, tentatively suggests that we observe ongoing acceleration of ions in the interaction zone between supernova remnant and cloud. These results shed light on the long-standing problem of the origin of galactic cosmic rays.

preprint2012arXiv

Mrk 421 active state in 2008: the MAGIC view, simultaneous multi-wavelength observations and SSC model constrained

Context: The blazar Markarian 421 is one of the brightest TeV gamma-ray sources of the northern sky. From December 2007 until June 2008 it was intensively observed in the very high energy (VHE, E > 100 GeV) band by the single-dish Major Atmospheric Gamma-ray Imaging Cherenkov telescope (MAGIC-I). Aims: We aimed to measure the physical parameters of the emitting region of the blazar jet during active states. Methods: We performed a dense monitoring of the source in VHE with MAGIC-I, and also collected complementary data in soft X-rays and optical-UV bands; then, we modeled the spectral energy distributions (SED) derived from simultaneous multi-wavelength data within the synchrotron self--compton (SSC) framework. Results: The source showed intense and prolonged gamma-ray activity during the whole period, with integral fluxes (E > 200 GeV) seldom below the level of the Crab Nebula, and up to 3.6 times this value. Eight datasets of simultaneous optical-UV (KVA, Swift/UVOT), soft X-ray (Swift/XRT) and MAGIC-I VHE data were obtained during different outburst phases. The data constrain the physical parameters of the jet, once the spectral energy distributions obtained are interpreted within the framework of a single-zone SSC leptonic model. Conclusions: The main outcome of the study is that within the homogeneous model high Doppler factors (40 <= delta <= 80) are needed to reproduce the observed SED; but this model cannot explain the observed short time-scale variability, while it can be argued that inhomogeneous models could allow for less extreme Doppler factors, more intense magnetic fields and shorter electron cooling times compatible with hour or sub-hour scale variability.

preprint2012arXiv

Observations of the magnetars 4U 0142+61 and 1E 2259+586 with the MAGIC telescopes

Magnetars are an extreme, highly magnetized class of isolated neutron stars whose large X-ray luminosity is believed to be driven by their high magnetic field. In this work we study for the first time the possible very high energy gamma-ray emission above 100 GeV from magnetars, observing the sources 4U 0142+61 and 1E 2259+586. We observed the two sources with atmospheric Cherenkov telescopes in the very high energy range (E > 100 GeV). 4U 0142+61 was observed with the MAGIC I telescope in 2008 for ~25 h and 1E 2259+586 was observed with the MAGIC stereoscopic system in 2010 for ~14 h. The data were analyzed with the standard MAGIC analysis software. Neither magnetar was detected. Upper limits to the differential and integral flux above 200 GeV were computed using the Rolke algorithm. We obtain integral upper limits to the flux of 1.52*10^-12cm^-2 s^-1 and 2.7*10^-12cm^-2 s^-1 with a confidence level of 95% for 4U 0142+61 and 1E 2259+586, respectively. The resulting differential upper limits are presented together with X-ray data and upper limits in the GeV energy range.

preprint2012arXiv

PG 1553+113: five years of observations with MAGIC

We present the results of five years (2005-2009) of MAGIC observations of the BL Lac object PG 1553+113 at very high energies (VHEs, E > 100 GeV). Power law fits of the individual years are compatible with a steady mean photon index Γ= 4.27 $\pm$ 0.14. In the last three years of data, the flux level above 150 GeV shows a clear variability (probability of constant flux < 0.001%). The flux variations are modest, lying in the range from 4% to 11% of the Crab Nebula flux. Simultaneous optical data also show only modest variability that seems to be correlated with VHE gamma ray variability. We also performed a temporal analysis of (all available) simultaneous Fermi/LAT data of PG 1553+113 above 1 GeV, which reveals hints of variability in the 2008-2009 sample. Finally, we present a combination of the mean spectrum measured at very high energies with archival data available for other wavelengths. The mean spectral energy distribution can be modeled with a one-zone Synchrotron Self Compton (SSC) model, which gives the main physical parameters governing the VHE emission in the blazar jet.

preprint2012arXiv

Phase-resolved energy spectra of the Crab Pulsar in the range of 50-400 GeV measured with the MAGIC Telescopes

We use 73 h of stereoscopic data taken with the MAGIC telescopes to investigate the very high-energy (VHE) gamma-ray emission of the Crab pulsar. Our data show a highly significant pulsed signal in the energy range from 50 to 400 GeV in both the main pulse (P1) and the interpulse (P2) phase regions. We provide the widest spectra to date of the VHE components of both peaks, and these spectra extend to the energy range of satellite-borne observatories. The good resolution and background rejection of the stereoscopic MAGIC system allows us to cross-check the correctness of each spectral point of the pulsar by comparison with the corresponding (strong and well-known) Crab nebula flux. The spectra of both P1 and P2 are compatible with power laws with photon indices of 4.0 \pm 0.8 (P1) and 3.42 \pm 0.26 (P2), respectively, and the ratio P1/P2 between the photon counts of the two pulses is 0.54 \pm 0.12. The VHE emission can be understood as an additional component produced by the inverse Compton scattering of secondary and tertiary e\pm pairs on IR-UV photons.

preprint2012arXiv

Quasi-stellar objects in the ALHAMBRA survey. I. Photometric redshift accuracy based on a 23 optical-NIR filter photometry

We characterize the ability of the ALHAMBRA survey to assign accurate photo-z's to BLAGN and QSOs based on their ALHAMBRA very-low-resolution optical-NIR spectroscopy. A sample of 170 spectroscopically identified BLAGN and QSOs have been used together with a library of templates (including SEDs from AGN, normal, starburst galaxies and stars) in order to fit the 23 photometric data points provided by ALHAMBRA in the optical and NIR (20 medium-band optical filters plus the standard JHKs). We find that the ALHAMBRA photometry is able to provide an accurate photo-z and spectral classification for ~88% of the spectroscopic sources over 2.5 deg^2 in different areas of the survey, all of them brighter than m678=23.5 (equivalent to r(SLOAN)~24.0). The derived photo-z accuracy is better than 1% and comparable to the most recent results in other cosmological fields. The fraction of outliers (~12%) is mainly caused by the larger photometric errors for the faintest sources and the intrinsic variability of the BLAGN/QSO population. A small fraction of outliers may have an incorrectly assigned spectroscopic redshift. The definition of the ALHAMBRA survey in terms of the number of filters, filter properties, area coverage and depth is able to provide photometric redshifts for BLAGN/QSOs with a precision similar to any previous survey that makes use of medium-band optical photometry. In agreement with previous literature results, our analysis also reveals that, in the 0<z<4 redshift interval, very accurate photo-z can be obtained without the use of near-IR broadband photometry at the expense of a slight increase of outliers. The NIR importance is expected to increase at higher redshifts (z>4). These results are relevant for the design of future optical follow-ups of surveys with a large fraction of BLAGN, as it is the case for X-rays or radio surveys.

preprint2012arXiv

The 2010 very high energy gamma-ray flare & 10 years of multi-wavelength observations of M 87

Abridged: The giant radio galaxy M 87 with its proximity, famous jet, and very massive black hole provides a unique opportunity to investigate the origin of very high energy (VHE; E>100 GeV) gamma-ray emission generated in relativistic outflows and the surroundings of super-massive black holes. M 87 has been established as a VHE gamma-ray emitter since 2006. The VHE gamma-ray emission displays strong variability on timescales as short as a day. In this paper, results from a joint VHE monitoring campaign on M 87 by the MAGIC and VERITAS instruments in 2010 are reported. During the campaign, a flare at VHE was detected triggering further observations at VHE (H.E.S.S.), X-rays (Chandra), and radio (43 GHz VLBA). The excellent sampling of the VHE gamma-ray light curve enables one to derive a precise temporal characterization of the flare: the single, isolated flare is well described by a two-sided exponential function with significantly different flux rise and decay times. While the overall variability pattern of the 2010 flare appears somewhat different from that of previous VHE flares in 2005 and 2008, they share very similar timescales (~day), peak fluxes (Phi(>0.35 TeV) ~= (1-3) x 10^-11 ph cm^-2 s^-1), and VHE spectra. 43 GHz VLBA radio observations of the inner jet regions indicate no enhanced flux in 2010 in contrast to observations in 2008, where an increase of the radio flux of the innermost core regions coincided with a VHE flare. On the other hand, Chandra X-ray observations taken ~3 days after the peak of the VHE gamma-ray emission reveal an enhanced flux from the core. The long-term (2001-2010) multi-wavelength light curve of M 87, spanning from radio to VHE and including data from HST, LT, VLA and EVN, is used to further investigate the origin of the VHE gamma-ray emission. No unique, common MWL signature of the three VHE flares has been identified.

preprint2012arXiv

The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: cosmological implications of the large-scale two-point correlation function

We obtain constraints on cosmological parameters from the spherically averaged redshift-space correlation function of the CMASS Data Release 9 (DR9) sample of the Baryonic Oscillation Spectroscopic Survey (BOSS). We combine this information with additional data from recent CMB, SN and BAO measurements. Our results show no significant evidence of deviations from the standard flat-Lambda CDM model, whose basic parameters can be specified by Omega_m = 0.285 +- 0.009, 100 Omega_b = 4.59 +- 0.09, n_s = 0.96 +- 0.009, H_0 = 69.4 +- 0.8 km/s/Mpc and sigma_8 = 0.80 +- 0.02. The CMB+CMASS combination sets tight constraints on the curvature of the Universe, with Omega_k = -0.0043 +- 0.0049, and the tensor-to-scalar amplitude ratio, for which we find r < 0.16 at the 95 per cent confidence level (CL). These data show a clear signature of a deviation from scale-invariance also in the presence of tensor modes, with n_s <1 at the 99.7 per cent CL. We derive constraints on the fraction of massive neutrinos of f_nu < 0.049 (95 per cent CL), implying a limit of sum m_nu < 0.51 eV. We find no signature of a deviation from a cosmological constant from the combination of all datasets, with a constraint of w_DE = -1.033 +- 0.073 when this parameter is assumed time-independent, and no evidence of a departure from this value when it is allowed to evolve as w_DE(a) = w_0 + w_a (1 - a). The achieved accuracy on our cosmological constraints is a clear demonstration of the constraining power of current cosmological observations.

preprint2011arXiv

A search for Very High Energy gamma-ray emission from Scorpius X-1 with the MAGIC telescopes

The acceleration of particles up to GeV or higher energies in microquasars has been the subject of considerable theoretical and observational efforts in the past few years. Sco X-1 is a microquasar from which evidence of highly energetic particles in the jet has been found when it is in the so-called Horizontal Branch (HB), a state when the radio and hard X-ray fluxes are higher and a powerful relativistic jet is present. Here we present the first very high energy gamma-ray observations of Sco X-1 obtained with the MAGIC telescopes. An analysis of the whole dataset does not yield a significant signal, with 95% CL flux upper limits above 300 GeV at the level of 2.4x10^{-12} ph/cm^2/s. Simultaneous RXTE observations were conducted to search for TeV emission during particular X-ray states of the source. A selection of the gamma-ray data obtained during the HB based on the X-ray colors did not yield a signal either, with an upper limit of 3.4x10^{-12} ph/cm^2/s. These upper limits place a constraint on the maximum TeV luminosity to non-thermal X-ray luminosity of L_{VHE}/L_{ntX}<0.02, that can be related to a maximum TeV luminosity to jet power ratio of L_{VHE}/L_{j}<10^{-3}. Our upper limits indicate that the underlying high-energy emission physics in Sco X-1 must be inherently different from that of the hitherto detected gamma-ray binaries.

preprint2011arXiv

Axion-like particle imprint in cosmological very-high-energy sources

Discoveries of very high energy (VHE) photons from distant blazars suggest that, after correction by extragalactic background light (EBL) absorption, there is a flatness or even a turn-up in their spectra at the highest energies that cannot be easily explained by the standard framework. Here, it is shown that a possible solution to this problem is achieved by assuming the existence of axion-like particles (ALPs) with masses ~1 neV. The ALP scenario is tested making use of observations of the highest redshift blazars known in the VHE energy regime, namely 3C 279, 3C 66A, PKS 1222+216 and PG 1553+113. In all cases, better fits to the observed spectra are found when including ALPs rather than considering EBL only. Interestingly, quite similar critical energies for photon/ALP conversions are also derived, independently of the source considered.

preprint2011arXiv

CLARA's view on the escape fraction of Lyman-Alpha photons in high redshift galaxies

Using CLARA (Code for Lyman Alpha Radiation Analysis) we constrain the escape fraction of Lyman-Alpha radiation in galaxies in the redshift range 5<z<7, based on the MareNostrum High-z Universe, a SPH cosmological simulation with more than 2 billion particles. We approximate Lyman-Alpha Emitters (LAEs) as dusty gaseous slabs with Lyman-Alpha radiation sources homogeneously mixed in the gas. Escape fractions for such a configuration and for different gas and dust contents are calculated using our newly developed radiative transfer code CLARA. The results are applied to the MareNostrum High-z Universe numerical galaxies. The model shows a weak redshift evolution and good agreement with estimations of the escape fraction as a function of reddening from observations at z \sim 2.2 and z \sim 3. We extend the slab model by including additional dust in a clumpy component in order to reproduce the UV con- tinuum luminosity function and UV colours at redshifts z>~5. The LAE Luminosity Function (LF) based on the extended clumpy model reproduces broadly the bright end of the LF derived from observations at z \sim 5 and z \sim 6. At z \sim 7 our model over-predicts the LF by roughly a factor of four, presumably because the effects of the neutral intergalactic medium are not taken into account. The remaining tension between the observed and simulated faint end of the LF, both in the UV-continuum and Lyman-Alpha at redshifts z \sim 5 and z \sim 6 points towards an overabundance of simulated LAEs hosted in haloes of masses 1.0x10^10h-1Msol < Mh < 4.0x10^10h-1Msol. Given the difficulties in explaining the observed overabundance by dust absorption, a probable origin of the mismatch are the high star formation rates in the simulated haloes around the quoted mass range. A more efficient supernova feedback should be able to regulate the star formation process in the shallow potential wells of these haloes.

preprint2011arXiv

Dark matter searches with Cherenkov telescopes: nearby dwarf galaxies or local galaxy clusters?

In this paper, we compare dwarf galaxies and galaxy clusters in order to elucidate which object class is the best target for gamma-ray DM searches with imaging atmospheric Cherenkov telescopes (IACTs). We have built a mixed dwarfs+clusters sample containing some of the most promising nearby dwarf galaxies (Draco, Ursa Minor, Wilman 1 and Segue 1) and local galaxy clusters (Perseus, Coma, Ophiuchus, Virgo, Fornax, NGC5813 and NGC5846), and then compute their DM annihilation flux profiles by making use of the latest modeling of their DM density profiles. We also include in our calculations the effect of DM substructure. Willman 1 appears as the best candidate in the sample. However, its mass modeling is still rather uncertain, so probably other candidates with less uncertainties and quite similar fluxes, namely Ursa Minor and Segue 1, might be better options. As for galaxy clusters, Virgo represents the one with the highest flux. However, its large spatial extension can be a serious handicap for IACT observations and posterior data analysis. Yet, other local galaxy cluster candidates with more moderate emission regions, such as Perseus, may represent good alternatives. After comparing dwarfs and clusters, we found that the former exhibit annihilation flux profiles that, at the center, are roughly one order of magnitude higher than those of clusters, although galaxy clusters can yield similar, or even higher, integrated fluxes for the whole object once substructure is taken into account. Even when any of these objects are strictly point-like according to the properties of their annihilation signals, we conclude that dwarf galaxies are best suited for observational strategies based on the search of point-like sources, while galaxy clusters represent best targets for analyses that can deal with rather extended emissions. Finally, we study the detection prospects for IACTs [ABRIDGED]

preprint2011arXiv

Detection of the gamma-ray binary LS I +61 303 in a low flux state at Very High Energy gamma-rays with the MAGIC Telescopes in 2009

We present very high energy (VHE, E > 100 GeV) γ-ray observations of the γ-ray binary system LS I+61 303 obtained with the MAGIC stereo system between 2009 October and 2010 January. We detect a 6.3σ γ-ray signal above 400 GeV in the combined data set. The integral flux above an energy of 300 GeV is F(E>300 GeV)=(1.4 +- 0.3stat +- 0.4syst) * 10^{-12} cm^{-2} s^{-1}, which corresponds to about 1.3% of the Crab Nebula flux in the same energy range. The orbit-averaged flux of LS I +61 303 in the orbital phase interval 0.6--0.7, where a maximum of the TeV flux is expected, is lower by almost an order of magnitude compared to our previous measurements between 2005 September and 2008 January. This provides evidence for a new low emission state in LS I +61 303. We find that the change to the low flux state cannot be solely explained by an increase of photon-photon absorption around the compact star.

preprint2011arXiv

Extragalactic background light inferred from AEGIS galaxy SED-type fractions

The extragalactic background light (EBL) is of fundamental importance both for understanding the entire process of galaxy evolution and for gamma-ray astronomy. However, the overall spectrum of the EBL between 0.1 and 1000 microns has never been determined directly neither from observed luminosity functions (LFs), over a wide redshift range, nor from any multiwavelength observation of galaxy spectral energy distributions (SEDs). The evolving, overall spectrum of the EBL is derived here utilizing a novel method based on observations only. The changing fractions of quiescent galaxies, star-forming galaxies, starburst galaxies and active galactic nucleus (AGN) galaxies from redshift 0.2 to 1 are estimated, and two alternative extrapolations of SED types to higher redshifts are considered. This allows calculation of the evolving EBL. The EBL uncertainties in our modelling based directly on the data are quantified, and their consequences for attenuation of very-high-energy gamma-rays due to pair production on the EBL are discussed. It is concluded that the EBL seems well constrained from the UV to the mid-IR at an intensity level roughly matching galaxy count data. Independent efforts from IR and gamma-ray astronomy are needed in order to reduce the uncertainties in the far-IR.

preprint2011arXiv

Gamma-ray excess from a stacked sample of high- and intermediate-frequency peaked blazars observed with the MAGIC telescope

Between 2004 and 2009 a sample of 28 X-ray selected high- and intermediate-frequency peaked blazars with a X-ray flux larger than 2 uJy at 1 keV in the redshift range from 0.018 to 0.361 was observed with the MAGIC telescope at energies above 100 GeV. Seven among them were detected and the results of these observations are discussed elsewhere. Here we concentrate on the remaining 21 blazars which were not detected during this observation campaign and present the 3 sigma (99.7 %) confidence upper limits on their flux. The individual flux upper limits lie between 1.6 % and 13.6 % of the integral flux from the Crab Nebula. Applying a stacking method to the sample of non-detections with a total of 394.1 hours exposure time, we find evidence for an excess with a cumulative significance of 4.9 standard deviations. It is not dominated by individual objects or flares, but increases linearly with the observation time as for a constant source with an integral flux level of ~1.5 % of that observed from the Crab Nebula above 150 GeV.

preprint2011arXiv

MAGIC contributions to the 32nd International Cosmic Ray Conference

Compilation of the papers contributed by the MAGIC collaboration to the 32nd International Cosmic Ray Conference (ICRC 2011), which took place between August 11 and 18, 2011 in Beijing, China. The papers are sorted in 6 categories: Overview and Highlight papers; Instrument, software and techniques; Galactic sources; Extragalactic sources; Multi-wavelength and joint campaigns; Fundamental physics, dark matter and cosmic rays.

preprint2011arXiv

MAGIC discovery of VHE Emission from the FSRQ PKS 1222+21

Very High Energy (VHE) gamma-ray emission from the flat spectrum radio quasar (FSRQ) PKS 1222+21 (4C 21.35, z=0.432) was detected with the MAGIC Cherenkov telescopes during a short observation (~0.5 hr) performed on 2010 June 17. The MAGIC detection coincides with high energy MeV/GeV gamma-ray activity measured by the Large Area Telescope (LAT) on board the Fermi satellite. The VHE spectrum measured by MAGIC extends from about 70 GeV up to at least 400 GeV and can be well described by a power law dN/dE \propto E^-Gamma with a photon index Gamma= 3.75+/-0.27stat +/-0.2syst. The averaged integral flux above 100 GeV is (4.56+/-0.46)x10^(-10) cm^-2 s^-1 (~1 Crab Nebula flux). The VHE flux measured by MAGIC varies significantly within the 30 min exposure implying a flux doubling time of about 10 min. The VHE and MeV/GeV spectra, corrected for the absorption by the extragalactic background light (EBL), can be described by a single power law with photon index 2.72+/-0.34 between 3 GeV and 400 GeV, and is consistent with emission belonging to a single component in the jet. The absence of a spectral cutoff constrains the gamma-ray emission region outside the Broad Line Region, which would otherwise absorb the VHE gamma-rays. Together with the detected fast variability, this challenges present emission models from jets in FSRQ. Moreover, the combined Fermi/LAT and MAGIC spectral data yield constraints on the density of the Extragalactic Background Light in the UV-optical to near-infrared range that are compatible with recent models.

preprint2011arXiv

MAGIC observations and multiwavelength properties of the quasar 3C279 in 2007 and 2009

Context. 3C 279, the first quasar discovered to emit VHE gamma-rays by the MAGIC telescope in 2006, was reobserved by MAGIC in January 2007 during a major optical flare and from December 2008 to April 2009 following an alert from the Fermi space telescope on an exceptionally high gamma -ray state. Aims. The January 2007 observations resulted in a detection on January 16 with significance 5.2 sigma, corresponding to a F(> 150 GeV) (3.8 \pm 0.8) \cdot 10^-11 ph cm^-2 s^-1 while the overall data sample does not show significant signal. The December 2008 - April 2009 observations did not detect the source. We study the multiwavelength behavior of the source at the epochs of MAGIC observations, collecting quasi-simultaneous data at optical and X-ray frequencies and for 2009 also gamma-ray data from Fermi. Methods. We study the light curves and spectral energy distribution of the source. The spectral energy distributions of three observing epochs (including the February 2006, which has been previously published in Albert et al. 2008a) are modeled with one-zone inverse Compton models and the emission on January 16, 2007 also with two zone model and with a lepto-hadronic model. Results. We find that the VHE gamma-ray emission detected in 2006 and 2007 challenges standard one-zone model, based on relativistic electrons in a jet scattering broad line region photons, while the other studied models fit the observed spectral energy distribution more satisfactorily.

preprint2011arXiv

MAGIC Upper Limits for two Milagro-detected, Bright Fermi Sources in the Region of SNR G65.1+0.6

We report on the observation of the region around supernova remnant G65.1+0.6 with the stand-alone MAGIC-I telescope. This region hosts the two bright GeV gamma-ray sources 1FGL J1954.3+2836 and 1FGL J1958.6+2845. They are identified as GeV pulsars and both have a possible counterpart detected at about 35 TeV by the Milagro observatory. MAGIC collected 25.5 hours of good quality data, and found no significant emission in the range around 1 TeV. We therefore report differential flux upper limits, assuming the emission to be point-like (<0.1 deg) or within a radius of 0.3 deg. In the point-like scenario, the flux limits around 1 TeV are at the level of 3 % and 2 % of the Crab Nebula flux, for the two sources respectively. This implies that the Milagro emission is either extended over a much larger area than our point spread function, or it must be peaked at energies beyond 1 TeV, resulting in a photon index harder than 2.2 in the TeV band.

preprint2011arXiv

Measuring equality horizon with the zero-crossing of the galaxy correlation function

The size of the horizon at the matter-radiation equality is a key scale of the Big Bang cosmology that is directly related to the energy-matter content of the Universe. In this letter, we argue that this scale can be accurately measured from the observed clustering of galaxies in new large scale surveys. We demonstrate that the zero-crossing, r_c, of the 2-point galaxy correlation function is closely related to the horizon size at matter-radiation equality for a large variety of flat LCDM models. Using large-volume cosmological simulations, we also show that the pristine zero-crossing is unaltered by non-linear evolution of density fluctuations, redshift distortions and galaxy biases. This makes r_c a very powerful standard ruler that can be accurately measured, at a percent level, in upcoming experiments that will collect redshifts of millions of galaxies and quasars.

preprint2011arXiv

Modelling the fraction of Lyman Break Galaxies with strong Lyman alpha emission at 5 < z < 7

We present theoretical results for the expected fraction of Lyman Break Galaxies (LBGs) to be detected as strong Lyman alpha emitters (LAEs) in the redshift range 5 < z < 7. We base our analysis on the 2-billion particle SPH simulation MareNostrum High-z Universe. We approximate galaxies as static dusty slabs with an additional clumpy dust distribution affecting stellar populations younger than 25 Myr. The model for the Lyman alpha escape fraction is based on the results of our Monte-Carlo radiative transfer code (CLARA) for a slab configuration. We also fix the transmission of Lyman alpha photons through the intergalactic medium to a constant value of 50% at all redshifts. From the results of this model we calculate xLya, the fraction of Lyman Break Galaxies with Lyman alpha equivalent width (EW) larger than 50 Angstrom. We find a remarkable agreement with observational data at 4.5 < z < 6. For bright (-22 < MUV < -20.5) and faint (-20.5 < MUV < -18.5) galaxies our model predicts xLya = 0.02 \pm 0.01 and xLya = 0.47 \pm 0.01 while observers report xLya = 0.08 \pm 0.02 and xLya = 0.47 \pm 0.16, respectively. Additional evolution of the extinction model at redshift z \sim 7, that decreases the intensity of transmitted Lyman alpha radiation by a factor of f_T = 0.4 as to match the LAE luminosity function at z \sim 6.5, naturally provides a good match for the recently reported xLya fractions at z > 6.3. Exploring different toy models for the Lyman alpha escape fraction, we show that a decreasing Lyman alpha escape fraction with increasing UV galaxy luminosity is a key element in our model to explain the of larger xLya fractions for fainter LBGs.

preprint2011arXiv

Observations of the Crab pulsar between 25 GeV and 100 GeV with the MAGIC I telescope

We report on the observation of $γ$-rays above 25\,GeV from the Crab pulsar (PSR B0532+21) using the MAGIC I telescope. Two data sets from observations during the winter period 2007/2008 and 2008/2009 are used. In order to discuss the spectral shape from 100\,MeV to 100\,GeV, one year of public {\it Fermi} Large Area Telescope ({\it Fermi}-LAT) data are also analyzed to complement the MAGIC data. The extrapolation of the exponential cutoff spectrum determined with the Fermi-LAT data is inconsistent with MAGIC measurements, which requires a modification of the standard pulsar emission models. In the energy region between 25 and 100\,GeV, the emission in the P1 phase (from -0.06 to 0.04, location of the main pulse) and the P2 phase (from 0.32 to 0.43, location of the interpulse) can be described by power laws with spectral indices of $-3.1 \pm 1.0_{stat} \pm 0.3_{syst}$ and $-3.5 \pm 0.5_{stat} \pm 0.3_{syst}$, respectively. Assuming an asymmetric Lorentzian for the pulse shape, the peak positions of the main pulse and the interpulse are estimated to be at phases $-0.009 \pm 0.007$ and $0.393 \pm 0.009$, while the full widths at half maximum are $0.025 \pm 0.008$ and $0.053 \pm 0.015$, respectively.

preprint2011arXiv

Performance of the MAGIC stereo system obtained with Crab Nebula data

MAGIC is a system of two Imaging Atmospheric Cherenkov Telescopes located in the Canary island of La Palma. Since autumn 2009 both telescopes have been working together in stereoscopic mode, providing a significant improvement with respect to the previous single-telescope observations. We use observations of the Crab Nebula taken at low zenith angles to assess the performance of the MAGIC stereo system. The trigger threshold of the MAGIC telescopes is 50-60 GeV. Advanced stereo analysis techniques allow MAGIC to achieve a sensitivity as good as (0.76 +/- 0.03)% of the Crab Nebula flux in 50 h of observations above 290 GeV. The angular resolution at those energies is better than ~0.07 degree. We also perform a detailed study of possible systematic effects which may influence the analysis of the data taken with the MAGIC telescopes.

preprint2011arXiv

Searches for Dark Matter annihilation signatures in the Segue 1 satellite galaxy with the MAGIC-I telescope

We report the results of the observation of the nearby satellite galaxy Segue 1 performed by the MAGIC-I ground-based gamma-ray telescope between November 2008 and March 2009 for a total of 43.2 hours. No significant gamma-ray emission was found above the background. Differential upper limits on the gamma-ray flux are derived assuming various power-law slopes for the possible emission spectrum. Integral upper limits are also calculated for several power-law spectra and for different energy thresholds. The values are of the order of 10^{-11} ph cm^{-2}$ s^{-1} above 100 GeV and 10^{-12} ph cm^{-2} s^{-1} above 200 GeV. Segue 1 is currently considered one of the most interesting targets for indirect dark matter searches. In these terms, the upper limits have been also interpreted in the context of annihilating dark matter particles. For such purpose, we performed a grid scan over a reasonable portion of the parameter space for the minimal SuperGravity model and computed the flux upper limit for each point separately, taking fully into account the peculiar spectral features of each model. We found that in order to match the experimental upper limits with the model predictions, a minimum flux boost of 10^{3} is required, and that the upper limits are quite dependent on the shape of the gamma-ray energy spectrum predicted by each specific model. Finally we compared the upper limits with the predictions of some dark matter models able to explain the PAMELA rise in the positron ratio, finding that Segue 1 data are in tension with the dark matter explanation of the PAMELA spectrum in the case of a dark matter candidate annihilating into tau+tau-. A complete exclusion however is not possible due to the uncertainties in the Segue 1 astrophysical factor.

preprint2011arXiv

The BigBOSS Experiment

BigBOSS is a Stage IV ground-based dark energy experiment to study baryon acoustic oscillations (BAO) and the growth of structure with a wide-area galaxy and quasar redshift survey over 14,000 square degrees. It has been conditionally accepted by NOAO in response to a call for major new instrumentation and a high-impact science program for the 4-m Mayall telescope at Kitt Peak. The BigBOSS instrument is a robotically-actuated, fiber-fed spectrograph capable of taking 5000 simultaneous spectra over a wavelength range from 340 nm to 1060 nm, with a resolution R = 3000-4800. Using data from imaging surveys that are already underway, spectroscopic targets are selected that trace the underlying dark matter distribution. In particular, targets include luminous red galaxies (LRGs) up to z = 1.0, extending the BOSS LRG survey in both redshift and survey area. To probe the universe out to even higher redshift, BigBOSS will target bright [OII] emission line galaxies (ELGs) up to z = 1.7. In total, 20 million galaxy redshifts are obtained to measure the BAO feature, trace the matter power spectrum at smaller scales, and detect redshift space distortions. BigBOSS will provide additional constraints on early dark energy and on the curvature of the universe by measuring the Ly-alpha forest in the spectra of over 600,000 2.2 < z < 3.5 quasars. BigBOSS galaxy BAO measurements combined with an analysis of the broadband power, including the Ly-alpha forest in BigBOSS quasar spectra, achieves a FOM of 395 with Planck plus Stage III priors. This FOM is based on conservative assumptions for the analysis of broad band power (kmax = 0.15), and could grow to over 600 if current work allows us to push the analysis to higher wave numbers (kmax = 0.3). BigBOSS will also place constraints on theories of modified gravity and inflation, and will measure the sum of neutrino masses to 0.024 eV accuracy.

preprint2010arXiv

Detection of very high energy gamma-ray emission from the Perseus cluster head-tail galaxy IC 310 by the MAGIC telescopes

We report on the detection with the MAGIC telescopes of very high energy gamma-rays from IC 310, a head-tail radio galaxy in the Perseus galaxy cluster, observed during the interval November 2008 to February 2010. The Fermi satellite has also detected this galaxy. The source is detected by MAGIC at a high statistical significance of 7.6sigma in 20.6 hr of stereo data. The observed spectral energy distribution is flat with a differential spectral index of -2.00 \pm 0.14. The mean flux above 300 GeV, between October 2009 and February 2010, (3.1 \pm 0.5)x10^{-12} cm^{-2} s^{-1}, corresponds to (2.5 \pm 0.4)% of Crab Nebula units. Only an upper limit, of 1.9% of Crab Nebula units above 300 GeV, was obtained with the 2008 data. This, together with strong hints (>3sigma) of flares in the middle of October and November 2009, implies that the emission is variable. The MAGIC results favour a scenario with the very high energy emission originating from the inner jet close to the central engine. More complicated models than a simple one-zone SSC scenario, e.g. multi-zone SSC, external Compton or hadronic, may be required to explain the very flat spectrum and its extension over more than three orders of magnitude in energy.

preprint2010arXiv

Extreme Multiplex Spectrograph: An efficient mechanical design for high-demanding requirements

XMS is a multi-channel wide-field spectrograph designed for the prime focus of the 3.5m Calar-Alto telescope. The instrument is composed by four quadrants, each of which contains a spectrograph channel. An innovative mechanical design -at concept/preliminary stage- has been implemented to: 1) Minimize the separation between the channels to achieve maximal filling factor; 2) Cope with the very constraining space and mass overall requirements; 3) Achieve very tight alignment tolerances; 4) Provide lens self-centering under large temperature excursions; 5) Provide masks including 4000 slits (edges thinner than 100μ). An overview of this extremely challenging mechanical design is here presented.

preprint2010arXiv

Herschel FIR counterparts of selected Ly-alpha emitters at z~2.2. Fast evolution since z~3 or missed obscured AGNs?

Ly-alpha emitters (LAEs) are seen everywhere in the redshift domain from local to z~7. Far-infrared (FIR) counterparts of LAEs at different epochs could provide direct clues on dust content, extinction, and spectral energy distribution (SED) for these galaxies. We search for FIR counterparts of LAEs that are optically detected in the GOODS-North field at redshift z~2.2 using data from the Herschel Space Telescope with the Photodetector Array Camera and Spectrometer (PACS). The LAE candidates were isolated via color-magnitude diagram using the medium-band photometry from the ALHAMBRA Survey, ancillary data on GOODS-North, and stellar population models. According to the fitting of these spectral synthesis models and FIR/optical diagnostics, most of them seem to be obscured galaxies whose spectra are AGN-dominated. From the analysis of the optical data, we have observed a fraction of AGN or composite over source total number of ~0.75 in the LAE population at z~2.2, which is marginally consistent with the fraction previously observed at z=2.25 and even at low redshift (0.2<z<0.45), but significantly different from the one observed at redshift ~3, which could be compatible either with a scenario of rapid change in the AGN fraction between the epochs involved or with a non detection of obscured AGN in other z=2-3 LAE samples due to lack of deep FIR observations. We found three robust FIR (PACS) counterparts at z~2.2 in GOODS-North. This demonstrates the possibility of finding dust emission in LAEs even at higher redshifts.

preprint2010arXiv

Impact of internal bremsstrahlung on the detection of gamma-rays from neutralinos

We present a detailed study of the effect of internal bremsstrahlung photons in the context of the minimal supersymmetric standard models and their impact on gamma-ray dark matter annihilation searches. We find that although this effect has to be included for the correct evaluation of fluxes of high energy photons from neutralino annihilation, its contribution is relevant only in models and at energies where the lines contribution is dominant over the secondary photons. Therefore, we find that the most optimistic supersymmetric scenarios for dark matter detection do not change significantly when including the internal bremsstrahlung. As an example, we review the gamma-ray dark matter detection prospects of the Draco dwarf spheroidal galaxy for the MAGIC stereoscopic system and the CTA project. Though the flux of high energy photons is enhanced by an order of magnitude in some regions of the parameter space, the expected fluxes are still much below the sensitivity of the instruments.

preprint2010arXiv

Magic constraints on Gamma-ray emission from Cygnus X-3

Cygnus X-3 is a microquasar consisting of an accreting compact object orbiting around a Wolf-Rayet star. It has been detected at radio frequencies and up to high-energy gamma rays (above 100 MeV). However, many models also predict a very high energy (VHE) emission (above hundreds of GeV) when the source displays relativistic persistent jets or transient ejections. Therefore, detecting such emission would improve the understanding of the jet physics. The imaging atmospheric Cherenkov telescope MAGIC observed Cygnus X-3 for about 70 hours between 2006 March and 2009 August in different X-ray/radio spectral states and also during a period of enhanced gamma-ray emission. MAGIC found no evidence for a VHE signal from the direction of the microquasar. An upper limit to the integral flux for energies higher than 250 GeV has been set to 2.2 x 10-12 photons cm-2 s-1 (95% confidence level). This is the best limit so far to the VHE emission from this source. The non-detection of a VHE signal during the period of activity in the high-energy band sheds light on the location of the possible VHE radiation favoring the emission from the innermost region of the jets, where absorption is significant. The current and future generations of Cherenkov telescopes may detect a signal under precise spectral conditions.

preprint2010arXiv

MAGIC Gamma-ray Telescope Observation of the Perseus Cluster of Galaxies: implications for cosmic rays, dark matter, and NGC1275

The Perseus galaxy cluster was observed by the MAGIC Cherenkov telescope for a total effective time of 24.4 hr during 2008 November and December. The resulting upper limits on the gamma-ray emission above 100 GeV are in the range of 4.6 to 7.5 x 10^{-12} cm^{-2} s^{-1} for spectral indices from -1.5 to -2.5, thereby constraining the emission produced by cosmic rays, dark matter annihilations, and the central radio galaxy NGC1275. Results are compatible with cosmological cluster simulations for the cosmic-ray-induced gamma-ray emission, constraining the average cosmic ray-to-thermal pressure to <4% for the cluster core region (<8% for the entire cluster). Using simplified assumptions adopted in earlier work (a power-law spectrum with an index of -2.1, constant cosmic ray-to-thermal pressure for the peripheral cluster regions while accounting for the adiabatic contraction during the cooling flow formation), we would limit the ratio of cosmic ray-to-thermal energy to E_CR/E_th<3%. The upper limit also translates into a level of gamma-ray emission from possible annihilations of the cluster dark matter (the dominant mass component) that is consistent with boost factors of ~10^4 for the typically expected dark matter annihilation-induced emission. Finally, the upper limits obtained for the gamma-ray emission of the central radio galaxy NGC1275 are consistent with the recent detection by the Fermi-LAT satellite. Due to the extremely large Doppler factors required for the jet, a one-zone synchrotron self-Compton model is implausible in this case. We reproduce the observed spectral energy density by using the structured jet (spine-layer) model which has previously been adopted to explain the high-energy emission of radio galaxies.

preprint2010arXiv

MAGIC TeV Gamma-Ray Observations of Markarian 421 during Multiwavelength Campaigns in 2006

The Major Atmospheric Gamma Imaging Cerenkov (MAGIC) telescope participated in three multiwavelength (MWL) campaigns, observing the blazar Markarian (Mkn) 421 during the nights of 2006 April 28, 29, and 2006 June 14. We analyzed the corresponding MAGIC very-high energy observations during 9 nights from 2006 April 22 to 30 and on 2006 June 14. We inferred light curves with sub-day resolution and night-by-night energy spectra. A strong gamma-ray signal was detected from Mkn 421 on all observation nights. The flux (E > 250 GeV) varied on night-by-night basis between (0.92+-0.11)10^-10 cm^-2 s^-1 (0.57 Crab units) and (3.21+-0.15)10^-10 cm^-2 s^-1 (2.0 Crab units) in 2006 April. There is a clear indication for intra-night variability with a doubling time of 36+-10(stat) minutes on the night of 2006 April 29, establishing once more rapid flux variability for this object. For all individual nights gamma-ray spectra could be inferred, with power-law indices ranging from 1.66 to 2.47. We did not find statistically significant correlations between the spectral index and the flux state for individual nights. During the June 2006 campaign, a flux substantially lower than the one measured by the Whipple 10-m telescope four days later was found. Using a log-parabolic power law fit we deduced for some data sets the location of the spectral peak in the very-high energy regime. Our results confirm the indications of rising peak energy with increasing flux, as expected in leptonic acceleration models.

preprint2010arXiv

Observations of the Blazar 3C 66A with the MAGIC Telescopes in Stereoscopic Mode

We report new observations of the intermediate-frequency peaked BL Lacertae object 3C 66A with the MAGIC telescopes. The data sample we use were taken in 2009 December and 2010 January, and comprises 2.3 hr of good quality data in stereoscopic mode. In this period, we find a significant signal from the direction of the blazar 3C 66A. The new MAGIC stereoscopic system is shown to play an essential role for the separation between 3C 66A and the nearby radio galaxy 3C 66B, which is at a distance of only $6^\prime$. The derived integral flux above $100\eh{GeV}$ is 8.3% of Crab Nebula flux and the energy spectrum is reproduced by a power law of photon index $3.64 \pm 0.39_{\rm stat} \pm 0.25_{\rm sys}$. Within errors, this is compatible with the one derived by VERITAS in 2009. From the spectra corrected for absorption by the extragalactic background light, we only find small differences between the four models that we applied, and constrain the redshift of the blazar to $z < 0.68$.

preprint2010arXiv

Search for an extended VHE gamma-ray emission from Mrk 421 and Mrk 501 with the MAGIC Telescope

Context: Part of the very high energy $γ$-ray radiation coming from extragalactic sources is absorbed through the pair production process on the extragalactic background light photons. Extragalactic magnetic fields alter the trajectories of these cascade pairs and, in turn, convert cosmic background photons to gamma-ray energies by inverse Compton scattering. These secondary photons can form an extended halo around bright VHE sources. Aims: We searched for an extended emission around the bright blazars Mrk 421 and Mrk 501 using the MAGIC telescope data. Methods: If extended emission is present, the angular distribution of reconstructed gamma-ray arrival directions around the source is broader than for a point-like source. In the analysis of a few tens of hours of observational data taken from Mrk 421 and Mrk 501 we used a newly developed method that provides better angular resolution. This method is based on the usage of multidimensional decision trees. Comparing the measured shapes of angular distributions with those expected from a point-like source one can detect or constrain possible extended emission around the source. We also studied the influence of different types of systematic errors on the shape of the distribution of reconstructed gamma-ray arrival directions for a point source. Results: We present upper limits for an extended emission calculated for both sources for various source extensions and emission profiles. We discuss possible constraints on the extragalactic magnetic fields strength. We obtain upper limits on the extended emission around the Mrk~421 (Mrk~501) on the level of < 5% (< 4%) of the Crab Nebula flux above the energy threshold of 300 GeV. Using these results we discuss possible constraints on the extragalactic magnetic fields strength around a few times $10^{-15}$~G.

preprint2010arXiv

Simultaneous multi-frequency observation of the unknown redshift blazar PG 1553+113 in March-April 2008

The blazar PG 1553+113 is a well known TeV gamma-ray emitter. In this paper, we determine its spectral energy distribution using simultaneous multi-frequency data in order to study its emission processes. An extensive campaign was carried out between March and April 2008, where optical, X-ray, high-energy (HE) gamma-ray, and very-high-energy (VHE) gamma-ray data were obtained with the KVA, Abastumani, REM, RossiXTE/ASM, AGILE and MAGIC telescopes, respectively. This is the first simultaneous broad-band (i.e., HE+VHE) gamma-ray observation, though AGILE did not detect the source. We combine data to derive source's spectral energy distribution and interpret its double peaked shape within the framework of a synchrotron self compton model

preprint2010arXiv

Spectral Energy Distribution of Markarian 501: Quiescent State vs. Extreme Outburst

The very high energy (VHE; E > 100 GeV) blazar Markarian 501 has a well-studied history of extreme spectral variability and is an excellent laboratory for studying the physical processes within the jets of active galactic nuclei. However, there are few detailed multiwavelength studies of Markarian 501 during its quiescent state, due to its low luminosity. A short-term multiwavelength study of Markarian 501 was coordinated in March 2009, focusing around a multi-day observation with the Suzaku X-ray satellite and including γ-ray data from VERITAS, MAGIC, and the Fermi Gamma-ray Space Telescope with the goal of providing a well-sampled multiwavelength baseline measurement of Markarian 501 in the quiescent state. The results of these quiescent-state observations are compared to the historically extreme outburst of April 16, 1997, with the goal of examining variability of the spectral energy distribution between the two states. The derived broadband spectral energy distribution shows the characteristic double-peaked profile. We find that the X-ray peak shifts by over two orders of magnitude in photon energy between the two flux states while the VHE peak varies little. The limited shift in the VHE peak can be explained by the transition to the Klein-Nishina regime. Synchrotron self-Compton models are matched to the data and the implied Klein-Nishina effects are explored.

preprint2010arXiv

The ALHAMBRA photometric system

This paper presents the characterization of the optical range of the ALHAMBRA photometric system, a 20 contiguous, equal-width, medium-band CCD system with wavelength coverage from 3500A to 9700A. The photometric description of the system is done by presenting the full response curve as a product of the filters, CCD and atmospheric transmission curves, and using some first and second order moments of this response function. We also introduce the set of standard stars that defines the system, formed by 31 classic spectrophotometric standard stars which have been used in the calibration of other known photometric systems, and 288 stars, flux calibrated homogeneously, from the Next Generation Spectral Library (NGSL). Based on the NGSL, we determine the transformation equations between Sloan Digital Sky Survey (SDSS) ugriz photometry and the ALHAMBRA photometric system, in order to establish some relations between both systems. Finally we develop and discuss a strategy to calculate the photometric zero points of the different pointings in the ALHAMBRA project.

preprint2010arXiv

The clustering of massive galaxies at z~0.5 from the first semester of BOSS data

We calculate the real- and redshift-space clustering of massive galaxies at z~0.5 using the first semester of data by the Baryon Oscillation Spectroscopic Survey (BOSS). We study the correlation functions of a sample of 44,000 massive galaxies in the redshift range 0.4<z<0.7. We present a halo-occupation distribution modeling of the clustering results and discuss the implications for the manner in which massive galaxies at z~0.5 occupy dark matter halos. The majority of our galaxies are central galaxies living in halos of mass 10^{13}Msun/h, but 10% are satellites living in halos 10 times more massive. These results are broadly in agreement with earlier investigations of massive galaxies at z~0.5. The inferred large-scale bias (b~2) and relatively high number density (nbar=3e-4 h^3 Mpc^{-3}) imply that BOSS galaxies are excellent tracers of large-scale structure, suggesting BOSS will enable a wide range of investigations on the distance scale, the growth of large-scale structure, massive galaxy evolution and other topics.

preprint2009arXiv

Near-IR Galaxy Counts and Evolution from the Wide-Field ALHAMBRA survey

The ALHAMBRA survey aims to cover 4 square degrees using a system of 20 contiguous, equal width, medium-band filters spanning the range 3500 A to 9700 A plus the standard JHKs filters. Here we analyze deep near-IR number counts of one of our fields (ALH08) for which we have a relatively large area (0.5 square degrees) and faint photometry (J=22.4, H=21.3 and K=20.0 at the 50% of recovery efficiency for point-like sources). We find that the logarithmic gradient of the galaxy counts undergoes a distinct change to a flatter slope in each band: from 0.44 at [17.0, 18.5] to 0.34 at [19.5, 22.0] for the J band; for the H band 0.46 at [15.5, 18.0] to 0.36 at [19.0, 21.0], and in Ks the change is from 0.53 in the range [15.0, 17.0] to 0.33 in the interval [18.0, 20.0]. These observations together with faint optical counts are used to constrain models that include density and luminosity evolution of the local type-dependent luminosity functions. Our models imply a decline in the space density of evolved early-type galaxies with increasing redshift, such that only 30% - 50% of the bulk of the present day red-ellipticals was already in place at z~1.

preprint2009arXiv

Simulated vs. observed UV emission at high redshift: a hint for a clumpy ISM?

We discuss the rest-frame UV emission between 5< z < 7 from the MareNostrum High-z Universe, a SPH simulation done with more than 2 billion particles. Cosmological simulations of galaxy formation generally overpredict the UV restframe luminosity function at high redshift, both at the bright and faint ends. In this Letter we explore a dust attenuation model where a larger extinction is applied to star populations younger than a given age, mimicking the effect of a clumpy interstellar medium. We show that this scenario fits reasonably well both the UV luminosity functions and the UV-continuum slopes derived from observations. The model assumes a large obscuration for stars younger than 25 Myr from the gas clouds where they should be embedded at their formation time. We find that the optical depth in these clouds should be between 30 and 100 times larger than the mean optical depth for the homogeneous part of the interstellar medium. These values are one order of magnitude larger than those estimated in local galaxies. Therefore, we conclude that LambdaCDM predictions for the high-z UV emission can accommodate the current observations if we consider a dust extinction model based on the assumption of a clumpy environment at high redshift.

preprint2005arXiv

GAW, Gamma Air Watch - A Large Field of View Imaging Atmospheric Cherenkov

GAW, acronym for Gamma Air Watch, is a path-finder experiment to test the feasibility of a new generation of imaging atmospheric Cherenkov telescopes that join high flux sensitivity with large field of view capability. GAW is conceived as an array of three identical imaging telescopes disposed at the vertexes of an equilateral triangle, about 80 m side. Two main features characterize GAW with respect to all the existing and presently planned ground-based Cherenkov telescopes. The first difference concerns the optics system: GAW uses a Fresnel refractive lens (2.13 m diameter) as light collector instead of classical reflective mirror. The second main difference is the detection working mode used: the detector at the focal surface operates in single photoelectron counting mode instead of the usual charge integration one. The GAW array is planned to be located at the Calar Alto Observatory site, Spain, 2150 m a.s.l. During its first phase, only 6x6 degrees of the focal plane detector will be implemented; moving it along the field of view, the sensitivity of the telescopes will be tested observing the Crab Nebula with on-axis and off-axis pointing up to 20 deg and with energy threshold of 300 GeV and energy peak of 700 GeV. As path-finder, GAW will also monitor the Very High Energy activity of some flaring Blazars as well as will follow-up GLAST detections at high energies. In a second phase, the focal plane will be enlarged to cover a field of view of 24 deg (full); pointing along different North-South directions, GAW would reach a survey of 360x60 degrees region of the sky. GAW is a collaboration effort of Research Institutes in Italy, Portugal and Spain.

preprint2004arXiv

Do the Unidentified EGRET Sources Trace Annihilating Dark Matter in the Local Group?

In a cold dark matter (CDM) framework of structure formation, the dark matter haloes around galaxies assemble through successive mergers with smaller haloes. This merging process is not completely efficient, and hundreds of surviving halo cores, or {\it subhaloes}, are expected to remain in orbit within the halo of a galaxy like the Milky Way. While the dozen visible satellites of the Milky Way may trace some of these subhaloes, the majority are currently undetected. A large number of high-velocity clouds (HVCs) of neutral hydrogen {\it are} observed around the Milky Way, and it is plausible that some of the HVCs may trace subhaloes undetected in the optical. Confirming the existence of concentrations of dark matter associated with even a few of the HVCs would represent a dramatic step forward in our attempts to understand the nature of dark matter. Supersymmetric (SUSY) extensions of the Standard Model of particle physics currently suggest neutralinos as a natural well-motivated candidate for the non-baryonic dark matter of the universe. If this is indeed the case, then it may be possible to detect dark matter indirectly as it annihilates into neutrinos, photons or positrons. In particular, the centres of subhaloes might show up as point sources in gamma-ray observations. In this work we consider the possibility that some of the unidentified EGRET $γ$-ray sources trace annihilating neutralino dark matter in the dark substructure of the Local Group. We compare the observed positions and fluxes of both the unidentified EGRET sources and the HVCs with the positions and fluxes predicted by a model of halo substructure, to determine to what extent any of these three populations could be associated.

preprint2004arXiv

Indirect searches for SUSY Dark Matter with the MAGIC Cherenkov Telescope

Neutralinos are the natural well-motivated candidates to provide the non-baryonic dark matter of the universe which may produce detectable signals through their annihilation into neutrinos, photons or positrons. Due to its high flux sensitivity and low energy threshold, the MAGIC Cherenkov telescope could potentially detect the neutralino annihilation high energy photon products. In the framework of minimal supergravity, the neutralino SUSY parameter space can be scanned in different benchmark scenarios defined after accelerator and cosmology constraints. Moreover, the neutralino density profiles in galaxy halos and sub-halos have to be understood to infer which is the optimal observation region to be explored by the MAGIC telescope for the detection of neutralino photon signatures within our Galaxy and our Local Group.

preprint2002arXiv

Near-IR spectroscopy of planetary nebulae precursors

We present near-IR spectroscopy of a sample of 30 IRAS sources recently identified as late AGB stars, post-AGB stars or early PNe. The spectra obtained are centered at various wavelengths covering the molecular hydrogen v=1-0 S(1) and v=2-1 S(1) emission lines, the recombination lines of hydrogen Br-gamma, Pf-gamma and Br-alpha, and the CO[v=2-0] first overtone bandhead at 2.294 microns. As a result of these observations we have increased from 4 to 13 the total number of proto-PNe detected in molecular hydrogen. When the molecular hydrogen is fluorescence-excited the detection rate is found to be directly correlated with the evolutionary stage of the central star, rather than with the nebular morphology. In contrast, shocked-excited molecular hydrogen is detected only in strongly bipolar proto-PNe, sometimes even at an early stage in the post-AGB phase. The strong correlation of shocked-excited molecular hydrogen emission with bipolarity found confirms the result previously reported by Kastner et al. (1996) in evolved PNe. However, our results show that this correlation does not exist in the case of fluorescence-excited molecular hydrogen.

preprint1999arXiv

A counterrotating central component in NGC 5728

We present a detailed study of the stellar kinematics in the barred galaxy NGC 5728 based on I-band photometry and long-slit spectroscopic observations in the region of the near-IR Ca II triplet. The analysis of the stellar line-of-sight velocity distribution (LOSVD) has revealed, in the central regions of the bar, the presence of a cold (v/sigma~2.5) prograde $S$-shaped velocity component that coexists in the central 4 kpc with a fainter and hotter (v/sigma~0.5) counterrotating component. Beyond 4 kpc from the nucleus the LOSVD shows the stellar bar kinematics. The comparison of the radial surface brightness profile of the velocity components with that obtained from an I-band image shows that the counterrotating core follows a r^(1/4) profile, while the S-shaped component does not follows the flat bar surface brightness profile. Several possible scenarios accounting for such kinematic signatures found in the center of the bar in NGC 5728 are discussed.

preprint1996arXiv

A Counter-rotating Bulge in the Sb Galaxy NGC 7331

We have found that the bulge of the large, nearby Sb galaxy NGC 7331 rotates retrograde to its disk. Analysis of spectra in the region of the near-IR Ca II triplet along the major axis shows that, in the radial range between 5'' and ~20'', the line of sight velocity distribution of the absorption lines has two distinct peaks, and can be decomposed into a fast-rotating component with v/sigma > 3, and a slower rotating, retrograde component with v/sigma between 1 -- 1.5. The radial surface brightness profile of the counter-rotating component follows that of the bulge, obtained from a 2-dimensional bulge-disk decomposition of a near-infrared K-band image, while the fast rotating component follows the disk. At the radius where the disk starts to dominate the isophotes change from being considerably boxy to very disky. Although a number of spiral galaxies have been found that contain cold, couter-rotating disks, this is the first galaxy known to have a boxy, probably triaxial, fairly warm, counter-rotating component, which is dominating in the central regions. If it is a bar seen end-on, this bar has to be thicker than the disk. We find that NGC 7331, even though it is a fairly early-type spiral, does not have a conventional, co-rotating bulge. The fact that the inner component is retrograde makes us believe that it was formed from infalling material, in either stellar or gaseous form (e.g. Balcells & Quinn 1990). Another possibility however is that the structure has been there since the formation of the galaxy. In this case it will be a challenge to explain the large change in orientation of the angular momentum when going outward radially.