The apparent superposition of galaxies with other astrophysical objects along the line of sight, a problem known as blending, will be a major challenge for upcoming, ground-based, deep, photometric galaxy surveys, such as the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST). Blending contributes to the systematic error budget of weak lensing studies by perturbing object detection and affecting photometric and shape measurements. Existing deblenders suffer from the lack of flexible yet accurate models of galaxy morphologies and therefore rely on assumptions (analytic profiles, symmetry, sparsity in a profile basis) to isolate galaxies within a blended scene. In this paper, we propose instead to use generative models based on deep neural networks, namely variational autoencoders (VAE), to learn a probabilistic model directly from data. Specifically, we train a VAE on images of centred, isolated galaxies, which we reuse in a second VAE-like neural network in charge of deblending galaxies. We train our networks on simulated images created with the GalSim software, including all six LSST bandpass filters as well as the visible and near-infrared bands of the Euclid satellite, as our method naturally generalises to multiple bands and data from multiple instruments. We validate our model and quantify deblending performance by measuring reconstruction errors in galaxy shapes and magnitudes. We obtain median errors on ellipticities between $\pm{0.01}$ and on $r$-band magnitude between $\pm{0.05}$ in most cases and shear multiplicative biases close to $10^{-2}$ in the optimal configuration. Finally, we discuss future challenges about training on real data and obtain encouraging results when applying transfer learning.

(Abridged) We describe here the most ambitious survey currently planned in the optical, the Large Synoptic Survey Telescope (LSST). A vast array of science will be enabled by a single wide-deep-fast sky survey, and LSST will have unique survey capability in the faint time domain. The LSST design is driven by four main science themes: probing dark energy and dark matter, taking an inventory of the Solar System, exploring the transient optical sky, and mapping the Milky Way. LSST will be a wide-field ground-based system sited at Cerro Pach\'{o}n in northern Chile. The telescope will have an 8.4 m (6.5 m effective) primary mirror, a 9.6 deg$^2$ field of view, and a 3.2 Gigapixel camera. The standard observing sequence will consist of pairs of 15-second exposures in a given field, with two such visits in each pointing in a given night. With these repeats, the LSST system is capable of imaging about 10,000 square degrees of sky in a single filter in three nights. The typical 5$\sigma$ point-source depth in a single visit in $r$ will be $\sim 24.5$ (AB). The project is in the construction phase and will begin regular survey operations by 2022. The survey area will be contained within 30,000 deg$^2$ with $\delta<+34.5^\circ$, and will be imaged multiple times in six bands, $ugrizy$, covering the wavelength range 320--1050 nm. About 90\% of the observing time will be devoted to a deep-wide-fast survey mode which will uniformly observe a 18,000 deg$^2$ region about 800 times (summed over all six bands) during the anticipated 10 years of operations, and yield a coadded map to $r\sim27.5$. The remaining 10\% of the observing time will be allocated to projects such as a Very Deep and Fast time domain survey. The goal is to make LSST data products, including a relational database of about 32 trillion observations of 40 billion objects, available to the public and scientists around the world.

The IceCube Neutrino Observatory is a cubic-kilometer-scale high-energy neutrino detector built into the ice at the South Pole. Construction of IceCube, the largest neutrino detector built to date, was completed in 2011 and enabled the discovery of high-energy astrophysical neutrinos. We describe here the design, production, and calibration of the IceCube digital optical module (DOM), the cable systems, computing hardware, and our methodology for drilling and deployment. We also describe the online triggering and data filtering systems that select candidate neutrino and cosmic ray events for analysis. Due to a rigorous pre-deployment protocol, 98.4% of the DOMs in the deep ice are operating and collecting data. IceCube routinely achieves a detector uptime of 99% by emphasizing software stability and monitoring. Detector operations have been stable since construction was completed, and the detector is expected to operate at least until the end of the next decade.

In cosmic ray air showers, the muon lateral separation from the center of the shower is a measure of the transverse momentum that the muon parent acquired in the cosmic ray interaction. IceCube has observed cosmic ray interactions that produce muons laterally separated by up to 400 m from the shower core, a factor of 6 larger distance than previous measurements. These muons originate in high pT (> 2 GeV/c) interactions from the incident cosmic ray, or high-energy secondary interactions. The separation distribution shows a transition to a power law at large values, indicating the presence of a hard pT component that can be described by perturbative quantum chromodynamics. However, the rates and the zenith angle distributions of these events are not well reproduced with the cosmic ray models tested here, even those that include charm interactions. This discrepancy may be explained by a larger fraction of kaons and charmed particles than is currently incorporated in the simulations.

