Source author record

A. Chiavassa

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

124works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

124 published item(s)

preprint2023arXiv

Primary Cosmic Rays Energy Spectrum and Mean Mass Composition by the Data of the TAIGA Astrophysical Complex

The corrected dependence of the mean depth of the EAS maximum $X_{max}$ on the energy was obtained from the data of the Tunka-133 array for 7 years and the TAIGA-HiSCORE array for 2 years. The parameter $\langle\ln A\rangle$, characterizing the mean mass compositon was derived from these results. The differential energy spectrum of primary cosmic rays in the energy range of $2\cdot 10^{14}$ - $2\cdot 10^{16}$\,eV was reconstructed using the new parameter $Q_{100}$ the Cherenkov light flux at the core distance 100 m.}

preprint2023arXiv

Results from recent analysis of KASCADE-Grande data

KASCADE and its extension array of KASCADE-Grande were devoted to measure individual air showers of cosmic rays in the primary energy range of 100 TeV to 1 EeV. The experiment has substantially contributed to investigate the energy spectrum and mass composition of cosmic rays in the transition region from galactic to extragalactic origin of cosmic rays as well as to quantify the characteristics of hadronic interaction models in the air shower development through validity tests using the multi-detector information from KASCADE-Grande. Although the data accumulation was completed in 2013, data analysis is still continuing. Recently, we investigated the reliability of the new hadronic interactions models of the SIBYLL version 2.3d only with the energy spectra from the KASCADE-Grande data. The evolution of the muon content of high energy air showers in the atmosphere is studied as well, using EPOS-LHC, SIBYLL 2.3, QGSJET-II-04 and SIBYLL 2.3c. In this talk, recent results from KASCADE-Grande and the update of the KASCADE Cosmic Ray Data Centre (KCDC) will be discussed.

preprint2022arXiv

A dusty veil shading Betelgeuse during its Great Dimming

Red supergiants are the most common final evolutionary stage of stars that have initial masses between 8 and 35 times that of the Sun. During this stage, which lasts roughly 100,000 years1, red supergiants experience substantial mass loss. However, the mechanism for this mass loss is unknown. Mass loss may affect the evolutionary path, collapse and future supernova light curve of a red supergiant, and its ultimate fate as either a neutron star or a black hole. From November 2019 to March 2020, Betelgeuse - the second-closest red supergiant to Earth (roughly 220 parsecs, or 724 light years, away) - experienced a historic dimming of its visible brightness. Usually having an apparent magnitude between 0.1 and 1.0, its visual brightness decreased to 1.614 +/- 0.008 magnitudes around 7-13 February 2020 - an event referred to as Betelgeuse's Great Dimming. Here we report high-angular-resolution observations showing that the southern hemisphere of Betelgeuse was ten times darker than usual in the visible spectrum during its Great Dimming. Observations and modelling support a scenario in which a dust clump formed recently in the vicinity of the star, owing to a local temperature decrease in a cool patch that appeared on the photosphere. The directly imaged brightness variations of Betelgeuse evolved on a timescale of weeks. Our findings suggest that a component of mass loss from red supergiants is inhomogeneous, linked to a very contrasted and rapidly changing photosphere

preprint2022arXiv

LOPES 3D -- studies on the benefits of EAS-radio measurements with vertically aligned antennas

The LOPES experiment was a radio interferometer built at the existing air shower array KASCADE-Grande in Karlsruhe, Germany. The last configuration of LOPES was called LOPES 3D and consisted of ten tripole antennas. Each of these antennas consisted of three crossed dipoles east-west, north-south, and vertically aligned. With this, LOPES 3D had the unique possibility to study the benefits of measurements with vertically aligned antennas in the environment of the well understood and calibrated particle detector array KASCADE-Grande. The measurements with three spatially coincident antennas allows a redundant reconstruction of the electric field vector. Several methods to exploit the redundancy were developed and tested. Furthermore, for the first time in LOPES, the background noise could be studied polarization- and direction dependent. With LOPES 3D it could be demonstrated that radio detection reaches a higher efficiency for inclined showers when including measurements with vertically aligned antennas and that the vertical component gets more important for the measurement of inclined showers. In this contribution we discuss a weighting scheme for the best combination of three redundant reconstructed electric field vectors. Furthermore, we discuss the influence of these weighting schemes on the ability to reconstruct air showers using the radio method. We show an estimate of the radio efficiency for inclined showers with focus on the benefits of measurements with vertically aligned antennas and we present the direction dependent noise in the different polarizations.

preprint2022arXiv

Probing Red Supergiant dynamics through photo-center displacements measured by Gaia

Red supergiant (RSGs) are cool massive stars in a late phase of their evolution when the stellar envelope becomes fully convective. They are the brightest stars in the universe at infrared light and can be detected in galaxies far beyond the Local Group, allowing for accurate determination of chemical composition of galaxies. The study of their physical properties is extremely important for various phenomena including the final fate of massive stars as type II supernovae and gravitational wave progenitors. We explore the well-studied nearby young stellar cluster chi Per. Using Gaia EDR3 data, we find the distance of the cluster (d = 2.260+-0.020 kpc). We then investigate the variability of the convection-related surface structure as a source for parallax measurement uncertainty. We use state-of-the-art 3D radiative hydrodynamics simulations with CO5BOLD and the post-processing radiative transfer code OPTIM3D to compute intensity maps in the Gaia G photometric system. We calculate the variabiltiy, as a function of time, of the intensity-weighted mean from the synthetic maps. We then select the RSG stars in the cluster and compare their uncertainty on parallaxes to the predictions of photocentre displacements. The synthetic maps of RSG show extremely irregular and temporal variable surfaces due to convection-related dynamics. Consequentially, the position of the photo-center varies during Gaia measurements between 0.033 and 0.130 AU (up to 5% of the corresponding simulation stellar radius). We argue that the variability of the convection-related surface structures accounts for a substantial part of the Gaia EDR3 parallax error of the RSG sample. We suggest that the variation of the uncertainty on Gaia parallax could be exploited quantitatively using appropriate 3D simulations to extract, in a unique way, important information about the stellar dynamics and parameters of RSG stars.

preprint2022arXiv

The Study of Atmosphere of Hot Jupiters and Their Host Stars

What makes the study of exoplanetary atmospheres so hard is the extraction of its tiny signal from observations, usually dominated by telluric absorption, stellar spectrum and instrumental noise. The High Resolution Spectroscopy has emerged as one of the leading techniques for detecting atomic and molecular species (Birkby 2018), but although it is particularly robust against contaminant absorption in the Earth's atmosphere, the non-stationary stellar spectrum -- in the form of either Doppler shift or distortion of the line profile during planetary transits -- creates a non-negligible source of noise that can alter or even prevent the detection. Recently, significant improvements have been achieved by using 3D, radiative hydrodynamical (RHD) simulations for the star and Global Circulation Models (GCM) for the planet (e.g., Chiavassa & Brogi 2019, Flowers et al. 2019). However, these numerical simulations have been computed independently so far, while acquired spectra are the result of the natural coupling at each phase along the planet orbit. With our work, we aim at generating emission spectra of G,F, and K-type stars and Hot Jupiters and coupling them at any phase of the orbit. This approach is expected to be particularly advantageous for those molecules that are present in both the atmospheres (e.g., CO) and form in the same region of the spectrum, resulting in mixed and overlapped spectral lines. We also present the analysis of transmission spectra of the Hot Saturn WASP-20b, observed in the K-band of the recently upgraded spectrograph CRIRES+ at a resolution ~92, 000 during the first night of the Science Verification of the instrument and that led to a tentative detection of H2O.

preprint2022arXiv

The Tunka-Grande scintillation array: current results

The Tunka-Grande experiment is a scintillation array with about 0.5 sq.km sensitive area at Tunka Valley, Siberia, for measuring charged particles and muons in extensive air showers (EASs). Tunka-Grande is optimized for cosmic ray studies in the energy range 10 PeV to about 1 EeV, where exploring the composition is of fundamental importance for understanding the transition from galactic to extragalactic origin of cosmic rays. This paper attempts to provide a synopsis of the current results of the experiment. In particular, the reconstruction of the all-particle energy spectrum in the range of 10 PeV to 1 EeV based on experimental data from four observation seasons is presented.

preprint2021arXiv

Final results of the LOPES radio interferometer for cosmic-ray air showers

LOPES, the LOFAR prototype station, was an antenna array for cosmic-ray air showers operating from 2003 - 2013 within the KASCADE-Grande experiment. Meanwhile, the analysis is finished and the data of air-shower events measured by LOPES are available with open access in the KASCADE Cosmic Ray Data Center (KCDC). This article intends to provide a summary of the achievements, results, and lessons learned from LOPES. By digital, interferometric beamforming the detection of air showers became possible in the radio-loud environment of the Karlsruhe Institute of Technology (KIT). As a prototype experiment, LOPES tested several antenna types, array configurations and calibration techniques, and pioneered analysis methods for the reconstruction of the most important shower parameters, i.e., the arrival direction, the energy, and mass-dependent observables such as the position of the shower maximum. In addition to a review and update of previously published results, we also present new results based on end-to-end simulations including all known instrumental properties. For this, we applied the detector response to radio signals simulated with the CoREAS extension of CORSIKA, and analyzed them in the same way as measured data. Thus, we were able to study the detector performance more accurately than before, including some previously inaccessible features such as the impact of noise on the interferometric cross-correlation beam. These results led to several improvements, which are documented in this paper and can provide useful input for the design of future cosmic-ray experiments based on the digital radio-detection technique.

preprint2021arXiv

Search for Astrophysical Nanosecond Optical Transients with TAIGA-HiSCORE Array

A wide-angle Cerenkov array TAIGA-HiSCORE (FOV $\sim$0.6 sr), was originally created as a part of TAIGA installation for high-energy gamma-ray astronomy and cosmic ray physics. Array now consist on nearly 100 optical stations on the area of 1 km$^2$. Due to high accuracy and stability ($\sim$1 ns) of time synchronization of the optical stations the accuracy of EAS arrival direction reconstruction is reached 0.1$^\mathrm{o}$. It was proven that the array can also be used to search for nanosecond events of the optical range. The report discusses the method of searching for optical transients using the HiSCORE array and demonstrates its performance on a real example of detecting signals from an artificial Earth satellite. The search for this short flares in the HiSCORE data of the winter season 2018--2019 is carried out. One candidate for double repeater has been detected, but the estimated probability of random simulation of such a transient by background EAS events is not less than 10%, which does not allow us to say that the detected candidate corresponds to a real astrophysical transient. An upper bound on the frequency of optical spikes with flux density of more than $10^{-4} \mathrm{erg/s/cm}^2$ and a duration of more than 5\,ns is established as $\sim 2 \times 10^{-3}$ events/sr/hour.

preprint2021arXiv

Sensitivity of the Cherenkov Telescope Array for probing cosmology and fundamental physics with gamma-ray propagation

The Cherenkov Telescope Array (CTA), the new-generation ground-based observatory for $γ$-ray astronomy, provides unique capabilities to address significant open questions in astrophysics, cosmology, and fundamental physics. We study some of the salient areas of $γ$-ray cosmology that can be explored as part of the Key Science Projects of CTA, through simulated observations of active galactic nuclei (AGN) and of their relativistic jets. Observations of AGN with CTA will enable a measurement of $γ$-ray absorption on the extragalactic background light with a statistical uncertainty below 15% up to a redshift $z=2$ and to constrain or detect $γ$-ray halos up to intergalactic-magnetic-field strengths of at least 0.3pG. Extragalactic observations with CTA also show promising potential to probe physics beyond the Standard Model. The best limits on Lorentz invariance violation from $γ$-ray astronomy will be improved by a factor of at least two to three. CTA will also probe the parameter space in which axion-like particles could constitute a significant fraction, if not all, of dark matter. We conclude on the synergies between CTA and other upcoming facilities that will foster the growth of $γ$-ray cosmology.

preprint2021arXiv

Sensitivity of the Cherenkov Telescope Array to a dark matter signal from the Galactic centre

We provide an updated assessment of the power of the Cherenkov Telescope Array (CTA) to search for thermally produced dark matter at the TeV scale, via the associated gamma-ray signal from pair-annihilating dark matter particles in the region around the Galactic centre. We find that CTA will open a new window of discovery potential, significantly extending the range of robustly testable models given a standard cuspy profile of the dark matter density distribution. Importantly, even for a cored profile, the projected sensitivity of CTA will be sufficient to probe various well-motivated models of thermally produced dark matter at the TeV scale. This is due to CTA's unprecedented sensitivity, angular and energy resolutions, and the planned observational strategy. The survey of the inner Galaxy will cover a much larger region than corresponding previous observational campaigns with imaging atmospheric Cherenkov telescopes. CTA will map with unprecedented precision the large-scale diffuse emission in high-energy gamma rays, constituting a background for dark matter searches for which we adopt state-of-the-art models based on current data. Throughout our analysis, we use up-to-date event reconstruction Monte Carlo tools developed by the CTA consortium, and pay special attention to quantifying the level of instrumental systematic uncertainties, as well as background template systematic errors, required to probe thermally produced dark matter at these energies. "Full likelihood tables complementing our analysis are provided here [ https://doi.org/10.5281/zenodo.4057987 ]"

preprint2021arXiv

The extended atmosphere and circumstellar environment of the cool evolved star VX Sagittarii as seen by MATISSE

Context. VX Sgr is a cool, evolved, and luminous red star whose stellar parameters are difficult to determine, which affects its classification. Aims. We aim to spatially resolve the photospheric extent as well as the circumstellar environment. Methods. We used interferometric observations obtained with the MATISSE instrument in the L (3 to 4 μm), M (4.5 to 5 μm), and N (8 to 13 μm) bands. We reconstructed monochromatic images using the MIRA software. We used 3D radiation-hydrodynamics (RHD) simulations carried out with CO5BOLD and a uniform disc model to estimate the apparent diameter and interpret the stellar surface structures. Moreover, we employed the radiative transfer codes Optim3D and Radmc3D to compute the spectral energy distribution for the L, M, and N bands, respectively. Results. MATISSE observations unveil, for the first time, the morphology of VX Sgr across the L, M, and N bands. The reconstructed images show a complex morphology with brighter areas whose characteristics depend on the wavelength probed. We measured the angular diameter as a function of the wavelength and showed that the photospheric extent in the L and M bands depends on the opacity through the atmosphere. In addition to this, we also concluded that the observed photospheric inhomogeneities can be interpreted as convection-related surface structures. The comparison in the N band yielded a qualitative agreement between the N band spectrum and simple dust radiative transfer simulations. However, it is not possible to firmly conclude on the interpretation of the current data because of the difficulty in constraing the model parameters using the limited accuracy of our absolute flux calibration. Conclusions. MATISSE observations and the derived reconstructed images unveil the appearance of the stellar surface and circumstellar environment across a very large spectral domain for the first time.

preprint2021arXiv

The NEVOD-EAS air-shower array

