This is the third paper of a series where we study the stellar population gradients (SP; ages, metallicities, $\alpha$-element abundance ratios and stellar initial mass functions) of early type galaxies (ETGs) at $z\le 0.08$ from the MaNGA-DR15 survey. In this work we focus on the S0 population and quantify how the SP varies across the population as well as with galactocentric distance. We do this by measuring Lick indices and comparing them to stellar population synthesis models. This requires spectra with high signal-to-noise which we achieve by stacking in bins of luminosity (L$_r$) and central velocity dispersion ($\sigma_0$). We find that: 1) There is a bimodality in the S0 population: S0s more massive than $3\times 10^{10}M_\odot$ show stronger velocity dispersion and age gradients (age and $\sigma_r$ decrease outwards) but little or no metallicity gradient, while the less massive ones present relatively flat age and velocity dispersion profiles, but a significant metallicity gradient (i.e. [M/H] decreases outwards). Above $2\times10^{11}M_\odot$ the number of S0s drops sharply. These two mass scales are also where global scaling relations of ETGs change slope. 2) S0s have steeper velocity dispersion profiles than fast rotating elliptical galaxies (E-FRs) of the same luminosity and velocity dispersion. The kinematic profiles and stellar population gradients of E-FRs are both more similar to those of slow rotating ellipticals (E-SRs) than to S0s, suggesting that E-FRs are not simply S0s viewed face-on. 3) At fixed $\sigma_0$, more luminous S0s and E-FRs are younger, more metal rich and less $\alpha$-enhanced. Evidently for these galaxies, the usual statement that 'massive galaxies are older' is not true if $\sigma_0$ is held fixed.

This paper documents the sixteenth data release (DR16) from the Sloan Digital Sky Surveys; the fourth and penultimate from the fourth phase (SDSS-IV). This is the first release of data from the southern hemisphere survey of the Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2); new data from APOGEE-2 North are also included. DR16 is also notable as the final data release for the main cosmological program of the Extended Baryon Oscillation Spectroscopic Survey (eBOSS), and all raw and reduced spectra from that project are released here. DR16 also includes all the data from the Time Domain Spectroscopic Survey (TDSS) and new data from the SPectroscopic IDentification of ERosita Survey (SPIDERS) programs, both of which were co-observed on eBOSS plates. DR16 has no new data from the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey (or the MaNGA Stellar Library "MaStar"). We also preview future SDSS-V operations (due to start in 2020), and summarize plans for the final SDSS-IV data release (DR17).

Early-type galaxies -- slow and fast rotating ellipticals (E-SRs and E-FRs) and S0s/lenticulars -- define a Fundamental Plane (FP) in the space of half-light radius $R_e$, enclosed surface brightness $I_e$ and velocity dispersion $\sigma_e$. Since $I_e$ and $\sigma_e$ are distance-independent measurements, the thickness of the FP is often expressed in terms of the accuracy with which $I_e$ and $\sigma_e$ can be used to estimate sizes $R_e$. We show that: 1) The thickness of the FP depends strongly on morphology. If the sample only includes E-SRs, then the observed scatter in $R_e$ is $\sim 16\%$, of which only $\sim 9\%$ is intrinsic. Removing galaxies with $M_*<10^{11}M_\odot$ further reduces the observed scatter to $\sim 13\%$ ($\sim 4\%$ intrinsic). The observed scatter increases to the $\sim 25\%$ usually quoted in the literature if E-FRs and S0s are added. If the FP is defined using the eigenvectors of the covariance matrix of the observables, then the E-SRs again define an exceptionally thin FP, with intrinsic scatter of only $5\%$ orthogonal to the plane. 2) The structure within the FP is most easily understood as arising from the fact that $I_e$ and $\sigma_e$ are nearly independent, whereas the $R_e-I_e$ and $R_e-\sigma_e$ correlations are nearly equal and opposite. 3) If the coefficients of the FP differ from those associated with the virial theorem the plane is said to be `tilted'. If we multiply $I_e$ by the global stellar mass-to-light ratio $M_*/L$ and we account for non-homology across the population by using S\'ersic photometry, then the resulting stellar mass FP is less tilted. Accounting self-consistently for $M_*/L$ gradients will change the tilt. The tilt we currently see suggests that the efficiency of turning baryons into stars increases and/or the dark matter fraction decreases as stellar surface brightness increases.

With the advent of future big-data surveys, automated tools for unsupervised discovery are becoming ever more necessary. In this work, we explore the ability of deep generative networks for detecting outliers in astronomical imaging datasets. The main advantage of such generative models is that they are able to learn complex representations directly from the pixel space. Therefore, these methods enable us to look for subtle morphological deviations which are typically missed by more traditional moment-based approaches. We use a generative model to learn a representation of expected data defined by the training set and then look for deviations from the learned representation by looking for the best reconstruction of a given object. In this first proof-of-concept work, we apply our method to two different test cases. We first show that from a set of simulated galaxies, we are able to detect $\sim90\%$ of merging galaxies if we train our network only with a sample of isolated ones. We then explore how the presented approach can be used to compare observations and hydrodynamic simulations by identifying observed galaxies not well represented in the models.

The masses of supermassive black holes at the centres of local galaxies appear to be tightly correlated with the mass and velocity dispersions of their galactic hosts. However, the local Mbh-Mstar relation inferred from dynamically measured inactive black holes is up to an order-of-magnitude higher than some estimates from active black holes, and recent work suggests that this discrepancy arises from selection bias on the sample of dynamical black hole mass measurements. In this work we combine X-ray measurements of the mean black hole accretion luminosity as a function of stellar mass and redshift with empirical models of galaxy stellar mass growth, integrating over time to predict the evolving Mbh-Mstar relation. The implied relation is nearly independent of redshift, indicating that stellar and black hole masses grow, on average, at similar rates. Matching the de-biased local Mbh-Mstar relation requires a mean radiative efficiency ~0.15, in line with theoretical expectations for accretion onto spinning black holes. However, matching the "raw" observed relation for inactive black holes requires a mean radiative efficiency around 0.02, far below theoretical expectations. This result provides independent evidence for selection bias in dynamically estimated black hole masses, a conclusion that is robust to uncertainties in bolometric corrections, obscured active black hole fractions, and kinetic accretion efficiency. For our fiducial assumptions, they favour moderate-to-rapid spins of typical supermassive black holes, to achieve a mean radiative efficiency ~0.12-0.20. Our approach has similarities to the classic Soltan analysis, but by using galaxy-based data instead of integrated quantities we are able to focus on regimes where observational uncertainties are minimized.

The mass and structural assembly of galaxies is a matter of intense debate. Current theoretical models predict the existence of a linear relationship between galaxy size ($R_e$) and the host dark matter halo virial radius ($R_h$).\\ By making use of semi-empirical models compared to the size distributions of central galaxies from the Sloan Digital Sky Survey, we provide robust constraints on the normalization and scatter of the $R_e-R_h$ relation. We explore the parameter space of models in which the $R_e-R_h$ relation is mediated by either the spin parameter or the concentration of the host halo, or a simple constant the nature of which is in principle unknown. We find that the data require extremely tight relations for both early-type and late-type galaxies (ETGs,LTGs), especially for more massive galaxies. These constraints challenge models based solely on angular momentum conservation, which predict significantly wider distributions of galaxy sizes and no trend with stellar mass, if taken at face value. We discuss physically-motivated alterations to the original models that bring the predictions into better agreement with the data. We argue that the measured tight size distributions of SDSS disk galaxies can be reproduced by semi-empirical models in which the $R_e-R_h$ connection is mediated by the \emph{stellar} specific angular momenta $j_{star}.$ We find that current cosmological models of galaxy formation broadly agree with our constraints for LTGs, and justify the strong link between $R_e$ and $j_{star}$ that we propose, however the tightness of the $R_e-R_h$ relation found in such ab-initio theoretical models for ETGs is in tension with our semi-empirical findings.

A supermassive black hole has been found at the centre of nearly every galaxy observed with sufficient sensitivity. The masses of these black holes are observed to increase with either the total mass or the mean (random) velocity of the stars in their host galaxies. The origin of these correlations remains elusive. Observational systematics and biases severely limit our knowledge of the local demography of supermassive black holes thus preventing accurate model comparisons and progress in this field. Here we show that the large-scale spatial distribution of local active galactic nuclei (AGN), believed to be accreting supermassive black holes, can constrain the shape and normalization of the black hole-stellar mass relation thus bypassing resolution-related observational biases. In turn, our results can set more stringent constraints on the so-called "radiative efficiency", a fundamental parameter describing the inner physics of supermassive black holes that is closely linked to their spin, geometry, and ability to release energy. The mean value of the radiative efficiency can be estimated by comparing the average total luminous output of AGN with the relic mass density locked up in quiescent supermassive black holes at galaxy centres today. For currently accepted values of the AGN obscured fractions and bolometric corrections, our newest estimates of the local supermassive black hole mass density favour mean radiative efficiencies of ~10-20\%, suggesting that the vast majority of supermassive black holes are spinning moderately to rapidly. With large-scale AGN surveys coming online, our novel methodology will enable even tighter constraints on the fundamental parameters that regulate the growth of supermassive black holes.

We estimate ages, metallicities, $\alpha$-element abundance ratios and stellar initial mass functions of elliptical (E) and S0 galaxies from the MaNGA-DR15 survey. We stack spectra and use a variety of single stellar population synthesis models to interpret the absorption line strengths in these spectra. We quantify how these properties vary across the population, as well as with galactocentric distance. This paper is the first of a series and is based on a sample of pure elliptical galaxies at z $\le$ 0.08. We show that the properties of the inner regions of Es with the largest luminosity (L$_r$) and central velocity dispersion ($\sigma_0$) are consistent with those associated with the commonly used Salpeter IMF, whereas a Kroupa-like IMF is a better description at $\sim$ 0.8R/Re (assuming [Ti/Fe] variations are limited). For these galaxies the stellar mass-to-light ratio decreases at most by a factor of 2 from the central regions to Re. In contrast, for lower L$_r$ and $\sigma_0$ galaxies, the IMF is shallower and M$_{*}$/L$_r$ in the central regions is similar to the outskirts. Although a factor of 2 is smaller than previous reports based on a handful of galaxies, it is still large enough to matter for dynamical mass estimates. Accounting self-consistently for these gradients when estimating both M$_{*}$ and M$_{dyn}$ brings the two into good agreement: gradients reduce M$_{dyn}$ by $\sim$ 0.2 dex while only slightly increasing the M$_{*}$ inferred using a Kroupa IMF. This is a different resolution of the M$_{*}$-M$_{dyn}$ discrepancy than has been followed in the recent literature where M$_{*}$ of massive galaxies is increased by adopting a Salpeter IMF while leaving Mdyn unchanged. A companion paper discusses how stellar population differences are even more pronounced if one separates slow from fast rotators.

We present estimates of stellar population (SP) gradients from stacked spectra of slow (SR) and fast (FR) rotator elliptical galaxies from the MaNGA-DR15 survey. We find that: 1) FRs are $\sim 5$ Gyrs younger, more metal rich, less $\alpha$-enhanced and smaller than SRs of the same luminosity $L_r$ and central velocity dispersion $\sigma_0$. This explains why when one combines SRs and FRs, objects which are small for their $L_r$ and $\sigma_0$ tend to be younger. Their SP gradients are also different. 2) Ignoring the FR/SR dichotomy leads one to conclude that compact galaxies are older than their larger counterparts of the same mass, even though almost the opposite is true for FRs and SRs individually. 3) SRs with $\sigma_0\le 250$ km s$^{-1}$ are remarkably homogeneous within $\sim R_e$: they are old, $\alpha$-enhanced and only slightly super-solar in metallicity. These SRs show no gradients in age and $M_*/L_r$, negative gradients in metallicity, and slightly positive gradients in [$\alpha$/Fe] (the latter are model dependent). SRs with $\sigma_0\ge 250$ km $s^{-1}$ are slightly younger and more metal rich, contradicting previous work suggesting that age increases with $\sigma_0$. They also show larger $M_*/L_r$ gradients. 4) Self-consistently accounting for $M_*/L$ gradients yields $M_{\rm dyn}\approx M_*$ because gradients reduce $M_{\rm dyn}$ by $\sim 0.2$ dex while only slightly increasing the $M_*$ inferred using a Kroupa (not Salpeter) IMF. 5) The FR population all but disappears above $M_*\ge 3\times 10^{11}M_\odot$; this is the same scale at which the size-mass correlation and other scaling relations change. Our results support the finding that this is an important mass scale which correlates with the environment and above which mergers matter.

The central regions of nearby elliptical galaxies are dominated by baryons (stars) and provide interesting laboratories for studying the radial acceleration relation (RAR). We carry out exploratory analyses and discuss the possibility of constraining the RAR in the super-critical acceleration range $(10^{-9.5},\hspace{1ex}10^{-8})$~${\rm m}~{\rm s}^{-2}$ using a sample of nearly round pure-bulge (spheroidal, dispersion-dominated) galaxies including 24 ATLAS$^{\rm 3D}$ galaxies and 4201 SDSS galaxies covering a wide range of masses, sizes and luminosity density profiles. We consider a range of current possibilities for the stellar mass-to-light ratio ($M_\star/L$), its gradient and dark or phantom matter (DM/PM) halo profiles. We obtain the probability density functions (PDFs) of the parameters of the considered models via Bayesian inference based on spherical Jeans Monte Carlo modeling of the observed velocity dispersions. We then constrain the DM/PM-to-baryon acceleration ratio $a_{\rm X}/a_{\rm B}$ from the PDFs. Unless we ignore observed radial gradients in $M_\star/L$, or assume unreasonably strong gradients, marginalization over nuisance factors suggests $a_{\rm X}/a_{\rm B} = 10^{p} (a_{\rm B}/a_{+1})^q$ with $p = -1.00 \pm 0.03$ (stat) $^{+0.11}_{-0.06}$ (sys) and $q=-1.02 \pm 0.09$ (stat) $^{+0.16}_{-0.00}$ (sys) around a super-critical acceleration $a_{+1}\equiv 1.2\times 10^{-9}~{\rm m}~{\rm s}^{-2}$. In the context of the $\Lambda$CDM paradigm, this RAR suggests that the NFW DM halo profile is a reasonable description of galactic halos even after the processes of galaxy formation and evolution. In the context of the MOND paradigm, this RAR favors the Simple interpolating function but is inconsistent with the vast majority of other theoretical proposals and fitting functions motivated mainly from sub-critical acceleration data.

We analyze the optical morphologies of galaxies in the IllustrisTNG simulation at $z\sim0$ with a Convolutional Neural Network trained on visual morphologies in the Sloan Digital Sky Survey. We generate mock SDSS images of a mass complete sample of $\sim12,000$ galaxies in the simulation using the radiative transfer code SKIRT and include PSF and noise to match the SDSS r-band properties. The images are then processed through the exact same neural network used to estimate SDSS morphologies to classify simulated galaxies in four morphological classes (E, S0/a, Sab, Scd). The CNN model finds that $\sim95\%$ of the simulated galaxies fall in one the four main classes with high confidence. The mass-size relations of the simulated galaxies divided by morphological type also reproduce well the slope and the normalization of observed relations which confirms the realism of optical morphologies in the TNG suite. However, the Stellar Mass Functions decomposed into different morphologies still show significant discrepancies with observations both at the low and high mass end. We find that the high mass end of the SMF is dominated in TNG by massive disk galaxies while early-type galaxies dominate in the observations according to the CNN classifications. The present work highlights the importance of detailed comparisons between observations and simulations in comparable conditions.

We investigate the anisotropy of the stellar velocity dispersions within the effective radius, $R_{\rm e}$, in 24 ATLAS$^{\rm 3D}$ pure-bulge galaxies, 16 of which are kinematic slow rotators (SRs). We allow the spherical anisotropy parameter $\beta$ to be radially varying and allow a radial gradient in the stellar mass-to-light ratio ($M_\star/L$) through the parameter $K$ introduced earlier. The median anisotropy for SRs depends on $K$ as follows: $\langle\beta_{\rm m}\rangle = a + b K$ with $a=0.19\pm 0.05$, $b=-0.13\pm 0.07$ ($\Lambda$CDM) or $a=0.21\pm 0.05$, $b=-0.26\pm 0.08$ (MOND), where $\beta_{\rm m}$ refers to the radially averaged quantity. Under the $\Lambda$CDM paradigm this scaling is tied to a scaling of $\langle f_{\rm DM}\rangle = (0.16\pm 0.03) +(0.31\pm 0.06) K$, where $f_{\rm DM}$ refers to the DM fraction within a sphere of $r=R_{\rm e}$. For $K=0$ (constant $M_\star/L$), we obtain radially biased results with $\langle\beta_{\rm m}\rangle \approx 0.2$ consistent with previous results. However, marginalizing over $0 < K < 1.5$ yields $\langle\beta_{\rm m}\rangle = 0.06 ^{+0.11}_{-0.14}$ with $\langle f_{\rm DM}\rangle = 0.35 \pm 0.08$: isotropy is preferred. This isotropy hides the fact that $\beta_{\rm m}$ is correlated with kinematic features such as counter rotating cores (CRCs), kinematically distinct cores (KDCs), and low-level velocities (LVs): SRs with LVs are likely to be radially biased while SRs with CRCs are likely to be tangentially biased, and SRs with KDCs are intermediate. Existing cosmological simulations allow us to understand these results qualitatively in terms of their dynamical structures and formation histories although there exist quantitative tensions. More realistic cosmological simulations, particularly allowing for $M_\star/L$ gradients, may be required to better understand SRs.

We study the rest-frame optical mass-size relation of bulges and discs from z~2 to z~0 for a complete sample of massive galaxies in the CANDELS fields using 2 component S\'ersic models (Dimauro et a. 2018). Discs and star forming galaxies follow similar mass-size relations. The mass-size relation of bulges is less steep than the one of quiescent galaxies (best fit slope of ~ 0.7 for quiescent galaxies against ~ 0.4 for bulges). We find little dependence of the structural properties of massive bulges and discs with the global morphology of galaxies (disc vs. bulge dominated) and the star formation activity (star-forming vs. quiescent). This result suggests similar bulge formation mechanisms for most massive galaxies and also that the formation of the bulge component does not significantly affect the disc structure. Our findings pose a challenge to models envisioning multiple channels for massive bulge growth, such as disc instabilities and mergers.

Local samples of quiescent galaxies with dynamically measured black hole masses (Mbh) may suffer from an angular resolution-related selection effect, which could bias the observed scaling relations between Mbh and host galaxy properties away from the intrinsic relations. In particular, previous work has shown that the observed Mbh-Mstar (stellar mass) relation is more strongly biased than the Mbh-sigma (velocity dispersion) relation. Local samples of active galactic nuclei (AGN) do not suffer from this selection effect, as in these samples Mbh is estimated from megamasers and/or reverberation mapping-based techniques. With the exception of megamasers, Mbh-estimates in these AGN samples are proportional to a virial coefficient fvir. Direct modelling of the broad line region suggests that fvir~3.5. However, this results in a Mbh-Mstar relation for AGN which lies below and is steeper than the one observed for quiescent black hole samples. A similar though milder trend is seen for the Mbh-sigma relation. Matching the high-mass end of the Mbh-Mstar and Mbh-sigma relations observed in quiescent samples requires fvir~15 and fvir~7, respectively. On the other hand, fvir~3.5 yields Mbh-sigma and Mbh-Mstar relations for AGN which are remarkably consistent with the expected `intrinsic' correlations for quiescent samples (i.e., once account has been made of the angular resolution-related selection effect), providing additional evidence that the sample of local quiescent black holes is biased. We also show that, as is the case for quiescent black holes, the Mbh-Mstar scaling relation of AGN is driven by velocity dispersion, thus providing additional key constraints to black hole-galaxy co-evolution models.

