Search

We introduce the Data Analysis Pipeline (DAP) for the Sloan Digital Sky Survey V (SDSS-V) Local Volume Mapper (LVM) project, referred to as the LVM-DAP. We outline our methods for recovering both stellar and emission line components from the optical integral field spectroscopy, highlighting the developments and changes implemented to address specific challenges of the data set. The observations from the LVM project are unique because they cover a wide range of physical resolutions, from approximately 0.05 pc to 100 pc, depending on the distance to the targets. This, along with the varying number of stars sampled in each aperture (ranging from zero, just one of a few, to thousands), presents challenges in using previous spectral synthesis methods and interpreting the spectral fits. We provide a detailed explanation of how we model the stellar content and separate it from the ionized gas emission lines. To assess the accuracy of our results, we compare them with both idealized and more realistic simulations, highlighting the limitations of our methods. We find that the DAP robustly correct for stellar continuum features and recover emission line parameters (e.g. flux, equivalent width, systemtic velocity and velocity dispersion) with a precision and accuracy that fulfill the requirements of the primary goal of the analysis. In addition, the recovered stellar parameters are reliable for single stars, the recovery of integrated populations is less precise. We conclude with a description of the data products we provide, instructions for downloading and using our software, and a showcase illustrating the quality of the data and the analysis on a deep exposure taken on the Huygens region at the center of the Orion Nebula., Comment: 35 pages, 17 figures, 8 tables, accepted for publication in AJ

Understanding the complex ionization structure and chemical composition of \hii\ regions poses a significant challenge in astrophysics. The abundance discrepancy problem, characterized by inconsistencies between abundances derived from recombination lines (RLs) and collisionally excited lines (CELs), has long been a puzzle in the field. In this theoretical study, we present novel photoionization models that incorporate temperature, density, and chemical inhomogeneities within a single cloud to comprehensively address this discrepancy. By accounting for the intricate interplay between ionization, excitation, and chemistry, our models successfully reproduce both observed RLs and CELs with with an average difference between our models and the observations of 25% -- within uncertainties inherent in Galactic archival long-slit and new SDSS-V Local Volume Mapper observations. Through comparisons between generic inhomogeneous model predictions and observations, demonstrating the ability of our theoretical framework to analyze the abundance discrepancy problem within \hii\ regions. Our results highlight the importance of incorporating spatially resolved temperature, density, and chemical structures when interpreting the physical processes governing emission line spectra in these astrophysical environments.

The SDSS-IV/MaNGA Survey data provide an unprecedented opportunity to study the internal motions of galaxies and, in particular, represent the largest sample of barred galaxy kinematic maps obtained to date. We present results from Nirvana, our non-axisymmetric kinematic modeling code built with a physically-motivated Bayesian forward modeling approach, which decomposes MaNGA velocity fields into first- and second-order radial and tangential rotational modes in a generalized and minimally-supervised fashion. We use Nirvana to produce models and rotation curves for 1263 unique barred MaNGA galaxies and a matched unbarred control sample We present our modeling approach, tests of its efficacy, and validation against existing visual bar classifications. Nirvana finds elevated non-circular motions in galaxies identified as bars in imaging, and bar position angles that agree well with visual measurements. The Nirvana-MaNGA barred and control samples provide a new opportunity for studying the influence of non-axisymmetric internal disk kinematics in a large statistical sample., Comment: 21 pages, 17 figures

The fifth Sloan Digital Sky Survey (SDSS-V) Local Volume Mapper (LVM) is a wide-field integral field unit (IFU) survey that uses an array of four 160 mm fixed telescopes with siderostats to minimize the number of moving parts. Individual telescope observes the science field or calibration field independently and is synchronized with the science exposure. We developed the LVM Acquisition and Guiding Package (LVMAGP) optimized telescope control software program for LVM observations, which can simultaneously control four focusers, three K-mirrors, one fiber selector, four mounts (siderostats), and seven guide cameras. This software is built on a hierarchical architecture and the SDSS framework and provides three key sequences: autofocus, field acquisition, and autoguide. We designed and fabricated a proto-model siderostat to test the telescope pointing model and LVMAGP software. The mirrors of the proto-model were designed as an isogrid open-back type, which reduced the weight by 46% and enabled reaching thermal equilibrium quickly. Additionally, deflection due to bolting torque, self-gravity, and thermal deformation was simulated, and the maximum scatter of the pointing model induced by the tilt of optomechanics was predicted to be $4'.4$, which can be compensated for by the field acquisition sequence. We performed a real sky test of LVMAGP with the proto-model siderostat and obtained field acquisition and autoguide accuracies of $0''.38$ and $1''.5$, respectively. It met all requirements except for the autoguide specification, which will be resolved by more precise alignment among the hardware components at Las Campanas Observatory., Comment: 29 pages, 15 figures, in press at Journal of Astronomical Telescopes, Instruments, and Systems

We measure the age-velocity relationship from the lag between ionized gas and stellar tangential speeds in ~500 nearby disk galaxies from MaNGA in SDSS-IV. Selected galaxies are kinematically axisymmetric. Velocity lags are asymmetric drift, seen in the Milky Way's (MW) solar neighborhood and other Local Group galaxies; their amplitude correlates with stellar population age. The trend is qualitatively consistent in rate (d(sigma)/dt) with a simple power-law model where sigma is proportional to t^b that explains the dynamical phase-space stratification in the solar neighborhood. The model is generalized based on disk dynamical times to other radii and other galaxies. We find in-plane radial stratification parameters sigma_(0,r} (dispersion of the youngest populations) in the range of 10-40 km/s and 0.2

We present the Sloan Digital Sky Survey V (SDSS-V) Local Volume Mapper (LVM). The LVM is an integral-field spectroscopic survey of the Milky Way, Magellanic Clouds, and of a sample of local volume galaxies, connecting resolved pc-scale individual sources of feedback to kpc-scale ionized interstellar medium (ISM) properties. The 4-year survey covers the southern Milky Way disk at spatial resolutions of 0.05 to 1 pc, the Magellanic Clouds at 10 pc resolution, and nearby large galaxies at larger scales totaling $>4300$ square degrees of sky, and more than 55M spectra. It utilizes a new facility of alt-alt mounted siderostats feeding 16 cm refractive telescopes, lenslet-coupled fiber-optics, and spectrographs covering 3600-9800A at R ~ 4000. The ultra-wide field IFU has a diameter of 0.5 degrees with 1801 hexagonally packed fibers of 35.3 arcsec apertures. The siderostats allow for a completely stationary fiber system, avoiding instability of the line spread function seen in traditional fiber feeds. Scientifically, LVM resolves the regions where energy, momentum, and chemical elements are injected into the ISM at the scale of gas clouds, while simultaneously charting where energy is being dissipated (via cooling, shocks, turbulence, bulk flows, etc.) to global scales. This combined local and global view enables us to constrain physical processes regulating how stellar feedback operates and couples to galactic kinematics and disk-scale structures, such as the bar and spiral arms, as well as gas in- and out-flows., Comment: 29 pages, 12 figures, accepted for publication in The Astronomical Journal

The Hobby-Eberly Telescope (HET) VIRUS Parallel Survey (HETVIPS) is a blind spectroscopic program that sparsely covers approximately two-thirds of the celestial sphere and consists of roughly 252 million fiber spectra. The spectra were taken in parallel mode with the Visible Integral-field Replicable Unit Spectrograph (VIRUS) instrument when the HET was observing a primary target with other HET facility instruments. VIRUS can simultaneously obtain approximately 35,000 spectra covering 3470A to 5540A at a spectral resolution of ~800. Although the vast majority of these spectra cover blank sky, we used the Pan-STARRS1 Data Release 2 Stacked Catalog to identify objects encompassed in the HETVIPS pointings and extract their spectra. This paper presents the first HETVIPS data release, containing 493,012 flux-calibrated spectra obtained through 31 March 2023, as well as a description of the data processing technique. Each of the object spectra were classified, resulting in a catalog of 74,196 galaxies, 4,087 quasars, 259,396 stars, and 154,543 unknown sources., Comment: Accepted in ApJ March 18, 2024. Published in ApJ May 1, 2024

