The initial-final mass relation (IFMR) links the birth mass of a star to the mass of the compact remnant left at its death. While the relevance of the IFMR across astrophysics is universally acknowledged, not all of its fine details have yet been resolved. A new analysis of a few carbon-oxygen white dwarfs in old open clusters of the Milky Way led us to identify a kink in the IFMR, located over a range of initial masses, $1.65 \lesssim M_{\rm i}/M_{\odot} \lesssim 2.10$. The kink's peak in WD mass of $\approx 0.70-0.75 \, M_{\odot}$ is produced by stars with $M_{\rm i} \simeq 1.8 - 1.9 \, M_{\odot}$, corresponding to ages of about $1.8 - 1.7 $ Gyr. Interestingly, this peak coincides with the initial mass limit between low-mass stars that develop a degenerate helium core after central hydrogen exhaustion, and intermediate-mass stars that avoid electron degeneracy. We interpret the IFMR kink as the signature of carbon star formation in the Milky Way. This finding is critical to constraining the evolution and chemical enrichment of low-mass stars, and their impact on the spectrophotometric properties of galaxies.

Parameter estimation with non-Gaussian stochastic fields is a common challenge in astrophysics and cosmology. In this paper we advocate performing this task using the scattering transform, a statistical tool rooted in the mathematical properties of convolutional neural nets. This estimator can characterize a complex field without explicitly computing higher-order statistics, thus avoiding the high variance and dimensionality problems. It generates a compact set of coefficients which can be used as robust summary statistics for non-Gaussian information. It is especially suited for fields presenting localized structures and hierarchical clustering, such as the cosmological density field. To demonstrate its power, we apply this estimator to the cosmological parameter inference problem in the context of weak lensing. Using simulated convergence maps with realistic noise, the scattering transform outperforms the power spectrum and peak counts, and is on par with the state-of-the-art CNN. It retains the advantages of traditional statistical descriptors (it does not require any training nor tuning), has provable stability properties, allows to check for systematics, and importantly, the scattering coefficients are interpretable. It is a powerful and attractive estimator for observational cosmology and, in general, the study of physically-motivated fields.

We present the status of on-going detector development efforts for our joint NASA/CNES balloon-borne UV multi-object spectrograph, the Faint Intergalactic Redshifted Emission Balloon (FIREBall-2; FB-2). FB-2 demonstrates a new UV detector technology, the delta-doped Electron Multiplying CCD (EMCCD), in a low risk suborbital environment, to prove the performance of EMCCDs for future space missions and Technology Readiness Level (TRL) advancement. EMCCDs can be used in photon counting (PC) mode to achieve extremely low readout noise ($<$1 electron). Our testing has focused on reducing clock-induced-charge (CIC) through wave shaping and well depth optimization with a \nuvu V2 CCCP Controller, measuring CIC at 0.001 e$^{-}$/pixel/frame. This optimization also includes methods for reducing dark current, via cooling, and substrate voltage levels. We discuss the challenges of removing cosmic rays, which are also amplified by these detectors, as well as a data reduction pipeline designed for our noise measurement objectives. FB-2 flew in 2018, providing the first time an EMCCD was used for UV observations in the stratosphere. FB-2 is currently being built up to fly again in 2020, and improvements are being made to the EMCCD to continue optimizing its performance for better noise control.

Double white dwarf (double-WD) binaries may merge within a Hubble time and produce high-mass WDs. Compared to other high-mass WDs, the double-WD merger products have higher velocity dispersion because they are older. With the power of Gaia data, we show strong evidence for double-WD merger products among high-mass WDs by analyzing the transverse-velocity distribution of more than a thousand high-mass WDs (0.8--1.3 $M_\odot$). We estimate that the fraction of double-WD merger products in our sample is about 20 %. We also obtain a precise double-WD merger rate and its mass dependence. Our merger rate estimates are close to binary population synthesis results and support the idea that double-WD mergers may contribute to a significant fraction of type Ia supernovae.

