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 build a volume-limited sample of galaxies derived from the SDSS-DR7 to characterize several physical properties of the dark matter halos where LSB galaxies reside. Using an observational proxy for the assembly time, we found that LSB galaxies assembly half of their total halo mass later than HSB ones, reinforcing the idea of them being unevolved systems. We use 5 different methods to estimate the total halo mass, finding that the total stellar-to-halo mass ratio is up to 22% lower in LSB galaxies. Finally, in order to estimate the spin parameter, we use a bulge+disk decomposition to obtain the specific angular momentum $j_*$ of the galaxy, a Tully-Fisher relation to estimate the rotation velocity of the disk, and the 5 different estimations of the halo mass to calculate the spin parameter. We found that the spin of LSB galaxies is 1.2 to 2 times higher than for HSB ones. We compare these results with a control sample that includes kinematic information, taken from the ALFALFA $\alpha$.100 galaxy catalog, allowing us to measure directly the rotation velocity of the disk. The trends in the values of $j_*$ and $\lambda$ are similar to the volume-limited sample.

We select a volume-limited sample of galaxies derived from the SDSS-DR7 to study the environment of low surface brightness (LSB) galaxies at different scales, as well as several physical properties of the dark matter haloes where the LSB galaxies of the sample are embedded. To characterize the environment we make use of a number of publicly available value-added galaxy catalogues. We find a slight preference for LSB galaxies to be found in filaments instead of clusters, with their mean distance to the nearest filament typically larger than for high surface brightness (HSB) galaxies. The fraction of isolated central LSB galaxies is higher than the same fraction for HSB ones, and the density of their local environment lower. The stellar-to-halo mass ratio using four different estimates is up to $\sim$20% for HSB galaxies. LSB central galaxies present more recent assembly times when compared with their HSB counterparts. Regarding the $\lambda$ spin parameter, using six different proxies for its estimation, we find that LSB galaxies present systematically larger values of $\lambda$ than the HSB galaxy sample, and constructing a control sample with direct kinematic information drawn from ALFALFA, we confirm that the spin parameter of LSB galaxies is 1.6 to 2 times larger than the one estimated for their HSB counterparts.

We present a Bayesian method to determine simultaneously the age, metallicity, distance modulus, and interstellar reddening by dust of any resolved stellar population, by comparing the observed and synthetic color magnitude diagrams on a star by star basis, with no need to bin the data into a carefully selected magnitude grid. We test the method with mock stellar populations, and show that it works correctly even for scarce stellar populations with only one or two hundred stars above the main sequence turn off. If the population is the result of two star formation bursts, we can infer the contribution of each event to the total stellar population. The code works automatically and has already been used to study massive amounts of Magellanic clouds photometric data. In this paper we analyze in detail three Large Magellanic Cloud star clusters and 6 Ultra Faint Dwarf Galaxies. For these galaxies we recover physical parameters in agreement with those quoted in the literature, age $\sim13.7$ Gyr and a very low metallicity $log\,Z\sim-4$. Searching for multiple populations in these galaxies, we find, at a very low significance level, signs of a double stellar population for Ursa Major I: a dominant old population and a younger one which contributes $\sim25$% of the stars, in agreement with independent results from other authors.

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 a study of the integrated properties of the 835 galaxies in the CALIFA survey. To derive the main physical parameters of the galaxies we have fitted their UV-to-IR spectral energy distributions (SED) with sets of theoretical models using CIGALE. We perform a comparison of the integrated galaxy parameters derived from multi-band SED fitting with those obtained from modelling the Integral Field Unit (IFU) spectra and show the clear advantage of using the SED-derived star formation rates (SFR). A detailed analysis of galaxies in the SFR/Mstar plane as a function of their properties reveals that quenching of star formation is caused by a combination of gas deficiency and the inefficiency of the existing gas to form new stars. Exploring the plausible mechanisms that could produce this effect, we find a strong correlation with galaxy morphology and the build-up of central bulge. On the other hand, the presence of AGN and/or a stellar bar, as well as the local environment have only temporal effects on the current star formation, a result also consistent with their model-derived star formation histories.