Gas rich galaxies are typically star-forming. We make use of HI-MaNGA, a program of HI follow-up for the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey of the Sloan Digital Sky Surveys to construct a sample of unusual neutral hydrogen (HI, 21cm) rich galaxies which have low Star Formation Rates (SFRs); using infra-red color from the Wide-field Infrared Survey Explorer (WISE) as a proxy for specific SFR. Out of a set of 1575 MaNGA galaxies with HI-MaNGA detections, we find 83 (5%) meet our selection criteria to be HI rich with low SFR. We construct two stellar mass-matched control samples: HI rich galaxies with typical SFR (High SF Control) and HI poor galaxies with low SFR (Low HI Control). We investigate the properties of each of these samples, comparing physical parameters such as ionization state maps, stellar and ionized gas velocity and dispersion, environment measures, metallicity, and morphology to search for the reasons why these unusual HI rich galaxies are not forming stars. We find evidence for recent external accretion of gas in some galaxies (via high counter-rotating fractions), along with some evidence for AGN feedback (from a high cLIER and/or red geyser fraction), and bar quenching (via an enhanced strong bar fraction). Some galaxies in the sample are consistent with simply having their HI in a high angular momentum, large radius, low density disc. We conclude that no single physical process can explain all HI rich, low SFR galaxies., Comment: 15 pages, in press MNRAS. v2 following corrections noticed in proofs

We present photometric and morphological analyses of nuclear star clusters (NSCs) -- very dense, massive star clusters present in the central regions of most galaxies -- in a sample of 33 massive disk galaxies within 20 Mpc, part of the "Composite Bulges Survey." We use data from the Hubble Space Telescope including optical (F475W and F814W) and near-IR (F160W) images from the Wide Field Camera 3. We fit the images in 2D to take into account the full complexity of the inner regions of these galaxies (including the contributions of nuclear disks and bars), isolating the nuclear star cluster and bulge components. We derive NSC radii and magnitudes in all 3 bands, which we then use to estimate NSC masses. Our sample significantly expands the sample of massive late-type galaxies with measured NSC properties. We clearly identify nuclear star clusters in nearly 80% of our galaxies, putting a lower limit on the nucleation fraction in these galaxies that is higher than previous estimates. We find that the NSCs in our massive disk galaxies are consistent with previous NSC mass-NSC radius and Galaxy Mass-NSC Mass relations. However, we also find a large spread in NSC masses, with a handful of galaxies hosting very low-mass, compact clusters. Our NSCs are aligned in PA with their host galaxy disks but are less flattened. They show no correlations with bar or bulge properties. Finally, we find the ratio of NSC to BH mass in our massive disk galaxy sample spans a factor of $\sim$300., Comment: Accepted to The Astrophysical Journal

Post-starburst galaxies (PSBs) are transition galaxies showing evidence of recent rapid star formation quenching. To understand the role of galaxy mergers in triggering quenching, we investigate the incidence of PSBs and resolved PSB properties in post-merger galaxies using both SDSS single-fiber spectra and MaNGA resolved IFU spectra. We find post-mergers have a PSB excess of 10 - 20 times that relative to their control galaxies using single-fiber PSB diagnostics. A similar excess of ~ 19 times is also found in the fraction of central (C)PSBs and ring-like (R)PSBs in post-mergers using the resolved PSB diagnostic. However, 60% of the CPSBs + RPSBs in both post-mergers and control galaxies are missed by the single-fiber data. By visually inspecting the resolved PSB distribution, we find that the fraction of outside-in quenching is 7 times higher than inside-out quenching in PSBs in post-mergers while PSBs in control galaxies do not show large differences in these quenching directions. In addition, we find a marginal deficit of HI gas in PSBs relative to non-PSBs in post-mergers using the MaNGA-HI data. The excesses of PSBs in post-mergers suggest that mergers play an important role in triggering quenching. Resolved IFU spectra are important to recover the PSBs missed by single-fiber spectra. The excess of outside-in quenching relative to inside-out quenching in post-mergers suggests that AGN are not the dominant quenching mechanism in these galaxies, but that processes from the disk (gas inflows/consumption and stellar feedback) play a more important role., Comment: Accepted in MNRAS on May 12 2023, 19 pages, 15 figures, 4 tables

We combine an unprecedented MaNGA sample of over 3,000 passive galaxies in the stellar mass range 10^{9}-10^{12} Msun with the Sloan Digital Sky Survey group catalog by Tinker to quantify how central and satellite formation, quantified by radial profiles in stellar age, [Fe/H], and [Mg/Fe], depends on the stellar mass of the galaxy (M*) and the mass of the host halo (Mh). After controlling for M* and Mh, the stacked spectra of centrals and satellites beyond the effective radius (r_e) show small, yet significant differences in multiple spectral features at the 1% level. According to spectral fitting with the code alf, a primary driver of these differences appears to be [Mg/Fe] variations, suggesting that stellar populations in the outskirts of satellites formed more rapidly than the outer populations of centrals. To probe the physical mechanisms that may be responsible for this signal, we examined how satellite stellar populations depend on Mh. We find that satellites in high-Mh halos show older stellar ages, lower [Fe/H], and higher [Mg/Fe] compared to satellites in low-Mh halos, especially for M*=10^{9.5}-10^{10.5} Msun. These signals lend support to environmentally driven processes that quench satellite galaxies, although variations in the merger histories of central and satellite galaxies also emerge as a viable explanation., Comment: Accepted for publication in ApJ. 19 pages and 6 figures

We study the bar pattern speeds and corotation radii of 225 barred galaxies, using IFU data from MaNGA and the Tremaine-Weinberg method. Our sample, which is divided between strongly and weakly barred galaxies identified via Galaxy Zoo, is the largest that this method has been applied to. We find lower pattern speeds for strongly barred galaxies than for weakly barred galaxies. As simulations show that the pattern speed decreases as the bar exchanges angular momentum with its host, these results suggest that strong bars are more evolved than weak bars. Interestingly, the corotation radius is not different between weakly and strongly barred galaxies, despite being proportional to bar length. We also find that the corotation radius is significantly different between quenching and star forming galaxies. Additionally, we find that strongly barred galaxies have significantly lower values for R, the ratio between the corotation radius and the bar radius, than weakly barred galaxies, despite a big overlap in both distributions. This ratio classifies bars into ultrafast bars (R < 1.0; 11% of our sample), fast bars (1.0 < R < 1.4; 27%) and slow bars (R > 1.4; 62%). Simulations show that R is correlated with the bar formation mechanism, so our results suggest that strong bars are more likely to be formed by different mechanisms than weak bars. Finally, we find a lower fraction of ultrafast bars than most other studies, which decreases the recently claimed tension with {\Lambda}CDM. However, the median value of R is still lower than what is predicted by simulations., Comment: 20 pages, 16 figures

