The "crosshatch" pattern is a recurring "feature" of HgCdTe arrays, specifically the Teledyne HAWAII family of near infrared detectors. It is a fixed pattern of high frequency QE variations along 3 directions generally thought to be related to the crystal structure of HgCdTe. The pattern is evident in detectors used in Hubble WFC3, WISE, JWST, and in candidate detectors for Euclid and WFIRST. Using undersampled point sources projected onto a HAWAII-2RG detector, we show that the pattern induces photometric variations that are not removed by a flat-field calibration, thus demonstrating that the QE variations occur on scales smaller than the 18 micron pixels. Our testbed is the Precision Projector Laboratory's astronomical scene generator, which can rapidly characterize the full detector by scanning thousands of undersampled spots. If not properly calibrated, detectors showing strong crosshatch may induce correlated errors in photometry, astrometry, spectroscopy, and shape measurements.

The "brighter-fatter" (BF) effect is a phenomenon (originally discovered in charge coupled devices) in which the size of the detector point spread function (PSF) increases with brightness. We present, for the first time, laboratory measurements demonstrating the existence of the effect in a Hawaii-2RG HgCdTe near infrared (NIR) detector. We use the Precision Projector Laboratory, a JPL facility for emulating astronomical observations with UV/VIS/NIR detectors, to project about 17,000 point sources onto the detector to stimulate the effect. After calibrating the detector for nonlinearity with flat-fields, we find evidence that charge is nonlinearly shifted from bright pixels to neighboring pixels during exposures of point sources, consistent with the existence of a BF-type effect. The Wide Field Infrared Survey Telescope (WFIRST) by NASA will use similar detectors to measure weak gravitational lensing from the shapes of hundreds of million of galaxies in the NIR. The WFIRST PSF size must be calibrated to approximately 0.1 percent to avoid biased inferences of dark matter and dark energy parameters; therefore further study and calibration of the BF effect in realistic images will be crucial.

Cosmic acceleration is the most surprising cosmological discovery in many decades. Testing and distinguishing among possible explanations requires cosmological measurements of extremely high precision probing the full history of cosmic expansion and structure growth and, ideally, compare and contrast matter and relativistic tracers of the gravity potential. This program is one of the defining objectives of the Wide-Field Infrared Survey Telescope (WFIRST), as set forth in the New Worlds, New Horizons report (NWNH) in 2010. The WFIRST mission has the ability to improve these measurements by 1-2 orders of magnitude compared to the current state of the art, while simultaneously extending their redshift grasp, greatly improving control of systematic effects, and taking a unified approach to multiple probes that provide complementary physical information and cross-checks of cosmological results. We describe in this annual report the activities of the Science Investigation Team (SIT) "Cosmology with the High Latitude Survey (HLS)" during the year 2017. This team was selected by NASA in December 2015 in order to address the stringent challenges of the WFIRST dark energy (DE) program through the Project's formulation phase. This SIT has elected to jointly address Galaxy Redshift Survey, Weak Lensing and Cluster Growth and thus fully embrace the fact that the imaging and spectroscopic elements of the HLS will be realized as an integrated observing program, and they jointly impose requirements on performance and operations. WFIRST is designed to be able to deliver a definitive result on the origin of cosmic acceleration. It is not optimized for Figure of Merit sensitivity but for control of systematic uncertainties and for having multiple techniques each with multiple cross-checks. Our SIT work focuses on understanding the potential systematics in the WFIRST DE measurements.

As astronomical observations from space benefit from improved sensitivity, the effectiveness of scientific programs is becoming limited by systematics that often originate in poorly understood image sensor behavior. Traditional, bottom-up detector characterization methods provide one way to model underlying detector physics, and generate ever more faithful numerical simulations, but this approach is vulnerable to preconceptions and over-simplification. The alternative top-down approach is laboratory emulation, which enables observation, calibration, and analysis scenarios to be tested without relying on a complete understanding of the underlying physics. This complements detector characterization and simulation efforts by testing their validity. We describe a laboratory facility and experimental testbed that supports the emulation of a wide range of mission concepts such as gravitational weak lensing measurements by WFIRST and high precision spectrophotometry of transiting exoplanets by JWST. An Offner relay projects readily customizable "scenes" (e.g. stars, galaxies, spectra) with very low optical aberration over the full area of a typical optical or near infrared image sensor. f/8 and slower focal ratios may be selected, spanning those of most proposed space missions and approximating the point spread function (PSF) size of seeing limited ground based surveys. Diffraction limited PSFs are projected over a wide field of view and wavelength range to deliver highly predictable image patterns down to sub-pixel scales with stable intensity and fine motion control. The testbed enables realistic validation of detector performance on science-like images, which aids mission design and survey strategy, as well as targeted investigations of various detector effects.