The measurement of muon energy is critical for many analyses in large Cherenkov detectors, particularly those that involve separating extraterrestrial neutrinos from the atmospheric neutrino background. Muon energy has traditionally been determined by measuring the specific energy loss (dE/dx) along the muon's path and relating the dE/dx to the muon energy. Because high-energy muons (E_mu > 1 TeV) lose energy randomly, the spread in dE/dx values is quite large, leading to a typical energy resolution of 0.29 in log10(E_mu) for a muon observed over a 1 km path length in the IceCube detector. In this paper, we present an improved method that uses a truncated mean and other techniques to determine the muon energy. The muon track is divided into separate segments with individual dE/dx values. The elimination of segments with the highest dE/dx results in an overall dE/dx that is more closely correlated to the muon energy. This method results in an energy resolution of 0.22 in log10(E_mu), which gives a 26% improvement. This technique is applicable to any large water or ice detector and potentially to large scintillator or liquid argon detectors.

IceTop, the surface component of the IceCube Neutrino Observatory at the South Pole, is an air shower array with an area of 1 km2. The detector allows a detailed exploration of the mass composition of primary cosmic rays in the energy range from about 100 TeV to 1 EeV by exploiting the correlation between the shower energy measured in IceTop and the energy deposited by muons in the deep ice. In this paper we report on the technical design, construction and installation, the trigger and data acquisition systems as well as the software framework for calibration, reconstruction and simulation. Finally the first experience from commissioning and operating the detector and the performance as an air shower detector will be discussed.

We report on the results of the search for extremely-high energy (EHE) neutrinos with energies above $10^7$ GeV obtained with the partially ($\sim$30%) constructed IceCube in 2007. From the absence of signal events in the sample of 242.1 days of effective livetime, we derive a 90% C.L. model independent differential upper limit based on the number of signal events per energy decade at $E^2 \phi_{\nu_e+\nu_\mu+\nu_\tau}\simeq 1.4 \times 10^{-6}$ GeV cm$^{-2}$ sec$^{-1}$ sr$^{-1}$ for neutrinos in the energy range from $3\times10^7$ to $3\times10^9$ GeV.

IceCube Collaboration Contributions to the 2009 International Cosmic Ray Conference

Over 5,000 PMTs are being deployed at the South Pole to compose the IceCube neutrino observatory. Many are placed deep in the ice to detect Cherenkov light emitted by the products of high-energy neutrino interactions, and others are frozen into tanks on the surface to detect particles from atmospheric cosmic ray showers. IceCube is using the 10-inch diameter R7081-02 made by Hamamatsu Photonics. This paper describes the laboratory characterization and calibration of these PMTs before deployment. PMTs were illuminated with pulses ranging from single photons to saturation level. Parameterizations are given for the single photoelectron charge spectrum and the saturation behavior. Time resolution, late pulses and afterpulses are characterized. Because the PMTs are relatively large, the cathode sensitivity uniformity was measured. The absolute photon detection efficiency was calibrated using Rayleigh-scattered photons from a nitrogen laser. Measured characteristics are discussed in the context of their relevance to IceCube event reconstruction and simulation efforts.

A photon induced shower at $E_{prim}\ge 10^{18}$ eV exhibits very specific features and is different from a hadronic one. At such energies, the LPM effect delays in average the first interactions of the photon in the atmosphere and hence slows down the whole shower development. They also have a smaller muonic content than hadronic ones. The response of a surface detector such as that of the Auger Observatory to these specific showers is thus different and has to be accounted for in order to enable potential photon candidates reconstruction correctly. The energy reconstruction in particular has to be adapted to the late development of photon showers. We propose in this article a method for the reconstruction of the energy of photon showers with a surface detector. The key feature of this method is to rely explicitly on the development stage of the shower. This approach leads to very satisfactory results ($\simeq 20%$). At even higher energies ($5.10^{19}$ eV and above) the probability for the photon to convert into a pair of e$^+$e$^-$ in the geomagnetic field becomes non negligible and requires a different function to evaluate the energy with the proposed method. We propose several approaches to deal with this issue in the scope of the establishment of an upper bound on the photon fraction in UHECR.