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

Proportional electroluminescence (EL) in noble gases is used in two-phase detectors for dark matter search to record (in the gas phase) the ionization signal induced by particle scattering in the liquid phase. The "standard" EL mechanism is considered to be due to noble gas excimer emission in the vacuum ultraviolet (VUV). In addition to it, there are two alternative mechanisms, producing light in the visible and near infrared (NIR) range. The first is due to bremsstrahlung of electrons scattered on neutral atoms (so-called neutral bremsstrahlung, NBrS). The second, responsible for electron avalanche scintillations in the NIR at higher electric fields, is due to transitions between excited atomic states. In this work, we have for the first time demonstrated two alternative techniques of the optical readout of two-phase argon detectors, in the visible and NIR range, using a SiPM matrix and electroluminescence due to either neutral bremsstrahlung or avalanche scintillations effects. The amplitude yield and position resolution were measured for these readout techniques, which allowed to assess the detection threshold for electron and nuclear recoils in two-phase argon detectors for dark matter search. It should be emphasized that this is the first practical application of the NBrS effect in detection science.

Large liquid argon detectors offer one of the best avenues for the detection of galactic weakly interacting massive particles (WIMPs) via their scattering on atomic nuclei. The liquid argon target allows exquisite discrimination between nuclear and electron recoil signals via pulse-shape discrimination of the scintillation signals. Atmospheric argon (AAr), however, has a naturally occurring radioactive isotope, $^{39}$Ar, a $\beta$ emitter of cosmogenic origin. For large detectors, the atmospheric $^{39}$Ar activity poses pile-up concerns. The use of argon extracted from underground wells, deprived of $^{39}$Ar, is key to the physics potential of these experiments. The DarkSide-20k dark matter search experiment will operate a dual-phase time projection chamber with 50 tonnes of radio-pure underground argon (UAr), that was shown to be depleted of $^{39}$Ar with respect to AAr by a factor larger than 1400. Assessing the $^{39}$Ar content of the UAr during extraction is crucial for the success of DarkSide-20k, as well as for future experiments of the Global Argon Dark Matter Collaboration (GADMC). This will be carried out by the DArT in ArDM experiment, a small chamber made with extremely radio-pure materials that will be placed at the centre of the ArDM detector, in the Canfranc Underground Laboratory (LSC) in Spain. The ArDM LAr volume acts as an active veto for background radioactivity, mostly $\gamma$-rays from the ArDM detector materials and the surrounding rock. This article describes the DArT in ArDM project, including the chamber design and construction, and reviews the background required to achieve the expected performance of the detector.

This white paper summarizes the activities of the Brazilian community concerning dark matter physics and highlights the importance of financial support to Brazilian groups that are deeply involved in experimental endeavours. The flagships of the Brazilian dark matter program are the Cherenkov Telescope Array, DARKSIDE, SBN and LHC experiments, but we emphasize that smaller experiments such as DAMIC and CONNIE constitute important probes to dark sectors as well and should receive special attention. Small experimental projects showing the potential to probe new regions of parameter space of dark matter models are encouraged. On the theoretical and phenomenological side, some groups are devoted to astrophysical aspects such as the dark matter density profile while others explore the signature of dark matter models at colliders, direct and indirect detection experiments. In summary, the Brazilian dark matter community that was born not long ago has grown tremendously in the past years and now plays an important role in the hunt for a dark matter particle.

We report measurements of the charged daughter fraction of $^{218}$Po as a result of the $^{222}$Rn alpha decay, and the mobility of $^{218}$Po$^+$ ions, using radon-polonium coincidences from the $^{238}$U chain identified in 532 live-days of DarkSide-50 WIMP-search data. The fraction of $^{218}$Po that is charged is found to be 0.37$\pm$0.03 and the mobility of $^{218}$Po$^+$ is (8.6$\pm$0.1)$\times$10$^{-4}$$\frac{\text{cm}^2}{\text{Vs}}$.