The Experimental complex NEVOD includes several different setups for studying various components of extensive air showers (EAS) in the energy range from 10^10 to 10^18 eV. The NEVOD-EAS array for detection of the EAS electron-photon component began its data taking in 2018. It is a distributed system of scintillation detectors installed over an area of about 10^4 m^2. A distinctive feature of this array is its cluster organization with different-altitude layout of the detecting elements. The main goal of the NEVOD-EAS array is to obtain an estimation of the primary particle energy for events measured by various detectors of the Experimental complex NEVOD. This paper describes the design, operation principles and data processing of the NEVOD-EAS array. The criteria for the event selection and the accuracy of the EAS parameters reconstruction obtained on the simulated events are discussed. The results of the preliminary analysis of experimental data obtained during a half-year operation are presented.

preprint2020arXiv

Optical interferometry and Gaia measurement uncertainties reveal the physics of asymptotic giant branch stars

Context. Asymptotic giant branch stars are cool luminous evolved stars that are well observable across the Galaxy and populating Gaia data. They have complex stellar surface dynamics Aims. On the AGB star CL Lac, it has been shown that the convection-related variability accounts for a substantial part of the Gaia DR2 parallax error. We observed this star with the MIRC-X beam combiner installed at the CHARA interferometer to detect the presence of stellar surface inhomogeneities. Methods. We performed the reconstruction of aperture synthesis images from the interferometric observations at different wavelengths. Then, we used 3D radiative hydrodynamics simulations of stellar convection with CO5BOLD and the post-processing radiative transfer code Optim3D to compute intensity maps in the spectral channels of MIRC-X observations. Then, we determined the stellar radius and compared the 3D synthetic maps to the reconstructed ones focusing on matching the intensity contrast, the morphology of stellar surface structures, and the photocentre position at two different spectral channels, 1.52 and 1.70 micron, simultaneously. Results. We measured the apparent diameter of CL Lac at two wavelengths and recovered the radius using a Gaia parallax. In addition to this, the reconstructed images are characterised by the presence of a brighter area that largely affects the position of the photocentre. The comparison with 3D simulation shows good agreement with the observations both in terms of contrast and surface structure morphology, meaning that our model is adequate for explaining the observed inhomogenities. Conclusions. This work confirms the presence of convection-related surface structures on an AGB star of Gaia DR2. Our result will help us to take a step forward in exploiting Gaia measurement uncertainties to extract the fundamental properties of AGB stars using appropriate RHD simulations.

preprint2020arXiv

Pulsating chromosphere of classical Cepheids. Calcium infrared triplet and H$α$ profile variations

It has been shown recently that the infrared emission of Cepheids, which is constant over the pulsation cycle, might be due to a pulsating shell of ionized gas of about 15\% of the stellar radius, which could be attributed to the chromospheric activity of Cepheids. The aim of this paper is to investigate the dynamical structure of the chromosphere of Cepheids along the pulsation cycle and quantify its size. We present H$α$ and Calcium Near InfraRed triplet (Ca IR) profile variations using high-resolution spectroscopy with the UVES spectrograph of a sample of 24 Cepheids with a good period coverage from $\approx$ 3 to 60 days. After a qualitative analysis of the spectral lines profiles, we quantify the Van Hoof effect (velocity gradient between the H$α$ and Ca IR) as a function of the period of the Cepheids. Then, we use the Schwarzschild mechanism (a line doubling due to a shock wave) to quantify the size of the chromosphere. We find a significant Van Hoof effect for Cepheids with period larger than $P=10$ days, in particular H$α$ lines are delayed with a velocity gradient up to $Δv \approx$30 km/s compared to Ca IR. We find that the size of the chromosphere of long-period Cepheids is of at least $\approx$ 50\% of the stellar radius, which is consistent at first order with the size of the shell made of ionized gas previously found from the analysis of infrared excess. Last, for most of the long-period Cepheids in the sample, we report a motionless absorption feature in the H$α$ line that we attribute to a circumstellar envelope that surrounds the chromosphere. Analyzing the Ca~IR lines of Cepheids is of importance to potentially unbias the period-luminosity relation from their infrared excess, particularly in the context of forthcoming observations from the Radial Velocity Spectrometer (RVS) on board \textit{Gaia}, that could be sensitive to their chromosphere.

preprint2020arXiv

SPHERE+: Imaging young Jupiters down to the snowline

SPHERE (Beuzit et al,. 2019) has now been in operation at the VLT for more than 5 years, demonstrating a high level of performance. SPHERE has produced outstanding results using a variety of operating modes, primarily in the field of direct imaging of exoplanetary systems, focusing on exoplanets as point sources and circumstellar disks as extended objects. The achievements obtained thus far with SPHERE (~200 refereed publications) in different areas (exoplanets, disks, solar system, stellar physics...) have motivated a large consortium to propose an even more ambitious set of science cases, and its corresponding technical implementation in the form of an upgrade. The SPHERE+ project capitalizes on the expertise and lessons learned from SPHERE to push high contrast imaging performance to its limits on the VLT 8m-telescope. The scientific program of SPHERE+ described in this document will open a new and compelling scientific window for the upcoming decade in strong synergy with ground-based facilities (VLT/I, ELT, ALMA, and SKA) and space missions (Gaia, JWST, PLATO and WFIRST). While SPHERE has sampled the outer parts of planetary systems beyond a few tens of AU, SPHERE+ will dig into the inner regions around stars to reveal and characterize by mean of spectroscopy the giant planet population down to the snow line. Building on SPHERE's scientific heritage and resounding success, SPHERE+ will be a dedicated survey instrument which will strengthen the leadership of ESO and the European community in the very competitive field of direct imaging of exoplanetary systems. With enhanced capabilities, it will enable an even broader diversity of science cases including the study of the solar system, the birth and death of stars and the exploration of the inner regions of active galactic nuclei.

preprint2020arXiv

VLTI-PIONIER imaging of the red supergiant V602 Carinae

Context. Red supergiant stars possess surface features and extended molecular atmospheres. Photospheric convection may be a crucial factor of the levitation of the outer atmospheric layers. However, the mechanism responsible is still poorly understood. Aims. We image the stellar surface of V602 Carinae (V602 Car) to constrain the morphology and contrast of the surface features and of the extended atmospheric layers. Methods. We observed V602 Car with the Very Large Telescope Interferometer (VLTI) PIONIER instrument (1.53-1.78 $\mathrmμ$m) between May and July 2016, and April and July 2019 with different telescope configurations. We compared the image reconstructions with 81 temporal snapshots of 3D radiative-hydrodynamics (RHD) CO$^5$BOLD simulations in terms of contrast and morphology, using the Structural Similarity Index. Results. The interferometric data are compatible with an overall spherical disk of angular diameter 4.4$\pm$0.2 mas, and an extended molecular layer. In 2016, the reconstructed image reveals a bright arc-like feature toward the northern rim of the photospheric surface. In 2019, an arc-like feature is seen at a different orientation and a new peak of emission is detected on the opposite side. The contrasts of the reconstructed surface images are 11%$\pm$2% and 9%$\pm$2% for 2016 and 2019, respectively. The morphology and contrast of the two images are consistent with 3D RHD simulations, within our achieved spatial resolution and dynamic range. The extended molecular layer contributes 10--13% of the total flux with an angular diameter of 6--8 mas. It is present but not clearly visible in the reconstructed images because it is close to the limits of the achieved dynamic range. The presence of the molecular layer is not reproduced by the 3D RHD simulations. Conclusions...

preprint2019arXiv

Surface imaging of cool evolved stars in the era of the ELT

Cool evolved stars are the main source of chemical enrichment of the interstellar medium. Understanding their mass loss offers a unique opportunity to study the cycle of matter. We discuss interferometric studies and their comparison to latest state-of-the-art dynamic model atmospheres. They show broad agreement for asymptotic giant branch stars. For red supergiants, however, current models cannot explain observed extensions by far, pointing to missing physical processes in their models, and uncertainties in our general understanding of mass loss. We present ongoing imaging and time-series observations that may provide the strongest constraint and may help to identify missing dynamic processes. VLTI studies will remain the highest spatial resolution observations at ESO into the ELT era, complemented by ALMA observations. We discuss crucial improvements in both instrumental and operational areas for surface imaging of cool evolved stars in the era of the ELT.

preprint2016arXiv

Azimuthal asymmetry in the risetime of the surface detector signals of the Pierre Auger Observatory

The azimuthal asymmetry in the risetime of signals in Auger surface detector stations is a source of information on shower development. The azimuthal asymmetry is due to a combination of the longitudinal evolution of the shower and geometrical effects related to the angles of incidence of the particles into the detectors. The magnitude of the effect depends upon the zenith angle and state of development of the shower and thus provides a novel observable, $(\sec θ)_\mathrm{max}$, sensitive to the mass composition of cosmic rays above $3 \times 10^{18}$ eV. By comparing measurements with predictions from shower simulations, we find for both of our adopted models of hadronic physics (QGSJETII-04 and EPOS-LHC) an indication that the mean cosmic-ray mass increases slowly with energy, as has been inferred from other studies. However, the mass estimates are dependent on the shower model and on the range of distance from the shower core selected. Thus the method has uncovered further deficiencies in our understanding of shower modelling that must be resolved before the mass composition can be inferred from $(\sec θ)_\mathrm{max}$.

preprint2016arXiv

Detection of extensive air showers with the NEVOD-EAS cluster type detector

A new cluster type shower array NEVOD-EAS is designed for estimating the size, axis position and arrival direction of extensive air showers registered by the Cherenkov water detector NEVOD and coordinate detector DECOR (Moscow, Russia). In 2015-2016, the central part of the array was deployed and started its operation. It includes 4 independent clusters of scintillation detector stations located around the NEVOD-DECOR experimental complex on the area of about $10^{4} m^{2}$. This article presents the results of studying amplitude and timing characteristics of the array clusters which are critical for EAS parameters reconstruction, as well as the examples of registered events.

preprint2016arXiv

Discovery of a complex linearly polarized spectrum of Betelgeuse dominated by depolarization of the continuum

Betelgeuse is an M supergiant that harbors spots and giant granules at its surface and presents linear polarization of its continuum. We have previously discovered linear polarization signatures associated with individual lines in the spectra of cool and evolved stars. Here, we investigate whether a similar linearly polarized spectrum exists for Betelgeuse. We used the spectropolarimeter Narval, combining multiple polarimetric sequences to obtain high signal-to-noise ratio spectra of individual lines, as well as the least-squares deconvolution (LSD) approach. We have discovered the existence of a linearly polarized spectrum for Betelgeuse, detecting a rather strong signal (at a few times 10$^{-4}$ of the continuum intensity level), both in individual lines and in the LSD profiles. Studying its properties and the signal observed for the resonant \ion{Na}{i}\,D lines, we conclude that we are mainly observing depolarization of the continuum by the absorption lines. The linear polarization of the Betelgeuse continuum is due to the anisotropy of the radiation field induced by brightness spots at the surface and Rayleigh scattering in the atmosphere. We have developed a geometrical model to interpret the observed polarization, from which we infer the presence of two brightness spots and their positions on the surface of Betelgeuse. We show that applying the model to each velocity bin along the Stokes Q and U profiles allows the derivation of a map of the bright spots. We use the Narval linear polarization observations of Betelgeuse obtained over a period of 1.4 years to study the evolution of the spots. Our study of the linearly polarized spectrum of Betelgeuse provides a novel method for studying the evolution of brightness spots at its surface and complements quasi-simultaneous observations obtained with PIONIER at the VLTI.

preprint2016arXiv

EELT-HIRES the high-resolution spectrograph for the E-ELT

The first generation of E-ELT instruments will include an optical-infrared High Resolution Spectrograph, conventionally indicated as EELT-HIRES, which will be capable of providing unique breakthroughs in the fields of exoplanets, star and planet formation, physics and evolution of stars and galaxies, cosmology and fundamental physics. A 2-year long phase A study for EELT-HIRES has just started and will be performed by a consortium composed of institutes and organisations from Brazil, Chile, Denmark, France, Germany, Italy, Poland, Portugal, Spain, Sweden, Switzerland and United Kingdom. In this paper we describe the science goals and the preliminary technical concept for EELT-HIRES which will be developed during the phase A, as well as its planned development and consortium organisation during the study.

preprint2016arXiv

Evidence for a mixed mass composition at the `ankle' in the cosmic-ray spectrum

We report a first measurement for ultra-high energy cosmic rays of the correlation between the depth of shower maximum and the signal in the water Cherenkov stations of air-showers registered simultaneously by the fluorescence and the surface detectors of the Pierre Auger Observatory. Such a correlation measurement is a unique feature of a hybrid air-shower observatory with sensitivity to both the electromagnetic and muonic components. It allows an accurate determination of the spread of primary masses in the cosmic-ray flux. Up till now, constraints on the spread of primary masses have been dominated by systematic uncertainties. The present correlation measurement is not affected by systematics in the measurement of the depth of shower maximum or the signal in the water Cherenkov stations. The analysis relies on general characteristics of air showers and is thus robust also with respect to uncertainties in hadronic event generators. The observed correlation in the energy range around the `ankle' at $\lg(E/{\rm eV})=18.5-19.0$ differs significantly from expectations for pure primary cosmic-ray compositions. A light composition made up of proton and helium only is equally inconsistent with observations. The data are explained well by a mixed composition including nuclei with mass $A > 4$. Scenarios such as the proton dip model, with almost pure compositions, are thus disfavoured as the sole explanation of the ultrahigh-energy cosmic-ray flux at Earth.

preprint2016arXiv

Nanosecond-level time synchronization of autonomous radio detector stations for extensive air showers

To exploit the full potential of radio measurements of cosmic-ray air showers at MHz frequencies, a detector timing synchronization within 1 ns is needed. Large distributed radio detector arrays such as the Auger Engineering Radio Array (AERA) rely on timing via the Global Positioning System (GPS) for the synchronization of individual detector station clocks. Unfortunately, GPS timing is expected to have an accuracy no better than about 5 ns. In practice, in particular in AERA, the GPS clocks exhibit drifts on the order of tens of ns. We developed a technique to correct for the GPS drifts, and an independent method is used for cross-checks that indeed we reach nanosecond-scale timing accuracy by this correction. First, we operate a "beacon transmitter" which emits defined sine waves detected by AERA antennas recorded within the physics data. The relative phasing of these sine waves can be used to correct for GPS clock drifts. In addition to this, we observe radio pulses emitted by commercial airplanes, the position of which we determine in real time from Automatic Dependent Surveillance Broadcasts intercepted with a software-defined radio. From the known source location and the measured arrival times of the pulses we determine relative timing offsets between radio detector stations. We demonstrate with a combined analysis that the two methods give a consistent timing calibration with an accuracy of 2 ns or better. Consequently, the beacon method alone can be used in the future to continuously determine and correct for GPS clock drifts in each individual event measured by AERA.

preprint2016arXiv

Near-infrared spectro-interferometry of Mira variables and comparisons to 1D dynamic model atmospheres and 3D convection simulations

