The Javalambre-Physics of the Accelerating Universe Astrophysical Survey (J-PAS) will soon start to scan thousands of square degrees of the northern extragalactic sky with a unique set of $56$ optical filters from a dedicated $2.55$m telescope, JST, at the Javalambre Astrophysical Observatory. Before the arrival of the final instrument (a 1.2 Gpixels, 4.2deg$^2$ field-of-view camera), the JST was equipped with an interim camera (JPAS-Pathfinder), composed of one CCD with a 0.3deg$^2$ field-of-view and resolution of 0.23 arcsec pixel$^{-1}$. To demonstrate the scientific potential of J-PAS, with the JPAS-Pathfinder camera we carried out a survey on the AEGIS field (along the Extended Groth Strip), dubbed miniJPAS. We observed a total of $\sim 1$ deg$^2$, with the $56$ J-PAS filters, which include $54$ narrow band (NB, $\rm{FWHM} \sim 145$Angstrom) and two broader filters extending to the UV and the near-infrared, complemented by the $u,g,r,i$ SDSS broad band (BB) filters. In this paper we present the miniJPAS data set, the details of the catalogues and data access, and illustrate the scientific potential of our multi-band data. The data surpass the target depths originally planned for J-PAS, reaching $\rm{mag}_{\rm {AB}}$ between $\sim 22$ and $23.5$ for the NB filters and up to $24$ for the BB filters ($5\sigma$ in a $3$~arcsec aperture). The miniJPAS primary catalogue contains more than $64,000$ sources extracted in the $r$ detection band with forced photometry in all other bands. We estimate the catalogue to be complete up to $r=23.6$ for point-like sources and up to $r=22.7$ for extended sources. Photometric redshifts reach subpercent precision for all sources up to $r=22.5$, and a precision of $\sim 0.3$% for about half of the sample. (Abridged)

We present a star/galaxy classification for the Southern Photometric Local Universe Survey (S-PLUS), based on a Machine Learning approach: the Random Forest algorithm. We train the algorithm using the S-PLUS optical photometry up to $r$=21, matched to SDSS/DR13, and morphological parameters. The metric of importance is defined as the relative decrease of the initial accuracy when all correlations related to a certain feature is vanished. In general, the broad photometric bands presented higher importance when compared to narrow ones. The influence of the morphological parameters has been evaluated training the RF with and without the inclusion of morphological parameters, presenting accuracy values of 95.0\% and 88.1\%, respectively. Particularly, the morphological parameter {\rm FWHM/PSF} performed the highest importance over all features to distinguish between stars and galaxies, indicating that it is crucial to classify objects into stars and galaxies. We investigate the misclassification of stars and galaxies in the broad-band colour-colour diagram $(g-r)$ versus $(r-i)$. The morphology can notably improve the classification of objects at regions in the diagram where the misclassification was relatively high. Consequently, it provides cleaner samples for statistical studies. The expected contamination rate of red galaxies as a function of the redshift is estimated, providing corrections for red galaxy samples. The classification of QSOs as extragalactic objects is slightly better using photometric-only case. An extragalactic point-source catalogue is provided using the classification without any morphology feature (only the SED information) with additional constraints on photometric redshifts and {\rm FWHM/PSF} values.

The Southern Photometric Local Universe Survey (S-PLUS) is imaging ~9300 deg^2 of the celestial sphere in twelve optical bands using a dedicated 0.8 m robotic telescope, the T80-South, at the Cerro Tololo Inter-American Observatory, Chile. The telescope is equipped with a 9.2k by 9.2k e2v detector with 10 um pixels, resulting in a field-of-view of 2 deg^2 with a plate scale of 0.55"/pixel. The survey consists of four main subfields, which include two non-contiguous fields at high Galactic latitudes (8000 deg^2 at |b| > 30 deg) and two areas of the Galactic plane and bulge (for an additional 1300 deg^2). S-PLUS uses the Javalambre 12-band magnitude system, which includes the 5 u, g, r, i, z broad-band filters and 7 narrow-band filters centered on prominent stellar spectral features: the Balmer jump/[OII], Ca H+K, H-delta, G-band, Mg b triplet, H-alpha, and the Ca triplet. S-PLUS delivers accurate photometric redshifts (delta_z/(1+z) = 0.02 or better) for galaxies with r < 20 AB mag and redshift < 0.5, thus producing a 3D map of the local Universe over a volume of more than 1 (Gpc/h)^3. The final S-PLUS catalogue will also enable the study of star formation and stellar populations in and around the Milky Way and nearby galaxies, as well as searches for quasars, variable sources, and low-metallicity stars. In this paper we introduce the main characteristics of the survey, illustrated with science verification data highlighting the unique capabilities of S-PLUS. We also present the first public data release of ~336 deg^2 of the Stripe-82 area, which is available at http://datalab.noao.edu/splus.

