We investigate the environmental dependence of the local gas-phase metallicity in a sample of star-forming galaxies from the MaNGA survey. Satellite galaxies with stellar masses in the range $9<\log(M_{*}/M_{\odot})<10$ are found to be $\sim 0.05 \, \mathrm{dex}$ higher in metallicity than centrals of similar stellar mass. Within the low-mass satellite population, we find that the interstellar medium (ISM) metallicity depends most strongly on the stellar mass of the galaxy that is central to the halo, though there is no obvious difference in the metallicity gradients. At fixed total stellar mass, the satellites of high mass ($M_{*}>10^{10.5} \, \mathrm{M_{\odot}}$) centrals are $\sim 0.1 \, \mathrm{dex}$ more metal rich than satellites of low-mass ($M_{*} < 10^{10} \, \mathrm{M_{\odot}}$) centrals, controlling for local stellar mass surface density and gas fraction. Fitting a gas-regulator model to the spaxel data, we are able to account for variations in the local gas fraction, stellar mass surface density and local escape velocity-dependent outflows. We find that the best explanation for the metallicity differences is the variation in the average metallicity of accreted gas between different environments that depends on the stellar mass of the dominant galaxies in each halo. This is interpreted as evidence for the exchange of enriched gas between galaxies in dense environments that is predicted by recent simulations.

Galaxy interaction is considered a key driver of galaxy evolution and star formation (SF) history. In this paper, we present an empirical picture of the radial extent of interaction-triggered SF along the merger sequence. The samples under study are drawn from the integral field spectroscopy (IFS) survey SDSS-IV MaNGA, including 205 star-forming galaxies in pairs/mergers and ~1350 control galaxies. For each galaxy in pairs, the merger stage is identified according to its morphological signatures: incoming phase, at first pericenter passage, at apocenter, in merging phase, and in final coalescence. The effect of interactions is quantified by the global and spatially resolved SF rate relative to the SF rate of a control sample selected for each individual galaxy ($\Delta$logSFR and $\Delta$logsSFR(r), respectively). Analysis of the radial $\Delta$logsSFR(r) distributions shows that galaxy interactions have no significant impact on the $\Delta$logsSFR(r) during the incoming phase. Right after the first pericenter passage, the radial $\Delta$logsSFR(r) profile decreases steeply from enhanced to suppressed activity for increasing galactocentric radius. Later on, SF is enhanced on a broad spatial scale out to the maximum radius we explore (~6.7 kpc) and the enhancement is in general centrally peaked. The extended SF enhancement is also observed for systems at their apocenters and in the coalescence phase, suggesting that interaction-triggered SF is not restricted to the central region of a galaxy. Further explorations of a wide range in parameter space of merger configurations (e.g., mass ratio) are required to constrain the whole picture of interaction-triggered SF.

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.

Asymmetries in HI in galaxies are often caused by the interaction with close companions, making isolated galaxies an ideal framework to study secular evolution. The AMIGA project has demonstrated that isolated galaxies show the lowest level of asymmetry in their HI integrated profiles, yet some present significant asymmetries. CIG 96 (NGC 864) is a representative case reaching a 16% level. Our aim is to investigate the HI asymmetries of this spiral galaxy and what processes have triggered the star-forming regions observed in the XUV pseudoring. We performed deep optical observations at CAHA 1.23m, 2.2m and VST telescopes. We reach surface brightness (SB) limits of mu_2.2m = 27.5 mag arcsec-2 (Cous R) and mu_VST = 28.7mag arcsec-2 (r) that show the XUV pseudoring of the galaxy in detail. Additionally, a wavelet filtering of the HI data cube from our deep observations with E/VLA telescope allowed us to reach a column density of N_HI = 8.9x10^18 cm -2 (5sigma) (28"x28" beam), lower than in any isolated galaxy. We confirm that the HI extends farther than 4xr_25 in all directions. Furthermore, we detect for the first time two gaseous structures (10^6 Msol) in the outskirts. The g-r colour index image from 1.23m shows extremely blue colours in certain regions of the pseudoring where N_HI>8.5x10^20 cm-2 , whereas the rest show red colours. Galactic cirrus contaminate the field, setting an unavoidable detection limit at 28.5mag arcsec-2 (r). We detect no stellar link within 1degx1deg or gaseous link within 40'x40' between CIG 96 and any companion. The isolation criteria rule out interactions with other similar-sized galaxies for at least 2.7Gyr. Using existing stellar evolution models, the age of the pseudoring is estimated at 1Gyr or older. Undetected previously accreted companions and cold gas accretion remain as the main hypothesis to explain the optical pseudoring and HI features of CIG 96.