We compare an observed Baryonic Tully-Fisher Relation (BTFR) from the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) and HI-MaNGA surveys to a simulated BTFR from the cosmological magnetohydrodynamical simulation IllustrisTNG. To do so, we calibrate the BTFR of the local universe using 377 galaxies from the MaNGA and HI-MaNGA surveys, and perform mock 21 cm observations of matching galaxies from IllustrisTNG. The mock observations are used to ensure that the comparison with the observed galaxies is fair since it has identical measurement algorithms, observational limitations, biases and uncertainties. For comparison, we also calculate the BTFR for the simulation without mock observations, and demonstrate how mock observations are necessary to fairly and consistently compare between observational and theoretical data. We report a MaNGA BTFR of log$_{10} (M_{ \rm Bary}/M_\odot)= (2.97 \pm 0.18)$ log$_{10} V_{ \rm Rot} + (4.04 \pm 0.41)\,\log_{10}{M_{\odot}}$ and an IllustrisTNG BTFR of log$_{10} (M_{ \rm Bary}/M_\odot) = (2.94 \pm 0.23$) log$_{10} V_{ \rm Rot} + (4.15 \pm 0.44)\,\log_{10}{M_{\odot}}$. Thus, MaNGA and IllustrisTNG produce BTFRs that agree within uncertainties, demonstrating that IllustrisTNG has created a galaxy population that obeys the observed relationship between mass and rotation velocity in the observed universe., Comment: 14 pages, 7 figures

We investigate the role of galaxy mergers in triggering AGN in the nearby Universe. Our analysis is based on a sample of 79 post-merger remnant galaxies with deep X-ray observations from Chandra/XMM-Newton capable of detecting a low-luminosity AGN of > 10^40.5 erg s^-1. This sample is derived from a visually classified, volume-limited sample of 807 post-mergers identified in the Sloan Digital Sky Survey Data Release 14 with log M*/M_sun > 10.5 and 0.02 < z < 0.06. We find the X-ray AGN fraction in this sample is 55.7% +\- 5.6% compared to 23.6% +\- 2.8% for a mass and redshift matched non-interacting control sample. The multi-wavelength AGN fraction (identified as an AGN in one of X-ray, IR, radio or optical diagnostics) for post-mergers is 76.6% +\- 4.8% compared to 39.1% +\- 3.2% for controls. Thus post-mergers exhibit a high overall AGN fraction with an excess between 2 - 4 depending on the AGN diagnostics used. In addition, we find most optical, IR, and radio AGN are also identified as X-ray AGN while a large fraction of X-ray AGN are not identified in any other diagnostic. This highlights the importance of deep X-ray imaging to identify AGN. We find the X-ray AGN fraction of post-mergers is independent of stellar mass above log M*/M_sun > 10.5 unlike the trend seen in control galaxies. Overall, our results show that post-merger galaxies are a good tracer of the merger-AGN connection and strongly support the theoretical expectations that mergers trigger AGN., Comment: Accepted for publication in ApJ, 27 pages, 16 figures, 3 tables

We present the first publicly released catalog of sources obtained from the Hobby-Eberly Telescope Dark Energy Experiment (HETDEX). HETDEX is an integral field spectroscopic survey designed to measure the Hubble expansion parameter and angular diameter distance at 1.88

In this paper, we place the Milky Way (MW) in the context of similar-looking galaxies in terms of their star-formation and chemical evolution histories. We select a sample of 138 Milky-Way analogues (MWAs) from the SDSS-IV/MaNGA survey based on their masses, Hubble types, and bulge-to-total ratios. To compare their chemical properties to the detailed spatially-resolved information available for the MW, we use a semi-analytic spectral fitting approach, which fits a self-consistent chemical-evolution and star-formation model directly to the MaNGA spectra. We model the galaxies' inner and outer regions assuming that some of the material lost in stellar winds falls inwards. We also incorporate chemical enrichment from type II and Ia supernovae to follow the alpha-element abundance at different metallicities and locations. We find some MWAs where the stellar properties closely reproduce the distribution of age, metallicity, and alpha enhancement at both small and large radii in the MW. In these systems, the match is driven by the longer timescale for star formation in the outer parts, and the inflow of enriched material to the central parts. However, other MWAs have very different histories. These divide into two categories: self-similar galaxies where the inner and outer parts evolve identically; and centrally-quenched galaxies where there is very little evidence of late-time central star formation driven by material accreted from the outer regions. We find that, although selected to be comparable, there are subtle morphological differences between galaxies in these different classes, and that the centrally-quenched galaxies formed their stars systematically earlier., Comment: 16 pages, 9 figures, MNRAS accepted version

We present a comparative study of active galactic nuclei (AGN) between galaxy pairs and isolated galaxies with the final data release of the MaNGA integral field spectroscopic survey. We build a sample of 391 kinematic galaxy pairs within the footprint of the survey and select AGN using the survey's spectra. We use the comoving volume densities of the AGN samples to quantify the effects that tidal interactions have on the triggering of nuclear accretion. Our hypothesis is that the pair sample contains AGN that are triggered by not only stochastic accretion but also tidally induced accretion and correlated accretion. With the level of stochastically triggered AGN fixed by the control sample, we model the strength of tidally induced accretion and correlated accretion as a function of projected separation (rp) and compare the model expectations with the observed volume densities of dual AGN and offset AGN (single AGN in a pair). At rp ~ 10 kpc, we find that tidal interactions induce ~30% more AGN than stochastic fueling and cause ~12% of the offset AGN to become dual AGN because of correlations. The strength of both these effects decreases with increasing rp. We also find that the OIII luminosities of the AGN in galaxy pairs are consistent with those found in isolated galaxies, likely because stochastically fed AGN dominate even among close pairs. Our results illustrates that while we can detect tidally induced effects statistically, it is challenging to separate tidally induced AGN and stochastically triggered AGN in interacting galaxies., Comment: To be published in ApJ. The pair sample is based on the catalog presented in the following work, https://iopscience.iop.org/article/10.3847/2515-5172./ac9282

We present a stellar parameter catalog built to accompany the MaStar Stellar Library, which is a comprehensive collection of empirical, medium-resolution stellar spectra. We constructed this parameter catalog by using a multicomponent $\chi^{2}$ fitting approach to match MaStar spectra to models generated by interpolating the ATLAS9-based BOSZ model spectra. The total $\chi^{2}$ for a given model is defined as the sum of components constructed to characterize narrow-band features of observed spectra (e.g., absorption lines) and the broadband continuum shape separately. Extinction and systematics due to flux calibration were taken into account in the fitting. The $\chi^{2}$ distribution for a given region of model space was sampled using a Markov Chain Monte Carlo (MCMC) algorithm, the data from which were then used to extract atmospheric parameter estimates ($T_{eff}$, log g, [Fe/H], and [$\alpha$/Fe]), their corresponding uncertainties, and direct extinction measurements. Two methods were used to extract parameters and uncertainties: one that accepts the MCMC's prescribed minimum-$\chi^{2}$ result, and one that uses Bayesian inference to compute a likelihood-weighted mean from the $\chi^{2}$ distribution sampled by the MCMC. Results were evaluated for internal consistency using repeat observations where available and by comparing them with external data sets (e.g., APOGEE-2 and Gaia DR2). Our spectral-fitting exercise reveals possible deficiencies in current theoretical model spectra, illustrating the potential power of MaStar spectra for helping to improve the models. This paper represents an update to the parameters that were originally presented with SDSS-IV DR17. The MaStar parameter catalog containing our BestFit results is available on the SDSS-IV DR17 website as part of version 2 of the MaStar stellar parameter value-added catalog., Comment: 22 pages, 21 figures. A&A in press