N-body simulations are essential tools in physical cosmology to understand the large-scale structure (LSS) formation of the Universe. Large-scale simulations with high resolution are important for exploring the substructure of universe and for determining fundamental physical parameters like neutrino mass. However, traditional particle-mesh (PM) based algorithms use considerable amounts of memory, which limits the scalability of simulations. Therefore, we designed a two-level PM algorithm CUBE towards optimal performance in memory consumption reduction. By using the fixed-point compression technique, CUBE reduces the memory consumption per N-body particle toward 6 bytes, an order of magnitude lower than the traditional PM-based algorithms. We scaled CUBE to 512 nodes (20,480 cores) on an Intel Cascade Lake based supercomputer with $\simeq$95\% weak-scaling efficiency. This scaling test was performed in "Cosmo-$\pi$" -- a cosmological LSS simulation using $\simeq$4.4 trillion particles, tracing the evolution of the universe over $\simeq$13.7 billion years. To our best knowledge, Cosmo-$\pi$ is the largest completed cosmological N-body simulation. We believe CUBE has a huge potential to scale on exascale supercomputers for larger simulations.

By comparing two age indicators of high-mass white dwarfs derived from Gaia data, two discoveries have been made recently: one is the existence of a cooling anomaly that produces the Q branch structure on the Hertzsprung--Russell diagram, the other is the existence of double-white-dwarf merger products. The former poses a challenge for white dwarf cooling models, and the latter has implications on binary evolution and type-Ia supernovae.

Recently, the power of Gaia data has revealed an enhancement of high-mass white dwarfs (WDs) on the Hertzsprung--Russell diagram, called the Q branch. This branch is located at the high-mass end of the recently identified crystallization branch. Investigating its properties, we find that the number density and velocity distribution on the Q branch cannot be explained by the cooling delay of crystallization alone, suggesting the existence of an extra cooling delay. To quantify this delay, we statistically compare two age indicators -- the dynamical age inferred from transverse velocity, and the photometric isochrone age -- for more than one thousand high-mass WDs (1.08--1.23 $M_\odot$) selected from Gaia Data Release 2. We show that about 6 % of the high-mass WDs must experience an 8 Gyr extra cooling delay on the Q branch, in addition to the crystallization and merger delays. This cooling anomaly is a challenge for WD cooling models. We point out that $^{22}$Ne settling in C/O-core WDs could account for this extra cooling delay.

The binary star Par 1802 in the Orion Nebula presents an interesting puzzle in the field of stellar dynamics and evolution. Binary systems such as Par 1802 are thought to form from the same natal material and thus the stellar members are expected to have very similar physical attributes. However, Par 1802's stars have significantly different temperatures despite their identical (within $3\%$) masses of about $0.39$ solar mass. The leading proof-of-concept idea is that a third companion gravitationally induced the two stars to orbit closer than their Roche-limit, which facilitated heating through tidal effects. Here we expand on this idea and study the three-body dynamical evolution of such a system, including tidal and pre-main-sequence evolution. We also include tidal heating and mass transfer at the onset of Roche-limit Crossing. We show, as a proof-of-concept, that mass transfer combined with tidal heating can naturally explain the observed temperature discrepancy. We also predict the orbital configuration of the possible tertiary companion. Finally we suggest that the dynamical evolution of such a system has pervasive consequences. We expect an abundance of systems to undergo mass transfer during their pre-main-sequence time, which can cause temperature differences.

We present lifetime measurements using beam-foil techniques for radiative transitions from the 3$p^4$($^1S$)4$s$ $^2S_{1/2}$, 3$p^4$($^3P$)5$s$ $^2P_{1/2,3/2}$, and 3$p^4$($^3P$)3$d$ $^2F_{5/2}$ levels in Cl I and the corresponding results of the oscillator strengths for transitions at 1004.68, 1079.88, 1090.73, and 1094.77 \AA, respectively. We compare our experimental results with available theoretical calculations and astronomical observations in an effort to resolve discrepancies among them.