Mergers of galaxies are an important mode for galaxy evolution because they serve as an efficient trigger of powerful starbursts. However, observational studies of the molecular gas properties during their early stages are scarce. We present interferometric CO(2-1) maps of two luminous infrared galaxies (LIRGs), NGC 3110 and NGC 232, obtained with the Submillimeter Array (SMA) with ~ 1 kpc resolution. While NGC 3110 is a spiral galaxy interacting with a minor (14:1 stellar mass) companion, NGC 232 is interacting with a similarly sized object. We find that such interactions have likely induced in these galaxies enhancements in the molecular gas content and central concentrations, partly at the expense of atomic gas. The obtained molecular gas surface densities in their circumnuclear regions are $\Sigma_{\rm mol}~\gtrsim10^{2.5}$ M$_\odot$ pc$^{-2}$, higher than in non-interacting objects by an order of magnitude. Gas depletion times of ~ 0.5 - 1 Gyr are found for the different regions, lying in between non-interacting disk galaxies and the starburst sequence. In the case of NGC 3110, the spiral arms show on average 0.5 dex shorter depletion times than in the circumnuclear regions if we assume a similar H$_2$-CO conversion factor. We show that even in the early stages of the interaction with a minor companion, a starburst is formed along the circumnuclear region and spiral arms, where a large population of SSCs is found (~350), and at the same time a large central gas concentration is building up which might be the fuel for an active galactic nucleus. The main morphological properties of the NGC 3110 system are reproduced by our numerical simulations and allow us to estimate that the current epoch of the interaction is at ~ 150 Myrs after closest approach.

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 investigate the star formation history and the dust attenuation in the galaxy merger Mrk848. Thanks to the multiwavelength photometry from the ultraviolet (UV) to the infrared (IR), and MaNGA's integral field spectroscopy, we are able to study this merger in a detailed way. We divide the whole merger into the core and tail regions, and fit both the optical spectrum and the multi-band spectral energy distribution (SED) to models to obtain the star formation properties for each region respectively. We find that the color excess of stars in the galaxy $E(B-V)_s^{SED}$ measured with the multi-band SED fitting is consistent with that estimated both from the infrared excess (the ratio of IR to UV flux) and from the slope of the UV continuum. Furthermore, the reliability of the $E(B-V)_s^{SED}$ is examined with a set of mock SEDs, showing that the dust attenuation of the stars can be well constrained by the UV-to-IR broadband SED fitting. The dust attenuation obtained from optical continuum $E(B-V)_s^{spec}$ is only about half of $E(B-V)_s^{SED}$. The ratio of the $E(B-V)_s^{spec}$ to the $E(B-V)_g$ obtained from the Balmer decrement is consistent with the local value (around 0.5). The difference between the results from the UV-to-IR data and the optical data is consistent with the picture that younger stellar populations are attenuated by an extra dust component from the birth clouds compared to older stellar populations which are only attenuated by the diffuse dust. Both with the UV-to-IR SED fitting and the spectral fitting, we find that there is a starburst younger than 100~Myr in one of the two core regions, consistent with the scenario that the interaction-induced gas inflow can enhance the star formation in the center of galaxies.

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.

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 study the internal gradients of stellar population properties within $1.5\;R_{\rm e}$ for a representative sample of 721 galaxies with stellar masses ranging between $10^{9}\;M_{\odot}$ to $10^{11.5}\;M_{\odot}$ from the SDSS-IV MaNGA IFU survey. Through the use of our full spectral fitting code FIREFLY, we derive light and mass-weighted stellar population properties and their radial gradients, as well as full star formation and metal enrichment histories. We also quanfify the impact that different stellar population models and full spectral fitting routines have on the derived stellar population properties, and the radial gradient measurements. In our analysis, we find that age gradients tend to be shallow for both early-type and late-type galaxies. {\em Mass-weighted} age gradients of early-types are positive ($\sim 0.09\; {\rm dex}/R_{\rm e}$) pointing to "outside-in" progression of star formation, while late-type galaxies have negative {\em light-weighted} age gradients ($\sim -0.11\; {\rm dex}/R_{\rm e}$), suggesting an "inside-out" formation of discs. We detect negative metallicity gradients in both early and late-type galaxies, but these are significantly steeper in late-types, suggesting that radial dependence of chemical enrichment processes and the effect of gas inflow and metal transport are far more pronounced in discs. Metallicity gradients of both morphological classes correlate with galaxy mass, with negative metallicity gradients becoming steeper with increasing galaxy mass. The correlation with mass is stronger for late-type galaxies, with a slope of $d(\nabla [Z/H])/d(\log M)\sim -0.2\pm 0.05\;$, compared to $d(\nabla [Z/H])/d(\log M)\sim -0.05\pm 0.05\;$ for early-types. This result suggests that the merger history plays a relatively small role in shaping metallicity gradients of galaxies.