Environment has long been known to have significant impact on the evolution of galaxies, but here we seek to quantify the subtler differences that might be found in disk galaxies, depending on whether they are isolated, the most massive galaxy in a group (centrals), or a lesser member (satellites). The MaNGA survey allows us to define a large mass-matched sample of 574 galaxies with high-quality integrated spectra in each category. Initial examination of their spectral indices indicates significant differences, particularly in low-mass galaxies. Semi-analytic spectral fitting of a full chemical evolution model to these spectra confirms these differences, with low-mass satellites having a shorter period of star formation and chemical enrichment typical of a closed box, while central galaxies have more extended histories, with evidence of on-going gas accretion over their lifetimes. The derived parameters for gas infall timescale and wind strength suggest that low-mass satellite galaxies have their hot halos of gas effectively removed, while central galaxies retain a larger fraction of gas than isolated galaxies due to the deeper group potential well in which they sit. S0 galaxies form a distinct subset within the sample, particularly at higher masses, but do not bias the inferred lower-mass environmental impact significantly. The consistent picture that emerges underlines the wealth of archaeological information that can be extracted from high-quality spectral data using techniques like semi-analytic spectral fitting., Comment: 9 pages, 3 figures, MNRAS accepted

A fundamental assumption adopted in nearly every extragalactic emission-line study is that the attenuation of different emission lines can be described by a single attenuation curve. Here we show this assumption fails in many cases with important implications for derived results. We developed a new method to measure the differential nebular attenuation among three kinds of transitions: the Balmer lines of hydrogen, high-ionization transitions, and low-ionization transitions. This method bins the observed data in a multidimensional space spanned by attenuation-insensitive line ratios. Within each small bin, the variations in line ratios are mainly driven by the variations in the nebular attenuation. This allows us to measure the nebular attenuation using both forbidden lines and Balmer lines. We applied this method to a sample of 2.4 million star-forming spaxels from SDSS-IV MaNGA. We found that the attenuation of high ionization lines and Balmer lines can be well described by a single Fitzpatrick (1999) extinction curve with $R_V=3.1$. However, no single attenuation curve can simultaneously account for all three transitions. This strongly suggests that different lines have different effective attenuations, likely because spectroscopy at kiloparsec resolutions mixes multiple regions with different intrinsic line ratios and different levels of attenuation. As a result, the assumption that different lines follow the same attenuation curve breaks down. Using a single attenuation curve determined by Balmer lines to correct attenuation-sensitive forbidden line ratios could bias the nebular parameters derived by 0.06--0.25 dex at $A_V = 1$, depending on the details of the dust attenuation model. Observations of a statistically large sample of H II regions with high spatial resolutions and large spectral coverage are vital for improved modeling and deriving accurate corrections for this effect., Comment: 28 pages, 17 figures, published in A&A

We consider the largest sample of 561 edge-on galaxies observed with integral field units by the MaNGA survey and find 300 galaxies where the ionised gas shows a negative vertical gradient (lag) in its rotational speed. We introduce the stop altitude as the distance to the galactic midplane at which the gas rotation should stop in the linear approximation. We find correlations between the lags, stop altitude and galactic mass, stellar velocity dispersion and overall Sersic index. We do not find any correlation of the lags or stop altitude with the star formation activity in the galaxies. We conclude that low mass galaxies (log(M*/Mo) < 10) with low Sersic index and with low stellar velocity dispersion posses a wider "zone of influence" in the extragalactic gas surrounding them with respect to higher mass galaxies that have a significant spherical component. We estimated the trend of the vertical rotational gradient with radius and find it flat for most of the galaxies in our sample. A small subsample of galaxies with negative radial gradients of lag has an enhanced fraction of objects with aged low surface brightness structures around them (e.g. faint shells), which indicates that noticeable accretion events in the past affected the extraplanar gas kinematics and might have contributed to negative radial lag gradients. We conclude that an isotropic accretion of gas from the circumgalactic medium plays a significant role in the formation of rotation velocity lags., Comment: 14 pages, 13 figures. Accepted by MNRAS

The Local Volume Mapper (LVM) project is one of three surveys that form the Sloan Digital Sky Survey V. It will map the interstellar gas emission in a large fraction of the southern sky using wide-field integral field spectroscopy. Four 16-cm telescopes in siderostat configuration feed the integral field units (IFUs). A reliable acquisition and guiding (A&G) strategy will help ensure that we meet our science goals. Each of the telescopes hosts commercial CMOS cameras used for A&G. In this work, we present our validation of the camera performance. Our tests show that the cameras have a readout noise of around 5.6e- and a dark current of 21e-/s, when operated at the ideal gain setting and at an ambient temperature of 20{\deg}C. To ensure their performance at a high-altitude observing site, such as the Las Campanas Observatory, we studied the thermal behaviour of the cameras at different ambient pressures and with different passive cooling solutions. Using the measured properties, we calculated the brightness limit for guiding exposures. With a 5 s exposure time, we reach a depth of around 16.5 Gaia gmag with a signal-to-noise ratio (SNR)>5. Using Gaia Early Data Release 3, we verified that there are sufficient guide stars for each of the around 25000 survey pointings. For accurate acquisition, we also need to know the focal plane geometry. We present an approach that combines on-chip astrometry and using a point source microscope to measure the relative positions of the IFU lenslets and the individual CMOS pixels to around 2 $\mu$m accuracy., Comment: Submitted to SPIE Astronomical Telescopes + Instrumentation 2022, Ground-based and Airborne Instrumentation for Astronomy IX, Paper 12184-256

The recent discovery of the extremely lensed Earendel object at $z=6.2$ is remarkable in that it is likely a single star or stellar multiple, observed within the first billion years of cosmic history. Depending on its mass, which is still uncertain but will soon be more tightly constrained with the James Webb Space Telescope, the Earendel star might even be a member of the first generation of stars, the so-called Population~III (Pop~III). By combining results from detailed cosmological simulations of the assembly of the first galaxies, including the enrichment of the pristine gas with heavy chemical elements, with assumptions on key stellar parameters, we quantify the probability that Earendel has indeed a Pop~III origin. We find that this probability is non-negligible throughout the mass range inferred for Earendel, specifically ranging from a few percent at the lower-mass end to near unity for some Pop~III initial mass function (IMF) models towards the high-mass end of the allowed range. For models that extend the metal-enriched IMF to $500$\,M$_\odot$, the likelihood of Earendel being a Pop~III star stays at the few to ten percent level. We discuss the implications of such a discovery for the overall endeavor to probe the hitherto so elusive first stars in the universe., Comment: 5 pages, 2 figures, resubmitted to ApJL

We spectroscopically detected candidate emission-lines of 8 likely, 17 probable, and 69 possible strong galaxy-galaxy gravitational lens candidates found within the spectra of ~10,000 galaxy targets contained within the completed Mapping of Nearby Galaxies at Apache Point Observatory (MaNGA) survey. This search is based upon the methodology of the Spectroscopic Identification of Lensing Objects (SILO) project, which extends the spectroscopic detection methods of the BOSS Emission-Line Lensing Survey (BELLS) and the Sloan Lens ACS Survey (SLACS). We scanned the co-added residuals that we constructed from stacks of foreground subtracted row-stacked-spectra (RSS) so a sigma-clipping method can be used to reject cosmic-rays and other forms of transients that impact only a small fraction of the combined exposures. We also constructed narrow-band images from the signal-to-noise of the co-added residuals to observe signs of lensed source images. We also use several methods to compute the probable strong lensing regime for each candidate lens to determine which candidate background galaxies may reside sufficiently near the galaxy centre for strong lensing to occur. We present the spectroscopic redshifts within a value-added catalogue (VAC) for data release 17 (DR17) of SDSS-IV. We also present the lens candidates, spectroscopic data, and narrow-band images within a VAC for DR17. High resolution follow-up imaging of these lens candidates are expected to yield a sample of confirmed grade-A lenses with sufficient angular size to probe possible discrepancies between the mass derived from a best-fitting lens model, and the dynamical mass derived from the observed stellar velocities., Comment: Accepted for publication in MNRAS, June 20, 2022. In press. 30 pages, 8 figures, 3 tables, 2 appendices with 15 tables

