Source author record

A. Herrero

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

ResearcherUnclaimed source record

astro-ph.SR astro-ph.HE astro-ph.GA astro-ph.CO astro-ph.IM astro-ph

Catalog footprint

What is connected

129works

6topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

preprint2025arXiv

Populations of evolved massive binary stars in the Small Magellanic Cloud II: Predictions from rapid binary evolution

Massive star evolution plays a crucial role in astrophysics but bares large uncertainties. This problem becomes more severe by the majority of massive stars being born in close binary systems, whose evolution is affected by the interaction of their components. We want to constrain major uncertainties in massive binary star evolution, in particular the efficiency and the stability of the first mass transfer phase. We use the rapid population synthesis code ComBinE to generate synthetic populations of post-interaction binaries, assuming constant mass-transfer efficiency. We employ a new merger criterion that adjusts self-consistently to any prescribed mass-transfer efficiency. We tailor our synthetic populations to be comparable to the expected binary populations in the Small Magellanic Cloud (SMC). We find that the observed populations of evolved massive binaries can not be reproduced with a single mass-transfer efficiency. Instead, a rather high efficiency (>50%) is needed to reproduce the number of Be stars and Be/X-ray binaries in the SMC, while a low efficiency (~10%) leads to a better agreement with the observed number of Wolf-Rayet stars. We construct a corresponding mass-dependent mass-transfer efficiency recipe to produce our fiducial synthetic SMC post-interaction binary population. It reproduces the observed number and properties of the Be/X-ray and WR-binaries rather well, and is not in stark disagreement with the observed OBe star population. It further predicts two large, yet unobserved populations of OB+BH binaries, that is ~100 OB+BH systems with rather small orbital periods (<20 days) and ~40 longer period OBe+BH systems.

preprint2022arXiv

The Gaia-ESO Survey: The analysis of the hot-star spectra