We obtained a total of 20 near-infrared K-band spectro-interferometric snapshot observations of the Mira variables o Cet, R Leo, R Aqr, X Hya, W Vel, and R Cnc with a spectral resolution of about 1500. We compared observed flux and visibility spectra with predictions by CODEX 1D dynamic model atmospheres and with azimuthally averaged intensities based on CO5BOLD 3D dynamic model atmospheres including convection. Our visibility data confirm the presence of spatially extended molecular atmospheres located above the continuum radii with large-scale inhomogeneities or clumps that contribute a few percent of the total flux. The detailed structure of the inhomogeneities or clumps show a variability on time scales of 3 months and above. Both modeling attempts provided satisfactory fits to our data. In particular, they are both consistent with the observed decrease in the visibility function at molecular bands of water vapor and CO, indicating a spatially extended molecular atmosphere. Observational variability phases are mostly consistent with those of the best-fit CODEX models, except for near-maximum phases, where data are better described by near-minimum models. Rosseland angular diameters derived from the model fits are broadly consistent between those based on the 1D and the 3D models and with earlier observations. We derived fundamental parameters including absolute radii, effective temperatures, and luminosities for our sources. Our results provide a first observational support for theoretical results that shocks induced by convection and pulsation in the 3D CO5BOLD models of AGB stars are roughly spherically expanding and of similar nature to those of self-excited pulsations in 1D CODEX models. Unlike for red supergiants, the pulsation- and shock-induced dynamics can levitate the molecular atmospheres of Mira variables to extensions that are consistent with observations.

preprint2016arXiv

New detectors of the Experimental complex NEVOD for multicomponent EAS detection

Experimental complex (EC) NEVOD includes a number of unique experimental facilities for studies of main components of cosmic rays on the Earth's surface. The complex is used for the basic research of CR flux characteristics and their interactions in the energy range 10^15 - 10^19 eV, and for applied investigations directed to the development of methods of the muon diagnostics of the atmosphere and the Earth's magnetosphere and near-terrestrial space. To extend the experimental capabilities and raising the status of the installation to the Mega Science level, nowadays new large-scale detectors: array for the EAS registration - NEVOD-EAS, detector of atmospheric neutrons - URAN, and large-area coordinate-tracking detector - TREK, are being deployed around EC NEVOD. The description of new detectors and a common trigger system to ensure the joint operation together with other detectors of EC NEVOD are presented.

preprint2016arXiv

Prototype muon detectors for the AMIGA component of the Pierre Auger Observatory

Auger Muons and Infill for the Ground Array) is an upgrade of the Pierre Auger Observatory to extend its range of detection and to directly measure the muon content of the particle showers. It consists of an infill of surface water-Cherenkov detectors accompanied by buried scintillator detectors used for muon counting. The main objectives of the AMIGA engineering array, referred to as the Unitary Cell, are to identify and resolve all engineering issues as well as to understand the muon-number counting uncertainties related to the design of the detector. The mechanical design, fabrication and deployment processes of the muon counters of the Unitary Cell are described in this document. These muon counters modules comprise sealed PVC casings containing plastic scintillation bars, wavelength-shifter optical fibers, 64 pixel photomultiplier tubes, and acquisition electronics. The modules are buried approximately 2.25 m below ground level in order to minimize contamination from electromagnetic shower particles. The mechanical setup, which allows access to the electronics for maintenance, is also described in addition to tests of the modules' response and integrity. The completed Unitary Cell has measured a number of air showers of which a first analysis of a sample event is included here.

preprint2016arXiv

Science with MATISSE

We present an overview of the scientific potential of MATISSE, the Multi Aperture mid-Infrared SpectroScopic Experiment for the Very Large Telescope Interferometer. For this purpose we outline selected case studies from various areas, such as star and planet formation, active galactic nuclei, evolved stars, extrasolar planets, and solar system minor bodies and discuss strategies for the planning and analysis of future MATISSE observations. Moreover, the importance of MATISSE observations in the context of complementary high-angular resolution observations at near-infrared and submillimeter/millimeter wavelengths is highlighted.

preprint2016arXiv

Search for correlations between the arrival directions of IceCube neutrino events and ultrahigh-energy cosmic rays detected by the Pierre Auger Observatory and the Telescope Array

This paper presents the results of different searches for correlations between very high-energy neutrino candidates detected by IceCube and the highest-energy cosmic rays measured by the Pierre Auger Observatory and the Telescope Array. We first consider samples of cascade neutrino events and of high-energy neutrino-induced muon tracks, which provided evidence for a neutrino flux of astrophysical origin, and study their cross-correlation with the ultrahigh-energy cosmic ray (UHECR) samples as a function of angular separation. We also study their possible directional correlations using a likelihood method stacking the neutrino arrival directions and adopting different assumptions on the size of the UHECR magnetic deflections. Finally, we perform another likelihood analysis stacking the UHECR directions and using a sample of through-going muon tracks optimized for neutrino point-source searches with sub-degree angular resolution. No indications of correlations at discovery level are obtained for any of the searches performed. The smallest of the p-values comes from the search for correlation between UHECRs with IceCube high-energy cascades, a result that should continue to be monitored.

preprint2016arXiv

Testing Hadronic Interactions at Ultrahigh Energies with Air Showers Measured by the Pierre Auger Observatory

Ultrahigh energy cosmic ray air showers probe particle physics at energies beyond the reach of accelerators. Here we introduce a new method to test hadronic interaction models without relying on the absolute energy calibration, and apply it to events with primary energy 6-16 EeV (E_CM = 110-170 TeV), whose longitudinal development and lateral distribution were simultaneously measured by the Pierre Auger Observatory. The average hadronic shower is 1.33 +- 0.16 (1.61 +- 0.21) times larger than predicted using the leading LHC-tuned models EPOS-LHC (QGSJetII-04), with a corresponding excess of muons.

preprint2016arXiv

The close circumstellar environment of Betelgeuse. IV. VLTI/PIONIER interferometric monitoring of the photosphere

Context. The mass-loss mechanism of cool massive evolved stars is poorly understood. The proximity of Betelgeuse makes it an appealing target to study its atmosphere, map the shape of its envelope, and follow the structure of its wind from the photosphere out to the interstellar medium. Aims. A link is suspected between the powerful convective motions in Betelgeuse and its mass loss. We aim to constrain the spatial structure and temporal evolution of the convective pattern on the photosphere and to search for evidence of this link. Methods. We report new interferometric observations in the infrared H band using the VLTI/PIONIER instrument. We monitored the photosphere of Betelgeuse between 2012 January and 2014 November to look for evolutions that may trigger the outflow. Results. Our interferometric observations at low spatial frequencies are compatible with the presence of a hot spot on the photosphere that has a characteristic width of one stellar radius. It appears to be superposed on the smaller scale convective pattern. In the higher spatial frequency domain, we observe a significant difference between the observations and the predictions of 3D hydrodynamical simulations. Conclusions. We bring new evidence for the presence of a convective pattern in the photosphere of red supergiants. The inferred hot spot is probably the top of a giant convection cell although an asymmetric extension of the star cannot be excluded by these interferometric observations alone. The properties of the observed surface features show a stronger contrast and inhomogeneity as predicted by 3D radiative hydrodynamical simulations. We propose that the large observed feature is modifying the signature of the convective pattern at the surface of the star in a way that simulations cannot reproduce.

preprint2016arXiv

The Pierre Auger Observatory Upgrade - Preliminary Design Report

The Pierre Auger Observatory has begun a major Upgrade of its already impressive capabilities, with an emphasis on improved mass composition determination using the surface detectors of the Observatory. Known as AugerPrime, the upgrade will include new 4 m$^2$ plastic scintillator detectors on top of all 1660 water-Cherenkov detectors, updated and more flexible surface detector electronics, a large array of buried muon detectors, and an extended duty cycle for operations of the fluorescence detectors. This Preliminary Design Report was produced by the Collaboration in April 2015 as an internal document and information for funding agencies. It outlines the scientific and technical case for AugerPrime. We now release it to the public via the arXiv server. We invite you to review the large number of fundamental results already achieved by the Observatory and our plans for the future.

preprint2016arXiv

VLTI/AMBER spectro-interferometry of the late-type supergiants V766 Cen (=HR 5171 A), sigma Oph, BM Sco, and HD 206859

We add four warmer late-type supergiants to our previous spectro-interferometric studies of red giants and supergiants. V766 Cen (=HR 5171 A) is found to be a high-luminosity log(L/L_sun)=5.8+-0.4 source of Teff 4290+-760 K and radius 1490+-540 Rsun located close to both the Hayashi and Eddington limits; this source is consistent with a 40 Msun evolutionary track without rotation and current mass 27-36 Msun. It exhibits NaI in emission arising from a shell of radius 1.5 Rphot and a photocenter displacement of about 0.1 Rphot. V766 Cen shows strong extended molecular (CO) layers and a dusty circumstellar background component. This suggest an optically thick pseudo-photosphere at about 1.5 Rphot at the onset of the wind. V766 Cen is a red supergiant located close to the Hayashi limit instead of a yellow hypergiant already evolving back toward warmer Teff as previously discussed. The stars sigma Oph, BM Sco, and HD 206859 are found to have lower luminosities of about log(L/Lsun)=3.4-3.5 and Teff of 3900-5300 K, corresponding to 5-9 Msun tracks. They do not show extended molecular layers as observed for higher luminosity red supergiants of our sample. BM Sco shows an unusually strong contribution by an over-resolved circumstellar dust component. These stars are more likely high-mass red giants instead of red supergiants. This leaves us with an unsampled locus in the HR diagram corresponding to luminosities log(L/Lsun)~3.8-4.8 or masses 10-13 Msun, possibly corresponding to the mass region where stars explode as type II-P supernovae during the RSG stage. Our previously found relation of increasing strength of extended molecular layers with increasing luminosities is now confirmed to extend to double our previous luminosities and up to the Eddington limit. This might further point to steadily increasing radiative winds with increasing luminosity. [Abridged]

preprint2015arXiv

A new view on exoplanet transits: Transit of Venus described using three-dimensional solar atmosphere Stagger-grid simulations

Stellar activity and, in particular, convection-related surface structures, potentially cause fluctuations that can affect the transit light curves. Surface convection simulations can help the interpretation of ToV. We used realistic three-dimensional radiative hydrodynamical simulation of the Sun from the Stagger-grid and synthetic images computed with the radiative transfer code Optim3D to provide predictions for the transit of Venus in 2004 observed by the satellite ACRIMSAT. We computed intensity maps from RHD simulation of the Sun and produced synthetic stellar disk image. We computed the light curve and compared it to the ACRIMSAT observations and also to the light curves obtained with solar surface representations carried out using radial profiles with different limb-darkening laws. We also applied the same spherical tile imaging method to the observations of center-to-limb Sun granulation with HINODE. We managed to explain ACRIMSAT observations of 2004 ToV and showed that the granulation pattern causes fluctuations in the transit light curve. We evaluated the contribution of the granulation to the ToV. We showed that the granulation pattern can partially explain the observed discrepancies between models and data. This confirms that the limb-darkening and the granulation pattern simulated in 3D RHD Sun represent well what is imaged by HINODE. In the end, we found that the Venus's aureole contribution during ToV is less intense than the solar photosphere, and thus negligible. Being able to explain consistently the data of 2004 ToV is a new step forward for 3D RHD simulations that are becoming essential for the detection and characterization of exoplanets. They show that the granulation have to be considered as an intrinsic incertitude, due to the stellar variability, on precise measurements of exoplanet transits of, most likely, planets with small diameters.

preprint2015arXiv

Comparing LOPES measurements of air-shower radio emission with REAS 3.11 and CoREAS simulations

Cosmic ray air showers emit radio pulses at MHz frequencies, which can be measured with radio antenna arrays - like LOPES at the Karlsruhe Institute of Technology in Germany. To improve the understanding of the radio emission, we test theoretical descriptions with measured data. The observables used for these tests are the absolute amplitude of the radio signal, and the shape of the radio lateral distribution. We compare lateral distributions of more than 500 LOPES events with two recent and public Monte Carlo simulation codes, REAS 3.11 and CoREAS (v 1.0). The absolute radio amplitudes predicted by REAS 3.11 are in good agreement with the LOPES measurements. The amplitudes predicted by CoREAS are lower by a factor of two, and marginally compatible with the LOPES measurements within the systematic scale uncertainties. In contrast to any previous versions of REAS, REAS 3.11 and CoREAS now reproduce the shape of the measured lateral distributions correctly. This reflects a remarkable progress compared to the situation a few years ago, and it seems that the main processes for the radio emission of air showers are now understood: The emission is mainly due to the geomagnetic deflection of the electrons and positrons in the shower. Less important but not negligible is the Askaryan effect (net charge variation). Moreover, we confirm that the refractive index of the air plays an important role, since it changes the coherence conditions for the emission: Only the new simulations including the refractive index can reproduce rising lateral distributions which we observe in a few LOPES events. Finally, we show that the lateral distribution is sensitive to the energy and the mass of the primary cosmic ray particles.

preprint2015arXiv

Improved absolute calibration of LOPES measurements and its impact on the comparison with REAS 3.11 and CoREAS simulations

LOPES was a digital antenna array detecting the radio emission of cosmic-ray air showers. The calibration of the absolute amplitude scale of the measurements was done using an external, commercial reference source, which emits a frequency comb with defined amplitudes. Recently, we obtained improved reference values by the manufacturer of the reference source, which significantly changed the absolute calibration of LOPES. We reanalyzed previously published LOPES measurements, studying the impact of the changed calibration. The main effect is an overall decrease of the LOPES amplitude scale by a factor of $2.6 \pm 0.2$, affecting all previously published values for measurements of the electric-field strength. This results in a major change in the conclusion of the paper 'Comparing LOPES measurements of air-shower radio emission with REAS 3.11 and CoREAS simulations' published in Astroparticle Physics 50-52 (2013) 76-91: With the revised calibration, LOPES measurements now are compatible with CoREAS simulations, but in tension with REAS 3.11 simulations. Since CoREAS is the latest version of the simulation code incorporating the current state of knowledge on the radio emission of air showers, this new result indicates that the absolute amplitude prediction of current simulations now is in agreement with experimental data.

preprint2015arXiv

Investigation of the radio wavefront of air showers with LOPES measurements and CoREAS simulations (ARENA 2014)

We investigated the radio wavefront of cosmic-ray air showers with LOPES measurements and CoREAS simulations: the wavefront is of approximately hyperbolic shape and its steepness is sensitive to the shower maximum. For this study we used 316 events with an energy above 0.1 EeV and zenith angles below $45^\circ$ measured by the LOPES experiment. LOPES was a digital radio interferometer consisting of up to 30 antennas on an area of approximately 200 m x 200 m at an altitude of 110 m above sea level. Triggered by KASCADE-Grande, LOPES measured the radio emission between 43 and 74 MHz, and our analysis might strictly hold only for such conditions. Moreover, we used CoREAS simulations made for each event, which show much clearer results than the measurements suffering from high background. A detailed description of our result is available in our recent paper published in JCAP09(2014)025. The present proceeding contains a summary and focuses on some additional aspects, e.g., the asymmetry of the wavefront: According to the CoREAS simulations the wavefront is slightly asymmetric, but on a much weaker level than the lateral distribution of the radio amplitude.