We analyze spatially resolved and co-added SDSS-IV MaNGA spectra with signal-to-noise ~100 from 2200 passive central galaxies (z~0.05) to understand how central galaxy assembly depends on stellar mass (M*) and halo mass (Mh). We control for systematic errors in Mh by employing a new group catalog from Tinker (2020a,b) and the widely-used Yang et al. (2007) catalog. At fixed M*, the strength of several stellar absorption features varies systematically with Mh. Completely model-free, this is one of the first indications that the stellar populations of centrals with identical M* are affected by the properties of their host halos. To interpret these variations, we applied full spectral fitting with the code alf. At fixed M*, centrals in more massive halos are older, show lower [Fe/H], and have higher [Mg/Fe] with 3.5 sigma confidence. We conclude that halos not only dictate how much M* galaxies assemble, but also modulate their chemical enrichment histories. Turning to our analysis at fixed Mh, high-M* centrals are older, show lower [Fe/H], and have higher [Mg/Fe] for Mh>10^{12}Msun/h with confidence > 4 sigma. While massive passive galaxies are thought to form early and rapidly, our results are among the first to distinguish these trends at fixed Mh. They suggest that high-M* centrals experienced unique early formation histories, either through enhanced collapse and gas fueling, or because their halos were early-forming and highly concentrated, a possible signal of galaxy assembly bias., Comment: Accepted for publication in ApJ. 28 pages and 12 figures

We introduce a novel non-parametric method to find solar neighborhood analogs (SNAs) in extragalactic IFS surveys. The main ansatz is that the physical properties of the solar neighborhood (SN) should be encoded in its optical stellar spectrum. We assume that our best estimate of such spectrum is the one extracted from the analysis performed by the Code for Stellar properties Heuristic Assignment (\textsc{CoSHA}) from the MaStar stellar library. It follows that finding SNAs in other galaxies consist in matching, in a $\chi^2$ sense, the SN reference spectrum across the optical extent of the observed galaxies. We apply this procedure to a selection of CALIFA galaxies, by requiring a close to face-on projection, relative isolation and non-AGN. We explore how the local and global properties of the SNAs (stellar age, metallicity, dust extinction, mass-to-light ratio, stellar surface mass and star-formation densities and galactocentric distance) and their corresponding host galaxies (morphological type, total stellar mass, star-formation rate, effective radius) compare with those of the SN and the Milky Way (MW). We find that SNAs are located preferentially in S(B)a~--~S(B)c galaxies, in a ring-like structure, which radii seems to scale with the galaxy size. Despite the known sources of systematics and errors, {most} properties present a considerable agreement with the literature on the SN. We conclude that the solar neighborhood is relatively common in our sample of SNAs. Our results warrants a systematic exploration of correlations among the physical properties of the SNAs and their host galaxies. We reckon that our method should inform current models of the galactic habitable zone in our MW and other galaxies., Comment: 5 pages, 2 figures, Accepted for publication in A&A

We present, for the first time, the relationship between local stellar mass surface density, $\mathrm{\Sigma_{*}}$, and N/O derived from SDSS-IV MaNGA data, using a sample of $792765$ high signal-to-noise ratio star-forming spaxels. Using a combination of phenomenological modelling and partial correlation analysis, we find that $\mathrm{\Sigma_{*}}$ alone is insufficient to predict the N/O in MaNGA spaxels, and that there is an additional dependence on the local star formation rate surface density, $\mathrm{\Sigma_{SFR}}$. This effect is a factor of $3$ stronger than the dependence of 12+log(O/H) on $\mathrm{\Sigma_{SFR}}$. Surprisingly, we find that the local N/O scaling relations also depend on the total galaxy stellar mass at fixed $\Sigma_{*}$ as well as the galaxy size at fixed stellar mass. We find that more compact galaxies are more nitrogen rich, even when $\mathrm{\Sigma_{*}}$ and $\mathrm{\Sigma_{SFR}}$ are controlled for. We show that $\sim50\%$ of the variance of N/O is explained by the total stellar mass and size. Thus, the evolution of nitrogen in galaxies is set by more than just local effects and does not simply track the build up of oxygen in galaxies. The precise form of the N/O-O/H relation is therefore sensitive to the sample of galaxies from which it is derived. This result casts doubt on the universal applicability of nitrogen-based strong-line metallicity indicators derived in the local universe., Comment: 22 pages, 14 figures. Accepted for publication in the Astrophysical Journal

Bars may induce morphological features, such as rings, through their resonances. Previous studies suggested that the presence of 'dark-gaps', or regions of a galaxy where the difference between the surface brightness along the bar major axis and along the bar minor axis are maximal, can be attributed to the location of bar corotation. Here, using GALAKOS, a high-resolution N-body simulation of a barred galaxy, we test this photometric method's ability to identify the bar corotation resonance. Contrary to previous work, our results indicate that 'dark-gaps' are a clear sign of the location of the 4:1 ultra-harmonic resonance instead of bar corotation. Measurements of the bar corotation can indirectly be inferred using kinematic information, e.g., by measuring the shape of the rotation curve. We demonstrate our concept on a sample of 578 face-on barred galaxies with both imaging and integral field observations and find the sample likely consists primarily of fast bars., Comment: 16 pages, 14 figures, Accepted for publication in ApJ

We show the results of a study using the spectral synthesis technique study for the full MaNGA sample showing their Chemical Enrichment History (ChEH) as well as the evolution of the stellar mass-metallicity relation (MZR) over cosmic time. We find that the more massive galaxies became enriched first and the lower mass galaxies did so later, producing a change in the MZR which becomes shallower in time. Separating the sample into morphology and star-forming status bins some particularly interesting results appear: The mass dependency of the MZR becomes less relevant for later morphological types, to the extent that it inverts for Sd/Irr galaxies, suggesting that morphology is at least as important a factor as mass in chemical evolution. The MZR for the full sample shows a flattening at the high-mass end and another at the low-mass range, but the former only appears for retired galaxies while the latter only appears for star-forming galaxies. We also find that the average metallicity gradient is currently negative for all mass bins but for low mass galaxies it was inverted at some point in the past, before which all galaxies had a positive gradient. We also compare how diverse the ChEHs are in the different bins considered as well as what primarily drives the diversity: How much galaxies become enriched or how quickly they do so., Comment: Currently in review at ApJ

We present the MaNGA FIREFLY Value-Added-Catalogue (VAC) - a catalogue of ~3.7 million spatially resolved stellar population properties across 10,010 nearby galaxies from the final data release of the MaNGA survey. The full spectral fitting code firefly is employed to derive parameters such as stellar ages, metallicities, stellar and remnant masses, star formation histories, star formation rates and dust attenuation. In addition to Voronoi-binned measurements, our VAC also provides global properties, such as central values and radial gradients. Two variants of the VAC are available: presenting the results from fits using the M11-MILES and the novel MaStar stellar population models. MaStar allows to constrain the fit over the whole MaNGA wavelength range, extends the age-metallicity parameter space, and uses empirical spectra from the same instrument as MaNGA. The fits employing MaStar models find on average slightly younger ages, higher mass-weighted metallicities and smaller colour excesses. These differences are reduced when matching wavelength range and converging template grids. We further report that FIREFLY stellar masses are systematically lower by ~0.3 dex than masses from the MaNGA PCA and Pipe3D VACs, but match masses from the NSA best with only ~0.1 dex difference. Finally, we show that FIREFLY stellar ages correlate with spectral index age indicators H$\delta_A$ and $D_n$(4000), though with a clear additional metallicity dependence., Comment: 20 pages, 16 figures (+appendix). Accepted for publication in MNRAS. The accepted version now also includes star formation rates and performance tests. The MaNGA FIREFLY VAC is publicly available at the SDSS webpage https://www.sdss.org/dr17/manga/manga-data/manga-firefly-value-added-catalog and at ICG Portsmouth's website http://www.icg.port.ac.uk/manga-firefly-vac