The Gaia-ESO Survey (GES) is a large public spectroscopic survey that has collected, over a period of 6 years, spectra of ~ 10^5 stars. This survey provides not only the reduced spectra, but also the stellar parameters and abundances resulting from the analysis of the spectra. The GES dataflow is organised in 19 working groups. Working group 13 (WG13) is responsible for the spectral analysis of the hottest stars (O, B and A type, with a formal cut-off of Teff > 7000 K) that were observed as part of GES. We present the procedures and techniques that have been applied to the reduced spectra, in order to determine the stellar parameters and abundances of these stars. The procedure used is similar to that of other working groups in GES. A number of groups (called `Nodes') each independently analyse the spectra, using their state-of-the-art techniques and codes. Specific for the analysis in WG13 is the large temperature range that is covered (Teff = 7000 - 50,000 K), requiring the use of different analysis codes. Most Nodes can therefore only handle part of the data. Quality checks are applied to the results of these Nodes by comparing them to benchmark stars, and by comparing them one to another. For each star the Node values are then homogenised into a single result: the recommended parameters and abundances. Eight Nodes each analysed (part of) the data. In total 17,693 spectra of 6462 stars were analysed, most of them in 37 open star clusters. The homogenisation led to stellar parameters for 5584 stars. Abundances were determined for a more limited number of stars. Elements studied are He, C, N, O, Ne, Mg, Al, Si and Sc. Abundances for at least one of those elements were determined for 292 stars. The hot-star data analysed here, as well as the Gaia-ESO Survey data in general, will be of considerable use in future studies of stellar evolution and open clusters.

preprint2022arXiv

The nature of the Cygnus extreme B-supergiant 2MASS J20395358+4222505

2MASS J20395358+4222505 is an obscured early B supergiant near the massive OB star association Cyg OB2. Despite its bright infrared magnitude (K$_{s}$=5.82) it has remained largely ignored because of its dim optical magnitude (B=16.63, V=13.68). In a previous paper we classified it as a highly reddened, potentially extremely luminous, early B-type supergiant. We obtained its spectrum in the U, B and R spectral bands during commissioning observations with the instrument MEGARA@GTC. It displays a particularly strong H$α$ emission for its spectral type, B1 Ia. The star seems to be in an intermediate phase between super- and hypergiant, a group that it will probably join in the near (astronomical) future. We observe a radial velocity difference between individual observations and determine the stellar parameters, obtaining T$_{eff}$ = 24000 K, logg$_{c}$= 2.88 $\pm$ 0.15. The rotational velocity found is large for a B-supergiant, vsini= 110 $\pm$ 25 km s$^{-1}$. The abundance pattern is consistent with solar, with a mild C underabundance (based on a single line). Assuming that J20395358+4222505 is at the distance of Cyg OB2 we derive the radius from infrared photometry, finding R= 41.2 $\pm$ 4.0 R$_{\odot}$, log(L/L$_{\odot}$)= 5.71 $\pm$ 0.04 and a spectroscopic mass of 46.5 $\pm$ 15.0 M$_{\odot}$. The clumped mass-loss rate (clumping factor 10) is very high for the spectral type, $\dot{M}$ = 2.4x10$^{-6}$ M$_{\odot}$ a$^{-1}$. The high rotational velocity and mass-loss rate place the star at the hot side of the bi-stability jump. Together with the nearly solar CNO abundance pattern, they may also point to evolution in a binary system, J20395358+4222505 being the initial secondary.

preprint2021arXiv

Near-IR narrow-band imaging with CIRCE at the Gran Telescopio Canarias: Searching for Ly$α$-emitters at $z \sim 9.3$

Identifying very high-redshift galaxies is crucial for understanding the formation and evolution of galaxies. However, many questions still remain, and the uncertainty on the epoch of reionization is large. In this approach, some models allow a double-reionization scenario, although the number of confirmed detections at very high $z$ is still too low to serve as observational proof. The main goal of this project is studying whether we can search for Lyman-$α$ emitters (LAEs) at $z \sim 9$ using a narrow-band (NB) filter that was specifically designed by our team and was built for this experiment. We used the NB technique to select candidates by measuring the flux excess due to the Ly$α$ emission. The observations were taken with an NB filter (full width at half minimum of 11 nm and central wavelength $λ_{c} = 1.257 μ$m) and the CIRCE near-infrared camera for the GTC. We describe a data reduction procedure that was especially optimized to minimize instrumental effects. With a total exposure time of 18.3 hours, the final NB image covers an area of $\sim 6.7$ arcmin$^{2}$, which corresponds to a comoving volume of $1.1 \times 10^{3}$ Mpc$^{3}$ at $z = 9.3$. We pushed the source detection to its limit, which allows us to analyze an initial sample of 97 objects. We detail the different criteria we applied to select the candidates. The criteria included visual verifications in different photometric bands. None of the objects resembled a reliable LAE, however, and we found no robust candidate down to an emission-line flux of $2.9 \times 10^{-16}$ erg s$^{-1} $cm$^{-2}$, which corresponds to a Ly$α$ luminosity limit of $3 \times 10^{44}$ erg s$^{-1}$. We derive an upper limit on the Ly$α$ luminosity function at $z \sim 9$ that agrees well with previous constraints. We conclude that deeper and wider surveys are needed to study the LAE population at the cosmic dawn.

preprint2020arXiv

A detailed non-LTE analysis of LB-1: Revised parameters and surface abundances

LB-1 has recently been proposed to be a binary system at 4 kpc consisting of a B star of 8 Msol and a massive stellar black hole of 70 Msol. This finding challenges our current theories of massive star evolution and formation of BHs at solar metallicity. Our objective is to derive the effective temperature, surface gravity and chemical composition of the B-type component in order to determine its nature and evolutionary status and, indirectly, to constrain the mass of the BH. We use the non-LTE stellar atmosphere code FASTWIND to analyse new and archival high resolution data. We determine (Teff, logg) values of (14000$\pm500$ K, 3.50$\pm0.15$ dex) that, combined with the Gaia parallax, implies a spectroscopic mass, from logg, of $3.2^{+2.1}_{-1.9}$ Msol and an evolutionary mass, assuming single star evolution, of $5.2^{+0.3}_{-0.6}$ Msol. We determine an upper limit of 8 km/s for the projected rotational velocity and derive the surface abundances, finding the star to have a silicon abundance below solar, to be significantly enhanced in nitrogen and iron, and depleted in carbon and magnesium. Complementary evidence derived from a photometric extinction analysis and Gaia yields similar results for Teff and logg and a consistent distance around 2~kpc. We propose that the B star is a slightly evolved main sequence star of 3-5 Msol with surface abundances reminiscent of diffusion in late B/A chemically peculiar stars with low rotational velocities. There is also evidence for CN-processed material in its atmosphere. These conclusions rely critically on the distance inferred from the Gaia parallax. The goodness of fit of the Gaia astrometry also favours a high-inclination orbit. If the orbit is edge-on and the B star has a mass of 3-5 Msol, the mass of the dark companion would be 4-5 Msol, which would be easier to explain with our current stellar evolutionary models.

preprint2020arXiv

An ionised bubble powered by a proto-cluster at z = 6.5

We show herein that a proto-cluster of Ly$α$ emitting galaxies, spectroscopically confirmed at redshift 6.5, produces a remarkable number of ionising continuum photons. We start from the Ly$α$ fluxes measured in the spectra of the sources detected spectroscopically. From these fluxes we derive the ionising emissivity of continuum photons of the proto-cluster, which we compare with the ionising emissivity required to reionise the proto-cluster volume. We find that the sources in the proto-cluster are capable of ionising a large bubble, indeed larger than the volume occupied by the proto-cluster. For various calculations we have used the model AMIGA, in particular to derive the emissivity of the Lyman continuum photons required to maintain the observed volume ionised. Besides, we have assumed the ionising photons escape fraction given by AMIGA at this redshift.

preprint2020arXiv

Analysis of absorption lines in the high resolution spectra of five hot post-AGB candidates

From an analysis of absorption lines in the high resolution spectra we have derived the radial velocities, stellar parameters (Teff, gravity, wind-strength parameter logQ and projected rotational velocity) and abundances (C, N, O, and Si ) of IRAS 17460-3114, IRAS 18131-3008, IRAS 19336-0400, LSE 45 and LSE 163. Abundances are found to be solar, except for a low Si abundance in IRAS 19336-0400 and a mild CNO pattern in LSE 163, that rotates at an unusual large rotational velocity for its spectral classification. Combining the stellar parameters information with Gaia DR2 data we are able to derive absolute magnitudes, radii and luminosities and clarify the possible post-AGB nature of the objects. IRAS 17460-3114 and IRAS 18131-3008 are found to be massive OB stars, whereas IRAS 19336-0400 is found to be a post-AGB star, already showing nebular lines in the spectrum. However, we could not confirm the nature of LSE 45 and LSE 163 as post-AGB stars, although their parameters are much more inconsistent with those of massive stars. In both cases, we find a discrepancy between the spectroscopic mass and that derived from the predictions of post-AGB evolutionary tracks. In addition, LSE 45 lacks nebular lines, that are present in IRAS 19336-0400 at a similar temperature. In the case of LSE 163 the rotational velocity (259+/-15 km/s) would be extremely large for a star evolving to CSPN. The combination of this rotational velocity, the high Galactic latitude, slightly large radial velocity and mild CNO enhancement suggests a history of binary interaction.

preprint2020arXiv

Multiplicity of the red supergiant population in the young massive cluster NGC 330

The multiplicity properties of massive stars are one of the important outstanding issues in stellar evolution. Quantifying the binary statistics of all evolutionary phases is essential to paint a complete picture of how and when massive stars interact with their companions, and to determine the consequences of these interactions. We investigate the multiplicity of an almost complete census of red supergiant stars (RSGs) in NGC 330, a young massive cluster in the SMC. Using a combination of multi-epoch HARPS and MUSE spectroscopy, we estimate radial velocities and assess the kinematic and multiplicity properties of 15 RSGs in NGC 330. Radial velocities are estimated to better than +/-100 m/s for the HARPS data. The line-of-sight velocity dispersion for the cluster is estimated as 3.20 +0.69-0.52 km/s. When virial equilibrium is assumed, the dynamical mass of the cluster is log (M{dyn} /M{sun}) = 5.20+/-0.17, in good agreement with previous upper limits. We detect significant radial velocity variability in our multi-epoch observations and distinguish between variations caused by atmospheric activity and those caused by binarity. The binary fraction of NGC 330 RSGs is estimated by comparisons with simulated observations of systems with a range of input binary fractions. In this way, we account for observational biases and estimate the intrinsic binary fraction for RSGs in NGC 330 as f{RSG} = 0.3+/-0.1 for orbital periods in the range 2.3< log P [days] <4.3, with q>0.1. Using the distribution of the luminosities of the RSG population, we estimate the age of NGC 330 to be 45+/-5 Myr and estimate a red straggler fraction of 50%. We estimate the binary fraction of RSGs in NGC 330 and conclude that it appears to be lower than that of main-sequence massive stars, which is expected because interactions between an RSG and a companion are assumed to effectively strip the RSG envelope.

preprint2020arXiv

Properties of OB star-black hole systems derived from detailed binary evolution models

The recent gravitational wave measurements have demonstrated the existence of stellar mass black hole binaries. It is essential for our understanding of massive star evolution to identify the contribution of binary evolution to the formation of double black holes. A promising way to progress is investigating the progenitors of double black hole systems and comparing predictions with local massive star samples such as the population in 30 Doradus in the Large Magellanic Cloud (LMC). Methods. To this purpose, we analyse a large grid of detailed binary evolution models at LMC metallicity with initial primary masses between 10 and 40 Msun, and identify which model systems potentially evolve into a binary consisting of a black hole and a massive main sequence star. We then derive the observable properties of such systems, as well as peculiarities of the OB star component. We find that about 3% of the LMC late O and early B stars in binaries are expected to possess a black hole companion, when assuming stars with a final helium core mass above 6.6 M to form black holes. While the vast majority of them may be X-ray quiet, our models suggest that these may be identified in spectroscopic binaries, either by large amplitude radial velocity variations ( > 50 km s ) and simultaneous nitrogen surface enrichment, or through a moderate radial velocity ( > 10 km/s ) and simultaneously rapid rotation of the OB star. The predicted mass ratios are such that main sequence companions could be excluded in most cases. A comparison to the observed OB+WR binaries in the LMC, Be/X-ray binaries, and known massive BH binaries supports our conclusion. We expect spectroscopic observations to be able to test key assumptions in our models, with important implications for massive star evolution in general, and for the formation of double-black hole mergers in particular.

preprint2020arXiv

The IACOB project. VI. On the elusive detection of massive O-type stars close to the ZAMS

The apparent lack of massive O-type stars near the zero-age main sequence (at ages < 2 Myr) is a topic widely discussed. Different explanations for this elusive detection have been proposed, but no firm conclusions have been reached yet. We reassess this empirical result benefiting from the high-quality spectroscopic observations of >400 Galactic O-type stars gathered by the IACOB and OWN surveys. We used temperatures and gravities from a iacob-gbat/fastwind spectroscopic analysis to locate our sample in the Kiel and spectroscopic HR diagrams. We evaluated the completeness of our sample of stars, observational biases using information from the Galactic O star catalog (GOSC), systematics of our methodology, and compare with other recent studies using smaller samples of Galactic O-type stars. We base our discussion on the spectroscopic HR diagram to avoid the use of uncertain distances. We performed a detailed study of the young cluster Trumpler-14 as an example of how Gaia cluster distances can help to construct the associated classical HR diagram. The apparent lack of massive O-type stars near the ZAMS with masses between 30 and 70 Msol persist even when spectroscopic results from a large, non-biased sample of stars are used. We do not find correlation between the dearth of stars and observational biases, limitations of our methodology, or the use of spectroscopic HR diagram instead of the classical one. Investigating the efficiency of mass accretion during the formation process we conclude that an adjustment of the accretion rate towards lower values could reconcile the hotter boundary of detected O-type stars and the theoretical birthline. Last, we discuss that the presence of a small sample of O2-O3.5 stars found closer to the ZAMS might be explained taking into account non-standard star evolution (e.g. binary interaction, mergers, or homogeneous evolution).

preprint2018arXiv

Response to comment on "An excess of massive stars in the local 30 Doradus starburst"

Farr and Mandel reanalyse our data, finding initial-mass-function slopes for high mass stars in 30 Doradus that agree with our results. However, their reanalysis appears to underpredict the observed number of massive stars. Their technique results in more precise slopes than in our work, strengthening our conclusion that there is an excess of massive stars above $30\,\mathrm{M}_\odot$ in 30 Doradus.

preprint2016arXiv

Contributions of the Cherenkov Telescope Array (CTA) to the 6th International Symposium on High-Energy Gamma-Ray Astronomy (Gamma 2016)

List of contributions from the Cherenkov Telescope Array (CTA) Consortium presented at the 6th International Symposium on High-Energy Gamma-Ray Astronomy (Gamma 2016), July 11-15, 2016, in Heidelberg, Germany.

preprint2016arXiv

Early-type massive stars in Carina Nebula within the Gaia-ESO Survey

The Gaia-ESO Survey (GES) is obtaining high quality spectra of $\sim$10$^{5}$ stars in our Galaxy, providing an homogeneous and unique overview of all the main components of the Milky Way, its formation history and the evolution of young, mature and ancient Galactic populations. Our group is in charge of the early-type massive stars that define the youngest population in the survey. In this contribution, we present the results of the quantitative spectroscopic analysis of O-type stars in the Carina Nebula within the Gaia-ESO Survey. For this aim, we have used FASTWIND and CMFGEN stellar atmosphere codes, providing stellar parameters for the current sample (GES data release iDR4).

preprint2016arXiv

Probing the Dragonfish star-forming complex: the ionizing population of the young massive cluster Mercer 30

The Dragonfish Nebula has been recently claimed to be powered by a superluminous but elusive OB association. Instead, systematic searches in near-infrared photometric surveys have found many other cluster candidates on this sky region. Among these, the first confirmed young massive cluster was Mercer 30, where Wolf-Rayet stars were found. We perform a new characterization of Mercer 30 with unprecedented accuracy, combining NICMOS/HST and VVV photometric data with multi-epoch ISAAC/VLT H- and K-band spectra. Stellar parameters for most of spectroscopically observed cluster members are found through precise non-LTE atmosphere modeling with the CMFGEN code. Our spectrophotometric study for this cluster yields a new, revised distance of d = (12.4 +- 1.7) kpc and a total of Q = 6.70 x 10^50 Lyman ionizing photons. A cluster age of (4.0 +- 0.8) Myr is found through isochrone fitting, and a total mass of (1.6 +- 0.6) x 10^4 Msol is estimated thanks to our extensive knowledge of the post-main-sequence population. As a consequence, membership of Mercer 30 to the Dragonfish star-forming complex is confirmed, allowing us to use this cluster as a probe for the whole complex, which turns out to be extremely large (400 pc across) and located at the outer edge of the Sagittarius-Carina spiral arm (11 kpc from the Galactic Center). The Dragonfish complex hosts 19 young clusters or cluster candidates (including Mercer 30 and a new candidate presented in this work) and an estimated minimum of 9 field Wolf-Rayet stars. The sum of all these contributions accounts for, at least, 73% of the Dragonfish Nebula ionization and leaves little or no room for the alleged superluminous OB association; alternative explanations are discussed.

preprint2016arXiv

The Galactic O-Star Spectroscopic Survey (GOSSS). III. 142 additional O-type systems

This is the third installment of GOSSS, a massive spectroscopic survey of Galactic O stars, based on new homogeneous, high signal-to-noise ratio, R~2500 digital observations selected from the Galactic O-Star Catalog (GOSC). In this paper we present 142 additional stellar systems with O stars from both hemispheres, bringing the total of O-type systems published within the project to 590. Among the new objects there are 20 new O stars. We also identify 11 new double-lined spectroscopic binaries (SB2s), of which 6 are of O+O type and 5 of O+B type, and an additional new tripled-lined spectroscopic binary (SB3) of O+O+B type. We also revise some of the previous GOSSS classifications, present some egregious examples of stars erroneously classified as O-type in the past, introduce the use of luminosity class IV at spectral types O4-O5.5, and adapt the classification scheme to the work of Arias et al. (2016).

preprint2016arXiv

The IACOB project: III. New observational clues to understand macroturbulent broadening in massive O- and B-type stars

We aim to provide new empirical clues about macroturbulent spectral line broadening in O- and B-type stars to evaluate its physical origin. We use high-resolution spectra of ~430 stars with spectral types in the range O4-B9 (all luminosity classes). We characterize the line-broadening of adequate diagnostic metal lines using a combined FT and GOF technique. We perform a quantitative spectroscopic analysis of the whole sample using automatic tools coupled with a huge grid of FASTWIND models. We also incorporate quantitative information about line asymmetries to our observational description of the characteristics of the line-profiles, and present a comparison of the shape and type of line-profile variability found in a small sample of O stars and B supergiants with still undefined pulsational properties and B main sequence stars with variable line-profiles. We present a homogeneous and statistically significant overview of the (single snapshot) line-broadening properties of stars in the whole O and B star domain. We find empirical evidence of the existence of various types of non-rotational broadening agents acting in the realm of massive stars. Even though all of them could be quoted and quantified as a macroturbulent broadening from a practical point of view, their physical origin can be different. Contrarily to the early- to late-B dwarfs/giants, which present a mixture of cases in terms of line-profile shape and variability, the whole O-type and B supergiant domain (or, roughly speaking, stars with M_ZAMS > 15 M_sol) is fully dominated by stars with a remarkable non-rotational broadening component and very similar profiles (including type of variability). We provide some examples illustrating how this observational dataset can be used to evaluate scenarios aimed at explaining the existence of sources of non-rotational broadening in massive stars.

preprint2016arXiv

The IACOB project: IV. New predictions for high-degree non-radial mode instability domains in massive stars and connection with macroturbulent broadening

Asteroseismology is a powerful tool to access the internal structure of stars. Apart from the important impact of theoretical developments, progress in this field has been commonly associated with the analysis of time-resolved observations. Recently, the so-called macroturbulent broadening has been proposed to be a complementary and less expensive way -- in terms of observational time -- to investigate pulsations in massive stars. We assess to what extent this ubiquitous non-rotational broadening component shaping the line profiles of O stars and B supergiants is a spectroscopic signature of pulsation modes driven by a heat mechanism. We compute stellar main sequence and post-main sequence models from 3 to 70Msun with the ATON stellar evolution code and determine the instability domains for heat-driven modes for degrees l=1-20 using the adiabatic and non-adiabatic codes LOSC and MAD. We use the observational material presented in Simón-Díaz et al. (2016) to investigate possible correlations between the single snapshot line-broadening properties of a sample of ~260 O and B-type stars and their location inside/outside the various predicted instability domains. We present an homogeneous prediction for the non-radial instability domains of massive stars for degree l up to 20. We provide a global picture of what to expect from an observational point of view in terms of frequency range of excited modes, and investigate the behavior of the instabilities with stellar evolution and increasing degree of the mode. Furthermore, our pulsational stability analysis, once compared to the empirical results of Simón-Díaz et al. (2016), indicates that stellar oscillations originated by a heat mechanism can not explain alone the occurrence of the large non-rotational line-broadening component commonly detected in the O star and B supergiant domain.

preprint2016arXiv

The R136 star cluster dissected with Hubble Space Telescope/STIS. I. Far-ultraviolet spectroscopic census and the origin of HeII 1640 in young star clusters

We introduce a HST/STIS stellar census of R136a, the central ionizing star cluster of 30 Doradus. We present low resolution far-ultraviolet STIS/MAMA spectroscopy of R136 using 17 contiguous 52x0.2 arcsec slits which together provide complete coverage of the central 0.85 parsec (3.4 arcsec). We provide spectral types of 90% of the 57 sources brighter than m_F555W = 16.0 mag within a radius of 0.5 parsec of R136a1, plus 8 additional nearby sources including R136b (O4\,If/WN8). We measure wind velocities for 52 early-type stars from CIV 1548-51, including 16 O2-3 stars. For the first time we spectroscopically classify all Weigelt & Baier members of R136a, which comprise three WN5 stars (a1-a3), two O supergiants (a5-a6) and three early O dwarfs (a4, a7, a8). A complete Hertzsprung-Russell diagram for the most massive O stars in R136 is provided, from which we obtain a cluster age of 1.5+0.3_-0.7 Myr. In addition, we discuss the integrated ultraviolet spectrum of R136, and highlight the central role played by the most luminous stars in producing the prominent HeII 1640 emission line. This emission is totally dominated by very massive stars with initial masses above ~100 Msun. The presence of strong HeII 1640 emission in the integrated light of very young star clusters (e.g A1 in NGC 3125) favours an initial mass function extending well beyond a conventional upper limit of 100 Msun. We include montages of ultraviolet spectroscopy for LMC O stars in the Appendix. Future studies in this series will focus on optical STIS/CCD medium resolution observations.

preprint2015arXiv

ANTARES Constrains a Blazar Origin of Two IceCube PeV Neutrino Events

The source(s) of the neutrino excess reported by the IceCube Collaboration is unknown. The TANAMI Collaboration recently reported on the multiwavelength emission of six bright, variable blazars which are positionally coincident with two of the most energetic IceCube events. Objects like these are prime candidates to be the source of the highest-energy cosmic rays, and thus of associated neutrino emission. We present an analysis of neutrino emission from the six blazars using observations with the ANTARES neutrino telescope.The standard methods of the ANTARES candidate list search are applied to six years of data to search for an excess of muons --- and hence their neutrino progenitors --- from the directions of the six blazars described by the TANAMI Collaboration, and which are possibly associated with two IceCube events. Monte Carlo simulations of the detector response to both signal and background particle fluxes are used to estimate the sensitivity of this analysis for different possible source neutrino spectra. A maximum-likelihood approach, using the reconstructed energies and arrival directions of through-going muons, is used to identify events with properties consistent with a blazar origin.Both blazars predicted to be the most neutrino-bright in the TANAMI sample (1653$-$329 and 1714$-$336) have a signal flux fitted by the likelihood analysis corresponding to approximately one event. This observation is consistent with the blazar-origin hypothesis of the IceCube event IC14 for a broad range of blazar spectra, although an atmospheric origin cannot be excluded. No ANTARES events are observed from any of the other four blazars, including the three associated with IceCube event IC20. This excludes at a 90\% confidence level the possibility that this event was produced by these blazars unless the neutrino spectrum is flatter than $-2.4$.

preprint2015arXiv

B fields in OB stars (BOB): Detection of a strong magnetic field in the O9.7 V star HD54879

The number of magnetic stars detected among massive stars is small; nevertheless, the role played by the magnetic field in stellar evolution cannot be disregarded. Links between line profile variability, enhancements/depletions of surface chemical abundances, and magnetic fields have been identified for low-mass B-stars, but for the O-type domain this is almost unexplored. Based on FORS2 and HARPS spectropolarimetric data, we present the first detection of a magnetic field in HD54879, a single slowly rotating O9.7 V star. Using two independent and different techniques we obtained the firm detection of a surface average longitudinal magnetic field with a maximum amplitude of about 600 G, in modulus. A quantitative spectroscopic analysis of the star with the stellar atmosphere code FASTWIND results in an effective temperature and a surface gravity of 33000$\pm1000$ K and 4.0$\pm0.1$ dex. The abundances of carbon, nitrogen, oxygen, silicon, and magnesium are found to be slightly lower than solar, but compatible within the errors. We investigate line-profile variability in HD54879 by complementing our spectra with spectroscopic data from other recent OB-star surveys. The photospheric lines remain constant in shape between 2009 and 2014, although H$α$ shows a variable emission. The H$α$ emission is too strong for a standard O9.7 V and is probably linked to the magnetic field and the presence of circumstellar material. Its normal chemical composition and the absence of photospheric line profile variations make HD54879 the most strongly magnetic, non-variable single O-star detected to date.

preprint2015arXiv

B fields in OB stars (BOB): low-resolution FORS2 spectropolarimetry of the first sample of 50 massive stars

Within the context of the collaboration "B fields in OB stars (BOB)", we used the FORS2 low-resolution spectropolarimeter to search for a magnetic field in 50 massive stars, including two reference magnetic massive stars. Because of the many controversies of magnetic field detections obtained with the FORS instruments, we derived the magnetic field values with two completely independent reduction and analysis pipelines. We compare and discuss the results obtained from the two pipelines. We obtained a general good agreement, indicating that most of the discrepancies on magnetic field detections reported in the literature are caused by the interpretation of the significance of the results (i.e., 3-4 sigma detections considered as genuine, or not), instead of by significant differences in the derived magnetic field values. By combining our results with past FORS1 measurements of HD46328, we improve the estimate of the stellar rotation period, obtaining P = 2.17950+/-0.00009 days. For HD125823, our FORS2 measurements do not fit the available magnetic field model, based on magnetic field values obtained 30 years ago. We repeatedly detect a magnetic field for the O9.7V star HD54879, the HD164492C massive binary, and the He-rich star CPD -57 3509. We obtain a magnetic field detection rate of 6+/-4%, while by considering only the apparently slow rotators we derive a detection rate of 8+/-5%, both comparable with what was previously reported by other similar surveys. We are left with the intriguing result that, although the large majority of magnetic massive stars is rotating slowly, our detection rate is not a strong function of the stellar rotational velocity.

preprint2015arXiv

CTA Contributions to the 34th International Cosmic Ray Conference (ICRC2015)

List of contributions from the CTA Consortium presented at the 34th International Cosmic Ray Conference, 30 July - 6 August 2015, The Hague, The Netherlands.

preprint2015arXiv

Discovery of the massive overcontact binary VFTS 352: Evidence for enhanced internal mixing

The contact phase expected to precede the coalescence of two massive stars is poorly characterized due to the paucity of observational constraints. Here we report on the discovery of VFTS 352, an O-type binary in the 30 Doradus region, as the most massive and earliest spectral type overcontact system known to date. We derived the 3D geometry of the system, its orbital period $P_{\rm orb}=1.1241452(4)$ d, components' effective temperatures -- $T_1=42\,540\pm280$ K and $T_2=41\,120\pm290$ K -- and dynamical masses -- $M_1=28.63\pm0.30 M_{\odot}$ and $M_2 = 28.85\pm0.30 M_{\odot}$. Compared to single-star evolutionary models, the VFTS 352 components are too hot for their dynamical masses by about 2700 and 1100 K, respectively. These results can be explained naturally as a result of enhanced mixing, theoretically predicted to occur in very short-period tidally-locked systems. The VFTS 352 components are two of the best candidates identified so far to undergo this so-called chemically homogeneous evolution. The future of VFTS 352 is uncertain. If the two stars merge, a very rapidly rotating star will be produced. Instead, if the stars continue to evolve homogeneously and keep shrinking within their Roche Lobes, coalescence can be avoided. In this case, tides may counteract the spin down by winds such that the VFTS 352 components may, at the end of their life, fulfill the requirements for long gamma-ray burst (GRB) progenitors in the collapsar scenario. Independently of whether the VFTS 352 components become GRB progenitors, this scenario makes VFTS 352 interesting as a progenitor of a black hole binary, hence as a potential gravitational wave source through black hole-black hole merger.

preprint2015arXiv

Massive Stars: Some open questions and the role of Multi-object spectroscopy

Massive stars are key ingredients in the evolution of the Universe. Yet, important uncertainties and limits persist in our understanding of these objects, even in their early phases, limiting their application as tools to interpret the Universe. Here we review some of these open questions and argue that large samples are needed to answer them, both in the Milky Way and nearby galaxies. Multiobject spectroscopy plays a crucial role in this process.

preprint2015arXiv

On the formation and evolution of the first Be star in a black hole binary MWC 656

We find that the formation of MWC 656 (the first Be binary containing a black hole) involves a common envelope phase and a supernova explosion. This result supports the idea that a rapidly rotating Be star can emerge out of a common envelope phase, which is very intriguing because this evolutionary stage is thought to be too fast to lead to significant accretion and spin up of the B star. We predict $\sim 10-100$ of B BH binaries to currently reside in the Galactic disk, among which around $1/3$ contain a Be star, but there is only a small chance to observe a system with parameters resembling MWC 656. If MWC 656 is representative of intrinsic Galactic Be BH binary population, it may indicate that standard evolutionary theory needs to be revised. This would pose another evolutionary problem in understanding BH binaries, with BH X-ray Novae formation issue being the prime example. The future evolution of MWC 656 with a $\sim 5$ M$_{\odot}$ black hole and with a $\sim 13$ M$_{\odot}$ main sequence companion on a $\sim 60$ day orbit may lead to the formation of a coalescing BH-NS system. The estimated Advanced LIGO/Virgo detection rate of such systems is up to $\sim 0.2$ yr$^{-1}$. This empirical estimate is a lower limit as it is obtained with only one particular evolutionary scenario, the MWC 656 binary. This is only a third such estimate available (after Cyg X-1 and Cyg X-3), and it lends additional support to the existence of so far undetected BH--NS binaries.

preprint2015arXiv

Optical and X-ray early follow-up of ANTARES neutrino alerts

High-energy neutrinos could be produced in the interaction of charged cosmic rays with matter or radiation surrounding astrophysical sources. Even with the recent detection of extraterrestrial high-energy neutrinos by the IceCube experiment, no astrophysical neutrino source has yet been discovered. Transient sources, such as gamma-ray bursts, core-collapse supernovae, or active galactic nuclei are promising candidates. Multi-messenger programs offer a unique opportunity to detect these transient sources. By combining the information provided by the ANTARES neutrino telescope with information coming from other observatories, the probability of detecting a source is enhanced, allowing the possibility of identifying a neutrino progenitor from a single detected event. A method based on optical and X-ray follow-ups of high-energy neutrino alerts has been developed within the ANTARES collaboration. This program, denoted as TAToO, triggers a network of robotic optical telescopes (TAROT and ROTSE) and the Swift-XRT with a delay of only a few seconds after a neutrino detection, and is therefore well-suited to search for fast transient sources. To identify an optical or X-ray counterpart to a neutrino signal, the images provided by the follow-up observations are analysed with dedicated pipelines. A total of 42 alerts with optical and 7 alerts with X-ray images taken with a maximum delay of 24 hours after the neutrino trigger have been analysed. No optical or X-ray counterparts associated to the neutrino triggers have been found, and upper limits on transient source magnitudes have been derived. The probability to reject the gamma-ray burst origin hypothesis has been computed for each alert.

preprint2015arXiv

Search for muon-neutrino emission from GeV and TeV gamma-ray flaring blazars using five years of data of the ANTARES telescope

The ANTARES telescope is well-suited for detecting astrophysical transient neutrino sources as it can observe a full hemisphere of the sky at all times with a high duty cycle. The background due to atmospheric particles can be drastically reduced, and the point-source sensitivity improved, by selecting a narrow time window around possible neutrino production periods. Blazars, being radio-loud active galactic nuclei with their jets pointing almost directly towards the observer, are particularly attractive potential neutrino point sources, since they are among the most likely sources of the very high-energy cosmic rays. Neutrinos and gamma rays may be produced in hadronic interactions with the surrounding medium. Moreover, blazars generally show high time variability in their light curves at different wavelengths and on various time scales. This paper presents a time-dependent analysis applied to a selection of flaring gamma-ray blazars observed by the FERMI/LAT experiment and by TeV Cherenkov telescopes using five years of ANTARES data taken from 2008 to 2012. The results are compatible with fluctuations of the background. Upper limits on the neutrino fluence have been produced and compared to the measured gamma-ray spectral energy distribution.

preprint2015arXiv

Search of Dark Matter Annihilation in the Galactic Centre using the ANTARES Neutrino Telescope

A search for high-energy neutrinos coming from the direction of the Galactic Centre is performed using the data recorded by the ANTARES neutrino telescope from 2007 to 2012. The event selection criteria are chosen to maximise the sensitivity to possible signals produced by the self-annihilation of weakly interacting massive particles accumulated around the centre of the Milky Way with respect to the atmospheric background. After data unblinding, the number of neutrinos observed in the line of sight of the Galactic Centre is found to be compatible with background expectations. The 90% C.L. upper limits in terms of the neutrino+anti-neutrino flux, $\rm Φ_{ν_μ+\barν_μ}$, and the velocity averaged annihilation cross-section, $\rm <σ_{A}v>$, are derived for the WIMP self-annihilation channels into $\rm b\bar{b},W^{+}W^{-},τ^{+}τ^{-},μ^{+}μ^{-},ν\barν$. The ANTARES limits for $\rm <σ_{A}v>$ are shown to be the most stringent for a neutrino telescope over the WIMP masses $\rm 25\,GeV < M_{WIMP} < 10\,TeV$.

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 Brightest Young Star Clusters in NGC 5253

The nearby dwarf starburst galaxy NGC5253 hosts a number of young, massive star clusters, the two youngest of which are centrally concentrated and surrounded by thermal radio emission (the `radio nebula'). To investigate the role of these clusters in the starburst energetics, we combine new and archival Hubble Space Telescope images of NGC5253 with wavelength coverage from 1500 Ang to 1.9 micron in 13 filters. These include H-alpha, P-beta, and P-alpha, and the imaging from the Hubble Treasury Program LEGUS (Legacy Extragalactic UV Survey). The extraordinarily well-sampled spectral energy distributions enable modeling with unprecedented accuracy the ages, masses, and extinctions of the 9 optically brightest clusters (M_V < -8.8) and the two young radio nebula clusters. The clusters have ages ~1-15 Myr and masses ~1x10^4 - 2.5x10^5 M_sun. The clusters' spatial location and ages indicate that star formation has become more concentrated towards the radio nebula over the last ~15 Myr. The most massive cluster is in the radio nebula; with a mass 2.5x10^5 M_sun and an age ~1 Myr, it is 2-4 times less massive and younger than previously estimated. It is within a dust cloud with A_V~50 mag, and shows a clear nearIR excess, likely from hot dust. The second radio nebula cluster is also ~1 Myr old, confirming the extreme youth of the starburst region. These two clusters account for about half of the ionizing photon rate in the radio nebula, and will eventually supply about 2/3 of the mechanical energy in present-day shocks. Additional sources are required to supply the remaining ionizing radiation, and may include very massive stars.

preprint2015arXiv

The IACOB spectroscopic database: recent updates and first data release

The IACOB project is an ambitious long-term project which is contributing to step forward in our knowledge about the physical properties and evolution of Galactic massive stars. The project aims at building a large database of high-resolution, multi-epoch, spectra of Galactic OB stars, and the scientific exploitation of the database using state-of-the-art models and techniques. In this proceeding, we summarize the latest updates of the IACOB spectroscopic database and highlight some of the first scientific results from the IACOB project; we also announce the first data release and the first public version of the iacob-broad tool for the line-broadening characterization of OB-type spectra.

preprint2015arXiv

The VLT-FLAMES Tarantula Survey XIX. B-type Supergiants - Atmospheric parameters and nitrogen abundances to investigate the role of binarity and the width of the main sequence

TLUSTY non-LTE model atmosphere calculations have been used to determine atmospheric parameters and nitrogen (N) abundances for 34 single and 18 binary B-type supergiants (BSGs). The effects of flux contribution from an unseen secondary were considered for the binary sample. We present the first systematic study of the incidence of binarity for a sample of BSGs across the theoretical terminal age main sequence (TAMS). To account for the distribution of effective temperatures of the BSGs it may be necessary to extend the TAMS to lower temperatures. This is consistent with the derived distribution of mass discrepancies, projected rotational velocities (vsini) and N abundances, provided that stars cooler than this temperature are post RSG objects. For the BSGs in the Tarantula and previous FLAMES surveys, most have small vsini. About 10% have larger vsini (>100 km/s) but surprisingly these show little or no N enhancement. All the cooler BSGs have low vsini of <70km/s and high N abundance estimates, implying that either bi-stability braking or evolution on a blue loop may be important. A lack of cool binaries, possibly reflects the small sample size. Single star evolutionary models, which include rotation, can account for the N enhancement in both the single and binary samples. The detailed distribution of N abundances in the single and binary samples may be different, possibly reflecting differences in their evolutionary history. The first comparative study of single and binary BSGs has revealed that the main sequence may be significantly wider than previously assumed, extending to Teff=20000K. Some marginal differences in single and binary atmospheric parameters and abundances have been identified, possibly implying non-standard evolution for some of the sample. This sample as a whole has implications for several aspects of our understanding of the evolution of BSGs. Full abstract in paper

preprint2015arXiv

Time calibration with atmospheric muon tracks in the ANTARES neutrino telescope

The ANTARES experiment consists of an array of photomultipliers distributed along 12 lines and located deep underwater in the Mediterranean Sea. It searches for astrophysical neutrinos collecting the Cherenkov light induced by the charged particles, mainly muons, produced in neutrino interactions around the detector. Since at energies of $\sim$10 TeV the muon and the incident neutrino are almost collinear, it is possible to use the ANTARES detector as a neutrino telescope and identify a source of neutrinos in the sky starting from a precise reconstruction of the muon trajectory. To get this result, the arrival times of the Cherenkov photons must be accurately measured. A to perform time calibrations with the precision required to have optimal performances of the instrument is described. The reconstructed tracks of the atmospheric muons in the ANTARES detector are used to determine the relative time offsets between photomultipliers. Currently, this method is used to obtain the time calibration constants for photomultipliers on different lines at a precision level of 0.5 ns. It has also been validated for calibrating photomultipliers on the same line, using a system of LEDs and laser light devices.

preprint2014arXiv

A Be-type star with a black hole companion

Stellar-mass black holes have all been discovered through X-ray emission, which arises from the accretion of gas from their binary companions (this gas is either stripped from low-mass stars or supplied as winds from massive ones). Binary evolution models also predict the existence of black holes accreting from the equatorial envelope of rapidly spinning Be-type stars (stars of the Be type are hot blue irregular variables showing characteristic spectral emission lines of hydrogen). Of the ~80 Be X-ray binaries known in the Galaxy, however, only pulsating neutron stars have been found as companions. A black hole was formally allowed as a solution for the companion to the Be star MWC 656 (also known as HD 215227), although that was based on a single radial velocity curve of the Be star, a mistaken spectral classification and rough estimates of the inclination angle. Here we report observations of an accretion disk line mirroring the orbit of the Be star. This, together with an improved radial velocity curve of the Be star through fitting sharp Fe II profiles from the equatorial disk, and a refined Be classification (to that of a B1.5-B2 III star), reveals a black hole of 3.8 to 6.9 solar masses orbiting MWC 656, the candidate counterpart of the gamma-ray source AGL J2241+4454. The black hole is X-ray quiescent and fed by a radiatively inefficient accretion flow giving a luminosity less than 1.6 x 10-7 times the Eddington luminosity. This implies that Be binaries with black-hole companions are difficult to detect by conventional X-ray surveys.

preprint2014arXiv

Constraining the neutrino emission of gravitationally lensed Flat-Spectrum Radio Quasars with ANTARES data

This paper proposes to exploit gravitational lensing effects to improve the sensitivity of neutrino telescopes to the intrinsic neutrino emission of distant blazar populations. This strategy is illustrated with a search for cosmic neutrinos in the direction of four distant and gravitationally lensed Flat-Spectrum Radio Quasars. The magnification factor is estimated for each system assuming a singular isothermal profile for the lens. Based on data collected from 2007 to 2012 by the ANTARES neutrino telescope, the strongest constraint is obtained from the lensed quasar B0218+357, providing a limit on the total neutrino luminosity of this source of $1.08\times 10^{46}\,\mathrm{erg}\,\mathrm{s}^{-1}$. This limit is about one order of magnitude lower than those previously obtained in the ANTARES standard point source searches with non-lensed Flat-Spectrum Radio Quasars.

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

Detection of bridge emission above 50 GeV from the Crab pulsar with the MAGIC telescopes

The Crab pulsar is the only astronomical pulsed source detected at very high energy (VHE, E>100GeV) gamma-rays. The emission mechanism of VHE pulsation is not yet fully understood, although several theoretical models have been proposed. In order to test the new models, we measured the light curve and the spectra of the Crab pulsar with high precision by means of deep observations. We analyzed 135 hours of selected MAGIC data taken between 2009 and 2013 in stereoscopic mode. In order to discuss the spectral shape in connection with lower energies, 4.6 years of {\it Fermi}-LAT data were also analyzed. The known two pulses per period were detected with a significance of $8.0 σ$ and $12.6 σ$. In addition, significant emission was found between the two pulses with $6.2 σ$. We discovered the bridge emission above 50 GeV between the two main pulses. This emission can not be explained with the existing theories. These data can be used for testing new theoretical models.

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

Legacy ExtraGalactic UV Survey (LEGUS) with The Hubble Space Telescope. I. Survey Description

The Legacy ExtraGalactic UV Survey (LEGUS) is a Cycle 21 Treasury program on the Hubble Space Telescope, aimed at the investigation of star formation and its relation with galactic environment in nearby galaxies, from the scales of individual stars to those of ~kpc-size clustered structures. Five-band imaging, from the near-ultraviolet to the I-band, with the Wide Field Camera 3, plus parallel optical imaging with the Advanced Camera for Surveys, is being collected for selected pointings of 50 galaxies within the local 12 Mpc. The filters used for the observations with the Wide Field Camera 3 are: F275W(2,704 A), F336W(3,355 A), F438W(4,325 A), F555W(5,308 A), and F814W(8,024 A); the parallel observations with the Advanced Camera for Surveys use the filters: F435W(4,328 A), F606W(5,921 A), and F814W(8,057 A). The multi-band images are yielding accurate recent (<~50 Myr) star formation histories from resolved massive stars and the extinction-corrected ages and masses of star clusters and associations. The extensive inventories of massive stars and clustered systems will be used to investigate the spatial and temporal evolution of star formation within galaxies. This will, in turn, inform theories of galaxy evolution and improve the understanding of the physical underpinning of the gas-star formation relation and the nature of star formation at high redshift. This paper describes the survey, its goals and observational strategy, and the initial science results. Because LEGUS will provide a reference survey and a foundation for future observations with JWST and with ALMA, a large number of data products are planned for delivery to the community.

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

Masgomas-4: Physical characterization of a double-core obscured cluster with a massive and very young stellar population

The discovery of new, obscured massive star clusters has changed our understanding of the Milky Way star-forming activity from a passive to a very active star-forming machine. The search for these obscured clusters is strongly supported by the use of all-sky, near-IR surveys. The main goal of the MASGOMAS project is to search for and study unknown, young, and massive star clusters in the Milky Way, using near-IR data. Here we try to determine the main physical parameters (distance, size, total mass, and age) of Masgomas-4, a new double-core obscured cluster. Using near-IR photometry ($J$, $H$, and $K_S$) we selected a total of 21 stars as OB-type star candidates. Multi-object, near-IR follow-up spectroscopy allowed us to carry out the spectral classification of the OB-type candidates. Of the 21 spectroscopically observed stars, ten are classified as OB-type stars, eight as F- to early G-type dwarf stars, and three as late-type giant stars. Spectroscopically estimated distances indicate that the OB-type stars belong to the same cluster, located at a distance of $1.90^{+1.28}_{-0.90}$ kpc. Our spectrophotometric data confirm a very young and massive stellar population, with a clear concentration of pre-main-sequence massive candidates (Herbig Ae/Be) around one of the cluster cores. The presence of a surrounding HII cloud and the Herbig Ae/Be candidates indicate an upper age limit of 5 Myr.

preprint2014arXiv

MGB and the new Galactic O-Star Spectroscopic Survey spectral classification standard grid

In this poster we present three developments related to the Galactic O-Star Spectroscopic Survey (GOSSS). First, we are making public the first version of MGB, an IDL code that allows the user to compare oberved spectra to a grid of spectroscopic standards to measure spectral types, luminosity classes, rotation indexes, and spectral qualifiers. Second, we present the associated grid of standard stars for the spectral types O2 to O9.7, with several improvements over the original GOSSS grid of Sota et al. (2011). Third, we present a list of egregious classification errors in SIMBAD: stars that are or have been listed there as being of O type but that in reality are late-type stars.

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

Rotational velocities of single and binary O-type stars in the Tarantula Nebula

Rotation is a key parameter in the evolution of massive stars, affecting their evolution, chemical yields, ionizing photon budget, and final fate. We determined the projected rotational velocity, $v_e\sin i$, of $\sim$330 O-type objects, i.e. $\sim$210 spectroscopic single stars and $\sim$110 primaries in binary systems, in the Tarantula nebula or 30 Doradus (30\,Dor) region. The observations were taken using VLT/FLAMES and constitute the largest homogeneous dataset of multi-epoch spectroscopy of O-type stars currently available. The most distinctive feature of the $v_e\sin i$ distributions of the presumed-single stars and primaries in 30 Dor is a low-velocity peak at around 100\,$\rm{km s^{-1}}$. Stellar winds are not expected to have spun-down the bulk of the stars significantly since their arrival on the main sequence and therefore the peak in the single star sample is likely to represent the outcome of the formation process. Whereas the spin distribution of presumed-single stars shows a well developed tail of stars rotating more rapidly than 300\,$\rm{km s^{-1}}$, the sample of primaries does not feature such a high-velocity tail. The tail of the presumed-single star distribution is attributed for the most part -- and could potentially be completely due -- to spun-up binary products that appear as single stars or that have merged. This would be consistent with the lack of such post-interaction products in the binary sample, that is expected to be dominated by pre-interaction systems. The peak in this distribution is broader and is shifted toward somewhat higher spin rates compared to the distribution of presumed-single stars. Systems displaying large radial velocity variations, typical for short period systems, appear mostly responsible for these differences.

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

Searches for clustering in the time integrated skymap of the ANTARES neutrino telescope

This paper reports a search for spatial clustering of the arrival directions of high energy muon neutrinos detected by the ANTARES neutrino telescope. An improved two-point correlation method is used to study the autocorrelation of 3058 neutrino candidate events as well as cross-correlations with other classes of astrophysical objects: sources of high energy gamma rays, massive black holes and nearby galaxies. No significant deviations from the isotropic distribution of arrival directions expected from atmospheric backgrounds are observed.

preprint2014arXiv

The B Fields in OB Stars (BOB) Survey

The B fields in OB stars (BOB) survey is an ESO large programme collecting spectropolarimetric observations for a large number of early-type stars in order to study the occurrence rate, properties, and ultimately the origin of magnetic fields in massive stars. As of July 2014, a total of 98 objects were observed over 20 nights with FORS2 and HARPSpol. Our preliminary results indicate that the fraction of magnetic OB stars with an organised, detectable field is low. This conclusion, now independently reached by two different surveys, has profound implications for any theoretical model attempting to explain the field formation in these objects. We discuss in this contribution some important issues addressed by our observations (e.g., the lower bound of the field strength) and the discovery of some remarkable objects.

preprint2014arXiv

The IACOB project: II. On the scatter of O-dwarf spectral type -- effective temperature calibrations

We are now in an era of large spectroscopic surveys of OB-type stars. Quantitative spectroscopic analysis of these modern datasets is enabling us to review the physical properties of blue massive stars with robust samples, not only revisiting mean properties and general trends, but also incorporating information about the effects of second-order parameters. We investigate the spectral type -- effective temperature (SpT - Teff) calibration for O-type dwarfs, and its claimed dependence on metallicity, using statistically-meaningful samples of stars extracted from the IACOB and VFTS surveys. We perform a homogeneous differential spectroscopic analysis of 33 Galactic and 53 LMC O~dwarfs (spanning spectral types of O4 -- O9.7) using the IACOB-GBAT package, a $χ^2$-fitting algorithm based on a large pre-computed grid of FASTWIND models, and standard techniques for the hydrogen/helium analysis of O-type stars. We compare the estimated effective temperatures and gravities as a function of (internally consistent) spectral classifications. While the general trend is that the temperature of a star increases with earlier spectral types and decreasing metallicity, we show that the large range of gravities found for O-type dwarfs -- spaning up to 0.45-0.50 dex in some spectral bins -- plays a critical role on the dependence of the effective temperature calibrations as a function of spectral type and metallicity. This result warns us about the use of SpT - Teff calibrations for O-dwarfs which ignore the effects of gravity, and highlights the risks of employing calibrations based on small samples. The effects of this scatter in gravities (evolutionary status) for O-type dwarfs should be included in future recipes which employ SpT - Teff calibrations.

preprint2014arXiv

The VLT-FLAMES Tarantula Survey XVI. The optical+NIR extinction laws in 30 Doradus and the photometric determination of the effective temperatures of OB stars

Context: The commonly used extinction laws of Cardelli et al. (1989) have limitations that, among other issues, hamper the determination of the effective temperatures of O and early B stars from optical+NIR photometry. Aims: We aim to develop a new family of extinction laws for 30 Doradus, check their general applicability within that region and elsewhere, and apply them to test the feasibility of using optical+NIR photometry to determine the effective temperature of OB stars. Methods: We use spectroscopy and NIR photometry from the VLT-FLAMES Tarantula Survey and optical photometry from HST/WFC3 of 30 Doradus and we analyze them with the software code CHORIZOS using different assumptions such as the family of extinction laws. Results: We derive a new family of optical+NIR extinction laws for 30 Doradus and confirm its applicability to extinguished Galactic O-type systems. We conclude that by using the new extinction laws it is possible to measure the effective temperatures of OB stars with moderate uncertainties and only a small bias, at least up to E(4405-5495) ~ 1.5 mag.

preprint2014arXiv

The VLT-FLAMES Tarantula Survey XVII. Physical and wind properties of massive stars at the top of the main sequence

The evolution and fate of very massive stars (VMS) is tightly connected to their mass-loss properties. Their initial and final masses differ significantly as a result of mass loss. VMS have strong stellar winds and extremely high ionising fluxes, which are thought to be critical sources of both mechanical and radiative feedback in giant Hii regions. However, how VMS mass-loss properties change during stellar evolution is poorly understood. In the framework of the VLT-Flames Tarantula Survey (VFTS), we explore the mass-loss transition region from optically thin O to denser WNh star winds, thereby testing theoretical predictions. To this purpose we select 62 O, Of, Of/WN, and WNh stars, an unprecedented sample of stars with the highest masses and luminosities known. We perform a spectral analysis of optical VFTS as well as near-infrared VLT/SINFONI data using the non-LTE radiative transfer code CMFGEN to obtain stellar and wind parameters. For the first time, we observationally resolve the transition between optically thin O and optically thick WNh star winds. Our results suggest the existence of a kink between both mass-loss regimes, in agreement with recent MC simulations. For the optically thick regime, we confirm the steep dependence on the Eddington factor from previous theoretical and observational studies. The transition occurs on the MS near a luminosity of 10^6.1Lsun, or a mass of 80...90Msun. Above this limit, we find that - even when accounting for moderate wind clumping (with f = 0.1) - wind mass-loss rates are enhanced with respect to standard prescriptions currently adopted in stellar evolution calculations. We also show that this results in substantial helium surface enrichment. Based on our spectroscopic analyses, we are able to provide the most accurate ionising fluxes for VMS known to date, confirming the pivotal role of VMS in ionising and shaping their environments.

preprint2013arXiv

Astrophysical parameters of the peculiar X-ray transient IGR J11215-5952

The hard X-ray source IGR J11215-5952 is a peculiar transient, displaying very short X-ray outbursts every 165 days. We obtained high-resolution spectra of the optical counterpart, HD 306414, at different epochs, spanning a total of three months, before and around the 2007 February outburst with the combined aims of deriving its astrophysical parameters and searching for orbital modulation. We fit model atmospheres generated with the fastwind code to the spectrum. We also cross-correlated each individual spectrum to the best-fit model to derive radial velocities. From its spectral features, we classify HD 306414 as B0.5 Ia. From the model fit, we find Teff = 24 700 K and log g = 2.7, in good agreement with the morphological classification. Using the interstellar lines in its spectrum, we estimate a distance to HD 306414 d > 7 kpc. Assuming this distance, we derive R* = 40 Rsol and Mspect = 30 Msol (consistent, within errors, with Mevol = 38 Msol). Radial velocity changes are not dominated by the orbital motion, and we find an upper limit on the semi-amplitude for the optical component Kopt < 11 +- 6 km/s. Large variations in the depth and shape of photospheric lines suggest the presence of strong pulsations, which may be the main cause of the radial velocity changes. Very significant variations, uncorrelated with those of the photospheric lines are seen in the shape and position of the Halpha emission feature around the time of the X-ray outburst, but large excursions are also observed at other times. HD 306414 is a normal B0.5 Ia supergiant. Its radial velocity curve is dominated by an effect that is different from binary motion, and is most likely stellar pulsations. The data available suggest that the X-ray outbursts are caused by the close passage of the neutron star in a very eccentric orbit, perhaps leading to localised mass outflow. (abridged).

preprint2013arXiv

CTA contributions to the 33rd International Cosmic Ray Conference (ICRC2013)

Compilation of CTA contributions to the proceedings of the 33rd International Cosmic Ray Conference (ICRC2013), which took place in 2-9 July, 2013, in Rio de Janeiro, Brazil

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

Rotational properties of the O-type star population in the Tarantula region

The 30 Doradus (30\,Dor) region in the Large Magellanic Cloud (also known as the Tarantula Nebula) is the nearest massive starburst region, containing the richest sample of massive stars in the Local Group. It is the best possible laboratory to investigate aspects of the formation and evolution of massive stars. Here, we focus on rotation which is a key parameter in the evolution of these objects. We establish the projected rotational velocity, $v_{e}\sin i$, distribution of an unprecedented sample of 216 radial velocity constant ($\rm{ΔRV\, \leq\, 20 \,km s^{-1}}$) O-type stars in 30\,Dor observed in the framework of the VLT-FLAMES Tarantula Survey (VFTS). The distribution of $v_{e}\sin i$ shows a two-component structure: a peak around 80 $\rm{km s^{-1}}$ and a high-velocity tail extending up to $\sim$600 $\rm{km s^{-1}}$. Around 75% of the sample has 0 $\leq\, v_{e}\sin i \leq$ 200 $\rm{km s^{-1}}$ with the other 25% distributed in the high-velocity tail. The presence of the low-velocity peak is consistent with that found in other studies of late-O and early-B stars. The high-velocity tail is compatible with expectations from binary interaction synthesis models and may be predominantly populated by post-binary interaction, spun-up, objects and mergers. This may have important implications for the nature of progenitors of long-duration gamma ray bursts.

preprint2013arXiv

Spectroscopic and physical parameters of Galactic O-type stars. II. Observational constraints on projected rotational and extra broadening velocities as a function of fundamental parameters and stellar evolution

Rotation is of key importance for the evolution of hot massive stars, however, the rotational velocities of these stars are difficult to determine. Based on our own data for 31 Galactic O stars and incorporating similar data for 86 OB supergiants from the literature, we aim at investigating the properties of rotational and extra line-broadening as a function of stellar parameters and at testing model predictions about the evolution of stellar rotation. Fundamental stellar parameters were determined by means of the code FASTWIND. Projected rotational and extra broadening velocities originate from a combined Ft + GOF method. Model calculations published previously were used to estimate the initial evolutionary masses. The sample O stars with Minit > 50 Msun rotate with less that 26% of their break-up velocity, and they also lack objects with v sin i < 50 km/s. For the stars with Minit > 35 Msun on the hotter side of the bi-stability jump, the observed and predicted rotational rates agree quite well; for those on the cooler side of the jump, the measured velocities are systematically higher than the predicted ones. In general, the derived extra broadening velocities decrease toward cooler Teff, whilst for later evolutionary phases they appear, at the same v sin i, higher for high-mass stars than for low-mass ones. None of the sample stars shows extra broadening velocities higher than 110 km/s. For the majority of the more massive stars, extra broadening either dominates or is in strong competition with rotation. Conclusions: For OB stars of solar metallicity, extra broadening is important and has to be accounted for in the analysis. When appearing at or close to the zero-age main sequence, most of the single and more massive stars rotate slower than previously thought. Model predictions for the evolution of rotation in hot massive stars may need to be updated.

preprint2013arXiv

Target Selection for the Apache Point Observatory Galactic Evolution Experiment (APOGEE)

The Apache Point Observatory Galactic Evolution Experiment (APOGEE) is a high-resolution infrared spectroscopic survey spanning all Galactic environments (i.e., bulge, disk, and halo), with the principal goal of constraining dynamical and chemical evolution models of the Milky Way. APOGEE takes advantage of the reduced effects of extinction at infrared wavelengths to observe the inner Galaxy and bulge at an unprecedented level of detail. The survey's broad spatial and wavelength coverage enables users of APOGEE data to address numerous Galactic structure and stellar populations issues. In this paper we describe the APOGEE targeting scheme and document its various target classes to provide the necessary background and reference information to analyze samples of APOGEE data with awareness of the imposed selection criteria and resulting sample properties. APOGEE's primary sample consists of ~100,000 red giant stars, selected to minimize observational biases in age and metallicity. We present the methodology and considerations that drive the selection of this sample and evaluate the accuracy, efficiency, and caveats of the selection and sampling algorithms. We also describe additional target classes that contribute to the APOGEE sample, including numerous ancillary science programs, and we outline the targeting data that will be included in the public data releases.

preprint2013arXiv

The IACOB project: I. Rotational velocities in Northern Galactic O and early B-type stars revisited. The impact of other sources of line-broadening

Stellar rotation is an important parameter in the evolution of massive stars. Accurate and reliable measurements of projected rotational velocities in large samples of OB stars are crucial to confront the predictions of stellar evolutionary models with observational constraints. We reassess previous determinations of projected rotational velocities (vsini) in Galactic OB stars using a large, high quality spectroscopic dataset, and a strategy which account for other sources of broadening appart from rotation affecting the diagnostic lines We present a versatile and user friendly IDL tool, based on a combined Fourier Transform (FT) + goodness of fit (GOF) methodology, for the line-broadening characterization in OB-type stars. We use this tool to (a) investigate the impact of macroturbulent and microturbulent broadenings on vsini measurements, and (b) determine vsini in a sample of 200 Galactic OB-type stars, also characterizing the amount of macroturbulent broadening (\vmacro) affecting the line profiles. We present observational evidence illustrating the strengths and limitations of the proposed FT+GOF methodology for the case of OB stars. We confirm previous statements (based on indirect arguments or smaller samples) that the macroturbulent broadening is ubiquitous in the massive star domain. We compare the newly derived vsini with previous determinations not accounting for this extra line-broadening contribution, and show that those cases with vsini< 120 km/s need to be systematically revised downwards by ~25 (+/-20) km/s. We suggest that microturbulence may impose an upper limit below which vsini and \vmacro\ could be incorrectly derived by means of the proposed methodology as presently used, and discuss the implications of this statement on the study of relatively narrow line massive stars.

preprint2013arXiv

The VLT-FLAMES Tarantula Survey - XI. A census of the hot luminous stars and their feedback in 30 Doradus

We compile the first comprehensive census of hot luminous stars in the 30 Doradus (30 Dor) star forming region of the LMC. The census extends to a radius of 10arcmin (150pc) from the central cluster, R136. Stars were selected photometrically and combined with the latest spectral types. 1145 candidate hot luminous stars were identified of which >700 were considered genuine early type stars that contribute to feedback. We assess the spectroscopic completeness to be 85% in outer regions (>5pc) but fall to 35% in the vicinity of R136, giving a total of 500 hot luminous stars with spectroscopy. Stellar calibrations and models were used to obtain their physical parameters before integrated values were compared to global observations and the population synthesis code, Starburst99. The 31 W-R and Of/WN stars made large contributions to the total ionising and wind luminosities of ~40% and ~50%, respectively. Stars with Minit>100Msun also showed high contributions to the global feedback, ~25% in both cases. Such massive stars are not accounted for by the current Starburst99 code, which underestimated the ionising and wind luminosities of R136 by factors of ~2 and ~9, respectively. The census inferred a SFR of 0.073+/-0.04Msun/yr for 30 Dor, typically higher than results from popular SFR calibrations. However, it remained consistent with a far-UV luminosity tracer and a combined Halpha and mid-infrared tracer, but only after correcting for Halpha extinction. The global ionising output exceeded measurements from the associated gas and dust, suggesting ~6(+55/-6)% of ionising photons escape the region. When studying very luminous star forming regions, it is therefore essential to include the most massive stars to ensure a reliable energy budget. If 30 Dor is typical of other large star forming regions, estimates of the SFR will be underpredicted if this escape fraction is not accounted for.(abridged)

preprint2013arXiv

The VLT-FLAMES Tarantula Survey XII. Rotational velocities of the single O-type stars

Aims. Using ground based multi-object optical spectroscopy obtained in the framework of the VLT-FLAMES Tarantula Survey (VFTS), we aim to establish the (projected) rotational velocity distribution for a sample of 216 presumably single O-type stars in 30 Doradus (30 Dor). Methods. We measured projected rotational velocities, \vrot, by means of a Fourier transform method and a profile fitting method applied on a set of isolated spectral lines. We also used an iterative deconvolution procedure to infer the probability density, $\rm{P(\veq)}$, of the equatorial rotational velocity, \veq. Results. The distribution of \vrot\ shows a two-component structure: a peak around 80 \kms\ and a high-velocity tail extending up to $\sim$600 \kms. This structure is also present in the inferred distribution $\rm{P(\veq)}$ with around 80% of the sample having 0 $<$ \veq\, $\leq\, 300$ \kms\ and the other 20% distributed in the high-velocity region. Conclusions. Most of the stars in our sample rotate with a rate less than 20%\ of their break-up velocity. For the bulk of the sample, mass-loss in a stellar wind and/or envelope expansion is not efficient enough to significantly spin down these stars within the first few Myr of evolution. The presence of a sizeable population of fast rotators is compatible with recent population synthesis computations that investigate the influence of binary evolution on the rotation rate of massive stars. Despite the fact that we have excluded stars that show significant radial velocity variations, our sample may have remained contaminated by post-interaction binary products. The fact that the high-velocity tail may be preferentially (and perhaps even exclusively), populated by post-binary interaction products, has important implications for the evolutionary origin of systems that produce gamma-ray bursts.

preprint2013arXiv

The VLT-FLAMES Tarantula Survey. XIII: On the nature of O Vz stars in 30 Doradus

AIMS. We test the hypothesis of O Vz stars (characterized by having HeII4686 stronger in absorption than other He lines in their blue-violet spectra) being at a younger evolutionary stage than are normal O-type dwarfs. METHODS. We have performed a quantitative spectroscopic analysis of a sample of 38 O Vz and 46 O V stars, identified by the VLT-FLAMES Tarantula Survey in the 30 Doradus region of the Large Magellanic Cloud (LMC). We obtained the stellar and wind parameters of both samples using the FASTWIND stellar atmosphere code and the IACOB-GBAT grid-based automatic tool. In the framework of a differential study, we compared the physical and evolutionary properties of both samples, regarding Teff, logg, logQ and logL. We also investigated the predictions of the FASTWIND code about the O Vz phenomenon. RESULTS. We find a differential distribution of objects in terms of effective temperature, with O Vz stars dominant at intermediate values. The O Vz stars in 30 Doradus tend to be younger and less luminous, and they have weaker winds than the O V stars, but we also find examples with ages of 2-4 Myr and with luminosities and winds that are similar to those of normal O dwarfs. Moreover, the O Vz stars do not appear to have higher gravities than the O V stars. In addition to effective temperature and wind strength, our FASTWIND predictions indicate how important it is to take other stellar parameters (gravity and projected rotational velocity) into account for correctly interpreting the O Vz phenomenon. CONCLUSIONS. In general, the O Vz stars appear to be on or very close to the ZAMS, but there are some examples where the Vz classification does not necessarily imply extreme youth. In particular, the presence of O Vz stars in our sample at more evolved phases than expected is likely a consequence of modest O-star winds owing to the low-metallicity environment of the LMC.

preprint2013arXiv

Very Massive Stars (VMS) in the Local Universe

Recent studies have claimed the existence of very massive stars (VMS) up to 300 solar masses in the local Universe. As this finding may represent a paradigm shift for the canonical stellar upper-mass limit of 150 Msun, it is timely to discuss the status of the data, as well as the far-reaching implications of such objects. We held a Joint Discussion at the General Assembly in Beijing to discuss (i) the determination of the current masses of the most massive stars, (ii) the formation of VMS, (iii) their mass loss, and (iv) their evolution and final fate. The prime aim was to reach broad consensus between observers and theorists on how to identify and quantify the dominant physical processes.

preprint2012arXiv

A peculiar Of star in the Local Group galaxy IC 1613

Context. Results from the theory of radiatively driven winds are incorporated in stellar evolutionary and population synthesis models, and used in our interpretation of the observations of the deep Universe. Yet, the theory has been confirmed only until Small Magellanic Cloud metallicities. Analyses of O-stars at lower metallicities are needed to prove the theory. Aims. We have observed GHV-62024, an O6.5 IIIf star in the low-metallicity galaxy IC1613. According to a previous preliminary analysis this star could challenge the radiatively driven wind theory at low metallicities. Methods. Our observations were obtained with VIMOS at VLT, at R~2000 and were analysed using the latest version of the model atmosphere code FASTWIND, which includes N III Results. We obtain the stellar parameters and conclude that the star follows the average wind momentum-luminosity relationship (WLR) expected for its metallicity, but with a high value for the exponent of the wind velocity law, beta. We suggest that this high value may be reached because GHV-62024 could be a fast rotator seen at a low inclination angle. While the derived beta value does not change by adopting a lower wind terminal velocity, a wrong $V_\infty$ has a clear impact on the position of the star in the WLR diagram. The N and He abundances are very high, consistent with strong CNO mixing that could have been caused by the fast rotation, although we cannot discard a different origin. We find again the well-known mass-discrepancy. Conclusions. We conclude that the star follows the WLR expected for its metallicity. The results are consistent with GHV-62024 being a fast rotator seen close to pole-on, contaminated at the surface with CNO products and with a wind structure altered by the fast rotation without modifying the global WLR. We suggest that this could be a general property of fast rotators.

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 a young and massive stellar cluster: Spectrophotometric near-infrared study of Masgomas-1

Context: Recent near-infrared data have contributed to the discovery of new (obscured) massive stellar clusters and massive stellar populations in previously known clusters in our Galaxy. These discoveries lead us to view the Milky Way as an active star-forming machine. Aims: The main purpose of this work is to determine physically the main parameters (distance, size, total mass and age) of Masgomas-1, the first massive cluster discovered by our systematic search programme. Methods: Using near-infrared (J, H, and Ks) photometry we selected 23 OB-type and five red supergiant candidates for multi-object H- and K-spectroscopy and spectral classification. Results: Of the 28 spectroscopically observed stars, 17 were classified as OB-type, four as supergiants, one as an A-type dwarf star, and six as late-type giant stars. The presence of a supergiant population implies a massive nature of Masgomas-1, supported by our estimate of the cluster initial total mass of (1.94\pm0.28)\cdot10^4 M_{sun}, obtained after integrating of the cluster mass function. The distance estimate of 3.53 kpc locates the cluster closer than the Scutum--Centaurus base but still within that Galactic arm. The presence of an O9V star and red supergiants in the same population indicates that the cluster age is in the range of 8 to 10 Myr.

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

On the nature of candidate luminous blue variables in M33

Luminous blue variables (LBVs) are expected to play an important role in massive stellar evolution as well as being the progenitors of some of the most luminous supernovae known. In this paper we provide a multiwavelength study of the population of (candidate) LBVs identified within M33. New spectra provide an observational baseline of >4yr with respect to published data, which is well suited to identifying LBV outbursts. Multi-epoch optical and mid-IR surveys of M33 further constrain the variability of the sample and permit a search for dusty circumstellar ejecta. Spectroscopic and photometric variability appears common amongst the sample, although in many cases further observations will be required to determine its physical origin. Nevertheless, we report a new outburst of M33 Var C, while the transition of the WNLh star B517 to a cooler B supergiant phase between 1993-2010 confirms an LBV classification. Proof-of-concept quantitative analysis is provided for Romano's star; the results being consistent with the finding that its bolometric luminosity varies during its LBV excursions. The combination of the temperature and luminosity of two stars, the B hypergiant [HS80] 110A and the cool hypergiant B324, appears to be in violation of the empirical Humphreys-Davidson limit. Mid-IR observations demonstrate that a number of candidates appear associated with hot circumstellar dust, although no objects as extreme as Eta Carinae are identified. The combined multiwavelength dataset suggests that the population of LBVs studied is contaminated by stars demonstrating the B[e] phenomenon. Of these, a subset of optically faint, low luminosity stars associated with hot dust are of particular interest since they appear similar to the likely progenitors of SN 2008S and the 2008 NGC300 transient, albeit suffering less intrinsic extinction. [ABRIDGED]

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

Summary of IAU GA SpS5 - I. Obscured and distant clusters

This first part of Special Session 5 explored the current status of infrared-based observations of obscured and distant stellar clusters in the Milky Way galaxy. Recent infrared surveys, either serendipitously or using targeted searches, have uncovered a rich population of young and massive clusters. However, cluster characterization is more challenging as it must be obtained often entirely in the infrared due to high line-of-sight extinction. Despite this, much is to be gained through the identification and careful analysis of these clusters, as they allow for the early evolution of massive stars to be better constrained. Further, they act as beacons delineating the Milky Way's structure and as nearby, resolved analogues to the distant unresolved massive clusters studied in distant galaxies.

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 ARAUCARIA project: Grid-Based Quantitative Spectroscopic Study of Massive Blue Stars in NGC55

The quantitative study of the physical properties and chemical abundances of large samples of massive blue stars at different metallicities is a powerful tool to understand the nature and evolution of these objects. Their analysis beyond the Milky Way is challenging, nonetheless it is doable and the best way to investigate their behavior in different environments. Fulfilling this task in an objective way requires the implementation of automatic analysis techniques that can perform the analyses systematically, minimizing at the same time any possible bias. As part of the ARAUCARIA project we carry out the first quantitative spectroscopic analysis of a sample of 12 B-type supergiants in the galaxy NGC55 at 1.94 Mpc away. By applying the methodology developed in this work, we derive their stellar parameters, chemical abundances and provide a characterization of the present-day metallicity of their host galaxy. Based on the characteristics of the stellar atmosphere/line formation code FASTWIND, we designed and created a grid of models for the analysis of massive blue supergiant stars. Along with this new grid, we implemented a spectral analysis algorithm. Both tools were specially developed to perform fully consistent quantitative spectroscopic analyses of low spectral resolution of B-type supergiants in a fast and objective way. We present the main characteristics of our FASTWIND model grid and perform a number of tests to investigate the reliability of our methodology. The automatic tool is applied afterward to a sample of 12 B-type supergiant stars in NGC55, deriving the stellar parameters and abundances. The results indicate that our stars are part of a young population evolving towards a red supergiant phase. The derived chemical composition hints to an average metallicity similar to the one of the Large Magellanic Cloud, with no indication of a spatial trend across the galaxy.

preprint2012arXiv

The VLT-FLAMES Tarantula Survey VIII. Multiplicity properties of the O-type star population

Aims. We analyze the multiplicity properties of the massive O-type star population. With 360 O-type stars, this is the largest homogeneous sample of massive stars analyzed to date. Methods. We use multi-epoch spectroscopy and variability analysis to identify spectroscopic binaries. We also use a Monte-Carlo method to correct for observational biases. Results. We observe a spectroscopic binary fraction of 0.35\pm0.03, which corresponds to the fraction of objects displaying statistically significant radial velocity variations with an amplitude of at least 20km/s. We compute the intrinsic binary fraction to be 0.51\pm0.04. We adopt power-laws to describe the intrinsic period and mass-ratio distributions: f_P ~ (log P)^π (with 0.15 < log P < 3.5) and f_q ~ q^κ with 0.1 < q=M_2/M_1 < 1.0. The power-law indexes that best reproduce the observed quantities are π= -0.45 +/- 0.30 and κ= -1.0\pm0.4. The obtained period distribution thus favours shorter period systems compared to an Oepik law. The mass ratio distribution is slightly skewed towards low mass ratio systems but remains incompatible with a random sampling of a classical mass function. The binary fraction seems mostly uniform across the field of view and independent of the spectral types and luminosity classes. The binary fraction in the outer region of the field of view (r > 7.8', i.e. approx117 pc) and among the O9.7 I/II objects are however significantly lower than expected from statistical fluctuations. Conclusions. Using simple evolutionary considerations, we estimate that over 50% of the current O star population in 30 Dor will exchange mass with its companion within a binary system. This shows that binary interaction is greatly affecting the evolution and fate of massive stars, and must be taken into account to correctly interpret unresolved populations of massive stars.

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

Blue luminous stars in nearby galaxies - UIT005: a possible link to the luminous blue variable stage

(abridged) A detailed study of the blue supergiant UIT005 (B2-2.5Ia+) in M33 is presented. The results of our quantitative spectral analysis indicate that the star is a very luminous, log(L/Lsun)~5.9 dex, and massive, M~50 Msun, object, showing a very high nitrogen-to-oxygen ratio in its surface (N/O~8, by mass). Based on the derived Mg and Si abundances, we argue that this high N/O ratio cannot be the result of an initial low O content due to its location on the disk of M33, known to present a steep metallicity gradient. In combination with the He abundance, the most plausible interpretation is that UIT005 is in an advanced stage of evolution, showing in its surface N enrichment and O depletion resulting from mixing with CNO processed material from the stellar interior. A comparison with the predictions of current stellar evolutionary models indicates that there are significant discrepancies, in particular with regard to the degree of chemical processing, with the models predicting a much lower degree of O depletion than observed. At the same time, the mass-loss rate derived in our analysis is an order of magnitude lower than the values considered in the evolutionary calculations. Based on a study of the surrounding stellar population and the nearby cluster NGC588, using WFPC2 photometry, we suggest that UIT005 could be in fact a runaway star from this cluster.

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

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

Macroturbulent broadening: a single snap-shot approach to investigate pulsations in massive stars?

We present a brief progress report of our project aimed at the investigation of the so-called macroturbulent broadening affecting the line-profiles of O and early B-type stars, and speculate about the possible use of this spectroscopic feature as a single snap-shot approach to investigate pulsations in massive stars.

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

Near-infrared study of the stellar population of Sh2-152

Context: The discovery of new massive star clusters and massive stellar populations in previously known clusters in our Galaxy by means of infrared studies has changed our view of the Milky Way from an inactive to an active star-forming machine. Within this scenario, we present a near-infrared spectrophotometric study of the stellar content of the compact H II region Sh2-152. Aims: We aim to determine the distance, extinction, age, and mass of Sh2-152, using for the first time near-infrared stellar classification for several sources in the region. Methods: Using our near-infrared (J, H, and Ks) photometry and the colour-magnitude diagram for the cluster field, we selected 13 bright stars, candidate members of the reddened cluster's main sequence, for H- and K- spectroscopy and spectral classification. This near-infrared information was complemented with an optical spectrum of the ionizing central star to confirm its spectral nature. Results: From the 13 spectroscopically observed stars, 5 were classified as B-type, 3 as G-type, 2 were young stellar objects (YSOs), and 3 remained unclassified (because of the poor data quality). The cluster's extinction varies from A_{Ks}=0.5 to 2.6 magnitudes (A_V=4.5 to 24 magnitudes) and the distance is estimated to be 3.21\pm0.21 kpc. The age of the cluster is younger than 9.4 Myr and the lower limit to the total mass of the cluster is (2.45\pm0.79)\cdot10^3 M_{Sun}. We compare the number of ionizing photons emitted from the OB-type stars with the Lyman continuum photons derived from the radio observations and conclude that both quantities are consistent for the central region of Sh2-152. In contrast, the main ionizing source of the lower region remains unidentified.

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

Spectroscopic and physical parameters of Galactic O-type stars. I. Effects of rotation and spectral resolving power in the spectral classification of dwarfs and giants

The modern-era spectral classification of O-stars relies on either the Walborn or the Conti-Mathys scheme. Since both of these approaches have been developed using low-quality photographic data, their application to high-quality digital data might not be straightforward and be hampered by problems and complications that have not yet been appreciated. Using high-resolution spectra obtained with the ESO/MPG 2.2\,m telescope in La Silla and following the premises of the Walborn and Conti classification schemes, we determined the spectral types and luminosity classes of 19 Galactic O-type stars and compared them to those attributed by Walborn and Mathys based on low-quality data. Our analysis reveals that the morphological spectral types assigned using high-resolution data are systematically later (by up to 1.5 subtypes) then those attributed by Walborn. By means of line-profile simulations, we show that part of this discrepancy is more likely caused by the combined effect of stellar rotation and high spectral resolution on the depth of helium lines used as spectral type indicators. In addition, we demonstrate that at least for narrow-lined stars the "rotational effect" does not disappear when the high-resolution spectra are degraded to the resolution of the Walborn standards. We also find evidence of a systematic difference between our high-resolution quantitative spectral types and those assigned by Mathys. Rotation and spectral resolution are important third parameters in the spectral classification of O-type stars. To obtain reliable spectral classes within the Walborn approach, the unknown and the standard spectra must be compared at the same resolution and \vsini. Owing to resolution effects, the Conti approach might also need to be updated.

preprint2011arXiv

The FLAMES Tarantula Survey

The Tarantula survey is an ESO Large Programme which has obtained multi-epochs spectroscopy of over 800 massive stars in the 30 Dor region in the Large Magelanic Cloud. Here we briefly describe the main drivers of the survey and the observational material derived.

preprint2011arXiv

The IACOB project: A grid-based automatic tool for the quantitative spectroscopic analysis of O-stars

We present the IACOB grid-based automatic tool for the quantitative spectroscopic analysis of O-stars. The tool consists of an extensive grid of FASTWIND models, and a variety of programs implemented in IDL to handle the observations, perform the automatic analysis, and visualize the results. The tool provides a fast and objective way to determine the stellar parameters and the associated uncertainties of large samples of O-type stars within a reasonable computational time.

preprint2011arXiv

The IACOB project: synergies for the Gaia era

The IACOB spectroscopic survey of Galactic OB stars is an ambitious observational project aimed at compiling a large, homogeneous, high-resolution database of optical spectra of massive stars observable from the Northern hemisphere. The quantitative spectroscopic analysis of this database, complemented by the invaluable information provided by Gaia (mainly regarding photometry and distances), will represent a major step forward in our knowledge of the fundamental physical characteristics of Galactic massive stars. In addition, results from this analysis will be of interest for other scientific questions to be investigated using Gaia observations. In this contribution we outline the present status of the IACOB spectroscopic database and indicate briefly some of the synergy links between the IACOB and Gaia scientific projects.

preprint2011arXiv

The VLT-FLAMES Tarantula Survey I: Introduction and observational overview

The VLT-FLAMES Tarantula Survey (VFTS) is an ESO Large Programme that has obtained multi-epoch optical spectroscopy of over 800 massive stars in the 30 Doradus region of the Large Magellanic Cloud (LMC). Here we introduce our scientific motivations and give an overview of the survey targets, including optical and near-infrared photometry and comprehensive details of the data reduction. One of the principal objectives was to detect massive binary systems via variations in their radial velocities, thus shaping the multi-epoch observing strategy. Spectral classifications are given for the massive emission-line stars observed by the survey, including the discovery of a new Wolf-Rayet star (VFTS 682, classified as WN5h), 2' to the northeast of R136. To illustrate the diversity of objects encompassed by the survey, we investigate the spectral properties of sixteen targets identified by Gruendl & Chu from Spitzer photometry as candidate young stellar objects or stars with notable mid-infrared excesses. Detailed spectral classification and quantitative analysis of the O- and B-type stars in the VFTS sample, paying particular attention to the effects of rotational mixing and binarity, will be presented in a series of future articles to address fundamental questions in both stellar and cluster evolution.

preprint2011arXiv

The VLT-FLAMES Tarantula Survey II: R139 revealed as a massive binary system

We report the discovery that R139 in 30 Doradus is a massive spectroscopic binary system. Multi-epoch optical spectroscopy of R139 was obtained as part of the VLT-FLAMES Tarantula Survey, revealing a double-lined system. The two components are of similar spectral types; the primary exhibits strong C III 4650 emission and is classified as an O6.5 Iafc supergiant, while the secondary is an O6 Iaf supergiant. The radial-velocity variations indicate a highly eccentric orbit with a period of 153.9 days. Photometry obtained with the Faulkes Telescope South shows no evidence for significant variability within an 18 month period. The orbital solution yields lower mass limits for the components of M1sin^3 i = 78 \pm 8 Msun and M2sin^3 i = 66 \pm 7 Msun. As R139 appears to be the most massive binary system known to contain two evolved Of supergiants, it will provide an excellent test for atmospheric and evolutionary models.

preprint2011arXiv

The VLT-FLAMES Tarantula Survey III: A very massive star in apparent isolation from the massive cluster R136

VFTS 682 is located in an active star-forming region, at a projected distance of 29 pc from the young massive cluster R136 in the Tarantula Nebula of the Large Magellanic Cloud. It was previously reported as a candidate young stellar object, and more recently spectroscopically revealed as a hydrogen-rich Wolf-Rayet (WN5h) star. Our aim is to obtain the stellar properties, such as its intrinsic luminosity, and to investigate the origin of VFTS 682. To this purpose, we model optical spectra from the VLT-FLAMES Tarantula Survey with the non-LTE stellar atmosphere code CMFGEN, as well as the spectral energy distribution from complementary optical and infrared photometry. We find the extinction properties to be highly peculiar (RV ~4.7), and obtain a surprisingly high luminosity log(L/Lsun) = 6.5 \pm 0.2, corresponding to a present-day mass of ~150Msun. The high effective temperature of 52.2 \pm 2.5kK might be explained by chemically homogeneous evolution - suggested to be the key process in the path towards long gamma-ray bursts. Lightcurves of the object show variability at the 10% level on a timescale of years. Such changes are unprecedented for classical Wolf-Rayet stars, and are more reminiscent of Luminous Blue Variables. Finally, we discuss two possibilities for the origin of VFTS 682: (i) the star either formed in situ, which would have profound implications for the formation mechanism of massive stars, or (ii) VFTS 682 is a slow runaway star that originated from the dense cluster R136, which would make it the most massive runaway known to date.

preprint2010arXiv

Are the Stellar Winds in IC1613 stronger than expected?

In this poster we present the results of our analyses of three early massive stars in IC 1613, whose spectra have been observed with VIMOS and analyzed with CMFGEN and FASTWIND. One of the targets resulted a possible LBV and the other two are Of stars with unexpectedly strong winds. The Of stars seem to be strongly contaminated by CNO products. Our preliminary results may represent a challenge for the theory of stellar atmospheres, but they still have to be confirmed by the analysis of more objects and a more complete coverage of the parameter space.

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

Is macroturbulence in OB Sgs related to pulsations?

As part of a long term observational project, we are investigating the macroturbulent broadening in O and B supergiants (Sgs) and its possible connection with spectroscopic variability phenomena and stellar oscillations. We present the first results of our project, namely firm observational evidence for a strong correlation between the extra broadening and photospheric line-profile variations in a sample of 13 Sgs with spectral types ranging from O9.5 to B8.

preprint2010arXiv

Is macroturbulent broadening in OB Supergiants related to pulsations?

The spectrum of O and B Supergiants is known to be affected by an important extra line-broadening (usually called macroturbulence) that adds to stellar rotation. Recent analysis of high resolution spectra has shown that the interpretation of this line-broadening as a consequence of large-scale turbulent motions would imply highly super-sonic velocity fields, making this scenario quite improbable. Stellar oscillations have been proposed as a likely alternative explanation. We present first encouraging results of an observational project aimed at investigating the $macroturbulent$ broadening in O and B Supergiants, and its possible connection with spectroscopic variability phenomena and stellar oscillations: a) all the studied B Supergiants show line profile variations, quantified by means of the first (< v >) and third velocity (< v^3 >) moments of the lines, b) there is a strong correlation between the peak-to-peak amplitudes of the < v > and < v^3 > variability and the size of the extra-broadening.

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