Deep Learning (DL) algorithms for morphological classification of galaxies have proven very successful, mimicking (or even improving) visual classifications. However, these algorithms rely on large training samples of labelled galaxies (typically thousands of them). A key question for using DL classifications in future Big Data surveys is how much of the knowledge acquired from an existing survey can be exported to a new dataset, i.e. if the features learned by the machines are meaningful for different data. We test the performance of DL models, trained with Sloan Digital Sky Survey (SDSS) data, on Dark Energy survey (DES) using images for a sample of $\sim$5000 galaxies with a similar redshift distribution to SDSS. Applying the models directly to DES data provides a reasonable global accuracy ($\sim$ 90%), but small completeness and purity values. A fast domain adaptation step, consisting in a further training with a small DES sample of galaxies ($\sim$500-300), is enough for obtaining an accuracy > 95% and a significant improvement in the completeness and purity values. This demonstrates that, once trained with a particular dataset, machines can quickly adapt to new instrument characteristics (e.g., PSF, seeing, depth), reducing by almost one order of magnitude the necessary training sample for morphological classification. Redshift evolution effects or significant depth differences are not taken into account in this study.

We describe the SDSS-IV MaNGA PyMorph Photometric (MPP-VAC) and MaNGA Deep Learning Morphology (MDLM-VAC) Value Added Catalogs. The MPP-VAC provides photometric parameters from S\'ersic and S\'ersic+Exponential fits to the 2D surface brightness profiles of the MaNGA DR15 galaxy sample. Compared to previous PyMorph analyses of SDSS imaging, our analysis of the MaNGA DR15 incorporates three improvements: the most recent SDSS images; modified criteria for determining bulge-to-disk decompositions; and the fits in MPP-VAC have been eye-balled, and re-fit if necessary, for additional reliability. A companion catalog, the MDLM-VAC, provides Deep Learning-based morphological classifications for the same galaxies. The MDLM-VAC includes a number of morphological properties (e.g., a TType, and a finer separation between elliptical and S0 galaxies). Combining the MPP- and MDLM-VACs allows to show that the MDLM morphological classifications are more reliable than previous work. It also shows that single-S\'ersic fits to late- and early-type galaxies are likely to return S\'ersic indices of $n \le 2$ and $\ge 4$, respectively, and this correlation between $n$ and morphology extends to the bulge component as well. While the former is well-known, the latter contradicts some recent work suggesting little correlation between $n$-bulge and morphology. Combining both VACs with MaNGA's spatially resolved spectroscopy allows us to study how the stellar angular momentum depends on morphological type. We find correlations between stellar kinematics, photometric properties, and morphological type even though the spectroscopic data played no role in the construction of the MPP- and MDLM-VACs.

Twenty years have passed since first light for the Sloan Digital Sky Survey (SDSS). Here, we release data taken by the fourth phase of SDSS (SDSS-IV) across its first three years of operation (July 2014-July 2017). This is the third data release for SDSS-IV, and the fifteenth from SDSS (Data Release Fifteen; DR15). New data come from MaNGA - we release 4824 datacubes, as well as the first stellar spectra in the MaNGA Stellar Library (MaStar), the first set of survey-supported analysis products (e.g. stellar and gas kinematics, emission line, and other maps) from the MaNGA Data Analysis Pipeline (DAP), and a new data visualisation and access tool we call "Marvin". The next data release, DR16, will include new data from both APOGEE-2 and eBOSS; those surveys release no new data here, but we document updates and corrections to their data processing pipelines. The release is cumulative; it also includes the most recent reductions and calibrations of all data taken by SDSS since first light. In this paper we describe the location and format of the data and tools and cite technical references describing how it was obtained and processed. The SDSS website (www.sdss.org) has also been updated, providing links to data downloads, tutorials and examples of data use. While SDSS-IV will continue to collect astronomical data until 2020, and will be followed by SDSS-V (2020-2025), we end this paper by describing plans to ensure the sustainability of the SDSS data archive for many years beyond the collection of data.

We present STEEL a STatistical sEmi-Empirical modeL designed to probe the distribution of satellite galaxies in groups and clusters. Our fast statistical methodology relies on tracing the abundances of central and satellite haloes via their mass functions at all cosmic epochs with virtually no limitation on cosmic volume and mass resolution. From mean halo accretion histories and subhalo mass functions the satellite mass function is progressively built in time via abundance matching techniques constrained by number densities of centrals in the local Universe. By enforcing dynamical merging timescales as predicted by high-resolution N-body simulations, we obtain satellite distributions as a function of stellar mass and halo mass consistent with current data. We show that stellar stripping, star formation, and quenching play all a secondary role in setting the number densities of massive satellites above $M_*\gtrsim 3\times 10^{10}\, M_{\odot}$. We further show that observed star formation rates used in our empirical model over predict low-mass satellites below $M_*\lesssim 3\times 10^{10}\, M_{\odot}$, whereas, star formation rates derived from a continuity equation approach yield the correct abundances similar to previous results for centrals.

Gas stripping of spiral galaxies or mergers are thought to be the formation mechanisms of lenticular galaxies. In order to determine the conditions in which each scenario dominates, we derive stellar populations of both the bulge and disk regions of 279 lenticular galaxies in the MaNGA survey. We find a clear bimodality in stellar age and metallicity within the population of S0s and this is strongly correlated with stellar mass. Old and metal-rich bulges and disks belong to massive galaxies, and young and metal-poor bulges and disks are hosted by low-mass galaxies. From this we conclude that the bulges and disks are co-evolving. When the bulge and disk stellar ages are compared, we find that the bulge is almost always older than the disk for massive galaxies ($\textrm{M}_{\star} > 10^{10}~\textrm{M}_{\odot}$). The opposite is true for lower mass galaxies. We conclude that we see two separate populations of lenticular galaxies. The old, massive, and metal-rich population possess bulges that are predominantly older than their disks, which we speculate may have been caused by morphological or inside-out quenching. In contrast, the less massive and more metal-poor population have bulges with more recent star formation than their disks. We postulate they may be undergoing bulge rejuvenation (or disk fading), or compaction. Environment doesn't play a distinct role in the properties of either population. Our findings give weight to the notion that while the faded spiral scenario likely formed low-mass S0s, other processes, such as mergers, may be responsible for high-mass S0s.