preprint2015arXiv

KCDC - The KASCADE Cosmic-ray Data Centre

KCDC, the KASCADE Cosmic-ray Data Centre, is a web portal, where data of astroparticle physics experiments will be made available for the interested public. The KASCADE experiment, financed by public money, was a large-area detector for the measurement of high-energy cosmic rays via the detection of air showers. KASCADE and its extension KASCADE-Grande stopped finally the active data acquisition of all its components including the radio EAS experiment LOPES end of 2012 after more than 20 years of data taking. In a first release, with KCDC we provide to the public the measured and reconstructed parameters of more than 160 million air showers. In addition, KCDC provides the conceptional design, how the data can be treated and processed so that they are also usable outside the community of experts in the research field. Detailed educational examples make a use also possible for high-school students and early stage researchers.

preprint2015arXiv

Pathways for Observing Stellar Surfaces Using 3D Hydrodynamical Simulations of Evolved Stars

Evolved stars are among the largest and brightest stars and they are ideal targets for the new generation of sensitive, high resolution instrumentation that provides spectrophotometric, interferometric, astrometric, and imaging observables. The interpretation of the complex stellar surface images requires numerical simulations of stellar convection that take into account multi-dimensional time-dependent radiation hydrodynamics with realistic input physics. We show how the evolved star simulations are obtained using the radiative hydrodynamics code CO5BOLD and how the accurate observables are computed with the post-processing radiative transfer code Optim3D. The synergy between observations and theoretical work is supported by a proper and quantitative analysis using these simulations, and by strong constraints from the observational side.

preprint2015arXiv

Revised absolute amplitude calibration of the LOPES experiment

One of the main aims of the LOPES experiment was the evaluation of the absolute amplitude of the radio signal of air showers. This is of special interest since the radio technique offers the possibility for an independent and highly precise determination of the energy scale of cosmic rays on the basis of signal predictions from Monte Carlo simulations. For the calibration of the amplitude measured by LOPES we used an external source. Previous comparisons of LOPES measurements and simulations of the radio signal amplitude predicted by CoREAS revealed a discrepancy of the order of a factor of two. A re-measurement of the reference calibration source, now performed for the free field, was recently performed by the manufacturer. The updated calibration values lead to a lowering of the reconstructed electric field measured by LOPES by a factor of $2.6 \pm 0.2$ and therefore to a significantly better agreement with CoREAS simulations. We discuss the updated calibration and its impact on the LOPES analysis results.

preprint2015arXiv

Stellar granulation and interferometry

Stars are not smooth. Their photosphere is covered by a granulation pattern associated with the heat transport by convection. The convection-related surface structures have different size, depth, and temporal variations with respect to the stellar type. The related activity (in addition to other phenomena such as magnetic spots, rotation, dust, etc.) potentially causes bias in stellar parameters determination, radial velocity, chemical abundances determinations, and exoplanet transit detections. The role of long-baseline interferometric observations in this astrophysical context is crucial to characterize the stellar surface dynamics and correct the potential biases. In this Chapter, we present how the granulation pattern is expected for different kind of stellar types ranging from main sequence to extremely evolved stars of different masses and how interferometric techniques help to study their photospheric dynamics.

preprint2015arXiv

The close circumstellar environment of Betelgeuse - III. SPHERE/ZIMPOL visible polarimetry of the inner envelope and photosphere

The physical mechanism through which the outgoing material of massive red supergiants is accelerated above the escape velocity is unclear. Thanks to the transparency of its circumstellar envelope, the nearby red supergiant Betelgeuse gives the opportunity to probe the innermost layers of the envelope of a typical red supergiant down to the photosphere, i.e. where the acceleration of the wind is expected to occur. We took advantage of the SPHERE/ZIMPOL adaptive optics imaging polarimeter to resolve the visible photosphere and close envelope of Betelgeuse. We detect an asymmetric gaseous envelope inside a radius of 2 to 3 times the near-infrared photospheric radius of the star (R*), and a significant Halpha emission mostly contained within 3 R*. From the polarimetric signal, we also identify the signature of dust scattering in an asymmetric and incomplete dust shell located at a similar radius. The presence of dust so close to the star may have a significant impact on the wind acceleration through radiative pressure on the grains. The 3 R* radius emerges as a major interface between the hot gaseous and dusty envelopes. The detected asymmetries strengthen previous indications that the mass loss of Betelgeuse is likely tied to the vigorous convective motions in its atmosphere.

preprint2015arXiv

The IceCube Neutrino Observatory, the Pierre Auger Observatory and the Telescope Array: Joint Contribution to the 34th International Cosmic Ray Conference (ICRC 2015)

We have conducted three searches for correlations between ultra-high energy cosmic rays detected by the Telescope Array and the Pierre Auger Observatory, and high-energy neutrino candidate events from IceCube. Two cross-correlation analyses with UHECRs are done: one with 39 cascades from the IceCube `high-energy starting events' sample and the other with 16 high-energy `track events'. The angular separation between the arrival directions of neutrinos and UHECRs is scanned over. The same events are also used in a separate search using a maximum likelihood approach, after the neutrino arrival directions are stacked. To estimate the significance we assume UHECR magnetic deflections to be inversely proportional to their energy, with values $3^\circ$, $6^\circ$ and $9^\circ$ at 100 EeV to allow for the uncertainties on the magnetic field strength and UHECR charge. A similar analysis is performed on stacked UHECR arrival directions and the IceCube sample of through-going muon track events which were optimized for neutrino point-source searches.

preprint2015arXiv

What causes the large extensions of red-supergiant atmospheres? Comparisons of interferometric observations with 1-D hydrostatic, 3-D convection, and 1-D pulsating model atmospheres

We present the atmospheric structure and the fundamental parameters of three red supergiants, increasing the sample of RSGs observed by near-infrared spectro-interferometry. Additionally, we test possible mechanisms that may explain the large observed atmospheric extensions of RSGs. We carried out spectro-interferometric observations of 3 RSGs in the near-infrared K-band with the VLTI/AMBER instrument at medium spectral resolution. To comprehend the extended atmospheres, we compared our observational results to predictions by available hydrostatic PHOENIX, available 3-D convection, and new 1-D self-excited pulsation models of RSGs. Our near-infrared flux spectra are well reproduced by the PHOENIX model atmospheres. The continuum visibility values are consistent with a limb-darkened disk as predicted by the PHOENIX models, allowing us to determine the angular diameter and the fundamental parameters of our sources. Nonetheless, in the case of V602 Car and HD 95686, the PHOENIX model visibilities do not predict the large observed extensions of molecular layers, most remarkably in the CO bands. Likewise, the 3-D convection models and the 1-D pulsation models with typical parameters of RSGs lead to compact atmospheric structures as well, which are similar to the structure of the hydrostatic PHOENIX models. They can also not explain the observed decreases in the visibilities and thus the large atmospheric molecular extensions. The full sample of our RSGs indicates increasing observed atmospheric extensions with increasing luminosity and decreasing surface gravity, and no correlation with effective temperature or variability amplitude, which supports a scenario of radiative acceleration on Doppler-shifted molecular lines.

preprint2014arXiv

3D hydrodynamical simulations of evolved stars and observations of stellar surfaces

Evolved stars are among the largest and brightest stars and they are ideal targets for the new generation of sensitive, high resolution instrumentation that pro- vides spectrophotometric, interferometric, astrometric, and imaging observables. The interpretation of the complex stellar surface images requires numerical simulations of stellar convection that take into account multi-dimensional time-dependent radiation hydrodynamics with realistic input physics. We show how the evolved star simulations are obtained using the radiative hydrodynamics code CO5BOLD and how the accurate observables are computed with the post-processing radiative transfer code Optim3D. The synergy between observations and theoretical work is supported by a proper and quantitative analysis using these simulations, and by strong constraints from the obser- vational side.

preprint2014arXiv

A search for point sources of EeV photons

Measurements of air showers made using the hybrid technique developed with the fluorescence and surface detectors of the Pierre Auger Observatory allow a sensitive search for point sources of EeV photons anywhere in the exposed sky. A multivariate analysis reduces the background of hadronic cosmic rays. The search is sensitive to a declination band from -85° to +20°, in an energy range from 10^17.3 eV to 10^18.5 eV. No photon point source has been detected. An upper limit on the photon flux has been derived for every direction. The mean value of the energy flux limit that results from this, assuming a photon spectral index of -2, is 0.06 eV cm^-2 s^-1, and no celestial direction exceeds 0.25 eV cm^-2 s^-1. These upper limits constrain scenarios in which EeV cosmic ray protons are emitted by non-transient sources in the Galaxy.

preprint2014arXiv

A Targeted Search for Point Sources of EeV Neutrons

A flux of neutrons from an astrophysical source in the Galaxy can be detected in the Pierre Auger Observatory as an excess of cosmic-ray air showers arriving from the direction of the source. To avoid the statistical penalty for making many trials, classes of objects are tested in combinations as nine "target sets", in addition to the search for a neutron flux from the Galactic Center or from the Galactic Plane. Within a target set, each candidate source is weighted in proportion to its electromagnetic flux, its exposure to the Auger Observatory, and its flux attenuation factor due to neutron decay. These searches do not find evidence for a neutron flux from any class of candidate sources. Tabulated results give the combined p-value for each class, with and without the weights, and also the flux upper limit for the most significant candidate source within each class. These limits on fluxes of neutrons significantly constrain models of EeV proton emission from non-transient discrete sources in the Galaxy.

preprint2014arXiv

First Experimental Characterization of Microwave Emission from Cosmic Ray Air Showers

We report the first direct measurement of the overall characteristics of microwave radio emission from extensive air showers. Using a trigger provided by the KASCADE-Grande air shower array, the signals of the microwave antennas of the CROME (Cosmic-Ray Observation via Microwave Emission) experiment have been read out and searched for signatures of radio emission by high-energy air showers in the GHz frequency range. Microwave signals have been detected for more than 30 showers with energies above 3*10^16 eV. The observations presented in this Letter are consistent with a mainly forward-directed and polarised emission process in the GHz frequency range. The measurements show that microwave radiation offers a new means of studying air showers at energies above 10^17 eV.

preprint2014arXiv

Reconstruction of inclined air showers detected with the Pierre Auger Observatory

We describe the method devised to reconstruct inclined cosmic-ray air showers with zenith angles greater than $60^\circ$ detected with the surface array of the Pierre Auger Observatory. The measured signals at the ground level are fitted to muon density distributions predicted with atmospheric cascade models to obtain the relative shower size as an overall normalization parameter. The method is evaluated using simulated showers to test its performance. The energy of the cosmic rays is calibrated using a sub-sample of events reconstructed with both the fluorescence and surface array techniques. The reconstruction method described here provides the basis of complementary analyses including an independent measurement of the energy spectrum of ultra-high energy cosmic rays using very inclined events collected by the Pierre Auger Observatory.

preprint2014arXiv

Reconstruction of the energy and depth of maximum of cosmic-ray air-showers from LOPES radio measurements

LOPES is a digital radio interferometer located at Karlsruhe Institute of Technology (KIT), Germany, which measures radio emission from extensive air showers at MHz frequencies in coincidence with KASCADE-Grande. In this article, we explore a method (slope method) which leverages the slope of the measured radio lateral distribution to reconstruct crucial attributes of primary cosmic rays. First, we present an investigation of the method on the basis of pure simulations. Second, we directly apply the slope method to LOPES measurements. Applying the slope method to simulations, we obtain uncertainties on the reconstruction of energy and depth of shower maximum Xmax of 13% and 50 g/cm^2, respectively. Applying it to LOPES measurements, we are able to reconstruct energy and Xmax of individual events with upper limits on the precision of 20-25% for the primary energy and 95 g/cm^2 for Xmax, despite strong human-made noise at the LOPES site.

preprint2014arXiv

Searches for Large-Scale Anisotropy in the Arrival Directions of Cosmic Rays Detected above Energy of $10^{19}$ eV at the Pierre Auger Observatory and the Telescope Array

Spherical harmonic moments are well-suited for capturing anisotropy at any scale in the flux of cosmic rays. An unambiguous measurement of the full set of spherical harmonic coefficients requires full-sky coverage. This can be achieved by combining data from observatories located in both the northern and southern hemispheres. To this end, a joint analysis using data recorded at the Telescope Array and the Pierre Auger Observatory above $10^{19}$ eV is presented in this work. The resulting multipolar expansion of the flux of cosmic rays allows us to perform a series of anisotropy searches, and in particular to report on the angular power spectrum of cosmic rays above $10^{19}$ eV. No significant deviation from isotropic expectations is found throughout the analyses performed. Upper limits on the amplitudes of the dipole and quadrupole moments are derived as a function of the direction in the sky, varying between 7% and 13% for the dipole and between 7% and 10% for a symmetric quadrupole.

preprint2014arXiv

The wavefront of the radio signal emitted by cosmic ray air showers

Analyzing measurements of the LOPES antenna array together with corresponding CoREAS simulations for more than 300 measured events with energy above $10^{17}\,$eV and zenith angles smaller than $45^\circ$, we find that the radio wavefront of cosmic-ray air showers is of approximately hyperbolic shape. The simulations predict a slightly steeper wavefront towards East than towards West, but this asymmetry is negligible against the measurement uncertainties of LOPES. At axis distances $\gtrsim 50\,$m, the wavefront can be approximated by a simple cone. According to the simulations, the cone angle is clearly correlated with the shower maximum. Thus, we confirm earlier predictions that arrival time measurements can be used to study the longitudinal shower development, but now using a realistic wavefront. Moreover, we show that the hyperbolic wavefront is compatible with our measurement, and we present several experimental indications that the cone angle is indeed sensitive to the shower development. Consequently, the wavefront can be used to statistically study the primary composition of ultra-high energy cosmic rays. At LOPES, the experimentally achieved precision for the shower maximum is limited by measurement uncertainties to approximately $140\,$g/cm$^2$. But the simulations indicate that under better conditions this method might yield an accuracy for the atmospheric depth of the shower maximum, $X_\mathrm{max}$, better than $30\,$g/cm$^2$. This would be competitive with the established air-fluorescence and air-Cherenkov techniques, where the radio technique offers the advantage of a significantly higher duty-cycle. Finally, the hyperbolic wavefront can be used to reconstruct the shower geometry more accurately, which potentially allows a better reconstruction of all other shower parameters, too.

preprint2014arXiv

Transiting exoplanets and magnetic spots characterized with optical interferometry