Pressure supported systems modeled under MONDian extended gravity are expected to show an outer flattening in their velocity dispersion profiles. A characteristic scaling between the amplitude of the asymptotic velocity dispersion and the radius at which the flattening occurs is also expected. By comprehensively analyzing the dynamical behavior of $\sim$300 extremely low rotating elliptical galaxies from the MaNGA survey, we show this type of pressure supported systems to be consistent with MONDian expectations, for a range of central velocity dispersion values of $60 km/s < \sigma_{central}< 280 km/s $ and asymptotic velocity dispersion values of $28km/s < \sigma_{\infty}<250 km/s$. We find that a universal velocity dispersion profile accurately describes the studied systems; the predicted kinematics of extended gravity are verified for all well observed galaxies.

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.

A large fraction of early-type galaxies (ETGs) host prominent dust features, and central dust rings are arguably the most interesting among them. We present here `Lord Of The Rings' (LOTR), a new methodology which allows to integrate the extinction by dust rings in a 2D fitting modelling of the surface brightness distribution. Our pipeline acts in two steps, first using the surface fitting software GALFIT to determine the unabsorbed stellar emission, and then adopting the radiative transfer code SKIRT to apply dust extinction. We apply our technique to NGC 4552 and NGC 4494, two nearby ETGs. We show that the extinction by a dust ring can mimic, in a surface brightness profile, a central point source (e.g. an unresolved nuclear stellar cluster or an active galactic nucleus; AGN) superimposed to a `core' (i.e. a central flattening of the stellar light commonly observed in massive ETGs). We discuss how properly accounting for dust features is of paramount importance to derive correct fluxes especially for low luminosity AGNs (LLAGNs). We suggest that the geometries of dust features are strictly connected with how relaxed is the gravitational potential, i.e. with the evolutionary stage of the host galaxy. Additionally, we find hints that the dust mass contained in the ring relates to the AGN activity.

The so-called jellyfish galaxies are objects exhibiting disturbed morphology, mostly in the form of tails of gas stripped from the main body of the galaxy. Several works have strongly suggested ram pressure stripping to be the mechanism driving this phenomenon. Here, we focus on one of these objects, drawn from a sample of optically selected jellyfish galaxies, and use it to validate SINOPSIS, the spectral fitting code that will be used for the analysis of the GASP (GAs Stripping Phenomena in galaxies with MUSE) survey, and study the spatial distribution and physical properties of gas and stellar populations in this galaxy. We compare the model spectra to those obtained with GANDALF, a code with similar features widely used to interpret the kinematic of stars and gas in galaxies from IFU data. We find that SINOPSIS can reproduce the pixel-by-pixel spectra of this galaxy at least as good as GANDALF does, providing reliable estimates of the underlying stellar absorption to properly correct the nebular gas emission. Using these results, we find strong evidences of a double effect of ram pressure exerted by the intracluster medium onto the gas of the galaxy. A moderate burst of star formation, dating between 20 and 500 Myr ago and involving the outer parts of the galaxy more strongly than the inner regions, was likely induced by a first interaction of the galaxy with the intracluster medium. Stripping by ram pressure, plus probable gas depletion due to star formation, contributed to create a truncated ionized gas disk. The presence of an extended stellar tail on only one side of the disk, points instead to another kind of process, likely a gravitational interaction by a fly-by or a close encounter with another galaxy in the cluster.

With the aim of assessing if low surface brightness galaxies host stellar bars, and study the dependence of the occurrence of bars as a function of surface brightness, we use the Galaxy Zoo 2 dataset to construct a large volume-limited sample of galaxies, and segregate the galaxies as low and high surface brightness in terms of their central surface brightness. We find that the fraction of low surface brightness galaxies hosting strong bars is systematically lower than the one found for high surface brightness galaxies. The dependence of the bar fraction on the central surface brightness is mostly driven by a correlation of the surface brightness with the spin and the gas-richness of the galaxies, showing only a minor dependence on the surface brightness. We also find that the length of the bars shows a strong dependence on the surface brightness, and although some of this dependence is attributed to the gas content, even at fixed gas-to-stellar mass ratio, high surface brightness galaxies host longer bars than their low surface brightness counterparts, which we attribute to an anticorrelation of the surface brightness with the spin.

