NuRadioMC is a Monte Carlo framework designed to simulate ultra-high energy neutrino detectors that rely on the radio detection method. This method exploits the radio emission generated in the electromagnetic component of a particle shower following a neutrino interaction. NuRadioMC simulates everything from the neutrino interaction in a medium, the subsequent Askaryan radio emission, the propagation of the radio signal to the detector and finally the detector response. NuRadioMC is designed as a modern, modular Python-based framework, combining flexibility in detector design with user-friendliness. It includes a state-of-the-art event generator, an improved modelling of the radio emission, a revisited approach to signal propagation and increased flexibility and precision in the detector simulation. This paper focuses on the implemented physics processes and their implications for detector design. A variety of models and parameterizations for the radio emission of neutrino-induced showers are compared and reviewed. Comprehensive examples are used to discuss the capabilities of the code and different aspects of instrumental design decisions.

The Askaryan Radio Array is an ultra-high energy (UHE, $>10^{17}$\,eV) neutrino detector designed to observe neutrinos by searching for the radio waves emitted by the relativistic byproducts of neutrino-nucleon interactions in Antarctic ice. In this paper, we present constraints on the diffuse flux of ultra-high energy neutrinos between $10^{16}-10^{21}$\,eV resulting from a search for neutrinos in two complementary analyses, both analyzing four years of data (2013-2016) from the two deep stations (A2, A3) operating at that time. We place a 90\% CL upper limit on the diffuse neutrino flux at $10^{18}$ eV of $EF(E)=3.9\times10^{-16}$~$\textrm{cm}^2$$\textrm{s}^{-1}$$\textrm{sr}^{-1}$. The analysis leverages more than quadruple the livetime of the previous ARA limit, and at this energy is the most sensitive reported by an in-ice radio neutrino detector by more than a factor of four. Looking forward, ARA has yet to analyze more than half of its archived data, and expects that dataset to double in size in the next three years as livetime accumulates in the full five-station array. With this additional livetime, we project ARA will have world-leading sensitivity to neutrinos above $10^{17}$ eV by 2022.

The past decade has welcomed the emergence of cosmic neutrinos as a new messenger to explore the most extreme environments of the universe. The discovery measurement of cosmic neutrinos, announced by IceCube in 2013, has opened a new window of observation that has already resulted in new fundamental information that holds the potential to answer key questions associated with the high-energy universe, including: what are the sources in the PeV sky and how do they drive particle acceleration; where are cosmic rays of extreme energies produced, and on which paths do they propagate through the universe; and are there signatures of new physics at TeV-PeV energies and above? The planned advancements in neutrino telescope arrays in the next decade, in conjunction with continued progress in broad multimessenger astrophysics, promise to elevate the cosmic neutrino field from the discovery to the precision era and to a survey of the sources in the neutrino sky. The planned detector upgrades to the IceCube Neutrino Observatory, culminating in IceCube-Gen2 (an envisaged $400M facility with anticipated operation in the next decade, described in this white paper) are the cornerstone that will drive the evolution of neutrino astrophysics measurements.

We report the observation of radar echoes from the ionization trails of high-energy particle cascades. These data were taken at the SLAC National Accelerator Laboratory, where the full electron beam ($\sim$10$^9$ e$^-$ at $\sim$10 GeV/e$^-$) was directed into a plastic target to simulate an ultra high-energy neutrino interaction. This target was interrogated with radio waves, and coherent radio reflections from the cascades were detected, with properties consistent with theoretical expectations. This is the first definitive observation of radar echoes from high-energy particle cascades, which may lead to a viable neutrino detection technology for energies $\gtrsim 10^{16}$ eV.