Near-infrared spectroscopy in NGC 7538

The characterisation of the stellar population toward young high-mass star-forming regions allows to constrain fundamental cluster properties like distance and age. These are essential when using high-mass clusters as probes to conduct Galactic studies. NGC 7538 is a star-forming region with an embedded stellar population only unearthed in the near-infrared. We present the first near-infrared spectro-photometric study of the candidate high-mass stellar content in NGC 7538. We obtained H and K spectra of 21 sources with both the multi-object and long-slit modes of LIRIS at the WHT, and complement these data with sub-arcsecond JHKs photometry of the region using the imaging mode of the same instrument. We find a wide variety of objects within the studied stellar population of NGC 7538. Our results discriminate between a stellar population associated to the HII region, but not contained within its extent, and several pockets of more recent star formation. We report the detection of CO bandhead emission toward several sources as well as other features indicative of a young stellar nature. We infer a spectro-photometric distance of 2.7+-0.5 kpc, an age spread in the range 0.5-2.2 Myr and a total mass ~1.7x10^3 Msun for the older population.

preprint2010arXiv

New Optical Results on gamma-ray Binaries

We present new optical spectroscopy of the gamma-ray binary LS 5039. Our data show evidence for sub-orbital modulation in the radial velocities with amplitude ~7 km/s and period ~Porb/4. This short-term oscillation is stable over at least 7 years and it is likely triggered by non-radial oscillations of the O6.5V optical star. We also present the results of a spectroscopic campaign on MWC 148, the optical counterpart of the new gamma-ray binary candidate HESS J0632+067. Long-term variations in the Halpha and Hbeta emission line parameters are clearly detected which, if modulated with the binary orbit, would imply a period >200 days.