We present maps of the ionized gas flux distributions, excitation, star-formation rate SFR, surface mass density $\Sigma_{H+}$, and obtain total values of SFR and ionized gas masses {\it M} for 62 Active Galactic Nuclei (AGN) observed with SDSS-IV MaNGA and compare them with those of a control sample of 112 non-active galaxies. The most luminous AGN -- with $L(\rm{[OIII]}\lambda 5007) \ge 3.8\times 10^{40}\,\mbox{erg}\,\mbox{s}^{-1}$, and those hosted by earlier-type galaxies are dominated by Seyfert excitation within 0.2 effective radius $R_e$ from the nucleus, surrounded by LINER excitation or transition regions, while the less luminous and hosted by later-type galaxies show equally frequent LINER and Seyfert excitation within $0.2\,R_e$. The extent $R$ of the region ionized by the AGN follows the relation $R\propto\,L(\rm{[OIII]})^{0.5}$ -- as in the case of the Broad-Line Region. The SFR distribution over the region ionized by hot stars is similar for AGN and controls, while the integrated SFR -- in the range $10^{-3}-10$\,M$_\odot$\,yr$^{-1}$ is also similar for the late-type sub-sample, but higher in the AGN for 75\% of the early-type sub-sample. We thus conclude that there is no signature of AGN quenching star formation in the body of the galaxy in our sample. We also find that 66\% of the AGN have higher ionized gas masses $M$ than the controls -- in the range 10$^5-3\times10^7$\,M$_\odot$ -- while 75\% of the AGN have higher $\Sigma_{H+}$ within $0.2\,R_e$ than the control galaxies.

Super-massive black holes (SMBH) are present at the center of most galaxies, with the related mass accretion processes giving origin to outflows in Active Galactic Nuclei (AGN). It has been presumed that only intense winds from luminous AGN were able to suppress star formation until the discovery of a new class of galaxies with no recent star formation and with the nucleus in a quiescent state showing kpc scale outflows. We used SDSS MaNGA and Gemini Integral Field Spectroscopy of the prototype Red Geyser Akira and found that the orientation of the outflow changes by about 50$^\circ$ from its nucleus to kpc scales. A possible interpretation is that the outflow is produced by a precessing accretion disk due to a misalignment between the orientation of the disk and the spin of the SMBH. The precession of the central source is also supported by a similar change in the orientation of the ionization pattern. Although similar behavior has commonly being reported for collimated relativistic jets, the precession of an AGN wide wind is reported here for the first time, implying on a larger work surface of the wind, which in turn increases the star formation suppression efficiency of the outflow.