We report on analysis of englacial radio-frequency (RF) pulser data received over horizontal baselines of 1--5 km, based on broadcasts from two sets of transmitters deployed to depths of up to 1500 meters at the South Pole. First, we analyze data collected usingtwo RF bicone transmitters 1400 meters below the ice surface, and frozen into boreholes drilled for the IceCube experiment in 2011. Additionally, in Dec., 2018, a fat-dipole antenna, fed by one of three high-voltage (~1 kV), fast (~(1-5 ns)) signal generators was lowered into the 1700-m deep icehole drilled for the South Pole Ice Core Experiment (SPICE), approximately 3 km from the geographic South Pole. Signals from transmitters were recorded on the five englacial multi-receiver ARA stations, with receiver depths between 60--200 m. We confirm the long, >1 km RF electric field attenuation length, test our observed signal arrival timing distributions against models, and measure birefringent asymmetries at the 0.15% level.

Testbeam experiment 576 (T576) at the SLAC National Accelerator Laboratory sought to make the first measurement of coherent radio reflections from the ionization produced in the wake of a high-energy particle shower. The >10 GeV electron beam at SLAC End Station A was directed into a large high-density polyethylene target to produce a shower analogous to that produced by an EeV neutrino interaction in ice. Continuous wave (CW) radio was transmitted into the target, and receiving antennas monitored for reflection of the transmitted signal from the ionization left in the wake of the shower. We detail the experiment and report on preliminary hints of a signal consistent with a radio reflection at a significance of 2.36$\sigma$. We recommend another test-beam measurement in order to verify the signal.

Background : The emergence of hyperon degrees of freedom in neutron star matter has been associated to first order phase transitions in some phenomenological models, but conclusions on the possible physical existence of an instability in the strangeness sector are strongly model dependent. Purpose : The purpose of the present study is to assess whether strangeness instabilities are related to specific values of the largely unconstrained hyperon interactions, and to study the effect of the strange meson couplings on phenomenological properties of neutron stars and supernova matter, once these latter are fixed to fulfill the constraints imposed by hypernuclear data. Method : We consider a phenomenological RMF model sufficiently simple to allow a complete exploration of the parameter space. Results : We show that no instability at supersaturation density exists for the RMF model, as long as the parameter space is constrained by basic physical requirements. This is at variance with a non-relativistic functional, with a functional behavior fitted through ab-initio calculations. Once the study is extended to include the full octet, we show that the parameter space allows reasonable radii for canonical neutron stars as well as massive stars above two-solar mass, together with an important strangeness content of the order of 30\%, slightly decreasing with increasing entropy, even in the absence of a strangeness driven phase transition. Conclusions : We conclude that the hyperon content of neutron stars and supernova matter cannot be established with present constraints, and is essentially governed by the unconstrained coupling to the strange isoscalar meson.

It has been argued that the small perturbations to the homogeneous and isotropic configurations of a canonical scalar field in an expanding universe do not grow. We show that this is not true in general, and clarify the root of the misunderstanding. We revisit a simple model in which the zero-mode of a free scalar field oscillates with high frequency around the minimum of the potential. Under this assumption the linear perturbations grow like those in the standard cold dark matter scenario, but with a Jeans length at the scale of the Compton wavelength of the scalar particle. Contrary to previous analyses in the literature our results do not rely on time-averages and/or fluid identifications, and instead we solve both analytically (in terms of a well-defined series expansion) and numerically the linearized Einstein-Klein-Gordon system. Also, we use gauge-invariant fields, which makes the physical analysis more transparent and simplifies the comparison with previous works carried out in different gauges. As a byproduct of this study we identify a time-dependent modulation of the different physical quantities associated to the background as well as the perturbations with potential observational consequences in dark matter models.