preprint2010arXiv

Observational evidence for a correlation between macroturbulent broadening and line-profile variations in OB Supergiants

The spectra of O and B supergiants are known to be affected by a significant form of extra line broadening (usually referred to as macroturbulence) in addition to that produced by stellar rotation. Recent analyses of high resolution spectra have shown that the interpretation of this line broadening as a consequence of large scale turbulent motions would imply highly supersonic velocity fields in photospheric regions, making this scenario quite improbable. Stellar oscillations have been proposed as a likely alternative explanation. As part of a long term observational project, we are investigating the macroturbulent broadening in O and B supergiants and its possible connection with spectroscopic variability phenomena and stellar oscillations. In this letter, we present the first encouraging results of our project, namely firm observational evidence for a strong correlation between the extra broadening and photospheric line-profile variations in a sample of 13 supergiants with spectral types ranging from O9.5 to B8.

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 IACOB spectroscopic database of Galactic OB stars

We present the IACOB spectroscopic database, the largest homogeneous database of high-resolution, high signal-to-noise ratio spectra of Northern Galactic OB-type stars compiled up to date. The spectra were obtained with the FIES spectrograph attached to the Nordic Optical Telescope. We briefly summarize the main characeristics and present status of the IACOB, first scientific results, and some future plans for its extension and scientific exploitation.