We investigate the effects of Active Galactic Nuclei (AGN) on the gas kinematics of their host galaxies, using MaNGA data for a sample of 62 AGN hosts and 109 control galaxies (inactive galaxies). We compare orientation of the line of nodes (kinematic Position Angle - PA) measured from the gas and stellar velocity fields for the two samples. We found that AGN hosts and control galaxies display similar kinematic PA offsets between gas and stars. However, we note that AGN have larger fractional velocity dispersion $\sigma$ differences between gas and stars [$\sigma_{frac}=(\sigma_{\rm gas}-\sigma_{stars})/\sigma_{\rm stars}$] when compared to their controls, as obtained from the velocity dispersion values of the central (nuclear) pixel (2.5" diameter). The AGN have a median value of $\sigma_{\rm frac}$ of $<\sigma_{frac}>_{\rm AGN}=0.04$, while the the median value for the control galaxies is $<\sigma_{frac}>_{\rm CTR}=-0.23$. 75% of the AGN show $\sigma_{frac}>-0.13$, while 75% of the normal galaxies show $\sigma_{\rm frac}<-0.04$, thus we suggest that the parameter $\sigma_{\rm frac}$ can be used as an indicative of AGN activity. We find a correlation between the [OIII]$\lambda$5007 luminosity and $\sigma_{frac}$ for our sample. Our main conclusion is that the AGN already observed with MaNGA are not powerful enough to produce important outflows at galactic scales, but at 1-2 kpc scales, AGN feedback signatures are always present on their host galaxies.

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.

The DarkSide-50 direct-detection dark matter experiment is a dual-phase argon time projection chamber operating at Laboratori Nazionali del Gran Sasso. This paper reports on the blind analysis of a (16,660+-270) kg d exposure using a target of low-radioactivity argon extracted from underground sources. We find no events in the dark matter selection box and set a 90% C.L. upper limit on the dark matter-nucleon spin-independent cross section of 1.14E-44 cm^2 (3.78E-44 cm^2, 3.43E-43 cm^2) for a WIMP mass of 100 GeV/c^2 (1 TeV/c^2, 10 TeV/c^2).

We present new constraints on sub-GeV dark matter particles scattering off electrons in argon based on an analysis of ionization signal data from the DarkSide-50 detector.

We present the results of a search for dark matter WIMPs in the mass range below 20 GeV/c^2 using a target of low-radioactivity argon. The data were obtained using the DarkSide-50 apparatus at Laboratori Nazionali del Gran Sasso (LNGS). The analysis is based on the ionization signal, for which the DarkSide-50 time projection chamber is fully efficient at 0.1 keVee. The observed rate in the detector at 0.5 keVee is about 1.5 events/keVee/kg/day and is almost entirely accounted for by known background sources. We obtain a 90% C.L. exclusion limit above 1.8 GeV/c^2 for the spin-independent cross section of dark matter WIMPs on nucleons, extending the exclusion region for dark matter below previous limits in the range 1.8-6 GeV/c^2.

We present spatially resolved stellar population age maps, average radial profiles and gradients for the first 62 Active Galactic Nuclei (AGN) observed with SDSS-IV MaNGA to study the effects of the active nuclei on the star formation history of the host galaxies. These results, derived using the STARLIGHT code, are compared with a control sample of non-active galaxies matching the properties of the AGN hosts. We find that the fraction of young stellar populations (SP) in high-luminosity AGN is higher in the inner ($R \leq 0.5\,R_e$) regions when compared with the control sample; low-luminosity AGN, on the other hand, present very similar fractions of young stars to the control sample hosts for the entire studied range ($1\,R_e$). The fraction of intermediate age SP of the AGN hosts increases outwards, with a clear enhancement when compared with the control sample. The inner region of the galaxies (AGN and control galaxies) presents a dominant old SP, whose fraction decreases outwards. We also compare our results (differences between AGN and control galaxies) for the early and late-type hosts and find no significant differences. In summary, our results suggest that the most luminous AGN seems to have been triggered by a recent supply of gas that has also triggered recent star formation ($t\,\leq\,40\,Myrs$) in the central region.

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 report the characterization of the first $62$ MaNGA Active Galactic Nuclei (AGN) hosts in the Fifth Product Launch (MPL-5) and the definition of a control sample of non-active galaxies. This control sample - comprising two galaxies for each AGN - was selected in order to match the AGN hosts in terms of stellar mass, redshift, visual morphology and inclination. The stellar masses are in the range $9.4<\log(M/M_\odot)<11.5$, and most objects have redshifts $\leq 0.08$. The AGN sample is mostly comprised of low-luminosity AGN, with only 17 nuclei with $L([OIII]\lambda 5007)\geq 3.8\times 10^{40}$ erg s$^{-1}$ (that we call "strong AGN"). The stellar population of the control sample galaxies within the inner $1$-$3$ kpc is dominated by the old ($\sim$ $4$ - $13$ Gyr) age component, with a small contribution of intermediate age ($\sim 640$-$940$ Myr) and young stars ($\leq 40$ Myr) to the total light at $5700\AA$. While the weaker AGN show a similar age distribution to that of the control galaxies, the strong AGN show an increased contribution of younger stars and a decreased contribution of older stars. Examining the relationship between the AGN stellar population properties and $L([OIII])$, we find that with increasing $L([OIII])$, the AGN exhibit a decreasing contribution from the oldest ($>4$ Gyr) stellar population relative to control galaxies, but have an increasing contribution from the younger components with ages $\sim 40$ Myr. We also find a correlation of the mean age differences (AGN - control) with $L([OIII])$, in the sense that more luminous AGN are younger than the control objects, while the low-luminosity AGN seem to be older. These results support a connection between the growth of the galaxy bulge via formation of new stars and the growth of the Supermassive Black Hole via matter accretion in the AGN phase.

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.

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.

The recent detection of the cosmic microwave background polarimeter experiment BICEP2 of tensor fluctuations in the B-mode power spectrum basically excludes all plausible axion models where its decay constant is above $10^{13}$ GeV. Moreover, there are strong theoretical, astrophysical, and cosmological motivations for models involving, in addition to the axion, also axion-like particles (ALPs), with decay constants in the intermediate scale range, between $10^9$ GeV and $10^{13}$ GeV. Here, we present a general analysis of models with an axion and further ALPs and derive bounds on the relative size of the axion and ALP photon (and electron) coupling. We discuss what we can learn from measurements of the axion and ALP photon couplings about the fundamental parameters of the underlying ultraviolet completion of the theory. For the latter we consider extensions of the Standard Model in which the axion and the ALP(s) appear as pseudo Nambu-Goldstone bosons from the breaking of global chiral $U(1)$ (Peccei-Quinn (PQ)) symmetries, occuring accidentally as low energy remnants from exact discrete symmetries. In such models, the axion and the further ALP are protected from disastrous explicit symmetry breaking effects due to Planck-scale suppressed operators. The scenarios considered exploit heavy right handed neutrinos getting their mass via PQ symmetry breaking and thus explain the small mass of the active neutrinos via a seesaw relation between the electroweak and an intermediate PQ symmetry breaking scale. We show some models that can accommodate simultaneously an axion dark matter candidate, an ALP explaining the anomalous transparency of the universe for $\gamma$-rays, and an ALP explaining the recently reported 3.55 keV gamma line from galaxies and clusters of galaxies, if the respective decay constants are of intermediate scale.

Aims. Recently it has been reported that the intrinsic dispersion at constant magnitude of the structural relations from early-type galaxies is a useful tool to study the universality of these structural relations, that is to say, to study whether the structural relations depend on luminosity, wavelength, redshift and/or environment. In this work we study the intrinsic dispersion of the Faber-Jackson relation as function of the luminosity, mass and redshift. Methods. We use a sample of approximately 90 000 early-type galaxies from the Sloan Digital Sky Survey (SDSS-DR7) spanning a magnitude range of 7 $mag$ in both $g$ and $r$ filters. We calculate the intrinsic dispersion of the Faber-Jackson relation at approximately constant magnitude and compare this at different luminosities, masses and redshifts. Results. The main results are the following: i) The intrinsic dispersion of the Faber-Jackson relation depends on the luminosity, mass and redshift. ii) The distribution for brighter and more massive galaxies has smaller intrinsic dispersion than that for fainter and less massive galaxies. iii) The distribution of bright and massive galaxies at higher redshift has smaller intrinsic dispersion than those similar galaxies at low redshift. Conclusions. Comparisons of the results found in this work with recent studies from the literature make us conclude that the intrinsic dispersion of the Faber-Jackson relation could depend on the history of galaxies, in other words, the intrinsic dispersion could depend on the number and nature of transformation events that have affected the galaxies along their life times, such as collapse, accretion, interaction and merging.

