LSST has chosen Kubernetes as the platform for deploying and operating the LSST Science Platform. We first present the background reasoning behind this decision, including both instrument-agnostic as well as LSST-specific requirements. We then discuss the basic principles of Kubernetes and Helm, and how they are used as the deployment base for the LSST Science Platform. Furthermore, we provide an example of how an external group may use these publicly available software resources to deploy their own instance of the LSST Science Platform, and customize it to their needs. Finally, we discuss how more astronomy software can follow these patterns to gain similar benefits.

(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 Astropy project supports and fosters the development of open-source and openly-developed Python packages that provide commonly-needed functionality to the astronomical community. A key element of the Astropy project is the core package Astropy, which serves as the foundation for more specialized projects and packages. In this article, we provide an overview of the organization of the Astropy project and summarize key features in the core package as of the recent major release, version 2.0. We then describe the project infrastructure designed to facilitate and support development for a broader ecosystem of inter-operable packages. We conclude with a future outlook of planned new features and directions for the broader Astropy project.

We present high-resolution Giant Metrewave Radio Telescope (GMRT) HI observations and deep Canada-France-Hawaii Telescope (CFHT) optical imaging of two galaxy groups: NGC 4725/47 and NGC 3166/9. These data are part of a multi-wavelength unbiased survey of the gas-rich dwarf galaxy populations in three nearby interacting galaxy groups. The NGC 4725/47 group hosts two tidal knots and one dIrr. Both tidal knots are located within a prominent HI tidal tail, appear to have sufficient mass (M_gas~10^8 M_sol) to evolve into long-lived tidal dwarf galaxies (TDGs) and are fairly young in age. The NGC 3166/9 group contains a TDG candidate, AGC 208457, at least three dIrrs and four HI knots. Deep CFHT imaging confirms that the optical component of AGC 208457 is bluer -- with a 0.28 mag g-r colour -- and a few Gyr younger than its purported parent galaxies. Combining the results for these groups with those from the NGC 871/6/7 group reported earlier, we find that the HI properties, estimated stellar ages and baryonic content of the gas-rich dwarfs clearly distinguish tidal features from their classical counterparts. We optimistically identify four potentially long-lived tidal objects associated to three separate pairs of interacting galaxies, implying that TDGs are not readily produced during interaction events as suggested by some recent simulations. The tidal objects examined in this survey also appear to have a wider variety of properties than TDGs of similar mass formed in current simulations of interacting galaxies, which could be the result of pre- or post-formation environmental influences.

The Large Synoptic Survey Telescope (LSST) is a large-aperture, wide-field, ground-based survey system that will image the sky in six optical bands from 320 to 1050 nm, uniformly covering approximately $18,000$deg$^2$ of the sky over 800 times. The LSST is currently under construction on Cerro Pach\'on in Chile, and expected to enter operations in 2022. Once operational, the LSST will explore a wide range of astrophysical questions, from discovering "killer" asteroids to examining the nature of Dark Energy. The LSST will generate on average 15 TB of data per night, and will require a comprehensive Data Management system to reduce the raw data to scientifically useful catalogs and images with minimum human intervention. These reductions will result in a real-time alert stream, and eleven data releases over the 10-year duration of LSST operations. To enable this processing, the LSST project is developing a new, general-purpose, high-performance, scalable, well documented, open source data processing software stack for O/IR surveys. Prototypes of this stack are already capable of processing data from existing cameras (e.g., SDSS, DECam, MegaCam), and form the basis of the Hyper-Suprime Cam (HSC) Survey data reduction pipeline.

Our proximity and external vantage point make M31 an ideal testbed for understanding the structure of spiral galaxies. The Andromeda Optical and Infrared Disk Survey (ANDROIDS) has mapped M31's bulge and disk out to R=40 kpc in $ugriJK_s$ bands with CFHT using a careful sky calibration. We use Bayesian modelling of the optical-infrared spectral energy distribution (SED) to estimate profiles of M31's stellar populations and mass along the major axis. This analysis provides evidence for inside-out disk formation and a declining metallicity gradient. M31's $i$-band mass-to-light ratio ($M/L_i^*$) decreases from 0.5 dex in the bulge to $\sim 0.2$ dex at 40 kpc. The best-constrained stellar population models use the full $ugriJK_s$ SED but are also consistent with optical-only fits. Therefore, while NIR data can be successfully modelled with modern stellar population synthesis, NIR data do not provide additional constraints in this application. Fits to the $gi$-SED alone yield $M/L_i^*$ that are systematically lower than the full SED fit by 0.1 dex. This is still smaller than the 0.3 dex scatter amongst different relations for $M/L_i$ via $g-i$ colour found in the literature. We advocate a stellar mass of $M_*(30\mathrm{kpc})=10.3^{+2.3}_{-1.7}\times 10^{10}\mathrm{M}_\odot$ for the M31 bulge and disk.

We present GMRT HI observations and deep CFHT MegaCam optical images of the gas-rich interacting galaxy group NGC 871/NGC 876/NGC 877 (hereafter NGC 871/6/7). Our high-resolution data sets provide a census of the HI and stellar properties of the detected gas-rich group members. In addition to a handful of spiral, irregular and dwarf galaxies, this group harbours an intriguing HI feature, AGC 749170, that has a gas mass of ~10^9.3 M_sol, a dynamical-to-gas mass ratio of ~1 (assuming the cloud is rotating and in dynamical equilibrium) and no optical counterpart in previous imaging. Our observations have revealed a faint feature in the CFHT g'- and r'-bands; if it is physically associated with AGC 749170, the latter has M/L_g > 1000 M_sol/L_sol as well as a higher metallicity (estimated using photometric colours) and a significantly younger stellar population than the other low-mass gas-rich group members. These properties, as well as its spectral and spatial location with respect to its suspected parent galaxies, strongly indicate a tidal origin for AGC 749170. Overall, the HI properties of AGC 749170 resemble those of other optically dark/dim clouds that have been found in groups. These clouds could represent a class of relatively long-lived HI-rich tidal remnants that survive in intermediate-density environments.

The Andromeda Optical and Infrared Disk Survey has mapped M31 in $u^* g^\prime r^\prime i^\prime J K_s$ wavelengths out to R=40 kpc using the MegaCam and WIRCam wide-field cameras on the Canada-France-Hawaii Telescope. Our survey is uniquely designed to simultaneously resolve stars while also carefully reproducing the surface brightness of M31, allowing us to study M31's global structure in the context of both resolved stellar populations and spectral energy distributions. We use the Elixir-LSB method to calibrate the optical $u^* g^\prime r^\prime i^\prime$ images by building real-time maps of the sky background with sky-target nodding. These maps are stable to $\mu_g \lesssim 28.5$ mag arcsec$^{-2}$ and reveal warps in the outer M31 disk in surface brightness. The equivalent WIRCam mapping in the near-infrared uses a combination of sky-target nodding and image-to-image sky offset optimization to produce stable surface brightnesses. This study enables a detailed analysis of the systematics of spectral energy distribution fitting with near-infrared bands where asymptotic giant branch stars impose a significant, but ill-constrained, contribution to the near-infrared light of a galaxy. Here we present panchromatic surface brightness maps and initial results from our near-infrared resolved stellar catalog.

We present wide-field near-infrared J and Ks images of the Andromeda Galaxy taken with WIRCam on the Canada-France-Hawaii Telescope (CFHT) as part of the Andromeda Optical and Infrared Disk Survey (ANDROIDS). This data set allows simultaneous observations of resolved stars and NIR surface brightness across M31's entire bulge and disk (within R=22 kpc). The primary concern of this work is the development of NIR observation and reduction methods to recover a uniform surface brightness map across the 3x1 degree disk of M31. This necessitates sky-target nodding across 27 WIRCam fields. Two sky-target nodding strategies were tested, and we find that strictly minimizing sky sampling latency does not maximize sky subtraction accuracy, which is at best 2% of the sky level. The mean surface brightness difference between blocks in our mosaic can be reduced from 1% to 0.1% of the sky brightness by introducing scalar sky offsets to each image. The true surface brightness of M31 can be known to within a statistical zeropoint of 0.15% of the sky level (0.2 mag arcsec sq. uncertainty at R=15 kpc). Surface brightness stability across individual WIRCam frames is limited by both WIRCam flat field evolution and residual sky background shapes. To overcome flat field variability of order 1% over 30 minutes, we find that WIRCam data should be calibrated with real-time sky flats. Due either to atmospheric or instrumental variations, the individual WIRCam frames have typical residual shapes with amplitudes of 0.2% of the sky after real-time flat fielding and median sky subtraction. We present our WIRCam reduction pipeline and performance analysis here as a template for future near-infrared observers needing wide-area surface brightness maps with sky-target nodding, and give specific recommendations for improving photometry of all CFHT/WIRCam programs. (Abridged)

Our GHOSTS survey measures the stellar envelope properties of 14 nearby disk galaxies by imaging their resolved stellar populations with HST/ACS&WFPC2. Most of the massive galaxies in the sample (Vrot >200 km/s) have very extended stellar envelopes with Sersic law profiles or mu(r)~r^-2.5 power law profiles in the outer regions. For these massive galaxies we can fit the central bulge light and the outer halo out to 30 kpc with one and the same Sersic profile and the stellar surface density of the profiles correlate with Hubble type and bulge- to-disk ratio. This suggests that the central bulges and inner halos are created in the same process. Smaller galaxies (Vrot ~100 km/s) have much smaller stellar envelopes, but depending on geometry, they could still be more luminous than expected from satellite remnants in hierarchical galaxy formation models. Alternatively, they could be created by disk heating through the bombardment of small dark matter sub-halos. The halos we fit are highly flattened, with minor- over-ma jor axis ratios c/a=0.6. The halos are somewhat more compact than hierarchical model predictions. The halos show small metallicity gradient out to 30 kpc and the massive galaxies have typical [Fe/H]~-0.8. We find indications of halo substructure in many galaxies, but some halos seem remarkable smooth.

In hierarchical galaxy formation the stellar halos of galaxies are formed by the accretion of minor satellites and therefore contain valuable information about the (early) assembly process of galaxies. Our GHOSTS survey measures the stellar envelope properties of 14 nearby disk galaxies by imaging their resolved stellar populations with HST/ACS&WFPC2. Most of the massive galaxies in the sample (Vrot>200 km/s) have very extended stellar envelopes with mu(r) ~ r^{-2.5} power law profiles in the outer regions. For these massive galaxies there is some evidence that the stellar surface density of the profiles correlates with Hubble type and bulge-to-disk ratio, begging the question whether these envelopes are more related to bulges than to a Milky Way-type stellar halo. Smaller galaxies (Vrot 100 km/s) have much smaller stellar envelopes, but depending on geometry, they could still be more luminous than expected from satellite remnants in hierarchical galaxy formation models. Alternatively, they could be created by disk heating through the bombardment of small dark matter sub-halos. We find that galaxies show varying amounts of halo substructure.

We use HST/ACS to study the resolved stellar populations of the nearby, nearly edge-on galaxy NGC4244 across its outer disk surface density break. The stellar photometry allows us to study the distribution of different stellar populations and reach very low equivalent surface brightnesses. We find that the break occurs at the same radius for young, intermediate age, and old stars. The stellar density beyond the break drops sharply by a factor of at least 600 in 5 kpc. The break occurs at the same radius independent of height above the disk, but is sharpest in the midplane and nearly disappears at large heights. These results make it unlikely that truncations are caused by a star formation threshold alone: the threshold would have to keep the same radial position from less than 100 Myr to 10 Gyr ago, in spite of potential disturbances such as infall and redistribution of gas by internal processes. A dynamical interpretation of truncation formation is more likely such as due to angular momentum redistribution by bars or density waves, or heating and stripping of stars caused by the bombardment of dark matter sub-halos. The latter explanation is also in quantitative agreement with the small diffuse component we see around the galaxy.