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

We investigate the Fe, C, N, O, Mg, Al, Si, K, Ca, Ce and Nd abundances of 2283 red giant stars in 31 globular clusters from high-resolution spectra observed in both the northern and southern hemisphere by the SDSS-IV APOGEE-2 survey. This unprecedented homogeneous dataset, largest to date, allows us to discuss the intrinsic Fe spread, the shape and statistics of Al-Mg and N-C anticorrelations as a function of cluster mass, luminosity, age and metallicity for all 31 clusters. We find that the Fe spread does not depend on these parameters within our uncertainties including cluster metallicity, contradicting earlier observations. We do not confirm the metallicity variations previously observed in M22 and NGC 1851. Some clusters show a bimodal Al distribution, while others exhibit a continuous distribution as has been previously reported in the literature. We confirm more than 2 populations in $\omega$ Cen and NGC 6752, and find new ones in M79. We discuss the scatter of Al by implementing a correction to the standard chemical evolution of Al in the Milky Way. After correction, its dependence on cluster mass is increased suggesting that the extent of Al enrichment as a function of mass was suppressed before the correction. We observe a turnover in the Mg-Al anticorrelation at very low Mg in $\omega$ Cen, similar to the pattern previously reported in M15 and M92. $\omega$ Cen may also have a weak K-Mg anticorrelation, and if confirmed, it would be only the third cluster known to show such a pattern.

We report the first APOGEE metallicities and alpha-element abundances measured for 3600 red giant stars spanning a large radial range of both the Large (LMC) and Small Magellanic Clouds (SMC), the largest Milky Way dwarf galaxies. Our sample is an order of magnitude larger than that of previous studies, and extends to much larger radial distances. These are the first results presented that make use of the newly installed Southern APOGEE instrument on the du Pont telescope at Las Campanas Observatory. Our unbiased sample of the LMC spans a large range in metallicity, from [Fe/H]=-0.2 to very metal-poor stars with [Fe/H]=-2.5, the most metal-poor Magellanic Clouds (MCs) stars detected to date. The LMC [alpha/Fe]-[Fe/H] distribution is very flat over a large metallicity range, but rises by ~0.1 dex at -1.0<[Fe/H]<-0.5. We interpret this as a sign of the known recent increase in MC star-formation activity, and are able to reproduce the pattern with a chemical evolution model that includes a recent "starburst". At the metal-poor end, we capture the increase of [alpha/Fe] with decreasing [Fe/H], and constrain the "alpha-knee" to [Fe/H]<-2.2 in both MCs, implying a low star-formation efficiency of ~0.01 Gyr^-1. The MC knees are more metal poor than those of less massive Milky Way (MW) dwarf galaxies such as Fornax, Sculptor, or Sagittarius. One possible interpretation is that the MCs formed in a lower-density environment than the MW, a hypothesis that is consistent with the paradigm that the MCs fell into the MW's gravitational potential only recently.

We study the spatially resolved properties of 333 elliptical galaxies with the MaNGA/SDSS-IV survey. The aim is to understand the fundamental processes of formation and quenching of elliptical galaxies. We used the DESI Legacy Imaging Surveys for accurate morphological classification. Based on integrated spectroscopic properties and colors, we classified galaxies into classical "red and dead", recently quenched, and blue star-forming ellipticals (CLEs, RQEs, and BSFs corresponding to 75%, 10%, and 4% of the sample, respectively). We inferred their stellar age and stellar metallicity gradients out to 1.5 effective radius, and reconstructed their global and radial histories of mass growth and star formation. We find the mass- and light-weighted age gradients of CLEs are nearly flat or mildly negative, with small differences between both ages. The respective metallicity gradients are negative, being flatter as less massive are the CLEs. The more massive CLEs assembled earlier and quenched faster than the less massive ones. The CLEs show a weak inside-out growth and a clear inside-out quenching. At masses < $10^{11}$ $M_{\bigodot}$, the age and $Z$ gradients of the RQEs and BSFs are flatter than those of the CLEs but with larger scatters. They show very weak inside-out growth and quenching, being the quenching slow and even not completed at $z\sim$ 0 for the BSFs. Instead, the massive RQEs show an outside-in quenching, and positive gradients in the light-weighted age and stellar metallicities. Our results are consistent with a scenario where the inner parts of CLEs formed by an early and coeval dissipative collapse with a consequent burst of star formation and further quenching, whereas the outer parts continued their assembly likely by dry mergers. We also discuss some evolutionary scenarios for the RQE and BSF galaxies that would agree with their generic results presented here.