We carry out spherical Jeans modeling of nearly round pure-bulge galaxies selected from the ATLAS$^{\rm 3D}$ sample. Our modeling allows for gradients in the stellar mass-to-light ratio ($M_\star/L$) through analytic prescriptions parameterized with a `gradient strength' $K$ introduced to accommodate any viable gradient. We use a generalized Osipkov-Merritt model for the velocity dispersion (VD) anisotropy. We produce Monte Carlo sets of models based on the stellar VD profiles under both the $\Lambda$CDM and MOND paradigms. Here, we describe the galaxy data, the empirical inputs, and the modeling procedures of obtaining the Monte Carlo sets. We then present the projected dynamical stellar mass, $M_{\rm \star e}$, within the effective radius $R_{\rm e}$, and the fundamental mass plane (FMP) as a function of $K$. We find the scaling of the $K$-dependent mass with respect to the ATLAS$^{\rm 3D}$ reported mass as: $\log_{10} \left[M_{\star{\rm e}}(K)/M_{\star{\rm e}}^{\rm A3D} \right] = a' + b' K$ with $a'=-0.019\pm 0.012$ and $b'=-0.18\pm 0.02$ ($\Lambda$CDM), or $a'=-0.023\pm 0.014$ and $b'=-0.23\pm 0.03$ (MOND), for $0\le K < 1.5$. The FMP has coefficients consistent with the virial expectation and only the zero-point scales with $K$. The median value of $K$ for the ATLAS$^{\rm 3D}$ galaxies is $\langle K\rangle =0.53^{+0.05}_{-0.04}$. We perform a similar analysis of the much larger SDSS DR7 spectroscopic sample. In this case, only the VD within a single aperture is available, so we impose the additional requirement that the VD slope be similar to that in the ATLAS$^{\rm 3D}$ galaxies. Our analysis of the SDSS galaxies suggests a positive correlation of $K$ with stellar mass.

The fourth generation of the Sloan Digital Sky Survey (SDSS-IV) has been in operation since July 2014. This paper describes the second data release from this phase, and the fourteenth from SDSS overall (making this, Data Release Fourteen or DR14). This release makes public data taken by SDSS-IV in its first two years of operation (July 2014-2016). Like all previous SDSS releases, DR14 is cumulative, including the most recent reductions and calibrations of all data taken by SDSS since the first phase began operations in 2000. New in DR14 is the first public release of data from the extended Baryon Oscillation Spectroscopic Survey (eBOSS); the first data from the second phase of the Apache Point Observatory (APO) Galactic Evolution Experiment (APOGEE-2), including stellar parameter estimates from an innovative data driven machine learning algorithm known as "The Cannon"; and almost twice as many data cubes from the Mapping Nearby Galaxies at APO (MaNGA) survey as were in the previous release (N = 2812 in total). This paper describes the location and format of the publicly available data from SDSS-IV surveys. We provide references to the important technical papers describing how these data have been taken (both targeting and observation details) and processed for scientific use. The SDSS website (www.sdss.org) has been updated for this release, and provides links to data downloads, as well as tutorials and examples of data use. SDSS-IV is planning to continue to collect astronomical data until 2020, and will be followed by SDSS-V.

We use machine learning to identify in color images of high-redshift galaxies an astrophysical phenomenon predicted by cosmological simulations. This phenomenon, called the blue nugget (BN) phase, is the compact star-forming phase in the central regions of many growing galaxies that follows an earlier phase of gas compaction and is followed by a central quenching phase. We train a Convolutional Neural Network (CNN) with mock "observed" images of simulated galaxies at three phases of evolution: pre-BN, BN and post-BN, and demonstrate that the CNN successfully retrieves the three phases in other simulated galaxies. We show that BNs are identified by the CNN within a time window of $\sim0.15$ Hubble times. When the trained CNN is applied to observed galaxies from the CANDELS survey at $z=1-3$, it successfully identifies galaxies at the three phases. We find that the observed BNs are preferentially found in galaxies at a characteristic stellar mass range, $10^{9.2-10.3} M_\odot$ at all redshifts. This is consistent with the characteristic galaxy mass for BNs as detected in the simulations, and is meaningful because it is revealed in the observations when the direct information concerning the total galaxy luminosity has been eliminated from the training set. This technique can be applied to the classification of other astrophysical phenomena for improved comparison of theory and observations in the era of large imaging surveys and cosmological simulations.

Understanding how bulges grow in galaxies is critical step towards unveiling the link between galaxy morphology and star-formation. To do so, it is necessary to decompose large sample of galaxies at different epochs into their main components (bulges and disks). This is particularly challenging, especially at high redshifts, where galaxies are poorly resolved. This work presents a catalog of bulge-disk decompositions of the surface brightness profiles of ~17.600 H-band selected galaxies in the CANDELS fields (F160W<23, 0<z<2) in 4 to 7 filters covering a spectral range of 430-1600nm. This is the largest available catalog of this kind up to z = 2. By using a novel approach based on deep-learning to select the best model to fit, we manage to control systematics arising from wrong model selection and obtain less contaminated samples than previous works. We show that the derived structural properties are within $\sim10-20\%$ of random uncertainties. We then fit stellar population models to the decomposed SEDs (Spectral Energy Distribution) of bulges and disks and derive stellar masses (and stellar mass bulge-to-total ratios) as well as rest-frame colors (U,V,J) for bulges and disks separately. All data products are publicly released with this paper and through the web page https://lerma.obspm.fr/huertas/form_CANDELS and will be used for scientific analysis in forthcoming works.

Within a galaxy the stellar mass-to-light ratio $\Upsilon_*$ is not constant. Spatially resolved kinematics of nearby early-type galaxies suggest that allowing for a variable initial mass function (IMF) returns significantly larger $\Upsilon_*$ gradients than if the IMF is held fixed. If $\Upsilon_*$ is greater in the central regions, then ignoring the IMF-driven gradient can overestimate $M_*^{\rm dyn}$ by as much as a factor of two for the most massive galaxies, though stellar population estimates $M_*^{\rm SP}$ are also affected. Large $\Upsilon_*$-gradients have four main consequences: First, $M_*^{\rm dyn}$ cannot be estimated independently of stellar population synthesis models. Second, if there is a lower limit to $\Upsilon_*$ and gradients are unknown, then requiring $M_*^{\rm dyn}=M_*^{\rm SP}$ constrains them. Third, if gradients are stronger in more massive galaxies, then $M_*^{\rm dyn}$ and $M_*^{\rm SP}$ can be brought into agreement, not by shifting $M_*^{\rm SP}$ upwards by invoking constant bottom-heavy IMFs, as advocated by a number of recent studies, but by revising $M_*^{\rm dyn}$ estimates in the literature downwards. Fourth, accounting for $\Upsilon_*$ gradients changes the high-mass slope of the stellar mass function $\phi(M_*^{\rm dyn})$, and reduces the associated stellar mass density. These conclusions potentially impact estimates of the need for feedback and adiabatic contraction, so our results highlight the importance of measuring $\Upsilon_*$ gradients in larger samples.

We present a morphological catalogue for $\sim$ 670,000 galaxies in the Sloan Digital Sky Survey in two flavours: T-Type, related to the Hubble sequence, and Galaxy Zoo 2 (GZ2 hereafter) classification scheme. By combining accurate existing visual classification catalogues with machine learning, we provide the largest and most accurate morphological catalogue up to date. The classifications are obtained with Deep Learning algorithms using Convolutional Neural Networks (CNNs). We use two visual classification catalogues, GZ2 and Nair & Abraham (2010), for training CNNs with colour images in order to obtain T-Types and a series of GZ2 type questions (disk/features, edge-on galaxies, bar signature, bulge prominence, roundness and mergers). We also provide an additional probability enabling a separation between pure elliptical (E) from S0, where the T-Type model is not so efficient. For the T-Type, our results show smaller offset and scatter than previous models trained with support vector machines. For the GZ2 type questions, our models have large accuracy (> 97\%), precision and recall values (> 90\%) when applied to a test sample with the same characteristics as the one used for training. The catalogue is publicly released with the paper.

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.

We carry out a systematic investigation of the total mass density profile of massive (Mstar~3e11 Msun) early-type galaxies and its dependence on redshift, specifically in the range 0<z<1. We start from a large sample of SDSS early-type galaxies with stellar masses and effective radii measured assuming two different profiles, de Vaucouleurs and S\'{e}rsic. We assign dark matter haloes to galaxies via abundance matching relations with standard LCDM profiles and concentrations. We then compute the total, mass-weighted density slope at the effective radius gamma', and study its redshift dependence at fixed stellar mass. We find that a necessary condition to induce an increasingly flatter gamma' at higher redshifts, as suggested by current strong lensing data, is to allow the intrinsic stellar profile of massive galaxies to be S\'{e}rsic and the input S\'{e}rsic index n to vary with redshift approximately as n(z)~(1+z)^(-1). This conclusion holds irrespective of the input Mstar-Mhalo relation, the assumed stellar initial mass function, or even the chosen level of adiabatic contraction in the model. Secondary contributors to the observed redshift evolution of gamma' may come from an increased contribution at higher redshifts of adiabatic contraction and/or bottom-light stellar initial mass functions. The strong lensing selection effects we have simulated seem not to contribute to this effect. A steadily increasing S\'{e}rsic index with cosmic time is supported by independent observations, though it is not yet clear whether cosmological hierarchical models (e.g., mergers) are capable of reproducing such a fast and sharp evolution.

We describe the Sloan Digital Sky Survey IV (SDSS-IV), a project encompassing three major spectroscopic programs. The Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2) is observing hundreds of thousands of Milky Way stars at high resolution and high signal-to-noise ratio in the near-infrared. The Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey is obtaining spatially-resolved spectroscopy for thousands of nearby galaxies (median redshift of z = 0.03). The extended Baryon Oscillation Spectroscopic Survey (eBOSS) is mapping the galaxy, quasar, and neutral gas distributions between redshifts z = 0.6 and 3.5 to constrain cosmology using baryon acoustic oscillations, redshift space distortions, and the shape of the power spectrum. Within eBOSS, we are conducting two major subprograms: the SPectroscopic IDentification of eROSITA Sources (SPIDERS), investigating X-ray AGN and galaxies in X-ray clusters, and the Time Domain Spectroscopic Survey (TDSS), obtaining spectra of variable sources. All programs use the 2.5-meter Sloan Foundation Telescope at Apache Point Observatory; observations there began in Summer 2014. APOGEE-2 also operates a second near-infrared spectrograph at the 2.5-meter du Pont Telescope at Las Campanas Observatory, with observations beginning in early 2017. Observations at both facilities are scheduled to continue through 2020. In keeping with previous SDSS policy, SDSS-IV provides regularly scheduled public data releases; the first one, Data Release 13, was made available in July 2016.

Recent work has confirmed that the masses of supermassive black holes, estimated from scaling relations with global properties such as the stellar masses of their host galaxies, may be biased high. Much of this may be caused by the requirement that the gravitational sphere of influence of the black hole must be resolved for the black-hole mass to be reliably estimated. We revisit this issue by using a comprehensive galaxy evolution semi-analytic model, which self-consistently evolves supermassive black holes from high-redshift seeds via gas accretion and mergers, and also includes AGN feedback. Once tuned to reproduce the (mean) correlation of black-hole mass with velocity dispersion, the model is unable to also account for the correlation with stellar mass. This behaviour is independent of the model's parameters, thus suggesting an internal inconsistency in the data. The predicted distributions, especially at the low-mass end, are also much broader than observed. However, if selection effects are included, the model's predictions tend to align with the observations. We also demonstrate that the correlations between the residuals of the local scaling relations are more effective than the scaling relations themselves at constraining AGN feedback models. In fact, we find that our semi-analytic model, while in apparent broad agreement with the scaling relations when accounting for selection biases, yields very weak correlations between their residuals at fixed stellar mass, in stark contrast with observations. This problem persists when changing the AGN feedback strength, and is also present in the $z\sim 0$ outputs of the hydrodynamic cosmological simulation Horizon-AGN, which includes state-of-the-art treatments of AGN feedback. This suggests that current AGN feedback models may be too weak or are simply not capturing the effect of the black hole on the stellar velocity dispersion.

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 carry out a systematic investigation of the total mass density profile of massive (Mstar>2e11 Msun) early-type galaxies and its dependence on galactic properties and host halo mass with the aid of a variety of lensing/dynamical data and large mock galaxy catalogs. The latter are produced via semi-empirical models that, by design, are based on just a few basic input assumptions. Galaxies, with measured stellar masses, effective radii and S\'{e}rsic indices, are assigned, via abundance matching relations, host dark matter halos characterized by a typical LCDM profile. Our main results are as follows: (i) In line with observational evidence, our semi-empirical models naturally predict that the total, mass-weighted density slope at the effective radius gamma' is not universal, steepening for more compact and/or massive galaxies, but flattening with increasing host halo mass. (ii) Models characterized by a Salpeter or variable initial mass function and uncontracted dark matter profiles are in good agreement with the data, while a Chabrier initial mass function and/or adiabatic contractions/expansions of the dark matter halos are highly disfavored. (iii) Currently available data on the mass density profiles of very massive galaxies (Mstar>1e12 Msun), with Mhalo>3e14 Msun, favor instead models with a stellar profile flatter than a S\'{e}rsic one in the very inner regions (r<3-5 kpc), and a cored NFW or Einasto dark matter profile with median halo concentration a factor of ~2 or <1.3, respectively, higher than those typically predicted by N-body numerical simulations.

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

We extend the comparison between the set of local galaxies having dynamically measured black holes with galaxies in the Sloan Digital Sky Survey (SDSS). We first show that the most up-to-date local black hole samples of early-type galaxies with measurements of effective radii, luminosities, and S\'ersic indices of the bulges of their host galaxies, have dynamical mass and S\'ersic index distributions consistent with those of SDSS early-type galaxies of similar bulge stellar mass. The host galaxies of local black hole samples thus do not appear structurally different from SDSS galaxies, sharing similar dynamical masses, light profiles and light distributions. Analysis of the residuals reveals that velocity dispersion is more fundamental than S\'ersic index n in the scaling relations between black holes and galaxies. Indeed, residuals with S\'ersic index could be ascribed to the (weak) correlation with bulge mass or even velocity dispersion. Finally, targetted Monte Carlo simulations that include the effects of the sphere of influence of the black hole, and tuned to reproduce the observed residuals and scaling relations in terms of velocity dispersion and stellar mass, show that, at least for galaxies with Mbulge > 1e10 Msun and n>5, the observed mean black hole mass at fixed S\'ersic index is biased significantly higher than the intrinsic value.

We compare the set of local galaxies having dynamically measured black holes with a large, unbiased sample of galaxies extracted from the Sloan Digital Sky Survey. We confirm earlier work showing that the majority of black hole hosts have significantly higher velocity dispersions sigma than local galaxies of similar stellar mass. We use Monte-Carlo simulations to illustrate the effect on black hole scaling relations if this bias arises from the requirement that the black hole sphere of influence must be resolved to measure black hole masses with spatially resolved kinematics. We find that this selection effect artificially increases the normalization of the Mbh-sigma relation by a factor of at least ~3; the bias for the Mbh-Mstar relation is even larger. Our Monte Carlo simulations and analysis of the residuals from scaling relations both indicate that sigma is more fundamental than Mstar or effective radius. In particular, the Mbh-Mstar relation is mostly a consequence of the Mbh-sigma and sigma-Mstar relations, and is heavily biased by up to a factor of 50 at small masses. This helps resolve the discrepancy between dynamically-based black hole-galaxy scaling relations versus those of active galaxies. Our simulations also disfavour broad distributions of black hole masses at fixed sigma. Correcting for this bias suggests that the calibration factor used to estimate black hole masses in active galaxies should be reduced to values of fvir~1. Black hole mass densities should also be proportionally smaller, perhaps implying significantly higher radiative efficiencies/black hole spins. Reducing black hole masses also reduces the gravitational wave signal expected from black hole mergers.

[abridged] We quantify the evolution of the stellar mass functions of star-forming and quiescent galaxies as a function of morphology from $z\sim 3$ to the present. Our sample consists of ~50,000 galaxies in the CANDELS fields ($\sim880$ $arcmin^2$), which we divide into four main morphological types, i.e. pure bulge dominated systems, pure spiral disk dominated, intermediate 2-component bulge+disk systems and irregular disturbed galaxies. Our main results are: Star-formation: At $z\sim 2$, 80\% of the stellar mass density of star-forming galaxies is in irregular systems. However, by $z\sim 0.5$, irregular objects only dominate at stellar masses below $10^9M\odot$. A majority of the star-forming irregulars present at $z\sim 2$ undergo a gradual transformation from disturbed to normal spiral disk morphologies by $z\sim 1$ without significant interruption to their star-formation. Rejuvenation after a quenching event does not seem to be common except perhaps for the most massive objects. Quenching: We confirm that galaxies reaching a stellar mass of $M_*\sim10^{10.8}M_\odot$ ($M^*$) tend to quench. Also, quenching implies the presence of a bulge: the abundance of massive red disks is negligible at all redshifts over 2~dex in stellar mass. However the dominant quenching mechanism evolves. At $z>2$, the SMF of quiescent galaxies above $M^*$ is dominated by compact spheroids. Quenching at this early epoch destroys the disk and produces a compact remnant unless the star-forming progenitors at even higher redshifts are significantly more dense. At $1<z<2$, the majority of newly quenched galaxies are disks with a significant central bulge. This suggests that mass-quenching at this epoch starts from the inner parts and preserves the disk. At $z<1$, the high mass end of the passive SMF is globally in place and the evolution mostly happens at stellar masses below $10^{10}M_\odot$.

Supermassive black hole -- host galaxy relations are key to the computation of the expected gravitational wave background (GWB) in the pulsar timing array (PTA) frequency band. It has been recently pointed out that standard relations adopted in GWB computations are in fact biased-high. We show that when this selection bias is taken into account, the expected GWB in the PTA band is a factor of about three smaller than previously estimated. Compared to other scaling relations recently published in the literature, the median amplitude of the signal at $f=1$yr$^{-1}$ drops from $1.3\times10^{-15}$ to $4\times10^{-16}$. Although this solves any potential tension between theoretical predictions and recent PTA limits without invoking other dynamical effects (such as stalling, eccentricity or strong coupling with the galactic environment), it also makes the GWB detection more challenging.

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 catalog of visual like H-band morphologies of $\sim50.000$ galaxies ($H_{f160w}<24.5$) in the 5 CANDELS fields (GOODS-N, GOODS-S, UDS, EGS and COSMOS). Morphologies are estimated with Convolutional Neural Networks (ConvNets). The median redshift of the sample is $<z>\sim1.25$. The algorithm is trained on GOODS-S for which visual classifications are publicly available and then applied to the other 4 fields. Following the CANDELS main morphology classification scheme, our model retrieves the probabilities for each galaxy of having a spheroid, a disk, presenting an irregularity, being compact or point source and being unclassifiable. ConvNets are able to predict the fractions of votes given a galaxy image with zero bias and $\sim10\%$ scatter. The fraction of miss-classifications is less than $1\%$. Our classification scheme represents a major improvement with respect to CAS (Concentration-Asymmetry-Smoothness)-based methods, which hit a $20-30\%$ contamination limit at high z. The catalog is released with the present paper via the $\href{http://rainbowx.fis.ucm.es/Rainbow_navigator_public}{Rainbow\,database}$

[abridged] We quantify the morphological evolution of z~0 massive galaxies ($M*/M_\odot\sim10^{11}$) from z~3 in the 5 CANDELS fields. The progenitors are selected using abundance matching techniques to account for the mass growth. The morphologies strongly evolve from z~3. At z<1, the population matches the massive end of the Hubble sequence, with 30% of spheroids, 50% of galaxies with equally dominant disk and bulge components and 20% of disks. At z~2-3 there is a majority of irregular systems (~60-70%) with still 30% of spheroids. We then analyze the SFRs, gas fractions and structural properties for the different morphologies independently. Our results suggest two distinct channels for the growth of bulges in massive galaxies. Around 30-40% were already bulges at z~2.5, with low average SFRs and gas-fractions (10-15%), high Sersic indices (n>3-4) and small effective radii ($R_e$~1 kpc) pointing towards an early formation through gas-rich mergers or VDI. Between z~ 2.5 and z~0, they rapidly increase their size by a factor of ~4-5, become all passive but their global morphology remains unaltered. The structural evolution is independent of the gas fractions, suggesting that it is driven by ex-situ events. The remaining 60% experience a gradual morphological transformation, from clumpy disks to more regular bulge+disks systems, essentially happening at z>1. It results in the growth of a significant bulge component (n~3) for 2/3 of the systems possibly through the migration of clumps while the remaining 1/3 keeps a rather small bulge (n~1.5-2). The transition phase between disturbed and relaxed systems and the emergence of the bulge is correlated with a decrease of the star formation activity and the gas fractions. The growth of the effective radii scales roughly with $H(z)^{-1}$ and it is therefore consistent with the expected growth of disks in galaxy haloes.

The mass and structural evolution of massive galaxies is one of the hottest topics in galaxy formation. This is because it may reveal invaluable insights into the still debated evolutionary processes governing the growth and assembly of spheroids. However, direct comparison between models and observations is usually prevented by the so-called "progenitor bias", i.e., new galaxies entering the observational selection at later epochs, thus eluding a precise study of how pre-existing galaxies actually evolve in size. To limit this effect, we here gather data on high-redshift brightest group and cluster galaxies, evolve their (mean) host halo masses down to z=0 along their main progenitors, and assign as their "descendants" local SDSS central galaxies matched in host halo mass. At face value, the comparison between high redshift and local data suggests a noticeable increase in stellar mass of a factor of >2 since z~1, and of >2.5 in mean effective radius. We then compare the inferred stellar mass and size growth with those predicted by hierarchical models for central galaxies, selected at high redshifts to closely match the halo and stellar mass bins as in the data. Only hierarchical models characterized by very limited satellite stellar stripping and parabolic orbits are capable of broadly reproducing the stellar mass and size increase of a factor ~2-4 observed in cluster galaxies since z ~1. The predicted, average (major) merger rate since z~1 is in good agreement with the latest observational estimates.

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.

We compare state-of-the-art semi-analytic models of galaxy formation as well as advanced sub-halo abundance matching models with a large sample of early-type galaxies from SDSS at z < 0.3. We focus our attention on the dependence of median sizes of central galaxies on host halo mass. The data do not show any difference in the structural properties of early-type galaxies with environment, at fixed stellar mass. All hierarchical models considered in this work instead tend to predict a moderate to strong environmental dependence, with the median size increasing by a factor of about 1.5-3 when moving from low to high mass host haloes. At face value the discrepancy with the data is highly significant, especially at the cluster scale, for haloes above log Mhalo > 14. The convolution with (correlated) observational errors reduces some of the tension. Despite the observational uncertainties, the data tend to disfavour hierarchical models characterized by a relevant contribution of disc instabilities to the formation of spheroids, strong gas dissipation in (major) mergers, short dynamical friction timescales, and very short quenching timescales in infalling satellites. We also discuss a variety of additional related issues, such as the slope and scatter in the local size-stellar mass relation, the fraction of gas in local early-type galaxies, and the general predictions on satellite galaxies.

We assemble a statistical set of global mass models for ~2,000 nearly spherical SDSS galaxies at a mean redshift of 0.12 based on their aperture velocity dispersions and newly derived luminosity profiles in conjunction with published velocity dispersion profiles and empirical properties and relations of galaxy and halo parameters. When two-component (i.e. stellar plus dark) mass models are fitted to the SDSS aperture velocity dispersions, the predicted velocity dispersion profile (VP) slopes within the effective radius R_eff match well the distribution in observed elliptical galaxies. In contrast, the single-component models cannot reproduce the VP slope distribution. From a number of input variations the models exhibit for the radial range 0.1 R_eff < r < R_eff a tight correlation <gamma_e>=(1.865+/-0.008)+(-4.93+/-0.15)<eta> where <gamma_e> is the mean slope absolute value of the total mass density and <eta> is the mean slope of the velocity dispersion profile, which leads to a super-isothermal <gamma_e> = 2.15+/-0.04 for <eta>=-0.058+/-0.008 in observed elliptical galaxies. Furthermore, the successful two-component models appear to imply a typical slope curvature pattern in the total mass profile because for the observed steep luminosity (stellar mass) profile and the weak lensing inferred halo profile at large radii a total mass profile with monotonically varying slope would require too high DM density in the optical region giving rise to too large aperture velocity dispersion and too shallow VP.

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.

[abridged] The mass-size relation of early-type galaxies (ETGs) has been largely studied in the last years to probe the mass assembly of the most massive objects in the Universe. In this paper, we focus on the mass-size relation of quiescent massive ETGs (Mstar/Msol > 3*10^10) living in massive clusters (M200 ~ 10^14 Mstar) at 0.8< z <1.5, as compared to those living in the field at the same epoch. Our sample contains ~ 400 ETGs in clusters and the same number in the field. Therefore, our sample is approximately an order of magnitude larger than previous studies in the same redshift range for galaxy clusters. We find that ETGs living in clusters are between ~30-50% larger than galaxies with the same stellar mass residing in the field. We parametrize the size using the mass-normalized size, gamma=Re/Mstar^0.57. The gamma distributions in both environments peak at the same position but the distributions in clusters are more skewed towards larger sizes. Since this size difference is not observed in the local Universe, the size evolution at fixed stellar mass from z~1.5 of cluster galaxies is less steep ((1+z)-0.53pm0.04) than the evolution of field galaxies ((1+z)-0.92pm0.04). The size difference seems to be essentially driven by the galaxies residing in the clusters cores (R<0.5*R200). If part of the size evolution is due to mergers, the difference we see between cluster and field galaxies could be due to higher merger rates in clusters at higher redshift, probably during the formation phase of the clusters when velocity dispersions are lower. We cannot exclude however that the difference is driven by newly quenched galaxies which are quenched more efficiently in clusters. The implications of these results for the hierarchical growth of ETGs will be discussed in a companion paper.

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.

We present the design and performance of the multi-object fiber spectrographs for the Sloan Digital Sky Survey (SDSS) and their upgrade for the Baryon Oscillation Spectroscopic Survey (BOSS). Originally commissioned in Fall 1999 on the 2.5-m aperture Sloan Telescope at Apache Point Observatory, the spectrographs produced more than 1.5 million spectra for the SDSS and SDSS-II surveys, enabling a wide variety of Galactic and extra-galactic science including the first observation of baryon acoustic oscillations in 2005. The spectrographs were upgraded in 2009 and are currently in use for BOSS, the flagship survey of the third-generation SDSS-III project. BOSS will measure redshifts of 1.35 million massive galaxies to redshift 0.7 and Lyman-alpha absorption of 160,000 high redshift quasars over 10,000 square degrees of sky, making percent level measurements of the absolute cosmic distance scale of the Universe and placing tight constraints on the equation of state of dark energy. The twin multi-object fiber spectrographs utilize a simple optical layout with reflective collimators, gratings, all-refractive cameras, and state-of-the-art CCD detectors to produce hundreds of spectra simultaneously in two channels over a bandpass covering the near ultraviolet to the near infrared, with a resolving power R = \lambda/FWHM ~ 2000. Building on proven heritage, the spectrographs were upgraded for BOSS with volume-phase holographic gratings and modern CCD detectors, improving the peak throughput by nearly a factor of two, extending the bandpass to cover 360 < \lambda < 1000 nm, and increasing the number of fibers from 640 to 1000 per exposure. In this paper we describe the original SDSS spectrograph design and the upgrades implemented for BOSS, and document the predicted and measured performances.

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.

With SDSS galaxy data and halo data from up-to-date N-body simulations we construct a semi-empirical catalog (SEC) of early-type systems by making a self-consistent bivariate statistical match of stellar mass (M_star) and velocity dispersion (sigma) with halo virial mass (M_vir). We then assign stellar mass profile and velocity dispersion profile parameters to each system in the SEC using their observed correlations with M_star and sigma. Simultaneously, we solve for dark matter density profile of each halo using the spherical Jeans equation. The resulting dark matter density profiles deviate in general from the dissipationless profile of NFW or Einasto and their mean inner density slope and concentration vary systematically with M_vir. Statistical tests of the distribution of profiles at fixed M_vir rule out the null hypothesis that it follows the distribution predicted by N-body simulations for M_vir ~< 10^{13.5-14.5} M_solar. These dark matter profiles imply that dark matter density is, on average, enhanced significantly in the inner region of halos with M_vir ~< 10^{13.5-14.5} M_solar supporting halo contraction. The main characteristics of halo contraction are: (1) the mean dark matter density within the effective radius has increased by a factor varying systematically up to ~ 3-4 at M_vir = 10^{12} M_solar, and (2) the inner density slope has a mean of <alpha> ~ 1.3 with rho(r) ~ r^{-alpha} and a halo-to-halo rms scatter of rms(alpha) ~ 0.4-0.5 for 10^{12} M_solar ~< M_vir ~< 10^{13-14} M_solar steeper than the NFW profile (alpha=1). Based on our results we predict that halos of nearby elliptical and lenticular galaxies can, in principle, be promising targets for gamma-ray emission from dark matter annihilation.

The coefficients a and b of the Fundamental Plane relation R ~ Sigma^a I^b depend on whether one minimizes the scatter in the R direction or orthogonal to the Plane. We provide explicit expressions for a and b (and confidence limits) in terms of the covariances between logR, logSigma and logI. Our analysis is more generally applicable to any other correlations between three variables: e.g., the color-magnitude-Sigma relation, the L-Sigma-Mbh relation, or the relation between the X-ray luminosity, Sunyaev-Zeldovich decrement and optical richness of a cluster, so we provide IDL code which implements these ideas, and we show how our analysis generalizes further to correlations between more than three variables. We show how to account for correlated errors and selection effects, and quantify the difference between the direct, inverse and orthogonal fit coefficients. We show that the three vectors associated with the Fundamental Plane can all be written as simple combinations of a and b because the distribution of I is much broader than that of Sigma, and Sigma and I are only weakly correlated. Why this should be so for galaxies is a fundamental open question about the physics of early-type galaxy formation. If luminosity evolution is differential, and Rs and Sigmas do not evolve, then this is just an accident: Sigma and I must have been correlated in the past. On the other hand, if the (lack of) correlation is similar to that at the present time, then differential luminosity evolution must have been accompanied by structural evolution. A model in which the luminosities of low-L galaxies evolve more rapidly than do those of higher-L galaxies is able to produce the observed decrease in a (by a factor of 2 at z~1) while having b decrease by only about 20 percent. In such a model, the Mdyn/L ratio is a steeper function of Mdyn at higher z.

There is mounting evidence that a significant fraction of Black Holes (BHs) today live in late-type galaxies, including bulge-less galaxies and those hosting pseudobulges, and are significantly undermassive with respect to the scaling relations followed by their counterpart BHs in classical bulges of similar stellar (or even bulge) mass. Here we discuss the predictions of two state-of-the-art hierarchical galaxy formation models in which BHs grow via mergers and, in one, also via disk instability. Our aim is to understand if the wealth of new data on local BH demography is consistent with standard models. We follow the merger trees of representative subsamples of BHs and compute the fractional contributions of different processes to the final BH mass. We show that the model in which BHs always closely follow the growth of their host bulges, also during late disk instabilities (i.e., bars), produces too narrow a distribution of BHs at fixed stellar mass to account for the numerous low-mass BHs now detected in later-type galaxies. Models with a looser connection between BH growth and bar instability instead predict the existence of a larger number of undermassive BHs, in better agreement with the observations. The scatter in the updated local BH-bulge mass relation (with no restriction on galaxy type) appears to be quite large when including later-type systems, but it can still be managed to be reproduced within current hierarchical models. However, the fuelling of BHs during the late bar-instability mode needs to be better quantified/improved to properly fit the data. We conclude discussing how the possibly large number of BHs in later type galaxies demands for an in-depth revision of the local BH mass function and its modelling.

For early-type galaxies, the correlations between stellar mass and size, velocity dispersion, surface brightness, color, axis ratio and color-gradient all indicate that two mass scales, M* = 3 x 10^10 Msun and M* = 2 x 10^11 Msun, are special. The smaller scale could mark the transition between wet and dry mergers, or it could be related to the interplay between SN and AGN feedback, although quantitative measures of this transition may be affected by morphological contamination. At the more massive scale, mean axis ratios and color gradients are maximal, and above it, the colors are redder, the sizes larger and the velocity dispersions smaller than expected based on the scaling at lower M*. In contrast, the color-sigma relation, and indeed, most scaling relations with sigma, are not curved: they are well-described by a single power law, or in some cases, are almost completely flat. When major dry mergers change masses, sizes, axis ratios and color gradients, they are expected to change the colors or velocity dispersions much less. Therefore, the fact that scaling relations at sigma > 150 km/s show no features, whereas the size-M*, b/a-M*, color-M* and color gradient-M* relations do, suggests that M* = 2 x 10^11 Msun is the scale above which major dry mergers dominate the assembly histories of early-type galaxies.

The color-magnitude relation of early-type galaxies differs slightly but significantly from a pure power-law, curving downwards at low and upwards at large luminosities (Mr>-20.5 and Mr<-22.5). This remains true of the color-size relation, and is even more apparent with stellar mass (M* < 3x10^10 Msun and M* > 2x10^11 Msun). The upwards curvature at the massive end does not appear to be due to stellar population effects. In contrast, the color-sigma relation is well-described by a single power law. Since major dry mergers change neither the colors nor sigma, but they do change masses and sizes, the clear features observed in the scaling relations with M*, but not with sigma > 150 km/s, suggest that M* > 2x10^11 Msun is the scale above which major dry mergers dominate the assembly history. We discuss three models of the merger histories since z ~ 1 which are compatible with our measurements. In all three models, dry mergers are responsible for the flattening of the color-M* relation at M* > 3x10^10 Msun - wet mergers only matter at smaller masses. At M* > 2 x 10^11 Msun, the merger histories in one model are dominated by major rather than minor dry mergers, as suggested by the axis ratio and color gradient trends. In another, although both major and minor mergers occur at the high mass end, the minor mergers contribute primarily to the formation of the ICL, rather than to the mass growth of the central massive galaxy. A final model assumes that the reddest objects were assembled by a mix of major and minor dry mergers.

We present an automated morphological classification in 4 types (E,S0,Sab,Scd) of ~700.000 galaxies from the SDSS DR7 spectroscopic sample based on support vector machines. The main new property of the classification is that we associate to each galaxy a probability of being in the four morphological classes instead of assigning a single class. The classification is therefore better adapted to nature where we expect a continuos transition between different morphological types. The algorithm is trained with a visual classification and then compared to several independent visual classifications including the Galaxy Zoo first release catalog. We find a very good correlation between the automated classification and classical visual ones. The compiled catalog is intended for use in different applications and can be downloaded at http://gepicom04.obspm.fr/sdss_morphology/Morphology_2010.html and soon from the CasJobs database.

We analyze a sample of 23 supermassive elliptical galaxies (central velocity dispersion larger than 330 km s-1), drawn from the SDSS. For each object, we estimate the dynamical mass from the light profile and central velocity dispersion, and compare it with the stellar mass derived from stellar population models. We show that these galaxies are dominated by luminous matter within the radius for which the velocity dispersion is measured. We find that the sizes and stellar masses are tightly correlated, with Re ~ M*^{1.1}$, making the mean density within the de Vaucouleurs radius a steeply declining function of M*: rho_e ~ M*^{-2.2}. These scalings are easily derived from the virial theorem if one recalls that this sample has essentially fixed (but large) sigma_0. In contrast, the mean density within 1 kpc is almost independent of M*, at a value that is in good agreement with recent studies of z ~ 2 galaxies. The fact that the mass within 1 kpc has remained approximately unchanged suggests assembly histories that were dominated by minor mergers -- but we discuss why this is not the unique way to achieve this. Moreover, the total stellar mass of the objects in our sample is typically a factor of ~ 5 larger than that in the high redshift (z ~ 2) sample, an amount which seems difficult to achieve. If our galaxies are the evolved objects of the recent high redshift studies, then we suggest that major mergers were required at z > 1.5, and that minor mergers become the dominant growth mechanism for massive galaxies at z < 1.5.

[abridged] Massive, passively evolving galaxies at redshifts z>1 exhibit on the average physical sizes smaller by factors ~3 than local early type galaxies (ETGs) endowed with the same stellar mass. Small sizes are in fact expected on theoretical grounds, if dissipative collapse occurs. Recent results show that the size evolution at z<1 is limited to less than 40%, while most of the evolution occurs at z>1, where both compact and already extended galaxies are observed and the scatter in size is remarkably larger than locally. The presence at high z of a significant number of ETGs with the same size as their local counterparts as well as of ETGs with quite small size, points to a timescale to reach the new, expanded equilibrium configuration of less than the Hubble time. We demonstrate that the projected mass of compact, high-z galaxies and that of local ETGs within the *same physical radius*, the nominal half-luminosity radius of high-z ETGs, differ substantially, in that the high-z ETGs are on the average significantly denser. We propose that quasar activity, which peaks at z~2, can remove large amounts of gas from central galaxy regions on a timescale shorter than of the dynamical one, triggering a puffing up of the stellar component at constant stellar mass; in this case the size increase goes together with a decrease of the central mass. The size evolution is expected to parallel that of the quasars and the inverse hierarchy, or downsizing, seen in the quasar evolution is mirrored in the size evolution.

We examine the colour-magnitude relation of approximately 5000 Brightest Cluster Galaxies (BCGs) in the Sloan Digital Sky Survey, and compare with non-BCG E/S0 galaxies. The colour-magnitude and colour-sigma (velocity dispersion) relations are flatter in slope (by a factor of about 2) for BCGs than for non-BCG E/S0s, and the BCGs also tend to be redder by 0.01 magnitudes in g-r. We investigate radial colour gradients in both samples, using the ratio of the de Vaucouleurs radii in the g and r bands. We find BCGs have significantly flatter (by 23%) mean colour gradients than other high luminosity E/S0s. In early-type galaxies, the colour gradients are strongest at intermediate luminosities of Mr=-22. Colour gradients in E/S0s increase with radius (up to 10kpc) and are negatively correlated with 10sigma + Mr (velocity dispersion relative to luminosity) and with mass density. The gradients also tend to decrease with increasing stellar age. These trends are weak or not seen in BCGs, in which the mean colour gradient is low whatever the other properties. We discuss possible explanations, which involve a greater amount of dry merging in the formation history of the BCGs.

We discuss how the effective radius Phi(Re) function (ERF) recently worked out by Bernardi et al. (2009) represents a new testbed to improve the current understanding of Semi-analytic Models of Galaxy formation. In particular, we here show that a detailed hierarchical model of structure formation can broadly reproduce the correct peak in the size distribution of local early-type galaxies, although it significantly overpredicts the number of very compact and very large galaxies. This in turn is reflected in the predicted size-mass relation, much flatter than the observed one, due to too large (~3 kpc) low-mass galaxies (<10^11 \msun), and to a non-negligible fraction of compact (< 0.5-1 kpc) and massive galaxies (> 10^11 \msun). We also find that the latter discrepancy is smaller than previously claimed, and limited to only ultracompact (Re < 0.5 kpc) galaxies when considering elliptical-dominated samples. We explore several causes behind these effects. We conclude that the former problem might be linked to the initial conditions, given that large and low-mass galaxies are present at all epochs in the model. The survival of compact and massive galaxies might instead be linked to their very old ages and peculiar merger histories. Overall, knowledge of the galactic stellar mass {\em and} size distributions allows a better understanding of where and how to improve models.

(Abridged) We present a detailed study of the stellar populations of a volume-limited sample of early-type galaxies from SDSS, across a range of environments -- defined as the mass of the host dark matter halo. The stellar populations are explored through the SDSS spectra, via projection onto a set of two spectral vectors determined from Principal Component Analysis. We find the velocity dispersion of the galaxy to be the main driver behind the different star formation histories of early-type galaxies. However, environmental effects are seen to play a role (although minor). Galaxies populating the lowest mass halos have stellar populations on average ~1Gyr younger than the rest of the sample. The fraction of galaxies with small amounts of recent star formation is also seen to be truncated when occupying halos more massive than 3E13Msun. The sample is split into satellite and central galaxies for a further analysis of environment. Satellites are younger than central galaxies of the same stellar mass. The younger satellite galaxies in 6E12Msun halos have stellar populations consistent with the central galaxies found in the lowest mass halos of our sample (i.e. 1E12Msun). This result is indicative of galaxies in lower mass halos being accreted into larger halos.

We provide fits to the distribution of galaxy luminosity, size, velocity dispersion and stellar mass as a function of concentration index C_r and morphological type in the SDSS. We also quantify how estimates of the fraction of `early' or `late' type galaxies depend on whether the samples were cut in color, concentration or light profile shape, and compare with similar estimates based on morphology. Our fits show that Es account for about 20% of the r-band luminosity density, rho_Lr, and 25% of the stellar mass density, rho_*; including S0s and Sas increases these numbers to 33% and 40%, and 50% and 60%, respectively. Summed over all galaxy types, we find rho_* ~ 3 * 10^8 M_Sun Mpc^{-3} at z ~ 0. This is in good agreement with expectations based on integrating the star formation history. However, compared to most previous work, we find an excess of objects at large masses, up to a factor of ~ 10 at M_* ~ 5*10^{11} M_Sun. The stellar mass density further increases at large masses if we assume different IMFs for Es and spiral galaxies, as suggested by some recent chemical evolution models, and results in a better agreement with the dynamical mass function. We also show that the trend for ellipticity to decrease with luminosity is primarily because the E/S0 ratio increases at large L. However, the most massive galaxies, M_* > 5 * 10^{11} M_Sun, are less concentrated and not as round as expected if one extrapolates from lower L, and they are not well-fit by pure deVaucouleur laws. This suggests formation histories with recent radial mergers. Finally, we show that the age-size relation is flat for Es of fixed dynamical mass, but, at fixed M_dyn, S0s and Sas with large sizes tend to be younger. Explaining this difference between E and S0 formation is a new challenge for models of early-type galaxy formation.