The Sloan Digital Sky Survey MaNGA program has now obtained integral field spectroscopy for over 10,000 galaxies in the nearby universe. We use the final MaNGA data release DR17 to study the correlation between ionized gas velocity dispersion and galactic star formation rate, finding a tight correlation in which sigma_Ha from galactic HII regions increases significantly from ~ 18-30 km/s broadly in keeping with previous studies. In contrast, sigma_Ha from diffuse ionized gas (DIG) increases more rapidly from 20-60 km/s. Using the statistical power of MaNGA, we investigate these correlations in greater detail using multiple emission lines and determine that the observed correlation of sigma_Ha with local star formation rate surface density is driven primarily by the global relation of increasing velocity dispersion at higher total SFR, as are apparent correlations with stellar mass. Assuming HII region models consistent with our finding that sigma_[O III] < sigma_Ha < sigma_[O I], we estimate the velocity dispersion of the molecular gas in which individual HII regions are embedded, finding values sigma_Mol = 5-30 km/s consistent with ALMA observations in a similar mass range. Finally, we use variations in the relation with inclination and disk azimuthal angle to constrain the velocity dispersion ellipsoid of the ionized gas sigma_z/sigma_r = 0.84 +- 0.03 and sigma_phi/sigma_r = 0.91 +- 0.03, similar to that of young stars in the Galactic disk. Our results are most consistent with theoretical models in which turbulence in modern galactic disks is driven primarily by star formation feedback., Comment: 23 pages, 18 figures. Revised version accepted for publication in ApJ

The MaNGA Stellar Library (MaStar) is a large collection of high-quality empirical stellar spectra designed to cover all spectral types and ideal for use in the stellar population analysis of galaxies observed in the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey. The library contains 59,266 spectra of 24,130 unique stars with spectral resolution $R\sim1800$ and covering a wavelength range of $3,622-10,354$ \r{A}. In this work, we derive five physical parameters for each spectrum in the library: effective temperature ($T_{eff}$), surface gravity ($\log g$), metalicity ($[Fe/H]$), micro-turbulent velocity ($\log(v_{micro})$), and alpha-element abundance ($[\alpha/Fe]$). These parameters are derived with a flexible data-driven algorithm that uses a neural network model. We train a neural network using the subset of 1,675 MaStar targets that have also been observed in the Apache Point Observatory Galactic Evolution Experiment (APOGEE), adopting the independently-derived APOGEE Stellar Parameter and Chemical Abundance Pipeline (ASPCAP) parameters for this reference set. For the regions of parameter space not well represented by the APOGEE training set ($7,000 \leq T \leq 30,000$ K), we supplement with theoretical model spectra. We present our derived parameters along with an analysis of the uncertainties and comparisons to other analyses from the literature., Comment: 16 figures, 3 tables, accepted to AJ

We calculate the fundamental stellar parameters effective temperature, surface gravity and iron abundance - T$_{\rm eff}$, log g, [Fe/H] - for the final release of the Mapping Nearby Galaxies at APO (MaNGA) Stellar Library (MaStar), containing 59,266 per-visit-spectra for 24,290 unique stars at intermediate resolution ($R\sim1800$) and high S/N (median = 96). We fit theoretical spectra from model atmospheres by both MARCS and BOSZ-ATLAS9 to the observed MaStar spectra, using the full spectral fitting code pPXF. We further employ a Bayesian approach, using a Markov Chain Monte Carlo (MCMC) technique to map the parameter space and obtain uncertainties. Originally in this paper, we cross match MaStar observations with Gaia photometry, which enable us to set reliable priors and identify outliers according to stellar evolution. In parallel to the parameter determination, we calculate corresponding stellar population models to test the reliability of the parameters for each stellar evolutionary phase. We further assess our procedure by determining parameters for standard stars such as the Sun and Vega and by comparing our parameters with those determined in the literature from high-resolution spectroscopy (APOGEE and SEGUE) and from lower-resolution matching template (LAMOST). The comparisons, considering the different methodologies and S/N of the literature surveys, are favourable in all cases. Our final parameter catalogue for MaStar cover the following ranges: $2592 \leq $ T$_{\rm eff} \leq 32983\;$K; $-0.7 \leq $ log g $ \leq 5.4\;$dex; $-2.9 \leq $ [Fe/H] $\leq 1.0\;$dex and will be available with the last SDSS-IV Data Release, in December 2021., Comment: 22 pages, 25 figures. MNRAS in press. Stellar parameters will be made available with the Data Release 17 of the SDSS, planned for December 2021

We describe the survey design, calibration, commissioning, and emission-line detection algorithms for the Hobby-Eberly Telescope Dark Energy Experiment (HETDEX). The goal of HETDEX is to measure the redshifts of over a million Ly$\alpha$ emitting galaxies between 1.88

The Hobby-Eberly Telescope (HET) Dark Energy Experiment (HETDEX) is undertaking a blind wide-field low-resolution spectroscopic survey of 540 square degrees of sky to identify and derive redshifts for a million Lyman-alpha emitting galaxies (LAEs) in the redshift range 1.9 < z < 3.5. The ultimate goal is to measure the expansion rate of the Universe at this epoch, to sharply constrain cosmological parameters and thus the nature of dark energy. A major multi-year wide field upgrade (WFU) of the HET was completed in 2016 that substantially increased the field of view to 22 arcminutes diameter and the pupil to 10 meters, by replacing the optical corrector, tracker, and prime focus instrument package and by developing a new telescope control system. The new, wide-field HET now feeds the Visible Integral-field Replicable Unit Spectrograph (VIRUS), a new low-resolution integral field spectrograph (LRS2), and the Habitable Zone Planet Finder (HPF), a precision near-infrared radial velocity spectrograph. VIRUS consists of 156 identical spectrographs fed by almost 35,000 fibers in 78 integral field units arrayed at the focus of the upgraded HET. VIRUS operates in a bandpass of 3500-5500 Angstroms with resolving power R~800. VIRUS is the first example of large scale replication applied to instrumentation in optical astronomy to achieve spectroscopic surveys of very large areas of sky. This paper presents technical details of the HET WFU and VIRUS, as flowed-down from the HETDEX science requirements, along with experience from commissioning this major telescope upgrade and the innovative instrumentation suite for HETDEX., Comment: 65 pages, 25 figures, published in the Astronomical Journal; replaced with final published version

The distribution of stellar metallicities within and across galaxies is an excellent relic of the chemical evolution across cosmic time. We present a detailed analysis of spatially resolved stellar populations based on $>2.6$ million spatial bins from 7439 nearby galaxies in the SDSS-IV MaNGA survey. To account for accurate inclination corrections, we derive an equation for morphology dependent determination of galaxy inclinations. Our study goes beyond the well-known global mass-metallicity relation and radial metallicity gradients by providing a statistically sound exploration of local relations between stellar metallicity $[Z/H]$, stellar surface mass density $\Sigma_\star$ and galactocentric distance in the global mass-morphology plane. We find a significant resolved mass density-metallicity relation $\rm r\Sigma_\star ZR$ for galaxies of all types and masses above $10^{9.8}\,\mathrm{M_\odot}$. Different radial distances make an important contribution to the spread of the relation. Particularly, in low and intermediate mass galaxies, we find that at fixed $\Sigma_\star$ metallicity increases with radius independently of morphology. For high masses, this radial dependence is only observed in high $\Sigma_\star$ regions of spiral galaxies. This result calls for a driver of metallicity, in addition to $\Sigma_\star$ that promotes chemical enrichment in the outer parts of galaxies more strongly than in the inner parts. We discuss gas accretion, outflows, recycling and radial migration as possible scenarios., Comment: 16 pages, 11 figures. MNRAS accepted, pending a few minor revisions

