Search

Observing planetary auroral radio emission is the most promising method to detect exoplanetary magnetic fields, the knowledge of which will provide valuable insights into the planet's interior structure, atmospheric escape, and habitability. We present LOFAR-LBA circularly polarized beamformed observations of the exoplanetary systems 55 Cancri, $\upsilon$ Andromedae, and $\tau$ Bo\"{o}tis. We tentatively detect circularly polarized bursty emission from the $\tau$ Bo\"{o}tis system in the range 14-21 MHz with a flux density of $\sim$890 mJy and with a significance of $\sim$3$\sigma$. For this detection, no signal is seen in the OFF-beams, and we do not find any potential causes which might cause false positives. We also tentatively detect slowly variable circularly polarized emission from $\tau$ Bo\"{o}tis in the range 21-30 MHz with a flux density of $\sim$400 mJy and with a statistical significance of $>$8$\sigma$. The slow emission is structured in the time-frequency plane and shows an excess in the ON-beam with respect to the two simultaneous OFF-beams. Close examination casts some doubts on the reality of the slowly varying signal. We discuss in detail all the arguments for and against an actual detection. Furthermore, a $\sim$2$\sigma$ marginal signal is found from the $\upsilon$ Andromedae system and no signal is detected from the 55 Cancri system. Assuming the detected signals are real, we discuss their potential origin. Their source probably is the $\tau$ Bootis planetary system, and a possible explanation is radio emission from the exoplanet $\tau$ Bootis b via the cyclotron maser mechanism. Assuming a planetary origin, we derived limits for the planetary polar surface magnetic field strength, finding values compatible with theoretical predictions. Further low-frequency observations are required to confirm this possible first detection of an exoplanetary radio signal. [Abridged], Comment: Accepted for publication in A&A (Oct. 22, 2020). 29 pages (15 pgs main text), 17 figures, 8 tables, 10 appendices

The discovery of ubiquitous habitable extrasolar planets, combined with revolutionary advances in instrumentation and observational capabilities, has ushered in a renaissance in the search for extra-terrestrial intelligence (SETI). Large scale SETI activities are now underway at numerous international facilities. The Five-hundred-meter Aperture Spherical radio Telescope (FAST) is the largest single-aperture radio telescope in the world, well positioned to conduct sensitive searches for radio emission indicative of exo-intelligence. SETI is one of the five key science goals specified in the original FAST project plan. A collaboration with the Breakthrough Listen Initiative has been initiated in 2016 with a joint statement signed both by Dr. Jun Yan, the then director of the National Astronomical Observatories, Chinese Academy of Sciences (NAOC), and Dr. Peter Worden, the Chairman of the Breakthrough Prize Foundation. In this paper, we highlight some of the unique features of FAST that will allow for novel SETI observations. We identify and describe three different signal types indicative of a technological source, namely, narrow-band, wide-band artificially dispersed, and modulated signals. We here propose observations with FAST to achieve sensitivities never before explored., Comment: 14 pages, 8 figures, Accepted by Research in Astron. Astrophys. (FAST special issue)

