(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 present a comprehensive methodology for the simulation of astronomical images from optical survey telescopes. We use a photon Monte Carlo approach to construct images by sampling photons from models of astronomical source populations, and then simulating those photons through the system as they interact with the atmosphere, telescope, and camera. We demonstrate that all physical effects for optical light that determine the shapes, locations, and brightnesses of individual stars and galaxies can be accurately represented in this formalism. By using large scale grid computing, modern processors, and an efficient implementation that can produce 400,000 photons/second, we demonstrate that even very large optical surveys can be now be simulated. We demonstrate that we are able to: 1) construct kilometer scale phase screens necessary for wide-field telescopes, 2) reproduce atmospheric point-spread-function moments using a fast novel hybrid geometric/Fourier technique for non-diffraction limited telescopes, 3) accurately reproduce the expected spot diagrams for complex aspheric optical designs, and 4) recover system effective area predicted from analytic photometry integrals. This new code, the photon simulator (PhoSim), is publicly available. We have implemented the Large Synoptic Survey Telescope (LSST) design, and it can be extended to other telescopes. We expect that because of the comprehensive physics implemented in PhoSim, it will be used by the community to plan future observations, interpret detailed existing observations, and quantify systematics related to various astronomical measurements. Future development and validation by comparisons with real data will continue to improve the fidelity and usability of the code.

We have measured spatial and temporal variability in the y band sky brightness over the course of four nights above Cerro Tololo near Cerro Pachon, Chile, the planned site for the Large Synoptic Survey Telescope (LSST). Our wide-angle camera lens provided a 41 deg field of view and a 145 arcsec pixel scale. We minimized potential system throughput differences by deploying a deep depletion CCD and a filter that matches the proposed LSST y_3 band (970 nm-1030 nm). Images of the sky exhibited coherent wave structure, attributable to atmospheric gravity waves at 90 km altitude, creating 3%-4% rms spatial sky flux variability on scales of about 2 degrees and larger. Over the course of a full night the y_3 band additionally showed highly coherent temporal variability of up to a factor of 2 in flux. We estimate the mean absolute sky level to be approximately y_3 = 17.8 mag (Vega), or y_3 = 18.3 mag (AB). While our observations were made through a y_3 filter, the relative sky brightness variability should hold for all proposed y bands, whereas the absolute levels should more strongly depend on spectral response. The spatial variability presents a challenge to wide-field cameras that require illumination correction strategies that make use of stacked sky flats. The temporal variability may warrant an adaptive y band imaging strategy for LSST, to take advantage of times when the sky is darkest.

(ABRIDGED) The ACS camera on board the Hubble Space Telescope has been used to obtain deep images of the giant elliptical galaxy NGC 3610, a well-established dissipative galaxy merger remnant. These observations supersede previous WFPC2 images which revealed the presence of a population of metal-rich globular clusters (GCs) of intermediate age (~1.5-4 Gyr). We detect a total of 580 GC candidates, 46% more than from the previous WFPC2 images. The new photometry strengthens the significance of the previously found bimodality of the color distribution of GCs. Peak colors in V-I are 0.93 +/-0.01 and 1.09 +/- 0.01 for the blue and red subpopulations, respectively. The luminosity function (LF) of the inner 50% of the metal-rich (`red') population of GCs differs markedly from that of the outer 50%. In particular, the LF of the inner 50% of the red GCs shows a flattening consistent with a turnover that is about 1.0 mag fainter than the turnover of the blue GC LF. This is consistent with predictions of recent models of GC disruption for the age range mentioned above and for metallicities that are consistent with the peak color of the red GCs as predicted by population synthesis models. We determine the specific frequency of GCs in NGC 3610 and find a present-day value of S_N = 1.4 +/- 0.6. We estimate that this value will increase to S_N = 3.8 +/- 1.7 at an age of 10 Gyr, which is consistent with typical S_N values for `normal' ellipticals. Our findings constitute further evidence in support of the notion that metal-rich GC populations formed during major mergers involving gas-rich galaxies can evolve dynamically (through disruption processes) into the red, metal-rich GC populations that are ubiquitous in `normal' giant ellipticals.

We present accurate integrated-light photometry in Johnson/Cousins V, R, and I for a sample of 28 globular clusters in the Magellanic Clouds. The majority of the clusters in our sample have reliable age and metallicity estimates available in the literature. The sample encompasses ages between 50 Myr and 7 Gyr, and metallicities ([Fe/H]) between -1.5 and 0.0 dex. The sample is dominated by clusters of ages between roughly 0.5 and 2 Gyr, an age range during which the bolometric luminosity of simple stellar populations is dominated by evolved red giant branch stars and thermally pulsing asymptotic giant branch (TP-AGB) stars whose theoretical colours are rather uncertain. The VRI colours presented in this paper have been used to calibrate stellar population synthesis model predictions.

We present high-resolution near-infrared spectra, obtained with the NIRSPEC spectrograph on the W. M. Keck II Telescope, of a collection of hot, massive stars within the central 25 arcseconds of the Galactic center. We have identified a total of twenty-one emission-line stars, seven of which are new radial velocity detections with five of those being classified as He I emission-line stars for the first time. These stars fall into two categories based on their spectral properties: 1) those with narrow 2.112, 2.113 micron He I doublet absorption lines, and 2) those with broad 2.058 micron He I emission lines. These data have the highest spectral resolution ever obtained for these sources and, as a result, both components of the absorption doublet are separately resolved for the first time. We use these spectral features to measure radial velocities. The majority of the measured radial velocities have relative errors of 20 kms, smaller than those previously obtained with proper-motion or radial velocity measurements for similar stellar samples in the Galactic center. The radial velocities estimated from the He I absorption doublet are more robust than those previously estimated from the 2.058 micron emission line, since they do not suffer from confusion due to emission from the surrounding ISM. Using this velocity information, we agree that the stars are orbiting in a somewhat coherent manner but are not as defined into a disk or disks as previously thought. Finally, multi-epoch radial velocity measurements for IRS 16NE show a change in its velocity presumably due to an unseen stellar companion.