preprint2010arXiv

The IACOB spectroscopic database of Northern Galactic OB stars

We present the IACOB spectroscopic database, an homogeneous set of high quality, high resolution spectra of Galactic O- and B-type stars obtained with the FIES spectrograph attached to the Nordic Optical Telescope. We also present some results from ongoing projects using the IACOB database.

preprint2010arXiv

The nature of V39: an LBV candidate or LBV impostor in the very low metallicity galaxy IC 1613?

[abridged] Context: Very few examples of luminous blue variable (LBV) stars or LBV candidates (LBVc) are known, particularly at metallicities below the SMC. The LBV phase is crucial for the evolution of massive stars, and its behavior with metallicity is poorly known. V39 in IC 1613 is a well-known photometric variable, with B-band changes larger than 1mag. over its period. The star, previously proposed to be a projection of a Galactic W Virginis and an IC 1613 red supergiant, shows features that render it a possible LBVc. Method: We investigate mid-resolution blue and red VLT-VIMOS spectra of V39, covering a time span of 40 days, and perform a quantitative analysis of the combined spectrum using the model atmosphere code CMFGEN. Results: We identify strong Balmer and FeII P-Cygni profiles, and a hybrid spectrum resembling a B-A supergiant in the blue and a G-star in the red. No significant Vrad variations are detected, and the spectral changes are small over the photometric period. Our analysis places V39 in the low-luminosity part of the LBV and LBVc region, but it is also consistent with a sgB[e] star. Conclusions: The radial velocity indicates that V39 belongs to IC 1613. The lack of Vrad changes and spectroscopic variations excludes binary scenarios. The features observed are not consistent with a W Virginis star, and this possibility is also discarded. We propose that the star is a B-A LBVc or sgB[e] star surrounded by a thick disk precessing around it. If confirmed, V39 would be the lowest metallicity resolved LBV candidate known to date. Alternatively, it could represent a new transient phase of massive star evolution, an LBV impostor.