In a sample of about 45,700 early-type galaxies extracted from SDSS, we find that the shape, normalization, and dispersion around the mean size-stellar mass relation is the same for young and old systems, provided the stellar mass is greater than 3*10^10 Msun. This is difficult to reproduce in pure passive evolution models, which generically predict older galaxies to be much more compact than younger ones of the same stellar mass. However, this aspect of our measurements is well reproduced by hierarchical models of galaxy formation. Whereas the models predict more compact galaxies at high redshifts, subsequent minor, dry mergers increase the sizes of the more massive objects, resulting in a flat size-age relation at the present time. At lower masses, the models predict that mergers are less frequent, so that the expected anti-correlation between age and size is not completely erased. This is in good agreement with our data: below 3*10^10 Msun, the effective radius R_e is a factor of ~2 lower for older galaxies. These successes of the models are offset by the fact that the predicted sizes have other serious problems, which we discuss.

We select a sample of 70378 E/S0 (early-type) galaxies at 0<z<0.36 from the Sloan Digital Sky Survey, excluding disk and star-forming galaxies. We estimate g and r magnitudes in the observer- and rest-frames directly from the SDSS DR6 spectra; this provides an object-by-object estimate of the k-correction. We use the k-corrections from the spectra to study the evolution of the rest-frame colour-magnitude relation (CMR) and colour-sigma (velocity dispersion) relation (CsigmaR). The evolution is very sensitive to the k-correction. Both the CMR and CsigmaR evolve blueward with increasing redshift, approximately in agreement with passive evolution models with age approx. 12 Gyr. The slope and zero-point of the CMR depend on whether colours were defined in fixed physical or angular apertures, a consequence of the fact that the centers of these objects tend to be redder: the relation is steeper for fixed angular apertures. One the other hand, the CsigmaR slope does not show this dependence on the aperture in which the colour was defined, suggesting that colour gradients are correlated with residuals from the sigma-M_r relation. As these residuals are age indicators, our findings suggest that colour gradients correlate with the age of the stellar population.