Long-slit spectra of several dozen young star clusters have been obtained at three positions in the Antennae galaxies with the Space Telescope Imaging Spectrograph (STIS) and its 52"x0.2" slit. Based on H-alpha emission-line measurements, the average cluster-to-cluster velocity dispersion in 7 different cluster aggregates ("knots") is <10 \kms. The fact that this upper limit is similar to the velocity dispersion of gas in the disks of typical spiral galaxies suggests that the triggering mechanism for the formation of young massive compact clusters ("super star clusters") is unlikely to be high velocity cloud-cloud collisions. On the other hand, models where preexisting giant molecular clouds in the disks of spiral galaxies are triggered into cluster formation are compatible with the observed low velocity dispersions. These conclusions are consistent with those reached by Zhang et al. (2001) based on comparisons between the positions of the clusters and the velocity and density structure of the nearby interstellar medium. We find evidence for systematically lower values of the line ratios [NII]/H-alpha and [SII]/H-alpha in the bright central regions of some of the knots, relative to their outer regions. This suggests that the harder ionizing photons are used up in the regions nearest the clusters, and the diffuse ionized gas farther out is photoionized by 'leakage' of the leftover low-energy photons. The low values of the [SII]/H-alpha line ratio, typically [SII]/H-alpha<0.4, indicates that the emission regions are photoionized rather than shock heated. The absence of evidence for shock-heated gas is an additional indication that high velocity cloud-cloud collisions are not playing a major role in the formation of the young clusters.

We acquired spectra of 24 LMC PNe in the 1150--3000 \AA range in order to determine carbon and other ionic abundances. The sample more than doubles the number of LMC PNe with good quality UV spectra in this wavelength range, and whose optical images are available in the {\it HST} archive. The {\it Space Telescope Imaging Spectrograph} was used with a very large aperture to obtain virtually slit-less spectra, thus the monochromatic images in the major nebulae emission lines are also available. The analysis of the data shows extremely good quality spectra. This paper presents the emission lines identified and measured, and the calculation of the ionic abundances of the emitting carbon and other ions, and total carbon abundance. P-Cygni profiles have been found in a fraction of the nebulae, and the limiting velocities of the stellar winds estimated. The total carbon abundance can be inferred reliably in most nebulae. We found that the average carbon abundance in round and elliptical PNe is one order of magnitude larger than that of the bipolar PNe, while elliptical and round PNe with a bipolar core have a bimodal behavior. This results confirm that bipolarity in LMC PNe is tightly correlated with high mass progenitors. When compared to predicted yields, we found that the observed abundance ratio show a shift toward higher carbon abundances, that may be due to initial conditions assumed in the models not appropriate for LMC PNe.