We study the star formation activity of nearby galaxies with bars using a sample of late-type galaxies at 0.02$\leq z \leq$ 0.05489 and $M_r <-19.5$ from the Sloan Digital Sky Survey. We compare the physical properties of strongly and weakly barred galaxies with those of non-barred galaxies that have stellar mass and redshift distributions similar to barred galaxies. We find that the star formation activity of strongly barred galaxies probed by starburstiness, $\it{g-r}$, NUV$-r$, and mid-infrared [3.4]$-$[12] colors is, on average, lower than that of non-barred galaxies. However, weakly barred galaxies do not show such a difference between barred and non-barred galaxies. The amounts of atomic and molecular gas in strongly barred galaxies are smaller than those of non-barred galaxies, and the gas metallicity is higher in strongly barred galaxies than in non-barred galaxies. The gas properties of weakly barred galaxies again show no difference from those of non-barred galaxies. We stack the optical spectra of barred and non-barred galaxies in several mass bins and fit to the stacked spectra with a spectral fitting code, STARLIGHT. We find no significant difference in stellar populations between barred and non-barred galaxies for both strongly and weakly barred galaxies. Our results are consistent with the idea that the star formation activity of barred galaxies is enhanced in the past along with significant gas consumption, and is currently lower than or similar to that of non-barred galaxies. The past star formation enhancement depends on the strength of bars.

For any MONDian extended theory of gravity where the rotation curves of spiral galaxies are explained through a change in physics rather than the hypothesis of dark matter, a generic dynamical behavior is expected for pressure supported systems: an outer flattening of the velocity dispersion profile occurring at a characteristic radius, where both the amplitude of this flat velocity dispersion and the radius at which it appears are predicted to show distinct scalings with the total mass of the system. By carefully analyzing the dynamics of globular clusters and elliptical galaxies, we are able to significantly extend the astronomical diversity of objects in which MONDian gravity has been tested, from spiral galaxies, to the much larger mass range covered by pressure supported systems. We show that a universal projected velocity dispersion profile accurately describes various classes of pressure supported systems, and further, that the expectations of extended gravity are met, across seven orders of magnitude in mass. These observed scalings are not expected under dark matter cosmology, and would require particular explanations tuned at the scales of each distinct astrophysical system.

We select a sample of galaxies from the Sloan Digital Sky Survey Data Release 7 (SDSS-DR7) where galaxies are classified, through visual inspection, as hosting strong bars, weak bars or as unbarred galaxies, and make use of HI mass and kinematic information from the Arecibo Legacy Fast ALFA (ALFALFA) survey catalog, to study the stellar, atomic gas and dark matter content of barred disk galaxies. We find, in agreement with previous studies, that the bar fraction increases with increasing stellar mass. A similar trend is found with total baryonic mass, although the dependence is not as strong as with stellar mass, this due to the contribution of gas. The bar fraction shows a decrease with increasing gas mass fraction. This anticorrelation between the likelihood of a galaxy hosting a bar with the gas richness of the galaxy results from the inhibiting effect the gas has in the formation of bars. We also find that for massive galaxies with stellar masses larger than 10$^{10} M_{\odot}$, at fixed stellar mass, the bar fraction decreases with increasing global halo mass (i.e. halo mass measured up to a radius of the order of the HI disk extent).