From a sample of ~50000 early-type galaxies from the SDSS, we measured the traditional Fundamental Plane in four bands. We then replaced luminosity with stellar mass, and measured the "stellar mass" FP. The FP steepens slightly as one moves from shorter to longer wavelengths: the orthogonal fit has slope 1.40 in g and 1.47 in z. The FP is thinner at longer wavelengths: scatter is 0.062 dex in g, 0.054 dex in z. The scatter is larger at small galaxy sizes/masses; at large masses measurement errors account for essentially all of the observed scatter. The FP steepens further when luminosity is replaced with stellar mass, to slope ~ 1.6. The intrinsic scatter also reduces further, to 0.048 dex. Since color and stellar mass-to-light ratio are closely related, this explains why color can be thought of as the fourth FP parameter. However, the slope of the stellar mass FP remains shallower than the value of 2 associated with the virial theorem. This is because the ratio of dynamical to stellar mass increases at large masses as M_d^0.17. The face-on view of the stellar mass kappa-space suggests that there is an upper limit to the stellar density for a given dynamical mass, and this decreases at large masses: M_*/R_e^3 ~ M_d^-4/3. We also study how the estimated coefficients a and b of the FP are affected by other selection effects (e.g. excluding small sigma biases a high; excluding fainter L biases a low). These biases are seen in FPs which have no intrinsic curvature, so the observation that a and b scale with L and sigma is not, by itself, evidence that the Plane is warped. We show that the FP appears to curve sharply downwards at the small mass end, and more gradually downwards towards larger masses. Whereas the drop at small sizes is real, most of the latter effect is due to correlated errors.

We use a sample of about 48,000 SDSS early-type galaxies to show that older galaxies have smaller half-light radii re and larger velocity dispersions sigma than younger ones of the same stellar mass Mstar. We use the age-corrected luminosity Lcorr as a proxy for Mstar to minimize biases: below Lcorr~1e11 Lsun, galaxies with age ~11 Gyrs have re smaller by 40% and sigma larger by 25%, compared to galaxies that are ~4 Gyrs younger. The sizes and velocity dispersions of more luminous galaxies vary by less than 15%, whatever their age, a challenge for current galaxy formation models. A closer check reveals that the lowering in the dispersion is caused by older galaxies that show a significant departure from the re--Lcorr and sigma--Lcorr relations at high Lcorr. Such features might find an explanation in models where more massive galaxies undergo more minor mergers than less massive galaxies at late times, thus causing a break in the homology. In terms of the Fundamental Plane of early-type galaxies, the data indicate that all galaxies show a significant and similar increase in the dynamical-to-stellar mass ratio with increasing mass, independent of their age. However, older galaxies have smaller Mdyn/Mstar ratios than objects which formed more recently. These findings may suggest that lower mass galaxies and, at fixed stellar mass, higher redshift galaxies, formed from gas-richer progenitors, thus underwent more dissipation and contraction.

At the present time, SDSS early-type BCGs have larger Re than early-type galaxies of similar L, whether these other objects are in the field, or are satellites in clusters (Re ~ L for BCGs). At fixed M* and formation time, BCGs at lower z are larger and have smaller velocity dispersions, i.e. Re increases and sigma decreases with age. As a result, at z~0.25, corresponding to a lookback time of order 3 Gyrs, BCGs are smaller than their lower z counterparts by as much as 70% for the brightest BCGs: Re evolves as (1+z)^{0.85(Mr+21)}. Qualitatively similar but weaker evolution in the sizes is also seen in the bulk of the early-type population: at Mr<-22 Re evolves as (1+z)^{0.7(Mr+21)}, while at Mr>-22 the evolution is approximately (1+z)^{-0.7}, independent of Mr. The sigma-L correlation also evolves: (1+z)^{-0.2(Mr+21)} at Mr < -22 (as for the BCGs) and (1+z)^{0.2} for fainter galaxies. The Re- and sigma-M* correlations yield consistent results. These trends are most easily understood if early-type BCGs grew from many dry minor mergers rather than a few major mergers. Only in such a scenario can BCGs be the descendents of the superdense galaxies seen at z~2; major dry mergers, which increase the size in proportion to the mass, cannot bring these galaxies onto the BCG Re-M* relation at z~0. We also compared the ages and sizes of our early-type BCGs with satellites. BCGs are larger than satellites of similar L or M* at the same redshift. Although both satellites and BCGs trace the same weak age-L or age-M* relation, this can be understood by noting that BCGs are typically about 1 Gyr older than the satellites in their group, and they are about 0.5 mags more luminous. Finally, we find that the mean satellite L is approximately independent of BCG L, in agreement with recent predictions based on the L-dependence of clustering.

We select a sample of about 50,000 early-type galaxies from the Sloan Digital Sky Survey (SDSS), calibrate fitting formulae which correct for known problems with photometric reductions of extended objects, apply these corrections, and then measure a number of pairwise scaling relations in the corrected sample. We show that, because they are not seeing corrected, the use of Petrosian-based quantities in magnitude limited surveys leads to biases, and suggest that this is one reason why Petrosian-based analyses of BCGs have failed to find significant differences from the bulk of the early-type population. These biases are not present when seeing-corrected parameters derived from deVaucouleur fits are used. Most of the scaling relations we study show evidence for curvature: the most luminous galaxies have smaller velocity dispersions, larger sizes, and fainter surface brightnesses than expected if there were no curva-ture. These statements remain true if we replace luminosities with stellar masses; they suggest that dissipation is less important at the massive end. There is curvature in the dynamical to stellar mass relation as well: the ratio of dynamical to stellar mass increases as stellar mass increases, but it curves upwards from this scaling both at small and large stellar masses. In all cases, the curvature at low masses becomes apparent when the sample becomes dominated by objects with stellar masses smaller than 3 x 10^10 M_Sun. We quantify all these trends using second order polynomials; these generally provide significantly better description of the data than linear fits, except at the least luminous end.

We used the Advanced Camera for Surveys on board the Hubble Space Telescope to obtain high resolution i-band images of the centers of 23 single galaxies, which were selected because they have SDSS velocity dispersions larger than 350 km/s. The surface brightness profiles of the most luminous of these objects (M_i<-24) have well-resolved `cores' on scales of 150-1000 pc, and share similar properties to BCGs. The total luminosity of the galaxy is a better predictor of the core size than is the velocity dispersion. The correlations of luminosity and velocity dispersion with core size agree with those seen in previous studies of galaxy cores. Because of high velocity dispersions, our sample of galaxies can be expected to harbor the most massive black holes, and thus have large cores with large amounts of mass ejection. The mass-deficits inferred from core-Sersic fits to the surface-brightness profiles are approximately double the black-hole masses inferred from the M_bh-sigma relation and the same as those inferred from the M_bh-L relation. The less luminous galaxies (M_i>-23) tend to have steeper `power-law' inner profiles, higher-ellipticity, diskier isophotes, and bulge-to-total ratios of order 0.5 -- all of which suggest that they are `fast-rotators' and rotational motions could have contaminated the velocity dispersion estimate. There are obvious dust features within about 300 pc of the center in about 35% of the sample, predominantly in power-law rather than core galaxies.

We study a sample of 43 early-type galaxies, selected from the SDSS because they appeared to have velocity dispersion > 350 km/s. High-resolution photometry in the SDSS i passband using HRC-ACS on board the HST shows that just less than half of the sample is made up of superpositions of two or three galaxies, so the reported velocity dispersion is incorrect. The other half of the sample is made up of single objects with genuinely large velocity dispersions. None of these objects has sigma larger than 426 +- 30 km/s. These objects define rather different relations than the bulk of the early-type galaxy population: for their luminosities, they are the smallest, most massive and densest galaxies in the Universe. Although the slopes of the scaling relations they define are rather different from those of the bulk of the population, they lie approximately parallel to those of the bulk "at fixed sigma". These objects appear to be of two distinct types: the less luminous (M_r>-23) objects are rather flattened and extremely dense for their luminosities -- their properties suggest some amount of rotational support and merger histories with abnormally large amounts of gaseous dissipation. The more luminous objects (M_r<-23) tend to be round and to lie in or at the centers of clusters. Their properties are consistent with the hypothesis that they are BCGs. Models in which BCGs form from predominantly radial mergers having little angular momentum predict that they should be prolate. If viewed along the major axis, such objects would appear to have abnormally large sigma for their sizes, and to be abnormally round for their luminosities. This is true of the objects in our sample once we account for the fact that the most luminous galaxies (M_r<-23.5), and BCGs, become slightly less round with increasing luminosity.

We present Hobby-Eberly Telescope (HET) observations for galaxies at redshift z < 0.3 from the Sloan Digital Sky Survey (SDSS) showing large velocity dispersions while appearing to be single galaxies in HST images. The high signal-to-noise HET spectra provide more definitive velocity dispersions. The maximum velocity dispersion we find is 444 km/s. Emission-line widths in QSOs indicate that black holes can exist with masses exceeding 5 billion solar masses, implying velocity dispersions greater than 500 km/s by the local black hole mass - velocity dispersion relationship. This suggests either that QSO black hole masses are overestimated or that the black hole - bulge relationship changes at high black hole mass. The latter option is consistent with evidence that the increase in velocity dispersion with luminosity levels off for the brightest elliptical galaxies.

We utilize the local velocity dispersion function (VDF) of spheroids, together with their inferred age--distributions, to predict the VDF at higher redshifts (0<z<6), under the assumption that (i) most of the stars in each nearby spheroid formed in a single episode, and (ii) the velocity dispersion sigma remained nearly constant afterward. We assume further that a supermassive black hole (BH) forms concurrently with the stars, and within ~1 Gyr of the formation of the potential well of the spheroid, and that the relation between the mass of the BH and host velocity dispersion maintains the form M_BH ~ sigma^{beta} with beta~4, but with the normalization allowed to evolve with redshift as ~(1+z)^{alpha}. We compute the BH mass function associated with the VDF at each redshift, and compare the accumulated total BH mass density with that inferred from the integrated quasar luminosity function (LF; the so--called Soltan argument). This comparison is insensitive to the assumed duty cycle or Eddington ratio of quasar activity, and we find that the match between the two BH mass densities favors a relatively mild redshift evolution, with alpha ~ 0.26, with a positive evolution as strong as alpha>1.3 excluded at the 99% confidence level. A direct match between the characteristic BH mass in the VDF--based and quasar LF--based BH mass functions also yields a mean Eddington ratio of lambda ~ 0.5-1 that is roughly constant within 0<z<3. A strong positive evolution in the M_BH-sigma relation is still allowed by the data if galaxies increase, on average, their velocity dispersions since the moment of formation, due to dissipative processes. If we assume that the mean velocity dispersion of the host galaxies evolves as sigma(z)=sigma(0)*(1+z)^{-gamma}, we find a lower limit of gamma>0.23 for alpha>1.5. abridged

(Abridged) We apply Principal Component Analysis (PCA) to a sample of early-type galaxies from the Sloan Digital Sky Survey (SDSS) in order to infer differences in their star formation histories. Out of the first few principal components (PC), we study four which give information about stellar populations and velocity dispersion. We construct two parameters (eta and zeta) as linear combinations of PC1 and PC2. We find zeta to be most sensitive to recent episodes of star formation, and eta to be strongly dependent on the average age of the stellar populations. The distribution of the eta component of the composites appear to be indistinguishable between high and low density regions, whereas the distribution of zeta parameters have a significant skew towards lower values for galaxies in low density regions. This result suggests that galaxies in lower density environments are less likely to present weak episodes of recent star formation. In contrast, a significant number of galaxies from our high density subsample -- which includes clusters (both outer regions and centres) and groups -- underwent small but detectable recent star formation at high metallicity, in agreement with recent estimates targeting elliptical galaxies in Hickson Compact Groups and in the field (Ferreras et al.).

Black hole masses are tightly correlated with the stellar velocity dispersions of the bulges which surround them, and slightly less-well correlated with the bulge luminosity. It is common to use these correlations to estimate the expected abundance of massive black holes. This is usually done by starting from an observed distribution of velocity dispersions or luminosities and then changing variables. This procedure neglects the fact that there is intrinsic scatter in these black hole mass--observable correlations. Accounting for this scatter results in estimates of black hole abundances which are larger by almost an order of magnitude at masses >10^9 M_sun. Including this scatter is particularly important for models which seek to infer quasar lifetimes and duty cycles from the local black hole mass function. However, even when scatter has been accounted for, the M_bh-sigma relation predicts fewer massive black holes than does the M_bh-L relation. This is because the sigma-L relation in the black hole samples currently available is inconsistent with that in the SDSS sample from which the distributions of L or sigma are based: the black hole samples have smaller L for a given sigma, or larger sigma for a given L. The sigma-L relation in the black hole samples is similarly discrepant with that in other samples of nearby early-type galaxies. This suggests that current black hole samples are biased: if this is a selection rather than physical effect, then the M_bh-sigma and M_bh-L relations currently in the literature are also biased from their true values.