Stellar activity causes difficulties in the characterization of transiting exoplanets. Studies have been performed to quantify its impact on infrared interferometry, but not in the visible domain, which however allows reaching better angular resolution and is also the one mostly used for spectroscopic and photometric measurements. We use a standard case to completely analyse the impact of an exoplanet and a spot on interferometric observables and relate it to current instrument capabilities, taking into account realistic achievable precisions. We built a numerical code called COMETS using analytical formulae to perform a simple comparison of exoplanet and spot signals. We explore instrumental specificities needed to detect them, like the baseline length required, the accuracy and S/N. We also discuss the impact of exoplanet and spot parameters on squared visibility and phase. We find that the main improvement to bring is the sensitivity of instruments. The accuracy on squared visibilities has to be improved by a factor 10 to detect an exoplanet of 0.1 mas, leading to $<0.5%$ precision, along with phase measurements of ~$5^°$ accuracy beyond the first null of visibility. For a 0.05 mas exoplanet, accuracies of ~$0.1%$ and ~$1^°$ from the first null are required on squared visibilities and phases, respectively. Magnetic spots can mimic these signals, leading to false exoplanet characterization. Phases measurements from the 3rd lobe is needed to distinguish between a spot and an exoplanet if they have the same radius. Increasing interferometer sensitivity, more objects will become common between interferometric targets and photometric ones. Furthermore, new missions like PLATO, CHEOPS or TESS will provide bright exoplanet host stars. Measurements will thus overlap and provide a better characterization of stellar activity and exoplanet.

preprint2013arXiv

A Community Science Case for E-ELT HIRES

Building on the experience of the high-resolution community with the suite of VLT high-resolution spectrographs, which has been tremendously successful, we outline here the (science) case for a high-fidelity, high-resolution spectrograph with wide wavelength coverage at the E-ELT. Flagship science drivers include: the study of exo-planetary atmospheres with the prospect of the detection of signatures of life on rocky planets; the chemical composition of planetary debris on the surface of white dwarfs; the spectroscopic study of protoplanetary and proto-stellar disks; the extension of Galactic archaeology to the Local Group and beyond; spectroscopic studies of the evolution of galaxies with samples that, unlike now, are no longer restricted to strongly star forming and/or very massive galaxies; the unraveling of the complex roles of stellar and AGN feedback; the study of the chemical signatures imprinted by population III stars on the IGM during the epoch of reionization; the exciting possibility of paradigm-changing contributions to fundamental physics. The requirements of these science cases can be met by a stable instrument with a spectral resolution of R~100,000 and broad, simultaneous spectral coverage extending from 370nm to 2500nm. Most science cases do not require spatially resolved information, and can be pursued in seeing-limited mode, although some of them would benefit by the E-ELT diffraction limited resolution. Some multiplexing would also be beneficial for some of the science cases. (Abridged)

preprint2013arXiv

Ankle-like Feature in the Energy Spectrum of Light Elements of Cosmic Rays Observed with KASCADE-Grande

Recent results of the KASCADE-Grande experiment provided evidence for a mild knee-like structure in the all-particle spectrum of cosmic rays at $E = 10^{16.92 \pm 0.10} \, \mathrm{eV}$, which was found to be due to a steepening in the flux of heavy primary particles. The spectrum of the combined components of light and intermediate masses was found to be compatible with a single power law in the energy range from $10^{16.3} \, \mathrm{eV}$ to $10^{18} \, \mathrm{eV}$. In this paper, we present an update of this analysis by using data with increased statistics, originating both from a larger data set including more recent measurements and by using a larger fiducial area. In addition, optimized selection criteria for enhancing light primaries are applied. We find a spectral feature for light elements, namely a hardening at $E = 10^{17.08 \pm 0.08} \, \mathrm{eV}$ with a change of the power law index from $-3.25 \pm 0.05$ to $-2.79 \pm 0.08$.

preprint2013arXiv

Comparison of LOPES data and CoREAS simulations using a full detector simulation (ICRC2013)

The LOPES experiment at the Karlsruhe Institute of Technology, Germany, has been measuring radio emission of air showers for almost 10 years. For a better understanding of the emission process a detailed comparison of data with simulations is necessary. This is possible using a newly developed detector simulation including all LOPES detector components. After propagating a simulated event through this full detector simulation a standard LOPES like event file is written. LOPES data and CoREAS simulations can then be treated equally and the same analysis software can be applied to both. This gives the opportunity to compare data and simulations directly. Furthermore, the standard analysis software can be used with simulations which provide the possibility to check the accuracy regarding reconstruction of air shower parameters. We point out the advantages and present first results using such a full LOPES detector simulation. A comparison of LOPES data and the Monte Carlo code CoREAS based on an analysis using this detector simulation is shown.

preprint2013arXiv

Cosmic Ray Measurements with LOPES: Status and Recent Results (ARENA 2012)

LOPES is a digital antenna array at the Karlsruhe Institute of Technology, Germany, for cosmic-ray air-shower measurements. Triggered by the co-located KASCADE-Grande air-shower array, LOPES detects the radio emission of air showers via digital radio interferometry. We summarize the status of LOPES and recent results. In particular, we present an update on the reconstruction of the primary-particle properties based on almost 500 events above 100 PeV. With LOPES, the arrival direction can be reconstructed with a precision of at least 0.65°, and the energy with a precision of at least 20 %, which, however, does not include systematic uncertainties on the absolute energy scale. For many particle and astrophysics questions the reconstruction of the atmospheric depth of the shower maximum, Xmax, is important, since it yields information on the type of the primary particle and its interaction with the atmosphere. Recently, we found experimental evidence that the slope of the radio lateral distribution is indeed sensitive to the longitudinal development of the air shower, but unfortunately, the Xmax precision at LOPES is limited by the high level of anthropogenic radio background. Nevertheless, the developed methods can be transferred to next generation experiments with lower background, which should provide an Xmax precision competitive to other detection technologies.

preprint2013arXiv

Energy spectra of KASCADE-Grande based on shower size measurements and different hadronic interaction models

KASCADE-Grande is dedicated for investigations of cosmic-ray air showers in the primary energy range from 10 PeV to 1 EeV. The multi-detector system allows us to reconstruct charged particles, electron and muon numbers for individual air showers with high accuracies. Based on the shower size ($N_{ch}$) spectra of the charged particle component, the all-particle energy spectrum of cosmic rays is reconstructed, where attenuation effects in the atmosphere are corrected by applying the constant intensity cut method. The energy calibration is performed by using CORSIKA simulations with high-energy interaction models QGSJET-II-2, QGSJET-II-4, EPOS 1.99 and SIBYLL 2.1, where FLUKA has been used as low-energy interaction model for all cases. In the different hadronic models, different abundances for shower particles are predicted. Such model differences in the observables will be compared and discussed in this contribution. Furthermore, by using data with increasing statistics, the updated energy spectra by means of different interaction models will be presented.

preprint2013arXiv

Fundamental parameters of 16 late-type stars derived from their angular diameter measured with VLTI/AMBER

Thanks to their large angular dimension and brightness, red giants and supergiants are privileged targets for optical long-baseline interferometers. Sixteen red giants and supergiants have been observed with the VLTI/AMBER facility over a two-years period, at medium spectral resolution (R=1500) in the K band. The limb-darkened angular diameters are derived from fits of stellar atmospheric models on the visibility and the triple product data. The angular diameters do not show any significant temporal variation, except for one target: TX Psc, which shows a variation of 4% using visibility data. For the eight targets previously measured by Long-Baseline Interferometry (LBI) in the same spectral range, the difference between our diameters and the literature values is less than 5%, except for TX Psc, which shows a difference of 11%. For the 8 other targets, the present angular diameters are the first measured from LBI. Angular diameters are then used to determine several fundamental stellar parameters, and to locate these targets in the Hertzsprung-Russell Diagram (HRD). Except for the enigmatic Tc-poor low-mass carbon star W Ori, the location of Tc-rich stars in the HRD matches remarkably well the thermally-pulsating AGB, as it is predicted by the stellar-evolution models. For pulsating stars with periods available, we compute the pulsation constant and locate the stars along the various sequences in the Period -- Luminosity diagram. We confirm the increase in mass along the pulsation sequences, as predicted by the theory, except for W Ori which, despite being less massive, appears to have a longer period than T Cet along the first-overtone sequence.

preprint2013arXiv

High-energy cosmic rays measured with KASCADE-Grande

The detection of high-energy cosmic rays above a few hundred TeV is realized by the observation of extensive air-showers. By using the multi-detector setup of KASCADE-Grande, energy spectrum, elemental composition, and anisotropies of high-energy cosmic rays in the energy range from below the knee up to 2 EeV are investigated. In addition, the large high-quality data set permits distinct tests of the validity of hadronic interaction models used in interpreting air-shower measurements. After more than 16 years, the KASCADE-Grande experiment terminated measurements end of 2012. This contribution will give an overview of the main results of the data analysis achieved so far, and will report about the status of KCDC, the KASCADE Cosmic-ray Data Center, where via a web-based interface the data will be made available for the interested public.

preprint2013arXiv

Investigation on the energy and mass composition of cosmic rays using LOPES radio data

The sensitivity to the mass composition as well as the reconstruction of the energy of the primary particle are explored here by leveraging the features of the radio lateral distribution function. For the purpose of this analysis, a set of events measured with the LOPES experiment is reproduced with the latest CoREAS radio simulation code. Based on simulation predictions, a method which exploits the slope of the radio lateral distribution function is developed (Slope Method) and directly applied on measurements. As a result, the possibility to reconstruct both the energy and the depth of the shower maximum of the cosmic ray air shower using radio data and achieving relatively small uncertainties is presented.

preprint2013arXiv

KASCADE-Grande measurements of energy spectra for elemental groups of cosmic rays

The KASCADE-Grande air shower experiment [W. Apel, et al. (KASCADE-Grande collaboration), Nucl. Instrum. Methods A 620 (2010) 202] consists of, among others, a large scintillator array for measurements of charged particles, Nch, and of an array of shielded scintillation counters used for muon counting, Nmu. KASCADE-Grande is optimized for cosmic ray measurements in the energy range 10 PeV to about 2000 PeV, where exploring the composition is of fundamental importance for understanding the transition from galactic to extragalactic origin of cosmic rays. Following earlier studies of the all-particle and the elemental spectra reconstructed in the knee energy range from KASCADE data [T. Antoni, et al. (KASCADE collaboration), Astropart. Phys. 24 (2005) 1], we have now extended these measurements to beyond 200 PeV. By analysing the two-dimensional shower size spectrum Nch vs. Nmu for nearly vertical events, we reconstruct the energy spectra of different mass groups by means of unfolding methods over an energy range where the detector is fully efficient. The procedure and its results, which are derived based on the hadronic interaction model QGSJET-II-02 and which yield a strong indication for a dominance of heavy mass groups in the covered energy range and for a knee-like structure in the iron spectrum at around 80 PeV, are presented. This confirms and further refines the results obtained by other analyses of KASCADE-Grande data, which already gave evidence for a knee-like structure in the heavy component of cosmic rays at about 80 PeV [W. Apel, et al. (KASCADE-Grande collaboration), Phys. Rev. Lett. 107 (2011) 171104].

preprint2013arXiv

KASCADE-Grande measurements of energy spectra for elemental groups of cosmic rays

The KASCADE-Grande experiment, located at KIT-Karlsruhe, Germany, consists of a large scintillator array for measurements of charged particles, N_ch, and of an array of shielded scintillation counters used for muon counting, N_mu. KASCADE-Grande is optimized for cosmic ray measurements in the energy range 10 PeV to 1000 PeV, thereby enabling the verification of a knee in the iron spectrum expected at approximately 100 PeV. Exploring the composition in this energy range is of fundamental importance for understanding the transition from galactic to extragalactic cosmic rays. Following earlier studies of elemental spectra reconstructed in the knee energy range from KASCADE data, we have now extended these measurements to beyond 100 PeV. By analysing the two-dimensional shower size spectrum N_ch vs. N_mu, we reconstruct the energy spectra of different mass groups by means of unfolding methods. The procedure and its results, giving evidence for a knee-like structure in the spectrum of iron nuclei, will be presented.

preprint2013arXiv

Long-term spectropolarimetric monitoring of the cool supergiant Betelgeuse

We report on a long-term monitoring of the cool supergiant Betelgeuse, using the NARVAL and ESPaDOnS high-resolution spectropolarimeters, respectively installed at Telescope Bernard Lyot (Pic du Midi Observatory, France) and at the Canada-France-Hawaii Telescope (Mauna Kea Observatory, Hawaii). The data set, constituted of circularly polarized (Stokes V) and intensity (Stokes I) spectra, was collected between 2010 and 2012. We investigate here the temporal evolution of magnetic field, convection and temperature at photospheric level, using simultaneous measurements of the longitudinal magnetic field component, the core emission of the Ca II infrared triplet, the line-depth ratio of selected photospheric lines and the radial velocity of the star.

preprint2013arXiv

LOPES 3D - vectorial measurements of radio emission from cosmic ray induced air showers

LOPES 3D is able to measure all three components of the electric field vector of the radio emission from air showers. This allows a better comparison with emission models. The measurement of the vertical component increases the sensitivity to inclined showers. By measuring all three components of the electric field vector LOPES 3D demonstrates by how much the reconstruction accuracy of primary cosmic ray parameters increases. Thus LOPES 3D evaluates the usefulness of vectorial measurements for large scale applications.

preprint2013arXiv

LOPES 3D reconfiguration and first measurements

The Radio detection technique of high-energy cosmic rays is based on the radio signal emitted by the charged particles in an air shower due to their deflection in the Earth's magnetic field. The LOPES experiment at Karlsruhe Institute of Technology, Germany with its simple dipoles made major contributions to the revival of this technique. LOPES is working in the frequency range from 40 to 80 MHz and was reconfigured several times to improve and further develop the radio detection technique. In the current setup LOPES consists of 10 tripole antennas which measure the complete electric field vector of the radio emission from cosmic rays. LOPES is the first experiment measuring all three vectorial components at once and thereby gaining the full information about the electric field vector and not only a two-dimensional projection. Such a setup including also measurements of the vertical electric field component is expected to increase the sensitivity to inclined showers and help to advance the understanding of the emission mechanism. We present the reconfiguration and calibration procedure of LOPES 3D and discuss first measurements.

preprint2013arXiv

LOPES-3D, an antenna array for full signal detection of air-shower radio emission

To better understand the radio signal emitted by extensive air-showers and to further develop the radio detection technique of high-energy cosmic rays, the LOPES experiment was reconfigured to LOPES-3D. LOPES-3D is able to measure all three vectorial components of the electric field of radio emission from cosmic ray air showers. The additional measurement of the vertical component ought to increase the reconstruction accuracy of primary cosmic ray parameters like direction and energy, provides an improved sensitivity to inclined showers, and will help to validate simulation of the emission mechanisms in the atmosphere. LOPES-3D will evaluate the feasibility of vectorial measurements for large scale applications. In order to measure all three electric field components directly, a tailor-made antenna type (tripoles) was deployed. The change of the antenna type necessitated new pre-amplifiers and an overall recalibration. The reconfiguration and the recalibration procedure are presented and the operationality of LOPES-3D is demonstrated.

preprint2013arXiv

Mass sensitivity in the radio lateral distribution function

