Spatially resolved kinematics of nearby galaxies has shown that the ratio of dynamical- to stellar population-based estimates of the mass of a galaxy ($M_*^{\rm JAM}/M_*$) correlates with $\sigma_e$, if $M_*$ is estimated using the same IMF for all galaxies and the stellar M/L ratio within each galaxy is constant. This correlation may indicate that, in fact, the IMF is more dwarf-rich for galaxies with large $\sigma$. We use this correlation to estimate a dynamical or IMF-corrected stellar mass, $M_*^{\rm \alpha_{JAM}}$, from $M_{*}$ and $\sigma_e$ for a sample of $6 \times 10^5$ SDSS galaxies for which spatially resolved kinematics is not available. We also compute the `virial' mass estimate $k(n,R)\,R_e\,\sigma_R^2/G$, where $n$ is the Sersic index, in the SDSS and ATLAS$^{\rm 3D}$ samples. We show that an $n$-dependent correction must be applied to the $k(n,R)$ values provided by Prugniel & Simien (1997). Our analysis also shows that the shape of the velocity dispersion profile in the ATLAS$^{\rm 3D}$ sample varies weakly with $n$: $(\sigma_R/\sigma_e) = (R/R_e)^{-\gamma(n)}$. The resulting stellar mass functions, based on $M_*^{\rm \alpha_{JAM}}$ and the recalibrated virial mass, are in good agreement. If the $M_*^{\rm \alpha_{JAM}}/M_* - \sigma_e$ correlation is indeed due to the IMF, and stellar M/L gradients can be ignored, then our $\phi(M_*^{\rm \alpha_{JAM}})$ is an estimate of the stellar mass function in which $\sigma_e$-dependent variations in the IMF across the population have been accounted for. Using a Fundamental Plane based observational proxy for $\sigma_e$ produces comparable results. By demonstrating that cheaper proxies are sufficiently accurate, our analysis should enable a more reliable census of the mass in stars for large galaxy samples, at a fraction of the cost. Our results are provided in tabular form.

The Sloan Digital Sky Survey pipeline photometry underestimates the brightnesses of the most luminous galaxies. This is mainly because (i) the SDSS overestimates the sky background and (ii) single or two-component Sersic-based models better fit the surface brightness profile of galaxies, especially at high luminosities, than does the de Vaucouleurs model used by the SDSS pipeline. We use the PyMorph photometric reductions to isolate effect (ii) and show that it is the same in the full sample as in small group environments, and for satellites in the most massive clusters as well. None of these are expected to be significantly affected by intracluster light (ICL). We only see an additional effect for centrals in the most massive halos, but we argue that even this is not dominated by ICL. Hence, for the vast majority of galaxies, the differences between PyMorph and SDSS pipeline photometry cannot be ascribed to the semantics of whether or not one includes the ICL when describing the stellar mass of massive galaxies. Rather, they likely reflect differences in star formation or assembly histories. Failure to account for the SDSS underestimate has significantly biased most previous estimates of the SDSS luminosity and stellar mass functions, and therefore Halo Model estimates of the z ~ 0.1 relation between the mass of a halo and that of the galaxy at its center. We also show that when one studies correlations, at fixed group mass, with a quantity which was not used to define the groups, then selection effects appear. We show why such effects arise, and should not be mistaken for physical effects.

We quantify the systematic effects on the stellar mass function which arise from assumptions about the stellar population, as well as how one fits the light profiles of the most luminous galaxies at z ~ 0.1. When comparing results from the literature, we are careful to separate out these effects. Our analysis shows that while systematics in the estimated comoving number density which arise from different treatments of the stellar population remain of order < 0.5 dex, systematics in photometry are now about 0.1 dex, despite recent claims in the literature. Compared to these more recent analyses, previous work based on Sloan Digital Sky Survey (SDSS) pipeline photometry leads to underestimates of rho_*(> M_*) by factors of 3-10 in the mass range 10^11 - 10^11.6 M_Sun, but up to a factor of 100 at higher stellar masses. This impacts studies which match massive galaxies to dark matter halos. Although systematics which arise from different treatments of the stellar population remain of order < 0.5 dex, our finding that systematics in photometry now amount to only about 0.1 dex in the stellar mass density is a significant improvement with respect to a decade ago. Our results highlight the importance of using the same stellar population and photometric models whenever low and high redshift samples are compared.

A number of recent estimates of the total luminosities of galaxies in the SDSS are significantly larger than those reported by the SDSS pipeline. This is because of a combination of three effects: one is simply a matter of defining the scale out to which one integrates the fit when defining the total luminosity, and amounts on average to < 0.1 mags even for the most luminous galaxies. The other two are less trivial and tend to be larger; they are due to differences in how the background sky is estimated and what model is fit to the surface brightness profile. We show that PyMorph sky estimates are fainter than those of the SDSS DR7 or DR9 pipelines, but are in excellent agreement with the estimates of Blanton et al. (2011). Using the SDSS sky biases luminosities by more than a few tenths of a magnitude for objects with half-light radii > 7 arcseconds. In the SDSS main galaxy sample these are typically luminous galaxies, so they are not necessarily nearby. This bias becomes worse when allowing the model more freedom to fit the surface brightness profile. When PyMorph sky values are used, then two component Sersic-Exponential fits to E+S0s return more light than single component deVaucouleurs fits (up to ~0.2 mag), but less light than single Sersic fits (0.1 mag). Finally, we show that PyMorph fits of Meert et al. (2015) to DR7 data remain valid for DR9 images. Our findings show that, especially at large luminosities, these PyMorph estimates should be preferred to the SDSS pipeline values.

We extend the catalogue of two-dimensional, PSF-corrected de Vacouleurs, Sersic, de Vacouleurs+Exponential, and Sersic+Exponential fits of ~7x10^5 galaxies presented in Meert, Vikram & Bernardi (2015) to include the g- and i-bands. Fits are analysed using the physically motivated flagging system presented in the original text, making adjustments for the differing signal-to-noise when necessary. We compare the fits in each of the g-, r-, and i-bands. Fixed aperture magnitudes and colours are also provided for all galaxies. The catalogues are available in electronic format.

We describe the luminosity function, based on Sersic fits to the light profiles, of CMASS galaxies at z ~ 0.55. Compared to previous estimates, our Sersic-based reductions imply more luminous, massive galaxies, consistent with the effects of Sersic- rather than Petrosian or de Vaucouleur-based photometry on the Sloan Digital Sky Survey (SDSS) main galaxy sample at z ~ 0.1. This implies a significant revision of the high mass end of the correlation between stellar and halo mass. Inferences about the evolution of the luminosity and stellar mass functions depend strongly on the assumed, and uncertain, k+e corrections. In turn, these depend on the assumed age of the population. Applying k+e corrections taken from fitting the models of Maraston et al. (2009) to the colors of both SDSS and CMASS galaxies, the evolution of the luminosity and stellar mass functions appears impressively passive, provided that the fits are required to return old ages. However, when matched in comoving number- or luminosity-density, the SDSS galaxies are less strongly clustered compared to their counterparts in CMASS. This rules out the passive evolution scenario, and, indeed, any minor merger scenarios which preserve the rank ordering in stellar mass of the population. Potential incompletenesses in the CMASS sample would further enhance this mismatch. Our analysis highlights the virtue of combining clustering measurements with number counts.

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 present a catalogue of two-dimensional, point spread function-corrected de Vacouleurs, S\'{e}rsic, de Vacouleurs+Exponential, and S\'{e}rsic+Exponential fits of $\sim7\times10^5$ spectroscopically selected galaxies drawn from the Sloan Digital Sky Survey (SDSS) Data Release 7. Fits are performed for the SDSS r band utilizing the fitting routine GALFIT and analysis pipeline PYMORPH. We compare these fits to prior catalogues. Fits are analysed using a physically motivated flagging system. The flags suggest that more than 90 per cent of two-component fits can be used for analysis. We show that the fits follow the expected behaviour for early and late galaxy types. The catalogues provide a robust set of structural and photometric parameters for future galaxy studies. We show that some biases remain in the measurements, e.g. the presence of bars significantly affect the bulge measurements although the bulge ellipticity may be used to separate barred and non-barred galaxies, and about 15 percent of bulges of two-component fits are also affected by resolution. The catalogues are available in electronic format. We also provide an interface for generating postage stamp images of the 2D model and residual as well as the 1D profile. These images can be generated for a user-uploaded list of galaxies on demand.

The stellar mass-halo mass relation is a key constraint in all semi-analytic, numerical, and semi-empirical models of galaxy formation and evolution. However, its exact shape and redshift dependence remain debated. Several recent works support a relation in the local Universe steeper than previously thought. Based on the comparisons with a variety of data on massive central galaxies, we show that this steepening holds up to z~1, for stellar masses Mstar>2e11 Msun. Specifically, we find significant evidence for a high-mass end slope of \beta>0.35-0.70, instead of the usual \beta~0.20-0.30 reported by a number of previous results. When including the independent constraints from the recent BOSS clustering measurements, the data, independent of any systematic errors in stellar masses, tend to favor a model with a very small scatter (< 0.15 dex) in stellar mass at fixed halo mass, in the redshift range z < 0.8 and for Mstar>3e11 Msun, suggesting a close connection between massive galaxies and host halos even at relatively recent epochs. We discuss the implications of our results with respect to the evolution of the most massive galaxies since z~1.

We quantify the systematics in the size-luminosity relation of galaxies in the SDSS main sample which arise from fitting different 1- and 2-component model profiles to the images. In objects brighter than L*, fitting a single Sersic profile to what is really a two-component SerExp system leads to biases: the half-light radius is increasingly overestimated as n of the fitted single component increases; it is also overestimated at B/T ~ 0.6. However, the net effect on the R-L relation is small, except for the most luminous tail, where it curves upwards towards larger sizes. We also study how this relation depends on morphological type. Our analysis is one of the first to use Bayesian-classifier derived weights, rather than hard cuts, to define morphology. Crudely, there appear to be only two relations: one for early-types (Es, S0s and Sa's) and another for late-types (Sbs and Scds). However, closer inspection shows that within the early-type sample S0s tend to be 15% smaller than Es of the same luminosity, and, among faint late types, Sbs are more than 25% smaller than Scds. Neither the early- nor the late-type relations are pure power-laws: both show significant curvature, which we quantify. However, the R-L relations of the bulges of early-types are almost pure power laws. Our analysis confirms that two mass scales are special for early-type galaxies: M* = 3e10 and 2e11 Msun. These same mass scales are also special for late types: there is almost no correlation between R and M* below the former, and almost no late-types above the latter. We also show that the intrinsic scatter around the relation decreases at large luminosity and/or stellar mass; this should provide additional constraints on models of how the most massive galaxies formed.

The early--type galaxy (ETG) mass--size relation has been largely studied to understand how these galaxies have assembled their mass. One key observational result of the last years is that massive galaxies increased their size by a factor of a few at fixed stellar mass from $z\sim2$. Hierarchical models favor minor mergers as a plausible driver of this size growth. Some of these models, predict a significant environmental dependence in the sense that galaxies residing in more massive halos tend to be larger than galaxies in lower mass halos, at fixed stellar mass and redshift. At present, observational results of this environmental dependence have been contradictory. In this paper we revisit this issue in the local Universe, by investigating how the sizes of massive ETGs depend on large-scale environment using an updated and accurate sample of ETGs in different environments - field, group, clusters - from the Sloan Digital Sky Survey DR7. Our analysis does not show any significant environmental dependence of the sizes of central and satellites ETGs at fixed stellar mass at $z\sim0$. The size-mass relation of early-type galaxies at $z\sim0$ seems to be universal, i.e., independent of the mass of the host halo and of the position of the galaxy in that halo (central or satellite). The result is robust to different galaxy selections based on star formation, morphology or central density. In fact, considering our observational errors and the size of the sample, any size ratio larger than $30-40%$ between massive galaxies ($log(M_*/M_\odot)>11$) living in clusters and in the field can be ruled out at $3\sigma$ level.

In addition to the large systematic differences arising from assumptions about the stellar mass-to-light ratio, the massive end of the stellar mass function is rather sensitive to how one fits the light profiles of the most luminous galaxies. We quantify this by comparing the luminosity and stellar mass functions based on SDSS cmodel magnitudes, and PyMorph single-Sersic and Sersic-Exponential fits to the surface brightness profiles of galaxies in the SDSS. The PyMorph fits return more light, so that the predicted masses are larger than when cmodel magnitudes are used. As a result, the total stellar mass density at z~0.1 is about 1.2x larger than in our previous analysis of the SDSS. The differences are most pronounced at the massive end, where the measured number density of objects having M* > 6 x 10^{11} Msun is ~5x larger. Alternatively, at number densities of 10^{-6} Mpc^{-3}, the limiting stellar mass is 2x larger. The differences with respect to fits by other authors, typically based on Petrosian-like magnitudes, are even more dramatic, although some of these differences are due to sky-subtraction problems, and are sometimes masked by large differences in the assumed $M_*/L$ (even after scaling to the same IMF). Our results impact studies of the growth and assembly of stellar mass in galaxies, and of the relation between stellar and halo mass, so we provide simple analytic fits to these new luminosity and stellar mass functions and quantify how they depend on morphology, as well as the binned counts in electronic format.

We present the results of fitting simulations of an unbiased sample of SDSS galaxies utilizing the fitting routine GALFIT and analysis pipeline PyMorph. These simulations are used to test the two-dimensional decompositions of SDSS galaxies. The simulations show that single S{\'e}rsic models of SDSS data are recovered with $\sigma_{\mathrm{mag}} \approx 0.025$ mag and $\sigma_{\mathrm{radius}} \approx 5\%$. The global values (half-light radius and magnitude) are equally well constrained when a two-component model is used. Sub-components of two-component models present more scatter. SDSS resolution is the primary source of error in the recovery of models. We use a simple statistical correction of the biases in fitted parameters, providing an example using the S{\'e}rsic index. Fitting a two-component S{\'e}rsic + Exponential model to a single S{\'e}rsic galaxy results in a noisier, but unbiased, recovery of the input parameters ($\sigma_{\mathrm{total mag}} \approx 0.075$ mag and $\sigma_{\mathrm{radius}} \approx 10\%$); fitting a single S{\'e}rsic profile to a two-component system results in biases of total magnitude and halflight radius of $\approx 0.05-0.10$ mag and 5\%-10\% in radius. Using an F-test to select the best fit model from among the single- and two-component models is sufficient to remove this bias. We recommend fitting a two-component model to all galaxies when attempting to measure global parameters such as total magnitude and halflight radius.

The statistics of peak counts in reconstructed shear maps contain information beyond the power spectrum, and can improve cosmological constraints from measurements of the power spectrum alone if systematic errors can be controlled. We study the effect of galaxy shape measurement errors on predicted cosmological constraints from the statistics of shear peak counts with the Large Synoptic Survey Telescope (LSST). We use the LSST image simulator in combination with cosmological N-body simulations to model realistic shear maps for different cosmological models. We include both galaxy shape noise and, for the first time, measurement errors on galaxy shapes. We find that the measurement errors considered have relatively little impact on the constraining power of shear peak counts for LSST.

A main science goal for the Large Synoptic Survey Telescope (LSST) is to measure the cosmic shear signal from weak lensing to extreme accuracy. One difficulty, however, is that with the short exposure time ($\simeq$15 seconds) proposed, the spatial variation of the Point Spread Function (PSF) shapes may be dominated by the atmosphere, in addition to optics errors. While optics errors mainly cause the PSF to vary on angular scales similar or larger than a single CCD sensor, the atmosphere generates stochastic structures on a wide range of angular scales. It thus becomes a challenge to infer the multi-scale, complex atmospheric PSF patterns by interpolating the sparsely sampled stars in the field. In this paper we present a new method, PSFent, for interpolating the PSF shape parameters, based on reconstructing underlying shape parameter maps with a multi-scale maximum entropy algorithm. We demonstrate, using images from the LSST Photon Simulator, the performance of our approach relative to a 5th-order polynomial fit (representing the current standard) and a simple boxcar smoothing technique. Quantitatively, PSFent predicts more accurate PSF models in all scenarios and the residual PSF errors are spatially less correlated. This improvement in PSF interpolation leads to a factor of 3.5 lower systematic errors in the shear power spectrum on scales smaller than $\sim13'$, compared to polynomial fitting. We estimate that with PSFent and for stellar densities greater than $\simeq1/{\rm arcmin}^{2}$, the spurious shear correlation from PSF interpolation, after combining a complete 10-year dataset from LSST, is lower than the corresponding statistical uncertainties on the cosmic shear power spectrum, even under a conservative scenario.

The complete 10-year survey from the Large Synoptic Survey Telescope (LSST) will image $\sim$ 20,000 square degrees of sky in six filter bands every few nights, bringing the final survey depth to $r\sim27.5$, with over 4 billion well measured galaxies. To take full advantage of this unprecedented statistical power, the systematic errors associated with weak lensing measurements need to be controlled to a level similar to the statistical errors. This work is the first attempt to quantitatively estimate the absolute level and statistical properties of the systematic errors on weak lensing shear measurements due to the most important physical effects in the LSST system via high fidelity ray-tracing simulations. We identify and isolate the different sources of algorithm-independent, \textit{additive} systematic errors on shear measurements for LSST and predict their impact on the final cosmic shear measurements using conventional weak lensing analysis techniques. We find that the main source of the errors comes from an inability to adequately characterise the atmospheric point spread function (PSF) due to its high frequency spatial variation on angular scales smaller than $\sim10'$ in the single short exposures, which propagates into a spurious shear correlation function at the $10^{-4}$--$10^{-3}$ level on these scales. With the large multi-epoch dataset that will be acquired by LSST, the stochastic errors average out, bringing the final spurious shear correlation function to a level very close to the statistical errors. Our results imply that the cosmological constraints from LSST will not be severely limited by these algorithm-independent, additive systematic effects.

We present basic predictions of an updated version of the Munich semi-analytic hierarchical galaxy formation model that grows bulges via mergers and disk instabilities. Overall, we find that while spheroids below Ms ~ 10^11 Msun grow their sizes via a mixture of disk instability and mergers, galaxies above it mainly evolve via mergers. Including gas dissipation in major mergers, efficiently shrinks galaxies, especially those with final mass Ms < 10^11 Msun that are the most gas-rich, improving the match with different observables. We find that the predicted scatter in sizes at fixed stellar mass is still larger than the observed one by up to <40%. Spheroids are, on average, more compact at higher redshifts at fixed stellar mass, and at fixed redshift and stellar mass larger galaxies tend to be more starforming. More specifically, while for bulge-dominated galaxies the model envisages a nearly mass-independent decrease in sizes, the predicted size evolution for intermediate-mass galaxies is more complex. The z=2 progenitors of massive galaxies with mass around Ms and B/T>0.7 at z=0, are found to be mostly disc-dominated galaxies with a median B/T ~ 0.3, with only ~20% remaining bulge-dominated. The model also predicts that central spheroids living in more massive haloes tend to have larger sizes at fixed stellar mass. Including host halo mass dependence in computing velocity dispersions, allows the model to properly reproduce the correlations with stellar mass. We also discuss the fundamental plane, the correlations with galaxy age, the structural properties of pseudobulges, and the correlations with central black holes.

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.