The Advanced Camera for Surveys on board the Hubble Space Telescope has been used to obtain deep high-resolution images of the giant early-type galaxy NGC 1316 which is an obvious merger remnant. These observations supersede previous, shallower observations which revealed the presence of a population of metal-rich globular clusters of intermediate age (~ 3 Gyr). We detect a total of 1496 cluster candidates, almost 4 times as many as from the previous WFPC2 images. We confirm the bimodality of the color distribution of clusters, even in V-I, with peak colors 0.93 and 1.06. The large number of detected clusters allows us to evaluate the globular cluster luminosity functions as a function of galactocentric radius. We find that the luminosity function of the inner 50% of the intermediate-age, metal-rich (`red') population of clusters differs markedly from that of the outer 50%. In particular, the luminosity function of the inner 50% of the red clusters shows a clear flattening consistent with a turnover that is about 1.0 mag fainter than the turnover of the blue clusters. This constitutes the first direct evidence that metal-rich cluster populations formed during major mergers of gas-rich galaxies can evolve dynamically (through disruption processes) into the red, metal-rich cluster populations that are ubiquitous in `normal' giant ellipticals.

We have calculated inverse-Compton Doppler factors for 3C 279 using the collection of VLBI data recently published by us, and the collection of multiwavelength spectra recently published by Hartman et al. From the Doppler factor and superluminal apparent speed, we then calculate the Lorentz factor and angle to the line-of-sight of the parsec-scale relativistic jet. We model the jet components as homogeneous spheres and the VLBI core as an unresolved inhomogeneous conical jet. The conical-jet model can be made to match both the observed X-ray emission and the VLBI properties of the core with a suitable choice of Doppler factor, implying the core makes a significant contribution to the X-ray emission. The parameters of the conical models indicate the jet is particle dominated at the radii that produce significant emission, and is not in equipartition. At the inner radius of the conical jet the magnetic field is of order 0.1 G and the relativistic-particle number density is of order 10 cm^{-3}. When all components are included in the calculation, then on average the core produces about half of the X-rays, with the other half being split between the long-lived component C4 and the brightest inner-jet component. We calculate an average speed and angle to the line-of-sight for the region of the jet interior to 1 mas of v=0.992c (gamma=8) and 4 degrees, and an average speed and angle to the line-of-sight for C4 (at a distance from the core of 3 mas) of v=0.997c (gamma=13) and 2 degrees. These values imply average Doppler factors of delta=12 for the inner jet, and delta=21 for C4.

As part of the Hubble Deep Field South program, a set of shorter 2-orbit observations were obtained of the area adjacent to the deep fields. The WFPC2 flanking fields cover a contiguous solid angle of 48 square arcminutes. Parallel observations with the STIS and NICMOS instruments produce a patchwork of additional fields with optical and near-infrared (1.6 micron) response. Deeper parallel exposures with WFPC2 and NICMOS were obtained when STIS observed the NICMOS deep field. These deeper fields are offset from the rest, and an extended low surface brightness object is visible in the deeper WFPC2 flanking field. In this data paper, which serves as an archival record of the project, we discuss the observations and data reduction, and present SExtractor source catalogs and number counts derived from the data. Number counts are broadly consistent with previous surveys from both ground and space. Among other things, these flanking field observations are useful for defining slit masks for spectroscopic follow-up over a wider area around the deep fields, for studying large-scale structure that extends beyond the deep fields, for future supernova searches, and for number counts and morphological studies, but their ultimate utility will be defined by the astronomical community.

We present and use new spectra and narrow-band images, along with previously published broad-band images, of stars in the Arches cluster to extract photometry, astrometry, equivalent width, and velocity information. The data are interpreted with a wind/atmosphere code to determine stellar temperatures, luminosities, mass-loss rates, and abundances. We have doubled the number of known emission-line stars, and we have also made the first spectroscopic identification of the main sequence for any population in the Galactic Center. We conclude that the most massive stars are bona-fide Wolf-Rayet (WR) stars and are some of the most massive stars known, having M_{initial} > 100 Msun, and prodigious winds, Mdot > 10^{-5} Msun yr^{-1}, that are enriched with helium and nitrogen; with these identifications, the Arches cluster contains about 5% of all known WR stars in the Galaxy. We find an upper limit to the velocity dispersion of 22 kms^{-1}, implying an upper limit to the cluster mass of 7(10^4) Msun within a radius of 0.23 pc; we also estimate the bulk heliocentric velocity of the cluster to be v_{cluster,odot} approximately +95 kms^{-1}.