Measuring the mass composition of ultra-high energy cosmic rays is one of the main tasks in the cosmic rays field. Here we are exploring the composition signature in the coherent electromagnetic emission from extensive air showers, detected in the MHz frequency range. One of the experiments that successfully detects radio events in the frequency band of 40-80 MHz is the LOPES experiment at KIT. It is a digital interferometric antenna array and has the important advantage of taking data in coincidence with the particle detector array KASCADE-Grande. A possible method to look at the composition signature in the radio data, predicted by simulations, concerns the radio lateral distribution function, since its slope is strongly correlated with Xmax. Recent comparison between REAS3 simulations and LOPES data showed a significantly improved agreement in the lateral distribution function and for this reason an analysis on a possible LOPES mass signature through the slope method is promising. Trying to reproduce a realistic case, proton and iron showers are simulated with REAS3 using the LOPES selection information as input parameters. The obtained radio lateral distribution slope is analyzed in detail. The lateral slope method to look at the composition signature in the radio data is shown here and a possible signature of mass composition in the LOPES data is discussed.

preprint2013arXiv

Observation of Polarised Microwave Emission from Cosmic Ray Air Showers

We report on the first direct measurement of the basic features of microwave radio emission from extensive air showers. Using a trigger provided by the KASCADE-Grande air shower array, the signals of the microwave antennas of the CROME (Cosmic-Ray Observation via Microwave Emission) experiment have been read out and searched for signatures of radio emission by high-energy air showers. Microwave signals have been detected for more than 30 showers with energies above $3\times10^{16}$\,eV. The observations presented in this Letter are consistent with a mainly forward-beamed, coherent and polarised emission process in the GHz frequency range. An isotropic, unpolarised radiation is disfavoured as the dominant emission model. The measurements show that microwave radiation offers a new means of studying air showers at very high energy.

preprint2013arXiv

Reconstructing energy and Xmax of cosmic ray air showers using the radio lateral distribution measured with LOPES

The LOPES experiment, a digital radio interferometer located at KIT (Karlsruhe Institute of Technology), obtained remarkable results for the detection of radio emission from extensive air showers at MHz frequencies. Features of the radio lateral distribution function (LDF) measured by LOPES are explored in this work for a precise reconstruction of two fundamental air shower parameters: the primary energy and the shower Xmax. The method presented here has been developed on (REAS3-)simulations, and is applied to LOPES measurements. Despite the high human-made noise at the LOPES site, it is possible to reconstruct both the energy and Xmax for individual events. On the one hand, the energy resolution is promising and comparable to the one of the co-located KASCADE-Grande experiment. On the other hand, Xmax values are reconstructed with the LOPES measurements with a resolution of 90 g/cm2 . A precision on Xmax better than 30 g/cm2 is predicted and achievable in a region with a lower human-made noise level.

preprint2013arXiv

Red supergiant star studies with CO5BOLD and Optim3D

We describe recent work focused towards a better understanding of red supergiant stars using 3D radiative-hydrodynamics (RHD) simulations with CO5BOLD. A small number of simulations now exist that span up to seven years of stellar time, at various numerical resolutions. Our discussion concentrates on interferometric and spectroscopic observations. We point out a number of problems, in particular the line depth and line width that are not well reproduced by simulations. The most recent introduction of a non-grey treatment of the radiation field dramatically improved the match with observations, without solving all difficulties. We also review the newly revived effective temperature scale controversy, and argue that it will only be solved using 3D RHD models.

preprint2013arXiv

Study of the muon content of very high-energy EAS measured with the KASCADE-Grande observatory

The KASCADE-Grande detector is an air-shower array devoted to the study of primary cosmic rays with very high-energies (E = 10^{16} - 10^{18} eV). The instrument is composed of different particle detector systems suitable for the detailed study of the properties of Extensive Air Showers (EAS) developed by cosmic rays in the atmosphere. Among the EAS observables studied with the detector, the charged number of particles, the muon content (at different energy thresholds), and the number of electrons are found. By comparing the measurements of these air-shower parameters with the expectations from MC simulations, different hadronic interaction models can be tested at the high-energy regime with the KASCADE-Grande experiment. In this work, the results of a study on the evolution of the muon content of EAS with zenith angle, performed with the KASCADE-Grande instrument, is presented. Measurements are compared with predictions from MC simulations based on the QGSJET II, QGSJET II-04, SIBYLL 2.1 and EPOS 1.99 hadronic interaction models. A mismatch between experiment and simulations is observed. A similar problem is found for the evolution of the lateral distribution function of muons in the atmosphere.

preprint2013arXiv

The <lnA> study in the primary energy range 10^{16} - 10^{17} eV with the Muon Tracking Detector in the KASCADE-Grande experiment

The KASCADE-Grande Muon Tracking Detector enables with high accuracy the measurement of directions of EAS muons with energy above 0.8 GeV and up to 700 m distance from the shower centre. Reconstructed muon tracks are used to investigate muon pseudorapidity (eta) distributions. These distributions are nearly identical to the pseudorapidity distributions of their parent mesons produced in hadronic interactions. Comparison of the eta distributions from measured and simulated showers can be used to test the quality of the high energy hadronic interaction models. In this context a comparison of the QGSJet-II-2 and QGSJet-II-4 model will be shown. The pseudorapidity distributions reflect the longitudinal development of EAS and, as such, are sensitive to the mass of the cosmic rays primary particles. With various parameters of the eta distribution, obtained from the MTD data, it is possible to calculate the mean logarithmic mass of CRs. The results of the <lnA> analysis in the primary energy range 10^{16} eV - 10^{17} eV with the 1st quartile (Q1) of eta distribution will be presented.

preprint2013arXiv

The energy-spectrum of light primaries in the range from 10^{16.6} to 10^{18.2} eV

Data of the Grande extension of the KASCADE experiment allows us to study extensive air showers induced by primary cosmic rays with energies above 10^{16} eV. The energy of an event is estimated in terms of the number of charged particles (Nch ) and the number of muons (Nμ ) measured at an altitude of 110 m a.s.l. While a combination of the two numbers is used for the energy, the ratio defines the primary mass (group). The spectrum of the combined light and medium mass components, recently measured with KASCADE-Grande, was found to be compatible with both a single power-law and a broken power-law in the energy range between 10^{16.3} and 10^{18} eV. In this contribution we will present the investigation of possible structures in the spectrum of light primaries with increased statistics both from a larger data set including more recent measurements and by using a larger fiducial area than in the previous study. With the better statistical accuracy and with optimized selection criteria for enhancing light primaries we have found evidence for a hardening (ankle) of the spectrum of the light component at an energy of 10^{17.08+-0.08} eV.

preprint2013arXiv

The Stagger-grid: A Grid of 3D Stellar Atmosphere Models - I. Methods and General Properties

We present the Stagger-grid, a comprehensive grid of time-dependent, 3D hydrodynamic model atmospheres for late-type stars with realistic treatment of radiative transfer, covering a wide range in stellar parameters. This grid of 3D models is intended for various applications like stellar spectroscopy, asteroseismology and the study of stellar convection. In this introductory paper, we describe the methods used for the computation of the grid and discuss the general properties of the 3D models as well as their temporal and spatial averages (<3D>). All our models were generated with the Stagger-code, using realistic input physics for the equation of state (EOS) and for continuous and line opacities. Our ~220 grid models range in Teff from 4000 to 7000K in steps of 500K, in log g from 1.5 to 5.0 in steps of 0.5 dex, and [Fe/H] from -4.0 to +0.5 in steps of 0.5 and 1.0 dex. We find a tight scaling relation between the vertical velocity and the surface entropy jump, which itself correlates with the constant entropy value of the adiabatic convection zone. The range in intensity contrast is enhanced at lower metallicity. The granule size correlates closely with the pressure scale height sampled at the depth of maximum velocity. We compare the <3D> models with widely applied 1D models, as well as with theoretical 1D hydrostatic models generated with the same EOS and opacity tables as the 3D models, in order to isolate the effects of using self-consistent and hydrodynamic modeling of convection, rather than the classical mixing length theory approach. For the first time, we are able to quantify systematically over a broad range of stellar parameters the uncertainties of 1D models arising from the simplified treatment of physics, in particular convective energy transport. In agreement with previous findings, we find that the differences can be significant, especially for metal-poor stars.

preprint2013arXiv

Thunderstorm Observations by Air-Shower Radio Antenna Arrays

Relativistic, charged particles present in extensive air showers lead to a coherent emission of radio pulses which are measured to identify the shower initiating high-energy cosmic rays. Especially during thunderstorms, there are additional strong electric fields in the atmosphere, which can lead to further multiplication and acceleration of the charged particles and thus have influence on the form and strength of the radio emission. For a reliable energy reconstruction of the primary cosmic ray by means of the measured radio signal it is very important to understand how electric fields affect the radio emission. In addition, lightning strikes are a prominent source of broadband radio emissions that are visible over very long distances. This, on the one hand, causes difficulties in the detection of the much lower signal of the air shower. On the other hand the recorded signals can be used to study features of the lightning development. The detection of cosmic rays via the radio emission and the influence of strong electric fields on this detection technique is investigated with the LOPES experiment in Karlsruhe, Germany. The important question if a lightning is initiated by the high electron density given at the maximum of a high-energy cosmic-ray air shower is also investigated, but could not be answered by LOPES. But, these investigations exhibit the capabilities of EAS radio antenna arrays for lightning studies. We report about the studies of LOPES measured radio signals of air showers taken during thunderstorms and give a short outlook to new measurements dedicated to search for correlations of lightning and cosmic rays.

preprint2013arXiv

Vectorial Radio Interferometry with LOPES 3D

One successful detection technique for high-energy cosmic rays is based on the radio signal emitted by the charged particles in an air shower. The LOPES experiment at Karlsruhe Institute of Technology, Germany, has made major contributions to the evolution of this technique. LOPES was reconfigured several times to improve and further develop the radio detection technique. In the latest setup LOPES consisted of 10 tripole antennas. With this, LOPES 3D was the first cosmic ray experiment measuring all three vectorial field components at once and thereby gaining the full information about the electric field vector. We present an analysis based on the data taken with special focus on the benefits of a direct measurement of the vertical polarization component. We demonstrate that by measuring all polarization components the detection and reconstruction efficiency is increased and noisy single channel data can be reconstructed by utilising the information from the other two channels of one antenna station.

preprint2012arXiv

A search for ultra-high energy neutrinos in highly inclined events at the Pierre Auger Observatory

The Surface Detector of the Pierre Auger Observatory is sensitive to neutrinos of all flavours above 0.1 EeV. These interact through charged and neutral currents in the atmosphere giving rise to extensive air showers. When interacting deeply in the atmosphere at nearly horizontal incidence, neutrinos can be distinguished from regular hadronic cosmic rays by the broad time structure of their shower signals in the water-Cherenkov detectors. In this paper we present for the first time an analysis based on down-going neutrinos. We describe the search procedure, the possible sources of background, the method to compute the exposure and the associated systematic uncertainties. No candidate neutrinos have been found in data collected from 1 January 2004 to 31 May 2010. Assuming an E^-2 differential energy spectrum the limit on the single flavour neutrino is (E^2 * dN/dE) < 1.74x10^-7 GeV cm^-2 s^-1 sr^-1 at 90% C.L. in the energy range 1x10^17 eV < E < 1x10^20 eV.

preprint2012arXiv

Antennas for the Detection of Radio Emission Pulses from Cosmic-Ray induced Air Showers at the Pierre Auger Observatory

The Pierre Auger Observatory is exploring the potential of the radio detection technique to study extensive air showers induced by ultra-high energy cosmic rays. The Auger Engineering Radio Array (AERA) addresses both technological and scientific aspects of the radio technique. A first phase of AERA has been operating since September 2010 with detector stations observing radio signals at frequencies between 30 and 80 MHz. In this paper we present comparative studies to identify and optimize the antenna design for the final configuration of AERA consisting of 160 individual radio detector stations. The transient nature of the air shower signal requires a detailed description of the antenna sensor. As the ultra-wideband reception of pulses is not widely discussed in antenna literature, we review the relevant antenna characteristics and enhance theoretical considerations towards the impulse response of antennas including polarization effects and multiple signal reflections. On the basis of the vector effective length we study the transient response characteristics of three candidate antennas in the time domain. Observing the variation of the continuous galactic background intensity we rank the antennas with respect to the noise level added to the galactic signal.

preprint2012arXiv

Description of Atmospheric Conditions at the Pierre Auger Observatory using the Global Data Assimilation System (GDAS)

Atmospheric conditions at the site of a cosmic ray observatory must be known for reconstructing observed extensive air showers. The Global Data Assimilation System (GDAS) is a global atmospheric model predicated on meteorological measurements and numerical weather predictions. GDAS provides altitude-dependent profiles of the main state variables of the atmosphere like temperature, pressure, and humidity. The original data and their application to the air shower reconstruction of the Pierre Auger Observatory are described. By comparisons with radiosonde and weather station measurements obtained on-site in Malargüe and averaged monthly models, the utility of the GDAS data is shown.

preprint2012arXiv

Experimental evidence for the sensitivity of the air-shower radio signal to the longitudinal shower development

We observe a correlation between the slope of radio lateral distributions, and the mean muon pseudorapidity of 59 individual cosmic-ray-air-shower events. The radio lateral distributions are measured with LOPES, a digital radio interferometer co-located with the multi-detector-air-shower array KASCADE-Grande, which includes a muon-tracking detector. The result proves experimentally that radio measurements are sensitive to the longitudinal development of cosmic-ray air-showers. This is one of the main prerequisites for using radio arrays for ultra-high-energy particle physics and astrophysics.

preprint2012arXiv

Fundamental properties of the Population II fiducial stars HD 122563 and Gmb 1830 from CHARA interferometric observations

We have determined the angular diameters of two metal-poor stars, HD 122563 and Gmb 1830, using CHARA and Palomar Testbed Interferometer observations. For the giant star HD 122563, we derive an angular diameter theta_3D = 0.940 +- 0.011 milliarcseconds (mas) using limb-darkening from 3D convection simulations and for the dwarf star Gmb 1830 (HD 103095) we obtain a 1D limb-darkened angular diameter theta_1D = 0.679 +- 0.007 mas. Coupling the angular diameters with photometry yields effective temperatures with precisions better than 55 K (Teff = 4598 +- 41 K and 4818 +- 54 K --- for the giant and the dwarf star, respectively). Including their distances results in very well-determined luminosities and radii (L = 230 +- 6 L_sun, R = 23.9 +- 1.9 R_sun and L = 0.213 +- 0.002 L_sun, R = 0.664 +- 0.015 R_sun, respectively). We used the CESAM2k stellar structure and evolution code in order to produce models that fit the observational data. We found values of the mixing-length parameter alpha (which describes 1D convection) that depend on the mass of the star. The masses were determined from the models with precisions of <3% and with the well-measured radii excellent constraints on the surface gravity are obtained (log g = 1.60 +- 0.04, 4.59 +- 0.02, respectively). The very small errors on both log g and Teff provide stringent constraints for spectroscopic analyses given the sensitivity of abundances to both of these values. The precise determination of Teff for the two stars brings into question the photometric scales for metal-poor stars.