In the axion monodromy inflation, the inflation is driven by the axion with super-Planckian field values in a monomial potential with superimposed sinusoidal modulations. The coupling of the axion to massless gauge fields can induce copious particle production during inflation, resulting in large non-Gaussian curvature perturbation that leads to the formation of primordial black holes. In this paper, we explore the parameter space in the axion monodromy inflation model that favors the formation of primordial black holes with masses ranging from $10^8$ grams to $20$ solar masses. We also study the associated gravitational waves and their detection in pulsar timing arrays and interferometry experiments.

We study the growth of superhorizon modes in the curvature perturbation during an ultra-slow-roll or a large-$\eta$ phase in single-field inflation. In a simple toy model, we derive the two-point correlation function of the curvature perturbation and show that the requirement for causality restricts the growth rate and hence puts a lower limit on the value of $\eta$. The toy model is then realized by considering an inflation potential with an inflection point. Our study is useful to assessing the growth of the curvature perturbation that seeds the formation of primordial black holes.

We report lifetimes, branching fractions, and the resulting oscillator strengths for transitions within the P II multiplet (3s23p2 3P - 3s3p3 3Po) at 1308 {\AA}. These comprehensive beam-foil measurements, which are the most precise set currently available experimentally, resolve discrepancies involving earlier experimental and theoretical results. Interstellar phosphorus abundances derived from {\lambda}1308 can now be interpreted with greater confidence. In the course of our measurements, we also obtained an experimental lifetime for the 3p4s 3P0o level of P IV. This lifetime agrees well with the available theoretical calculation.

Trapped inflation has been proposed to provide a successful inflation with a steep potential. We discuss the formation of primordial black holes in the trapped inflationary scenario. We show that primordial black holes are naturally produced during inflation with a steep trapping potential. In particular, we have given a recipe for an inflaton potential with which particle production can induce large non-Gaussian curvature perturbation that leads to the formation of high stellar-mass primordial black holes. These primordial black holes could be dark matter observed by the LIGO detectors through a binary black-hole merger. At the end, we have given an attempt to realize the required inflaton potential in the axion monodromy inflation, and discussed the gravitational waves sourced by the particle production.

A number of exciting concepts are under development for Flagship, Probe class, Explorer class, and Suborbital class NASA missions in the ultraviolet/optical spectral ranges. These missions will depend on high performance silicon detector arrays being delivered affordably and in high numbers. In a focused effort we have advanced delta-doping technology to high throughput and high yield wafer-scale processing, encompassing a multitude of state-of-the-art silicon-based detector formats and designs. As part of this technology advancement and in preparation for upcoming missions, we have embarked on a number of field observations, instrument integrations, and independent evaluations of delta-doped arrays. In this paper, we present recent data and innovations from the Advanced Detectors and Systems program at JPL, including two-dimensional doping technology; our end-to-end post-fabrication processing of high performance UV/Optical/NIR arrays; and advanced coatings for detectors and optical elements. Additionally, we present examples of past, in-progress, and planned observations and deployments of delta-doped arrays.

We report on multilayer high efficiency antireflection coating (ARC) design and development for use at UV wavelengths on CCDs and other Si-based detectors. We have previously demonstrated a set of single-layer coatings, which achieve >50% quantum efficiency (QE) in four bands from 130 to 300 nm. We now present multilayer coating designs that significantly outperform our previous work between 195 and 215 nm. Using up to 11 layers, we present several model designs to reach QE above 80%. We also demonstrate the successful performance of 5 and 11 layer ARCs on silicon and fused silica substrates. Finally, we present a five-layer coat- ing deposited onto a thinned, delta-doped CCD and demonstrate external QE greater than 60% between 202 and 208 nm, with a peak of 67.6% at 206 nm.

Lifetime measurements using beam foil techniques for the $5s5p^{2}$ $^{2}D_{3/2,5/2}$ levels in Sn II are presented. The resulting oscillator strengths for transitions at 1699.4, 1831.8 and 1811.2 \AA{} are reported. We also studied these levels with multi-configuration Dirac Hartree-Fock (MCDHF) calculations using a development version of the GRASP2K package. Our experimental and theoretical results are compared with other available studies.

