In this paper we investigate the cosmological dynamics of an up to cubic curvature correction to General Relativity (GR) known as Cosmological Einsteinian Cubic Gravity (CECG), whose vacuum spectrum consists of the graviton exclusively and its cosmology is well-posed as an initial value problem. We are able to uncover the global asymptotic structure of the phase space of this theory. It is revealed that an inflationary matter-dominated bigbang is the global past attractor which means that inflation is the starting point of any physically meaningful cosmic history. Given that higher order curvature corrections to GR are assumed to influence the cosmological dynamics at early times -- high energies/large curvature limit -- the late-time inflation can not be a consequence of the up to cubic order curvature modifications. We confirm this assumption by showing that late-time acceleration of the expansion in the CECG model is possible only if add a cosmological constant term.

The production of very high energy muons inside an extensive air shower is observable at $\nu$ telescopes and sensitive to the composition of the primary cosmic ray. Here we discuss five different sources of these muons: pion and kaon decays; charmed hadron decays; rare decays of unflavored mesons; photon conversion into a muon pair; and photon conversion into a $J/\psi$ vector meson decaying into muons. We solve the cascade equations for a $10^{10.5}$ GeV proton primary and find that unflavored mesons and gamma conversions are the two main sources of $E\ge 10^{8.5}$ GeV muons, while charm decays dominate at $10^{5.5}\,{\rm GeV}< E< 10^{8.5}\,{\rm GeV}$. In inclined events one of these muons may deposite a large fraction of its energy near the surface, implying fluctuations in the longitudinal profile of the shower and in the muon to electron count at the ground level. In particular, we show that 1 out of 6 proton showers of $10^{10.5}$ GeV include an $E>10^6$ GeV deposition within 500 g/cm$^2$, while only in 1 out of 330 showers it is above $10^7$ GeV. We also show that the production of high energy muons is very different in proton, iron or photon showers ({e.g., conversions $\gamma\to \mu^+ \mu^-$ are the main source of $E\ge 10^4$ GeV muons in photon showers). Finally, we use Monte Carlo simulations to discuss the validity of our results.

The PTOLEMY project aims to develop a scalable design for a Cosmic Neutrino Background (CNB) detector, the first of its kind and the only one conceived that can look directly at the image of the Universe encoded in neutrino background produced in the first second after the Big Bang. The scope of the work for the next three years is to complete the conceptual design of this detector and to validate with direct measurements that the non-neutrino backgrounds are below the expected cosmological signal. In this paper we discuss in details the theoretical aspects of the experiment and its physics goals. In particular, we mainly address three issues. First we discuss the sensitivity of PTOLEMY to the standard neutrino mass scale. We then study the perspectives of the experiment to detect the CNB via neutrino capture on tritium as a function of the neutrino mass scale and the energy resolution of the apparatus. Finally, we consider an extra sterile neutrino with mass in the eV range, coupled to the active states via oscillations, which has been advocated in view of neutrino oscillation anomalies. This extra state would contribute to the tritium decay spectrum, and its properties, mass and mixing angle, could be studied by analyzing the features in the beta decay electron spectrum.

We report on the serendipitous observations of Solar System objects imaged during the High cadence Transient Survey (HiTS) 2014 observation campaign. Data from this high cadence, wide field survey was originally analyzed for finding variable static sources using Machine Learning to select the most-likely candidates. In this work we search for moving transients consistent with Solar System objects and derive their orbital parameters. We use a simple, custom detection algorithm to link trajectories and assume Keplerian motion to derive the asteroid's orbital parameters. We use known asteroids from the Minor Planet Center (MPC) database to assess the detection efficiency of the survey and our search algorithm. Trajectories have an average of nine detections spread over 2 days, and our fit yields typical errors of $\sigma_a\sim 0.07 ~{\rm AU}$, $\sigma_{\rm e} \sim 0.07 $ and $\sigma_i\sim 0.^{\circ}5~ {\rm deg}$ in semi-major axis, eccentricity, and inclination respectively for known asteroids in our sample. We extract 7,700 orbits from our trajectories, identifying 19 near Earth objects, 6,687 asteroids, 14 Centaurs, and 15 trans-Neptunian objects. This highlights the complementarity of supernova wide field surveys for Solar System research and the significance of machine learning to clean data of false detections. It is a good example of the data--driven science that LSST will deliver.

We report the discovery of distant RR Lyrae stars, including the most distant known in the Milky Way, using data taken in the $g-$band with the Dark Energy Camera as part of the High cadence Transient Survey (HiTS; 2014 campaign). We detect a total of 173 RR Lyrae stars over a ~120 deg^2 area, including both known RR Lyrae and new detections. The heliocentric distances d_H of the full sample range from 9 to >200 kpc, with 18 of them beyond 90 kpc. We identify three sub-groups of RR Lyrae as members of known systems: the Sextans dwarf spheroidal galaxy, for which we report 46 new discoveries, and the ultra-faint dwarf galaxies Leo IV and Leo V. Following an MCMC methodology, we fit spherical and ellipsoidal profiles of the form rho(R) ~ R^n to the radial density distribution of RR Lyrae in the Galactic halo. The best fit corresponds to the spherical case, for which we obtain a simple power-law index of n=-4.17^{+0.18}_{-0.20}, consistent with recent studies made with samples covering shorter distances. The pulsational properties of the outermost RR Lyrae in the sample (d_H>90 kpc) differ from the ones in the halo population at closer distances. The distribution of the stars in a Period-Amplitude diagram suggest they belong to Oosterhoff-intermediate or Oosterhoff II groups, similar to what is found in the ultra-faint dwarf satellites around the Milky Way. The new distant stars discovered represent an important addition to the few existing tracers of the Milky Way potential in the outer halo.

We present ultraviolet, optical and infrared photometry and optical spectroscopy of the type Ic superluminous supernova (SLSN) Gaia16apd (= SN 2016eay), covering its evolution from 26 d before the $g$-band peak to 234.1 d after the peak. Gaia16apd was followed as a part of the NOT Unbiased Transient Survey (NUTS). It is one of the closest SLSNe known ($z = 0.102\pm0.001$), with detailed optical and ultraviolet (UV) observations covering the peak. Gaia16apd is a spectroscopically typical type Ic SLSN, exhibiting the characteristic blue early spectra with O II absorption, and reaches a peak $M_{g} = -21.8 \pm 0.1$ mag. However, photometrically it exhibits an evolution intermediate between the fast- and slowly-declining type Ic SLSNe, with an early evolution closer to the fast-declining events. Together with LSQ12dlf, another SLSN with similar properties, it demonstrates a possible continuum between fast- and slowly-declining events. It is unusually UV-bright even for a SLSN, reaching a non-$K$-corrected $M_{uvm2} \simeq -23.3$ mag, the only other type Ic SLSN with similar UV brightness being SN 2010gx. Assuming that Gaia16apd was powered by magnetar spin-down, we derive a period of $P = 1.9\pm0.2$ ms and a magnetic field of $B = 1.9\pm0.2 \times 10^{14}$ G for the magnetar. The estimated ejecta mass is between 8 and 16 $\mathrm{M}_{\odot}$ and the kinetic energy between 1.3 and $2.5 \times 10^{52}$ erg, depending on opacity and assuming that the entire ejecta is swept up into a thin shell. Despite the early photometric differences, the spectra at late times are similar to slowly-declining type Ic SLSNe, implying that the two subclasses originate from similar progenitors.

We present the first results of the High cadence Transient Survey (HiTS), a survey whose objective is to detect and follow up optical transients with characteristic timescales from hours to days, especially the earliest hours of supernova (SN) explosions. HiTS uses the Dark Energy Camera (DECam) and a custom made pipeline for image subtraction, candidate filtering and candidate visualization, which runs in real-time to be able to react rapidly to the new transients. We discuss the survey design, the technical challenges associated with the real-time analysis of these large volumes of data and our first results. In our 2013, 2014 and 2015 campaigns we have detected more than 120 young SN candidates, but we did not find a clear signature from the short-lived SN shock breakouts (SBOs) originating after the core collapse of red supergiant stars, which was the initial science aim of this survey. Using the empirical distribution of limiting-magnitudes from our observational campaigns we measured the expected recovery fraction of randomly injected SN light curves which included SBO optical peaks produced with models from Tominaga et al. (2011) and Nakar & Sari (2010). From this analysis we cannot rule out the models from Tominaga et al. (2011) under any reasonable distributions of progenitor masses, but we can marginally rule out the brighter and longer-lived SBO models from Nakar & Sari (2010) under our best-guess distribution of progenitor masses. Finally, we highlight the implications of this work for future massive datasets produced by astronomical observatories such as LSST.

The investigation of the non-linearity of the Leavitt law is a topic that began more than seven decades ago, when some of the studies in this field found that the Leavitt law has a break at about ten days. The goal of this work is to investigate a possible statistical cause of this non-linearity. By applying linear regressions to OGLE-II and OGLE-IV data, we find that, in order to obtain the Leavitt law by using linear regression, robust techniques to deal with influential points and/or outliers are needed instead of the ordinary least-squares regression traditionally used. In particular, by using $M$- and $MM$-regressions we establish firmly and without doubts the linearity of the Leavitt law in the Large Magellanic Cloud, without rejecting or excluding Cepheid data from the analysis. This implies that light curves of Cepheids suggesting blending, bumps, eclipses or period changes, do not affect the Leavitt law for this galaxy. For the SMC, including this kind of Cepheids, it is not possible to find an adequate model, probably due to the geometry of the galaxy. In that case, a possible influence of these stars could exist.

We present the sensitivity of HAWC to Gamma Ray Bursts (GRBs). HAWC is a very high-energy gamma-ray observatory currently under construction in Mexico at an altitude of 4100 m. It will observe atmospheric air showers via the water Cherenkov method. HAWC will consist of 300 large water tanks instrumented with 4 photomultipliers each. HAWC has two data acquisition (DAQ) systems. The main DAQ system reads out coincident signals in the tanks and reconstructs the direction and energy of individual atmospheric showers. The scaler DAQ counts the hits in each photomultiplier tube (PMT) in the detector and searches for a statistical excess over the noise of all PMTs. We show that HAWC has a realistic opportunity to observe the high-energy power law components of GRBs that extend at least up to 30 GeV, as it has been observed by Fermi LAT. The two DAQ systems have an energy threshold that is low enough to observe events similar to GRB 090510 and GRB 090902b with the characteristics observed by Fermi LAT. HAWC will provide information about the high-energy spectra of GRBs which in turn could help to understanding about e-pair attenuation in GRB jets, extragalactic background light absorption, as well as establishing the highest energy to which GRBs accelerate particles.

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.