We take a sample of early-type galaxies from the Sloan Digital Sky Survey (SDSS-DR7, $\sim$ 90 000 galaxies) spanning a range of approximately 7 $mag$ in both $g$ and $r$ filters and analyse the behaviour of the Faber-Jackson relation parameters as functions of the magnitude range. We calculate the parameters in two ways: i) We consider the faintest (brightest) galaxies in each sample and we progressively increase the width of the magnitude interval by inclusion of the brighter (fainter) galaxies (increasing-magnitude-intervals), and ii) we consider narrow-magnitude intervals of the same width ($\Delta M = 1.0$ $mag$) over the whole magnitude range available (narrow-magnitude-intervals). Our main results are that: i) in both increasing and narrow-magnitude-intervals the Faber-Jackson relation parameters change systematically, ii) non-parametric tests show that the fluctuations in the values of the slope of the Faber-Jackson relation are not products of chance variations. We conclude that the values of the Faber-Jackson relation parameters depend on the width of the magnitude range and the luminosity of galaxies within the magnitude range. This dependence is caused, to a great extent by the selection effects and because the geometrical shape of the distribution of galaxies on the $M - \log (\sigma_{0})$ plane depends on luminosity. We therefore emphasize that if the luminosity of galaxies or the width of the magnitude range or both are not taken into consideration when comparing the structural relations of galaxy samples for different wavelengths, environments, redshifts and luminosities, any differences found may be misinterpreted.