We study the internal radial gradients of stellar population properties within $1.5\;R_{\rm e}$ and analyse the impact of galaxy environment. We use a representative sample of 721 galaxies with masses ranging between $10^{9}\;M_{\odot}$ to $10^{11.5}\;M_{\odot}$ from the SDSS-IV survey MaNGA. We split this sample by morphology into early-type and late-type galaxies. Using the full spectral fitting code FIREFLY, we derive the light and mass-weighted stellar population properties age and metallicity, and calculate the gradients of these properties. We use three independent methods to quantify galaxy environment, namely the $N^{th}$ nearest neighbour, the tidal strength parameter $Q$ and distinguish between central and satellite galaxies. In our analysis, we find that early-type galaxies generally exhibit shallow light-weighted age gradients in agreement with the literature and mass-weighted median age gradients tend to be slightly positive. Late-type galaxies, instead, have negative light-weighted age gradients. We detect negative metallicity gradients in both early and late-type galaxies that correlate with galaxy mass, with the gradients being steeper and the correlation with mass being stronger in late-types. We find, however, that stellar population gradients, for both morphological classifications, have no significant correlation with galaxy environment for all three characterisations of environment. Our results suggest that galaxy mass is the main driver of stellar population gradients in both early and late-type galaxies, and any environmental dependence, if present at all, must be very subtle.

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 study the properties of 66 galaxies with kinematically misaligned gas and stars from MaNGA survey. The fraction of kinematically misaligned galaxies varies with galaxy physical parameters, i.e. M*, SFR and sSFR. According to their sSFR, we further classify these 66 galaxies into three categories, 10 star-forming, 26 "Green Valley" and 30 quiescent ones. The properties of different types of kinematically misaligned galaxies are different in that the star-forming ones have positive gradient in D4000 and higher gas-phase metallicity, while the green valley/quiescent ones have negative D4000 gradients and lower gas-phase metallicity on average. There is evidence that all types of the kinematically misaligned galaxies tend to live in more isolated environment. Based on all these observational results, we propose a scenario for the formation of star forming galaxies with kinematically misaligned gas and stars - the progenitor accretes misaligned gas from a gas-rich dwarf or cosmic web, the cancellation of angular momentum from gas-gas collisions between the pre-existing gas and the accreted gas largely accelerates gas inflow, leading to fast centrally-concentrated star-formation. The higher metallicity is due to enrichment from this star formation. For the kinematically misaligned green valley and quiescent galaxies, they might be formed through gas-poor progenitors accreting kinematically misaligned gas from satellites which are smaller in mass.

