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 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)

A survey that can cover the sky in optical bands over wide fields to faint magnitudes with a fast cadence will enable many of the exciting science opportunities of the next decade. The Large Synoptic Survey Telescope (LSST) will have an effective aperture of 6.7 meters and an imaging camera with field of view of 9.6 deg^2, and will be devoted to a ten-year imaging survey over 20,000 deg^2 south of +15 deg. Each pointing will be imaged 2000 times with fifteen second exposures in six broad bands from 0.35 to 1.1 microns, to a total point-source depth of r~27.5. The LSST Science Book describes the basic parameters of the LSST hardware, software, and observing plans. The book discusses educational and outreach opportunities, then goes on to describe a broad range of science that LSST will revolutionize: mapping the inner and outer Solar System, stellar populations in the Milky Way and nearby galaxies, the structure of the Milky Way disk and halo and other objects in the Local Volume, transient and variable objects both at low and high redshift, and the properties of normal and active galaxies at low and high redshift. It then turns to far-field cosmological topics, exploring properties of supernovae to z~1, strong and weak lensing, the large-scale distribution of galaxies and baryon oscillations, and how these different probes may be combined to constrain cosmological models and the physics of dark energy.

We analyze warps in the nearby edge-on spiral galaxies observed in the {\em Spitzer/IRAC} 4.5 micron band. In our sample of 24 galaxies we find evidence of warp in 14 galaxies. We estimate the observed onset radii for the warps in a subsample of 10 galaxies. The dark matter distribution in each of these galaxies are calculated using the mass distribution derived from the observed light distribution and the observed rotation curves. The theoretical predictions of the onset radii for the warps are then derived by applying a self-consistent linear response theory to the obtained mass models for 6 galaxies with rotation curves in the literature. By comparing the observed onset radii to the theoretical ones, we find that discs with constant thickness can not explain the observations; moderately flaring discs are needed. The required flaring is consistent with the observations. Our analysis shows that the onset of warp is not symmetric in our sample of galaxies. We define a new quantity called the onset-asymmetry index and study its dependence on galaxy properties. The onset asymmetries in warps tend to be larger in galaxies with smaller disc scale lengths. We also define and quantify the global asymmetry in the stellar light distribution, that we call the edge-on asymmetry in edge-on galaxies. It is shown that in most cases the onset asymmetry in warp is actually anti-correlated with the measured edge-on asymmetry in our sample of edge-on galaxies and this could plausibly indicate that the surrounding dark matter distribution is asymmetric.

We present a summary of our recent work on the vertical distribution of the resolved stellar populations in six low mass, edge-on, spiral galaxies observed with the Hubble Space Telescope Advanced Camera for Surveys (HST/ACS). In each galaxy we find evidence for an extraplanar stellar component extending up to 15 scale heights (~3.5 kpc) above the plane. We analyze the vertical distribution as a function of stellar age by tracking changes in the color-magnitude diagram. The young stellar component (< 10^8 yrs) is found to have a scale height larger than the young component in the Milky Way, suggesting that stars in these low mass galaxies form in a thicker disk. We also find that the scale height of a stellar population increases with age, with young main sequence stars, intermediate age asymptotic giant branch stars, and old red giant branch stars having succesively larger scale heights in each galaxy. The color distributions of the red giant branch (RGB) stars well above the midplane indicate that the extended stellar components we see are moderately metal-poor, with peak metallicities of [Fe/H] ~ -1 and with little or no metallicity gradient with height. The lack of metallicity gradient can be explained if a majority of extraplanar RGB stars were formed at early times and are not dominated by a younger heated population. Our observations suggest that, like the Milky Way, low mass disk galaxies also have multiple stellar components. We examine our results in light of disk heating and merger scenarios and conclude that both mechanisms likely played a role in forming the disks of our sample galaxies.

This paper addresses available constraints on mass models fitted to rotation curves. Mass models of disk galaxies have well-known degeneracies, that prevent a unique mass decomposition. The most notable is due to the unknown value of the stellar mass-to-light ratio (the disk-halo degeneracy); even with this known, degeneracies between the halo parameters themselves may prevent an unambiguous determination of the shape of the dark halo profile, which includes the inner density slope of the dark matter halo. The latter is often referred to as the ``cusp-core degeneracy''. We explore constraints on the disk and halo parameters and apply these to four mock and six observed disk galaxies with high resolution and extended rotation curves. Our full set of constraints consists of mass-to-light (M/L) ratios from stellar population synthesis models based on B-R colors, constraints on halo parameters from N-body simulations, and constraining the halo virial velocity to be less than the maximum observed velocity. These constraints are only partially successful in lifting the cusp-core degeneracy. The effect of adiabatic contraction of the halo by the disk is to steepen cores into cusps and reduce the best-fit halo concentration and M/L values (often significantly). We also discuss the effect of disk thickness, halo flattening, distance errors, and rotation curve error values on mass modeling. Increasing the imposed minimum rotation curve error from typically low, underestimated values to more realistic estimates decreases the chi-square substantially and makes distinguishing between a cuspy or cored halo profile even more difficult. In spite of the degeneracies and uncertainties present, our constrained mass modeling favors sub-maximal disks (i.e., a dominant halo) at 2.2 disk scale lengths, with Vdisk/Vtot < 0.6.