We introduce \cosha{}: a Code for Stellar properties Heuristic Assignment. In order to estimate the stellar properties, \cosha{} implements a Gradient Tree Boosting algorithm to label each star across the parameter space ($T_\mathrm{eff}$, $\log{g}$, $[\mathrm{Fe}/\mathrm{H}]$, and $[\alpha/\mathrm{Fe}]$). We use \cosha{} to estimate these stellar atmospheric parameters of $22\,$k unique stars in the MaNGA Stellar Library (MaStar). To quantify the reliability of our approach, we run both internal tests using the G\"ottingen Stellar Library (GSL, a theoretical library) and the first data release of MaStar, and external tests by comparing the resulting distributions in the parameter space with the APOGEE estimates of the same properties. In summary, our parameter estimates span in the ranges: $T_\mathrm{eff}=[2900,12000]\,$K, $\log{g}=[-0.5,5.6]$, $[\mathrm{Fe}/\mathrm{H}]=[-3.74,0.81]$, $[\alpha/\mathrm{Fe}]=[-0.22,1.17]$. {We report internal (external) uncertainties of the properties of $\sigma_{T_\mathrm{eff}}\sim43\,(240)\,$K, $\sigma_{\log{g}}\sim0.2\,(0.4)$, $\sigma_{[\mathrm{Fe}/\mathrm{H}]}\sim0.16\,(0.24)$, $\sigma_{[\alpha/\mathrm{Fe}]}\sim0.09\,(0.08)$.} These uncertainties are comparable to those of other methods with similar objectives. Despite the fact that \cosha{} is not aware of the spatial distribution of these physical properties in the Milky Way, we are able to recover the main trends known in the literature. The catalog of physical properties for MaStar can be accessed in \url{http://ifs.astroscu.unam.mx/MaStar}., Comment: 25 pages, 12 figures, published in ApJS

The fundamental nature and extent of the coronal line region (CLR), which may serve as a vital tracer for Active Galactic Nucleus (AGN) activity, remain unresolved. Previous studies suggest that the CLR is produced by AGN-driven outflows and occupies a distinct region between the broad line region and the narrow line region, which places it tens to hundreds of parsecs from the galactic center. Here, we investigate 10 coronal line (CL; ionization potential $\ge$ 100 eV) emitting galaxies from the SDSS-IV MaNGA catalog with emission from one or more CLs detected at $\ge$ $5{\sigma}$ above the continuum in at least 10 spaxels - the largest such MaNGA catalog. We find that the CLR is far more extended, reaching out to 1.3 - 23 kpc from the galactic center. We cross-match our sample of 10 CL galaxies with the largest existing MaNGA AGN catalog and identify 7 in it; two of the remaining three are galaxy mergers and the final one is an AGN candidate. Further, we measure the average CLR electron temperatures to range between 12,331 K - 22,530 K, slightly above the typical threshold for pure AGN photoionization ($\sim$ 20,000 K) and indicative of shocks (e.g., merger-induced or from supernova remnants) in the CLR. We reason that ionizing photons emitted by the central continuum source (i.e. AGN photoionization) primarily generate the CLs, and that energetic shocks are an additional ionization mechanism that likely produce the most extended CLRs we measure., Comment: 25 pages, 19 figures

We study cool neutral gas traced by NaD absorption in 140 local ($\rm z<0.1)$ early-type ``red geyser'' galaxies. These galaxies show unique signatures in spatially-resolved strong-line emission maps that have been interpreted as large-scale active galactic nuclei driven ionized winds. To investigate the possible fuel source for these winds, we examine the abundance and kinematics of cool gas ($\rm T \sim 100-1000 K$) inferred from Na I D absorption in red geysers and matched control samples drawn from SDSS-IV MaNGA. We find that red geysers host greater amounts of NaD-associated material. Substantial cool gas components are detected in more than $\rm 50 \%$ of red geysers (compared to 25\% of the control sample) going up to 78$\%$ for radio-detected red geysers. Our key result is that cool gas in red geysers is predominantly infalling. Among our 30 radio-detected red geysers, 86$\%$ show receding NaD absorption velocities (with respect to the systemic velocity) between $\rm 40 - 50~km~s^{-1}$. We verify this result by stacking NaD profiles across each sample which confirms the presence of infalling NaD velocities within red geysers ( $\sim\rm 40~km~s^{-1}$) with no velocity offsets detected in the control samples. Interpreting our observations as signatures of inflowing cool neutral clouds, we derive an approximate mass inflow rate of $\rm \dot{M}_{in} \sim 0.1 M_{\odot} yr^{-1}$, similar to that expected from minor merging and internal recycling. Some red geysers show much higher rates ($\rm \dot{M}_{in} \sim 5 M_{\odot} yr^{-1}$) that may indicate an ongoing accretion event., Comment: 16 pages, 14 figures, 2 tables. Accepted for publication in the Astrophysical journal

We assemble a sample of 17 low metallicity (7.45 < log(O/H)+12 < 8.12) galaxies with z < 0.1 found spectroscopically, without photometric pre-selection, in early data from the Hobby Eberly Telescope Dark Energy Experiment (HETDEX). Star forming galaxies that occupy the lowest mass and metallicity end of the mass-metallicity relation tend to be under sampled in continuum-based surveys as their spectra are typically dominated by emission from newly forming stars. We search for galaxies with high [OIII]$\lambda$5007 / [OII]$\lambda$3727, implying highly ionized nebular emission often indicative of low metallicity systems. With the Second Generation Low Resolution Spectrograph on the Hobby Eberly Telescope we acquired follow-up spectra, with higher resolution and broader wavelength coverage, of each low-metallicity candidate in order to confirm the redshift, measure the H$\alpha$ and [NII] line strengths and, in many cases, obtain deeper spectra of the blue lines. We find our galaxies are consistent with the mass-metallicity relation of typical low mass galaxies. However, galaxies in our sample tend to have similar specific star formation rates (sSFRs) as the incredibly rare "blueberry" galaxies found in (Yang et. al. 2017). We illustrate the power of spectroscopic surveys for finding low mass and metallicity galaxies and reveal that we find a sample of galaxies that are a hybrid between the properties of typical dwarf galaxies and the more extreme blueberry galaxies., Comment: 32 pages, 12 figures, 8 tables, Accepted for publication in ApJ

As available data sets grow in size and complexity, advanced visualization tools enabling their exploration and analysis become more important. In modern astronomy, integral field spectroscopic galaxy surveys are a clear example of increasing dimensionality and complexity of datasets, which challenge the traditional methods used to extract the physical information they contain. We present the use of a novel self-supervised Machine Learning method to visualize the multi-dimensional information on stellar population and kinematics in the MaNGA survey in a two dimensional plane. Our framework is insensitive to non-physical properties such as the size of integral field unit (IFU) and is therefore able to order galaxies according to their resolved physical properties. Using the extracted representations, we study how galaxies distribute based on their resolved and global physical properties. We show that even when using exclusively information about the internal structure, galaxies naturally cluster into two well-known categories from a purely data driven perspective: rotating main-sequence disks and massive slow rotators, hence confirming distinct assembly channels. Low-mass rotation-dominated quenched galaxies appear only as a third cluster if information about the integrated physical properties is preserved, suggesting a mixture of assembly processes for these galaxies without any particular signature in their internal kinematics that distinguishes them from the two main groups. The framework for data exploration is publicly released with this publication, ready to be used with the MaNGA or other integral field data sets., Comment: 17 pages, 9 figures. Accepted for publication in The Astrophysical Journal

