We present a detailed method to simulating sensor distortions using a photon and electron Monte Carlo method. We use three dimensional electrostatic simulations to parameterize the perturbed electric field profile for non-ideal sensor details. We follow the conversion of simulated photons, and the subsequent response of the converted electrons to the electric field pattern. These non-ideal sensor details can be implemented efficiently in a Monte Carlo approach. We demonstrate that the non-ideal sensor distortions have a variety of observable consequence including the modification of the astrometric pattern, the distortion of the electron diffusion size and shape, and the distortion of flats. We show analytic validation of the diffusion physics, reproduce two kinds of edge distortion, and show qualitative validation of field-free regions, lithography errors, and fringing. We also demonstrate that there are two related effects of doping variation having different observable consequences. We show that field distortions from accumulated electrons lead to intensity-dependent point-spread-functions and the sub-linear variance in flats. The method is implemented in the Photon Simulator (PhoSim) and the code is publically available.

Development of the hardware, data analysis, and simulation techniques for large compact radio arrays dedicated to mapping the 21 cm line of neutral hydrogen gas has proven to be more difficult than imagined twenty years ago when such telescopes were first proposed. Despite tremendous technical and methodological advances, there are several outstanding questions on how to optimally calibrate and analyze such data. On the positive side, it has become clear that the outstanding issues are purely technical in nature and can be solved with sufficient development activity. Such activity will enable science across redshifts, from early galaxy evolution in the pre-reionization era to dark energy evolution at low redshift.

While Type Ia Supernovae (SNe Ia) are one of the most mature cosmological probes, the next era promises to be extremely exciting in the number of different ways SNe Ia are used to measure various cosmological parameters. Here we review the experiments in the 2020s that will yield orders of magnitudes more SNe Ia, and the new understandings and capabilities to constrain systematic uncertainties at a level to match these statistics. We then discuss five different cosmological probes with SNe Ia: the conventional Hubble diagram for measuring dark energy properties, the distance ladder for measuring the Hubble constant, peculiar velocities and weak lensing for measuring sigma8 and strong-lens measurements of H0 and other cosmological parameters. For each of these probes, we discuss the experiments that will provide the best measurements and also the SN Ia-related systematics that affect each one.

Our current understanding of the Universe is established through the pristine measurements of structure in the cosmic microwave background (CMB) and the distribution and shapes of galaxies tracing the large scale structure (LSS) of the Universe. One key ingredient that underlies cosmological observables is that the field that sources the observed structure is assumed to be initially Gaussian with high precision. Nevertheless, a minimal deviation from Gaussianityis perhaps the most robust theoretical prediction of models that explain the observed Universe; itis necessarily present even in the simplest scenarios. In addition, most inflationary models produce far higher levels of non-Gaussianity. Since non-Gaussianity directly probes the dynamics in the early Universe, a detection would present a monumental discovery in cosmology, providing clues about physics at energy scales as high as the GUT scale.

The hot dense environment of the early universe is known to have produced large numbers of baryons, photons, and neutrinos. These extreme conditions may have also produced other long-lived species, including new light particles (such as axions or sterile neutrinos) or gravitational waves. The gravitational effects of any such light relics can be observed through their unique imprint in the cosmic microwave background (CMB), the large-scale structure, and the primordial light element abundances, and are important in determining the initial conditions of the universe. We argue that future cosmological observations, in particular improved maps of the CMB on small angular scales, can be orders of magnitude more sensitive for probing the thermal history of the early universe than current experiments. These observations offer a unique and broad discovery space for new physics in the dark sector and beyond, even when its effects would not be visible in terrestrial experiments or in astrophysical environments. A detection of an excess light relic abundance would be a clear indication of new physics and would provide the first direct information about the universe between the times of reheating and neutrino decoupling one second later.

(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.

We describe a system to measure the Quantum Efficiency in the wavelength range of 300 nm to 1100 nm of 40x40 mm n-channel CCD sensors for the construction of the 3.2 gigapixel LSST focal plane. The technique uses a series of instrument to create a very uniform flux of photons of controllable intensity in the wavelength range of interest across the face the sensor. This allows the absolute Quantum Efficiency to be measured with an accuracy in the 7% range. This system will be part of a production facility at Brookhaven National Lab for the basic component of the LSST camera.

We employ electrostatic conversion drift calculations to match CCD pixel signal covariances observed in flat field exposures acquired using candidate sensor devices for the LSST Camera. We thus constrain pixel geometry distortions present at the end of integration, based on signal images recorded. We use available data from several operational voltage parameter settings to validate our understanding. Our primary goal is to optimize flux point-spread function (FPSF) estimation quantitatively, and thereby minimize sensor-induced errors which may limit performance in precision astronomy applications. We consider alternative compensation scenarios that will take maximum advantage of our understanding of this underlying mechanism in data processing pipelines currently under development. To quantitatively capture the pixel response in high-contrast/high dynamic range operational extrema, we propose herein some straightforward laboratory tests that involve altering the time order of source illumination on sensors, within individual test exposures. Hence the word {\it hysteretic} in the title of this paper.

Near-future astronomical survey experiments, such as LSST, possess system requirements of unprecedented fidelity that span photometry, astrometry and shape transfer. Some of these requirements flow directly to the array of science imaging sensors at the focal plane. Availability of high quality characterization data acquired in the course of our sensor development program has given us an opportunity to develop and test a framework for simulation and modeling that is based on a limited set of physical and geometric effects. In this paper we describe those models, provide quantitative comparisons between data and modeled response, and extrapolate the response model to predict imaging array response to astronomical exposure. The emergent picture departs from the notion of a fixed, rectilinear grid that maps photo-conversions to the potential well of the channel. In place of that, we have a situation where structures from device fabrication, local silicon bulk resistivity variations and photo-converted carrier patterns still accumulating at the channel, together influence and distort positions within the photosensitive volume that map to pixel boundaries. Strategies for efficient extraction of modeling parameters from routinely acquired characterization data are described. Methods for high fidelity illumination/image distribution parameter retrieval, in the presence of such distortions, are also discussed.

We study the potential of International Linear Collider (ILC) at sqrt(s) = 500 GeV to probe new dark matter motivated scenario where the bottom squark (sbottom) is the next to lightest supersymmetric particle. For this scenario, which is virtually impossible for the LHC to test, the ILC has a potential to cover a large fraction of the parameter space. The challenge is due to a very low energy of jets, below 20-30 GeV, which pushes the jet clustering and flavour tagging algorithms to their limits. The process of sbottom pair production was studied within the SiD detector concept. We demonstrate that ILC offers a unique opportunity to test the SUSY parameter space motivated by the sbottom neutralino co-annihilation scenario in cases when the sbottom production is kinematically accessible. The study was done with the full SiD simulation and reconstruction chain including all Standard Model and beam backgrounds.