We use the improved IllustrisTNG300 magneto-hydrodynamical cosmological simulation to revisit the effect that secondary halo bias has on the clustering of the central galaxy population. With a side length of 205 $h^{-1}$Mpc and significant improvements on the sub-grid model with respect to the previous Illustris boxes, IllustrisTNG300 allows us to explore the dependencies of galaxy clustering over a large cosmological volume and wide halo-mass range. We show, at high statistical significance, that the halo assembly bias signal (i.e., the secondary dependence of halo bias on halo formation redshift) manifests itself on the clustering of the central galaxy population when this is split by stellar mass, colour, specific star formation rate, and surface density. A significant detection is also obtained for galaxy size: at fixed halo mass, larger central galaxies are more tightly clustered than smaller central galaxies in haloes of mass M$_{\rm vir} \lesssim 10^{12.5}$ $h^{-1}$M$_{\odot}$. This effect, however, seems to be uncorrelated with halo formation time, unlike the rest of the secondary dependencies analysed. We also explore the transmission of the halo spin bias signal, i.e., the secondary dependence of halo bias on halo spin. Although galaxy spin retains little information about the total spin of the halo, the correlation is enough to produce a significant galaxy spin bias signal. We discuss possible ways to probe the spin bias effects with observations.

Using 3D positions and kinematics of stars relative to the Sagittarius (Sgr) orbital plane and angular momentum, we identify 166 Sgr stream members observed by the Apache Point Observatory Galactic Evolution Experiment (APOGEE) that also have Gaia DR2 astrometry. This sample of 63/103 stars in the Sgr trailing/leading arm are combined with an APOGEE sample of 710 members of the Sgr dwarf spheroidal core (385 of them newly presented here) to establish differences of 0.6 dex in median metallicity and 0.1 dex in [$\alpha$/Fe] between our Sgr core and dynamically older stream samples. Mild chemical gradients are found internally along each arm, but these steepen when anchored by core stars. With a model of Sgr tidal disruption providing estimated dynamical ages (i.e., stripping times) for each stream star, we find a mean metallicity gradient of 0.12 +/- 0.03 dex/Gyr for stars stripped from Sgr over time. For the first time, an [$\alpha$/Fe] gradient is also measured within the stream, at 0.02 +/- 0.01 dex/Gyr using magnesium abundances and 0.04 +/- 0.01 dex/Gyr using silicon, which imply that the Sgr progenitor had significant radial abundance gradients. We discuss the magnitude of those inferred gradients and their implication for the nature of the Sgr progenitor within the context of the current family of Milky Way satellite galaxies, and suggest that more sophisticated Sgr models are needed to properly interpret the growing chemodynamical detail we have on the Sgr system.

We use a sample of ~3000 galaxies from the MaNGA MPL-7 internal data release to study the alpha abundance distribution within low-redshift galaxies. We use the Lick index ratio Mgb/<Fe> as an alpha abundance indicator to study relationships between the alpha abundance distribution and galaxy properties such as effective stellar velocity dispersion within 0.3 effective radii (sigma_*), galaxy environment, and dark matter halo formation time (z_f). We find that (1) all galaxies show a tight correlation between Mgb/<Fe> and sigma_*; (2) `old' (H_beta < 3) low-sigma_* galaxies in high local density environment and inner regions within galaxy groups are enhanced in Mgb/<Fe>, while `young' (H_beta>3) galaxies and high-mass galaxies show no or less environmental dependence; (3) `old' galaxies with high-z_f show enhanced Mgb/<Fe> over low- and medium-z_f; (4) Mgb/<Fe> gradients are close to zero and show dependence on sigma_* but no obvious dependence on the environment or z_f. Our study indicates that stellar velocity dispersion or galaxy mass is the main parameter driving the Mgb/<Fe> enhancement, although environments appear to have modest effects, particularly for low- and medium-mass galaxies.

We examine the evolution of assembly bias using a semi-analytical model of galaxy formation implemented in the Millennium-WMAP7 N-body simulation. We consider fixed number density galaxy samples ranked by stellar mass or star formation rate. We investigate how the clustering of haloes and their galaxy content depend on halo formation time and concentration, and how these relationships evolve with redshift. At $z=0$ the dependences of halo clustering on halo concentration and formation time are similar. However, at higher redshift, halo assembly bias weakens for haloes selected by age, and reverses and increases for haloes selected by concentration. The variation of the halo occupation with concentration and formation time is also similar at $z=0$ and changes at higher redshifts. In this case, the occupancy variation with halo age stays mostly constant with redshift but decreases for concentration. Finally, we look at the evolution of assembly bias reflected in the galaxy distribution by examining the galaxy correlation functions relative to those of shuffled galaxy samples which remove the occupancy variation. This correlation functions ratio monotonically decreases with larger redshift and for lower number density samples, going below unity in some cases, leading to reduced galaxy clustering. While the halo occupation functions themselves vary, the assembly bias trends are similar whether selecting galaxies by stellar mass or star formation rate. Our results provide further insight into the origin and evolution of assembly bias. Our extensive occupation function measurements and fits are publicly available and can be used to create realistic mock catalogues.