preprint2009arXiv

The stellar population of the star forming region G61.48+0.09

- Context: We present the results of a near-infrared photometric and spectroscopic study of the star forming region G61.48+0.09. - Aims: The purpose of this study is to characterize the stellar content of the cluster and to determine its distance, extinction, age and mass. - Methods: The stellar population was studied by using color-magnitude diagrams to select twenty promising cluster members, for which follow up spectroscopy was done. The observed spectra allowed a spectral classification of the stars. - Results: Two stars have emission lines, twelve are G-type stars, and six are late-O or early-B stars. - Conclusions: The cluster's extinction varies from A_{K_S} = 0.9 to A_{K_S} = 2.6, (or A_{V}~8 to A_{V}~23). G61.48+0.09 is a star forming region located at 2.5+/-0.4 Kpc. The cluster is younger than 10 Myr and has a minimum stellar mass of 1500+/-500 Solar masses. However, the actual total mass of the cluster remains undetermined, as we cannot see its whole stellar content.

preprint2009arXiv

The Ultraviolet Spectrum and Physical Properties of the Mass Donor Star in HD 226868 = Cygnus X-1

We present an examination of high resolution, ultraviolet spectroscopy from Hubble Space Telescope of the photospheric spectrum of the O-supergiant in the massive X-ray binary HD 226868 = Cyg X-1. We analyzed this and ground-based optical spectra to determine the effective temperature and gravity of the O9.7 Iab supergiant. Using non-local thermodynamic equilibrium (non-LTE), line blanketed, plane parallel models from the TLUSTY grid, we obtain T_eff = 28.0 +/- 2.5kK and log g > 3.00 +/- 0.25, both lower than in previous studies. The optical spectrum is best fit with models that have enriched He and N abundances. We fit the model spectral energy distribution for this temperature and gravity to the UV, optical, and IR fluxes to determine the angular size of and extinction towards the binary. The angular size then yields relations for the stellar radius and luminosity as a function of distance. By assuming that the supergiant rotates synchronously with the orbit, we can use the radius - distance relation to find mass estimates for both the supergiant and black hole as a function of the distance and the ratio of stellar to Roche radius. Fits of the orbital light curve yield an additional constraint that limits the solutions in the mass plane. Our results indicate masses of 23^{+8}_{-6} M_sun for the supergiant and 11^{+5}_{-3} M_sun for the black hole.