We report on CO (J = 2 - 1) mapping with the IRAM 30-m HERA receiver array of CGCG 97-079, an irregular galaxy in the merging galaxy cluster Abell 1367 (z = 0.022). We find that $\sim$ 80% of the detected CO (J = 2 - 1) is projected within a 16 arcsec$^{2}$ (6.5 kpc$^{2}$) region to the north and west of the optical/NIR centre, with the intensity maximum offset $\sim 10$ arcsec (4 kpc) NW of the optical/NIR centre and $\sim$ 7 arcsec (3 kpc) south-east of the HI intensity maximum. Evolutionary synthesis models indicate CGCG 97-079 experienced a burst of star formation $\sim$ 10$^8$ yr ago, most likely triggered by a tidal interaction with CGCG 97-073. For CGCG 97-079 we deduce an infall velocity to the cluster of $\sim$ 1000 km s$^{-1}$ and moderate ram pressure (P$_\mathrm{ram} \sim 10^{-11}$ dyn cm$^{-2}$). The observed offset in CGCG 97-079 of the highest density HI and CO (J = 2 - 1) from the stellar components has not previously been observed in galaxies currently undergoing ram pressure stripping, although previous detailed studies of gas morphology and kinematics during ram pressure stripping were restricted to significantly more massive galaxies with deeper gravitational potential wells. We conclude the observed cold gas density maxima offsets are most likely the result of ram pressure and/or the high-speed tidal interaction with CGCG 97-073. However ram pressure stripping is likely to be playing a major role in the perturbation of lower density gas.

HI kinematic asymmetries are common in late-type galaxies irrespective of environment, although the amplitudes are strikingly low in isolated galaxies. As part of our studies of the HI morphology and kinematics in isolated late-type galaxies we have chosen several very isolated galaxies from the AMIGA sample for HI mapping. Here we present GMRT 21-cm HI line mapping of CIG 340 which was selected because its integrated HI spectrum has a very symmetric profile, Aflux = 1.03. Optical images of the galaxy hinted at a warped disk in contrast to the symmetric integrated HI spectrum profile. Our aim is to determine the extent to which the optical asymmetry is reflected in the resolved HI morphology and kinematics. GMRT observations reveal significant HI morphological asymmetries in CIG 340 despite it's overall symmetric optical form and highly symmetric HI spectrum. The most notable HI features are: 1) a warp in the HI disk (with an optical counterpart), 2) the HI north/south flux ratio = 1.32 is much larger than expected from the integrated HI spectrum profile and 3) a ~ 45" (12 kpc) HI extension, containing ~ 6% of the detected HI mass on the northern side of the disk. We conclude that in isolated galaxies a highly symmetric HI spectrum can mask significant HI morphological asymmetries. The northern HI extension appears to be the result of a recent perturbation (10^8 yr), possibly by a satellite which is now disrupted or projected within the disk. This study provides an important step in our ongoing program to determine the predominant source of HI asymmetries in isolated galaxies. For CIG 340 the isolation from major companions, symmetric HI spectrum, optical morphology and interaction timescales have allowed us to narrow the possible causes the HI asymmetries and identify tests to further constrain the source of the asymmetries.

Multi-molecule observations towards an increasing variety of galaxies have been showing that the relative molecular abundances are affected by the type of activity. However, these studies are biased towards bright active galaxies, which are typically in interaction. We study the molecular composition of one of the most isolated galaxies in the local Universe where the physical and chemical properties of their molecular clouds have been determined by intrinsic mechanisms. We present 3 mm broad band observations of the galaxy CIG 638, extracted from the AMIGA sample of isolated galaxies. The emission of the J=1-0 transitions of CCH, HCN, HCO+, and HNC are detected. Integrated intensity ratios between these line are compared with similar observations from the literature towards active galaxies including starburst galaxies (SB), active galactic nuclei (AGN), luminous infrared galaxies (LIRG), and GMCs in M33. A significantly high ratio of CCH with respect to HCN, HCO+, and HNC is found towards CIG 638 when compared with all other galaxies where these species have been detected. This points to either an overabundance of CCH or to a relative lack of dense molecular gas as supported by the low HCN/CO ratio, or both. The data suggest that the CIG 638 is naturally a less perturbed galaxy where a lower fraction of dense molecular gas, as well as a more even distribution could explain the measured ratios. In this scenario the dense gas tracers would be naturally dimmer, while the UV enhanced CCH, would be overproduced in a less shielded medium.