We present an analytical model to identify thin discs in galaxies, and apply this model to a sample of SDSS MaNGA galaxies. This model fits the velocity and velocity dispersion fields of galaxies with regular kinematics. By introducing two parameters $\zeta$ related to the comparison of the model's asymmetric drift correction to the observed gas kinematics and $\eta$ related to the dominant component of a galaxy, we classify the galaxies in the sample as "disc-dominated", "non-disc-dominated", or "disc-free" indicating galaxies with a dominating thin disc, a non-dominating thin disc, or no thin disc detection with our method, respectively. The dynamical mass resulting from our model correlates with stellar mass, and we investigate discrepancies by including gas mass and variation of the initial mass function. As expected, most spiral galaxies in the sample are disc-dominated, while ellipticals are predominantly disc-free. Lenticular galaxies show a dichotomy in their kinematic classification, which is related to their different star formation rates and gas fractions. We propose two possible scenarios to explain these results. In the first scenario, disc-free lenticulars formed in more violent processes than disc-dominated ones, while in the second scenario, the quenching processes in lenticulars lead to a change in their kinematic structures as disc-dominated lenticulars evolve to disc-free ones., Comment: 13 pages, 16 figures, accepted by MNRAS

Spatially resolved spectroscopy from SDSS-IV MaNGA survey has revealed a class of quiescent, relatively common early-type galaxies, termed "red geysers", that possibly host large scale active galactic nuclei driven winds. Given their potential importance in maintaining low level of star formation at late times, additional evidence confirming that winds are responsible for the red geyser phenomenon is critical. In this work, we present follow-up observations with the Echellette Spectrograph and Imager (ESI) at the Keck telescope of two red geysers (z$<$0.1) using multiple long slit positions to sample different regions of each galaxy. Our ESI data with a spectral resolution (R) $\sim$ 8000 improves upon MaNGA's resolution by a factor of four, allowing us to resolve the ionized gas velocity profiles along the putative wind cone with an instrumental resolution of $\rm \sigma = 16~km~s^{-1}$. The line profiles of H$\alpha$ and [NII]$\rm \lambda 6584$ show asymmetric shapes that depend systematically on location $-$ extended blue wings on the red-shifted side of the galaxy and red wings on the opposite side. We construct a simple wind model and show that our results are consistent with geometric projections through an outflowing conical wind oriented at an angle towards the line of sight. An alternative hypothesis that assigns the asymmetric pattern to "beam-smearing" of a rotating, ionized gas disk does a poor job matching the line asymmetry profiles. While our study features just two sources, it lends further support to the notion that red geysers are the result of galaxy-scale winds., Comment: 20 pages, 23 figures, accepted for publication in the Astrophysical journal

The Habitable zone Planet Finder (HPF) is a fiber fed precise radial velocity spectrograph at the 10 m Hobby Eberly Telescope (HET). Due to its fixed altitude design, the HET pupil changes appreciably across a track, leading to significant changes of the fiber far-field illumination. HPF's fiber scrambler is designed to suppress the impact of these illumination changes on the radial velocities -- but the residual impact on the radial velocity measurements has yet to be probed on sky. We use GJ 411, a bright early type (M2) M dwarf to probe the effects of far-field input trends due to these pupil variations on HPF radial velocities (RVs). These large changes ($\sim$ 2x) in pupil area and centroid present a harsh test of HPF's far-field scrambling. Our results show that the RVs are effectively decoupled from these extreme far-field input changes due to pupil centroid offsets, attesting to the effectiveness of the scrambler design. This experiment allows us to test the impact of these changes with large pupil variation on-sky, something we would not easily be able to do at a conventional optical telescope. While the pupil and illumination changes expected at these other telescopes are small, scaling from our results enables us to estimate and bound these effects, and show that they are controllable even for the new and next generation of RV instruments in their quest to beat down instrumental noise sources towards the goal of a few cm/s., Comment: 16 pages. Published in The Astrophysical Journal

We investigate the dust attenuation in both stellar populations and ionized gas in kpc-scale regions in nearby galaxies, using integral field spectroscopy data from MaNGA MPL-9. We identify star-forming (HII) and diffuse ionized gas (DIG) regions from MaNGA datacubes. From the stacked spectrum of each region, we measure the stellar attenuation, $E(B-V)_{\rm star}$, using the technique developed by Li et al.(2020), as well as the gas attenuation, $E(B-V)_{\rm gas}$, from the Balmer decrement. We then examine the correlation of $E(B-V)_{\rm star}$, $E(B-V)_{\rm gas}$, $E(B-V)_{\rm gas}-E(B-V)_{\rm star}$ and $E(B-V)_{\rm star}/E(B-V)_{\rm gas}$ with 16 regional/global properties, and for regions with different $\rm H{\alpha}$ surface brightnesses ($\Sigma_{\rm H\alpha}$). We find a stronger correlation between $E(B-V)_{\rm star}$ and $E(B-V)_{\rm gas}$ in regions of higher $\Sigma_{\rm H\alpha}$. Luminosity-weighted age ($t_L$) is found to be the property that is the most strongly correlated with $E(B-V)_{\rm star}$, and consequently with $E(B-V)_{\rm gas}-E(B-V)_{\rm star}$ and $E(B-V)_{\rm star}/E(B-V)_{\rm gas}$. At fixed $\Sigma_{\rm H\alpha}$, $\log_{10}t_L$ is linearly and negatively correlated with $E(B-V)_{\rm star}/E(B-V)_{\rm gas}$ at all ages. Gas-phase metallicity and ionization level are important for the attenuation in the gas. Our results indicate that the ionizing source for DIG regions is likely distributed in the outer-skirt of galaxies, while for HII regions our results can be well explained by the two-component dust model of Charlot & Fall (2000)., Comment: accepted by ApJ

We derive ages, metallicities, and individual element abundances of early- and late-type galaxies (ETGs and LTGs) out to 1.5 R$_e$. We study a large sample of 1900 galaxies spanning $8.6 - 11.3 \log M/M_{\odot}$ in stellar mass, through key absorption features in stacked spectra from the SDSS-IV/MaNGA survey. We use mock galaxy spectra with extended star formation histories to validate our method for LTGs and use corrections to convert the derived ages into luminosity- and mass-weighted quantities. We find flat age and negative metallicity gradients for ETGs and negative age and negative metallicity gradients for LTGs. Age gradients in LTGs steepen with increasing galaxy mass, from $-0.05\pm0.11~\log$ Gyr/R$_e$ for the lowest mass galaxies to $-0.82\pm0.08~\log$ Gyr/R$_e$ for the highest mass ones. This strong gradient-mass relation has a slope of $-0.70\pm0.18$. Comparing local age and metallicity gradients with the velocity dispersion $\sigma$ within galaxies against the global relation with $\sigma$ shows that internal processes regulate metallicity in ETGs but not age, and vice versa for LTGs. We further find that metallicity gradients with respect to local $\sigma$ show a much stronger dependence on galaxy mass than radial metallicity gradients. Both galaxy types display flat [C/Fe] and [Mg/Fe], and negative [Na/Fe] gradients, whereas only LTGs display gradients in [Ca/Fe] and [Ti/Fe]. ETGs have increasingly steep [Na/Fe] gradients with local $\sigma$ reaching $6.50\pm0.78$ dex/$\log$ km/s for the highest masses. [Na/Fe] ratios are correlated with metallicity for both galaxy types across the entire mass range in our sample, providing support for metallicity dependent supernova yields., Comment: 21 pages, 21 figures, 4 tables + Appendix