preprint2009arXiv

The VLT-FLAMES Tarantula Survey

The Tarantula Survey is an ambitious ESO Large Programme that has obtained multi-epoch spectroscopy of over 1,000 massive stars in the 30 Doradus region of the Large Magellanic Cloud. Here we introduce the scientific motivations of the survey and give an overview of the observational sample. Ultimately, quantitative analysis of every star, paying particular attention to the effects of rotational mixing and binarity, will be used to address fundamental questions in both stellar and cluster evolution.

preprint2009arXiv

The Young Stellar Population of IC1613. I. A New Catalogue of OB Associations

Context: Determining the parameters of massive stars is crucial to understand many processes in galaxies and the Universe, since these objects are important sources of ionization, chemical enrichment and momentum. 10m class telescopes enable us to perform detailed quantitative spectroscopic analyses of massive stars in other galaxies, sampling areas of different metallicity. Relating the stars to their environment is crucial to understand the physical processes ruling their formation and evolution. Aims: In preparation for the GTC, our goal is to build a catalogue of massive star candidates in the metal-poor irregular galaxy IC1613 with high astrometric accuracy, apt for the current generation of multi-object spectrographs. A census of OB associations in this galaxy is also needed, to provide important additional information about age and environment of the candidate OB stars. Methods: From INT-WFC observations, we have built an astrometric and photometric catalogue of stars in IC1613. Candidate blue massive stars are preselected from their colors. A friends-of-friends algorithm is developed to find their clustering in the galaxy. While a common physical origin for all the members of the associations cannot be ensured, this is a necessary first step to place candidate OB stars in a population context. Results: We have produced a deep catalogue of targets in IC1613 that covers a large field of view. To achieve high astrometric accuracy a new astrometric procedure is developed for the INT-WFC data. We have also built a catalogue of OB associations in IC1613. We have found that they concentrate in the central regions, specially in the HII bubbles. The study of extinction confirms that it is patchy, with local values of color-excess above the foreground value.