The large Integral Field Spectroscopy (IFS) surveys have allowed the classification of ionizing sources of emission lines on sub-kpc scales. In this work, we define two non-parametric parameters, quiescence (f$_{q}$) and its concentration (c$_{q}$), to quantify the strength and the spatial distribution of the quenched areas, respectively, traced by the LI(N)ER regions with low EW(H$\alpha$). With these two measurements, we classify MaNGA galaxies into inside-out and outside-in quenching types according to their locations on the f$_{q}$ vs. c$_{q}$ plane and we measure the fraction of inside-out (outside-in) quenching galaxies as a function of halo mass. We find that the fraction of galaxies showing inside-out quenching increases with halo mass, irrespective of stellar mass or galaxy type (satellites vs. centrals). In addition, high stellar mass galaxies exhibit a greater fraction of inside-out quenching compared to low stellar mass ones in all environments. In contrast, the fraction of outside-in quenching does not depend on halo mass. Our results suggest that morphological quenching may be responsible for the inside-out quenching seen in all environments. On the other hand, the flat dependence of the outside-in quenching on halo mass could be a mixed result of ram-pressure stripping and galaxy mergers. Nevertheless, at a given environment and stellar mass, the fraction of inside-out quenching is systematically greater than that of outside-in quenching, suggesting that inside-out quenching is the dominant quenching mode in all environments.

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 investigate the relationship of kinematically misaligned galaxies with their large-scale environment, in the context of halo assembly bias. According to numerical simulations, halo age at fixed halo mass is intrinsically linked to the large-scale tidal environment created by the cosmic web. We investigate the relationship between distances to various cosmic web features and present-time gas accretion rate. We select a sub-sample of ~900 central galaxies from the MaNGA survey with defined global position angles (PA; angle at which velocity change is greatest) for their stellar and H$\alpha$ gas components up to a minimum of 1.5 effective radii ($R_e$). We split the sample by misalignment between the gas and stars as defined by the difference in their PA. For each central galaxy we find its distance to nodes and filaments within the cosmic web, and estimate the host halo's age using the central stellar mass to total halo mass ratio $M_{*}/M_{h}$. We also construct halo occupation distributions using a background subtraction technique for galaxy groups split using the central galaxy's kinematic misalignment. We find, at fixed halo mass, no statistical difference in these properties between our kinematically aligned and misaligned galaxies. We suggest that the lack of correlation could be indicative of cooling flows from the hot halo playing a far larger role than `cold mode' accretion from the cosmic web or a demonstration that the spatial extent of current large-scale integral field unit (IFU) surveys hold little information about large-scale environment extractable through this method.

Galaxy mass and environment play a major role in the evolution of galaxies. In the transition from star-forming to quenched galaxies, Active galactic nuclei (AGN) have also a principal action. However, the connections between these three actors are still uncertain. In this work we investigate the effects of stellar mass and the large-scale environment (LSS), on the fraction of optical nuclear activity in a population of isolated galaxies, where AGN would not be triggered by recent galaxy interactions or mergers. As a continuation of a previous work, we focus on isolated galaxies to study the effect of stellar mass and the LSS in terms of morphology (early- and late-type), colour (red and blue), and specific star formation rate (quenched and star-forming). To explore where AGN activity is affected by the LSS we fix the stellar mass into low- and high-mass galaxies. We use the tidal strength parameter to quantify their effects. We found that AGN is strongly affected by stellar mass in 'active' galaxies (namely late-type, blue, and star-forming), however it has no influence for 'quiescent' galaxies (namely early-type, red, and quenched), at least for masses down to $\rm 10^{10}\,[M_\odot]$. In relation to the LSS, we found an increment on the fraction of SFN with denser LSS in low-mass star forming and red isolated galaxies. Regarding AGN, we find a clear increment of the fraction of AGN with denser environment in quenched and red isolated galaxies, independently of the stellar mass. AGN activity would be 'mass triggered' in 'active' isolated galaxies. This means that AGN is independent of the intrinsic property of the galaxies, but on its stellar mass. On the other hand, AGN would be 'environment triggered' in 'quiescent' isolated galaxies, where the fraction of AGN in terms of sSFR and colour increases from void regions to denser LSS, independently of its stellar mass.

The formation and evolution of elliptical galaxies in low-density environments are less understood than classical elliptical galaxies in high-density environments. Isolated galaxies are defined as galaxies without massive neighbors within scales of galaxy groups. The effect of the environment at several Mpc scales on their properties has been barely explored. Here we study the role of large-scale environment in some physical properties of 573 isolated elliptical galaxies out to z=0.08. We use three environmental estimators of the large-scale structure within a projected radius of 5 Mpc around isolated galaxies: the tidal strength parameter, the projected density eta_k, and the distance to the fifth nearest neighbor galaxy. We find 80% of galaxies at lower densities correspond to 'red and dead' elliptical galaxies. Blue and red galaxies do not tend to be located in different environments according to eta_k. Almost all the isolated ellipticals in the densest large-scale environments are red or quenched, where a third of them are low-mass galaxies. The percentage of isolated elliptical galaxies located in the AGN region of the BPT diagram is 64%. We have identified 33 blue, star-forming isolated ellipticals using both color and sSFR. Half of them are star-forming nuclei in the BPT diagram, which is 5% of the galaxies in this diagram. The large-scale environment is not playing the primary role to determine the color or sSFR of isolated elliptical galaxies. The large-scale environment seems to be negligible from a stellar mass scale around 10^10.6 Msun, probably because of the dominant presence of AGN at higher masses. For lower masses, the processes of cooling and infall of gas from large scales are very inefficient in ellipticals. AGN might also be an essential ingredient to keep most of the low-mass isolated elliptical galaxies quenched.

Understanding the formation and evolution of early-type, spheroid-dominated galaxies is an open question within the context of the hierarchical clustering scenario, particularly, in low-density environments. Our goal is to study the main structural, dynamical, and stellar population properties and assembly histories of field spheroid-dominated galaxies formed in a LCDM scenario to assess to what extend they are consistent with observations. We selected spheroid-dominated systems from a LCDM simulation that includes star formation, chemical evolution and Supernova feedback. A sample of 18 field systems with Mstar <= 6x10^10 Msun that are dominated by the spheroid component. For this sample we estimate the fundamental relations of ellipticals and then compared with current observations. The simulated spheroid galaxies have sizes in good agreement with observations. The bulges follow a Sersic law with Sersic indexes that correlate with the bulge-to-total mass ratios. The structural-dynamical properties of the simulated galaxies are consistent with observed Faber-Jackson, Fundamental Plane, and Tully-Fisher relations. However, the simulated galaxies are bluer and with higher star formation rates than observed isolated early-type galaxies. The archaeological mass growth histories show a slightly delayed formation and more prominent inside-out growth mode than observational inferences based on the fossil record method. The main structural and dynamical properties of the simulated spheroid-dominated galaxies are consistent with observations. This is remarkable since none of them has been tuned to be reproduced. However, the simulated galaxies are blue and star-forming, and with later stellar mass growth histories as compared to observational inferences. This is mainly due to the persistence of extended discs in the simulations. Abridged

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.

Galaxy integrated H{\alpha} star formation rate-stellar mass relation, or SFR(global)-M*(global) relation, is crucial for understanding star formation history and evolution of galaxies. However, many studies have dealt with SFR using unresolved measurements, which makes it difficult to separate out the contamination from other ionizing sources, such as active galactic nuclei and evolved stars. Using the integral field spectroscopic observations from SDSS-IV MaNGA, we spatially disentangle the contribution from different H{\alpha} powering sources for ~1000 galaxies. We find that, when including regions dominated by all ionizing sources in galaxies, the spatially-resolved relation between H{\alpha} surface density ({\Sigma}H{\alpha}(all)) and stellar mass surface density ({\Sigma}*(all)) progressively turns over at high {\Sigma}*(all) end for increasing M*(global) and bulge dominance (bulge-to-total light ratio, B/T). This in turn leads to the flattening of the integrated H{\alpha}(global)-M*(global) relation in the literature. By contrast, there is no noticeable flattening in both integrated H{\alpha}(HII)-M*(HII) and spatially-resolved {\Sigma}H{\alpha}(HII)-{\Sigma}*(HII) relations when only regions where star formation dominates the ionization are considered. In other words, the flattening can be attributed to the increasing regions powered by non-star-formation sources, which generally have lower ionizing ability than star formation. Analysis of the fractional contribution of non-star-formation sources to total H{\alpha} luminosity of a galaxy suggests a decreasing role of star formation as an ionizing source toward high-mass, high-B/T galaxies and bulge regions. This result indicates that the appearance of the galaxy integrated SFR-M* relation critically depends on their global properties (M*(global) and B/T) and relative abundances of various ionizing sources within the galaxies.

We study the correlation between the specific star formation rate of central galaxies and neighbour galaxies, also known as 'galactic conformity', out to 20 Mpc/h using three semi-analytic models (SAMs, one from L-GALAXIES and other two from GALFORM). The aim is to establish whether SAMs are able to show galactic conformity using different models and selection criteria. In all the models, when the selection of primary galaxies is based on an isolation criterion in real space, the mean fraction of quenched galaxies around quenched primary galaxies is higher than that around star-forming primary galaxies of the same stellar mass. The overall signal of conformity decreases when we remove satellites selected as primary galaxies, but the effect is much stronger in GALFORM models compared with the L-GALAXIES model. We find this difference is partially explained by the fact that in GALFORM once a galaxy becomes a satellite remains as such, whereas satellites can become centrals at a later time in L-GALAXIES. The signal of conformity decreases down to 60% in the L-GALAXIES model after removing central galaxies that were ejected from their host halo in the past. Galactic conformity is also influenced by primary galaxies at fixed stellar mass that reside in dark matter haloes of different masses. Finally, we explore a proxy of conformity between distinct haloes. In this case the conformity is weak beyond ~ 3 Mpc/h (<3% in L-GALAXIES, <1-2% in GALFORM models). Therefore, it seems difficult that conformity is directly related with a long-range effect.

Regions of disc galaxies with widespread star formation tend to be both gravitationally unstable and self-shielded against ionizing radiation, whereas extended outer discs with little or no star formation tend to be stable and unshielded on average. We explore what drives the transition between these two regimes, specifically whether discs first meet the conditions for self-shielding (parameterized by dust optical depth, $\tau$) or gravitational instability (parameterized by a modified version of Toomre's instability parameters, $Q_{\rm thermal}$, which quantifies the stability of a gas disc that is thermally supported at $T=10^4$ K). We first introduce a new metric formed by the product of these quantities, $Q_{\rm thermal}\tau$, which indicates whether the conditions for disk instability or self-shielding are easier to meet in a given region of a galaxy, and we discuss how $Q_{\rm thermal}\tau$ can be constrained even in the absence of direct gas information. We then analyse a sample of 13 galaxies with resolved gas measurements and find that on average galaxies will reach the threshold for disk instabilities ($Q_{\rm thermal}<1$) before reaching the threshold for self-shielding ($\tau>1$). Using integral field spectroscopic observations of a sample of 236 galaxies from the MaNGA survey, we find that the value of $Q_{\rm thermal}\tau$ in star-forming discs is consistent with similar behavior. These results support a scenario where disc fragmentation and collapse occurs before self-shielding, suggesting that gravitational instabilities are the primary condition for widespread star formation in galaxy discs. Our results support similar conclusions based on recent galaxy simulations.

We study the role of cold gas in quenching star formation in the green valley by analysing ALMA $^{12}$CO (1-0) observations of three galaxies with resolved optical spectroscopy from the MaNGA survey. We present resolution-matched maps of the star formation rate and molecular gas mass. These data are used to calculate the star formation efficiency (SFE) and gas fraction ($f_{\rm~gas}$) for these galaxies separately in the central `bulge' regions and outer disks. We find that, for the two galaxies whose global specific star formation rate (sSFR) deviates most from the star formation main sequence, the gas fraction in the bulges is significantly lower than that in their disks, supporting an `inside-out' model of galaxy quenching. For the two galaxies where SFE can be reliably determined in the central regions, the bulges and disks share similar SFEs. This suggests that a decline in $f_{\rm~gas}$ is the main driver of lowered sSFR in bulges compared to disks in green valley galaxies. Within the disks, there exist common correlations between the sSFR and SFE and between sSFR and $f_{\rm~gas}$ on kpc scales -- the local SFE or $f_{\rm~gas}$ in the disks declines with local sSFR. Our results support a picture in which the sSFR in bulges is primarily controlled by $f_{\rm~gas}$, whereas both SFE and $f_{\rm~gas}$ play a role in lowering the sSFR in disks. A larger sample is required to confirm if the trend established in this work is representative of green valley as a whole.

