The Simons Observatory (SO) is a ground-based cosmic microwave background (CMB) experiment sited on Cerro Toco in the Atacama Desert in Chile that promises to provide breakthrough discoveries in fundamental physics, cosmology, and astrophysics. Supported by the Simons Foundation, the Heising-Simons Foundation, and with contributions from collaborating institutions, SO will see first light in 2021 and start a five year survey in 2022. SO has 287 collaborators from 12 countries and 53 institutions, including 85 students and 90 postdocs. The SO experiment in its currently funded form ('SO-Nominal') consists of three 0.4 m Small Aperture Telescopes (SATs) and one 6 m Large Aperture Telescope (LAT). Optimized for minimizing systematic errors in polarization measurements at large angular scales, the SATs will perform a deep, degree-scale survey of 10% of the sky to search for the signature of primordial gravitational waves. The LAT will survey 40% of the sky with arc-minute resolution. These observations will measure (or limit) the sum of neutrino masses, search for light relics, measure the early behavior of Dark Energy, and refine our understanding of the intergalactic medium, clusters and the role of feedback in galaxy formation. With up to ten times the sensitivity and five times the angular resolution of the Planck satellite, and roughly an order of magnitude increase in mapping speed over currently operating ("Stage 3") experiments, SO will measure the CMB temperature and polarization fluctuations to exquisite precision in six frequency bands from 27 to 280 GHz. SO will rapidly advance CMB science while informing the design of future observatories such as CMB-S4.

The Simons Observatory (SO) is a new cosmic microwave background experiment being built on Cerro Toco in Chile, due to begin observations in the early 2020s. We describe the scientific goals of the experiment, motivate the design, and forecast its performance. SO will measure the temperature and polarization anisotropy of the cosmic microwave background in six frequency bands: 27, 39, 93, 145, 225 and 280 GHz. The initial configuration of SO will have three small-aperture 0.5-m telescopes (SATs) and one large-aperture 6-m telescope (LAT), with a total of 60,000 cryogenic bolometers. Our key science goals are to characterize the primordial perturbations, measure the number of relativistic species and the mass of neutrinos, test for deviations from a cosmological constant, improve our understanding of galaxy evolution, and constrain the duration of reionization. The SATs will target the largest angular scales observable from Chile, mapping ~10% of the sky to a white noise level of 2 $\mu$K-arcmin in combined 93 and 145 GHz bands, to measure the primordial tensor-to-scalar ratio, $r$, at a target level of $\sigma(r)=0.003$. The LAT will map ~40% of the sky at arcminute angular resolution to an expected white noise level of 6 $\mu$K-arcmin in combined 93 and 145 GHz bands, overlapping with the majority of the LSST sky region and partially with DESI. With up to an order of magnitude lower polarization noise than maps from the Planck satellite, the high-resolution sky maps will constrain cosmological parameters derived from the damping tail, gravitational lensing of the microwave background, the primordial bispectrum, and the thermal and kinematic Sunyaev-Zel'dovich effects, and will aid in delensing the large-angle polarization signal to measure the tensor-to-scalar ratio. The survey will also provide a legacy catalog of 16,000 galaxy clusters and more than 20,000 extragalactic sources.

In recent years, machine learning (ML) methods have remarkably improved how cosmologists can interpret data. The next decade will bring new opportunities for data-driven cosmological discovery, but will also present new challenges for adopting ML methodologies and understanding the results. ML could transform our field, but this transformation will require the astronomy community to both foster and promote interdisciplinary research endeavors.

We train a deep residual convolutional neural network (CNN) to predict the gas-phase metallicity ($Z$) of galaxies derived from spectroscopic information ($Z \equiv 12 + \log(\rm O/H)$) using only three-band $gri$ images from the Sloan Digital Sky Survey. When trained and tested on $128 \times 128$-pixel images, the root mean squared error (RMSE) of $Z_{\rm pred} - Z_{\rm true}$ is only 0.085 dex, vastly outperforming a trained random forest algorithm on the same data set (RMSE $=0.130$ dex). The amount of scatter in $Z_{\rm pred} - Z_{\rm true}$ decreases with increasing image resolution in an intuitive manner. We are able to use CNN-predicted $Z_{\rm pred}$ and independently measured stellar masses to recover a mass-metallicity relation with $0.10$ dex scatter. Because our predicted MZR shows no more scatter than the empirical MZR, the difference between $Z_{\rm pred}$ and $Z_{\rm true}$ can not be due to purely random error. This suggests that the CNN has learned a representation of the gas-phase metallicity, from the optical imaging, beyond what is accessible with oxygen spectral lines.

We report on newly identified galaxy clusters from the high signal-to-noise ($>5$-$\sigma$) end of the second all-sky \textit{Planck} Sunyaev-Zel'dovich (SZ) catalog (PSZ2). The clusters are identified in deep, optical imaging from the Kitt Peak National Observatory 4m Mayall telescope taken between 2014 and 2017. Here we focus on the highest richness systems, and identify galaxy clusters through a combination of the maxBCG algorithm and visual image inspection. Galaxy clusters are considered to be confirmed if they are both rich and spatially coincident ($\lesssim 6'$) with the reported PSZ2 position. Of the 85 fields containing unconfirmed PSZ2 candidates observed, we find 15 (17.6% of the observed sample) corresponding galaxy clusters ($ 0.13 < z < 0.78$), 12 of which are previously unrecognized as counterparts. To explain this low identification fraction, we consider three possible scenarios: that clusters are (1) mostly at low-z, (2) mostly at high-z, or (3) located in fields with high object density. None of these scenarios alone can account for the low purity of rich galaxy clusters among the high signal-to-noise PSZ2 unconfirmed candidates.