A numerical study of a pseudoscalar inflation having an axion-photon-like coupling is performed by solving numerically the coupled differential equations of motion for inflaton and photon mode functions from the onset of inflation to the end of reheating. The backreaction due to particle production is also included self-consistently. We find that this particular inflation model realizes the idea of a warm inflation in which a steady thermal bath is established by the particle production. In most cases this thermal bath exceeds the amount of radiation released in the reheating process. In the strong coupling regime, the transition from the inflationary to the radiation-dominated phase does not involve either a preheating or reheating process. In addition, energy density peaks produced near the end of inflation may lead to the formation of primordial black holes.

We present the results of lifetime measurements made using beam-foil techniques on levels of astrophysical interest in Pb II producing lines at 1203.6 A ($6s6p^{2}$ $^{2}D_{3/2}$) and 1433.9 A ($6s^{2}6d$ $^{2}D_{3/2}$). We also report the first detection of the Pb II $\lambda1203$ line in the interstellar medium (ISM) from an analysis of archival spectra acquired by the Space Telescope Imaging Spectrograph (STIS) onboard the Hubble Space Telescope (HST). The oscillator strengths derived from our experimental lifetimes for Pb II $\lambda\lambda1203$, $1433$ are generally consistent with recent theoretical results, including our own relativistic calculations. Our analysis of high-resolution HST/STIS spectra helps to confirm the relative strengths of the Pb II $\lambda\lambda 1203$, $1433$ lines. However, the oscillator strength that we obtain for Pb II $\lambda1433$ ($0.321\pm0.034$) is significantly smaller than earlier theoretical values, which have been used to study the abundance of Pb in the ISM. Our revised oscillator strength for $\lambda1433$ yields an increase in the interstellar abundance of Pb of 0.43 dex over determinations based on the value given by Morton, indicating that the depletion of Pb onto interstellar dust grains is less severe than previously thought.

We discuss the generation of primordial magnetic fields during inflation in the dilaton-axion electromagnetism, in which the dilaton and axion dynamics are introduced in terms of two time dependent functions of the cosmic scale factor, $I(a) F^2/4$ and $J(a) F\tilde{F}/4$, respectively, where $F$ is the electromagnetic field strength and $\tilde{F}$ is its dual. We study the form of $J(a)$ that can generate a large seed magnetic field in the weak coupling regime, $I(a)<1$. Although the $J(a)$ function is model dependent, the axion-photon coupling indeed opens up a new window for a successful inflationary magnetogenesis.

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.

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.

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{\deg} to +20{\deg}, 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.

The Pierre Auger Observatory is the world's largest cosmic ray observatory. Our current exposure reaches nearly 40,000 km$^2$ str and provides us with an unprecedented quality data set. The performance and stability of the detectors and their enhancements are described. Data analyses have led to a number of major breakthroughs. Among these we discuss the energy spectrum and the searches for large-scale anisotropies. We present analyses of our X$_{max}$ data and show how it can be interpreted in terms of mass composition. We also describe some new analyses that extract mass sensitive parameters from the 100% duty cycle SD data. A coherent interpretation of all these recent results opens new directions. The consequences regarding the cosmic ray composition and the properties of UHECR sources are briefly discussed.

Joint contributions of the Pierre Auger and Telescope Array Collaborations to the 33rd International Cosmic Ray Conference, Rio de Janeiro, Brazil, July 2013: cross-calibration of the fluorescence telescopes, large scale anisotropies and mass composition.

Contributions of the Pierre Auger Collaboration to the 33rd International Cosmic Ray Conference, Rio de Janeiro, Brazil, July 2013

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.

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.

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.

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\"ue and averaged monthly models, the utility of the GDAS data is shown.

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.

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.

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.

Technical reports on operations and monitoring of the Pierre Auger Observatory