The fourth generation of the Sloan Digital Sky Survey (SDSS-IV) began observations in July 2014. It pursues three core programs: APOGEE-2, MaNGA, and eBOSS. In addition, eBOSS contains two major subprograms: TDSS and SPIDERS. This paper describes the first data release from SDSS-IV, Data Release 13 (DR13), which contains new data, reanalysis of existing data sets and, like all SDSS data releases, is inclusive of previously released data. DR13 makes publicly available 1390 spatially resolved integral field unit observations of nearby galaxies from MaNGA, the first data released from this survey. It includes new observations from eBOSS, completing SEQUELS. In addition to targeting galaxies and quasars, SEQUELS also targeted variability-selected objects from TDSS and X-ray selected objects from SPIDERS. DR13 includes new reductions of the SDSS-III BOSS data, improving the spectrophotometric calibration and redshift classification. DR13 releases new reductions of the APOGEE-1 data from SDSS-III, with abundances of elements not previously included and improved stellar parameters for dwarf stars and cooler stars. For the SDSS imaging data, DR13 provides new, more robust and precise photometric calibrations. Several value-added catalogs are being released in tandem with DR13, in particular target catalogs relevant for eBOSS, TDSS, and SPIDERS, and an updated red-clump catalog for APOGEE. This paper describes the location and format of the data now publicly available, as well as providing references to the important technical papers that describe the targeting, observing, and data reduction. The SDSS website, http://www.sdss.org, provides links to the data, tutorials and examples of data access, and extensive documentation of the reduction and analysis procedures. DR13 is the first of a scheduled set that will contain new data and analyses from the planned ~6-year operations of SDSS-IV.