We study the nonlinearity (NL) in the conversion from charge to voltage in infrared detectors (HXRG) for use in precision astronomy. We present laboratory measurements of the NL function of a H2RG detector and discuss the accuracy to which it would need to be calibrated in future space missions to perform cosmological measurements through the weak gravitational lensing technique. In addition, we present an analysis of archival data from the infrared H1RG detector of the Wide Field Camera 3 in the Hubble Space Telescope that provides evidence consistent with the existence of a sensor effect analogous to the brighter-fatter effect found in Charge-Coupled Devices. We propose a model in which this effect could be understood as shifts in the effective pixel boundaries, and discuss prospects of laboratory measurements to fully characterize this effect.

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.

Unlike optical CCDs, near-infrared detectors, which are based on CMOS hybrid readout technology, typically suffer from electrical crosstalk between the pixels. The interpixel capacitance (IPC) responsible for the crosstalk affects the point-spread function (PSF) of the telescope, increasing the size and modifying the shape of all objects in the images while correlating the Poisson noise. Upcoming weak lensing surveys that use these detectors, such as WFIRST, place stringent requirements on the PSF size and shape (and the level at which these are known), which in turn must be translated into requirements on IPC. To facilitate this process, we present a first study of the effect of IPC on WFIRST PSF sizes and shapes. Realistic PSFs are forward-simulated from physical principles for each WFIRST bandpass. We explore how the PSF size and shape depends on the range of IPC coupling with pixels that are connected along an edge or corner; for the expected level of IPC in WFIRST, IPC increases the PSF sizes by $\sim$5\%. We present a linear fitting formula that describes the uncertainty in the PSF size or shape due to uncertainty in the IPC, which could arise for example due to unknown time evolution of IPC as the detectors age or due to spatial variation of IPC across the detector. We also study of the effect of a small anisotropy in the IPC, which further modifies the PSF shapes. Our results are a first, critical step in determining the hardware and characterization requirements for the detectors used in the WFIRST survey.

Weak gravitational lensing (WL) is one of the most powerful techniques to learn about the dark sector of the universe. To extract the WL signal from astronomical observations, galaxy shapes must be measured and corrected for the point spread function (PSF) of the imaging system with extreme accuracy. Future WL missions (such as the Wide-Field Infrared Survey Telescope, WFIRST) will use a family of hybrid nearinfrared CMOS detectors (HAWAII-4RG) that are untested for accurate WL measurements. Like all image sensors, these devices are subject to conversion gain nonlinearities (voltage response to collected photo-charge) that bias the shape and size of bright objects such as reference stars that are used in PSF determination. We study this type of detector nonlinearity (NL) and show how to derive requirements on it from WFIRST PSF size and ellipticity requirements. We simulate the PSF optical profiles expected for WFIRST and measure the fractional error in the PSF size and the absolute error in the PSF ellipticity as a function of star magnitude and the NL model. For our nominal NL model (a quadratic correction), we find that, uncalibrated, NL can induce an error of 0.01 (fractional size) and 0.00175 (absolute ellipticity error) in the H158 bandpass for the brightest unsaturated stars in WFIRST. In addition, our simulations show that to limit the bias of the size and ellipticity errors in the H158 band to approximately 10% of the estimated WFIRST error budget, the parameter of our quadratic NL model must be calibrated to about 1% and 2.4%, respectively. We present a fitting formula that can be used to estimate WFIRST detector NL requirements once a true PSF error budget is established.

