We report studies on mitigation of optical effects of bright Low-Earth-Orbit (LEO) satellites on Vera C. Rubin Observatory and its Legacy Survey of Space and Time (LSST). These include options for pointing the telescope to avoid satellites, laboratory investigations of bright trails on the Rubin Observatory LSST Camera sensors, algorithms for correcting image artifacts caused by bright trails, experiments on darkening SpaceX Starlink satellites, and ground-based follow-up observations. Starlink satellites with no brightness mitigation are presently g ~ 5.1 mag, and an initial experiment "DarkSat" is g ~ 6.1 mag. Future Starlink darkening plans may reach g ~ 7 mag, a brightness level that enables non-linear image artifact correction to the same level as background noise. However, the satellite trails will still exist at S/N ~ 100, generating systematics that may impact data analysis and limiting some science. LEO satellite trails are wider than a point-spread function because satellites are slightly out of focus due to their finite distance; for Rubin Observatory's 8.4-m mirror and a satellite at 550 km, this widening effect is about 3 arcsec, which helps avoid saturation by decreasing the trail's peak surface brightness.

A dedicated simulator, Poisson_CCD, has been constructed which models astronomical CCDs by solving Poisson's equation numerically and simulating charge transport within the CCD. The potentials and free carrier densities within the CCD are self-consistently solved for, giving realistic results for the charge distribution within the CCD storage wells. The simulator has been used to model the CCDs which are being used to construct the LSST digital camera. The simulator output has been validated by comparing its predictions with several different types of CCD measurements, including astrometric shifts, brighter-fatter induced pixel-pixel covariances, saturation effects, and diffusion spreading. The code is open source and freely available.

We explore the weak lensing E- and B-mode shear signals of a field of galaxy clusters using both large scale structure N-body simulations and multi-color Suprime-cam & Hubble Space Telescope observations. Using the ray-traced and observed shears along with photometric redshift catalogs, we generate mass maps of the foreground overdensities by optimally filtering the tangential shear that they induce on background galaxies. We then develop and test a method to approximate the foreground structure as a superposition of NFW-like halos by locating these overdensities and determining their mass and redshift, thereby modeling the background correlated shear field as a sum of lensings induced by the foreground clusters. We demonstrate that the B-mode maps and shear correlation functions, which are generated by similarly filtering the cross shear in this method, are in agreement with observations and are related to the estimated cluster masses and locations as well as the distribution of background sources. Using the foreground mass model, we identify several sources of weak lensing B-modes including leakage and edge effects, source clustering, and multiple lensing which can be observed in deep cosmic shear surveys.

The 3.2 gigapixel LSST camera, an array of 189 thick fully-depleted CCDs, will repeatedly image the southern sky and accomplish a wide variety of science goals. However, its trove of tens of billions of object images implies stringent requirements on systematic biases imprinted during shift-and-stare CCD observation. In order to correct for these biases which, without correction, violate requirements on weak lensing precision, we investigate CCD systematics using both simulations of charge transport as well as with a unique bench-top optical system matched to the LSST's fast f/1.2 beam. By illuminating single CCDs with realistic scenes of stars and galaxies and then analyzing these images with the LSST data management pipelines, we can characterize the survey's imaging performance well before the camera's first light. We present measurements of several CCD systematics under varying conditions in the laboratory, including the brightness-dependent broadening of star and galaxy images, charge transport anomalies in the silicon bulk as well as the edges, and serial deferred charge. Alongside these measurements, we also present the development and testing of physics-based models which inform corrections or mitigation strategies for these systematics. Optimization of the CCD survey operation under a variety of realistic observational conditions, including systematic effects from the optics, clocking, sky brightness, and image analysis, will be critical to achieve the LSST's goals of precision astronomy and cosmology.

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

Reduction of images and science analysis from ground-based telescopes such as the LSST requires detailed knowledge of the PSF of the image, which includes components attributable to the instrument as well as components attributable to the atmosphere. Because the atmospheric component is constantly changing, the PSF is typically extracted from each image by measuring the size and shape of star images across the CCD, then building a fitting function over the focal plane which is used to model the PSF for analysis of extended sources such as galaxies. Since the stars in each CCD field have a range of brightnesses, accurate knowledge of the PSF for point sources of varying brightness is essential. It has been found that in thick, fully-depleted CCDs, the electrostatic repulsion of charge stored in the collecting wells gives rise to a larger and slightly more elliptical PSF for brighter stars. This "brighter-fatter" effect has been reported in some detail in the literature. In this work, we report direct and indirect measurements of this effect in prototype LSST sensors, and describe a detailed physics-based model of the electrostatics and charge transport within the CCD.

We characterize the astrometric distortion at the edges of thick, fully-depleted CCDs in the lab using a bench-top simulation of LSST observing. By illuminating an array of forty thousand pinholes (30mu m diameter) at the object plane of a f/1.2 optical reimager, thousands of PSFs can be imaged over a 4Kx4K pixel CCD. Each high purity silicon pixel, 10mu m square by 100mu m deep, can then be individually characterized through a series of sub-pixel dithers in the X/Y plane. The unique character [response, position, shape] of each pixel as a function of flux, wavelength, back side bias, etc. can be investigated. We measure the magnitude and onset of astrometric error at the edges of the detector as a test of the experimental setup, using a LSST prototype CCD. We show that this astrometric error at the edge is sourced from non-uniformities in the electric field lines that define pixel boundaries. This edge distortion must be corrected in order to optimize the science output of weak gravitational lensing and large scale structure measurements for the LSST.

We describe a camera beam simulator for the LSST which is capable of illuminating a 60mm field at f/1.2 with realistic astronomical scenes, enabling studies of CCD astrometric and photometric performance. The goal is to fully simulate LSST observing, in order to characterize charge transport and other features in the thick fully depleted CCDs and to probe low level systematics under realistic conditions. The automated system simulates the centrally obscured LSST beam and sky scenes, including the spectral shape of the night sky. The doubly telecentric design uses a nearly unit magnification design consisting of a spherical mirror, three BK7 lenses, and one beam-splitter window. To achieve the relatively large field the beam-splitter window is used twice. The motivation for this LSST beam test facility was driven by the need to fully characterize a new generation of thick fully-depleted CCDs, and assess their suitability for the broad range of science which is planned for LSST. Due to the fast beam illumination and the thick silicon design [each pixel is 10 microns wide and over 100 microns deep] at long wavelengths there can be effects of photon transport and charge transport in the high purity silicon. The focal surface covers a field more than sufficient for a 40x40 mm LSST CCD. Delivered optical quality meets design goals, with 50% energy within a 5 micron circle. The tests of CCD performance are briefly described.

We are conducting a long-term photometric survey of the nearby galaxy M33 to discover Cepheids, eclipsing binaries, and long-period variables. The dataset combines previously-obtained optical images from the DIRECT project with new observations acquired at the WIYN 3.5m telescope. The entire data set spans over 7 years with excellent synoptic coverage which will enable the discovery and characterization of stars displaying variability over a wide range of timescales (days, weeks, months, years). In this preliminary work we show representative light curves of different variables we found so far in two fields, color-magnitude diagrams, and optical Cepheid Period-Luminosity relations for M33. The ultimate goal of the project is to provide an absolute calibration of the Cepheid Period-Luminosity relation, and to study its metallicity dependence at optical wavelengths.