The Large Synoptic Survey Telescope (LSST) will enable revolutionary studies of galaxies, dark matter, and black holes over cosmic time. The LSST Galaxies Science Collaboration has identified a host of preparatory research tasks required to leverage fully the LSST dataset for extragalactic science beyond the study of dark energy. This Galaxies Science Roadmap provides a brief introduction to critical extragalactic science to be conducted ahead of LSST operations, and a detailed list of preparatory science tasks including the motivation, activities, and deliverables associated with each. The Galaxies Science Roadmap will serve as a guiding document for researchers interested in conducting extragalactic science in anticipation of the forthcoming LSST era.

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.

We have studied a sample of 89 very isolated, elliptical galaxies at z < 0.08 and compared their properties with elliptical galaxies located in a high-density environment such as the Coma supercluster. Our aim is to probe the role of environment on the morphological transformation and quenching of elliptical galaxies as a function of mass. In addition, we elucidate the nature of a particular set of blue and star-forming isolated ellipticals identified here. We study physical properties of ellipticals such as color, specific star formation rate, galaxy size, and stellar age, as a function of stellar mass and environment based on SDSS data. We analyze the blue star-forming isolated ellipticals in more detail, through photometric characterization using GALFIT, and infer their star formation history using STARLIGHT. Among the isolated ellipticals ~ 20% are blue, 8% are star forming, and ~ 10% are recently quenched, while among the Coma ellipticals ~ 8% are blue and just <= 1% are star forming or recently quenched. There are four isolated galaxies (~ 4.5%) that are blue and star forming at the same time. These galaxies, with masses between 7 x 10^9 and 2 x 10^10 h-2 M_sun, are also the youngest galaxies with light-weighted stellar ages <= 1 Gyr and exhibit bluer colors toward the galaxy center. Around 30-60% of their present-day luminosity, but only < 5% of their present-day mass, is due to star formation in the last 1 Gyr. The processes of morphological transformation and quenching seem to be in general independent of environment since most of elliptical galaxies are 'red and dead', although the transition to the red sequence should be faster for isolated ellipticals. In some cases, the isolated environment seems to propitiate the rejuvenation of ellipticals by recent (< 1 Gyr) cold gas accretion.