preprint2009arXiv

Time-scales of Line-broadening Variability in OB Supergiants

Several works have recently shown that there is an important extra line-broadening (usually called $macroturbulence$) affecting the spectra of O and B Supergiants that adds to stellar rotation. So far, the only (very recent) physical explanation for the appearance of $macroturbulence$ relates to oscillations. This is a plausible explanation, but no direct evidence confirming its validity has been presented yet. We recently started an observational project to obtain constraints on the time-scales of variability associated to this extra line-broadening and its possible origin. Our observational strategy consists of the study of a well selected group of O and B stars, for which we obtain time series of high-quality spectra. We present some preliminary results from our first campaign with FIES@NOT2.5m.

preprint2007arXiv

Fourier method in the determination of rotational velocities in OB stars

We present a comprehensive study that applies the Fourier transform to a sample of O and early B-type stars (either dwarfs, giants, or supergiants) to determine their projected rotational velocities, compare with previous values obtained with other methods, and seek for evidence of extra broadening in the spectral lines The Fourier technique, extensively used in the study of cooler stars, has remained only marginally applied for the case of early-type stars. The comparison of \vsini values obtained through the \ft and \fwhm methods shows that the \fwhm technique must be used with care in the analysis of OB giants and supergiants, and when it is applied to \ion{He}{i} lines. Contrarily, the \ft method appears to be a powerful tool to derive reliable projected rotational velocities, and separate the effect of rotation from other broadening mechanisms present in these stars. The analysis of the sample of OB stars shows that while dwarfs and giants display a broad range of projected rotational velocities, from less than 30 up to 450 \kms, supergiants have in general values close to or below 100 \kms. The analysis has also definitely shown that while the effect of extra broadening is negligible in OB dwarfs, it is clearly present in supergiants. When examining the behavior of the projected rotational velocities with the stellar parameters and across the HR diagram, we conclude, in agreement with previous researchers, that the rotational velocity should decrease when the stars evolve. On the contrary, macroturbulence may be constant, resulting therefore in an increasing importance as compared to rotation when the stars evolve.

preprint2005arXiv

A possible black hole in the gamma-ray microquasar LS 5039

The population of high energy and very high energy gamma-ray sources, detected with EGRET and the new generation of ground-based Cherenkov telescopes, conforms a reduced but physically important sample. Most of these sources are extragalactic (e.g., blazars), while among the galactic ones there are pulsars and SN remnants. The microquasar LS 5039, previously proposed to be associated with an EGRET source by Paredes et al. (2000), has recently been detected at TeV energies, confirming that microquasars should be regarded as a class of high energy gamma-ray sources. To model and understand how the energetic photons are produced and escape from LS 5039 it is crucial to unveil the nature of the compact object, which remains unknown. Here we present new intermediate-dispersion spectroscopy of this source which, combined with values reported in the literature, provides an orbital period of 3.90603+/-0.00017 d, a mass function f(M)=0.0053+/-0.0009 M_sun, and an eccentricity e=0.35+/-0.04. Atmosphere model fitting to the spectrum of the optical companion, together with our new distance estimate of d=2.5+/-0.1 kpc, yields R_opt=9.3+0.7-0.6 R_sun, log (L_opt/L_sun)=5.26+/-0.06, and M_opt=22.9+3.4-2.9 M_sun. These, combined with our dynamical solution and the assumption of pseudo-synchronization, yield an inclination i=24.9+/-2.8 degree and a compact object mass M_X=3.7+1.3-1.0 M_sun. This is above neutron star masses for most of the standard equations of state and, therefore, we propose that the compact object in LS 5039 is a black hole. We finally discuss about the implications of our orbital solution and new parameters of the binary system on the CNO products, the accretion/ejection energetic balance, the SN explosion scenario, and the behaviour of the TeV emission with the new orbital period.

A. Herrero

What is connected

Connect this record

See the researcher in context

Building this map preview

129 published item(s)

Populations of evolved massive binary stars in the Small Magellanic Cloud II: Predictions from rapid binary evolution

The Gaia-ESO Survey: The analysis of the hot-star spectra

The nature of the Cygnus extreme B-supergiant 2MASS J20395358+4222505

Near-IR narrow-band imaging with CIRCE at the Gran Telescopio Canarias: Searching for Ly$α$-emitters at $z \sim 9.3$

A detailed non-LTE analysis of LB-1: Revised parameters and surface abundances

An ionised bubble powered by a proto-cluster at z = 6.5

Analysis of absorption lines in the high resolution spectra of five hot post-AGB candidates

Multiplicity of the red supergiant population in the young massive cluster NGC 330

Properties of OB star-black hole systems derived from detailed binary evolution models

The IACOB project. VI. On the elusive detection of massive O-type stars close to the ZAMS

Response to comment on "An excess of massive stars in the local 30 Doradus starburst"

Contributions of the Cherenkov Telescope Array (CTA) to the 6th International Symposium on High-Energy Gamma-Ray Astronomy (Gamma 2016)

Early-type massive stars in Carina Nebula within the Gaia-ESO Survey

Probing the Dragonfish star-forming complex: the ionizing population of the young massive cluster Mercer 30

The Galactic O-Star Spectroscopic Survey (GOSSS). III. 142 additional O-type systems

The IACOB project: III. New observational clues to understand macroturbulent broadening in massive O- and B-type stars

The IACOB project: IV. New predictions for high-degree non-radial mode instability domains in massive stars and connection with macroturbulent broadening

The R136 star cluster dissected with Hubble Space Telescope/STIS. I. Far-ultraviolet spectroscopic census and the origin of HeII 1640 in young star clusters

ANTARES Constrains a Blazar Origin of Two IceCube PeV Neutrino Events

B fields in OB stars (BOB): Detection of a strong magnetic field in the O9.7 V star HD54879

B fields in OB stars (BOB): low-resolution FORS2 spectropolarimetry of the first sample of 50 massive stars

CTA Contributions to the 34th International Cosmic Ray Conference (ICRC2015)

Discovery of the massive overcontact binary VFTS 352: Evidence for enhanced internal mixing

Massive Stars: Some open questions and the role of Multi-object spectroscopy

On the formation and evolution of the first Be star in a black hole binary MWC 656

Optical and X-ray early follow-up of ANTARES neutrino alerts

Search for muon-neutrino emission from GeV and TeV gamma-ray flaring blazars using five years of data of the ANTARES telescope

Search of Dark Matter Annihilation in the Galactic Centre using the ANTARES Neutrino Telescope

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

The Brightest Young Star Clusters in NGC 5253

The IACOB spectroscopic database: recent updates and first data release

The VLT-FLAMES Tarantula Survey XIX. B-type Supergiants - Atmospheric parameters and nitrogen abundances to investigate the role of binarity and the width of the main sequence

Time calibration with atmospheric muon tracks in the ANTARES neutrino telescope

A Be-type star with a black hole companion

Constraining the neutrino emission of gravitationally lensed Flat-Spectrum Radio Quasars with ANTARES data

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

Detection of bridge emission above 50 GeV from the Crab pulsar with the MAGIC telescopes

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

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

Legacy ExtraGalactic UV Survey (LEGUS) with The Hubble Space Telescope. I. Survey Description

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

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

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

Masgomas-4: Physical characterization of a double-core obscured cluster with a massive and very young stellar population

MGB and the new Galactic O-Star Spectroscopic Survey spectral classification standard grid

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

Multiwavelength observations of Mrk 501 in 2008

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

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

Rotational velocities of single and binary O-type stars in the Tarantula Nebula

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

Searches for clustering in the time integrated skymap of the ANTARES neutrino telescope

The B Fields in OB Stars (BOB) Survey

The IACOB project: II. On the scatter of O-dwarf spectral type -- effective temperature calibrations

The VLT-FLAMES Tarantula Survey XVI. The optical+NIR extinction laws in 30 Doradus and the photometric determination of the effective temperatures of OB stars

The VLT-FLAMES Tarantula Survey XVII. Physical and wind properties of massive stars at the top of the main sequence

Astrophysical parameters of the peculiar X-ray transient IGR J11215-5952

CTA contributions to the 33rd International Cosmic Ray Conference (ICRC2013)

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

Rotational properties of the O-type star population in the Tarantula region

Spectroscopic and physical parameters of Galactic O-type stars. II. Observational constraints on projected rotational and extra broadening velocities as a function of fundamental parameters and stellar evolution

Target Selection for the Apache Point Observatory Galactic Evolution Experiment (APOGEE)

The IACOB project: I. Rotational velocities in Northern Galactic O and early B-type stars revisited. The impact of other sources of line-broadening

The VLT-FLAMES Tarantula Survey - XI. A census of the hot luminous stars and their feedback in 30 Doradus

The VLT-FLAMES Tarantula Survey XII. Rotational velocities of the single O-type stars

The VLT-FLAMES Tarantula Survey. XIII: On the nature of O Vz stars in 30 Doradus

Very Massive Stars (VMS) in the Local Universe

A peculiar Of star in the Local Group galaxy IC 1613

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

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

Discovery of a young and massive stellar cluster: Spectrophotometric near-infrared study of Masgomas-1

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

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

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