We show how the Radon transform (defined as a series of line integrals through an image at different orientations and offsets from the origin) can be used as a simple, non-parametric tool to characterize galaxy velocity fields, specifically their global kinematic position angles (PA_k) and any radial variation or asymmetry in PA_k. This method is fast and easily automated, making it particularly beneficial in an era where IFU and interferometric surveys are yielding samples of thousands of galaxies. We demonstrate the Radon transform by applying it to gas and stellar velocity fields from the first ~2800 galaxies of the SDSS-IV MaNGA IFU survey. We separately classify gas and stellar velocity fields into five categories based on the shape of their radial PA_k profiles. At least half of stellar velocity fields and two-thirds of gas velocity fields are found to show detectable deviations from uniform coplanar circular motion, although most of these variations are symmetric about the center of the galaxy. The behavior of gas and stellar velocity fields is largely independent, even when PA_k profiles for both components are measured over the same radii. We present evidence that one class of symmetric PA_k variations is likely associated with bars and/or oval distortions, while another class is more consistent with warped disks. This analysis sets the stage for more in-depth future studies which explore the origin of diverse kinematic behavior in the galaxy population.

We perform full spectrum fitting stellar population analysis and Jeans Anisotropic modelling (JAM) of the stellar kinematics for about 2000 early-type galaxies (ETGs) and spiral galaxies from the MaNGA DR14 sample. Galaxies with different morphologies are found to be located on a remarkably tight mass plane which is close to the prediction of the virial theorem, extending previous results for ETGs. By examining an inclined projection (`the mass-size' plane), we find that spiral and early-type galaxies occupy different regions on the plane, and their stellar population properties (i.e. age, metallicity and stellar mass-to-light ratio) vary systematically along roughly the direction of velocity dispersion, which is a proxy for the bulge fraction. Galaxies with higher velocity dispersions have typically older ages, larger stellar mass-to-light ratios and are more metal rich, which indicates that galaxies increase their bulge fractions as their stellar populations age and become enriched chemically. The age and stellar mass-to-light ratio gradients for low-mass galaxies in our sample tend to be positive ($\rm centre<outer$), while the gradients for most massive galaxies are negative. The metallicity gradients show a clear peak around velocity dispersion $\log_{10} \sigma_{\rm e}\approx 2.0$, which corresponds to the critical mass $\sim 3\times 10^{10}M_{\odot}$ of the break in the mass-size relation. Spiral galaxies with large mass and size have the steepest gradients, while the most massive ETGs, especially above the critical mass $M_{\rm crit}\ge 2\times 10^{11} M_{\odot}$, where slow rotator ETGs start dominating, have much flatter gradients. This may be due to differences in their evolution histories, e.g. mergers.

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 perform Jeans anisotropic modeling (JAM) on elliptical and spiral galaxies from the MaNGA DR13 sample. By comparing the stellar mass-to-light ratios estimated from stellar population synthesis (SPS) and from JAM, we find a similar systematic variation of the initial mass function (IMF) as in the earlier $\rm ATLAS^{3D}$ results. Early type galaxies (elliptical and lenticular) with lower velocity dispersions within one effective radius are consistent with a Chabrier-like IMF while galaxies with higher velocity dispersions are consistent with a more bottom heavy IMF such as the Salpeter IMF. Spiral galaxies have similar systematic IMF variations, but with slightly different slopes and larger scatters, due to the uncertainties caused by higher gas fractions and extinctions for these galaxies. Furthermore, we examine the effects of stellar mass-to-light ratio gradients on our JAM modeling, and find that the trends from our results becomes stronger after considering the gradients.

Diffuse Ionized Gas (DIG) is prevalent in star-forming galaxies. Using a sample of 365 nearly face-on star-forming galaxies observed by MaNGA, we demonstrate how DIG in star-forming galaxies impacts the measurements of emission line ratios, hence the interpretation of diagnostic diagrams and gas-phase metallicity measurements. At fixed metallicity, DIG-dominated low H\alpha\ surface brightness regions display enhanced [SII]/H\alpha, [NII]/H\alpha, [OII]/H\beta, and [OI]/H\alpha. The gradients in these line ratios are determined by metallicity gradients and H\alpha\ surface brightness. In line ratio diagnostic diagrams, contamination by DIG moves HII regions towards composite or LI(N)ER-like regions. A harder ionizing spectrum is needed to explain DIG line ratios. Leaky HII region models can only shift line ratios slightly relative to HII region models, and thus fail to explain the composite/LI(N)ER line ratios displayed by DIG. Our result favors ionization by evolved stars as a major ionization source for DIG with LI(N)ER-like emission. DIG can significantly bias the measurement of gas metallicity and metallicity gradients derived using strong-line methods. Metallicities derived using N2O2 are optimal because they exhibit the smallest bias and error. Using O3N2, R23, N2=[NII]/H\alpha, and N2S2H\alpha\ (Dopita et al. 2016) to derive metallicities introduces bias in the derived metallicity gradients as large as the gradient itself. The strong-line method of Blanc et al. (2015; IZI hereafter) cannot be applied to DIG to get an accurate metallicity because it currently contains only HII region models which fail to describe the DIG.

We present a study on the stellar age and metallicity distributions for 1105 galaxies using the STARLIGHT software on MaNGA integral field spectra. We derive age and metallicity gradients by fitting straight lines to the radial profiles, and explore their correlations with total stellar mass M*, NUV-r colour and environments, as identified by both the large scale structure (LSS) type and the local density. We find that the mean age and metallicity gradients are close to zero but slightly negative, which is consistent with the inside-out formation scenario. Within our sample, we find that both the age and metallicity gradients show weak or no correlation with either the LSS type or local density environment. In addition, we also study the environmental dependence of age and metallicity values at the effective radii. The age and metallicity values are highly correlated with M* and NUV-r and are also dependent on LSS type as well as local density. Low-mass galaxies tend to be younger and have lower metallicity in low-density environments while high-mass galaxies are less affected by environment.

With the availability of large integral-field unit (IFU) spectral surveys of nearby galaxies, there is now the potential to extract spectral information from across the bulges and discs of galaxies in a systematic way. This information can address questions such as how these components built up with time, how galaxies evolve and whether their evolution depends on other properties of the galaxy such as its mass or environment. We present BUDDI, a new approach to fit the two-dimensional light profiles of galaxies as a function of wavelength to extract the spectral properties of these galaxies' discs and bulges. The fitting is carried out using GalfitM, a modified form of Galfit which can fit multi-waveband images simultaneously. The benefit of this technique over traditional multi-waveband fits is that the stellar populations of each component can be constrained using knowledge over the whole image and spectrum available. The decomposition has been developed using commissioning data from the SDSS-IV Mapping Nearby Galaxies at APO (MaNGA) survey with redshifts z<0.14 and coverage of at least 1.5 effective radii for a spatial resolution of 2.5 arcsec FWHM and field of view of >22 arcsec, but can be applied to any IFU data of a nearby galaxy with similar or better spatial resolution and coverage. We present an overview of the fitting process, the results from our tests, and we finish with example stellar population analyses of early-type galaxies from the MaNGA survey to give an indication of the scientific potential of applying bulge-disc decomposition to IFU data.

Galaxies grow through both internal and external processes. In about 10% of nearby red galaxies with little star formation, gas and stars are counter-rotating, demonstrating the importance of external gas acquisition in these galaxies. However, systematic studies of such phenomena in blue, star-forming galaxies are rare, leaving uncertain the role of external gas acquisition in driving evolution of blue galaxies. Based on new measurements with integral field spectroscopy of a large representative galaxy sample, we find an appreciable fraction of counter-rotators among blue galaxies (9 out of 489 galaxies). The central regions of blue counter-rotators show younger stellar populations and more intense, ongoing star formation than their outer parts, indicating ongoing growth of the central regions. The result offers observational evidence that the acquisition of external gas in blue galaxies is possible; the interaction with pre-existing gas funnels the gas into nuclear regions (< 1 kpc) to form new stars.

We present the stellar surface mass density {\it vs.} gas metallicity ($\Sigma_*-Z$) relation for more than 500,000 spatially-resolved star-forming resolution elements (spaxels) from a sample of 653 disk galaxies included in the SDSS IV MaNGA survey. We find a tight relation between these local properties, with higher metallicities as the surface density increases. This relation extends over three orders of magnitude in the surface mass density and a factor of four in metallicity. We show that this local relationship can simultaneously reproduce two well-known properties of disk galaxies: their global mass-metallicity relationship {\it and} their radial metallicity gradients. We also find that the $\Sigma_* - Z$ relation is largely independent of the galaxy's total stellar mass and specific star-formation rate (sSFR), except at low stellar mass and high sSFR. These results suggest that in the present-day universe local properties play a key role in determining the gas-phase metallicity in typical disk galaxies.

We study the galaxy cosmological mass function (GCMF) in a semi-empirical relativistic approach using observational data provided by galaxy redshift surveys. Starting from the theory of Ribeiro & Stoeger (2003, arXiv:astro-ph/0304094) between the mass-to-light ratio, the selection function obtained from the luminosity function (LF) data and the luminosity density, the average luminosity $L$ and the average galactic mass $\mathcal{M}_g$ are computed in terms of the redshift. $\mathcal{M}_g$ is also alternatively estimated by a method that uses the galaxy stellar mass function (GSMF). Comparison of these two forms of deriving the average galactic mass allows us to infer a possible bias introduced by the selection criteria of the survey. We used the FORS Deep Field galaxy survey sample of 5558 galaxies in the redshift range $0.5 < z < 5.0$ and its LF Schechter parameters in the B-band, as well as this sample's stellar mass-to-light ratio and its GSMF data. Assuming ${\mathcal{M}_{g_0}} \approx 10^{11} \mathcal{M}_\odot$ as the local value of the average galactic mass, the LF approach results in $L_{B} \propto (1+z)^{(2.40 \pm 0.03)}$ and $\mathcal{M}_g \propto (1+z)^{(1.1\pm0.2)}$. However, using the GSMF results produces $\mathcal{M}_g \propto (1+z)^{(-0.58 \pm 0.22)}$. We chose the latter result as it is less biased. We then obtained the theoretical quantities of interest, such as the differential number counts, to calculate the GCMF, which can be fitted by a Schechter function. The derived GCMF follows theoretical predictions in which the less massive objects form first, being followed later by more massive ones. In the range $0.5 < z < 2.0$ the GCMF has a strong variation that can be interpreted as a higher rate of galaxy mergers or as a strong evolution in the star formation history of these galaxies.

The goal of this work is to compute the number density of far-IR selected galaxies in the comoving frame and along the past lightcone of observationally constrained LTB models and to compare those results with their standard model counterparts. We derived integral number densities and differential number densities using different cosmological distance definitions in the LTB models. Then, we computed selection functions and consistency functions for the luminosity functions in the combined fields of the Herschel/PACS evolutionary probe (PEP) survey in both standard and void cosmologies, from which we derived the observed values of the above-mentioned densities. We used the Kolmogorov-Smirnov statistics to study both the evolution of the consistency functions and its connection to the evolution of the comoving density of sources. Finally, we fitted the power-law behaviour of the densities along the observer's past lightcone. The analysis of the comoving number density shows that the increased flexibility of the LTB models is not enough to fit the observed redshift evolution of the number counts, if it is specialised to a recent best-fit giant void parametrisation. The results for the power-law fits of the densities along the observer's past lightcone show general agreement across both cosmological models studied here around a slope of -2.5 $\pm$ 0.1 for the integral number density on the luminosity-distance volumes. The differential number densities show much bigger slope discrepancies. We conclude that the differential number densities on the observer's past lightcone were still rendered dependent on the cosmological model by the flux limits of the PEP survey. In addition, we show that an intrinsic evolution of the sources must be assumed to fit the comoving number-density redshift evolution in the giant void parametrisation for the LTB models used in this work.

This paper studies the connection between the relativistic number density of galaxies down the past light cone in a Friedmann-Lemaitre-Robertson-Walker spacetime with non-vanishing cosmological constant and the galaxy luminosity function (LF) data. It extends the redshift range of previous results presented in Albani et al. (2007:astro-ph/0611032) where the galaxy distribution was studied out to z=1. Observational inhomogeneities were detected at this range. This research also searches for LF evolution in the context of the framework advanced by Ribeiro and Stoeger (2003:astro-ph/0304094), further developing the theory linking relativistic cosmology theory and LF data. Selection functions are obtained using the Schechter parameters and redshift parametrization of the galaxy luminosity functions obtained from an I-band selected dataset of the FORS Deep Field galaxy survey in the redshift range 0.5<z<5.0 for its blue bands and 0.75<z<3.0 for its red ones. Differential number counts, densities and other related observables are obtained, and then used with the calculated selection functions to study the empirical radial distribution of the galaxies in a fully relativistic framework. The redshift range of the dataset used in this work, which is up to five times larger than the one used in previous studies, shows an increased relevance of the relativistic effects of expansion when compared to the evolution of the LF at the higher redshifts. The results also agree with the preliminary ones presented in Albani et al. (2007:astro-ph/0611032), suggesting a power-law behavior of relativistic densities at high redshifts when they are defined in terms of the luminosity distance.

We report the discovery of a young massive stellar cluster and infrared nebula in the direction of the CS molecular cloud associated to the IRAS point source 16132-5039. The analysis of the mid-infrared images from the more accurate MSX catalog, reveled that there are two independent components associated with the IRAS source. The integral of the spectral energy distribution for these components, between 8.28 $\mu$m and 100 $\mu$m, gave lower limits for the bolometric luminosity of the embedded objects of $8.7 \times 10^4 L_\odot$ and $9 \times 10^3 L_\odot$, which corresponds to ZAMS O8 and B0.5 stars, respectively. The number of Lyman continuum photons expected from the stars that lie along the reddening line for early-type stars is about 1.7 $ \times$ 10$^{49}$ s$^{-1}$, enough to produce the detected flux densities at 5 GHz. The NIR spectrum of the nebula increases with frequency, implying that free-free emission cannot be the main source of the extended luminosity, from which we conclude that the observed emission must be mainly dust scattered light. A comparison of the cluster described in this paper with the young stellar cluster associated with the IRAS source 16177-5018, which is located at the same distance and direction, shows that the mean visual absorption of the newly discovered cluster is about 10 magnitudes smaller and it contains less massive stars, suggesting that it was formed from a less massive molecular cloud.

Wide-separation lensed QSOs measure the mass function and evolution of massive galaxy clusters, in a similar way to the cluster mass function deduced from X-ray-selected samples or statistical measurements of the Sunyaev-Zeldovich effect. We compute probabilities of strong lensing of QSOs by galaxy clusters in dark energy cosmologies using semianalytical modelling and explore the sensitivity of the method to various input parameters and assumptions. We highlight the importance of considering both the variation of halo properties with mass, redshift and cosmology and the effect of cosmic scatter in halo concentration. We then investigate the extent to which observational surveys for wide-separation lensed QSOs may be used to measure cosmological parameters such as the fractional matter density Omega_M, the rms linear density fluctuation in spheres of 8 Mpc/h, sigma_8, and the dark energy equation of state parameter w. We find that wide-separation lensed QSOs can measure sigma_8 and Omega_M in an equivalent manner to other methods such as cluster abundance studies and cosmic shear measurements. In assessing whether lensing statistics can distinguish between values of w, we conclude that at present the uncertainty in the calibration of sigma_8 in quintessence models dominates the conclusions reached. Nonetheless, lensing searches based on current QSO surveys such as the Two-degree Field and the Sloan Digital Sky Survey with 10^4-10^5 QSOs should detect systems with angular separations greater than 5'' and hence can provide an important test of the standard cosmological model that is complementary to measurements of cosmic microwave background anisotropies.

We report the possible discovery of multiple gravitationally-lensed images of QSOs with angular separations on arcminute scales. The QSOs were selected from the completed 2dF QSO survey as having redshifts and optical colours which are statistically consistent. In this paper we present higher-quality optical spectra of the candidates and discuss the likelihood of their genuinely being systems lensed by massive clusters of galaxies. From a comparison of the spectra it appears that up to six pairs of QSOs may be lensed multiple images, although the true number may be less than that and further observations should be undertaken to amass more evidence and to detect the lensing clusters. Two of the candidates may be associated with low redshift clusters of galaxies.