preprint2012arXiv

Imaging the heart of astrophysical objects with optical long-baseline interferometry

The number of publications of aperture-synthesis images based on optical long-baseline interferometry measurements has recently increased due to easier access to visible and infrared interferometers. The interferometry technique has now reached a technical maturity level that opens new avenues for numerous astrophysical topics requiring milli-arcsecond model-independent imaging. In writing this paper our motivation was twofold: 1) review and publicize emblematic excerpts of the impressive corpus accumulated in the field of optical interferometry image reconstruction; 2) discuss future prospects for this technique by selecting four representative astrophysical science cases in order to review the potential benefits of using optical long baseline interferometers. For this second goal we have simulated interferometric data from those selected astrophysical environments and used state-of-the-art codes to provide the reconstructed images that are reachable with current or soon-to-be facilities. The image reconstruction process was "blind" in the sense that reconstructors had no knowledge of the input brightness distributions. We discuss the impact of optical interferometry in those four astrophysical fields. We show that image reconstruction software successfully provides accurate morphological information on a variety of astrophysical topics and review the current strengths and weaknesses of such reconstructions. We investigate how to improve image reconstruction and the quality of the image possibly by upgrading the current facilities. We finally argue that optical interferometers and their corresponding instrumentation, existing or to come, with 6 to 10 telescopes, should be well suited to provide images of complex sceneries.

preprint2012arXiv

Radar reflection off extensive air showers

We investigate the possibility of detecting extensive air showers by the radar technique. Considering a bistatic radar system and different shower geometries, we simulate reflection of radio waves off the static plasma produced by the shower in the air. Using the Thomson cross-section for radio wave reflection, we obtain the time evolution of the signal received by the antennas. The frequency upshift of the radar echo and the power received are studied to verify the feasibility of the radar detection technique.

preprint2012arXiv

The Rapid Atmospheric Monitoring System of the Pierre Auger Observatory

The Pierre Auger Observatory is a facility built to detect air showers produced by cosmic rays above 10^17 eV. During clear nights with a low illuminated moon fraction, the UV fluorescence light produced by air showers is recorded by optical telescopes at the Observatory. To correct the observations for variations in atmospheric conditions, atmospheric monitoring is performed at regular intervals ranging from several minutes (for cloud identification) to several hours (for aerosol conditions) to several days (for vertical profiles of temperature, pressure, and humidity). In 2009, the monitoring program was upgraded to allow for additional targeted measurements of atmospheric conditions shortly after the detection of air showers of special interest, e.g., showers produced by very high-energy cosmic rays or showers with atypical longitudinal profiles. The former events are of particular importance for the determination of the energy scale of the Observatory, and the latter are characteristic of unusual air shower physics or exotic primary particle types. The purpose of targeted (or "rapid") monitoring is to improve the resolution of the atmospheric measurements for such events. In this paper, we report on the implementation of the rapid monitoring program and its current status. The rapid monitoring data have been analyzed and applied to the reconstruction of air showers of high interest, and indicate that the air fluorescence measurements affected by clouds and aerosols are effectively corrected using measurements from the regular atmospheric monitoring program. We find that the rapid monitoring program has potential for supporting dedicated physics analyses beyond the standard event reconstruction.

preprint2012arXiv

The spectrum of high-energy cosmic rays measured with KASCADE-Grande

The energy spectrum of cosmic rays between 10**16 eV and 10**18 eV, derived from measurements of the shower size (total number of charged particles) and the total muon number of extensive air showers by the KASCADE-Grande experiment, is described. The resulting all-particle energy spectrum exhibits strong hints for a hardening of the spectrum at approximately 2x10**16 eV and a significant steepening at c. 8x10**16 eV. These observations challenge the view that the spectrum is a single power law between knee and ankle. Possible scenarios generating such features are discussed in terms of astrophysical processes that may explain the transition region from galactic to extragalactic origin of cosmic rays.

preprint2012arXiv

The Tunka-133 EAS Cherenkov light array: status of 2011

A new EAS Cherenkov light array, Tunka-133, with ~1 km^2 geometrical area has been installed at the Tunka Valley (50 km from Lake Baikal) in 2009. The array permits a detailed study of cosmic ray energy spectrum and mass composition in the energy range 10^16 - 10^18 eV with a uniform method. We describe the array construction, DAQ and methods of the array calibration.The method of energy reconstruction and absolute calibration of measurements are discussed. The analysis of spatial and time structure of EAS Cherenkov light allows to estimate the depth of the EAS maximum X_max. The results on the all particles energy spectrum and the mean depth of the EAS maximum X_max vs. primary energy derived from the data of two winter seasons (2009 -- 2011), are presented. Preliminary results of joint operation of the Cherenkov array with antennas for detection of EAS radio signals are shown. Plans for future upgrades -- deployment of remote clusters, radioantennas and a scintillator detector network and a prototype of the HiSCORE gamma-telescope -- are discussed.

preprint2012arXiv

Three-dimensional interferometric, spectrometric, and planetary views of Procyon

We used a new realistic 3D radiative-hydrodynamical model atmosphere of Procyon generated with the Stagger Code and synthetic spectra computed with the radiative transfer code Optim3D to re-analyze interferometric and spectroscopic data from the optical to the infrared of Procyon. We compute intensity maps in two optical filters centered at 500 and 800 nm (MARK III) and one infrared filter centered at 2200 nm (VINCI). We constructed stellar disk images accounting for the center-to-limb variations and used them to derive visibility amplitudes and closure phases. We provide 3D limb-darkening coefficients in the optical as well as in the infrared. We show that visibility curves and closure phases show clear deviations from circular symmetry from the 3rd lobe on. These deviations are detectable with current interferometers using closure phases. We derive new angular diameters at different wavelengths with two independent methods based on 3D simulations. We find a diameter_Vinci = 5.390 \pm 0.03 mas that this is confirmed by an independent asteroseismic estimation. The resulting Teff is 6591 K, which is consistent with the infrared flux method determinations. We find also a value of the surface gravity log g = 4.01 \pm 0.03 that is larger by 0.05 dex from literature values. Spectrophotometric comparisons with observations provide very good agreement with the spectral energy distribution and photometric colors, allowing us to conclude that the thermal gradient of the simulation matches fairly well Procyon. Finally, we show that the granulation pattern of a planet hosting Procyon-like star has a non-negligible impact on the detection of hot Jupiters in the infrared using interferometry closure phases. It is then crucial to have a comprehensive knowledge of the host star to directly detect and characterize hot Jupiters. In this respect, RHD simulations are very important to reach this aim.

preprint2011arXiv

3D hydrodynamical model atmospheres: a tool to correct radial velocities and parallaxes for Gaia

Convection plays an essential role in the emerging intensity for many stars that will be observed by Gaia. Convective-related surface structures affect the shape, shift, and asymmetry of absorption lines, the phocentric and photometric variability causing bias in Gaia measurements. Regarding the importance of Gaia mission and its goals, it is mandatory to have the best models of the observed stars. 3D time-dependent hydrodynamical simulations of surface convection are crucial to model the photosphere of late type stars in a very realistic way. These simulations are an important tool to correct the radial velocities and better estimates the parallaxes and photometric variability.

preprint2011arXiv

Advanced functionality for radio analysis in the Offline software framework of the Pierre Auger Observatory

The advent of the Auger Engineering Radio Array (AERA) necessitates the development of a powerful framework for the analysis of radio measurements of cosmic ray air showers. As AERA performs "radio-hybrid" measurements of air shower radio emission in coincidence with the surface particle detectors and fluorescence telescopes of the Pierre Auger Observatory, the radio analysis functionality had to be incorporated in the existing hybrid analysis solutions for fluoresence and surface detector data. This goal has been achieved in a natural way by extending the existing Auger Offline software framework with radio functionality. In this article, we lay out the design, highlights and features of the radio extension implemented in the Auger Offline framework. Its functionality has achieved a high degree of sophistication and offers advanced features such as vectorial reconstruction of the electric field, advanced signal processing algorithms, a transparent and efficient handling of FFTs, a very detailed simulation of detector effects, and the read-in of multiple data formats including data from various radio simulation codes. The source code of this radio functionality can be made available to interested parties on request.

preprint2011arXiv

Anisotropy and chemical composition of ultra-high energy cosmic rays using arrival directions measured by the Pierre Auger Observatory

The Pierre Auger Collaboration has reported evidence for anisotropy in the distribution of arrival directions of the cosmic rays with energies $E>E_{th}=5.5\times 10^{19}$ eV. These show a correlation with the distribution of nearby extragalactic objects, including an apparent excess around the direction of Centaurus A. If the particles responsible for these excesses at $E>E_{th}$ are heavy nuclei with charge $Z$, the proton component of the sources should lead to excesses in the same regions at energies $E/Z$. We here report the lack of anisotropies in these directions at energies above $E_{th}/Z$ (for illustrative values of $Z=6,\ 13,\ 26$). If the anisotropies above $E_{th}$ are due to nuclei with charge $Z$, and under reasonable assumptions about the acceleration process, these observations imply stringent constraints on the allowed proton fraction at the lower energies.

preprint2011arXiv

Kneelike structure in the spectrum of the heavy component of cosmic rays observed with KASCADE-Grande

We report the observation of a steepening in the cosmic ray energy spectrum of heavy primary particles at about $8 \times 10^{16}$ eV. This structure is also seen in the all-particle energy spectrum, but is less significant. Whereas the `knee' of the cosmic ray spectrum at $3-5 \times 10^{15}$ eV was assigned to light primary masses by the KASCADE experiment, the new structure found by the KASCADE-Grande experiment is caused by heavy primaries. The result is obtained by independent measurements of the charged particle and muon components of the secondary particles of extensive air showers in the primary energy range of $10^{16}$ to $10^{18}$ eV. The data are analyzed on a single-event basis taking into account also the correlation of the two observables.

preprint2011arXiv

Radiative hydrodynamics simulations of red supergiant stars. IV gray versus non-gray opacities

Red supergiants are massive evolved stars that contribute extensively to the chemical enrichment of our Galaxy. It has been shown that convection in those stars gives rise to large granules that cause surface inhomogeneities and shock waves in the photosphere. The understanding of their dynamics is crucial to unveil the unknown mass-loss mechanism, their chemical composition and stellar parameters. We present a new generation of red supergiants simulations with a more sophisticated opacity treatment done with 3D radiative- hydrodynamics CO5BOLD. In the code, the coupled equations of compressible hydrodynamics and non-local radiation transport are solved in the presence of a spherical potential. The stellar core is replaced by a special spherical inner boundary condition, where the gravitational potential is smoothed and the energy production by fusion is mimicked by a simply producing heat corresponding to the stellar luminosity. The post-processing radiative transfer code OPTIM3D is used to extract spectroscopic and interferometric observables. We show that the relaxation of the assumption of frequency-independent opacities shows a steeper mean thermal gradient in the optical thin region that affect strongly the atomic strengths and the spectral energy distribution. Moreover, the weaker temperature fluctuations reduce the incertitude on the radius determination with interferometry. We show that 1D models of red supergiants must include a turbulent velocity calibrated on 3D simulations to obtain the effective surface gravity that mimic the effect of turbulent pressure on the stellar atmosphere. We provide an empirical calibration of the ad-hoc micro- and macroturbulence parameters for 1D models using the 3D simulations: we find that there is not a clear distinction between the different macroturbulent profiles needed in 1D models to fit 3D synthetic lines.

preprint2011arXiv

Search for First Harmonic Modulation in the Right Ascension Distribution of Cosmic Rays Detected at the Pierre Auger Observatory

We present the results of searches for dipolar-type anisotropies in different energy ranges above $2.5\times 10^{17}$ eV with the surface detector array of the Pierre Auger Observatory, reporting on both the phase and the amplitude measurements of the first harmonic modulation in the right-ascension distribution. Upper limits on the amplitudes are obtained, which provide the most stringent bounds at present, being below 2% at 99% $C.L.$ for EeV energies. We also compare our results to those of previous experiments as well as with some theoretical expectations.

preprint2011arXiv

Search for signatures of magnetically-induced alignment in the arrival directions measured by the Pierre Auger Observatory

We present the results of an analysis of data recorded at the Pierre Auger Observatory in which we search for groups of directionally-aligned events (or `multiplets') which exhibit a correlation between arrival direction and the inverse of the energy. These signatures are expected from sets of events coming from the same source after having been deflected by intervening coherent magnetic fields. The observation of several events from the same source would open the possibility to accurately reconstruct the position of the source and also measure the integral of the component of the magnetic field orthogonal to the trajectory of the cosmic rays. We describe the largest multiplets found and compute the probability that they appeared by chance from an isotropic distribution. We find no statistically significant evidence for the presence of multiplets arising from magnetic deflections in the present data.

preprint2011arXiv

The effect of the geomagnetic field on cosmic ray energy estimates and large scale anisotropy searches on data from the Pierre Auger Observatory

We present a comprehensive study of the influence of the geomagnetic field on the energy estimation of extensive air showers with a zenith angle smaller than $60^\circ$, detected at the Pierre Auger Observatory. The geomagnetic field induces an azimuthal modulation of the estimated energy of cosmic rays up to the ~2% level at large zenith angles. We present a method to account for this modulation of the reconstructed energy. We analyse the effect of the modulation on large scale anisotropy searches in the arrival direction distributions of cosmic rays. At a given energy, the geomagnetic effect is shown to induce a pseudo-dipolar pattern at the percent level in the declination distribution that needs to be accounted for.

preprint2011arXiv

The Lateral Trigger Probability function for the Ultra-High Energy Cosmic Ray Showers detected by the Pierre Auger Observatory

In this paper we introduce the concept of Lateral Trigger Probability (LTP) function, i.e., the probability for an extensive air shower (EAS) to trigger an individual detector of a ground based array as a function of distance to the shower axis, taking into account energy, mass and direction of the primary cosmic ray. We apply this concept to the surface array of the Pierre Auger Observatory consisting of a 1.5 km spaced grid of about 1600 water Cherenkov stations. Using Monte Carlo simulations of ultra-high energy showers the LTP functions are derived for energies in the range between 10^{17} and 10^{19} eV and zenith angles up to 65 degs. A parametrization combining a step function with an exponential is found to reproduce them very well in the considered range of energies and zenith angles. The LTP functions can also be obtained from data using events simultaneously observed by the fluorescence and the surface detector of the Pierre Auger Observatory (hybrid events). We validate the Monte-Carlo results showing how LTP functions from data are in good agreement with simulations.

preprint2011arXiv

Trigger and Aperture of the Surface Detector Array of the Pierre Auger Observatory

The surface detector array of the Pierre Auger Observatory consists of 1600 water-Cherenkov detectors, for the study of extensive air showers (EAS) generated by ultra-high-energy cosmic rays. We describe the trigger hierarchy, from the identification of candidate showers at the level of a single detector, amongst a large background (mainly random single cosmic ray muons), up to the selection of real events and the rejection of random coincidences. Such trigger makes the surface detector array fully efficient for the detection of EAS with energy above $3\times 10^{18}$ eV, for all zenith angles between 0$^\circ$ and 60$^\circ$, independently of the position of the impact point and of the mass of the primary particle. In these range of energies and angles, the exposure of the surface array can be determined purely on the basis of the geometrical acceptance.