We present a study of the environment of barred galaxies using galaxies drawn from the SDSS. We use several different statistics to quantify the environment: the projected two-point cross-correlation function, the background-subtracted number count of neighbor galaxies, the overdensity of the local environment, the membership of our galaxies to galaxy groups to segregate central and satellite systems, and for central galaxies we estimate the stellar to halo mass ratio (M$_{\mathrm{*}}/$M$_{\mathrm{h}}$) . When we split our sample into early- and late-type galaxies, we see a weak but significant trend for early-type galaxies with a bar to be more strongly clustered on scales from a few 100 kpc to 1 Mpc when compared to unbarred early-type galaxies. This indicates that the presence of a bar in early-type galaxies depends on the location within their host dark matter halos. This is confirmed by the group catalog in the sense that for early-types, the fraction of central galaxies is smaller if they have a bar. For late-type galaxies, we find fewer neighbors within $\sim50$ kpc around the barred galaxies when compared to unbarred galaxies from the control sample, suggesting that tidal forces from close companions suppress the formation/growth of bars. For central late-type galaxies, bars are more common on galaxies with high M$_{\mathrm{*}}/$M$_{\mathrm{h}}$ values, as expected from early theoretical works which showed that systems with massive dark matter halos are more stable against bar instabilities. Finally, we find no obvious correlation between overdensity and the bars in our sample, showing that galactic bars are not obviously linked to the large-scale structure of the universe.

We use a large volume-limited sample of disk galaxies drawn from the Sloan Digital Sky Survey Data Release 7 to study the dependence of the bar fraction on the stellar-to-halo mass ratio, making use of a group catalog, we identify central and satellite galaxies in our sample. For the central galaxies in the sample we estimate the stellar-to-halo mass ratio (M$_{\mathrm{*}}/$M$_{\mathrm{h}}$) and find that the fraction of barred galaxies is a strong function of this ratio, especially for the case of strong bars. Bars are more common in galaxies with high M$_{\mathrm{*}}/$M$_{\mathrm{h}}$ values, as expected from early theoretical works that showed that systems with massive dark matter halos are more stable against bar instabilities. We find that the change of the bar fraction with M$_{\mathrm{h}}$ and M$_{\mathrm{*}}$ is stronger if we consider a relation with the form $f_{\mathrm{bar}}=f_{\mathrm{bar}}$(M$_{\mathrm{*}}^{\alpha}/$M$_{\mathrm{h}}$) with $\alpha=1.5$, and that the bar fraction is largely independent of other physical properties such as color and spin parameter when M$_{\mathrm{*}}^{3/2}/$M$_{\mathrm{h}}$ is fixed. With our sample of galaxies segregated into centrals and satellites, we also compare the fraction of barred galaxies in each group, finding a slightly higher bar fraction for satellites when compared with centrals at fixed stellar mass, but at fixed color this difference becomes very weak. This result, in agreement with previous studies, confirms that the bar fraction does not directly depend on the group/cluster environment, but the dependence exists through its dependence on internal morphology.

We present a study of the environment of barred galaxies using a volume-limited sample of over 30,000 galaxies drawn from the Sloan Digital Sky Survey. We use four different statistics to quantify the environment: the projected two-point cross-correlation function, the background-subtracted number count of neighbor galaxies, the overdensity of the local environment, and the membership of our galaxies to galaxy groups to segregate central and satellite systems. For barred galaxies as a whole, we find a very weak difference in all the quantities compared to unbarred galaxies of the control sample. When we split our sample into early- and late-type galaxies, we see a weak but significant trend for early-type galaxies with a bar to be more strongly clustered on scales from a few 100 kpc to 1 Mpc when compared to unbarred early-type galaxies. This indicates that the presence of a bar in early-type galaxies depends on the location within their host dark matter halos. This is confirmed by the group catalog in the sense that for early-types, the fraction of central galaxies is smaller if they have a bar. For late-type galaxies, we find fewer neighbors within $\sim$50 kpc around the barred galaxies when compared to unbarred galaxies form the control sample, suggesting that tidal forces from close companions suppress the formation/growth of bars. Finally, we find no obvious correlation between overdensity and the bars in our sample, showing that galactic bars are not obviously linked to the large-scale structure of the universe.