The discovery of almost 2000 exoplanets has revealed an unexpectedly diverse planet population. Observations to date have shown that our Solar System is certainly not representative of the general population of planets in our Milky Way. The key science questions that urgently need addressing are therefore: What are exoplanets made of? Why are planets as they are? What causes the exceptional diversity observed as compared to the Solar System? EChO (Exoplanet Characterisation Observatory) has been designed as a dedicated survey mission for transit and eclipse spectroscopy capable of observing a large and diverse planet sample within its four-year mission lifetime. EChO can target the atmospheres of super-Earths, Neptune-like, and Jupiter-like planets, in the very hot to temperate zones (planet temperatures of 300K-3000K) of F to M-type host stars. Over the next ten years, several new ground- and space-based transit surveys will come on-line (e.g. NGTS, CHEOPS, TESS, PLATO), which will specifically focus on finding bright, nearby systems. The current rapid rate of discovery would allow the target list to be further optimised in the years prior to EChO's launch and enable the atmospheric characterisation of hundreds of planets. Placing the satellite at L2 provides a cold and stable thermal environment, as well as a large field of regard to allow efficient time-critical observation of targets randomly distributed over the sky. A 1m class telescope is sufficiently large to achieve the necessary spectro-photometric precision. The spectral coverage (0.5-11 micron, goal 16 micron) and SNR to be achieved by EChO, thanks to its high stability and dedicated design, would enable a very accurate measurement of the atmospheric composition and structure of hundreds of exoplanets.

We perform numerical evolutions of cosmological scenarios using a standard general relativistic code in spherical symmetry. We concentrate on two different situations: initial matter distributions that are homogeneous and isotropic, and perturbations to those that respect the spherical symmetry. As matter models we consider the case of a pressureless perfect fluid, i.e. dust, and the case of a real massive scalar field oscillating around the minimum of the potential. Both types of matter have been considered as possible dark matter candidates in the cosmology literature, dust being closely related to the standard cold dark matter paradigm. We confirm that in the linear regime the perturbations associated with these types of matter grow in essentially the same way, the main difference being that in the case of a scalar field the dynamics introduce a cut-off in the power spectrum of the density perturbations at scales comparable with the Compton wavelength of the field. We also follow the evolutions well beyond the linear regime showing that both models are able to form structure. In particular we find that, once in the nonlinear regime, perturbations collapse faster in a universe dominated by dust. This is expected to delay the formation of the first structures in the scalar field dark matter scenario with respect to the standard cold dark matter one.

We investigate the existence of possible stable strange matter and related stability windows at finite temperature for different models that are generally applied to describe quark stars, namely, the quark-mass density dependent model, the MIT bag model and the Nambu-Jona-Lasinio model. We emphasize that, although the limits for stable strange matter depend on a comparison with the ground state of 56Fe, which is a zero temperature state, the quantity that has to be used in the search for strange matter in proto-quark stars is the free energy and we analyze stability windows up to temperatures of the order of 40 MeV. The effects of strong magnetic fields on stability windows are computed and the resulting mass-radius relations for different stages of the proto-quark star are analyzed.

The assumption underlying the existence of quark stars is based on the Bodmer-Witten conjecture. These authors have claimed that it is possible that the interior of a neutron-like star does not consist primarily of hadrons, but rather of the strange matter (SM). Strange matter is composed of deconfined quarks, including up, down and strange quarks, plus the leptons necessary to ensure charge neutrality and \beta-equilibrium. This possibility arises because a phase transition from hadronic to quark phase is possible at densities present in the interior of neutron stars. It has been argued that strange matter is the true ground state of all matter. If this is the case, as soon as the core of the star converts to the quark phase, the entire star converts. SM was first considered in calculations obtained within the MIT bag model framework. More sophisticated treatments for SM, based on the quark-mass density dependent model is considered in this work. We discuss next the following models for proto-quark stars: the QMDD model, and the MIT bag model.

We imaged the error box of a gamma-ray burst of the short (0.5 s), hard type (GRB 000313), with the BOOTES-1 experiment in southern Spain, starting 4 min after the gamma-ray event, in the I-band. A bright optical transient (OT 000313) with I = 9.4 +/- 0.1 was found in the BOOTES-1 image, close to the error box (3-sigma) provided by BATSE. Late time VRIK'-band deep observations failed to reveal an underlying host galaxy. If the OT 000313 is related to the short, hard GRB 000313, this would be the first optical counterpart ever found for this kind of events (all counterparts to date have been found for bursts of the long, soft type). The fact that only prompt optical emission has been detected (but no afterglow emission at all, as supported by theoretical models) might explain why no optical counterparts have ever been found for short, hard GRBs.This fact suggests that most short bursts might occur in a low-density medium and favours the models that relate them to binary mergers in very low-density enviroments.