preprint2010arXiv

Longitudinal EAS-Development Studies in the Air-Shower Experiment KASCADE-Grande

A large area (128 m^2) Muon Tracking Detector (MTD), located within the KASCADE experiment, has been built with the aim to identify muons (E_mu > 0.8 GeV) and their directions in extensive air showers by track measurements under more than 18 r.l. shielding. The orientation of the muon track with respect to the shower axis is expressed in terms of the radial- and tangential angles. By means of triangulation the muon production height H_mu is determined. By means of H_mu, a transition from light to heavy cosmic ray primary particle with increasing shower energy Eo from 1-10 PeV is observed. Muon pseudorapidity distributions for the first interactions above 15 km are studied and compared to Monte Carlo simulations.

preprint2010arXiv

On noise treatment in radio measurements of cosmic ray air showers

Precise measurements of the radio emission by cosmic ray air showers require an adequate treatment of noise. Unlike to usual experiments in particle physics, where noise always adds to the signal, radio noise can in principle decrease or increase the signal if it interferes by chance destructively or constructively. Consequently, noise cannot simply be subtracted from the signal, and its influence on amplitude and time measurement of radio pulses must be studied with care. First, noise has to be determined consistently with the definition of the radio signal which typically is the maximum field strength of the radio pulse. Second, the average impact of noise on radio pulse measurements at individual antennas is studied for LOPES. It is shown that a correct treatment of noise is especially important at low signal-to-noise ratios: noise can be the dominant source of uncertainty for pulse height and time measurements, and it can systematically flatten the slope of lateral distributions. The presented method can also be transfered to other experiments in radio and acoustic detection of cosmic rays and neutrinos.

preprint2010arXiv

Radiative hydrodynamics simulations of red supergiant stars: II. simulations of convection on Betelgeuse match interferometric observations

Context. The red supergiant (RSG) Betelgeuse is an irregular variable star. Convection may play an important role in understanding this variability. Interferometric observations can be interpreted using sophisticated simulations of stellar convection. Aims. We compare the visibility curves and closure phases obtained from our 3D simulation of RSG convection with CO5BOLD to various interferometric observations of Betelgeuse from the optical to the H band in order to characterize and measure the convection pattern on this star. Methods. We use 3D radiative-hydrodynamics (RHD) simulation to compute intensity maps in different filters and we thus derive interferometric observables using the post-processing radiative transfer code OPTIM3D. The synthetic visibility curves and closure phases are compared to observations. Results. We provide a robust detection of the granulation pattern on the surface of Betelgeuse in the optical and in the H band based on excellent fits to the observed visibility points and closure phases. Moreover, we determine that the Betelgeuse surface in the H band is covered by small to medium scale (5-15 mas) convection-related surface structures and a large (30 mas) convective cell. In this spectral region, H2O molecules are the main absorbers and contribute to the small structures and to the position of the first null of the visibility curve (i.e. the apparent stellar radius).

preprint2010arXiv

Radiative hydrodynamics simulations of red supergiant stars. III. Spectro-photocentric variability, photometric variability, and consequences on Gaia measurements

Context. It has been shown that convection in red supergiant stars gives rise to large granules causing surface inhomogeneities together with shock waves in the photosphere. The resulting motion of the photocenter (on time scales ranging from months to years) could possibly have adverse effects on the parallax determination with Gaia. Aims. We explore the impact of the granulation on the photocentric and photometric variability. We quantify these effects in order to better characterize the error possibly altering the parallax. Methods. We use 3D radiative-hydrodynamics simulations of convection with CO5BOLD and the post-processing radiative transfer code OPTIM3D to compute intensity maps and spectra in the Gaia G band [325-1030 nm]. Results. We provide astrometric and photometric predictions from 3D simulations of RSGs that are used to evaluate the degradation of the astrometric parameters of evolved stars derived by Gaia. We show from RHD simulations that a supergiant like Betelgeuse exhibits a photocentric noise characterised by a standard deviation of the order of 0.1 AU. The number of bright giant and supergiant stars whose Gaia parallaxes will be altered by the photocentric noise ranges from a few tens to several thousandths. The degradation of the astrometric fit due to the presence of this photocentric noise will be noticeable up to about 5 kpc for the brightest supergiants. Moreover, parallaxes of supergiants are affected by a error of the order of a few percents. We show that the photocentric noise, as predicted by the 3D simulation, does account for a substantial part of the supplementary 'cosmic noise' that affects Hipparcos measurements of Betelgeuse and Antares.

preprint2010arXiv

Sensitivity of KASCADE-Grande data to hadronic interaction models

KASCADE-Grande is a large detector array dedicated for studies of high-energy cosmic rays in the primary energy range from 100 TeV to 1 EeV. The multi-detector concept of the experimental set-up offers the possibility to measure simultaneously various observables related to the electromagnetic, muonic, and hadronic air shower components. The experimental data are compared to predictions of CORSIKA simulations using high-energy hadronic interaction models (e.g. QGSJET or EPOS), as well as low-energy interaction models (e.g. FLUKA or GHEISHA). This contribution will summarize the results of such investigations. In particular, the validity of the new EPOS version 1.99 for EAS with energy around 100 PeV will be discussed.

preprint2010arXiv

The KASCADE-Grande experiment: measurements of the all-particle energy spectrum of cosmic rays

The all-particle energy spectrum as measured by the KASCADE-Grande experiment for E = 10^{16} - 10^{18} eV is presented within the framework of the QGSJET II/FLUKA hadronic interaction models. Three different methods were applied based on the muon size and the total number of charged particles individually and in combination. From the study it is found that the spectrum cannot be completely described by a smooth power law due to the presence of characteristic features.

preprint2010arXiv

The proton-air inelastic cross-section measurement at $\sqrt{s} \approx $ 2 TeV from EAS-TOP experiment

The proton-air inelastic cross section value \sigmapairin=338$\pm$21({\it stat})$\pm$19({\it syst})-28({\it syst}) mb at $\sqrt{s} \approx $ 2 TeV has been measured by the EAS-TOP Extensive Air Shower experiment. The absorption length of cosmic ray proton primaries cascades reaching the maximum development at the observation level is obtained from the flux attenuation for different zenith angles (i.e. atmospheric depths). The analysis, including the effects of the heavier primaries contribution and systematic uncertainties, is described. The experimental result is compared with different high energy interaction models and the relationships with the {\it pp} ($\bar pp$) total cross section measurements are discussed.

preprint2010arXiv

Three-dimensional hydrodynamical simulations of red giant stars: semi-global models for the interpretation of interferometric observations

Context. Theoretical predictions from models of red giant branch stars are a valuable tool for various applications in astrophysics ranging from galactic chemical evolution to studies of exoplanetary systems. Aims. We use the radiative transfer code OPTIM3D and realistic 3D radiative-hydrodynamical (RHD) surface convection simulations of red giants to explore the impact of granulation on interferometric observables. Methods. We compute intensity maps for the 3D simulation snapshots in two filters: in the optical at 5000 \pm 300 Å and in the K band 2.14 $\pm$ 0.26 μm FLUOR filter, corresponding to the wavelength-range of instruments mounted on the CHARA interferometer. From the intensity maps, we construct images of the stellar disks, accounting for center-to-limb variations. We then derive interferometric visibility amplitudes and phases. We study their behavior with position angle and wavelength. Results. We provide average limb-darkening coefficients for different metallicities and wavelength-ranges. We detail the prospects for the detection and characterization of granulation and center-to-limb variations of red giant stars with today's interferometers. We find that the effect of convective-related surface structures depends on metallicity and surface gravity. We provided theoretical closure phases that should be incorporated into the analysis of red giant planet companion closure phase signals. We estimate 3D-1D corrections to stellar radii determination: 3D models are ~ 3.5% smaller to ~ 1% larger in the optical with respect to 1D, and roughly 0.5 to 1.5% smaller in the infrared. Even if these corrections are small, they are important to properly set the zero point of effective temperature scale derived by interferometry and to strengthen the confidence of existing red giant catalogues of calibrating stars for interferometry.

preprint2010arXiv

Update on the correlation of the highest energy cosmic rays with nearby extragalactic matter

Data collected by the Pierre Auger Observatory through 31 August 2007 showed evidence for anisotropy in the arrival directions of cosmic rays above the Greisen-Zatsepin-Kuz'min energy threshold, \nobreak{$6\times 10^{19}$eV}. The anisotropy was measured by the fraction of arrival directions that are less than $3.1^\circ$ from the position of an active galactic nucleus within 75 Mpc (using the Véron-Cetty and Véron $12^{\rm th}$ catalog). An updated measurement of this fraction is reported here using the arrival directions of cosmic rays recorded above the same energy threshold through 31 December 2009. The number of arrival directions has increased from 27 to 69, allowing a more precise measurement. The correlating fraction is $(38^{+7}_{-6})%$, compared with $21%$ expected for isotropic cosmic rays. This is down from the early estimate of $(69^{+11}_{-13})%$. The enlarged set of arrival directions is examined also in relation to other populations of nearby extragalactic objects: galaxies in the 2 Microns All Sky Survey and active galactic nuclei detected in hard X-rays by the Swift Burst Alert Telescope. A celestial region around the position of the radiogalaxy Cen A has the largest excess of arrival directions relative to isotropic expectations. The 2-point autocorrelation function is shown for the enlarged set of arrival directions and compared to the isotropic expectation.

preprint2009arXiv

Lateral Distribution of the Radio Signal in Extensive Air Showers Measured with LOPES

The antenna array LOPES is set up at the location of the KASCADE-Grande extensive air shower experiment in Karlsruhe, Germany and aims to measure and investigate radio pulses from Extensive Air Showers. The coincident measurements allow us to reconstruct the electric field strength at observation level in dependence of general EAS parameters. In the present work, the lateral distribution of the radio signal in air showers is studied in detail. It is found that the lateral distributions of the electric field strengths in individual EAS can be described by an exponential function. For about 20% of the events a flattening towards the shower axis is observed, preferentially for showers with large inclination angle. The estimated scale parameters R0 describing the slope of the lateral profiles range between 100 and 200 m. No evidence for a direct correlation of R0 with shower parameters like azimuth angle, geomagnetic angle, or primary energy can be found. This indicates that the lateral profile is an intrinsic property of the radio emission during the shower development which makes the radio detection technique suitable for large scale applications.

preprint2009arXiv

Radiative hydrodynamics simulations of red supergiant stars: I. interpretation of interferometric observations

Context. It has been suggested that convection in Red Supergiant (RSG) stars gives rise to large-scale granules causing observable surface inhomogeneities. This convection is also extremely vigorous, and suspected to be one of the causes of mass-loss in RSGs. It must thus be understood in details. Evidence has been accumulated that there are asymmetries in the photospheres of RSGs, but detailedstudies of granulation are still lacking. Interferometric observations offer an exciting possibility to tackle this question, but they are still often interpreted using smooth symmetrical limb-darkened intensity distributions, or very simple spotted ad hoc models. Aims. We explore the impact of the granulation on visibility curves and closure phases using the radiative transfer code OPTIM3D. We simultaneously assess how 3D simulations of convection in RSG with CO5BOLD can be tested against these observations. Methods. We use 3D radiative-hydrodynamics (RHD) simulations of convection to compute intensity maps at various wavelengths and time, from which we derive interferometric visibility amplitudes and phases. We study their behaviour with time, position angle, and wavelength, and compare them to observations of the RSG alpha Ori Results. We provide average limb-darkening coefficients for RSGs. We detail the prospects for the detection and characterization of granulation (contrast, size) on RSGs. We demonstrate that our RHD simulations provide an excellent fit to existing interferometric observation of alpha Ori, contrary to limb darkened disks. This confirms the existence of large convective cells on the surface of Betelgeuse.

preprint2009arXiv

The Fluorescence Detector of the Pierre Auger Observatory

The Pierre Auger Observatory is a hybrid detector for ultra-high energy cosmic rays. It combines a surface array to measure secondary particles at ground level together with a fluorescence detector to measure the development of air showers in the atmosphere above the array. The fluorescence detector comprises 24 large telescopes specialized for measuring the nitrogen fluorescence caused by charged particles of cosmic ray air showers. In this paper we describe the components of the fluorescence detector including its optical system, the design of the camera, the electronics, and the systems for relative and absolute calibration. We also discuss the operation and the monitoring of the detector. Finally, we evaluate the detector performance and precision of shower reconstructions.

preprint2009arXiv

Upper limit on the cosmic-ray photon fraction at EeV energies from the Pierre Auger Observatory

From direct observations of the longitudinal development of ultra-high energy air showers performed with the Pierre Auger Observatory, upper limits of 3.8%, 2.4%, 3.5% and 11.7% (at 95% c.l.) are obtained on the fraction of cosmic-ray photons above 2, 3, 5 and 10 EeV (1 EeV = 10^18 eV) respectively. These are the first experimental limits on ultra-high energy photons at energies below 10 EeV. The results complement previous constraints on top-down models from array data and they reduce systematic uncertainties in the interpretation of shower data in terms of primary flux, nuclear composition and proton-air cross-section.

preprint2009arXiv

VLTI/AMBER spectro-interferometric imaging of VX Sgr's inhomogenous outer atmosphere

Aims. We aim to explore the photosphere of the very cool late-type star VX Sgr and in particular the existence and characterization of molecular layers above the continuum forming photosphere. Methods. We obtained interferometric observations with the VLTI/AMBER interferometer using the fringe tracker FINITO in the spectral domain 1.45-2.50 micron with a spectral resolution of about 35 and baselines ranging from 15 to 88 meters.We perform independent image reconstruction for different wavelength bins and fit the interferometric data with a geometrical toy model.We also compare the data to 1D dynamical models of Miras atmosphere and to 3D hydrodynamical simulations of red supergiant (RSG) and asymptotic giant branch (AGB) stars. Results. Reconstructed images and visibilities show a strong wavelength dependence. The H-band images display two bright spots whose positions are confirmed by the geometrical toy model. The inhomogeneities are qualitatively predicted by 3D simulations. At about 2,00 micron and in the region 2,35 - 2,50 micron, the photosphere appears extended and the radius is larger than in the H band. In this spectral region, the geometrical toy model locates a third bright spot outside the photosphere that can be a feature of the molecular layers. The wavelength dependence of the visibility can be qualitatively explained by 1D dynamical models of Mira atmospheres. The best-fitting photospheric models show a good match with the observed visibilities and give a photospheric diameter of theta = 8,82+-0,50 mas. The H2O molecule seems to be the dominant absorber in the molecular layers. Conclusions. We show that the atmosphere of VX Sgr rather resembles Mira/AGB star model atmospheres than RSG model atmospheres. In particular, we see molecular (water) layers that are typical for Mira stars.