We report on a search for artificial narrowband signals of 20 stars within the restricted Earth Transit Zone as a part of the ten-year Breakthrough Listen (BL) search for extraterrestrial intelligence. The restricted Earth Transit Zone is the region of the sky from which an observer would see the Earth transit the Sun with an impact parameter of less than 0.5. This region of the sky is geometrically unique, providing a potential way for an extraterrestrial intelligence to discover the Solar System. The targets were nearby (7-143 pc) and the search covered an electromagnetic frequency range of 3.95-8.00 GHz. We used the Robert C. Byrd Green Bank Telescope to perform these observations with the standard BL data recorder. We searched these data for artificial narrowband ($\sim$Hz) signals with Doppler drift rates of $\pm 20$ Hz s$^{-1}$. We found one set of potential candidate signals on the target HIP 109656 which was then found to be consistent with known properties of anthropogenic radio frequency interference. We find no evidence for radio technosignatures from extraterrestrial intelligence in our observations. The observing campaign achieved a minimum detectable flux which would have allowed detections of emissions that were $10^{-3}$ to $0.88$ times as powerful as the signaling capability of the Arecibo radar transmitter, for the nearest and furthest stars respectively. We conclude that at least $8\%$ of the systems in the restricted Earth Transit Zone within 150 pc do not possess the type of transmitters searched in this survey. To our knowledge, this is the first targeted search for extraterrestrial intelligence of the restricted Earth Transit Zone. All data used in this paper are publicly available via the Breakthrough Listen Public Data Archive (http://seti.berkeley.edu/bldr2)., Comment: 17 pages, 8 figures, submitted to ApJ

In this paper we report the current status of a new research program. The primary goal of the "Vanishing & Appearing Sources during a Century of Observations" (VASCO) project is to search for vanishing and appearing sources using existing survey data to find examples of exceptional astrophysical transients. The implications of finding such objects extend from traditional astrophysics fields to the more exotic searches for evidence of technologically advanced civilizations. In this first paper we present new, deeper observations of the tentative candidate discovered by Villarroel et al. (2016). We then perform the first searches for vanishing objects throughout the sky by comparing 600 million objects from the US Naval Observatory Catalogue (USNO) B1.0 down to a limiting magnitude of $\sim 20 - 21$ with the recent Pan-STARRS Data Release-1 (DR1) with a limiting magnitude of $\sim$ 23.4. We find about 150,000 preliminary candidates that do not have any Pan-STARRS counterpart within a 30 arcsec radius. We show that these objects are redder and have larger proper motions than typical USNO objects. We visually examine the images for a subset of about 24,000 candidates, superseding the 2016 study with a sample ten times larger. We find about $\sim$ 100 point sources visible in only one epoch in the red band of the USNO which may be of interest in searches for strong M dwarf flares, high-redshift supernovae or other catagories of unidentified red transients., Comment: SETI meets time domain astronomy. Accepted into the Astronomical Journal

The star EPIC 249706694 (HD 139139) was found to exhibit 28 transit-like events over an 87 day period during the Kepler mission's K2 Campaign 15 (Rappaport et al. 2019). These events did not fall into an identifiable pattern and could not be explained by a multitude of transit scenarios explored by the authors. We conduct follow-up observations at C-band frequencies with the Green Bank Telescope as part of the ongoing Breakthrough Listen search for technosignatures. We search for narrow band signals above a signal-to-noise threshold of 10 and with Doppler drift rates within +-5 Hz/s. We detect no evidence of technosignatures from EPIC 249706694 and derive an upper limit for the EIRP (Equivalent Isotropic Radiated Power) of putative transmissions to be 10 TW., Comment: 3 pages, 1 table. Accepted for publication in RNAAS

A radio transmitter which is accelerating with a non-zero radial component with respect to a receiver will produce a signal that appears to change its frequency over time. This effect, commonly produced in astrophysical situations where orbital and rotational motions are ubiquitous, is called a drift rate. In radio SETI (Search for Extraterrestrial Intelligence) research, it is unknown a priori which frequency a signal is being sent at, or even if there will be any drift rate at all besides motions within the solar system. Therefore a range of potential drift rates need to be individually searched, and a maximum drift rate needs to be chosen. The middle of this range is zero, indicating no acceleration, but the absolute value for the limits remains unconstrained. A balance must be struck between computational time and the possibility of excluding a signal from an ETI. In this work, we examine physical considerations that constrain a maximum drift rate and highlight the importance of this problem in any narrowband SETI search. We determine that a normalized drift rate of 200 nHz (eg. 200 Hz/s at 1 GHz) is a generous, physically motivated guideline for the maximum drift rate that should be applied to future narrowband SETI projects if computational capabilities permit., Comment: Published in the Astrophysical Journal, 16 pages, 2 figures

The discovery of the ubiquity of habitable extrasolar planets, combined with revolutionary advances in instrumentation and observational capabilities, have ushered in a renaissance in the millenia-old quest to answer our most profound question about the Universe and our place within it - Are we alone? The Breakthrough Listen Initiative, announced in July 2015 as a 10-year 100M USD program, is the most comprehensive effort in history to quantify the distribution of advanced, technologically capable life in the universe. In this white paper, we outline the status of the on-going observing campaign with our primary observing facilities, as well as planned activities with these instruments over the next few years. We also list collaborative facilities which will conduct searches for technosignatures in either primary observing mode, or commensally. We highlight some of the novel analysis techniques we are bringing to bear on multi-petabyte data sets, including machine learning tools we are deploying to search for a broader range of technosignatures than was previously possible., Comment: Astronomy and Astrophysics Decadal Survey 2020, APC white paper, To be published in BAAS

Breakthrough Listen (BL) is a ten-year initiative to search for signatures of technologically capable life beyond Earth via radio and optical observations of the local Universe. A core part of the BL program is a comprehensive survey of 1702 nearby stars at radio wavelengths (1-10 GHz). Here, we report on observations with the 64-m CSIRO Parkes radio telescope in New South Wales, Australia, and the 100-m Robert C. Byrd Green Bank radio telescope in West Virginia, USA. Over 2016 January to 2019 March, a sample of 1138 stars was observed at Green Bank using the 1.10-1.90 GHz and 1.80-2.80 GHz receivers, and 189 stars were observed with Parkes over 2.60-3.45 GHz. We searched these data for the presence of engineered signals with Doppler-acceleration drift rates between -4 to 4 Hz/s. Here, we detail our data analysis techniques and provide examples of detected events. After excluding events with characteristics consistent with terrestrial radio interference, we are left with zero candidates. Given the sensitivity of our observations, we can put an upper limit on the power of potential radio transmitters at these frequencies at 2x10^12 W, and 9x10^12 W for GBT and Parkes respectively. These observations constitute the most comprehensive search over 1.10-3.45 GHz for technosignatures to date. All data products, totalling ~219 TB, are available for download as part of the first BL data release (DR1), as described in a companion paper (Lebofsky et. al., 2019), Comment: Published 2020-02-05 (AJ)

Breakthrough Listen is the most comprehensive and sensitive search for extraterrestrial intelligence (SETI) to date, employing a collection of international observational facilities including both radio and optical telescopes. During the first three years of the Listen program, thousands of targets have been observed with the Green Bank Telescope (GBT), Parkes Telescope and Automated Planet Finder. At GBT and Parkes, observations have been performed ranging from 700 MHz to 26 GHz, with raw data volumes averaging over 1PB / day. A pseudo-real time software spectroscopy suite is used to produce multi-resolution spectrograms amounting to approximately 400 GB hr^-1 GHz^-1 beam^-1. For certain targets, raw baseband voltage data is also preserved. Observations with the Automated Planet Finder produce both 2-dimensional and 1-dimensional high resolution (R~10^5) echelle spectral data. Although the primary purpose of Listen data acquisition is for SETI, a range of secondary science has also been performed with these data, including studies of fast radio bursts. Other current and potential research topics include spectral line studies, searches for certain kinds of dark matter, probes of interstellar scattering, pulsar searches, radio transient searches and investigations of stellar activity. Listen data are also being used in the development of algorithms, including machine learning approaches to modulation scheme classification and outlier detection, that have wide applicability not just for astronomical research but for a broad range of science and engineering. In this paper, we describe the hardware and software pipeline used for collection, reduction, archival, and public dissemination of Listen data. We describe the data formats and tools, and present Breakthrough Listen Data Release 1.0 (BLDR 1.0), a defined set of publicly-available raw and reduced data totalling 1 PB., Comment: 28 pages, 15 figures, PASP accepted

Here, we report on the detection and verification of Fast Radio Burst FRB 180301, which occurred on UTC 2018 March 1 during the Breakthrough Listen observations with the Parkes telescope. Full-polarization voltage data of the detection were captured--a first for non-repeating FRBs--allowing for coherent de-dispersion and additional verification tests. The coherently de-dispersed dynamic spectrum of FRB 180301 shows complex, polarized frequency structure over a small fractional bandwidth. As FRB 180301 was detected close to the geosynchronous satellite band during a time of known 1-2 GHz satellite transmissions, we consider whether the burst was due to radio interference emitted or reflected from an orbiting object. Based on the preponderance of our verification tests, we find that FRB 180301 is likely of astrophysical origin, but caution that anthropogenic sources cannot conclusively be ruled out., Comment: Submitted to MNRAS

Boyajian' s Star (KIC 8462852) has received significant attention due to its unusual periodic brightness fluctuations detected by the Kepler Spacecraft and subsequent ground based observations. Possible explanations for the dips in the photometric measurements include interstellar or circumstellar dust, and it has been speculated that an artificial megastructure could be responsible. We analyze 177 high-resolution spectra of Boyajian's Star in an effort to detect potential laser signals from extraterrestrial civilizations. The spectra were obtained by the Lick Observatory's Automated Planet Finder telescope as part of the Breakthrough Listen Project, and cover the wavelength range of visible light from 374 to 970 nm. We calculate that the APF would be capable of detecting lasers of power greater than approximately 24 MW at the distance of Boyajian's Star, d = 1470 ly. The top candidates from the analysis can all be explained as either cosmic ray hits, stellar emission lines or atmospheric air glow emission lines.

We report on high-resolution spectra obtained by the Automated Planet Finder and high resolution optical Levy Spectrometer and the search for periodic spectral modulations, such as those reported in Borra (2016). In their analysis of 2.5 million spectra from the Sloan Digital Sky Survey, Borra (2016) report periodic spectral modulations in 234 stars, and suggest that these signals may be evidence of extra-terrestrial civilizations. To further evaluate this claim, we observed a total of three of the 234 stars with the Automated Planet Finder Telescope and Levy Spectrometer including all stars brighter than a visual magnitude of 14. Fourier analysis of the resultant spectra of these three sources does not reveal any periodic spectral modulations at the reported period, nor at any other period., Comment: 8 pages, 5 figures

Breakthrough Listen is a ten-year initiative to search for signatures of technologies created by extraterrestrial civilizations at radio and optical wavelengths. Here, we detail the digital data recording system deployed for Breakthrough Listen observations at the 64-m aperture CSIRO Parkes Telescope in New South Wales, Australia. The recording system currently implements two recording modes: a dual-polarization, 1.125 GHz bandwidth mode for single beam observations, and a 26-input, 308-MHz bandwidth mode for the 21-cm multibeam receiver. The system is also designed to support a 3 GHz single-beam mode for the forthcoming Parkes ultra-wideband feed. In this paper, we present details of the system architecture, provide an overview of hardware and software, and present initial performance results., Comment: v2, accepted to PASA

We report the first detections of the repeating fast radio burst source FRB 121102 above 5.2 GHz. Observations were performed using the 4$-$8 GHz receiver of the Robert C. Byrd Green Bank Telescope with the Breakthrough Listen digital backend. We present the spectral, temporal and polarization properties of 21 bursts detected within the first 60 minutes of a total 6-hour observations. These observations comprise the highest burst density yet reported in the literature, with 18 bursts being detected in the first 30 minutes. A few bursts clearly show temporal sub-structures with distinct spectral properties. These sub-structures superimpose to provide enhanced peak signal-to-noise ratio at higher trial dispersion measures. Broad features occur in $\sim 1$ GHz wide subbands that typically differ in peak frequency between bursts within the band. Finer-scale structures ($\sim 10-50$ MHz) within these bursts are consistent with that expected from Galactic diffractive interstellar scintillation. The bursts exhibit nearly 100% linear polarization, and a large average rotation measure of 9.359$\pm$0.012 $\times$ 10$^{\rm 4}$ rad m$^{\rm -2}$ (in the observer's frame). No circular polarization was found for any burst. We measure an approximately constant polarization position angle in the 13 brightest bursts. The peak flux densities of the reported bursts have average values (0.2$\pm$0.1 Jy), similar to those seen at lower frequencies ($<3$ GHz), while the average burst widths (0.64$\pm$0.46 ms) are relatively narrower., Comment: Accepted for publication in ApJ

We have conducted a search for radio emission consistent with an artificial source targeting 1I/'Oumuamua with the Robert C. Byrd Green Bank Telescope (GBT) between 1.1 and 11.6 GHz. We searched the data for narrowband signals and found none. Given the close proximity to this interstellar object, we can place limits to putative transmitters with extremely low power (0.08 W)., Comment: 1 table, Submitted to RN of AAS

We undertook observations with the Green Bank Telescope, simultaneously with the 300m telescope in Arecibo, as a follow-up of a possible flare of radio emission from Ross 128. We report here the non-detections from the GBT observations in C band (4-8 GHz), as well as non-detections in archival data at L band (1.1-1.9 GHz). We suggest that a likely scenario is that the emission comes from one or more satellites passing through the same region of the sky., Comment: 6 pages, 2 figures. Accepted to International Journal of Astrobiology

We report on a search for engineered signals from a sample of 692 nearby stars using the Robert C. Byrd Green Bank Telescope, undertaken as part of the $Breakthrough~Listen~Initiative$ search for extraterrestrial intelligence. Observations were made over 1.1$-$1.9 GHz (L band), with three sets of five-minute observations of the 692 primary targets, interspersed with five-minute observations of secondary targets. By comparing the "ON" and "OFF" observations we are able to identify terrestrial interference and place limits on the presence of engineered signals from putative extraterrestrial civilizations inhabiting the environs of the target stars. During the analysis, eleven events passed our thresholding algorithm, but a detailed analysis of their properties indicates they are consistent with known examples of anthropogenic radio frequency interference. We conclude that, at the time of our observations, none of the observed systems host high-duty-cycle radio transmitters emitting between 1.1 and 1.9 GHz with an Equivalent Isotropic Radiated Power of $\sim10^{13}$ W, which is readily achievable by our own civilization. Our results suggest that fewer than $\sim$ 0.1$\%$ of the stellar systems within 50 pc possess the type of transmitters searched in this survey., Comment: ApJ. 13 pages, 7 figures, 4 tables. This version has a typo correction to Eq. (4)

The Breakthrough Listen Initiative is undertaking a comprehensive search for radio and optical signatures from extraterrestrial civilizations. An integral component of the project is the design and implementation of wide-bandwidth data recorder and signal processing systems. The capabilities of these systems, particularly at radio frequencies, directly determine survey speed; further, given a fixed observing time and spectral coverage, they determine sensitivity as well. Here, we detail the Breakthrough Listen wide-bandwidth data recording system deployed at the 100-m aperture Robert C. Byrd Green Bank Telescope. The system digitizes up to 6 GHz of bandwidth at 8 bits for both polarizations, storing the resultant 24 GB/s of data to disk. This system is among the highest data rate baseband recording systems in use in radio astronomy. A future system expansion will double recording capacity, to achieve a total Nyquist bandwidth of 12 GHz in two polarizations. In this paper, we present details of the system architecture, along with salient configuration and disk-write optimizations used to achieve high-throughput data capture on commodity compute servers and consumer-class hard disk drives.

We observe a large fraction of circular polarization in radio emission from extensive air showers recorded during thunderstorms, much higher than in the emission from air showers measured during fair-weather circumstances. We show that the circular polarization of the air showers measured during thunderstorms can be explained by the change in the direction of the transverse current as a function of altitude induced by atmospheric electric fields. Thus by using the full set of Stokes parameters for these events, we obtain a good characterization of the electric fields in thunderclouds. We also measure a large horizontal component of the electric fields in the two events that we have analysed.

We present the target selection for the Breakthrough Listen search for extraterrestrial intelligence during the first year of observations at the Green Bank Telescope, Parkes Telescope and Automated Planet Finder. On the way to observing 1,000,000 nearby stars in search of technological signals, we present three main sets of objects we plan to observe in addition to a smaller sample of exotica. We choose the 60 nearest stars, all within 5.1 pc from the sun. Such nearby stars offer the potential to observe faint radio signals from transmitters having a power similar to those on Earth. We add a list of 1649 stars drawn from the Hipparcos catalog that span the Hertzprung-Russell diagram, including all spectral types along the main sequence, subgiants, and giant stars. This sample offers diversity and inclusion of all stellar types, but with thoughtful limits and due attention to main sequence stars. Our targets also include 123 nearby galaxies composed of a "morphological-type-complete" sample of the nearest spirals, ellipticals, dwarf spherioidals, and irregulars. While their great distances hamper the detection of technological electromagnetic radiation, galaxies offer the opportunity to observe billions of stars simultaneously and to sample the bright end of the technological luminosity function. We will also use the Green Bank and Parkes telescopes to survey the plane and central bulge of the Milky Way. Finally, the complete target list includes several classes of exotica, including white dwarfs, brown dwarfs, black holes, neutron stars, and asteroids in our Solar System.

The low flux of the ultra-high energy cosmic rays (UHECR) at the highest energies provides a challenge to answer the long standing question about their origin and nature. Even lower fluxes of neutrinos with energies above $10^{22}$ eV are predicted in certain Grand-Unifying-Theories (GUTs) and e.g.\ models for super-heavy dark matter (SHDM). The significant increase in detector volume required to detect these particles can be achieved by searching for the nano-second radio pulses that are emitted when a particle interacts in Earth's moon with current and future radio telescopes. In this contribution we present the design of an online analysis and trigger pipeline for the detection of nano-second pulses with the LOFAR radio telescope. The most important steps of the processing pipeline are digital focusing of the antennas towards the Moon, correction of the signal for ionospheric dispersion, and synthesis of the time-domain signal from the polyphased-filtered signal in frequency domain. The implementation of the pipeline on a GPU/CPU cluster will be discussed together with the computing performance of the prototype., Comment: Proceedings of the 22nd International Conference on Computing in High Energy and Nuclear Physics (CHEP2016), USA

The low flux of the ultra-high energy cosmic rays (UHECR) at the highest energies provides a challenge to answer the long standing question about their origin and nature. A significant increase in the number of detected UHECR is expected to be achieved by employing Earth's moon as detector, and search for short radio pulses that are emitted when a particle interacts in the lunar rock. Observation of these short pulses with current and future radio telescopes also allows to search for the even lower fluxes of neutrinos with energies above $10^{22}$ eV, that are predicted in certain Grand-Unifying-Theories (GUTs), and e.g. models for super-heavy dark matter (SHDM). In this contribution we present the initial design for such a search with the LOFAR radio telescope., Comment: To be published in the Proceedings of the ARENA2016 conference, Groningen, The Netherlands

Photometric lightcurves were also obtained for 14 of them during the Spitzer observations to provide the context of the observations and reliable estimates of their absolute magnitudes. The extracted mid-IR spectra were analyzed using a modified standard thermal model (STM) and a thermophysical model (TPM) that takes into account the shape and geometry of the large primary at the time of the Spitzer observation. We derived a reliable estimate of the size, albedo, and beaming factor for each of these asteroids, representing three main taxonomic groups: C, S, and X. For large (volume-equivalent system diameter Deq $\lt$ 130 km) binary asteroids, the TPM analysis indicates a low thermal inertia ($\Gamma$ < $\sim$100 J s-1/2K-1m-2) and their emissivity spectra display strong mineral features, implying that they are covered with a thick layer of thermally insulating regolith. The smaller (surface-equivalent system diameter Deff $\lt$17 km) asteroids also show some emission lines of minerals, but they are significantly weaker, consistent with regoliths with coarser grains, than those of the large binary asteroids. The average bulk densities of these multiple asteroids vary from 0.7-1.7 g/cm3 (P-, C- type) to $\sim$2 g/cm3 (S-type). The highest density is estimated for the M-type (22) Kalliope (3.2 $\pm$ 0.9 g/cm3). The spectral energy distributions (SED) and emissivity spectra, made available as a supplement document, could help to constrain the surface compositions of these asteroids., Comment: 69 pages, 5 Figures, 8 Tables

We describe a method for spectral cleaning and timing calibration of short voltage time series data from individual radio interferometer receivers. It makes use of the phase differences in Fast Fourier Transform (FFT) spectra across antenna pairs. For strong, localized terrestrial sources these are stable over time, while being approximately uniform-random for a sum over many sources or for noise. Using only milliseconds-long datasets, the method finds the strongest interfering transmitters, a first-order solution for relative timing calibrations, and faulty data channels. No knowledge of gain response or quiescent noise levels of the receivers is required. With relatively small data volumes, this approach is suitable for use in an online system monitoring setup for interferometric arrays. We have applied the method to our cosmic-ray data collection, a collection of measurements of short pulses from extensive air showers, recorded by the LOFAR radio telescope. Per air shower, we have collected 2 ms of raw time series data for each receiver. The spectral cleaning has a calculated optimal sensitivity corresponding to a power signal-to-noise ratio of 0.08 (or -11 dB) in a spectral window of 25 kHz, for 2 ms of data in 48 antennas. This is well sufficient for our application. Timing calibration across individual antenna pairs has been performed at 0.4 ns precision; for calibration of signal clocks across stations of 48 antennas the precision is 0.1 ns. Monitoring differences in timing calibration per antenna pair over the course of the period 2011 to 2015 shows a precision of 0.08 ns, which is useful for monitoring and correcting drifts in signal path synchronizations. A cross-check method for timing calibration is presented, using a pulse transmitter carried by a drone flying over the array. Timing precision is similar, 0.3 ns., Comment: Accepted for publication in Astronomy & Astrophysics; 10 pages, 7 figures

Cosmic rays are the highest energy particles found in nature. Measurements of the mass composition of cosmic rays between 10^{17} eV and 10^{18} eV are essential to understand whether this energy range is dominated by Galactic or extragalactic sources. It has also been proposed that the astrophysical neutrino signal comes from accelerators capable of producing cosmic rays of these energies. Cosmic rays initiate cascades of secondary particles (air showers) in the atmosphere and their masses are inferred from measurements of the atmospheric depth of the shower maximum, Xmax, or the composition of shower particles reaching the ground. Current measurements suffer from either low precision, or a low duty cycle and a high energy threshold. Radio detection of cosmic rays is a rapidly developing technique, suitable for determination of Xmax with a duty cycle of in principle nearly 100%. The radiation is generated by the separation of relativistic charged particles in the geomagnetic field and a negative charge excess in the shower front. Here we report radio measurements of Xmax with a mean precision of 16 g/cm^2 between 10^{17}-10^{17.5} eV. Because of the high resolution in $Xmax we can determine the mass spectrum and find a mixed composition, containing a light mass fraction of ~80%. Unless the extragalactic component becomes significant already below 10^{17.5} eV, our measurements indicate an additional Galactic component dominating at this energy range., Comment: 35 pages, 11 figures, updated version: Pierre Auger Observatory data ICRC 2015 added to Fig 2

The atmospheric electric fields in thunderclouds have been shown to significantly modify the intensity and polarization patterns of the radio footprint of cosmic-ray-induced extensive air showers. Simulations indicated a very non-linear dependence of the signal strength in the frequency window of 30-80 MHz on the magnitude of the atmospheric electric field. In this work we present an explanation of this dependence based on Monte-Carlo simulations, supported by arguments based on electron dynamics in air showers and expressed in terms of a simplified model. We show that by extending the frequency window to lower frequencies additional sensitivity to the atmospheric electric field is obtained.

Air showers induced by cosmic rays create nanosecond pulses detectable at radio frequencies. These pulses have been measured successfully in the past few years at the LOw Frequency ARray (LOFAR) and are used to study the properties of cosmic rays. For a complete understanding of this phenomenon and the underlying physical processes, an absolute calibration of the detecting antenna system is needed. We present three approaches that were used to check and improve the antenna model of LOFAR and to provide an absolute calibration of the whole system for air shower measurements. Two methods are based on calibrated reference sources and one on a calibration approach using the diffuse radio emission of the Galaxy, optimized for short data-sets. An accuracy of 19% in amplitude is reached. The absolute calibration is also compared to predictions from air shower simulations. These results are used to set an absolute energy scale for air shower measurements and can be used as a basis for an absolute scale for the measurement of astronomical transients with LOFAR., Comment: 34 pages, 10 figures

The energy reconstruction of extensive air showers measured with the LOFAR Radboud Air Shower Array (LORA) is presented in detail. LORA is a particle detector array located in the center of the LOFAR radio telescope in the Netherlands. The aim of this work is to provide an accurate and independent energy measurement for the air showers measured through their radio signal with the LOFAR antennas. The energy reconstruction is performed using a parameterized relation between the measured shower size and the cosmic-ray energy obtained from air shower simulations. In order to illustrate the capabilities of LORA, the all-particle cosmic-ray energy spectrum has been reconstructed, assuming that cosmic rays are composed only of protons or iron nuclei in the energy range between $\sim2\times10^{16}$ and $2\times10^{18}$ eV. The results are compatible with literature values and a changing mass composition in the transition region from a galactic to an extragalactic origin of cosmic rays., Comment: 12 pages, Accepted by Astroparticle Physics, updated authors list

We present measurements of radio emission from cosmic ray air showers that took place during thunderstorms. The intensity and polarization patterns of these air showers are radically different from those measured during fair-weather conditions. With the use of a simple two-layer model for the atmospheric electric field, these patterns can be well reproduced by state-of-the-art simulation codes. This in turn provides a novel way to study atmospheric electric fields., Comment: 6 pages, 3 figures, accepted for publication in Physical Review Letters

The pattern of the radio emission of air showers is finely sampled with the Low-Frequency ARray (LOFAR). A set of 382 measured air showers is used to test a fast, analytic parameterization of the distribution of pulse powers. Using this parameterization we are able to reconstruct the shower axis and give estimators for the energy of the air shower as well as the distance to the shower maximum., Comment: 15 pages, 10 figures, accepted for publication in JCAP

Measuring radio emission from air showers offers a novel way to determine properties of the primary cosmic rays such as their mass and energy. Theory predicts that relativistic time compression effects lead to a ring of amplified emission which starts to dominate the emission pattern for frequencies above ~100 MHz. In this article we present the first detailed measurements of this structure. Ring structures in the radio emission of air showers are measured with the LOFAR radio telescope in the frequency range of 110 - 190 MHz. These data are well described by CoREAS simulations. They clearly confirm the importance of including the index of refraction of air as a function of height. Furthermore, the presence of the Cherenkov ring offers the possibility for a geometrical measurement of the depth of shower maximum, which in turn depends on the mass of the primary particle., Comment: 20 pages, 10 figures, accpeted for publication in Astroparticle Physics

We provide a taxonomic and compositional characterization of Multiple Asteroid Systems (MASs) located in the main belt (MB) using visible and near-infrared (0.45-2.5 um) spectral data of 42 MB MASs. The mineralogical analysis is applied to determine meteorite analogs for the MASs, which, in turn, are applied to the MAS density measurements of Marchis et al. (2012) to estimate the system porosity. The macroporosities are used to evaluate the primary MAS formation hypotheses. The visible observing campaign includes 25 MASs obtained using the SOAR telescope with the Goodman High Throughput Spectrometer. The infrared observing campaign includes 34 MASs obtained using the NASA IRTF with the SpeX spectragraph. The MASs are classified using the Bus-DeMeo taxonomic system. We perform a NIR spectral band parameter analysis using a new analysis routine, the Spectral Analysis Routine for Asteroids (SARA). The SARA routine determines band centers, areas, and depths by utilizing the diagnostic absorption features near 1- and 2-um. The band parameter analysis provides the Gaffey subtype for the S-type MASs; the relative abundance olivine-to-pyroxene ratio; and olivine and pyroxene modal abundances for S-complex and V-type MASs. This mineralogical information is applied to determine meteorite analogs. We determine the H, L, and LL meteorite analogs for 15 MASs with ordinary chondrite-like (OC) mineralogies. We observe an excess (10/15) of LL-like mineralogies. Of the MASs with LL-like mineralogies, seven are consistent with Flora family membership, supporting the hypothesis that the Flora family is a source of LL-like NEAs and LL chondrites on Earth. Using the measured densities of the meteorite analog and the MAS densities from Marchis et al. (2012), we estimate the macroporosity for 13 MASs and find that all estimated macroporosities are in agreement with formation hypotheses., Comment: 45 pages, 8 figures Accepted to Icarus

The mass composition of cosmic rays contains important clues about their origin. Accurate measurements are needed to resolve long-standing issues such as the transition from Galactic to extragalactic origin, and the nature of the cutoff observed at the highest energies. Composition can be studied by measuring the atmospheric depth of the shower maximum Xmax of air showers generated by high-energy cosmic rays hitting the Earth's atmosphere. We present a new method to reconstruct Xmax based on radio measurements. The radio emission mechanism of air showers is a complex process that creates an asymmetric intensity pattern on the ground. The shape of this pattern strongly depends on the longitudinal development of the shower. We reconstruct Xmax by fitting two-dimensional intensity profiles, simulated with CoREAS, to data from the LOFAR radio telescope. In the dense LOFAR core, air showers are detected by hundreds of antennas simultaneously. The simulations fit the data very well, indicating that the radiation mechanism is now well-understood. The typical uncertainty on the reconstruction of Xmax for LOFAR showers is 17 g/cm^2., Comment: 12 pages, 10 figures, submitted to Phys. Rev. D

The measurement of the radio emission from extensive air showers, induced by high-energy cosmic rays is one of the key science projects of the LOFAR radio telescope. The LOfar Radboud air shower Array (LORA) has been installed in the core of LOFAR in the Netherlands. The main purpose of LORA is to measure the properties of air showers and to trigger the read-out of the LOFAR radio antennas to register extensive air showers. The experimental set-up of the array of scintillation detectors and its performance are described., Comment: 10 pages, Accepted for publication in Nuclear Instruments and Methods A

We present LOFAR measurements of radio emission from extensive air showers. We find that this emission is strongly polarized, with a median degree of polarization of nearly $99\%$, and that the angle between the polarization direction of the electric field and the Lorentz force acting on the particles, depends on the observer location in the shower plane. This can be understood as a superposition of the radially polarized charge-excess emission mechanism, first proposed by Askaryan and the geomagnetic emission mechanism proposed by Kahn and Lerche. We calculate the relative strengths of both contributions, as quantified by the charge-excess fraction, for $163$ individual air showers. We find that the measured charge-excess fraction is higher for air showers arriving from closer to the zenith. Furthermore, the measured charge-excess fraction also increases with increasing observer distance from the air shower symmetry axis. The measured values range from $(3.3\pm 1.0)\%$ for very inclined air showers at $25\, \mathrm{m}$ to $(20.3\pm 1.3)\%$ for almost vertical showers at $225\, \mathrm{m}$. Both dependencies are in qualitative agreement with theoretical predictions., Comment: 22 pages, 14 figures, accepted for publication in JCAP

Extensive air showers, induced by high energy cosmic rays impinging on the Earth's atmosphere, produce radio emission that is measured with the LOFAR radio telescope. As the emission comes from a finite distance of a few kilometers, the incident wavefront is non-planar. A spherical, conical or hyperbolic shape of the wavefront has been proposed, but measurements of individual air showers have been inconclusive so far. For a selected high-quality sample of 161 measured extensive air showers, we have reconstructed the wavefront by measuring pulse arrival times to sub-nanosecond precision in 200 to 350 individual antennas. For each measured air shower, we have fitted a conical, spherical, and hyperboloid shape to the arrival times. The fit quality and a likelihood analysis show that a hyperboloid is the best parametrization. Using a non-planar wavefront shape gives an improved angular resolution, when reconstructing the shower arrival direction. Furthermore, a dependence of the wavefront shape on the shower geometry can be seen. This suggests that it will be possible to use a wavefront shape analysis to get an additional handle on the atmospheric depth of the shower maximum, which is sensitive to the mass of the primary particle., Comment: Accepted for publication in Astroparticle Physics