We present a new color-selection technique, based on the Bruzual & Charlot models convolved with the bands of the ALHAMBRA survey, and the redshifted position of the Balmer jump to select star-forming galaxies in the redshift range 0.5 < z < 1.5. These galaxies are dubbed Balmer jump Galaxies BJGs. We apply the iSEDfit Bayesian approach to fit each detailed SED and determine star-formation rate (SFR), stellar mass, age and absolute magnitudes. The mass of the haloes where these samples reside are found via a clustering analysis. Five volume-limited BJG sub-samples with different mean redshifts are found to reside in haloes of median masses $\sim 10^{12.5 \pm 0.2} M_\odot$ slightly increasing toward z=0.5. This increment is similar to numerical simulations results which suggests that we are tracing the evolution of an evolving population of haloes as they grow to reach a mass of $\sim 10^{12.7 \pm 0.1} M_\odot$ at z=0.5. The likely progenitors of our samples at z$\sim$3 are Lyman Break Galaxies, which at z$\sim$2 would evolve into star-forming BzK galaxies, and their descendants in the local Universe are elliptical galaxies.Hence, this allows us to follow the putative evolution of the SFR, stellar mass and age of these galaxies. From z$\sim$1.0 to z$\sim$0.5, the stellar mass of the volume limited BJG samples nearly does not change with redshift, suggesting that major mergers play a minor role on the evolution of these galaxies. The SFR evolution accounts for the small variations of stellar mass, suggesting that star formation and possible minor mergers are the main channels of mass assembly.

We analyse galaxies in groups in the Sloan Digital Sky Survey (SDSS) and find a weak but significant assembly-type bias, where old central galaxies have a higher clustering amplitude (61 $\pm$ 9 per cent) at scales > 1 Mpc than young central galaxies of equal host halo mass ($M_{h} \sim 10^{11.8} h^{-1}$ $M_{\odot}$). The observational sample is volume-limited out to z=0.1 with $M_r -$ 5 log$(h) \le -19.6$. We construct a mock catalogue of galaxies that shows a similar signal of assembly bias (46 $\pm$ 9 per cent) at the same halo mass. We then adapt the model presented by Lacerna & Padilla (Paper I) to redefine the overdensity peak height, which traces the assembly bias such that galaxies in equal density peaks show the same clustering regardless of their stellar age, but this time using observational features such as a flux limit. The proxy for peak height, which is proposed as a better alternative than the virial mass, consists in the total mass given by the mass of neighbour host haloes in cylinders centred at each central galaxy. The radius of the cylinder is parametrized as a function of stellar age and virial mass. The best-fitting set of parameters that make the assembly bias signal lower than 5$-$15 per cent for both SDSS and mock central galaxies are similar. The idea behind the parametrization is not to minimize the bias, but it is to use this method to understand the physical features that produce the assembly bias effect. Even though the tracers of the density field used here differ significantly from those used in paper I, our analysis of the simulated catalogue indicates that the different tracers produce correlated proxies, and therefore the reason behind this assembly bias is the crowding of peaks in both simulations and the SDSS.