Using a sample of 608 Type Ia supernovae from the SDSS-II and BOSS surveys, combined with a sample of foreground galaxies from SDSS-II, we estimate the weak lensing convergence for each supernova line-of-sight. We find that the correlation between this measurement and the Hubble residuals is consistent with the prediction from lensing (at a significance of 1.7sigma. Strong correlations are also found between the residuals and supernova nuisance parameters after a linear correction is applied. When these other correlations are taken into account, the lensing signal is detected at 1.4sigma. We show for the first time that distance estimates from supernovae can be improved when lensing is incorporated by including a new parameter in the SALT2 methodology for determining distance moduli. The recovered value of the new parameter is consistent with the lensing prediction. Using CMB data from WMAP7, H0 data from HST and SDSS BAO measurements, we find the best-fit value of the new lensing parameter and show that the central values and uncertainties on Omega_m and w are unaffected. The lensing of supernovae, while only seen at marginal significance in this low redshift sample, will be of vital importance for the next generation of surveys, such as DES and LSST, which will be systematics dominated.

Extremely accurate shape measurements of galaxy images are needed to probe dark energy properties with weak gravitational lensing surveys. To increase survey area with a fixed observing time and pixel count, images from surveys such as the Wide Field Infrared Survey Telescope (WFIRST) or Euclid will necessarily be undersampled and therefore distorted by aliasing. Oversampled, unaliased images can be obtained by combining multiple, dithered exposures of the same source with a suitable reconstruction algorithm. Any such reconstruction must minimally distort the reconstructed images for weak lensing analyses to be unbiased. In this paper, we use the IMage COMbination (IMCOM) algorithm of Rowe, Hirata, and Rhodes to investigate the effect of image combination on shape measurements (size and ellipticity). We simulate dithered images of sources with varying amounts of ellipticity and undersampling, reconstruct oversampled output images from them using IMCOM, and measure shape distortions in the output. Our simulations show that IMCOM creates no significant distortions when the relative offsets between dithered images are precisely known. Distortions increase with the uncertainty in those offsets but become problematic only with relatively poor astrometric precision. E.g. for images similar to those from the Astrophysics Focused Telescope Asset (AFTA) implementation of WFIRST, combining eight undersampled images (sampling ratio Q=1) with highly pessimistic uncertainty in astrometric registration (\sigma_d~10^{-3} pixels) yields an RMS shear error of O(10^{-4}). Our analysis pipeline is adapted from that of the Precision Projector Laboratory -- a joint project between NASA Jet Propulsion Laboratory and Caltech which characterizes image sensors using laboratory emulations of astronomical data.

Weak gravitational lensing observations are a key science driver for the NASA Wide Field Infrared Survey Telescope (WFIRST). To validate the performance of the WFIRST infrared detectors, we have performed a laboratory emulation of weak gravitational lensing measurements. Our experiments used a custom precision projector system to image a target mask composed of a grid of pinholes, emulating stellar point sources, onto a 1.7 micron cut-off Teledyne HgCdTe/H2RG detector. We used a 880nm LED illumination source and f/22 pupil stop to produce undersampled point spread functions similar to those expected from WFIRST. We also emulated the WFIRST image reconstruction strategy, using the IMage COMbination (IMCOM) algorithm to derive oversampled images from dithered, undersampled input images. We created shear maps for this data and computed shear correlation functions to mimic a real weak lensing analysis. After removing only 2nd order polynomial fits to the shear maps, we found that the correlation functions could be reduced to O(10^-6). This places a conservative upper limit on the detector-induced bias to the correlation function (under our test conditions). This bias is two orders of magnitude lower than the expected weak lensing signal. Restricted to scales relevant to dark energy analyses (sky separations > 0.5 arcmin), the bias is O(10^-7): comparable to the requirement for future weak lensing missions to avoid biasing cosmological parameter estimates. Our experiment will need to be upgraded and repeated under different configurations to fully characterize the shape measurement performance of WFIRST IR detectors.

Assuming General Relativity is correct on large-scales, Redshift-Space Distortions (RSDs) and the Integrated Sachs-Wolfe effect (ISW) are both sensitive to the time derivative of the linear growth function. We investigate the extent to which these probes provide complementary or redundant information when they are combined to constrain the evolution of the linear velocity power spectrum, often quantified by the function $f(z)\sigma_8(z)$, where $f$ is the logarithmic derivative of $\sigma_8$ with respect to $(1+z)$. Using a spherical Fourier-Bessel (SFB) expansion for galaxy number counts and a spherical harmonic expansion for the CMB anisotropy, we compute the covariance matrices of the signals for a large galaxy redshift survey combined with a CMB survey like Planck. The SFB basis allows accurate ISW estimates by avoiding the plane-parallel approximation, and it retains RSD information that is otherwise lost when projecting angular clustering onto redshift shells. It also allows straightforward calculations of covariance with the CMB. We find that the correlation between the ISW and RSD signals are low since the probes are sensitive to different modes. For our default surveys, on large scales ($k<0.05 \Mpc/h$), the ISW can improve constraints on $f\sigma_8$ by more than 10% compared to using RSDs alone. In the future, when precision RSD measurements are available on smaller scales, the cosmological constraints from ISW measurements will not be competitive; however, they will remain a useful consistency test for possible systematic contamination and alternative models of gravity.

One of the most pressing issues in cosmology is whether general relativity (GR) plus a dark sector is the underlying physical theory or whether a modified gravity model is needed. Upcoming dark energy experiments designed to probe dark energy with multiple methods can address this question by comparing the results of the different methods in constraining dark energy parameters. Disagreement would signal the breakdown of the assumed model (GR plus dark energy). We study the power of this consistency test by projecting constraints in the $w_0-w_a$ plane from the four different techniques of the Dark Energy Survey in the event that the underlying true model is modified gravity. We find that the standard technique of looking for overlap has some shortcomings, and we propose an alternative, more powerful Multi-dimensional Consistency Test. We introduce the methodology for projecting whether a given experiment will be able to use this test to distinguish a modified gravity model from GR.

We present an analysis of the luminosity distances of Type Ia Supernovae from the Sloan Digital Sky Survey-II (SDSS-II) Supernova Survey in conjunction with other intermediate redshift (z<0.4) cosmological measurements including redshift-space distortions from the Two-degree Field Galaxy Redshift Survey (2dFGRS), the Integrated Sachs-Wolfe (ISW) effect seen by the SDSS, and the latest Baryon Acoustic Oscillation (BAO) distance scale from both the SDSS and 2dFGRS. We have analysed the SDSS-II SN data alone using a variety of "model-independent" methods and find evidence for an accelerating universe at >97% level from this single dataset. We find good agreement between the supernova and BAO distance measurements, both consistent with a Lambda-dominated CDM cosmology, as demonstrated through an analysis of the distance duality relationship between the luminosity (d_L) and angular diameter (d_A) distance measures. We then use these data to estimate w within this restricted redshift range (z<0.4). Our most stringent result comes from the combination of all our intermediate-redshift data (SDSS-II SNe, BAO, ISW and redshift-space distortions), giving w = -0.81 +0.16 -0.18(stat) +/- 0.15(sys) and Omega_M=0.22 +0.09 -0.08 assuming a flat universe. This value of w, and associated errors, only change slightly if curvature is allowed to vary, consistent with constraints from the Cosmic Microwave Background. We also consider more limited combinations of the geometrical (SN, BAO) and dynamical (ISW, redshift-space distortions) probes.

Supermassive black hole binary systems (SMBHB) are standard sirens -- the gravitational wave analogue of standard candles -- and if discovered by gravitational wave detectors, they could be used as precise distance indicators. Unfortunately, gravitational lensing will randomly magnify SMBHB signals, seriously degrading any distance measurements. Using a weak lensing map of the SMBHB line of sight, we can estimate its magnification and thereby remove some uncertainty in its distance, a procedure we call "delensing." We find that delensing is significantly improved when galaxy shears are combined with flexion measurements, which reduce small-scale noise in reconstructed magnification maps. Under a Gaussian approximation, we estimate that delensing with a 2D mosaic image from an Extremely Large Telescope (ELT) could reduce distance errors by about 30-40% for a SMBHB at z=2. Including an additional wide shear map from a space survey telescope could reduce distance errors by 50%. Such improvement would make SMBHBs considerably more valuable as cosmological distance probes or as a fully independent check on existing probes.

The weak gravitational lensing of distant galaxies by large-scale structure is expected to become a powerful probe of dark energy. By measuring the ellipticities of large numbers of background galaxies, the subtle gravitational distortion called "cosmic shear" can be measured and used to constrain dark energy parameters. The observed galaxy ellipticities, however, are induced not by shear but by reduced shear, which also accounts for slight magnifications of the images. This distinction is negligible for present weak lensing surveys, but it will become more important as we improve our ability to measure and understand small-angle cosmic shear modes. I calculate the discrepancy between shear and reduced shear in the context of power spectra and cross spectra, finding the difference could be as high as 10% on the smallest accessible angular scales. I estimate how this difference will bias dark energy parameters obtained from two weak lensing methods: weak lensing tomography and the shear ratio method known as offset-linear scaling. For weak lensing tomography, ignoring the effects of reduced shear will cause future surveys considered by the Dark Energy Task Force to measure dark energy parameters that are biased by amounts comparable to their error bars. I advocate that reduced shear be properly accounted for in such surveys, and I provide a semi-analytic formula for doing so. Since reduced shear cross spectra do not follow an offset-linear scaling relation, the shear ratio method is similarly biased but with smaller significance.

Two methods of constraining the properties of dark energy are weak lensing tomography and cluster counting. Uncertainties in mass calibration of clusters can be reduced by using the properties of halo clustering (the clustering of clusters). However, within a single survey, weak lensing and halo clustering probe the same density fluctuations. We explore the question of whether this information can be used twice -- once in weak lensing and then again in halo clustering to calibrate cluster masses -- or whether the combined dark energy constraints are weaker than the sum of the individual constraints. For a survey like the Dark Energy Survey (DES), we find that the cosmic shearing of source galaxies at high redshifts is indeed highly correlated with halo clustering at lower redshifts. Surprisingly, this correlation does not degrade cosmological constraints for a DES-like survey, and in fact, constraints are marginally improved since the correlations themselves act as additional observables. This considerably simplifies the analysis for a DES-like survey: when weak lensing and halo clustering are treated as independent experiments, the combined dark energy constraints (cluster counts included) are accurate if not slightly conservative. Our findings mirror those of Takada and Bridle, who investigated correlations between the cosmic shear and cluster counts.

We revisit the estimation of higher order corrections to the angular power spectra of weak gravitational lensing. Extending a previous calculation of Cooray and Hu, we find two additional terms to the fourth order in potential perturbations of large-scale structure corresponding to corrections associated with the Born approximation and the neglect of line-of-sight coupling of two foreground lenses in the standard first order result. These terms alter the convergence ($\kappa\kappa$), the lensing shear E-mode ($\epsilon\epsilon$), and their cross-correlation ($\kappa\epsilon$) power spectra on large angular scales, but leave the power spectra of the lensing shear B-mode ($\beta\beta$) and rotational ($\omega\omega$) component unchanged as compared to previous estimates. The new terms complete the calculation of corrections to weak lensing angular power spectra associated with both the Born approximation and the lens-lens coupling to an order in which the contributions are most significant. Taking these features together, we find that these corrections are unimportant for any weak lensing survey, including for a full sky survey limited by cosmic variance.

Essentially all of our knowledge of the acceleration history of the Universe - including the acceleration itself - is predicated upon the validity of general relativity. Without recourse to this assumption, we use SNeIa to analyze the expansion history and find (i) very strong (5 sigma) evidence for a period of acceleration, (ii) strong evidence that the acceleration has not been constant, (iii) evidence for an earlier period of deceleration and (iv) only weak evidence that the Universe has not been decelerating since z~0.3.

Measurements of ellipticities of background galaxies are sensitive to the reduced shear, the cosmic shear divided by $(1-\kappa)$ where $\kappa$ is the projected density field. We compute the difference between shear and reduced shear both analytically and with simulations. The difference becomes more important on smaller scales, and will impact cosmological parameter estimation from upcoming experiments. A simple recipe is presented to carry out the required correction.