Early-type galaxy velocity dispersions and luminosities are correlated. The correlation estimated in local samples (< 100 Mpc) differs from that measured more recently in the SDSS. This is true even when systematics in the SDSS photometric and spectroscopic parameters have been accounted-for. We show that this is also true for the ENEAR sample if galaxy luminosities are estimated using distances which have been corrected for peculiar motions. We then show that, because the estimate of the `true' distance is derived from a correlation with velocity dispersion, in this case the D_n-sigma relation, using it in the sigma-L relation leads to an artificially tight relation with a biased slope. Making no correction for peculiar velocities results in a sigma-L relation which is very similar to that of the SDSS, although with larger scatter. We also measure the sigma-L correlation in a mock ENEAR catalog, in which the underlying galaxy sample has the same sigma-L correlation as seen in the SDSS. The mock catalog produces the same D_n-sigma relation as the data, the same biased slope when D_n-sigma distances are used to estimate luminosities, and good agreement with the input sigma-L relation when redshift is used as the distance indicator. This provides further evidence that the true sigma-L relation of ENEAR galaxies is indeed very similar to that of SDSS early-types. Our results suggest that local sigma-L relations which are based on Fundamental Plane distances should also be re-evaluated. Our findings also have important implications for black hole demographics; the best direct estimates of the masses of supermassive black holes come from local galaxies, so estimates of the black hole mass function are more safely made by working with the Mbh-sigma correlation than with Mbh-L.

It is common to estimate black hole abundances by using a measured correlation between black hole mass and another more easily measured observable such as the velocity dispersion or luminosity of the surrounding bulge. The correlation is used to transform the distribution of the observable into an estimate of the distribution of black hole masses. However, different observables provide different estimates: the Mbh-sigma relation predicts fewer massive black holes than does the Mbh-L relation. This is because the sigma-L relation in black hole samples currently available is inconsistent with that in the SDSS sample, from which the distributions of L or sigma are based: the black hole samples have smaller L for a given sigma or have larger sigma for a given L. This is true whether L is estimated in the optical or in the NIR. If this is a selection rather than physical effect, then the Mbh-sigma and Mbh-L relations currently in the literature are also biased from their true values. We provide a framework for describing the effect of this bias. We then combine it with a model of the bias to make an estimate of the true intrinsic relations. While we do not claim to have understood the source of the bias, our simple model is able to reproduce the observed trends. If we have correctly modeled the selection effect, then our analysis suggests that the bias in the <Mbh|sigma> relation is likely to be small, whereas the <Mbh|L> relation is biased towards predicting more massive black holes for a given luminosity. In addition, it is likely that the Mbh-L relation is entirely a consequence of more fundamental relations between Mbh and sigma, and between sigma and L. The intrinsic relation we find suggests that at fixed luminosity, older galaxies tend to host more massive black holes.

The size-luminosity relation of early-type Brightest Cluster Galaxies (BCGs), R_e ~ L^0.88, is steeper than that for the bulk of the early-type galaxy population, for which R_e ~ L^0.68. In addition, although BCGs are hardly offset from the Fundamental Plane defined by the bulk of the early-type population, they show considerably smaller scatter. The larger than expected sizes of BCGs, and the increased homogeneity, are qualitatively consistent with models which seek to explain the colors of the most massive galaxies by invoking dry dissipationless mergers, since dissipation tends to reduce the sizes of galaxies, and wet mergers which result in star formation would tend to increase the scatter in luminosity at fixed size and velocity dispersion. Furthermore, BCGs define the same g-r color-magnitude relation as the bulk of the early-type population. If BCGs formed from dry mergers, then BCG progenitors must have been red for their magnitudes, suggesting that they hosted older stellar populations than typical for their luminosities. Our findings have two other consequences. First, the R_e-L relation of the early-type galaxy population as a whole (i.e., normal plus BCG) exhibits some curvature. Some of this curvature must be a consequence of the fact that an increasing fraction of the most luminous galaxies are BCGs. The second consequence is suggested by the fact that, despite following a steeper size-luminosity relation, BCGs tend to define a tight relation between dynamical mass R_e sigma^2/G and luminosity. As consequence, we find that BCGs define a shallower sigma-L relation than the bulk of the early-type galaxy population.

We use the spectra of ~ 22,000 early-type galaxies, selected from the Sloan Digital Sky Survey, to infer the ages, metallicities and star formation histories of these galaxies. We find clear evidence of "downsizing", i.e. galaxies with larger velocity dispersion have older stellar populations. In particular, most early-type galaxies with velocity dispersion exceeding 200 km s-1 formed more than 90% of their current stellar mass at redshift z > 2.5. Therefore, star formation was suppressed around this redshift. We also show that chemical enrichment was rapid, lasting 1-2 Gyr and find evidence that [Fe/H] is sub-solar. We study the robustness of these results by comparing three different approaches: using (i) Lick absorption line indices; (ii) fitting a single-burst stellar population model to the whole spectrum (lines+continuum); and (iii) reconstructing the star formation and metallicity histories in multiple age-bins, providing a method to measure mass-weighted ages and metallicities. We find good agreement between the luminosity-weighted ages and metallicities computed with these three methods.

We use the spectra of 22,000 nearby early-type galaxies from the Sloan Digital Sky Survey (SDSS) to determine the age distribution of these galaxies as a function of their velocity dispersion sigma_v in the range 100 km/s < sigma_v < 280 km/s. We then combine the inferred age-distributions with the local abundance of spheroids, including early-type galaxies and late-type bulges, to predict the evolution of the quasar luminosity function (LF) in the redshift range 0<z<6. We make the following simple assumptions: (i) the formation of stars in each galaxy, at the epoch identified with the mean mass-weighted stellar age, is accompanied by the prompt assembly of the nuclear supermassive black hole (SMBH); (ii) the mass of the SMBH obeys the M_bh-sigma_v correlation observed in nearby galaxies; (iii) the SMBH radiates at a fraction f_Edd of the Eddington limit for a fixed duration t_Q, and is identified as a luminous quasar during this epoch, (iv) the intrinsic dispersions in the Eddington ratio and the M_bh-sigma_v relation produce a combined scatter of Delta(log L_Q) around the mean logarithmic quasar luminosity <log L_Q> at fixed sigma_v. These assumptions require that the SMBH remnants of quasars with bolometric luminosity below L_bol=10^{12.5} f_Edd L_sun reside predominantly in bulges of late type galaxies. We find that evolution of the observed quasar LF can be fit over the entire redshift range in this simple model, 0<z<6 with the choices of Delta(log L_Q)=0.6-0.9, t_Q= (6-8)x10^7 yr, and <f_Edd>=0.3-0.5. We find no evidence that any of the model parameters evolves with redshift, supporting the strong connection between the formation of stars and nuclear SMBHs in spheroids.

We study the photometric and spectral properties of 39320 early-type galaxies within the SDSS, as a function of both local environment and redshift. The distance to the nearest cluster of galaxies and the distance to the 10th nearest luminous neighbor were used to define two extremes in environment. The properties of early-type galaxies are weakly but significantly different in these two extremes: objects in dense environments are less than ~1 Gyr older and alpha-enhanced by ~0.02 relative to their counterparts of the same velocity dispersion in less dense regions, but the metallicities show no dependence on environment. This suggests that, in dense regions, the stars in early-type galaxies formed at slightly earlier times, and on a slightly shorter timescale, than in less dense regions. In addition, objects at lower redshifts are older, but have the same metallicities and alpha-enhancements as their counterparts of the same sigma at higher redshifts, as expected if the low redshift sample is a passively-aged version of the sample at higher redshifts. These significant conclusions are derived from several pieces of evidence discussed in this paper, including the FP of early-types, the observed evolution of absorption linestrengths (measured from a large ensemble of composite spectra which we provide in tabular form) over small lookback times, and comparison with the predictions of single burst (non-solar abundance) stellar population synthesis models (which also indicate that age, metallicity and alpha-enhancement all increase with velocity dispersion).

We describe the results of a search for galaxies with large (> 350 km/s) velocity dispersions. The largest systems we have found appear to be the extremes of the early-type galaxy population: compared to other galaxies with similar luminosities, they have the largest velocity dispersions and the smallest sizes. However, they are not distant outliers from the Fundamental Plane and mass-to-light scaling relations defined by the bulk of the early-type galaxy population. They may host the most massive black holes in the Universe, and their abundance and properties can be used to constrain galaxy formation models. Clear outliers from the scaling relations tend to be objects in superposition (angular separations smaller than 1 arcsec), evidence for which comes sometimes from the spectra, sometimes from the images, and sometimes from both. The statistical properties of the superposed pairs, e.g., the distribution of pair separations and velocity dispersions, can be used to provide useful information about the expected distribution of image multiplicities, separations and flux ratios due to gravitational lensing by multiple lenses, and may also constrain models of their interaction rates.

We present the "C4 Cluster Catalog", a new sample of 748 clusters of galaxies identified in the spectroscopic sample of the Second Data Release (DR2) of the Sloan Digital Sky Survey (SDSS). The C4 cluster--finding algorithm identifies clusters as overdensities in a seven-dimensional position and color space, thus minimizing projection effects which plagued previous optical clusters selection. The present C4 catalog covers ~2600 square degrees of sky with groups containing 10 members to massive clusters having over 200 cluster members with redshifts. We provide cluster properties like sky location, mean redshift, galaxy membership, summed r--band optical luminosity (L_r), velocity dispersion, and measures of substructure. We use new mock galaxy catalogs to investigate the sensitivity to the various algorithm parameters, as well as to quantify purity and completeness. These mock catalogs indicate that the C4 catalog is ~90% complete and 95% pure above M_200 = 1x10^14 solar masses and within 0.03 <=z <= 0.12. The C4 algorithm finds 98% of X-ray identified clusters and 90% of Abell clusters within 0.03 <= z <= 0.12. We show that the L_r of a cluster is a more robust estimator of the halo mass (M_200) than the line-of-sight velocity dispersion or the richness of the cluster. L_r. The final SDSS data will provide ~2500 C4 clusters and will represent one of the largest and most homogeneous samples of local clusters.

We present an analysis of the color-magnitude-velocity dispersion relation for a sample of 39320 early-type galaxies within the Sloan Digital Sky Survey. We demonstrate that the color-magnitude relation is entirely a consequence of the fact that both the luminosities and colors of these galaxies are correlated with stellar velocity dispersions. Previous studies of the color-magnitude relation over a range of redshifts suggest that the luminosity of an early-type galaxy is an indicator of its metallicity, whereas residuals in color from the relation are indicators of the luminosity-weighted age of its stars. We show that this, when combined with our finding that velocity dispersion plays a crucial role, has a number of interesting implications. First, galaxies with large velocity dispersions tend to be older (i.e., they scatter redward of the color-magnitude relation). Similarly, galaxies with large dynamical mass estimates also tend to be older. In addition, at fixed luminosity, galaxies which are smaller, or have larger velocity dispersions, or are more massive, tend to be older. Second, models in which galaxies with the largest velocity dispersions are also the most metal poor are difficult to reconcile with our data. However, at fixed velocity dispersion, galaxies have a range of ages and metallicities: the older galaxies have smaller metallicities, and vice-versa. Finally, a plot of velocity dispersion versus luminosity can be used as an age indicator: lines of constant age run parallel to the correlation between velocity dispersion and luminosity.

We have studied the morphology-density relation and morphology--cluster-centric-radius relation using a volume limited sample (0.05$<z<$0.1, $Mr^*<-$20.5) of the Sloan Digital Sky Survey (SDSS) data. We found there are two characteristic changes in both the morphology-density and the morphology-radius relations, suggesting two different mechanisms are responsible for the relations. In the sparsest regions (below 1 Mpc$^{-2}$ or outside of 1 virial radius), both relations become less noticeable, suggesting the responsible physical mechanisms for galaxy morphological change require denser environment. In the intermediate density regions (density between 1 and 6 Mpc$^{-2}$ or virial radius between 0.3 and 1), intermediate-type fractions increase toward denser regions, whereas late-disc fractions decrease. In the densest regions (above 6 Mpc$^{-2}$ or inside of 0.3 virial radius), intermediate-type fractions decrease radically and early-type fractions increase in turn. This is a contrasting result to that in intermediate regions and it suggests that yet another mechanism is more responsible for the morphological change in these regions. We also compared the morphology-density relation from the SDSS (0.01$<z<$0.054) with that of the MORPHS data ($z\sim$0.5). Two relations lie on top of each other, suggesting that the morphology-density relation was already established at $z\sim$0.5 as is in the present universe.

Using mean relative peculiar velocity measurements for pairs of galaxies, we estimate the cosmological density parameter $\Omega_m$ and the amplitude of density fluctuations $\sigma_8$. Our results suggest that our statistic is a robust and reproducible measure of the mean pairwise velocity and thereby the $\Omega_m$ parameter. We get $\Omega_m = 0.30^{+0.17}_{-0.07}$ and $\sigma_8 = 1.13^{+0.22}_{-0.23}$. These estimates do not depend on prior assumptions on the adiabaticity of the initial density fluctuations, the ionization history, or the values of other cosmological parameters.

The Sloan Digital Sky Survey (SDSS) first data release provides a database of 106000 unique galaxies in the main galaxy sample with measured spectra. A sample of star-forming (SF) galaxies are identified from among the 3079 of these having 1.4 GHz luminosities from FIRST, by using optical spectral diagnostics. Using 1.4 GHz luminosities as a reference star formation rate (SFR) estimator insensitive to obscuration effects, the SFRs derived from the measured SDSS Halpha, [OII] and u-band luminosities, as well as far-infrared luminosities from IRAS, are compared. It is established that straightforward corrections for obscuration and aperture effects reliably bring the SDSS emission line and photometric SFR estimates into agreement with those at 1.4 GHz, although considerable scatter (~60%) remains in the relations. It thus appears feasible to perform detailed investigations of star formation for large and varied samples of SF galaxies through the available spectroscopic and photometric measurements from the SDSS. We provide herein exact prescriptions for determining the SFR for SDSS galaxies. The expected strong correlation between [OII] and Halpha line fluxes for SF galaxies is seen, but with a median line flux ratio F_[OII]/F_Halpha=0.23, about a factor of two smaller than that found in the sample of Kennicutt (1992). This correlation, used in deriving the [OII] SFRs, is consistent with the luminosity-dependent relation found by Jansen et al. (2001). The median obscuration for the SDSS SF systems is found to be A_Halpha=1.2 mag, while for the radio detected sample the median obscuration is notably higher, 1.6 mag, and with a broader distribution.

We present central velocity dispersions and Mg$_2$ line indices for an all-sky sample of $\sim$ 1178 elliptical and S0 galaxies, of which 984 had no previous measures. This sample contains the largest set of homogeneous spectroscopic data for a uniform sample of elliptical galaxies in the nearby universe. These galaxies were observed as part of the ENEAR project, designed to study the peculiar motions and internal properties of the local early-type galaxies. Using 523 repeated observations of 317 galaxies obtained during different runs, the data are brought to a common zero-point. These multiple observations, taken during the many runs and different instrumental setups employed for this project, are used to derive statistical corrections to the data and are found to be relatively small, typically $\lsim$ 5% of the velocity dispersion and 0.01 mag in the \mg line-strength. Typical errors are about 8% in velocity dispersion and 0.01 mag in \mg, in good agreement with values published elsewhere.

In previous work on galaxy clusters, several authors reported a discovery of an unusual population of galaxies, which have spiral morphologies, but do not show any star formation activity. These galaxies are called ``passive spirals'', and have been interesting since it has been difficult to understand the existence of such galaxies. Using a volume limited sample (0.05<z<0.1 and Mr<-20.5; 25813 galaxies) of the Sloan Digital Sky Survey data, we have found 73 (0.28$\pm$0.03%) passive spiral galaxies and studied their environments. It is found that passive spiral galaxies live in local galaxy density 1-2 Mpc$^{-2}$ and 1-10 cluster-centric virial radius. Thus the origins of passive spiral galaxies are likely to be cluster related. These characteristic environments coincide with the previously reported environment where galaxy star formation rate suddenly declines and the so-called morphology-density relation turns. It is likely that the same physical mechanism is responsible for all of these observational results. The existence of passive spiral galaxies suggests that a physical mechanism that works calmly is preferred to dynamical origins such as major merger/interaction since such a mechanism can destroy spiral arm structures. Compared with observed cluster galaxy evolution such as the Butcher-Oemler effect and the morphological Butcher-Oemler effect, passive spiral galaxies are likely to be a galaxy population in transition between red, elliptical/S0 galaxies in low redshift clusters and blue, spiral galaxies more numerous in higher redshift clusters.

We report the discovery of a new gravitationally lensed quasar from the Sloan Digital Sky Survey, SDSS J090334.92+502819.2. This object was targeted for SDSS spectroscopy as a Luminous Red Galaxy (LRG), but manual examination of the spectrum showed the presence of a quasar at z= 3.6 in addition to a red galaxy at z=0.388, and the SDSS image showed a second possible quasar image nearby. Follow-up imaging and spectroscopy confirmed the lensing hypothesis. In images taken at the ARC 3.5-meter telescope, two quasars are separated by 2.8 arc-seconds; the lensing galaxy is clearly seen and is blended with one of the quasar images. Spectroscopy taken at the Keck II telescope shows that the quasars have identical redshifts of z=3.6 and both show the presence of the same broad absorption line-like troughs. We present simple lens models which account for the geometry and magnifications. The lens galaxy lies near two groups of galaxies and may be a part of them. The models suggest that the groups may contribute considerable shear and may have a strong effect on the lens configuration.

The distribution of early-type galaxy velocity dispersions, phi(sigma), is measured using a sample drawn from the SDSS database. Its shape differs significantly from that which one obtains by simply using the mean correlation between luminosity, L, and velocity dispersion, sigma, to transform the luminosity function into a velocity function: ignoring the scatter around the mean sigma-L relation is a bad approximation. An estimate of the contribution from late-type galaxies is also made, which suggests that phi(sigma) is dominated by early-type galaxies at velocities larger than ~ 200 km/s.

We present the observed fraction of galaxies with an Active Galactic Nucleus (AGN) as a function of environment in the Early Data Release of the Sloan Digital Sky Survey (SDSS). Using 4921 galaxies between 0.05 <= z <= 0.095, and brighter than M_r* = -20.0 (or M* +1.45), we find at least ~ 20% of these galaxies possess an unambiguous detection of an AGN, but this fraction could be as high as ~40% after we model the ambiguous emission line galaxies in our sample. We have studied the environmental dependence of galaxies using the the distance to the 10^th nearest neighbor. As expected, we observe that the fraction of star--forming galaxies decreases with density, while the fraction of passive galaxies increases with density. In contrast, the fraction of galaxies with an AGN remains constant from the cores of galaxy clusters to the rarefied field population. We conclude that the presence of an AGN is independent of the disk component of a galaxy. Our analyses are robust against measurement error, definition of an AGN, aperture bias, stellar absorption, survey geometry and signal--to--noise. Our observations are consistent with the hypothesis that a supermassive black hole resides in the bulge of all massive galaxies and ~40% of these black holes are seen as AGNs in our sample. A high fraction of local galaxies with an AGN suggests that either the mean lifetime of these AGNs is longer than previously thought (>10^8 years), or that the AGN burst more often than expected; ~40 times over the redshift range of our sample.

We study the variation of the dark matter mass fraction of elliptical galaxies as a function of their luminosity, stellar mass, and size using a sample of 29,469 elliptical galaxies culled from the Sloan Digital Sky Survey. We model ellipticals as a stellar Hernquist profile embedded in an adiabatically compressed dark matter halo. This model allows us to estimate a dynamical mass ($M_{dynm}$) at the half-light radius from the velocity dispersion of the spectra, and to compare these to the stellar mass estimates ($M_{*}$) from Kauffmann et al (2003). We find that $M_{*}/L$ is independent of luminosity, while $M_{dynm}/L$ increases with luminosity, implying that the dark matter fraction increases with luminosity. We also observe that at a fixed luminosity or stellar mass, the dark matter fraction increases with increasing galaxy size or, equivalently, increases with decreasing surface brightness: high surface brightness galaxies show almost no evidence for dark matter, while in low surface brightness galaxies, the dark matter exceeds the stellar mass at the half light radius. We relate this to the fundamental plane of elliptical galaxies, suggesting that the tilt of this plane from simple virial predictions is due to the dark matter in galaxies. We find that a simple model where galaxies are embedded in dark matter halos and have a star formation efficiency independent of their surface brightness explains these trends. We estimate the virial mass of ellipticals as being approximately 7-30 times their stellar mass, with the lower limit suggesting almost all of the gas within the virial radius is converted into stars.

We investigate the evolution of the fractions of late type cluster galaxies as a function of redshift, using one of the largest, most uniform cluster samples available. The sample consists of 514 clusters of galaxies in the range 0.02<z<0.3 from the Sloan Digital Sky Survey Cut & Enhance galaxy cluster catalog. This catalog was created using a single automated cluster finding algorithm on uniform data from a single telescope, with accurate CCD photometry, thus, minimizing selection biases. We use four independent methods to analyze the evolution of the late type galaxy fraction. Specifically, we select late type galaxies based on: restframe g-r color, u-r color, galaxy profile fitting and concentration index. The first criterion corresponds to the one used in the classical Butcher-Oemler analyses. The last three criteria are more sensitive to the morphological type of the galaxies. In all four cases, we find an increase in the fraction of late type galaxies with increasing redshift, significant at the 99.9% level. The results confirm that cluster galaxies do change colors with redshift (the Butcher-Oemler effect) and, in addition, they change their morphology to later-type toward higher redshift -- indicating a morphological equivalent of the Butcher-Oemler effect. We also find a tendency of richer clusters to have lower fractions of late type galaxies. The trend is consistent with a ram pressure stripping model, where richer clusters have more effective ram pressure due to their higher temperature.

We present here results on the composite luminosity functions of galaxies in the clusters of galaxies selected from the SDSS Cut and Enhance cluster catalog (CE; Goto et al. 2001). We construct the composite luminosity function in the five SDSS bands, u,g,r,i and z, using 204 CE clusters ranging from z=0.02 to z=0.25. We use photometric redshifts to construct composite luminosity functions. The robustness of the weighting scheme was tested using Monte Carlo simulation. We find the slope of composite LFs become flatter toward redder color band. Comparing with field LFs of SDSS (Blanton et al. 2001), cluster LFs have brighter characteristic magnitude and flatter slopes in g, r, i and z bands. These results are consistent with the hypothesis that the cluster LF has two distinct underlying populations i.e. the bright end of the LF is dominated by bright early types that follow a gaussian--like luminosity distribution, while the faint--end of the cluster LF is a steep power--law like function dominated by star--forming (bluer) galaxies. We also study the composite luminosity functions for early--type and late--type galaxies using profile fits, a concentration parameter and u-r color to classify galaxy morphology. The strong dependence of LF on galaxy morphology is found.The faint end slope of the LF is always flatter for early--type galaxies than late--type regardless of passband and methodology. These results are consistent with the hypothesis that the cluster region are dominated by bright elliptical galaxies.

We investigate the hypothesis that quasars formed together with the stellar populations of early-type galaxies. This hypothesis, in conjunction with the stellar ages of early-type galaxies from population synthesis models, the relation of black hole mass to bulge velocity dispersion, and the velocity dispersion distribution from the Sloan Digital Sky Survey, completely determines the cosmic accretion hystory of supermassive black holes and the redshift evolution of the characteristic luminosity. On the other hand the precise shape of the luminosity function of quasars depends on the light curve of quasars and, in the optical, but not so much in X-rays, on the covering factor of the dust surrounding the AGN. We find a plausible set of assumptions for which the coeval formation of supermassive black holes and elliptical galaxies is in good agreement with the observed B-band and X-ray luminosity functions of quasars.

The primary observational goals of the Sloan Digital Sky Survey are to obtain CCD imaging of 10,000 deg^2 of the north Galactic cap in five passbands, with a limiting magnitude in the r-band of 22.5, to obtain spectroscopic redshifts of 10^6 galaxies and 10^5 quasars, and to obtain similar data for three ~ 200 deg^2 stripes in the south Galactic cap, with repeated imaging to allow co-addition and variability studies in at least one of these stripes. The resulting photometric and spectroscopic galaxy datasets allow one to map the large scale structure traced by optical galaxies over a wide range of scales to unprecedented precision. Results relevant to the large scale structure of our Universe include: a flat model with a cosmological constant Omega_Lambda=0.7 provides a good description of the data; the galaxy-galaxy correlation function shows departures from a power law which are statistically significant; and galaxy clustering is a strong function of galaxy type.

We present R-band CCD photometry for 1332 early-type galaxies, observed as part of the ENEAR survey of peculiar motions using early-type galaxies in the nearby Universe. Circular apertures are used to trace the surface brightness profiles, which are then fit by a two-component bulge-disk model. From the fits we obtain the structural parameters required to estimate galaxy distances using the D_n-sigma and Fundamental Plane relations. We find that about 12% of the galaxies are well represented by a pure r^{1/4} law while 87% are best fit by a two component model. There are 356 repeated observations of 257 galaxies obtained during different runs that are used to derive statistical corrections and bring the data to a common system. We also use these repeated observations to estimate our internal errors. The accuracy of our measurements are tested by the comparison of 354 galaxies in common with other authors. Typical errors in our measurements are 0.011 dex for log{D_n}, 0.064 dex for log{r_e}, 0.086 mag arcsec^{-2} for <{mu_e}> and 0.09 for m_{R_C}, comparable to those estimated by other authors. The photometric data reported here represent one of the largest high-quality and uniform all-sky samples currently available for early-type galaxies in the nearby universe, especially suitable for peculiar motion studies.

A sample of nearly 9000 early-type galaxies, in the redshift range 0.01 < z < 0.3, was selected from the Sloan Digital Sky Survey using morphological and spectral criteria. This paper describes how the sample was selected, presents examples of images and seeing corrected fits to the observed surface brightness profiles, describes our method for estimating K-corrections, and shows that the SDSS spectra are of sufficiently high quality to measure velocity dispersions accurately. It also provides catalogs of the measured photometric and spectroscopic parameters. In related papers, these data are used to study how early-type galaxy observables, including luminosity, effective radius, surface brightness, color, and velocity dispersion, are correlated with one another.

A magnitude limited sample of nearly 9000 early-type galaxies, in the redshift range 0.01 < z < 0.3, was selected from the Sloan Digital Sky Survey using morphological and spectral criteria. The sample was used to study how early-type galaxy observables, including luminosity L, effective radius R_o, surface brightness I_o, color, and velocity dispersion sigma, are correlated with one another. Measurement biases are understood with mock catalogs which reproduce all of the observed scaling relations and their dependences on fitting technique. At any given redshift, the intrinsic distribution of luminosities, sizes and velocity dispersions in our sample are all approximately Gaussian. A maximum likelihood analysis shows that sigma ~ L^{0.25\pm 0.012}, R_o ~ L^{0.63\pm 0.025}, and R_o ~ I^{-0.75\pm 0.02} in the r* band. In addition, the mass-to-light ratio within the effective radius scales as M_o/L ~ L^{0.14\pm 0.02} or M_o/L ~ M_o^{0.22\pm 0.05}, and galaxies with larger effective masses have smaller effective densities: Delta_o ~ M_o^{-0.52\pm 0.03}. These relations are approximately the same in the g*, i* and z* bands. Relative to the population at the median redshift in the sample, galaxies at lower and higher redshifts have evolved only little, with more evolution in the bluer bands. The luminosity function is consistent with weak passive luminosity evolution and a formation time of about 9 Gyrs ago.

A magnitude limited sample of nearly 9000 early-type galaxies, in the redshift range 0.01 < z < 0.3, was selected from the Sloan Digital Sky Survey using morphological and spectral criteria. The Fundamental Plane relation in this sample is R_o ~ sigma^{1.49\pm 0.05} I_o^{-0.75\pm 0.01} in the r* band. It is approximately the same in the g*, i* and z* bands. Relative to the population at the median redshift in the sample, galaxies at lower and higher redshifts have evolved only little. If the Fundamental Plane is used to quantify this evolution then the apparent magnitude limit can masquerade as evolution; once this selection effect has been accounted for, the evolution is consistent with that of a passively evolving population which formed the bulk of its stars about 9 Gyrs ago. One of the principal advangtages of the SDSS sample over previous samples is that the galaxies in it lie in environments ranging from isolation in the field to the dense cores of clusters. The Fundamental Plane shows that galaxies in dense regions are slightly different from galaxies in less dense regions.

The colors and chemical abundances of early-type galaxies at redshifts z<0.3 are studied using a sample of nearly 9000 galaxies, selected from the Sloan Digital Sky Survey using morphological and spectral criteria. In this sample, redder galaxies have larger velocity dispersions. Color also correlates with magnitude and size, but these correlations are entirely due to the L-sigma and R_o-sigma relations: the primary correlation is color-sigma. The red light in early-type galaxies is, on average, slightly more centrally concentrated than the blue. Because of these color gradients, the strength of the color--magnitude relation depends on whether or not the colors are defined using a fixed metric aperture. Chemical evolution and star formation histories of early-type galaxies are investigated using co-added spectra of similar objects. At fixed sigma, the population at z ~ 0.2 had weaker Mg_2 and stronger H_beta absorption compared to the population at z ~ 0. It was also bluer. Comparison of these colors and line-strengths, and their evolution, with single-burst stellar population models suggests a formation time of 9 Gyrs ago, consistent with a Fundamental Plane analysis of this sample.

[Abridged] We present here a new and homogeneous sample of 3340 galaxies selected from the Sloan Digital Sky Survey (SDSS) based solely on the observed strength of their Hdelta absorption line. These galaxies are commonly known as ``post-starburst'' or ``E+A'' galaxies, and the study of these galaxies has been severely hampered by the lack of a large, statistical sample of such galaxies. In this paper, we rectify this problem by selecting a sample of galaxies which possess an absorption Hdelta equivalent width of EW(Hdelta_max) - Delta EW(Hdelta_max) > 4A from 106682 galaxies in the SDSS. We have performed extensive tests on our catalog including comparing different methodologies of measuring the Hdelta absorption and studying the effects of stellar absorption, dust extinction, emission-filling and measurement error. The measured abundance of our Hdelta-selected (HDS) galaxies is 2.6 +/- 0.1% of all galaxies within a volume-limited sample of 0.05<z<0.1 and M(r*)<-20.5, which is consistent with previous studies of such galaxies in the literature. We find that only 25 of our HDS galaxies in this volume-limited sample (3.5+/-0.7%) show no evidence for OII and Halpha emission, thus indicating that true E+A (or k+a) galaxies are extremely rare objects at low redshift, i.e., only 0.09+/-0.02% of all galaxies in this volume-limited sample are true E+A galaxies. In contrast, 89+/-5% of our HDS galaxies in the volume-limited sample have significant detections of the OII and Halpha emission lines. We find 27 galaxies in our volume-limited HDS sample that possess no detectable OII emission, but do however possess detectable Halpha emission. These galaxies may be dusty star-forming galaxies. We provide the community with this new catalog of Hdelta-selected galaxies to aid in the understanding of these galaxies.

We combine Sloan Digital Sky Survey spectra of 22,000 luminous, red, bulge-dominated galaxies to get high S/N average spectra in the rest-frame optical and ultraviolet (2600A to 7000A). The average spectra of these massive, quiescent galaxies are early-type with weak emission lines and with absorption lines indicating an apparent excess of alpha elements over solar abundance ratios. We make average spectra of subsamples selected by luminosity, environment and redshift. The average spectra are remarkable in their similarity. What variations do exist in the average spectra as a function of luminosity and environment are found to form a nearly one-parameter family in spectrum space. We present a high signal-to-noise ratio spectrum of the variation. We measure the properties of the variation with a modified version of the Lick index system and compare to model spectra from stellar population syntheses. The variation may be a combination of age and chemical abundance differences, but the conservative conclusion is that the quality of the data considerably exceeds the current state of the models.

We develop a new method to constrain the star formation histories, dust attenuation and stellar masses of galaxies. It is based on two stellar absorption line indices, the 4000 Angstrom break strength and the Balmer absorption line index Hdelta_A. Together, these indices allow us to constrain the mean stellar ages of galaxies and the fractional stellar mass formed in bursts over the past few Gyr. A comparison with broad band photometry then yields estimates of dust attenuation and of stellar mass. We generate a large library of Monte Carlo realizations of different star formation histories, including starbursts of varying strength and a range of metallicities. We use this library to generate median likelihood estimates of burst mass fractions, dust attenuation strengths, stellar masses and stellar mass-to-light ratios for a sample of 122,208 galaxies drawn from the Sloan Digital Sky Survey. The typical 95% confidence range in our estimated stellar masses is +-40%. We study how the stellar mass-to-light ratios of galaxies vary as a function of absolute magnitude, concentration index and photometric pass-band and how dust attenuation varies as a function of absolute magnitude and 4000 Angstrom break strength. We also calculate how the total stellar mass of the present Universe is distributed over galaxies as a function of their mass, size, concentration, colour, burst mass fraction and surface mass density. We find that most of the stellar mass in the local Universe resides in galaxies that have stellar masses \~5\times 10^10 M_sol, half light radii ~3 kpc, and half- light surface mass densities ~10^9 M_sol/kpc^2. The distribution of D(4000) is strongly bimodal, showing a clear division between galaxies dominated by old stellar populations and galaxies with more recent star formation.

(Abridged) We present in this paper a detailed analysis of the effect of environment on the star-formation activity of galaxies within the EDR of the SDSS. We have used the Halpha emission line to derive the star-formation rate (SFR) for each galaxy within a volume-limited sample of 8598 galaxies with 0.05 < z < 0.095 and M(r)<= -20.45. We find that the SFR of galaxies is strongly correlated with the local (projected) galaxy density and thus we present here the density-SFR relation that is analogous to the density-morphology relation. The effect of density on the SFR of galaxies is seen in three ways. First, the overall distribution of SFRs is shifted to lower values in dense environments compared with the field population. Second, the effect is most noticeable for the strongly star-forming galaxies in the 75th percentile of the SFR distribution. Third, there is a ``break'' (or characteristic density) in the density-SFR relation at a local galaxy density of 1h-2 Mpc-2. To understand this break further, we have studied the SFR of galaxies as a function of clustercentric radius from 17 clusters and groups objectively selected from the SDSS EDR data. The distribution of SFRs of cluster galaxies begins to change, compared with the field population, at a clustercentric radius of 3-4 virial radii, which is consistent with the characteristic break in density that we observe in the density-SFR relation. Our tests suggest that the density-morphology relation alone is unlikely to explain the density-SFR relation we observe. Taken all together, these works demonstrate that the decrease in SFR of galaxies in dense environments is a universal phenomenon over a wide range in density (from 0.08 to 10h-2 Mpc-2) and redshift (out to z = 0.5).

The effective optical depth in the Ly-alpha forest region of 1061 low-resolution QSO spectra drawn from the SDSS database decreases with decreasing redshift over the range 2.5 < z < 4. Although the evolution is relatively smooth, tau_eff ~ (1+z)^{3.8 pm 0.2}, at z ~ 3.2 the effective optical depth decreases suddenly, by about ten percent with respect to this smoother evolution. It climbs back to the original smooth scaling again by z ~ 2.9. We describe two techniques, one of which is new, for quantifying this evolution which give consistent results. A variety of tests show that the feature is not likely to be a consequence of how the QSO sample was selected, nor the result of flux calibration or other systematic effects. Other authors have argued that, at this same epoch, the temperature of the IGM also shows a departure from an otherwise smooth decrease with time. These features in the evolution of the temperature and the optical depth are signatures of the reionization of He II.

The expansion of the Universe leads to a rapid drop in the mean hydrogen Ly-alpha opacity tau of the intergalactic medium (IGM), tau (1+z)^{3.8}, between redshifts 4 and 3. Measurements of the temperature evolution of the IGM and of the He II opacity both suggest that He II reionizes in this redshift range. We use hydrodynamical simulations to show that the temperature increase associated with He II reionization leads to a relatively sudden decrease in tau around the reionization epoch of ~ 10 per cent. We find clear evidence for such a feature in the evolution of tau determined from a sample of ~ 1100 quasars obtained from the SDSS. He II reionization starts at redshift ~ 3.4, and lasts for (Delta z) ~ 0.4. The increase in the IGM temperature also explains the widths of hydrogen absorption lines as measured in high-resolution spectra.

This paper presents data on the ENEARc subsample of the larger ENEAR survey of the nearby early-type galaxies. The ENEARc galaxies belong to clusters and were specifically chosen to be used for the construction of a D_n-sigma template. The ENEARc sample includes new measurements of spectroscopic and photometric parameters (redshift, velocity dispersion, line index Mg_2, and the angular diameter d_n) as well as data from the literature. New spectroscopic data are given for 229 cluster early-type galaxies and new photometry is presented for 348 objects. Repeat and overlap observations with external data sets are used to construct a final merged catalog consisting of 640 early-type galaxies in 28 clusters. Objective criteria, based on catalogs of groups of galaxies derived from complete redshift surveys of the nearby universe, are used to assign galaxies to clusters. In a companion paper these data are used to construct the template D_n-sigma distance relation for early-type galaxies which has been used to estimate galaxy distances and derive peculiar velocities for the ENEAR all-sky sample.

In this paper R-band photometric and velocity dispersion measurements for a sample of 452 elliptical and S0 galaxies in 28 clusters are used to construct a template D_n-sigma relation. This template relation is constructed by combining the data from the 28 clusters, under the assumption that galaxies in different clusters have similar properties. The photometric and spectroscopic data used consist of new as well as published measurements, converted to a common system, as presented in a accompanying paper. The resulting direct relation, corrected for incompleteness bias, is log{D_n} =1.203 log{sigma} + 1.406; the zero-point has been defined by requiring distant clusters to be at rest relative to the CMB. This zero-point is consistent with the value obtained by using the distance to Virgo as determined by the Cepheid period-luminosity relation. This new D_n-sigma relation leads to a peculiar velocity of -72 (\pm 189) km/s for the Coma cluster. The scatter in the distance relation corresponds to a distance error of about 20%, comparable to the values obtained for the Fundamental Plane relation. Correlations between the scatter and residuals of the D_n-sigma relation with other parameters that characterize the cluster and/or the galaxy stellar population are also analyzed. The direct and inverse relations presented here have been used in recent studies of the peculiar velocity field mapped by the ENEAR all-sky sample.

A sample of nearly 9000 early-type galaxies, in the redshift range 0.01<z<0.3, was selected from the Sloan Digital Sky Survey using morphological and spectral criteria. The sample was used to study how early-type galaxy observables, including luminosity L, effective radius R, surface brightness I, color, and velocity dispersion V, are correlated with one another. Measurement biases are understood with mock catalogs which reproduce all of the observed scaling relations. At any given redshift, the intrinsic distribution of luminosities, sizes and velocity dispersions in our sample are all approximately Gaussian. In the r* band L ~ V^3.91, L ~ R^1.58, R ~ I^(-0.75), and the Fundamental Plane relation is R ~ V^(1.49) I^(-0.75). These relations are approximately the same in the g*, i* and z* bands. At fixed luminosity, the mass-to-light ratio scales as M/L ~ L^0.14. The g*-r* color scales as V^0.25. Color also correlates with magnitude and size, but these correlations are entirely due to the L-V and R-V relations. Chemical evolution and star formation histories are investigated using co-added spectra of similar objects in our sample. Chemical abundances correlate primarily with velocity dispersion. At fixed V, the higher redshift population is bluer, is weaker in Mg2, and is stronger in Hbeta than the population nearby. In addition, the population at higher redshifts is slightly more luminous. These differences are consistent with that of a passively evolving population which formed the bulk of its stars about 9 Gyrs ago. The Fundamental Plane suggests that galaxies in dense regions are slightly different from those in less dense regions, but the co-added spectra and color--magnitude relations show no statistically significant dependence on environment.

We analyze new measurements of the Mg_2 central line strength index and velocity dispersion (sigma) for the galaxies of the ENEAR survey. The observations are now complete (da Costa et al. 2000) and the sample contains 1223 early-type galaxies. We also analyze the line strength indices for a sample of 95 spiral bulges (from Sa to Sbc). For the early-type galaxies we find: i) that the Mg_2-sigma relation for Es and S0s are nearly the same, with both populations showing comparable scatter, and ii) a marginal difference in the slope of the Mg_2-sigma relation for cluster and field early-type galaxies. However, we suggest that before interpreting such a difference in the framework of a mass-metallicity relation, it is important to take into account the effects of rotation in the Mg_2-sigma relation. Our preliminary results indicate that once the rotation effects are minimized by choosing a sample containing only slow rotators, the Mg_2-sigma relation is similar both for isolated and clustered galaxies. More data on rotational velocities of early-type galaxies are needed to confirm this result. For spiral bulges, we find that their locus in the Mg_2-sigma plane lies always below the one occupied by early-type galaxies.

We describe the target selection and resulting properties of a spectroscopic sample of luminous, red galaxies (LRG) from the imaging data of the Sloan Digital Sky Survey (SDSS). These galaxies are selected on the basis of color and magnitude to yield a sample of luminous, intrinsically red galaxies that extends fainter and further than the main flux-limited portion of the SDSS galaxy spectroscopic sample. The sample is designed to impose a passively-evolving luminosity and rest-frame color cut to a redshift of 0.38. Additional, yet more luminous, red galaxies are included to a redshift of 0.5. Approximately 12 of these galaxies per square degree are targeted for spectroscopy, so the sample will number over 100,000 with the full survey. SDSS commissioning data indicate that the algorithm efficiently selects luminous (M_g=-21.4), red galaxies, that the spectroscopic success rate is very high, and that the resulting set of galaxies is approximately volume-limited out to z=0.38. When the SDSS is complete, the LRG spectroscopic sample will fill over 1h^-3 Gpc^3 with an approximately homogeneous population of galaxies and will therefore be well suited to studies of large-scale structure and clusters out to z=0.5.

We use the recently completed redshift-distance survey of nearby early-type galaxies (ENEAR) to measure the dipole component of the peculiar velocity field to a depth of cz ~ 6000 km/s. The sample consists of 1145 galaxies brighter than m_B=14.5 and cz < 7000 km/s, uniformly distributed over the whole sky, and 129 fainter cluster galaxies within the same volume. Most of the Dn-sigma distances were obtained from new spectroscopic and photometric observations conducted by this project, ensuring the homogeneity of the data over the whole sky. These 1274 galaxies are objectively assigned to 696 objects -- 282 groups/clusters and 414 isolated galaxies. We find that within a volume of radius ~ 6000 km/s, the best-fitting bulk flow has an amplitude of |vbulk| =220 +/- 42 km/s in the CMB restframe, pointing towards l=304 +/- 16 degrees, b=25 +/- 11 degrees. The error in the amplitude includes statistical, sampling and possible systematic errors. This solution is in excellent agreement with that obtained by the SFI Tully-Fisher survey. Our results suggest that most of the motion of the Local Group is due to fluctuations within 6000 km/s, in contrast to recent claims of large amplitude bulk motions on larger scales.

We present a comparison between the peculiar velocity field measured from the ENEAR all-sky $D_n-\sigma$ catalog and that derived from the galaxy distribution of the IRAS PSCz redshift survey. The analysis is based on a modal expansion of these data in redshift space by means of spherical harmonics and Bessel functions. The effective smoothing scale of the expansion is almost linear with redshift reaching 1500km/s at 3000km/s. The general flow patterns in the filtered ENEAR and PSCz velocity fields agree well within 6000km/s, assuming a linear biasing relation between the mass and the PSCz galaxies. The comparison allows us to determine the parameter $\beta=\Omega^{0.6}/b$, where $\Omega$ is the cosmological density parameter and $b$ is the linear biasing factor. A likelihood analysis of the ENEAR and PSCz modes yields $\beta=0.5 +- 0.1$, in good agreement with values obtained from Tully-Fisher surveys.

We estimate the mass density fluctuations power spectrum (PS) on large scales by applying a maximum likelihood technique to the peculiarvelocity data of the recently completed redshift-distance survey of early-type galaxies (ENEAR). The general results are in agreement with the high amplitude power spectra found from similar analysis of other independent all-sky catalogs of peculiar velocity data such as MARK III and SFI. For Lambda & Open CDM COBE normalized PS models, the best-fit parameters are confined by a contour approximately defined by Omega h^{1.3}=0.377+-0.08 and Omega h^{0.88}=0.517+-0.083, respectively. Gamma-shape models, free of COBE normalization, resultsin the weak constraint of $\Gamma \geq 0.17$ and in the rather stringent constraint of sigma_8 Omega^{0.6}=1.0+-0.25. All quoted uncertainties refer to 3-sigma confidence-level. The calculated PS is used as a prior for Wiener reconstruction of the density field at different resolutions and the three-dimensional velocity field within a volume of radius ~80 Mpc/h. All major structures in the nearby universe are recovered and are well matched to those predicted from all-sky redshift surveys.

We present an analysis of the ENEAR sample of peculiar velocities of elliptical galaxies, obtained with D_n-\sigma distances. We use the velocity correlation function to analyze the statistics of the field-object's velocities, while the analysis of the cluster data is based on the estimate of their rms peculiar velocity, Vrms. The statistics of the model velocity field is parameterized by the amplitude, \eta_8=\sigma_8 \Omega_m^{0.6}, and by the shape parameter, \Gamma. From the velocity correlation statistics we obtain \eta_8=0.51{-0.09}{+0.24} for \Gamma=0.25 at the 2\sigma level. Even though less constraining, a consistent result is obtained by comparing the measured Vrms of clusters to linear theory predictions. For \Gamma=0.25 we find \eta_8=0.63{-0.19}{+0.22}$ at 1\sigma. Overall, our results point toward a statistical concordance of the cosmic flows traced by spirals and early-type galaxies, with galaxy distances estimated using TF and D_n-\sigma distance indicators, respectively.

This is the first in a series of papers describing the recently completed all-sky redshift-distance survey of nearby early-type galaxies (ENEAR) carried out for peculiar velocity analysis. The sample is divided into two parts and consists of 1607 elliptical and lenticular galaxies with cz < 7000 km/s and with blue magnitudes brighter than m_B=14.5 (ENEARm), and of galaxies in clusters (ENEARc). Galaxy distances based on the Dn-sigma and Fundamental Plane (FP) relations are now available for 1359 and 1107 ENEARm galaxies, respectively, with roughly 80% based on new data gathered by our group. The Dn-sigma and FP template distance relations are derived by combining 569 and 431 galaxies in 28 clusters, respectively, of which about 60% are based on our new measurements. The ENEARm redshift-distance survey extends the earlier work of the 7S and the recent Tully-Fisher surveys sampling a comparable volume. In subsequent papers of this series we intend to use the ENEAR sample by itself or in combination with the SFI Tully-Fisher survey to analyze the properties of the local peculiar velocity field and to test how sensitive the results are to different sampling and to the distance indicators. We also anticipate that the homogeneous database assembled will be used for a variety of other applications and serve as a benchmark for similar studies at high-redshift.

The ENEAR project is an all-sky survey of nearby early-type galaxies. About 2200 new spectra and R-band images of over 1500 galaxies together with published data were assembled into a large catalog containing over 2000 objects with cz, sigma, photometric data, and line strengths on a uniform system. From this extensive database a magnitude-limited sample has been drawn comprising ~ 1400 galaxies brighter than m_B = 14.5, cz < 7000 km/s, and type T -2 or less with measured distances (ENEARf) and about 500 early-type galaxies in 28 clusters/groups (ENEARc) to derive an internally consistent D_n-sigma relation to estimate galaxy distances. In this contribution we discuss some general properties of ENEAR and briefly describe some preliminary results.

Using new observations for a sample of 931 early-type galaxies we investigate whether the \mg2--\so relation shows any dependence on the local environment. The galaxies have been assigned to three different environments depending on the local overdensity: clusters, groups, and field, having used our completeredshift database to guide the assignment of galaxies. It is found that cluster, group and field early-type galaxies follow almost identical \mg2--\so\ relations, with the largest \mg2 zero-point difference (clusters minus field) being only $0.007\pm 0.002$ mag. No correlation of the residuals is found with the morphological type or the bulge to disk ratio. Using stellar population models in a differential fashion, this small zero-point difference implies a luminosity-weighted age difference of only $\sim 1$ Gyr between the corresponding stellar populations, with field galaxies being younger. The mass-weighted age difference could be significantly smaller, if minor events of late star formation took place preferentially in field galaxies. We combine these results with the existing evidence for the bulk of stars in cluster early-type galaxies having formed at very high redshift, and conclude that the bulk of stars in galactic spheroids had to form at high redshifts ($z\gsim 3$), no matter whether such spheroids now reside in low or high density regions. The cosmological implications of these findings are briefly discussed.