We perform an exhaustive comparison among central galaxies from SDSS catalogs in different local environments at 0.01<=z<=0.08. The central galaxies are separated into two categories: group centrals (host halos containing satellites) and field centrals (host halos without satellites). From the latter, we select other two subsamples: isolated centrals and bright field centrals, both with the same magnitude limit. The stellar mass (Ms) distributions of the field and group central galaxies are different, which explains why in general the field central galaxies are mainly located in the blue cloud/star forming regions, whereas the group central galaxies are strongly biased to the red sequence/passive regions. The isolated centrals occupy the same regions as the bright field centrals since both populations have similar Ms distributions. At parity of Ms, the color and specific star formation rate (sSFR) distributions of the samples are similar, specially between field and group centrals. Furthermore, we find that the stellar-to-halo mass (Ms-Mh) relation of isolated galaxies does not depend on the color, sSFR and morphological type. For systems without satellites, the Ms-Mh relation steepens at high halo masses compared to group centrals, which is a consequence of assuming a one-to-one relation between group total stellar mass and halo mass. Under the same assumption, the scatter around the Ms-Mh relation of centrals with satellites increases with halo mass. Our results suggest that the mass growth of central galaxies is mostly driven by the halo mass, with environment and mergers playing a secondary role.

The mass aggregation and merger histories of present-day distinct haloes selected from the cosmological Millennium Simulations I and II are mapped into stellar mass aggregation and galaxy merger histories of central galaxies by using empirical stellar-to-halo and stellar-to-gas mass relations. The growth of bulges driven by the galaxy mergers/interactions is calculated using dynamical prescriptions. The predicted bulge demographics at redshift z~0 is consistent with observations (Zavala+2012). Here we present the evolution of the morphological mix (traced by the bulge-to-total mass ratio, B/T) as a function of mass up to z=3. This mix remains qualitatively the same up to z~1: B/T<0.1 galaxies dominate at low masses, 0.1<B/T<0.45 at intermediate masses, and B/T>0.45 at large masses. At z>1, the fractions of disc-dominated and bulgeless galaxies increase strongly, and by z~2 the era of pure disc galaxies is reached. Bulge-dominated galaxies acquire such a morphology, and most of their mass, following a downsizing trend. Since our results are consistent with most of the recent observational studies of the morphological mix at different redshifts, a LCDM-based scenario of merger-driven bulge assembly does not seem to face critical issues. However, if the stellar-to-halo mass relation changes too little with redshift, then some tensions with observations appear.

We present Gemini/GMOS-S multi-object spectroscopy of 31 galaxy cluster candidates at redshifts between 0.2 and 1.0 and centered on QSO sight-lines taken from Lopez et al. (2008). The targets were selected based on the presence of a intervening MgII absorption system at a similar redshift to that of a galaxy cluster candidate lying at a projected distance < 2 h^{-1}Mpc from the QSO sight-line (a 'photometric-hit'). The absorption systems span rest-frame equivalent widths between 0.015 and 2.028 angstroms. Our aim was 3-fold: 1) identify the absorbing galaxies and determine their impact parameters, 2) confirm the galaxy cluster candidates in the vicinity of each quasar sightline, and 3) determine whether the absorbing galaxies reside in galaxy clusters. Our main findings are: 1) the identification of 10 out of 24 absorbing galaxies with redshifts up to 1.0955. 2) The spectroscopic confirmation of 20 out of 31 cluster/group candidates, with most of the confirmed clusters/groups at z < 0.7. 3) Following from the results above, the spectroscopic confirmation of 10 out of 14 photometric hits within ~ 650 km/s from galaxy clusters/groups, in addition to 2 new ones related to galaxy group environments. These numbers imply efficiencies of 71% in finding such systems with MOS spectroscopy. This is a remarkable result since we defined a photometric hit as those cluster-absorber pairs having a redshift difference dz = 0.1. Absorbing cluster-galaxies hosting weak absorbers are consistent with lower star formation activity than the rest, which produce strong absorption and agree with typical MgII absorbing galaxies found in the literature. Our spectroscopic confirmations lend support to the selection of photometric hits made in Lopez et al. (2008).

We study an assembly-type bias parametrized by the dimensionless spin parameter that affects massive structures. In numerical simulations higher spin haloes are more strongly clustered than lower spin haloes of equal mass. We detect a difference of over a 30 per cent in the clustering strength for dark matter haloes of 10^13-10^14 Msun, which is similar to the result of Bett et al. We explore whether the dependence of clustering strength on halo spin is removed if we apply the redefinition of overdensity peak height proposed by Lacerna & Padilla (Paper I) obtained using assembly ages. We find that this is not the case due to two reasons. Firstly, only a few objects of low-virial mass are moved into the mass range where the spin introduces an assembly bias after using this redefinition. Secondly, this formalism does not alter the mass of massive objects. We then repeat the process of finding the redefined peak height of Paper I but using the spin. In this case, the new masses show no spin-related assembly bias but they introduce a previously absent assembly bias with respect to relative age. From this result, we conclude that the assembly-type bias with respect to the halo spin has a different origin than with respect to assembly age. The former may be due to the material from filaments, which is accreted by massive haloes, that is enhanced in high-density environments, thus causing more extreme spin values without significantly changing the formation age of the halo. In addition, high-mass objects may correspond, in some cases, to a different peak height than that suggested by their mass in numerical simulations, providing a possible explanation for the assembly bias with respect to spin. (abridged)