Meeting the science goals for many current and future ground-based optical large-area sky surveys requires that the calibrated broadband photometry is stable in time and uniform over the sky to 1% precision or better. Past surveys have achieved photometric precision of 1-2% by calibrating the survey's stellar photometry with repeated measurements of a large number of stars observed in multiple epochs. The calibration techniques employed by these surveys only consider the relative frame-by-frame photometric zeropoint offset and the focal plane position-dependent illumination corrections, which are independent of the source color. However, variations in the wavelength dependence of the atmospheric transmission and the instrumental throughput induce source color-dependent systematic errors. These systematic errors must also be considered to achieve the most precise photometric measurements. In this paper, we examine such systematic chromatic errors using photometry from the Dark Energy Survey (DES) as an example. We define a natural magnitude system for DES and calculate the systematic errors on stellar magnitudes, when the atmospheric transmission and instrumental throughput deviate from the natural system. We conclude that the systematic chromatic errors caused by the change of airmass in each exposure, the change of the precipitable water vapor and aerosol in the atmosphere over time, and the non-uniformity of instrumental throughput over the focal plane, can be up to 2% in some bandpasses. We compare the calculated systematic chromatic errors with the observed DES data. For the test sample data, we correct these errors using measurements of the atmospheric transmission and instrumental throughput. The residual after correction is less than 0.3%. We also find that the errors for non-stellar objects are redshift-dependent and can be larger than those for stars at certain redshifts.

We use data from the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey to study how the spatial variation in the stellar populations of galaxies relate to the formation of galaxies at $1.5 < z < 3.5$. We use the Internal Color Dispersion (ICD), measured between the rest-frame UV and optical bands, which is sensitive to age (and dust attenuation) variations in stellar populations. The ICD shows a relation with the stellar masses and morphologies of the galaxies. Galaxies with the largest variation in their stellar populations as evidenced by high ICD have disk-dominated morphologies (with S\'{e}rsic indexes $< 2$) and stellar masses between $10 < \mathrm{Log~M/ M_\odot}< 11$. There is a marked decrease in the ICD as the stellar mass and/or the S\'ersic index increases. By studying the relations between the ICD and other galaxy properties including sizes, total colors, star-formation rate, and dust attenuation, we conclude that the largest variations in stellar populations occur in galaxies where the light from newly, high star-forming clumps contrasts older stellar disk populations. This phase reaches a peak for galaxies only with a specific stellar mass range, $10 < \mathrm{Log~M/ M_\odot} < 11$, and prior to the formation of a substantial bulge/spheroid. In contrast, galaxies at higher or lower stellar masses, and/or higher S\'{e}rsic index ($n > 2$) show reduced ICD values, implying a greater homogeneity of their stellar populations. This indicates that if a galaxy is to have both a quiescent bulge along with a star forming disk, typical of Hubble Sequence galaxies, this is most common for stellar masses $10 < \mathrm{Log~M/M_\odot} < 11$ and when the bulge component remains relatively small ($n<2$).

We have built an Atmospheric Transmission Monitoring Camera (aTmCam), which consists of four telescopes and detectors each with a narrow-band filter that monitors the brightness of suitable standard stars. Each narrowband filter is selected to monitor a different wavelength region of the atmospheric transmission, including regions dominated by the precipitable water vapor and aerosol optical depth. The colors of the stars are measured by this multi narrow-band imager system simultaneously. The measured colors, a model of the observed star, and the measured throughput of the system can be used to derive the atmospheric transmission of a site on sub-minute time scales. We deployed such a system to the Cerro Tololo Inter-American Observatory (CTIO) and executed two one-month-long observing campaigns in Oct-Nov 2012 and Sept-Oct 2013. We have determined the time and angular scales of variations in the atmospheric transmission above CTIO during these observing runs. We also compared our results with those from a GPS Water Vapor Monitoring System and find general agreement. The information for the atmospheric transmission can be used to improve photometric precision of large imaging surveys such as the Dark Energy Survey and the Large Synoptic Survey Telescope.

Traditional color and airmass corrections can typically achieve ~0.02 mag precision in photometric observing conditions. A major limiting factor is the variability in atmospheric throughput, which changes on timescales of less than a night. We present preliminary results for a system to monitor the throughput of the atmosphere, which should enable photometric precision when coupled to more traditional techniques of less than 1% in photometric conditions. The system, aTmCam, consists of a set of imagers each with a narrow-band filter that monitors the brightness of suitable standard stars. Each narrowband filter is selected to monitor a different wavelength region of the atmospheric transmission, including regions dominated by the precipitable water absorption and aerosol scattering. We have built a prototype system to test the notion that an atmospheric model derived from a few color indices measurements can be an accurate representation of the true atmospheric transmission. We have measured the atmospheric transmission with both narrowband photometric measurements and spec- troscopic measurements; we show that the narrowband imaging approach can predict the changes in the throughput of the atmosphere to better than ~10% across a broad wavelength range, so as to achieve photometric precision less than 0.01 mag.