We present a new proxy for the overdensity peak height for which the large-scale clustering of haloes of a given mass does not vary significantly with the assembly history. The peak height, usually taken to be well represented by the virial mass, can instead be approximated by the mass inside spheres of different radii, which in some cases can be larger than the virial radius and therefore include mass outside the individual host halo. The sphere radii are defined as r = $a$ delta_t + $b$ log_10(M_vir/M_nl), where delta_t is the age relative to the typical age of galaxies hosted by haloes with virial mass M_vir, M_nl is the non-linear mass, and $a$=0.2 and $b$=-0.02 are the free parameters adjusted to trace the assembly bias effect. Note that $r$ depends on both halo mass and age. In this new approach, some of the objects which were initially considered low-mass peaks belong to regions with higher overdensities. At large scales, i.e. in the two-halo regime, this model properly recovers the simple prescription where the bias responds to the height of the mass peak alone, in contrast to the usual definition (virial mass) that shows a strong dependence on additional halo properties such as formation time. The dependence on the age in the one-halo term is also remarkably reduced. The population of galaxies whose "peak height" changes with this new definition consists mainly of old stellar populations and are preferentially hosted by low-mass haloes located near more massive objects. The latter is in agreement with recent results which indicate that old, low-mass haloes would suffer truncation of mass accretion by nearby larger haloes or simply due to the high density of their surroundings, thus showing an assembly bias effect. The change in mass is small enough that the Sheth et al. (2001) mass function is still a good fit to the resulting distribution of new masses.

We present the results of a search for variable stars in the globular cluster NGC 5286, which has recently been suggested to be associated with the Canis Major dwarf spheroidal galaxy. 57 variable stars were detected, only 19 of which had previously been known. Among our detections one finds 52 RR Lyrae (22 RRc and 30 RRab), 4 LPV's, and 1 type II Cepheid of the BL Herculis type. Periods are derived for all of the RR Lyrae as well as the Cepheid, and BV light curves are provided for all the variables. The mean period of the RRab variables is <Pab> = 0.656 days, and the number fraction of RRc stars is N(c)/N(RR) = 0.42, both consistent with an Oosterhoff II (OoII) type -- thus making NGC 5286 one of the most metal-rich ([Fe/H] = -1.67; Harris 1996) OoII globulars known to date. The minimum period of the \RRab's, namely Pab,min = 0.513 d, while still consistent with an OoII classification, falls towards the short end of the observed Pab,min distribution for OoII globular clusters. As was recently found in the case of the prototypical OoII globular cluster M15 (NGC 7078), the distribution of stars in the Bailey diagram does not strictly conform to the previously reported locus for OoII stars. We provide Fourier decomposition parameters for all of the RR Lyrae stars detected in our survey, and discuss the physical parameters derived therefrom. The values derived for the RRc's are not consistent with those typically found for OoII clusters, which may be due to the cluster's relatively high metallicity -- the latter being confirmed by our Fourier analysis of the ab-type RR Lyrae light curves. We derive for the cluster a revised distance modulus of (m-M)V = 16.04 mag. (ABRIDGED)

Based on recent results on the frequency of MgII absorption line systems in the "QSO behind RCS clusters" survey (QbC), we analyse the effects of the cluster environment on the sizes of baryonic haloes around galaxies. We use two independent models, i) an empirical halo occupation model which fits current measurements of the clustering and luminosity function of galaxies at low and high redshifts, and ii) the GALFORM semi-analytic model of galaxy formation, which follows the evolution of the galaxy population from first principles, adjusted to match the statistics of low and high redshift galaxies. In both models we constrain the MgII halo sizes of field and cluster galaxies using observational results on the observed MgII statistics. Our results for the field are in good agreement with previous works, indicating a typical \mgii\ halo size of $r_MgII ~ 50h_71^-1kpc in the semi-analytic model, and slightly lower in the halo occupation number approach. For the cluster environment, we find that both models require a median MgII halo size of r_MgII< 10h_71^-1kpc in order to reproduce the observed statistics on absorption line systems in clusters of galaxies. Based on the Chen & Tinker (2008) result that stronger systems occur closer to the MgII halo centre, we find that strong absorption systems in clusters of galaxies occur at roughly a fixed fraction of the cold-warm halo size out to 1h_71^-1Mpc from the cluster centres. In contrast, weaker absorption systems appear to occur at progressively shorter relative fractions of this halo as the distance to the cluster centre decreases.