The launch of the James Webb Space Telescope will open up a new window for observations at the highest redshifts, reaching out to $z\approx 15$. However, even with this new facility, the first stars will remain out of reach, as they are born in small minihalos with luminosities too faint to be detected even by the longest exposure times. In this Letter, we investigate the basic properties of the "Ultimately Large Telescope", a facility that can detect Population III stars. Observations will take place in the near-infrared and therefore a moon-based facility is proposed. An instrument needs to reach magnitudes as faint as 39mag$_\mathrm{AB}$, corresponding to a primary mirror size of about 100m in diameter. Assuming JWST NIRCam filters, we estimate that Pop~III sources will have unique signatures in a colour-colour space and can be identified unambiguously.

Gas interior to the bar of the Milky Way has recently been shown as the closest example of a Low Ionization (Nuclear) Emission Region--LI(N)ER--in the universe. To better understand the nature of this gas, a sample of face-on galaxies with integral field spectroscopy are used to study the ionized gas conditions of 240 barred and 250 nonbarred galaxies, focusing on those that are most similar to the Milky Way. Strong optical line emission of $[NII]$ $\lambda 6584$, H$\alpha$, $[OIII]$ $\lambda 5007$, and H$\beta$ are used to diagnose the dominant ionization mechanisms of gas across galaxies and the Galaxy via Baldwin-Phillips-Terlevich (BPT) Diagrams. Barred galaxies show a strong suppression of star formation and an increase in composite and LI(N)ER like spectra in their inner regions when compared with similar nonbarred counterparts. This effect is lessened in galaxies of very low ($\log_{10}(M_\star/M_\odot) \lesssim 10.4$) or very high ($\log_{10}(M_\star/M_\odot) \gtrsim 11.1$) total stellar mass. Bar masks from Galaxy Zoo:3D show the bar's non-axisymmetric effect on the ionized gas and help predict the face-on distribution of ionized gas conditions near the bar of the Milky Way.

We report the stellar atmospheric parameters for 7503 spectra contained in the first release of the MaNGA stellar library (MaStar) in SDSS DR15. The first release of MaStar contains 8646 spectra measured from 3321 unique stars, each covering the wavelength range 3622 \AA\ to 10354 \AA\ with a resolving power of $R \sim$ 1800. In this work, we first determined the basic stellar parameters: effective temperature ($\rm T_{eff}$), surface gravity ($\log g$), and metallicity ($\rm[Fe/H]$), which best fit the data using an empirical interpolator based on the Medium-resolution Isaac Newton Telescope library of empirical spectra (MILES), as implemented by the University of Lyon Spectroscopic analysis Software (Koleva et al. 2008, ULySS) package. While we analyzed all 8646 spectra from the first release of MaStar, since MaStar has a wider parameter-space coverage than MILES, not all of these fits are robust. In addition, not all parameter regions covered by MILES yield robust results, likely due to the non-uniform coverage of the parameter space by MILES. We tested the robustness of the method using the MILES spectra itself and identified a proxy based on the local density of the training set. With this proxy, we identified 7503 MaStar spectra with robust fitting results. They cover the range from 3179K to 20,517K in effective temperature ($\rm T_{eff}$), from 0.40 to 5.0 in surface gravity ($\log g$), and from $-$2.49 to $+$0.73 in metallicity ($\rm[Fe/H]$).

The distribution of galaxies in optical diagnostic diagrams can provide information about their physical parameters when compared with ionization models under proper assumptions. By using a sample of central emitting regions from the MaNGA survey, we find evidence of the existence of upper boundaries for narrow-line regions (NLRs) of active galactic nuclei (AGN) in optical BPT diagrams, especially in the diagrams involving [S II]$\lambda \lambda$6716, 6731/H$\alpha$. Photoionization models can well reproduce the boundaries as a consequence of the decrease of [S II]$\lambda \lambda$6716, 6731/H$\alpha$ and [O III]$\lambda$5007/H$\beta$ ratios at very high metallicity. Whilst the exact location of the upper boundary in the [S II] BPT diagram only weakly depends on the electron density of the ionized cloud and the secondary nitrogen prescription, its dependence on the shapes of the input SEDs is much stronger. This allows us to constrain the power-law index of the AGN SED between 1 Ryd and $\sim100$ Ryd to be less than or equal to $-1.40\pm 0.05$. The coverage of the photoionization models in the [N II] BPT diagram has a stronger dependence on the electron density and the secondary nitrogen prescription. With the density constrained by the [S II] doublet ratio and the input SED constrained by the [S II] BPT diagram, we find that the extent of the data in the [N II] BPT diagram favors those prescriptions with high N/O ratios. Although shock-ionized cloud can produce similar line ratios as those by photoionization, the resulting shapes of the upper boundaries, if exist, would likely be different from those of a photoionizing origin.

We use the first release of the SDSS/MaStar stellar library comprising ~9000, high S/N spectra, to calculate integrated spectra of stellar population models. The models extend over the wavelength range 0.36-1.03 micron and share the same spectral resolution (R~1800) and flux calibration as the SDSS-IV/MaNGA galaxy data. The parameter space covered by the stellar spectra collected thus far allows the calculation of models with ages and chemical composition in the range t>200 Myr, -2 <=[Z/H]<= + 0.35, which will be extended as MaStar proceeds. Notably, the models include spectra for dwarf Main Sequence stars close to the core H-burning limit, as well as spectra for cold, metal-rich giants. Both stellar types are crucial for modelling lambda>0.7 micron absorption spectra. Moreover, a better parameter coverage at low metallicity allows the calculation of models as young as 500 Myr and the full account of the Blue Horizontal Branch phase of old populations. We present models adopting two independent sets of stellar parameters (T_eff, logg, [Z/H]). In a novel approach, their reliability is tested 'on the fly' using the stellar population models themselves. We perform tests with Milky Way and Magellanic Clouds globular clusters, finding that the new models recover their ages and metallicities remarkably well, with systematics as low as a few per cent for homogeneous calibration sets. We also fit a MaNGA galaxy spectrum, finding residuals of the order of a few per cent comparable to the state-of-art models, but now over a wider wavelength range.

We present the first asymmetric drift (AD) measurements for unresolved stellar populations of different characteristic ages above and below 1.5 Gyr. These measurements sample the age-velocity relation (AVR) in galaxy disks. In this first paper we develop two efficient algorithms to extract AD on a spaxel-by-spaxel basis from optical integral-field spectroscopic (IFS) data-cubes. The algorithms apply different spectral templates, one using simple stellar populations and the other a stellar library; their comparison allows us to assess systematic errors in derived multi-component velocities, such as template-mismatch. We test algorithm reliability using mock spectra and Monte Carlo Markov Chains on real data from the MaNGA survey in SDSS-IV. We quantify random and systematic errors in AD as a function of signal-to-noise and stellar population properties with the aim of applying this technique to large subsets of the MaNGA galaxy sample. As a demonstration of our methods, we apply them to an initial sample of seven galaxies with comparable stellar mass and color to the Milky Way. We find a wide range of distinct AD radial profiles for young and old stellar populations.

Bars inhabit the majority of local-Universe disk galaxies and may be important drivers of galaxy evolution through the redistribution of gas and angular momentum within disks. We investigate the star formation and gas properties of bars in galaxies spanning a wide range of masses, environments, and star formation rates using the MaNGA galaxy survey. Using a robustly-defined sample of 684 barred galaxies, we find that fractional (or scaled) bar length correlates with the host's offset from the star-formation main sequence. Considering the morphology of the H$\alpha$ emission we separate barred galaxies into different categories, including barred, ringed, and central configurations, together with H$\alpha$ detected at the ends of a bar. We find that only low-mass galaxies host star formation along their bars, and that this is located predominantly at the leading edge of the bar itself. Our results are supported by recent simulations of massive galaxies, which show that the position of star formation within a bar is regulated by a combination of shear forces, turbulence and gas flows. We conclude that the physical properties of a bar are mostly governed by the existing stellar mass of the host galaxy, but that they also play an important role in the galaxy's ongoing star formation.

This paper documents the sixteenth data release (DR16) from the Sloan Digital Sky Surveys; the fourth and penultimate from the fourth phase (SDSS-IV). This is the first release of data from the southern hemisphere survey of the Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2); new data from APOGEE-2 North are also included. DR16 is also notable as the final data release for the main cosmological program of the Extended Baryon Oscillation Spectroscopic Survey (eBOSS), and all raw and reduced spectra from that project are released here. DR16 also includes all the data from the Time Domain Spectroscopic Survey (TDSS) and new data from the SPectroscopic IDentification of ERosita Survey (SPIDERS) programs, both of which were co-observed on eBOSS plates. DR16 has no new data from the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey (or the MaNGA Stellar Library "MaStar"). We also preview future SDSS-V operations (due to start in 2020), and summarize plans for the final SDSS-IV data release (DR17).

We perform a "fossil record" analysis for ~800 low-redshift spiral galaxies, using STARLIGHT applied to integral field spectroscopic observations from the SDSS-IV MaNGA survey to obtain fully spatially-resolved high-resolution star formation histories (SFHs). From the SFHs, we are able to build maps indicating the present-day distribution of stellar populations of different ages in each galaxy. We find small negative mean age gradients in most spiral galaxies, especially at high stellar mass, which reflects the formation times of stellar populations at different galactocentric radii. We show that the youngest (<10^{8.5} years) populations exhibit significantly more extended distributions than the oldest (>10^{9.5} years), again with a strong dependence on stellar mass. By interpreting the radial profiles of "time slices" as indicative of the size of the galaxy at the time those populations had formed, we are able to trace the simultaneous growth in mass and size of the spiral galaxies over the last 10 Gyr. Despite finding that the evolution of the measured light-weighted radius is consistent with inside-out growth in the majority of spiral galaxies, the evolution of an equivalent mass-weighted radius has changed little over the same time period. Since radial migration effects are likely to be small, we conclude that the growth of disks in spiral galaxies has occurred predominantly through an inside-out mode (with the effect greatest in high-mass galaxies), but this has not had anywhere near as much impact on the distribution of mass within spiral galaxies.

Dust attenuation in star-forming spiral galaxies affects stars and gas in different ways due to local variations in dust geometry. We present spatially resolved measurements of dust attenuation for a sample of 232 such star-forming spiral galaxies, derived from spectra acquired by the SDSS-IV MaNGA survey. The dust attenuation affecting the stellar populations of these galaxies (obtained using full spectrum stellar population fitting methods) is compared with the dust attenuation in the gas (derived from the Balmer decrement). Both of these attenuation measures increase for local regions of galaxies with higher star formation rates; the dust attenuation affecting the stellar populations increases more so than the dust attenuation in the gas, causing the ratio of the dust attenuation affecting the stellar populations to the dust attenuation in the gas to decrease for local regions of galaxies with higher star formation rate densities. No systematic difference is discernible in any of these dust attenuation quantities between the spiral arm and inter-arm regions of the galaxies. While both the dust attenuation in the gas and the dust attenuation affecting the stellar populations decrease with galactocentric radius, the ratio of the two quantities does not vary with radius. This ratio does, however, decrease systematically as the stellar mass of the galaxy increases. Analysis of the radial profiles of the two dust attenuation measures suggests that there is a disproportionately high concentration of birth clouds (incorporating gas, young stars and clumpy dust) nearer to the centres of star-forming spiral galaxies.

We study the stellar populations and assembly of the nearby spiral galaxy NGC 2903's bulge, bar, and outer disc using the VIRUS-P Exploration of Nearby Galaxies IFS survey. We observe NGC 2903 with a spatial resolution of 185 pc using the Mitchell Spectrograph's 4.25 arcsec fibres at the 2.7 Harlan J. Smith telescope. Bulge-bar-disc decomposition on the 2MASS Ks-band image of NGC 2903 shows that it has ~6%, 6%, and 88%, of its stellar mass in the bulge, bar, and outer disc, respectively, and its bulge has a low Sersic index of ~0.27, suggestive of a disky bulge. We perform stellar population synthesis and find that the outer disc has 46% of its mass in stars >5 Gyr, 48% in stars between 1 and 5 Gyr, and <10% in younger stars. Its stellar bar has 65% of its mass in ages 1-5 Gyr and has metallicities similar to the outer disc, suggestive of the evolutionary picture where the bar forms from disc material. Its bulge is mainly composed of old high-metallicity stars though it also has a small fraction of young stars. We find enhanced metallicity in the spiral arms and central region, tracing areas of high star formation as seen in the Halpha map. These results are consistent with the idea that galaxies of low bulge-to-total mass ratio and low bulge Sersic index like NGC 2903 has not had a recent major merger event, but has instead grown mostly through minor mergers and secular processes.

We measured gas-phase metallicity, ionisation parameter and dust extinction for 1795 representative local star-forming galaxies using integral field spectroscopy from the SDSS-IV MaNGA survey. We self-consistently derive these quantities by comparing observed line fluxes with photoionisation models using a Bayesian framework. We also present the first comprehensive study of the [SIII]$\lambda\lambda$9069,9532 nebular lines, which have long been predicted to be ideal tracers of the ionisation parameter. Unfortunately, we find that current photoionisation models substantially over-predict [SIII] lines intensity, while broadly reproducing other observed optical line ratios. We discuss how to nonetheless make use of the information provided by [SIII] lines by setting a prior on the ionisation parameter. Following this approach, we derive spatially-resolved maps and radial profiles of metallicity and ionisation parameter. The metallicity radial profiles are comparable with previous works, with metallicity declining toward the outer parts and a flattening in the central regions, in agreement with infall models of galaxy formation, that predict that spiral discs build up through accretion of material, which leads to an inside-out growth. On the other hand, ionisation parameter radial profiles are flat for low-mass galaxies, while their slope becomes positive as galaxy mass increases. However, the ionisation parameter maps we obtain are clumpy, especially for low-mass galaxies. Ionisation parameter is tightly correlated with the H$\alpha$ equivalent width [EW(H$\alpha$)], following a nearly universal relation, which we attribute to the change of the spectral shape of ionising sources due to ageing of HII regions. We derive a positive correlation between ionisation parameter and metallicity at fixed EW(H$\alpha$), in disagreement with previous theoretical works expecting an anti-correlation.

The Milky Way provides an ideal laboratory to test our understanding of galaxy evolution, owing to our ability to observe our Galaxy over fine scales. However, connecting the Galaxy to the wider galaxy population remains difficult, due to the challenges posed by our internal perspective and to the different observational techniques employed. Here, we present a sample of galaxies identified as Milky Way Analogs (MWAs) on the basis of their stellar masses and bulge-to-total ratios, observed as part of the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey. We analyse the galaxies in terms of their stellar kinematics and populations as well as their ionised gas contents. We find our sample to contain generally young stellar populations in their outskirts. However, we find a wide range of stellar ages in their central regions, and we detect central AGN-like or composite-like activity in roughly half of the sample galaxies, with the other half consisting of galaxies with central star-forming emission or emission consistent with old stars. We measure gradients in gas metallicity and stellar metallicity that are generally flatter in physical units than those measured for the Milky Way; however, we find far better agreement with the Milky Way when scaling gradients by galaxies' disc scale lengths. From this, we argue much of the discrepancy in metallicity gradients to be due to the relative compactness of the Milky Way, with differences in observing perspective also likely to be a factor.

In this paper we use strong line calibrations of N/O and O/H in MaNGA spaxel data to explore the systematics introduced by variations in N/O on various strong-line metallicity diagnostics. We find radial variations in N/O at fixed O/H which correlate with total galaxy stellar-mass; and which can induce $\sim 40 \%$ systematic uncertainties in oxygen abundance gradients when nitrogen-dependent abundance calibrations are used. Empirically, we find that these differences are associated with variation in the local star formation efficiency, as predicted by recent chemical evolution models for galaxies, but we cannot rule out other processes such as radial migration and the accretion of passive dwarf galaxies also playing a role.

We present the first release of the MaNGA Stellar Library (MaStar), which is a large, well-calibrated, high-quality empirical library covering the wavelength range of 3,622-10,354A at a resolving power of R~1800. The spectra were obtained using the same instrument as used by the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) project, by piggybacking on the SDSS-IV/APOGEE-2N observations. Compared to previous empirical libraries, the MaStar library will have a higher number of stars and a more comprehensive stellar-parameter coverage, especially of cool dwarfs, low-metallicity stars, and stars with different [alpha/Fe], achieved by a sophisticated target selection strategy that takes advantage of stellar-parameter catalogs from the literature. This empirical library will provide a new basis for stellar population synthesis, and is particularly well-suited for stellar-population analysis of MaNGA galaxies. The first version of the library contains 8646 high-quality per-visit spectra for 3321 unique stars. Compared to photometry, the relative flux calibration of the library is accurate to 3.9% in g-r, 2.7% in r-i, and 2.2% in i-z. The data are released as part of Sloan Digital Sky Survey Data Release 15. We expect the final release of the library to contain more than 10,000 stars.

Galaxy interactions leave imprints in the motions of their stars, and so observing the two-dimensional stellar kinematics allows us to uncover their formation process. Slow rotators, which have stellar orbits dominated by random motions, are thought to be the fossil relics of a sequence of multiple gas-poor mergers, in an environment where the cold gas required to form new stars is nearly absent. Indeed, observations of a handful of nearby galaxy clusters have indicated that slow rotators are preferentially found in the gas-poor, dense cores of clusters, which themselves must form by merging of smaller groups. However, the generality of this result and connection between kinematics and environment is currently unclear, as recent studies have suggested that, at given stellar mass, the environment does not influence the formation of slow rotators. Here we address this issue by combining a careful quality-assessed sample selection with two-dimensional stellar kinematics from a large galaxy survey and a novel photometric classification approach where kinematics are unavailable. We obtain the first complete census of the location of massive slow rotators in ten large clusters: in all cases, slow rotators are extremely rare and generally trace the clusters density peaks. This result unambiguously establishes that massive slow rotators are the relics of violent hierarchical cluster formation.

Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) is acquiring integral-field spectroscopy for the largest sample of galaxies to date. By 2020, the MaNGA Survey --- one of three core programs in the fourth-generation Sloan Digital Sky Survey (SDSS-IV) --- will have observed a statistically representative sample of 10$^4$ galaxies in the local Universe ($z\lesssim0.15$). In addition to a robust data-reduction pipeline (DRP), MaNGA has developed a data-analysis pipeline (DAP) that provides higher-level data products. To accompany the first public release of its code base and data products, we provide an overview of the MaNGA DAP, including its software design, workflow, measurement procedures and algorithms, performance, and output data model. In conjunction with our companion paper Belfiore et al., we also assess the DAP output provided for 4718 observations of 4648 unique galaxies in the recent SDSS Data Release 15 (DR15). These analysis products focus on measurements that are close to the data and require minimal model-based assumptions. Namely, we provide stellar kinematics (velocity and velocity dispersion), emission-line properties (kinematics, fluxes, and equivalent widths), and spectral indices (e.g., D4000 and the Lick indices). We find that the DAP provides robust measurements and errors for the vast majority ($>$99%) of analyzed spectra. We summarize assessments of the precision and accuracy of our measurements as a function of signal-to-noise, and provide specific guidance to users regarding the limitations of the data. The MaNGA DAP software is publicly available and we encourage community involvement in its development.

It has been observed that low redshift early-type galaxies can be separated into slow and fast rotators according to a proxy of specific stellar angular momentum, $\lambda_{R_e}$. Detailed studies of a handful of nearby clusters have shown that slow rotators are generally found at the centres of clusters where the number density is highest, whereas the fast rotators trace the trend followed by early-type galaxies of increasing in number with local density. In this paper, we study the environmental distribution of slow and fast rotators using the stellar kinematics of about 3900 galaxies from the Sloan Digital Sky Survey's Mapping Nearby Galaxies at Apache Point Observatory survey. For galaxies in groups closer than $z=0.08$ that are not observed with MaNGA but satisfy the necessary conditions for slow rotators, we visually assign slow/fast rotator classifications to obtain a complete sample. Our final catalogue contains about 900 groups of five or more members. We observe the kinematic morphology-density (kT-$\Sigma$) relation for each group and find an increasing fraction of massive slow rotators with increasing number density. We provide evidence suggesting that the observed lack of trends in angular momentum with environment at fixed stellar mass is partly because the maximum density varies between clusters, and that the locations of massive slow rotators are strongly correlated with peak densities in galaxy groups and clusters. We conclude that the (projected) number density relative to the cluster peak is more fundamental than the absolute number density in influencing the abundance of slow rotators. We find that the kT-$\Sigma$ relation does exist at fixed stellar mass, and we rule out the hypothesis that the kT-$\Sigma$ relation is a result of dynamical friction alone, instead arguing that massive slow rotators grow hierarchically in tandem with their host clusters.

Galaxy morphology is inextricably linked to environment. The morphology-density relation quantifies this relationship. However, optical morphology is only loosely related to the kinematic structure of galaxies, and about two thirds of elliptical ("spheroidal") galaxies are actually misclassified face-on disks, and would appear flattened if view edge on. A more robust classification is the slow/fast rotator classification which describes the kinematic structure. Slow and fast rotators form a bimodality in galaxy properties and are thought to follow distinct evolutionary paths, and so a kinematic morphology-density (kT-$\Sigma$) relation is more meaningful. To date the kT-$\Sigma$ relation has only been studied for a handful of nearby clusters,or across large numbers of clusters but with incomplete coverage. Here, we combine stellar kinematics obtained with the Sloan Digital Sky Survey's (SDSS) Mapping Nearby Galaxies at Apache Point Observatory survey with classifications using a novel visual method to obtain the largest complete census of slow and fast rotators in nearby galaxy groups and clusters. To account for incompleteness in the SDSS spectroscopic catalogue, we combine the catalogue with the photometric catalogue, which we clean using empirical criteria. We test our visual classification method and find our false negative rate to be only $\sim$6% (0% for $M \gtrsim 10^{11.7}M_{\odot}$). In contrast, our false positive rate is about 50% implying half of all slow rotator candidates will be confirmed as fast rotators if stellar kinematics become available. Hence, our misclassification is essentially random with only a weak dependence on stellar mass, and as slow rotators are intrinsically rare, the absolute number of misclassifications will be small compared to the sample size. The result is the largest complete census of massive slow rotator candidates and their environments to date.

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.

SDSS-IV MaNGA (Mapping Nearby Galaxies at Apache Point Observatory) is the largest integral-field spectroscopy survey to date, aiming to observe a statistically representative sample of 10,000 low-redshift galaxies. In this paper we study the reliability of the emission-line fluxes and kinematic properties derived by the MaNGA Data Analysis Pipeline (DAP). We describe the algorithmic choices made in the DAP with regards to measuring emission-line properties, and the effect of our adopted strategy of simultaneously fitting the continuum and line emission. The effect of random errors are quantified by studying various fit-quality metrics, idealized recovery simulations and repeat observations. This analysis demonstrates that the emission lines are well-fit in the vast majority of the MaNGA dataset and the derived fluxes and errors are statistically robust. The systematic uncertainty on emission-line properties introduced by the choice of continuum templates is also discussed. In particular, we test the effect of using different stellar libraries and simple stellar-population models on the derived emission-line fluxes and the effect of introducing different tying prescriptions for the emission-line kinematics. We show that these effects can generate large ($>$ 0.2 dex) discrepancies at low signal-to-noise and for lines with low equivalent width (EW); however, the combined effect is noticeable even for H$\alpha$ EW $>$ 6~\AA. We provide suggestions for optimal use of the data provided by SDSS data release 15 and propose refinements on the \DAP\ for future MaNGA data releases.

We assemble an unbiased sample of 29 galaxies with [O II] $\lambda 3727$ and/or [O III] $\lambda 5007$ detections at $z < 0.15$ from the Hobby-Eberly Telescope Dark Energy Experiment (HETDEX) Pilot Survey (HPS). HPS finds galaxies without pre-selection based on their detected emission lines via integral field spectroscopy. Sixteen of these objects were followed up with the second-generation, low resolution spectrograph (LRS2) on the upgraded Hobby-Eberly Telescope. Oxygen abundances were then derived via strong emission lines using a Bayesian approach. We find most of the galaxies fall along the mass-metallicity relation derived from photometrically selected star forming galaxies in the Sloan Digital Sky Survey (SDSS). However, two of these galaxies have low metallicity (similar to the very rare green pea galaxies in mass-metallicity space). The star formation rates of this sample fall in an intermediate space between the SDSS star forming main sequence and the extreme green pea galaxies. We conclude that spectroscopic selection fills part of the mass-metallicity-SFR phase space that is missed in photometric surveys with pre-selection like SDSS, i.e., we find galaxies that are actively forming stars but are faint in continuum. We use the results of this pilot investigation to make predictions for the upcoming unbiased, large spectroscopic sample of local line emitters from HETDEX. With the larger HETDEX survey we will determine if galaxies selected spectroscopically without continuum brightness pre-selection have metallicities that fall on a continuum that bridges typical star forming and rarer, more extreme systems like green peas.

We present estimates of stellar population (SP) gradients from stacked spectra of slow (SR) and fast (FR) rotator elliptical galaxies from the MaNGA-DR15 survey. We find that: 1) FRs are $\sim 5$ Gyrs younger, more metal rich, less $\alpha$-enhanced and smaller than SRs of the same luminosity $L_r$ and central velocity dispersion $\sigma_0$. This explains why when one combines SRs and FRs, objects which are small for their $L_r$ and $\sigma_0$ tend to be younger. Their SP gradients are also different. 2) Ignoring the FR/SR dichotomy leads one to conclude that compact galaxies are older than their larger counterparts of the same mass, even though almost the opposite is true for FRs and SRs individually. 3) SRs with $\sigma_0\le 250$ km s$^{-1}$ are remarkably homogeneous within $\sim R_e$: they are old, $\alpha$-enhanced and only slightly super-solar in metallicity. These SRs show no gradients in age and $M_*/L_r$, negative gradients in metallicity, and slightly positive gradients in [$\alpha$/Fe] (the latter are model dependent). SRs with $\sigma_0\ge 250$ km $s^{-1}$ are slightly younger and more metal rich, contradicting previous work suggesting that age increases with $\sigma_0$. They also show larger $M_*/L_r$ gradients. 4) Self-consistently accounting for $M_*/L$ gradients yields $M_{\rm dyn}\approx M_*$ because gradients reduce $M_{\rm dyn}$ by $\sim 0.2$ dex while only slightly increasing the $M_*$ inferred using a Kroupa (not Salpeter) IMF. 5) The FR population all but disappears above $M_*\ge 3\times 10^{11}M_\odot$; this is the same scale at which the size-mass correlation and other scaling relations change. Our results support the finding that this is an important mass scale which correlates with the environment and above which mergers matter.

We estimate ages, metallicities, $\alpha$-element abundance ratios and stellar initial mass functions of elliptical (E) and S0 galaxies from the MaNGA-DR15 survey. We stack spectra and use a variety of single stellar population synthesis models to interpret the absorption line strengths in these spectra. We quantify how these properties vary across the population, as well as with galactocentric distance. This paper is the first of a series and is based on a sample of pure elliptical galaxies at z $\le$ 0.08. We show that the properties of the inner regions of Es with the largest luminosity (L$_r$) and central velocity dispersion ($\sigma_0$) are consistent with those associated with the commonly used Salpeter IMF, whereas a Kroupa-like IMF is a better description at $\sim$ 0.8R/Re (assuming [Ti/Fe] variations are limited). For these galaxies the stellar mass-to-light ratio decreases at most by a factor of 2 from the central regions to Re. In contrast, for lower L$_r$ and $\sigma_0$ galaxies, the IMF is shallower and M$_{*}$/L$_r$ in the central regions is similar to the outskirts. Although a factor of 2 is smaller than previous reports based on a handful of galaxies, it is still large enough to matter for dynamical mass estimates. Accounting self-consistently for these gradients when estimating both M$_{*}$ and M$_{dyn}$ brings the two into good agreement: gradients reduce M$_{dyn}$ by $\sim$ 0.2 dex while only slightly increasing the M$_{*}$ inferred using a Kroupa IMF. This is a different resolution of the M$_{*}$-M$_{dyn}$ discrepancy than has been followed in the recent literature where M$_{*}$ of massive galaxies is increased by adopting a Salpeter IMF while leaving Mdyn unchanged. A companion paper discusses how stellar population differences are even more pronounced if one separates slow from fast rotators.

Spectra of galaxies contain a wealth of information about the stellar populations from which they are made. With integral-field unit (IFU) surveys, such data can be used to map out stellar population properties across the face of a galaxy, allowing one to go beyond simple radial profiles and study details of non-axisymmetric structure. To-date, however, such studies have been limited by the quality of available data and the power of spectral analysis tools. We now take the next step and study the barred spiral galaxy MCG+07-28-064 from observations obtained as part of the SDSS-IV MaNGA project. We find that we can decompose this galaxy into "time slices," which reveal the varying contributions that stars of differing ages make to its bar and spiral structure, offering new insight into the evolution of these features. We find evidence for the ongoing growth of the bar, including the most recent star formation on its leading edge, and for the underlying density wave responsible for spiral structure. This pilot study indicates that there is a wealth of untapped information on the spatial distribution of SFHs available in the current generation of IFU galaxy surveys.

A galaxy's stellar mass is one of its most fundamental properties, but it remains challenging to measure reliably. With the advent of very large optical spectroscopic surveys, efficient methods that can make use of low signal-to-noise spectra are needed. With this in mind, we created a new software package for estimating effective stellar mass-to-light ratios $\log \Upsilon^*$ that uses principal component analysis(PCA) basis set to optimize the comparison between observed spectra and a large library of stellar population synthesis models. In Paper I, we showed that a with a set of six PCA basis vectors we could faithfully represent most optical spectra from the Mapping Nearby Galaxies at APO (MaNGA) survey;and we tested the accuracy of our M/L estimates using synthetic spectra. Here, we explore sources of systematic error in our mass measurements by comparing our new measurements to data from the literature. We compare our stellar mass surface density estimates to kinematics-derived dynamical mass surface density measurements from the DiskMass Survey and find some tension between the two which could be resolved if the disk scale-heights used in the kinematic analysis were overestimated by a factor of $\sim 1.5$. We formulate an aperture-corrected stellar mass catalog for the MaNGA survey, and compare to previous stellar mass estimates based on multi-band optical photometry, finding typical discrepancies of 0.1 dex. Using the spatially resolved MaNGA data, we evaluate the impact of estimating total stellar masses from spatially unresolved spectra, and we explore how the biases that result from unresolved spectra depend upon the galaxy's dust extinction and star formation rate. Finally, we describe a SDSS Value-Added Catalog which will include both spatially resolved and total (aperture-corrected) stellar masses for MaNGA galaxies.

We present a method of fitting optical spectra of galaxies using a basis set of six vectors obtained from principal component analysis (PCA) of a library of synthetic spectra of 40000 star formation histories (SFHs). Using this library, we provide estimates of resolved effective stellar mass-to-light ratio ($\log \Upsilon^*$) for thousands of galaxies from the SDSS-IV/MaNGA integral-field spectroscopic survey. Using a testing framework built on additional synthetic SFHs, we show that the estimates of stellar mass-to-light ratio are reliable (as are their uncertainties) at a variety of signal-to-noise ratios, stellar metallicities, and dust attenuation conditions. Finally, we describe the future release of the resolved stellar mass-to-light ratios as a SDSS-IV/MaNGA Value-Added Catalog (VAC) and provide a link to the software used to conduct this analysis.

The role of gas accretion in galaxy evolution is still a matter of debate. The presence of inflows of metal-poor gas that trigger star formation bursts of low metallicity has been proposed as an explanation for the local anti-correlation between star formation rate (SFR) and gas-phase metallicity ($\rm Z_g$) found in the literature. In the present study, we show how the anti-correlation is also present as part of a diversified range of behaviours for a sample of more than 700 nearby spiral galaxies from the SDSS IV MaNGA survey. We have characterized the local relation between SFR and $\rm Z_g$ after subtracting the azimuthally averaged radial profiles of both quantities. $60\%$ of the analyzed galaxies display a $\rm SFR-Z_g$ anti-correlation, with the remaining $40\%$ showing no correlation ($19\%$) or positive correlation ($21\%$). Applying a Random Forest machine-learning algorithm, we obtain that the slope of the correlation is mainly determined by the average gas-phase metallicity of the galaxy. Galaxy mass, $g-r$ colors, stellar age, and mass density seem to play a less significant role. This result is supported by the performed 2nd-order polynomial regression analysis. Thus, the local $\rm SFR-Z_g$ slope varies with the average metallicity, with the more metal-poor galaxies presenting the lowest slopes (i.e., the strongest $\rm SFR-Z_g$ anti-correlations), and reversing the relation for more metal-rich systems. Our results suggest that external gas accretion fuels star-formation in metal-poor galaxies, whereas in metal-rich systems the gas comes from previous star formation episodes.

The late assembly of massive galaxies is thought to be dominated by stellar accretion in their outskirts (beyond 2 effective radii Re) due to dry, minor galaxy mergers. We use observations of 1010 passive early-type galaxies (ETGs) within z<0.15 from SDSS IV MaNGA to search for evidence of this accretion. The outputs from the stellar population fitting codes FIREFLY, pPXF, and Prospector are compared to control for systematic errors in stellar metallicity (Z) estimation. We find that the average radial logZ/Zsun profiles of ETGs in various stellar mass (M) bins are not linear. As a result, these profiles are poorly characterized by a single gradient value, explaining why weak trends reported in previous work can be difficult to interpret. Instead, we examine the full radial extent of stellar metallicity profiles and find them to flatten in the outskirts of M>10^{11}Msun ETGs. This is a signature of stellar accretion. Based on a toy model for stellar metallicity profiles, we infer the ex-situ stellar mass fraction in ETGs as a function of M and galactocentric radius. We find that ex-situ stars at 2Re make up 20% of the projected stellar mass of M<10^{10.5}Msun ETGs, rising up to 80% for M>10^{11.5}Msun ETGs.

Bars in galaxies are thought to stimulate both inflow of material and radial mixing along them. Observational evidence for this mixing has been inconclusive so far however, limiting the evaluation of the impact of bars on galaxy evolution. We now use results from the MaNGA integral field spectroscopic survey to characterise radial stellar age and metallicity gradients along the bar and outside the bar in 128 strongly barred galaxies. We find that age and metallicity gradients are flatter in the barred regions of almost all barred galaxies when compared to corresponding disk regions at the same radii. Our results re-emphasize the key fact that by azimuthally averaging integral field spectroscopic data one loses important information from non-axisymmetric galaxy components such as bars and spiral arms. We interpret our results as observational evidence that bars are radially mixing material in galaxies of all stellar masses, and for all bar morphologies and evolutionary stages.

The exact nature of the arms of spiral galaxies is still an open question. It has been widely assumed that spiral arms in galaxies with two distinct symmetrical arms are the products of density waves that propagate around the disk, with the spiral arms being visibly enhanced by the star formation that is triggered as the passing wave compresses gas in the galaxy disk. Such a persistent wave would propagate with an approximately constant angular speed, its pattern speed Op. The quasi-stationary density wave theory can be tested by measuring this quantity and showing that it does not vary with radius in the galaxy. Unfortunately, this measurement is difficult because Op is only indirectly connected to observables such as the stellar rotation speed. Here, we use the detailed information on stellar populations of the grand-design spiral galaxy UGC 3825, extracted from spectral mapping, to measure the offset between young stars of a known age and the spiral arm in which they formed, allowing the first direct measure of Op at a range of radii. The offset in this galaxy is found to be as expected for a pattern speed that varies little with radius, indicating consistency with a quasi-stationary density wave, and lending credence to this new method.

Crucial progress in our understanding of star formation and feedback will depend on the ability to obtain spatially resolved spectroscopic observations of \ion{H}{ii} regions, from which reliable instantaneous measurements of their physical conditions can be obtained. Acquiring these datasets across full galactic systems will prove crucial for obtaining population samples that enable us to understand the time evolution of similar regions, and the variability of conditions among coeval regions. Separating the spatial and temporal dependencies in such way for different physical processes involved in star formation and the injection of feedback is crucial to overcome the inherit degeneracies associated with observing instantaneous snapshots of a dynamic ISM at any given time. Emission line diagnostics are at the core of measuring the physical condition in \ion{H}{ii} regions (e.g. dynamics, SFR, chemical abundances, dust extinction, ionization and excitation, etc.). These measurements require high spatial resolution, contiguous coverage across full galactic systems, and sensitivities significantly deeper than past efforts. The spatial scale required to resolve the \ion{H}{ii} regions of a few pc is only attainable in the Local Group where very large sky coverage is necessary.

Several models have predicted that stars could form inside galactic outflows and that this would be a new major mode of galaxy evolution. Observations of galactic outflows have revealed that they host large amounts of dense and clumpy molecular gas, which provide conditions suitable for star formation. We have investigated the properties of the outflows in a large sample of galaxies by exploiting the integral field spectroscopic data of the large MaNGA-SDSS4 galaxy survey. We find that star formation occurs inside at least half of the galactic outflows in our sample. We also show that even if star formation is prominent inside many other galactic outflows, this may have not been revealed as the diagnostics are easily dominated by the presence of even faint AGN and shocks. If very massive outflows typical of distant galaxies and quasars follow the same scaling relations observed locally, then the star formation inside high-z outflows can be up to several 100 Msun/yr and could contribute substantially to the early formation of the spheroidal component of galaxies. Star formation in outflows can also potentially contribute to establishing the scaling relations between black holes and their host spheroids. Moreover, supernovae exploding on large orbits can chemically enrich in-situ and heat the circumgalactic and intergalactic medium. Finally, young stars ejected on large orbits may also contribute to the reionization of the Universe.

By applying spectroscopic decomposition methods to a sample of MaNGA early-type galaxies, we separate out spatially and kinematically distinct stellar populations, allowing us to explore the similarities and differences between galaxy bulges and discs, and how they affect the global properties of the galaxy. We find that the components have interesting variations in their stellar populations, and display different kinematics. Bulges tend to be consistently more metal rich than their disc counterparts, and while the ages of both components are comparable, there is an interesting tail of younger, more metal poor discs. Bulges and discs follow their own distinct kinematic relationships, both on the plane of the stellar spin parameter, lambda_R, and ellipticity, and in the relation between stellar mass and specific angular momentum, j, with the location of the galaxy as a whole on these planes being determined by how much bulge and disc it contains. As a check of the physical significance of the kinematic decompositions, we also dynamically model the individual galaxy components within the global potential of the galaxy. The resulting components exhibit kinematic parameters consistent with those from the spectroscopic decomposition, and though the dynamical modelling suffers from some degeneracies, the bulges and discs display systematically different intrinsic dynamical properties. This work demonstrates the value in considering the individual components of galaxies rather than treating them as a single entity, which neglects information that may be crucial in understanding where, when and how galaxies evolve into the systems we see today.

Deriving a simple, analytic galaxy star formation history (SFH) using observational data is a complex task without the proper tool to hand. We therefore present SNITCH, an open source code written in Python, developed to quickly (~2 minutes) infer the parameters describing an analytic SFH model from the emission and absorption features of a galaxy spectrum dominated by star formation gas ionisation. SNITCH uses the Flexible Stellar Population Synthesis models of Conroy et al. (2009), the MaNGA Data Analysis Pipeline and a Markov Chain Monte Carlo method in order to infer three parameters (time of quenching, rate of quenching and model metallicity) which best describe an exponentially declining quenching history. This code was written for use on the MaNGA spectral data cubes but is customisable by a user so that it can be used for any scenario where a galaxy spectrum has been obtained, and adapted to infer a user-defined analytic SFH model for specific science cases. Herein we outline the rigorous testing applied to SNITCH and show that it is both accurate and precise at deriving the SFH of a galaxy spectra. The tests suggest that SNITCH is sensitive to the most recent epoch of star formation but can also trace the quenching of star formation even if the true decline does not occur at an exponential rate. With the use of both an analytical SFH and only five spectral features, we advocate that this code be used as a comparative tool across a large population of spectra, either for integral field unit data cubes or across a population of galaxy spectra.

The large Integral Field Spectroscopy (IFS) surveys have allowed the classification of ionizing sources of emission lines on sub-kpc scales. In this work, we define two non-parametric parameters, quiescence (f$_{q}$) and its concentration (c$_{q}$), to quantify the strength and the spatial distribution of the quenched areas, respectively, traced by the LI(N)ER regions with low EW(H$\alpha$). With these two measurements, we classify MaNGA galaxies into inside-out and outside-in quenching types according to their locations on the f$_{q}$ vs. c$_{q}$ plane and we measure the fraction of inside-out (outside-in) quenching galaxies as a function of halo mass. We find that the fraction of galaxies showing inside-out quenching increases with halo mass, irrespective of stellar mass or galaxy type (satellites vs. centrals). In addition, high stellar mass galaxies exhibit a greater fraction of inside-out quenching compared to low stellar mass ones in all environments. In contrast, the fraction of outside-in quenching does not depend on halo mass. Our results suggest that morphological quenching may be responsible for the inside-out quenching seen in all environments. On the other hand, the flat dependence of the outside-in quenching on halo mass could be a mixed result of ram-pressure stripping and galaxy mergers. Nevertheless, at a given environment and stellar mass, the fraction of inside-out quenching is systematically greater than that of outside-in quenching, suggesting that inside-out quenching is the dominant quenching mode in all environments.

We use data from 1222 late-type star-forming galaxies in the SDSS IV MaNGA survey to identify regions in which the gas-phase metallicity is anomalously-low compared to expectations from the tight empirical relation between metallicity and stellar surface mass-density at a given stellar mass. We find anomalously low metallicity (ALM) gas in 10% of the star-forming spaxels, and in 25% of the galaxies in the sample. The incidence rate of ALM gas increases strongly with both global and local measures of the specific star-formation rate, and is higher in lower mass galaxies and in the outer regions of galaxies. The incidence rate is also significantly higher in morphologically disturbed galaxies. We estimate that the lifetimes of the ALM regions are a few hundred Myr. We argue that the ALM gas has been delivered to its present location by a combination of interactions, mergers, and accretion from the halo, and that this infusion of gas stimulates star-formation. Given the estimated lifetime and duty cycle of such events, we estimate that the time-averaged accretion rate of ALM gas is similar to the star-formation rate in late type galaxies over the mass-range M$_* \sim10^9$ to 10$^{10}$ M$_{\odot}$.

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.

Chemical enrichment signatures strongly constrain galaxy formation and evolution, and a detailed understanding of abundance patterns provides clues regarding the nucleosynthetic production pathways of elements. Using the SDSS-IV MaNGA IFU survey, we study radial gradients of chemical element abundances in detail. We use stacked spectra out to 1 Re of 366 early-type galaxies with masses 9.9 - 10.8 log $M/M_{\odot}$ to probe the abundances of the elements C, N, Na, Mg, Ca, and Ti, relative to the abundance of Fe, by fitting stellar population models to a combination of Lick absorption indices. We find that C, Mg, and Ti trace each other both as a function of galaxy radius and galaxy mass. These similar C and Mg abundances within and across galaxies set a lower limit for star-formation timescales. Conversely, N and Ca are generally offset to lower abundances. The under-abundance of Ca compared to Mg implies delayed enrichment of Ca through Type Ia supernovae, whereas the correlated behaviour of Ti and the lighter $\alpha$ elements, C and Mg, suggest contributions to Ti from Type II supernovae. We obtain shallow radial gradients in [Mg/Fe], [C/Fe], and [Ti/Fe], meaning that these inferences are independent of radius. However, we measure strong negative radial gradients for [N/Fe] and [Na/Fe], of up to $-0.25\pm0.05$ and $-0.29\pm0.02$ dex/Re respectively. These gradients become shallower with decreasing galaxy mass. We find that N and Na abundances increase more steeply with velocity dispersion within galaxies than globally, while the other elements show the same relation locally and globally. This implies that the high Na and N abundances found in massive early type galaxies are generated by internal processes within galaxies. These are strongly correlated with the total metallicity, suggesting metallicity-dependent Na enrichment, and secondary N production in massive early-type galaxies.

We investigate the relationship of kinematically misaligned galaxies with their large-scale environment, in the context of halo assembly bias. According to numerical simulations, halo age at fixed halo mass is intrinsically linked to the large-scale tidal environment created by the cosmic web. We investigate the relationship between distances to various cosmic web features and present-time gas accretion rate. We select a sub-sample of ~900 central galaxies from the MaNGA survey with defined global position angles (PA; angle at which velocity change is greatest) for their stellar and H$\alpha$ gas components up to a minimum of 1.5 effective radii ($R_e$). We split the sample by misalignment between the gas and stars as defined by the difference in their PA. For each central galaxy we find its distance to nodes and filaments within the cosmic web, and estimate the host halo's age using the central stellar mass to total halo mass ratio $M_{*}/M_{h}$. We also construct halo occupation distributions using a background subtraction technique for galaxy groups split using the central galaxy's kinematic misalignment. We find, at fixed halo mass, no statistical difference in these properties between our kinematically aligned and misaligned galaxies. We suggest that the lack of correlation could be indicative of cooling flows from the hot halo playing a far larger role than `cold mode' accretion from the cosmic web or a demonstration that the spatial extent of current large-scale integral field unit (IFU) surveys hold little information about large-scale environment extractable through this method.

A new class of quiescent galaxies harboring possible AGN-driven winds has been discovered using spatially resolved optical spectroscopy from the ongoing SDSS-IV MaNGA survey. These galaxies, termed "red geysers", constitute $5-10\%$ of the local quiescent population and are characterized by narrow bisymmetric patterns in ionized gas emission features. Cheung et al. argued that these galaxies host large-scale AGN-driven winds that may play a role in suppressing star formation at late times. In this work, we test the hypothesis that AGN activity is ultimately responsible for the red geyser phenomenon. We compare the nuclear radio activity of the red geysers to a matched control sample with similar stellar mass, redshift, rest frame $NUV-r$ color, axis ratio and presence of ionized gas. We have used the 1.4 GHz radio continuum data from VLA FIRST survey to stack the radio flux from the red geyser and control samples. In addition to a 3 times higher FIRST detection rate, we find that red geysers have a 5$\sigma$ higher level of average radio flux than control galaxies. After restricting to rest-frame $NUV-r$ color $>$ 5 and checking mid-IR WISE photometry, we rule out star formation contamination and conclude that red geysers are associated with more active AGN. Red geysers and a possibly-related class with disturbed H$\alpha$ emission account for 40\% of all radio-detected red galaxies with $\rm log~(M_\star/M_\odot) < 11$. Our results support a picture in which episodic AGN activity drives large-scale-relatively weak ionized winds that may provide a feedback mechanism for many early-type galaxies.

Negative feedback from accretion onto super-massive black holes (SMBHs), that is to remove gas and suppress star formation in galaxies, has been widely suggested. However, for Seyfert galaxies which harbor less active, moderately accreting SMBHs in the local universe, the feedback capability of their black hole activity is elusive. We present spatially-resolved H$\alpha$ measurements to trace ongoing star formation in Seyfert galaxies and compare their specific star formation rate with a sample of star-forming galaxies whose global galaxy properties are controlled to be the same as the Seyferts. From the comparison we find that the star formation rates within central kpc of Seyfert galaxies are mildly suppressed as compared to the matched normal star forming galaxies. This suggests that the feedback of moderate SMBH accretion could, to some extent, regulate the ongoing star formation in these intermediate to late type galaxies under secular evolution.

Gas stripping of spiral galaxies or mergers are thought to be the formation mechanisms of lenticular galaxies. In order to determine the conditions in which each scenario dominates, we derive stellar populations of both the bulge and disk regions of 279 lenticular galaxies in the MaNGA survey. We find a clear bimodality in stellar age and metallicity within the population of S0s and this is strongly correlated with stellar mass. Old and metal-rich bulges and disks belong to massive galaxies, and young and metal-poor bulges and disks are hosted by low-mass galaxies. From this we conclude that the bulges and disks are co-evolving. When the bulge and disk stellar ages are compared, we find that the bulge is almost always older than the disk for massive galaxies ($\textrm{M}_{\star} > 10^{10}~\textrm{M}_{\odot}$). The opposite is true for lower mass galaxies. We conclude that we see two separate populations of lenticular galaxies. The old, massive, and metal-rich population possess bulges that are predominantly older than their disks, which we speculate may have been caused by morphological or inside-out quenching. In contrast, the less massive and more metal-poor population have bulges with more recent star formation than their disks. We postulate they may be undergoing bulge rejuvenation (or disk fading), or compaction. Environment doesn't play a distinct role in the properties of either population. Our findings give weight to the notion that while the faded spiral scenario likely formed low-mass S0s, other processes, such as mergers, may be responsible for high-mass S0s.

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

A large population of ultra-diffuse galaxies (UDGs) was recently discovered in the Coma cluster. Here we present optical spectra of three such UDGs, DF7, DF44 and DF17, which have central surface brightnesses of $\mu_g \approx 24.4-25.1$ mag arcsec$^{-2}$. The spectra were acquired as part of an ancillary program within the SDSS-IV MaNGA Survey. We stacked 19 fibers in the central regions from larger integral field units (IFUs) per source. With over 13.5 hours of on-source integration we achieved a mean signal-to-noise ratio (S/N) in the optical of $9.5$\AA$^{-1}$, $7.9$\AA$^{-1}$ and $5.0$\AA$^{-1}$, respectively, for DF7, DF44 and DF17. Stellar population models applied to these spectra enable measurements of recession velocities, ages and metallicities. The recession velocities of DF7, DF44 and DF17 are $6599^{+40}_{-25}$km/s, $6402^{+41}_{-39}$km/s and $8315^{+43}_{-43}$km/s, spectroscopically confirming that all of them reside in the Coma cluster. The stellar populations of these three galaxies are old and metal-poor, with ages of $7.9^{+3.6}_{-2.5}$Gyr, $8.9^{+4.3}_{-3.3}$Gyr and $9.1^{+3.9}_{-5.5}$Gyr, and iron abundances of $\mathrm{[Fe/H]}$ $-1.0^{+0.3}_{-0.4}$, $-1.3^{+0.4}_{-0.4}$ and $-0.8^{+0.5}_{-0.5}$, respectively. Their stellar masses are $3$-$6\times10^8 M_\odot$. The UDGs in our sample are as old or older than galaxies at similar stellar mass or velocity dispersion (only DF44 has an independently measured dispersion). They all follow the well-established stellar mass$-$stellar metallicity relation, while DF44 lies below the velocity dispersion-metallicity relation. These results, combined with the fact that UDGs are unusually large for their stellar mass, suggest that stellar mass plays a more important role in setting stellar population properties for these galaxies than either size or surface brightness.

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 study radial profiles in H$\alpha$ equivalent width and specific star formation rate (sSFR) derived from spatially-resolved SDSS-IV MaNGA spectroscopy to gain insight on the physical mechanisms that suppress star formation and determine a galaxy's location in the SFR-$\rm M_\star$ diagram. Even within the star-forming `main sequence', the measured sSFR decreases with stellar mass, both in an integrated and spatially-resolved sense. Flat sSFR radial profiles are observed for $\rm log(M_\star/ M_\odot) < 10.5$, while star-forming galaxies of higher mass show a significant decrease in sSFR in the central regions, a likely consequence of both larger bulges and an inside-out growth history. Our primary focus is the green valley, constituted by galaxies lying below the star formation main sequence, but not fully passive. In the green valley we find sSFR profiles that are suppressed with respect to star-forming galaxies of the same mass at all galactocentric distances out to 2 effective radii. The responsible quenching mechanism therefore appears to affect the entire galaxy, not simply an expanding central region. The majority of green valley galaxies of $\rm log(M_\star/ M_\odot) > 10.0$ are classified spectroscopically as central low-ionisation emission-line regions (cLIERs). Despite displaying a higher central stellar mass concentration, the sSFR suppression observed in cLIER galaxies is not simply due to the larger mass of the bulge. Drawing a comparison sample of star forming galaxies with the same $\rm M_\star$ and $\rm \Sigma_{1~kpc}$ (the mass surface density within 1 kpc), we show that a high $\rm \Sigma_{1~kpc}$ is not a sufficient condition for determining central quiescence.

MaNGA provides the opportunity to make precise spatially resolved measurements of the IMF slope in galaxies owing to its unique combination of spatial resolution, wavelength coverage and sample size. We derive radial gradients in age, element abundances and IMF slope analysing optical and near-infrared absorption features from stacked spectra out to the half-light radius of 366 early-type galaxies with masses $9.9 - 10.8\;\log M/M_{\odot}$. We find flat gradients in age and [$\alpha$/Fe] ratio, as well as negative gradients in metallicity, consistent with the literature. We further derive significant negative gradients in the [Na/Fe] ratio with galaxy centres being well enhanced in Na abundance by up to 0.5 dex. Finally, we find a gradient in IMF slope with a bottom-heavy IMF in the centre (typical mass excess factor of 1.5) and a Milky Way-type IMF at the half-light radius. This pattern is mass-dependent with the lowest mass galaxies in our sample featuring only a shallow gradient around a Milky Way IMF. Our results imply the local IMF-$\sigma$ relation within galaxies to be even steeper than the global relation and hint towards the local metallicity being the dominating factor behind the IMF variations. We also employ different stellar population models in our analysis and show that a radial IMF gradient is found independently of the stellar population model used. A similar analysis of the Wing-Ford band provides inconsistent results and further evidence of the difficulty in measuring and modelling this particular feature.

We study the kinematics and the stellar populations of the bulge and disc of the spiral galaxy NGC 3521. At each position in the field of view, we separate the contributions of the bulge and the disc from the total observed spectrum and study their kinematics, age, and metallicities independently. Their properties are clearly distinct: the bulge rotates more slowly, has a higher velocity dispersion, and is less luminous than the disc. We identify three main populations of stars in NGC 3521: old ($\geq7$ Gyr), intermediate ($\approx$ 3 Gyr), and young ($\leq$1 Gyr). The mass and light of NGC 3521 are dominated by the intermediate stellar population. The youngest population contributes mostly to the disc component and its contribution increases with radius. We also study the luminosity-weighed properties of the stars in NGC 3521. Along the photometric major axis, we find: i) no age gradient for the stars in the bulge, and a negative age gradient for the stars in the disc; ii) negative metallicity gradients and sub-solar $\alpha$-enhancement for both the bulge and the disc. We propose the following picture for the formation of NGC 3521: initial formation a long time ago ($\geq 7$ Gyr), followed by a second burst of star formation or a merger ($\approx$ 3 Gyrs ago), which contributed predominantly to the mass-build up of the bulge. Recently ($\leq 1$ Gyr), the disc of NGC 3521 experienced an additional episode of star formation that started in the innermost regions.

We measure $\lambda_{R_e}$, a proxy for galaxy specific stellar angular momentum within one effective radius, and the ellipticity, $\epsilon$, for about 2300 galaxies of all morphological types observed with integral field spectroscopy as part of the MaNGA survey, the largest such sample to date. We use the $(\lambda_{R_e}, \epsilon)$ diagram to separate early-type galaxies into fast and slow rotators. We also visually classify each galaxy according to its optical morphology and two-dimensional stellar velocity field. Comparing these classifications to quantitative $\lambda_{R_e}$ measurements reveals tight relationships between angular momentum and galaxy structure. In order to account for atmospheric seeing, we use realistic models of galaxy kinematics to derive a general approximate analytic correction for $\lambda_{R_e}$. Thanks to the size of the sample and the large number of massive galaxies, we unambiguously detect a clear bimodality in the $(\lambda_{R_e}, \epsilon)$ diagram which may result from fundamental differences in galaxy assembly history. There is a sharp secondary density peak inside the region of the diagram with low $\lambda_{R_e}$ and $\epsilon < 0.4$, previously suggested as the definition for slow rotators. Most of these galaxies are visually classified as non-regular rotators and have high velocity dispersion. The intrinsic bimodality must be stronger, as it tends to be smoothed by noise and inclination. The large sample of slow rotators allows us for the first time to unveil a secondary peak at +/-90 degrees in their distribution of the misalignments between the photometric and kinematic position angles. We confirm that genuine slow rotators start appearing above a stellar mass of 2\times10^{11} M_{\odot}$ where a significant number of high-mass fast rotators also exist.

We present new evidence for AGN feedback in a subset of 69 quenched low-mass galaxies ($M_{\star} \lesssim 5\times10^{9}$ M$_{\odot}$, $M_{\rm{r}} > -19$) selected from the first two years of the SDSS-IV MaNGA survey. The majority (85 per cent) of these quenched galaxies appear to reside in a group environment. We find 6 galaxies in our sample that appear to have an active AGN that is preventing on-going star-formation; this is the first time such a feedback mechanism has been observed in this mass range. Interestingly, five of these six galaxies have an ionised gas component that is kinematically offset from their stellar component, suggesting the gas is either recently accreted or outflowing. We hypothesise these six galaxies are low-mass equivalents to the "red geysers" observed in more massive galaxies. Of the other 63 galaxies in the sample, we find 8 do appear for have some low-level, residual star formation, or emission from hot, evolved stars. The remaining galaxies in our sample have no detectable ionised gas emission throughout their structures, consistent with them being quenched. This work shows the potential for understanding the detailed physical properties of dwarf galaxies through spatially resolved spectroscopy.

We study the spatially resolved star formation of 1494 galaxies in the SDSSIV-MaNGA Survey. SFRs are calculated using a two-step process, using $H_\alpha$ in star forming regions and $D_n4000$ in regions identified as AGN/LI(N)ER or lineless. The roles of secular and environmental quenching processes are investigated by studying the dependence of the radial profiles of specific star formation rate on stellar mass, galaxy structure and environment. We report on the existence of `Centrally Suppressed' galaxies, which have suppressed SSFR in their cores compared to their disks. The profiles of centrally suppressed and unsuppressed galaxies are distibuted in a bimodal way. Galaxies with high stellar mass and core velocity dispersion are found to be much more likely to be centrally suppressed than low mass galaxies, and we show that this is related to morphology and the presence of AGN/LI(N)ER like emission. Centrally suppressed galaxies also display lower star formation at all radii compared to unsuppressed galaxies. The profiles of central and satellite galaxies are also compared, and we find that satellite galaxies experience lower specific star formation rates at all radii than central galaxies. This uniform suppression could be a signal of the stripping of hot halo gas in the process known as strangulation. We find that satellites are not more likely to be suppressed in their cores than centrals, indicating that the core suppression is an entirely internal process. We find no correlation between the local environment density and the profiles of star formation rate surface density.

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

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.

We determine the local metallicity of the ionized gas for more than $9.2\times 10^5$ star forming regions (spaxels) located in 1023 nearby galaxies included in the SDSS-IV MaNGA IFU survey. We use the dust extinction derived from the Balmer decrement and stellar template fitting in each spaxel to estimate the local gas and stellar mass densities, respectively. We also use the measured rotation curves to determine the local escape velocity ($\mathrm{V_{esc}}$). We have then analyze the relationships between the local metallicity and both the local gas fraction ($\mu$) and $\mathrm{V_{esc}}$. We find that metallicity decreases with both increasing $\mu$ and decreasing $\mathrm{V_{esc}}$. By examining the residuals in these relations we show that the gas fraction plays a more primary role in the local chemical enrichment than $\mathrm{V_{esc}}$. We show that the gas-regulator model of chemical evolution provides a reasonable explanation of the metallicity on local scales. The best-fit parameters for this model are consistent with metal loss caused by momentum-driven galactic outflows. We also argue that both the gas fraction and local escape velocity are connected to the local stellar surface density, which in turn is a tracer of the epoch at which the dominant local stellar population formed

Regions of disc galaxies with widespread star formation tend to be both gravitationally unstable and self-shielded against ionizing radiation, whereas extended outer discs with little or no star formation tend to be stable and unshielded on average. We explore what drives the transition between these two regimes, specifically whether discs first meet the conditions for self-shielding (parameterized by dust optical depth, $\tau$) or gravitational instability (parameterized by a modified version of Toomre's instability parameters, $Q_{\rm thermal}$, which quantifies the stability of a gas disc that is thermally supported at $T=10^4$ K). We first introduce a new metric formed by the product of these quantities, $Q_{\rm thermal}\tau$, which indicates whether the conditions for disk instability or self-shielding are easier to meet in a given region of a galaxy, and we discuss how $Q_{\rm thermal}\tau$ can be constrained even in the absence of direct gas information. We then analyse a sample of 13 galaxies with resolved gas measurements and find that on average galaxies will reach the threshold for disk instabilities ($Q_{\rm thermal}<1$) before reaching the threshold for self-shielding ($\tau>1$). Using integral field spectroscopic observations of a sample of 236 galaxies from the MaNGA survey, we find that the value of $Q_{\rm thermal}\tau$ in star-forming discs is consistent with similar behavior. These results support a scenario where disc fragmentation and collapse occurs before self-shielding, suggesting that gravitational instabilities are the primary condition for widespread star formation in galaxy discs. Our results support similar conclusions based on recent galaxy simulations.

SDSS-V will be an all-sky, multi-epoch spectroscopic survey of over six million objects. It is designed to decode the history of the Milky Way, trace the emergence of the chemical elements, reveal the inner workings of stars, and investigate the origin of planets. It will also create an integral-field spectroscopic map of the gas in the Galaxy and the Local Group that is 1,000x larger than the current state of the art and at high enough spatial resolution to reveal the self-regulation mechanisms of galactic ecosystems. SDSS-V will pioneer systematic, spectroscopic monitoring across the whole sky, revealing changes on timescales from 20 minutes to 20 years. The survey will thus track the flickers, flares, and radical transformations of the most luminous persistent objects in the universe: massive black holes growing at the centers of galaxies. The scope and flexibility of SDSS-V will be unique among extant and future spectroscopic surveys: it is all-sky, with matched survey infrastructures in both hemispheres; it provides near-IR and optical multi-object fiber spectroscopy that is rapidly reconfigurable to serve high target densities, targets of opportunity, and time-domain monitoring; and it provides optical, ultra-wide-field integral field spectroscopy. SDSS-V, with its programs anticipated to start in 2020, will be well-timed to multiply the scientific output from major space missions (e.g., TESS, Gaia, eROSITA) and ground-based projects. SDSS-V builds on the 25-year heritage of SDSS's advances in data analysis, collaboration infrastructure, and product deliverables. The project is now refining its science scope, optimizing the survey strategies, and developing new hardware that builds on the SDSS-IV infrastructure. We present here an overview of the current state of these developments as we seek to build our worldwide consortium of institutional and individual members.

In this paper, we investigate 2727 galaxies observed by MaNGA as of June 2016 to develop spatially resolved techniques for identifying signatures of active galactic nuclei (AGN). We identify 303 AGN candidates. The additional spatial dimension imposes challenges in identifying AGN due to contamination from diffuse ionized gas, extra-planar gas and photoionization by hot stars. We show that the combination of spatially-resolved line diagnostic diagrams and additional cuts on H$\alpha$ surface brighness and H$\alpha$ equivalent width can distinguish between AGN-like signatures and high-metallicity galaxies with LINER-like spectra. Low mass galaxies with high specific star formation rates are particularly difficult to diagnose and routinely show diagnostic line ratios outside of the standard star-formation locus. We develop a new diagnostic -- the distance from the standard diagnostic line in the line-ratios space -- to evaluate the significance of the deviation from the star-formation locus. We find 173 galaxies that would not have been selected as AGN candidates based on single-fibre spectral measurements but exhibit photoionization signatures suggestive of AGN activity in the MaNGA resolved observations, underscoring the power of large integral field unit (IFU) surveys. A complete census of these new AGN candidates is necessary to understand their nature and probe the complex co-evolution of supermassive black holes and their hosts.

We study the role of cold gas in quenching star formation in the green valley by analysing ALMA $^{12}$CO (1-0) observations of three galaxies with resolved optical spectroscopy from the MaNGA survey. We present resolution-matched maps of the star formation rate and molecular gas mass. These data are used to calculate the star formation efficiency (SFE) and gas fraction ($f_{\rm~gas}$) for these galaxies separately in the central `bulge' regions and outer disks. We find that, for the two galaxies whose global specific star formation rate (sSFR) deviates most from the star formation main sequence, the gas fraction in the bulges is significantly lower than that in their disks, supporting an `inside-out' model of galaxy quenching. For the two galaxies where SFE can be reliably determined in the central regions, the bulges and disks share similar SFEs. This suggests that a decline in $f_{\rm~gas}$ is the main driver of lowered sSFR in bulges compared to disks in green valley galaxies. Within the disks, there exist common correlations between the sSFR and SFE and between sSFR and $f_{\rm~gas}$ on kpc scales -- the local SFE or $f_{\rm~gas}$ in the disks declines with local sSFR. Our results support a picture in which the sSFR in bulges is primarily controlled by $f_{\rm~gas}$, whereas both SFE and $f_{\rm~gas}$ play a role in lowering the sSFR in disks. A larger sample is required to confirm if the trend established in this work is representative of green valley as a whole.

Do the theorised different formation mechanisms of fast and slow rotators produce an observable difference in their star formation histories? To study this we identify quenching slow rotators in the MaNGA sample by selecting those which lie below the star forming sequence and identify a sample of quenching fast rotators which were matched in stellar mass. This results in a total sample of 194 kinematically classified galaxies, which is agnostic to visual morphology. We use u-r and NUV-u colours from SDSS and GALEX and an existing inference package, STARPY, to conduct a first look at the onset time and exponentially declining rate of quenching of these galaxies. An Anderson-Darling test on the distribution of the inferred quenching rates across the two kinematic populations reveals they are statistically distinguishable ($3.2\sigma$). We find that fast rotators quench at a much wider range of rates than slow rotators, consistent with a wide variety of physical processes such as secular evolution, minor mergers, gas accretion and environmentally driven mechanisms. Quenching is more likely to occur at rapid rates ($\tau \lesssim 1~\rm{Gyr}$) for slow rotators, in agreement with theories suggesting slow rotators are formed in dynamically fast processes, such as major mergers. Interestingly, we also find that a subset of the fast rotators quench at these same rapid rates as the bulk of the slow rotator sample. We therefore discuss how the total gas mass of a merger, rather than the merger mass ratio, may decide a galaxy's ultimate kinematic fate.

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.

The "kinematic" morphology-density relation for early-type galaxies posits that those galaxies with low angular momentum are preferentially found in the highest-density regions of the universe. We use a large sample of galaxy groups with halo masses 10^12.5 < M_halo < 10^14.5 M_sun/h observed with the Mapping Nearby Galaxies at APO (MaNGA) survey to examine whether there is a correlation between local environment and rotational support that is independent of stellar mass. We find no compelling evidence for a relationship between the angular momentum content of early-type galaxies and either local overdensity or radial position within the group at fixed stellar mass.

We describe the sample design for the SDSS-IV MaNGA survey and present the final properties of the main samples along with important considerations for using these samples for science. Our target selection criteria were developed while simultaneously optimizing the size distribution of the MaNGA integral field units (IFUs), the IFU allocation strategy, and the target density to produce a survey defined in terms of maximizing S/N, spatial resolution, and sample size. Our selection strategy makes use of redshift limits that only depend on i-band absolute magnitude ($M_i$), or, for a small subset of our sample, $M_i$ and color (NUV-i). Such a strategy ensures that all galaxies span the same range in angular size irrespective of luminosity and are therefore covered evenly by the adopted range of IFU sizes. We define three samples: the Primary and Secondary samples are selected to have a flat number density with respect to $M_i$ and are targeted to have spectroscopic coverage to 1.5 and 2.5 effective radii (Re), respectively. The Color-Enhanced supplement increases the number of galaxies in the low-density regions of color-magnitude space by extending the redshift limits of the Primary sample in the appropriate color bins. The samples cover the stellar mass range $5\times10^8 \leq M_* \leq 3\times10^{11} M_{\odot}$ and are sampled at median physical resolutions of 1.37 kpc and 2.5 kpc for the Primary and Secondary samples respectively. We provide weights that will statistically correct for our luminosity and color-dependent selection function and IFU allocation strategy, thus correcting the observed sample to a volume limited sample.

We describe the Sloan Digital Sky Survey IV (SDSS-IV), a project encompassing three major spectroscopic programs. The Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2) is observing hundreds of thousands of Milky Way stars at high resolution and high signal-to-noise ratio in the near-infrared. The Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey is obtaining spatially-resolved spectroscopy for thousands of nearby galaxies (median redshift of z = 0.03). The extended Baryon Oscillation Spectroscopic Survey (eBOSS) is mapping the galaxy, quasar, and neutral gas distributions between redshifts z = 0.6 and 3.5 to constrain cosmology using baryon acoustic oscillations, redshift space distortions, and the shape of the power spectrum. Within eBOSS, we are conducting two major subprograms: the SPectroscopic IDentification of eROSITA Sources (SPIDERS), investigating X-ray AGN and galaxies in X-ray clusters, and the Time Domain Spectroscopic Survey (TDSS), obtaining spectra of variable sources. All programs use the 2.5-meter Sloan Foundation Telescope at Apache Point Observatory; observations there began in Summer 2014. APOGEE-2 also operates a second near-infrared spectrograph at the 2.5-meter du Pont Telescope at Las Campanas Observatory, with observations beginning in early 2017. Observations at both facilities are scheduled to continue through 2020. In keeping with previous SDSS policy, SDSS-IV provides regularly scheduled public data releases; the first one, Data Release 13, was made available in July 2016.

We study the gas phase metallicity (O/H) and nitrogen abundance gradients traced by star forming regions in a representative sample of 550 nearby galaxies in the stellar mass range $\rm 10^9-10^{11.5} M_\odot$ with resolved spectroscopic data from the SDSS-IV MaNGA survey. Using strong-line ratio diagnostics (R23 and O3N2 for metallicity and N2O2 for N/O) and referencing to the effective (half-light) radius ($\rm R_e$), we find that the metallicity gradient steepens with stellar mass, lying roughly flat among galaxies with $\rm log(M_\star/M_\odot) = 9.0$ but exhibiting slopes as steep as -0.14 dex $\rm R_e^{-1}$ at $\rm log(M_\star/M_\odot) = 10.5$ (using R23, but equivalent results are obtained using O3N2). At higher masses, these slopes remain typical in the outer regions of our sample ($\rm R > 1.5 ~R_e$), but a flattening is observed in the central regions ($\rm R < 1~ R_e$). In the outer regions ($\rm R > 2.0 ~R_e$) we detect a mild flattening of the metallicity gradient in stacked profiles, although with low significance. The N/O ratio gradient provides complementary constraints on the average chemical enrichment history. Unlike the oxygen abundance, the average N/O profiles do not flatten out in the central regions of massive galaxies. The metallicity and N/O profiles both depart significantly from an exponential form, suggesting a disconnect between chemical enrichment and stellar mass surface density on local scales. In the context of inside-out growth of discs, our findings suggest that central regions of massive galaxies today have evolved to an equilibrium metallicity, while the nitrogen abundance continues to increase as a consequence of delayed secondary nucleosynthetic production.

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

We report the discovery of a mysterious giant $H_{\alpha}$ blob that is $\sim 8$ kpc away from the main MaNGA target 1-24145, one component of a dry galaxy merger, identified in the first-year SDSS-IV MaNGA data. The size of the $H_{\alpha}$ blob is $\sim$ 3-4 kpc in radius, and the $H_{\alpha}$ distribution is centrally concentrated. However, there is no optical continuum counterpart in deep broadband images reaching $\sim$26.9 mag arcsec$^{-2}$ in surface brightness. We estimate that the masses of ionized and cold gases are $3.3 \times 10^{5}$ $\rm M_{\odot}$ and $< 1.3 \times 10^{9}$ $\rm M_{\odot}$, respectively. The emission-line ratios indicate that the $H_{\alpha}$ blob is photoionized by a combination of massive young stars and AGN. Furthermore, the ionization line ratio decreases from MaNGA 1-24145 to the $H_{\alpha}$ blob, suggesting that the primary ionizing source may come from MaNGA 1-24145, likely a low-activity AGN. Possible explanations of this $H_{\alpha}$ blob include AGN outflow, the gas remnant being tidally or ram-pressure stripped from MaNGA 1-24145, or an extremely low surface brightness (LSB) galaxy. However, the stripping scenario is less favoured according to galaxy merger simulations and the morphology of the $H_{\alpha}$ blob. With the current data, we can not distinguish whether this $H_{\alpha}$ blob is ejected gas due to a past AGN outburst, or a special category of `ultra-diffuse galaxy' (UDG) interacting with MaNGA 1-24145 that further induces the gas inflow to fuel the AGN in MaNGA 1-24145.

A minimum in stellar velocity dispersion is often observed in the central regions of disc galaxies. To investigate the origin of this feature, known as a {\sigma}-drop, we analyse the stellar kinematics of a high-resolution N-body + smooth particle hydrodynamical simulation, which models the secular evolution of an unbarred disc galaxy. We compared the intrinsic mass-weighted kinematics to the recovered luminosity-weighted ones. The latter were obtained by analysing synthetic spectra produced by a new code, SYNTRA, that generates synthetic spectra by assigning a stellar population synthesis model to each star particle based on its age and metallicity. The kinematics were derived from the synthetic spectra as in real spectra to mimic the kinematic analysis of real galaxies. We found that the recovered luminosity-weighted kinematics in the centre of the simulated galaxy are biased to higher rotation velocities and lower velocity dispersions due to the presence of young stars in a thin and kinematically cool disc, and are ultimately responsible for the {\sigma}-drop.

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

Using spatially resolved spectroscopy from SDSS-IV MaNGA we have demonstrated that low ionisation emission line regions (LIERs) in local galaxies result from photoionisation by hot evolved stars, not active galactic nuclei. LIERs are ubiquitous in both quiescent galaxies and in the central regions of galaxies where star formation takes place at larger radii. We refer to these two classes of galaxies as extended LIER (eLIER) and central LIER (cLIER) galaxies respectively. cLIERs are late type galaxies located around the green valley, in the transition region between the star formation main sequence and quiescent galaxies. These galaxies display regular disc rotation in both stars and gas, although featuring a higher central stellar velocity dispersion than star forming galaxies of the same mass. cLIERs are consistent with being slowly quenched inside-out; the transformation is associated with massive bulges, pointing towards the importance of bulge growth via secular evolution. eLIERs are morphologically early types and are indistinguishable from passive galaxies devoid of line emission in terms of their stellar populations, morphology and central stellar velocity dispersion. Ionised gas in eLIERs shows both disturbed and disc-like kinematics. When a large-scale flow/rotation is observed in the gas, it is often misaligned relative to the stellar component. These features indicate that eLIERs are passive galaxies harbouring a residual cold gas component, acquired mostly via external accretion. Importantly, quiescent galaxies devoid of line emission reside in denser environments and have significantly higher satellite fraction than eLIERs. Environmental effects thus represent the likely cause for the existence of line-less galaxies on the red sequence.

We have conducted a study of extra-planar diffuse ionized gas using the first year data from the MaNGA IFU survey. We have stacked spectra from 49 edge-on, late-type galaxies as a function of distance from the midplane of the galaxy. With this technique we can detect the bright emission lines Halpha, Hbeta, [OII]3726, 3729, [OIII]5007, [NII]6549, 6584, and [SII]6717, 6731 out to about 4 kpc above the midplane. With 16 galaxies we can extend this analysis out to about 9 kpc, i.e. a distance of ~2R_e, vertically from the midplane. In the halo, the surface brightnesses of the [OII] and Halpha emission lines are comparable, unlike in the disk where Halpha dominates. When we split the sample by specific star formation rate, concentration index, and stellar mass, each subsample's emission line surface brightness profiles and ratios differ, indicating that extra-planar gas properties can vary. The emission line surface brightnesses of the gas around high specific star formation rate galaxies are higher at all distances, and the line ratios are closer to ratios characteristic of HII regions compared with low specific star formation rate galaxies. The less concentrated and lower stellar mass samples exhibit line ratios that are more like HII regions at larger distances than their more concentrated and higher stellar mass counterparts. The largest difference between different subsamples occurs when the galaxies are split by stellar mass. We additionally infer that gas far from the midplane in more massive galaxies has the highest temperatures and steepest radial temperature gradients based on their [NII]/Halpha and [OII]/Halpha ratios between the disk and the halo.

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.

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

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.

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

The nature of warm, ionized gas outside of galaxies may illuminate several key galaxy evolutionary processes. A serendipitous observation by the MaNGA survey has revealed a large, asymmetric H$\alpha$ complex with no optical counterpart that extends $\approx8"$ ($\approx6.3$ kpc) beyond the effective radius of a dusty, starbursting galaxy. This H$\alpha$ extension is approximately three times the effective radius of the host galaxy and displays a tail-like morphology. We analyze its gas-phase metallicities, gaseous kinematics, and emission-line ratios, and discuss whether this H$\alpha$ extension could be diffuse ionized gas, a gas accretion event, or something else. We find that this warm, ionized gas structure is most consistent with gas accretion through recycled wind material, which could be an important process that regulates the low-mass end of the galaxy stellar mass function.

By means of the fossil record method implemented through Pipe3D, we reconstruct the global and radial stellar mass growth histories (MGHs) of an unprecedentedly large sample of galaxies, ranging from dwarf to giant objects, from the "Mapping Nearby Galaxies at the Apache Point Observatory" survey. We confirm that the main driver of the global MGHs is mass, with more massive galaxies assembling their masses earlier (downsizing), though for a given mass, the global MGHs segregate by color, specific star formation rate (sSFR), and morphological type. From the inferred radial mean MGHs, we find that at the late evolutionary stages (or for fractions of assembled mass larger than ~ 80%), the innermost regions formed stars on average earlier than the outermost ones (inside-out). At earlier epochs, when the age resolution of the method becomes poor, the mass assembly seems to be spatially homogeneous or even in the outside-in mode, specially for the red/quiescent/early-type galaxies. The innermost MGHs are in general more regular (less scatter around the mean) than the outermost ones. For dwarf and low-mass galaxies, we do not find evidence of an outside-in formation mode; instead their radial MGHs are very diverse most of the time, with periods of outside- in and inside-out modes (or strong radial migration), suggesting this an episodic SF history. Blue/star-forming/late-type galaxies present on average a significantly more pronounced inside-out formation mode than red/quiescent/early-type galaxies, independently of mass. We discuss our results in the light of the processes of galaxy formation, quenching, and radial migration. We discuss also on the uncertainties and biases of the fossil record method and how they could affect our results.

Using kinematic maps from the Sloan Digital Sky Survey (SDSS) Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey, we reveal that the majority of low-mass quenched galaxies exhibit coherent rotation in their stellar kinematics. Our sample includes all 39 quenched low-mass galaxies observed in the first year of MaNGA. The galaxies are selected with $M_{r} > -19.1$, stellar masses $10^{9}$ M$_{\odot} < M_{\star} < 5\times10^{9}$ M$_{\odot}$, EW$_{H\alpha} <2$ \AA, and all have red colours $(u-r)>1.9$. They lie on the size-magnitude and $\sigma$-luminosity relations for previously studied dwarf galaxies. Just six ($15\pm5.7$ per cent) are found to have rotation speeds $v_{e,rot} < 15$ km s$^{-1}$ at $\sim1$ $R_{e}$, and may be dominated by pressure support at all radii. Two galaxies in our sample have kinematically distinct cores in their stellar component, likely the result of accretion. Six contain ionised gas despite not hosting ongoing star formation, and this gas is typically kinematically misaligned from their stellar component. This is the first large-scale Integral Field Unit (IFU) study of low mass galaxies selected without bias against low-density environments. Nevertheless, we find the majority of these galaxies are within $\sim1.5$ Mpc of a bright neighbour ($M_{K} < -23$; or M$_{\star} > 5\times10^{10}$ M$_{\odot}$), supporting the hypothesis that galaxy-galaxy or galaxy-group interactions quench star formation in low-mass galaxies. The local bright galaxy density for our sample is $\rho_{proj} = 8.2\pm2.0$ Mpc$^{-2}$, compared to $\rho_{proj} = 2.1\pm0.4$ Mpc$^{-2}$ for a star forming comparison sample, confirming that the quenched low mass galaxies are preferentially found in higher density environments.

We present a study of the excitation conditions and metallicity of ionized gas (Z_gas) in eight nearby barred and unbarred spiral galaxies from the VIRUS-P Exploration of Nearby Galaxies (VENGA) survey, which provides high spatial sampling and resolution (median ~ 387 pc), large coverage from the bulge to outer disc, broad wavelength range (3600-6800 Ang.), and medium spectral resolution (~ 120 km s^-1 at 5000 Ang.). Our results are: (1) We present high resolution gas excitation maps to differentiate between regions with excitation typical of Seyfert, LINER, or recent star formation. We find LINER-type excitation at large distances (3-10 kpc) from the centre, and associate this excitation with diffuse ionized gas (DIG). (2) After excluding spaxels dominated by Seyfert, LINER, or DIG, we produce maps with the best spatial resolution and sampling to date of the ionization parameter q, star formation rate, and Z_gas using common strong line diagnostics. We find that isolated barred and unbarred spirals exhibit similarly shallow Z_gas profiles from the inner kpc out to large radii (7-10 kpc or 0.5-1.0 R_25). This implies that if profiles had steeper gradients at earlier epochs, then the present-day bar is not the primary driver flattening gradients over time. This result contradicts earlier claims, but agrees with recent IFU studies. (3) The Z_gas gradients in our z ~ 0 massive spirals are markedly shallower, by ~ 0.2 dex kpc^-1, than published gradients for lensed lower mass galaxies at z ~ 1.5-2.0. Cosmologically-motivated hydrodynamical simulations best match this inferred evolution, but the match is sensitive to adopted stellar feedback prescriptions.

The MaNGA Survey (Mapping Nearby Galaxies at Apache Point Observatory) is one of three core programs in the Sloan Digital Sky Survey IV. It is obtaining integral field spectroscopy (IFS) for 10K nearby galaxies at a spectral resolution of R~2000 from 3,622-10,354A. The design of the survey is driven by a set of science requirements on the precision of estimates of the following properties: star formation rate surface density, gas metallicity, stellar population age, metallicity, and abundance ratio, and their gradients; stellar and gas kinematics; and enclosed gravitational mass as a function of radius. We describe how these science requirements set the depth of the observations and dictate sample selection. The majority of targeted galaxies are selected to ensure uniform spatial coverage in units of effective radius (Re) while maximizing spatial resolution. About 2/3 of the sample is covered out to 1.5Re (Primary sample), and 1/3 of the sample is covered to 2.5Re (Secondary sample). We describe the survey execution with details that would be useful in the design of similar future surveys. We also present statistics on the achieved data quality, specifically, the point spread function, sampling uniformity, spectral resolution, sky subtraction, and flux calibration. For our Primary sample, the median r-band signal-to-noise ratio is ~73 per 1.4A pixel for spectra stacked between 1-1.5 Re. Measurements of various galaxy properties from the first year data show that we are meeting or exceeding the defined requirements for the majority of our science goals.

Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) is an optical fiber-bundle integral-field unit (IFU) spectroscopic survey that is one of three core programs in the fourth-generation Sloan Digital Sky Survey (SDSS-IV). With a spectral coverage of 3622 - 10,354 Angstroms and an average footprint of ~ 500 arcsec^2 per IFU the scientific data products derived from MaNGA will permit exploration of the internal structure of a statistically large sample of 10,000 low redshift galaxies in unprecedented detail. Comprising 174 individually pluggable science and calibration IFUs with a near-constant data stream, MaNGA is expected to obtain ~ 100 million raw-frame spectra and ~ 10 million reduced galaxy spectra over the six-year lifetime of the survey. In this contribution, we describe the MaNGA Data Reduction Pipeline (DRP) algorithms and centralized metadata framework that produces sky-subtracted, spectrophotometrically calibrated spectra and rectified 3-D data cubes that combine individual dithered observations. For the 1390 galaxy data cubes released in Summer 2016 as part of SDSS-IV Data Release 13 (DR13), we demonstrate that the MaNGA data have nearly Poisson-limited sky subtraction shortward of ~ 8500 Angstroms and reach a typical 10-sigma limiting continuum surface brightness mu = 23.5 AB/arcsec^2 in a five arcsec diameter aperture in the g band. The wavelength calibration of the MaNGA data is accurate to 5 km/s rms, with a median spatial resolution of 2.54 arcsec FWHM (1.8 kpc at the median redshift of 0.037) and a median spectral resolution of sigma = 72 km/s.

Quiescent galaxies with little or no ongoing star formation dominate the galaxy population above $M_{*}\sim 2 \times 10^{10}~M_{\odot}$, where their numbers have increased by a factor of $\sim25$ since $z\sim2$. Once star formation is initially shut down, perhaps during the quasar phase of rapid accretion onto a supermassive black hole, an unknown mechanism must remove or heat subsequently accreted gas from stellar mass loss or mergers that would otherwise cool to form stars. Energy output from a black hole accreting at a low rate has been proposed, but observational evidence for this in the form of expanding hot gas shells is indirect and limited to radio galaxies at the centers of clusters, which are too rare to explain the vast majority of the quiescent population. Here we report bisymmetric emission features co-aligned with strong ionized gas velocity gradients from which we infer the presence of centrally-driven winds in typical quiescent galaxies that host low-luminosity active nuclei. These galaxies are surprisingly common, accounting for as much as $10\%$ of the population at $M_* \sim 2 \times 10^{10}~ M_{\odot}$. In a prototypical example, we calculate that the energy input from the galaxy's low-level active nucleus is capable of driving the observed wind, which contains sufficient mechanical energy to heat ambient, cooler gas (also detected) and thereby suppress star formation.

We study the spatially resolved excitation properties of the ionised gas in a sample of 646 galaxies using integral field spectroscopy data from SDSS-IV MaNGA. Making use of Baldwin-Philips-Terlevich diagnostic diagrams we demonstrate the ubiquitous presence of extended (kpc scale) low ionisation emission-line regions (LIERs) in both star forming and quiescent galaxies. In star forming galaxies LIER emission can be associated with diffuse ionised gas, most evident as extra-planar emission in edge-on systems. In addition, we identify two main classes of galaxies displaying LIER emission: `central LIER' (cLIER) galaxies, where central LIER emission is spatially extended, but accompanied by star formation at larger galactocentric distances, and `extended LIER' (eLIER) galaxies, where LIER emission is extended throughout the whole galaxy. In eLIER and cLIER galaxies, LIER emission is associated with radially flat, low H$\alpha$ equivalent width of line emission ($<$ 3 \AA) and stellar population indices demonstrating the lack of young stellar populations, implying that line emission follows tightly the continuum due to the underlying old stellar population. The H$\alpha$ surface brightness radial profiles are always shallower than $\rm 1/r^{2}$ and the line ratio [OIII]$\lambda$5007/[OII]$\lambda$3727,29 (a tracer of the ionisation parameter of the gas) shows a flat gradient. This combined evidence strongly supports the scenario in which LIER emission is not due to a central point source but to diffuse stellar sources, the most likely candidates being hot, evolved (post-asymptotic giant branch) stars. Shocks are observed to play a significant role in the ionisation of the gas only in rare merging and interacting systems.

We present a Bayesian approach to the redshift classification of emission-line galaxies when only a single emission line is detected spectroscopically. We consider the case of surveys for high-redshift Lyman-alpha-emitting galaxies (LAEs), which have traditionally been classified via an inferred rest-frame equivalent width (EW) greater than 20 angstrom. Our Bayesian method relies on known prior probabilities in measured emission-line luminosity functions and equivalent width distributions for the galaxy populations, and returns the probability that an object in question is an LAE given the characteristics observed. This approach will be directly relevant for the Hobby-Eberly Telescope Dark Energy Experiment (HETDEX), which seeks to classify ~10^6 emission-line galaxies into LAEs and low-redshift [O II] emitters. For a simulated HETDEX catalog with realistic measurement noise, our Bayesian method recovers 86% of LAEs missed by the traditional EW > 20 angstrom cutoff over 2 < z < 3, outperforming the EW cut in both contamination and incompleteness. This is due to the method's ability to trade off between the two types of binary classification error by adjusting the stringency of the probability requirement for classifying an observed object as an LAE. In our simulations of HETDEX, this method reduces the uncertainty in cosmological distance measurements by 14% with respect to the EW cut, equivalent to recovering 29% more cosmological information. Rather than using binary object labels, this method enables the use of classification probabilities in large-scale structure analyses. It can be applied to narrowband emission-line surveys as well as upcoming large spectroscopic surveys including Euclid and WFIRST.

We present post-cryogenic Spitzer imaging at 3.6 and 4.5 micron with the Infrared Array Camera (IRAC) of the Spitzer/HETDEX Exploratory Large-Area (SHELA) survey. SHELA covers $\sim$deg$^2$ of the Sloan Digital Sky Survey "Stripe 82" region, and falls within the footprints of the Hobby-Eberly Telescope Dark Energy Experiment (HETDEX) and the Dark Energy Survey. The HETDEX blind R $\sim$ 800 spectroscopy will produce $\sim$ 200,000 redshifts from the Lyman-$\alpha$ emission for galaxies in the range 1.9 < z < 3.5, and an additional $\sim$200,000 redshifts from the [OII] emission for galaxies at z < 0.5. When combined with deep ugriz images from the Dark Energy Camera, K-band images from NEWFIRM, and other ancillary data, the IRAC photometry from Spitzer will enable a broad range of scientific studies of the relationship between structure formation, galaxy stellar mass, halo mass, AGN, and environment over a co-moving volume of $\sim$0.5 Gpc$^3$ at 1.9 < z < 3.5. Here, we discuss the properties of the SHELA IRAC dataset, including the data acquisition, reduction, validation, and source catalogs. Our tests show the images and catalogs are 80% (50%) complete to limiting magnitudes of 22.0 (22.6) AB mag in the detection image, which is constructed from the weighted sum of the IRAC 3.6 and 4.5 micron images. The catalogs reach limiting sensitivities of 1.1 $\mu$Jy at both 3.6 and 4.5 micron (1$\sigma$, for R=2 arcsec circular apertures). As a demonstration of science, we present IRAC number counts, examples of highly temporally variable sources, and galaxy surface density profiles of rich galaxy clusters. In the spirit of Spitzer Exploratory programs we provide all images and catalogs as part of the publication.

In this review our aim is to summarize the observed properties of pseudobulges and classical bulges. We utilize an empirical approach to studying the properties of bulges in disk galaxies, and restrict our analysis to statistical proper- ties. A clear bimodality is observed in a number of properties including morphology, structural properties, star formation, gas content & stellar population, and kinematics. As well as summarizing known methods to identify pseudobulges and classical bulges we also show new results, including absorption line indices that can be used to identify different bulge types. We conclude by summarizing those properties that isolate pseudobulges from classical bulges. Our intention is to describe a practical, easy to use, list of criteria for identifying bulge types.

Mapping Nearby Galaxies at Apache Point Observatory (MaNGA), one of three core programs in the Sloan Digital Sky Survey-IV (SDSS-IV), is an integral-field spectroscopic (IFS) survey of roughly 10,000 nearby galaxies. It employs dithered observations using 17 hexagonal bundles of 2 arcsec fibers to obtain resolved spectroscopy over a wide wavelength range of 3,600-10,300A. To map the internal variations within each galaxy, we need to perform accurate {\it spectral surface photometry}, which is to calibrate the specific intensity at every spatial location sampled by each individual aperture element of the integral field unit. The calibration must correct only for the flux loss due to atmospheric throughput and the instrument response, but not for losses due to the finite geometry of the fiber aperture. This requires the use of standard star measurements to strictly separate these two flux loss factors (throughput versus geometry), a difficult challenge with standard single-fiber spectroscopy techniques due to various practical limitations. Therefore, we developed a technique for spectral surface photometry using multiple small fiber-bundles targeting standard stars simultaneously with galaxy observations. We discuss the principles of our approach and how they compare to previous efforts, and we demonstrate the precision and accuracy achieved. MaNGA's relative calibration between the wavelengths of H$\alpha$ and H$\beta$ has a root-mean-square (RMS) of 1.7%, while that between [NII] $\lambda$6583A and [OII] $\lambda$3727A has an RMS of 4.7%. Using extinction-corrected star formation rates and gas-phase metallicities as an illustration, this level of precision guarantees that flux calibration errors will be sub-dominant when estimating these quantities. The absolute calibration is better than 5% for more than 89% of MaNGA's wavelength range.

Galaxy proto-clusters at z >~ 2 provide a direct probe of the rapid mass assembly and galaxy growth of present day massive clusters. Because of the need of precise galaxy redshifts for density mapping and the prevalence of star formation before quenching, nearly all the proto-clusters known to date were confirmed by spectroscopy of galaxies with strong emission lines. Therefore, large emission-line galaxy surveys provide an efficient way to identify proto-clusters directly. Here we report the discovery of a large-scale structure at z = 2.44 in the HETDEX Pilot Survey. On a scale of a few tens of Mpc comoving, this structure shows a complex overdensity of Lya emitters (LAE), which coincides with broad-band selected galaxies in the COSMOS/UltraVISTA photometric and zCOSMOS spectroscopic catalogs, as well as overdensities of intergalactic gas revealed in the Lya absorption maps of Lee et al. (2014). We construct mock LAE catalogs to predict the cosmic evolution of this structure. We find that such an overdensity should have already broken away from the Hubble flow, and part of the structure will collapse to form a galaxy cluster with 10^14.5 +- 0.4 M_sun by z = 0. The structure contains a higher median stellar mass of broad-band selected galaxies, a boost of extended Lya nebulae, and a marginal excess of active galactic nuclei relative to the field, supporting a scenario of accelerated galaxy evolution in cluster progenitors. Based on the correlation between galaxy overdensity and the z = 0 descendant halo mass calibrated in the simulation, we predict that several hundred 1.9 < z < 3.5 proto-clusters with z = 0 mass of > 10^14.5 M_sun will be discovered in the 8.5 Gpc^3 of space surveyed by the Hobby Eberly Telescope Dark Energy Experiment.

MaNGA (Mapping Nearby Galaxies at Apache Point Observatory) is an integral-field spectroscopic survey of 10,000 nearby galaxies that is one of three core programs in the fourth-generation Sloan Digital Sky Survey (SDSS-IV). MaNGA's 17 pluggable optical fiber-bundle integral field units (IFUs) are deployed across a 3 deg field, they yield spectral coverage 3600-10,300 Ang at a typical resolution R ~ 2000, and sample the sky with 2" diameter fiber apertures with a total bundle fill factor of 56%. Observing over such a large field and range of wavelengths is particularly challenging for obtaining uniform and integral spatial coverage and resolution at all wavelengths and across each entire fiber array. Data quality is affected by the IFU construction technique, chromatic and field differential refraction, the adopted dithering strategy, and many other effects. We use numerical simulations to constrain the hardware design and observing strategy for the survey with the aim of ensuring consistent data quality that meets the survey science requirements while permitting maximum observational flexibility. We find that MaNGA science goals are best achieved with IFUs composed of a regular hexagonal grid of optical fibers with rms displacement of 5 microns or less from their nominal packing position, this goal is met by the MaNGA hardware, which achieves 3 microns rms fiber placement. We further show that MaNGA observations are best obtained in sets of three 15-minute exposures dithered along the vertices of a 1.44 arcsec equilateral triangle, these sets form the minimum observational unit, and are repeated as needed to achieve a combined signal-to-noise ratio of 5 per Angstrom per fiber in the r-band continuum at a surface brightness of 23 AB/arcsec^2. (abbrev.)

We construct a set of model spectra specifically designed to match the colours of the BOSS CMASS galaxies and to be used with photometric redshift template fitting techniques. As a basis we use a set of spectral energy distributions (SEDs) of single and composite stellar population models. These models cannot describe well the whole colour range populated by the CMASS galaxies at all redshifts, wherefore we modify them by multiplying the SEDs with $\lambda^{-\beta}$ for $\lambda>\lambda_i$ for different values of $\lambda_i$ and $\beta$. When fitting these SEDs to the colours of the CMASS sample, with a burst and dust components in superposition, we can recreate the location in colour spaces inhabited by the CMASS galaxies. From the best fitting models we select a small subset in a two-dimensional plane, whereto the galaxies were mapped by a self-organizing map. These models are used for the estimation of photometric redshifts with a Bayesian template fitting code. The photometric redshifts with the novel templates have a very small outlier rate of $0.22\,\%$, a low bias $\langle\Delta z/(1+z)\rangle=2.0\cdot10^{-3}$, and scatter of $\sigma_{68}=0.026$ in the restframe. Using our models, the galaxy colours are reproduced to a better extent with the photometric redshifts of this work than with photometric redshifts of SDSS.

MaNGA (Mapping Nearby Galaxies at Apache Point Observatory) is a 6-year SDSS-IV survey that will obtain resolved spectroscopy from 3600 $\AA$ to 10300 $\AA$ for a representative sample of over 10,000 nearby galaxies. In this paper, we derive spatially resolved stellar population properties and radial gradients by performing full spectral fitting of observed galaxy spectra from P-MaNGA, a prototype of the MaNGA instrument. These data include spectra for eighteen galaxies, covering a large range of morphological type. We derive age, metallicity, dust and stellar mass maps, and their radial gradients, using high spectral-resolution stellar population models, and assess the impact of varying the stellar library input to the models. We introduce a method to determine dust extinction which is able to give smooth stellar mass maps even in cases of high and spatially non-uniform dust attenuation. With the spectral fitting we produce detailed maps of stellar population properties which allow us to identify galactic features among this diverse sample such as spiral structure, smooth radial profiles with little azimuthal structure in spheroidal galaxies, and spatially distinct galaxy sub-components. In agreement with the literature, we find the gradients for galaxies identified as early-type to be on average flat in age, and negative (- 0.15 dex / R$_e$ ) in metallicity, whereas the gradients for late-type galaxies are on average negative in age (- 0.39 dex / R$_e$ ) and flat in metallicity. We demonstrate how different levels of data quality change the precision with which radial gradients can be measured. We show how this analysis, extended to the large numbers of MaNGA galaxies, will have the potential to shed light on galaxy structure and evolution.

We present an analysis of the data produced by the MaNGA prototype run (P-MaNGA), aiming to test how the radial gradients in recent star formation histories, as indicated by the 4000AA-break (D4000), Hdelta absorption (EW(Hd_A)) and Halpha emission (EW(Ha)) indices, can be useful for understanding disk growth and star formation cessation in local galaxies. We classify 12 galaxies observed on two P-MaNGA plates as either centrally quiescent (CQ) or centrally star-forming (CSF), according to whether D4000 measured in the central spaxel of each datacube exceeds 1.6. For each galaxy we generate both 2D maps and radial profiles of D4000, EW(Hd_A) and EW(Ha). We find that CSF galaxies generally show very weak or no radial variation in these diagnostics. In contrast, CQ galaxies present significant radial gradients, in the sense that D4000 decreases, while both EW(Hd_A) and EW(Ha) increase from the galactic center outward. The outer regions of the galaxies show greater scatter on diagrams relating the three parameters than their central parts. In particular, the clear separation between centrally-measured quiescent and star-forming galaxies in these diagnostic planes is largely filled in by the outer parts of galaxies whose global colors place them in the green valley, supporting the idea that the green valley represents a transition between blue-cloud and red-sequence phases, at least in our small sample. These results are consistent with a picture in which the cessation of star formation propagates from the center of a galaxy outwards as it moves to the red sequence.

MaNGA (Mapping Nearby Galaxies at Apache Point Observatory) is a SDSS-IV survey that will obtain spatially resolved spectroscopy from 3600 \AA\ to 10300 \AA\ for a representative sample of over 10000 nearby galaxies. In this paper we present the analysis of nebular emission line properties in 14 galaxies obtained with P-MaNGA, a prototype of the MaNGA instrument. Using spatially resolved diagnostic diagrams we find extended star formation in galaxies that are centrally dominated by Seyfert/LINER-like emission, which illustrates that galaxy characterisations based on single fibre spectra are necessarily incomplete. We observe extended LINER-like emission (up to $\rm 1 R_{e}$) in three galaxies. We make use of the $\rm EW(H \alpha)$ to argue that the observed emission is consistent with ionisation from hot evolved stars. We derive stellar population indices and demonstrate a clear correlation between $\rm D_n(4000)$ and $\rm EW(H \delta_A)$ and the position in the ionisation diagnostic diagram: resolved galactic regions which are ionised by a Seyfert/LINER-like radiation field are also devoid of recent star formation and host older and/or more metal rich stellar populations. We also detect extraplanar LINER-like emission in two highly inclined galaxies, and identify it with diffuse ionised gas. We investigate spatially resolved metallicities and find a positive correlation between metallicity and star formation rate (SFR) surface density. We study the relation between N/O vs O/H on resolved scales. We find that, at given N/O, regions within individual galaxies are spread towards lower metallicities, deviating from the sequence defined by galactic central regions as traced by Sloan $3''$ fibre spectra. We suggest that the observed dispersion can be a tracer for gas flows in galaxies: infalls of pristine gas and/or the effect of a galactic fountain.

We derive the stellar-to-halo mass relations, SHMR, of local blue and red central galaxies separately, as well as the fraction of halos hosting blue/red central galaxies. We find that: 1) the SHMR of central galaxies is segregated by color, with blue centrals having a SHMR above the one of red centrals; at logMh~12, the Ms/Mh ratio of the blue centrals is ~0.05, which is ~1.7 times larger than the value of red centrals. 2) The intrinsic scatters of the SHMRs of red and blue centrals are ~0.14 and ~0.11dex, respectively. The intrinsic scatter of the average SHMR of all central galaxies changes from ~0.20dex to ~0.14dex in the 11.3<logMh<15 range. 3) The fraction of halos hosting blue centrals at Mh=1E11Msun is 87%, but at 2x1E12Msun decays to ~20%, approaching to a few per cents at higher masses. The characteristic mass at which this fraction is the same for blue and red galaxies is Mh~7x1E11Msun. Our results suggest that the SHMR of central galaxies at large masses is shaped by halo mass quenching (likely through shock virial heating and AGN feedback), but group richness also plays an important role: central galaxies living in less dense environments quenched their star formation later or did not quench it yet. At low masses, processes that delay star formation without invoking too strong supernova-driven outflows could explain the high Ms/Mh ratios of blue centrals as compared to those of the scarce red centrals.

We describe the design, manufacture, and performance of bare-fiber integral field units (IFUs) for the SDSS-IV survey MaNGA (Mapping Nearby Galaxies at APO) on the the Sloan 2.5 m telescope at Apache Point Observatory (APO). MaNGA is a luminosity-selected integral-field spectroscopic survey of 10,000 local galaxies covering 360-1030 nm at R ~ 2200. The IFUs have hexagonal dense packing of fibers with packing regularity of 3 um (RMS), and throughput of 96+/-0.5% from 350 nm to 1 um in the lab. Their sizes range from 19 to 127 fibers (3-7 hexagonal layers) using Polymicro FBP 120:132:150 um core:clad:buffer fibers to reach a fill fraction of 56%. High throughput (and low focal-ratio degradation) is achieved by maintaining the fiber cladding and buffer intact, ensuring excellent surface polish, and applying a multi-layer AR coating of the input and output surfaces. In operations on-sky, the IFUs show only an additional 2.3% FRD-related variability in throughput despite repeated mechanical stressing during plate plugging (however other losses are present). The IFUs achieve on-sky throughput 5% above the single-fiber feeds used in SDSS-III/BOSS, attributable to equivalent performance compared to single fibers and additional gains from the AR coating. The manufacturing process is geared toward mass-production of high-multiplex systems. The low-stress process involves a precision ferrule with hexagonal inner shape designed to lead inserted fibers to settle in a dense hexagonal pattern. The ferrule inner diameter is tapered at progressively shallower angles toward its tip and the final 2 mm are straight and only a few um larger than necessary to hold the desired number of fibers. This process scales to accommodate other fiber sizes and to IFUs with substantially larger fiber count. (Abridged)

We present an overview of a new integral field spectroscopic survey called MaNGA (Mapping Nearby Galaxies at Apache Point Observatory), one of three core programs in the fourth-generation Sloan Digital Sky Survey (SDSS-IV) that began on 2014 July 1. MaNGA will investigate the internal kinematic structure and composition of gas and stars in an unprecedented sample of 10,000 nearby galaxies. We summarize essential characteristics of the instrument and survey design in the context of MaNGA's key science goals and present prototype observations to demonstrate MaNGA's scientific potential. MaNGA employs dithered observations with 17 fiber-bundle integral field units that vary in diameter from 12" (19 fibers) to 32" (127 fibers). Two dual-channel spectrographs provide simultaneous wavelength coverage over 3600-10300 A at R~2000. With a typical integration time of 3 hr, MaNGA reaches a target r-band signal-to-noise ratio of 4-8 (per A, per 2" fiber) at 23 AB mag per sq. arcsec, which is typical for the outskirts of MaNGA galaxies. Targets are selected with stellar mass greater than 1e9 Msun using SDSS-I redshifts and i-band luminosity to achieve uniform radial coverage in terms of the effective radius, an approximately flat distribution in stellar mass, and a sample spanning a wide range of environments. Analysis of our prototype observations demonstrates MaNGA's ability to probe gas ionization, shed light on recent star formation and quenching, enable dynamical modeling, decompose constituent components, and map the composition of stellar populations. MaNGA's spatially resolved spectra will enable an unprecedented study of the astrophysics of nearby galaxies in the coming 6 yr.

The Hobby-Eberly Dark Energy Experiment pilot survey identified 284 [O II] 3727 emitting galaxies in a 169 square-arcminute field of sky in the redshift range 0 < z < 0.57. This line flux limited sample provides a bridge between studies in the local universe and higher-redshift [O II] surveys. We present an analysis of the star formation rates (SFRs) of these galaxies as a function of stellar mass as determined via spectral energy distribution fitting. The [O II] emitters fall on the "main sequence" of star-forming galaxies with SFR decreasing at lower masses and redshifts. However, the slope of our relation is flatter than that found for most other samples, a result of the metallicity dependence of the [O II] star formation rate indicator. The mass specific SFR is higher for lower mass objects, supporting the idea that massive galaxies formed more quickly and efficiently than their lower mass counterparts. This is confirmed by the fact that the equivalent widths of the [O II] emission lines trend smaller with larger stellar mass. Examination of the morphologies of the [O II] emitters reveals that their star formation is not a result of mergers, and the galaxies' half-light radii do not indicate evolution of physical sizes.

We compare the H-beta line strengths of 1.90 < z < 2.35 star-forming galaxies observed with the near-IR grism of the Hubble Space Telescope with ground-based measurements of Ly-alpha from the HETDEX Pilot Survey and narrow-band imaging. By examining the line ratios of 73 galaxies, we show that most star-forming systems at this epoch have a Ly-alpha escape fraction below ~6%. We confirm this result by using stellar reddening to estimate the effective logarithmic extinction of the H-beta emission line (c_Hbeta = 0.5) and measuring both the H-beta and Ly-alpha luminosity functions in a ~ 100,000 cubic Mpc volume of space. We show that in our redshift window, the volumetric Ly-alpha escape fraction is at most 4.4+/-2.1(1.2)%, with an additional systematic ~25% uncertainty associated with our estimate of extinction. Finally, we demonstrate that the bulk of the epoch's star-forming galaxies have Ly-alpha emission line optical depths that are significantly greater than that for the underlying UV continuum. In our predominantly [O~III] 5007-selected sample of galaxies, resonant scattering must be important for the escape of Ly-alpha photons.

We present the results from a VLT/SINFONI and Keck/NIRSPEC near-infrared spectroscopic survey of 16 Lyman-alpha emitters (LAEs) at $z$ = 2.1 - 2.5 in the COSMOS and GOODS-N fields discovered from the HETDEX Pilot Survey. We detect rest-frame optical nebular lines (H$\alpha$ and/or [OIII]$\lambda$5007) for 10 of the LAEs and measure physical properties, including the star formation rate (SFR), gas-phase metallicity, gas-mass fraction, and Ly$\alpha$ velocity offset. We find that LAEs may lie below the mass-metallicity relation for continuum-selected star-forming galaxies at the same redshift. The LAEs all show velocity shifts of Ly$\alpha$ relative to the systemic redshift ranging between +85 and +296 km s$^{-1}$ with a mean of +180 km s$^{-1}$. This value is smaller than measured for continuum-selected star-forming galaxies at similar redshifts. The Ly$\alpha$ velocity offsets show a moderate correlation with the measured star formation rate (2.5$\sigma$), but no significant correlations are seen with the SFR surface density, specific SFR, stellar mass, or dynamical mass ($\lesssim$ 1.5$\sigma$). Exploring the role of dust, kinematics of the interstellar medium (ISM), and geometry on the escape of Ly$\alpha$ photons, we find no signature of selective quenching of resonantly scattered Ly$\alpha$ photons. However, we also find no evidence that a clumpy ISM is enhancing the Ly$\alpha$ equivalent width. Our results suggest that the low metallicity in LAEs may be responsible for yielding an environment with a low neutral hydrogen column density as well as less dust, easing the escape of Ly$\alpha$ photons over that in continuum-selected star-forming galaxies.

Previous studies have reported the existence of two counter-rotating stellar disks in the early-type spiral galaxy NGC7217. We have obtained high-resolution optical spectroscopic data (R ~ 9000) with the new fiber-based Integral Field Unit instrument VIRUS-W at the 2.7m telescope of the McDonald Observatory in Texas. Our analysis confirms the existence of two components. However, we find them to be co-rotating. The first component is the more luminous (~ 77% of the total light), has the higher velocity dispersion (~ 170 km/s) and rotates relatively slowly (projected $v_{max}$ = 50 km/s). The lower luminosity second component, (~ 23% of the total light), has a low velocity dispersion (~ 20 km/s) and rotates quickly (projected $v_{max}$ = 150 km/s). The difference in the kinematics of the two stellar components allows us to perform a kinematic decomposition and to measure the strengths of their Mg and Fe Lick indices separately. The rotational velocities and dispersions of the less luminous and faster component are very similar to those of the interstellar gas as measured from the [OIII] emission. Morphological evidence of active star formation in this component further suggests that NGC7217 may be in the process of (re)growing a disk inside a more massive and higher dispersion stellar halo. The kinematically cold and regular structure of the gas disk in combination with the central almost dust-free morphology allows us to compare the dynamical mass inside of the central 500pc with predictions from a stellar population analysis. We find agreement between the two if a Kroupa stellar initial mass function is assumed.

We use broadband photometry extending from the rest-frame UV to the near-IR to fit the individual spectral energy distributions (SEDs) of 63 bright (L(Ly-alpha) > 10^43 ergs/s) Ly-alpha emitting galaxies (LAEs) in the redshift range 1.9 < z < 3.6. We find that these LAEs are quite heterogeneous, with stellar masses that span over three orders of magnitude, from 7.5 < log M < 10.5. Moreover, although most LAEs have small amounts of extinction, some high-mass objects have stellar reddenings as large as E(B-V) ~0.4. Interestingly, in dusty objects the optical depths for Ly-alpha and the UV continuum are always similar, indicating that Ly-alpha photons are not undergoing many scatters before escaping their galaxy. In contrast, the ratio of optical depths in low-reddening systems can vary widely, illustrating the diverse nature of the systems. Finally, we show that in the star formation rate (SFR)-log mass diagram, our LAEs fall above the "main-sequence" defined by z ~ 3 continuum selected star-forming galaxies. In this respect, they are similar to sub-mm-selected galaxies, although most LAEs have much lower mass.

We present a quintuply lensed z ~ 6 candidate discovered in the field of the galaxy cluster RXC J2248.7-4431 (z ~ 0.348) targeted within the Cluster Lensing and Supernova survey with Hubble (CLASH) and selected in the deep HST Frontier Fields survey. Thanks to the CLASH 16-band HST imaging, we identify the quintuply lensed z ~ 6 candidate as an optical dropout in the inner region of the cluster, the brightest image having magAB=24.81+-0.02 in the f105w filter. We perform a detailed photometric analysis to verify its high-z and lensed nature. We get as photometric redshift z_phot ~ 5.9, and given the extended nature and NIR colours of the lensed images, we rule out low-z early type and galactic star contaminants. We perform a strong lensing analysis of the cluster, using 13 families of multiple lensed images identified in the HST images. Our final best model predicts the high-z quintuply lensed system with a position accuracy of 0.8''. The magnifications of the five images are between 2.2 and 8.3, which leads to a delensed UV luminosity of L_1600 ~ 0.5L*_1600 at z=6. We also estimate the UV slope from the observed NIR colours, finding a steep beta=-2.89+-0.38. We use singular and composite stellar population SEDs to fit the photometry of the hiz candidate, and we conclude that it is a young (age <300 Myr) galaxy with mass of M ~ 10^8Msol, subsolar metallicity (Z<0.2Zsol) and low dust content (AV ~ 0.2-0.4).

Survey observations of the three-dimensional locations of galaxies are a powerful approach to measure the distribution of matter in the universe, which can be used to learn about the nature of dark energy, physics of inflation, neutrino masses, etc. A competitive survey, however, requires a large volume (e.g., Vsurvey is roughly 10 Gpc3) to be covered, and thus tends to be expensive. A "sparse sampling" method offers a more affordable solution to this problem: within a survey footprint covering a given survey volume, Vsurvey, we observe only a fraction of the volume. The distribution of observed regions should be chosen such that their separation is smaller than the length scale corresponding to the wavenumber of interest. Then one can recover the power spectrum of galaxies with precision expected for a survey covering a volume of Vsurvey (rather than the volume of the sum of observed regions) with the number density of galaxies given by the total number of observed galaxies divided by Vsurvey (rather than the number density of galaxies within an observed region). We find that regularly-spaced sampling yields an unbiased power spectrum with no window function effect, and deviations from regularly-spaced sampling, which are unavoidable in realistic surveys, introduce calculable window function effects and increase the uncertainties of the recovered power spectrum. While we discuss the sparse sampling method within the context of the forthcoming Hobby-Eberly Telescope Dark Energy Experiment, the method is general and can be applied to other galaxy surveys.

We present new results on the spectrally resolved Lyman-alpha (LyA) emission of three LyA emitting field galaxies at z~2.4 with high LyA equivalent width (>100 Angstroms) and LyA luminosity (~10^43 erg/s). At 120 km/s (FWHM) spectral resolution, the prominent double-peaked LyA profile straddles the systemic velocity, where the velocity zero-point is determined from spectroscopy of the galaxies' rest-frame optical nebular emission lines. The average velocity offset from systemic of the stronger redshifted emission component for our sample is 176 km/s while the average total separation between the redshifted and main blueshifted emission components is 380 km/s. These measurements are a factor of ~2 smaller than for UV continuum-selected galaxies that show LyA in emission with lower LyA equivalent width. We compare our LyA spectra to the predicted line profiles of a spherical "expanding shell" LyA radiative transfer grid that models large-scale galaxy outflows. Specifically blueward of the systemic velocity where two galaxies show a weak, highly blueshifted (by ~1000 km/s) tertiary emission peak, the model line profiles are a relatively poor representation of the observed spectra. Since the neutral gas column density has a dominant influence over the shape of the LyA line profile, we caution against equating the observed LyA velocity offset with a physical outflow velocity, especially at lower spectral resolution where the unresolved LyA velocity offset is a convoluted function of several degenerate parameters. Referring to rest-frame ultraviolet and optical Hubble Space Telescope imaging, we find that galaxy-galaxy interactions may play an important role in inducing a starburst that results in copious LyA emission, as well as perturbing the gas distribution and velocity field which have strong influence over the LyA emission line profile.

Several occupational distributions for satellite galaxies more massive than ms~4E7 Msun around MW-sized hosts are presented and used to predict their internal dynamics. For this, a galaxy group mock catalog is constructed on the basis of (sub)halo-to-stellar mass relations constrained with the observed Galaxy Stellar Mass Function decomposed into centrals and satellites, and the 2-point correlation functions. For the MW-sized hosts we find that: 6.6% have 2 satellites in the SMC-LMC mass range; their probabilities to have 1, 2 or 3 satellites equal or larger than the LMC, the SMC or Sgr, respectively, are 0. 26, 0. 14, and 0. 14; those hosts with 3 satellites >= Sgr (as the MW) are among the most common cases, though the most and 2nd most massive satellites in these systems are SMALLER than the LMC and SMC by 0.7 and 0.8 dex, respectively. The cumulative satellite mass function is quite broad, the particular configuration of the MW being of low frequency but NOT an outlier (it is within the upper 1sig yet). For the MW satellite configuration, we constrain the halo mass to be Mh>1.4E12 Msun (1sig level). We find that the abundance of massive subhalos should agree with the abundance of massive satellites in all MW-sized hosts (there is not a massive satellite missing problem for LCDM). However, the vmax of the subhalos of satellites with ms<1E8 Msun is systematically larger than the vmax inferred from current observational studies of the MW bright dwarfs; at difference of previous works, this conclusion is based on an analysis of the overall population of MW-sized galaxies. Some pieces of evidence suggest that the issue could refer only to satellite dwarfs but not to central dwarfs; then, environmental processes associated to dwarfs inside host halos combined with SN-driven core expansion could be at the basis of the lowering of vmax.[Abridged]

We present an analysis of the luminosities and equivalent widths of the 284 z < 0.56 [O II]-emitting galaxies found in the 169 square arcmin pilot survey for the Hobby-Eberly Telescope Dark Energy Experiment (HETDEX). By combining emission-line fluxes obtained from the Mitchell spectrograph on the McDonald 2.7-m telescope with deep broadband photometry from archival data, we derive each galaxy's de-reddened [O II] 3727 luminosity and calculate its total star formation rate. We show that over the last ~5 Gyr of cosmic time there has been substantial evolution in the [O II] emission-line luminosity function, with L* decreasing by ~0.6 +/-0.2 dex in the observed function, and by ~0.9 +/-0.2 dex in the de-reddened relation. Accompanying this decline is a significant shift in the distribution of [O II] equivalent widths, with the fraction of high equivalent-width emitters declining dramatically with time. Overall, the data imply that the relative intensity of star formation within galaxies has decreased over the past ~5 Gyr, and that the star formation rate density of the universe has declined by a factor of ~2.5 between z ~ 0.5 and z ~ 0. These observations represent the first [O II]-based star formation rate density measurements in this redshift range, and foreshadow the advancements which will be generated by the main HETDEX survey.

We describe the construction of a template set of spectral energy distributions (SEDs) for the estimation of photometric redshifts of luminous red galaxies (LRGs) with a Bayesian template fitting method. By examining the color properties of several publicly available SED sets within a redshift range of 0<z<0.5 and comparing them to SDSS DR7 data, we show that only some of the investigated SEDs approximately match the colors of the LRG data throughout the redshift range, however not at the quantitative level required for precise photometric redshifts. We generate new SEDs by superposing model SEDs of composite stellar populations with a burst model, allowing both components to be reddened by dust, in order to match the data in five different redshift bins. We select a set of SEDs which represents the LRG data in color space within five redshift bins, thus defining our new SED template set for photometric redshift estimates. The results we get with the new template set and our Bayesian template fitting photometric redshift code (PhotoZ) are nearly unbiased, with a scatter of \sigma(\Delta z)=0.027 (including outliers), and a fraction of catastrophic outliers (|z_phot-z_spec|/(1+z_spec)>0.15) of 0.12%. We show that templates that optimally describe the brightest galaxies (-24.5<M_R<-22.7) indeed vary from z=0.1 to z=0.5, consistent with aging of the stellar population. Furthermore, we find that templates that optimally describe galaxies at z<0.1 strongly differ as a function of the absolute magnitude of the galaxies, indicating an increase in star formation activity for less luminous galaxies. Our findings based on the photometry of the SDSS LRGs and our SED template fitting are supported by comparison to the average SDSS LRG spectra in different luminosity and redshift bins.

We infer the local stellar-to-halo/subhalo mass relations (MRs) for central and satellite galaxies separately. We constraint this relations by using several combinations of observational data, consisting of the total galaxy stellar mass function, its decomposition into centrals and satellites, and the projected two-point correlation functions measured in different stellar mass bins. The differences among the resulting MRs are within the model-fit uncertainties (which are very small, smaller than the intrinsic scatter between galaxy and halo mass), no matter what combination of data are used. We also constrain the values of the intrinsic scatter around the central-halo (CH) and satellite-subhalo (SS) MRs assuming them to be constant: sigma_c= 0.168+-0.051 dex and sigma_s = 0.172+-0.057 dex, respectively. The CH and SS MRs are actually different, in particular when we take the subhalo mass at the present-day epoch instead of at their accretion time. We calculate several population statistics. We find that the central galaxy is not on average within the mass distribution of the most-massive satellite, even for cluster-sized halos, i.e., centrals are not a mere realization of the high-end of the satellite mass function; however for > 3x10^13 msun halos, ~15% of centrals could be. We also find that the probabilities of MW-sized halos of having N Magellanic-Clouds (MCs)-sized satellites agree well with observational measures; for a halo mass of 2x10^12 msun, the probability to have 2 MCs is 5.4%, but if we exclude those systems with satellites larger than the MCs, then the probability decreases to <2.2%.

We present the survey design, data reduction, and spectral fitting pipeline for the VIRUS-P Exploration of Nearby Galaxies (VENGA). VENGA is an integral field spectroscopic survey, which maps the disks of 30 nearby spiral galaxies. Targets span a wide range in Hubble type, star formation activity, morphology, and inclination. The VENGA data-cubes have 5.6'' FWHM spatial resolution, ~5A FWHM spectral resolution, sample the 3600A-6800A range, and cover large areas typically sampling galaxies out to ~0.7 R_25. These data-cubes can be used to produce 2D maps of the star formation rate, dust extinction, electron density, stellar population parameters, the kinematics and chemical abundances of both stars and ionized gas, and other physical quantities derived from the fitting of the stellar spectrum and the measurement of nebular emission lines. To exemplify our methods and the quality of the data, we present the VENGA data-cube on the face-on Sc galaxy NGC 628 (a.k.a. M 74). The VENGA observations of NGC 628 are described, as well as the construction of the data-cube, our spectral fitting method, and the fitting of the stellar and ionized gas velocity fields. We also propose a new method to measure the inclination of nearly face-on systems based on the matching of the stellar and gas rotation curves using asymmetric drift corrections. VENGA will measure relevant physical parameters across different environments within these galaxies, allowing a series of studies on star formation, structure assembly, stellar populations, chemical evolution, galactic feedback, nuclear activity, and the properties of the interstellar medium in massive disk galaxies.

We measure the radial profile of the 12CO(1-0) to H_2 conversion factor (Xco) in NGC 628. The H\alpha emission from the VENGA integral field spectroscopy is used to map the star formation rate surface density (\Sigma_{SFR}). We estimate the molecular gas surface density (\Sigma_{H2}) from \Sigma_{SFR} by inverting the molecular star formation law (SFL), and compare it to the CO intensity to measure Xco. We study the impact of systematic uncertainties by changing the slope of the SFL, using different SFR tracers (H\alpha vs. far-UV plus 24\mu m), and CO maps from different telescopes (single-dish and interferometers). The observed Xco profile is robust against these systematics, drops by a factor of 2 from R~7 kpc to the center of the galaxy, and is well fit by a gradient \Delta log(Xco)=0.06\pm0.02 dex kpc^-1. We study how changes in Xco follow changes in metallicity, gas density, and ionization parameter. Theoretical models show that the gradient in Xco can be explained by a combination of decreasing metallicity, and decreasing \Sigma_{H2} with radius. Photoelectric heating from the local UV radiation field appears to contribute to the decrease of Xco in higher density regions. Our results show that galactic environment plays an important role at setting the physical conditions in star forming regions, in particular the chemistry of carbon in molecular complexes, and the radiative transfer of CO emission. We caution against adopting a single Xco value when large changes in gas surface density or metallicity are present.

In this paper we present gas density, star formation rate, stellar masses, and bulge disk decompositions for a sample of 60 galaxies. Our sample is the combined sample of BIMA SONG, CARMA STING, and PdBI NUGA surveys. We study the effect of using CO-to-H_2 conversion factors that depend on the CO surface brightness, and also that of correcting star formation rates for diffuse emission from old stellar populations. We estimate that star formation rates in bulges are typically lower by 20% when correcting for diffuse emission. We find that over half of the galaxies in our sample have molecular gas surface density >100 M_sun pc^-2. We find a trend between gas density of bulges and bulge Sersic index; bulges with lower Sersic index have higher gas density. Those bulges with low Sersic index (pseudobulges) have gas fractions that are similar to that of disks. We also find that there is a strong correlation between bulges with the highest gas surface density and the galaxy being barred. However, we also find that classical bulges with low gas surface density can be barred as well. Our results suggest that understanding the connection between the central surface density of gas in disk galaxies and the presence of bars should also take into account the total gas content of the galaxy and/or bulge Sersic index. Indeed, we find that high bulge Sersic index is the best predictor of low gas density inside the bulge (not barredness of the disk). Finally, we show that when using the corrected star formation rates and gas densities, the correlation between star formation rate surface density and gas surface density of bulges is similar to that of disks.

We introduce a novel implementation of orbit-based (or Schwarzschild) modeling that allows dark matter density profiles to be calculated non-parametrically in nearby galaxies. Our models require no assumptions to be made about velocity anisotropy or the dark matter profile. The technique can be applied to any dispersion-supported stellar system, and we demonstrate its use by studying the Local Group dwarf spheroidal (dSph) galaxy Draco. We use existing kinematic data at larger radii and also present 12 new radial velocities within the central 13 pc obtained with the VIRUS-W integral field spectrograph on the 2.7m telescope at McDonald Observatory. Our non-parametric Schwarzschild models find strong evidence that the dark matter profile in Draco is cuspy for 20 < r < 700 pc. The profile for r > 20 pc is well-fit by a power law with slope \alpha=-1.0 +/- 0.2, consistent with predictions from Cold Dark Matter (CDM) simulations. Our models confirm that, despite its low baryon content relative to other dSphs, Draco lives in a massive halo.

We extend the abundance matching technique (AMT) to infer the satellite-subhalo and central-halo mass relations (MRs) of galaxies, as well as the corresponding satellite conditional mass functions (CMFs). We use the observed galaxy stellar mass function (GSMF) decomposed into centrals and satellites and the LCDM halo/subhalo mass functions as inputs. We explore the effects of defining the subhalo mass at the time of accretion (m_acc) vs. at the time of observation (m_obs). We test the standard assumption that centrals and satellites follow the same MRs, showing that this assumption leads to predictions in disagreement with observations, specially for m_obs. Instead, when the satellite-subhalo MRs are constrained following our AMT, they are always different from the central-halo MR: the smaller the stellar mass (Ms), the less massive is the subhalo of satellites as compared to the halo of centrals of the same Ms. On average, for Ms<2x10^11Msol, the dark mass of satellites decreased by 60-65% with respect to their masses at accretion time. The resulting MRs for both definitions of subhalo mass yield satellite CMFs in agreement with observations. Also, when these MRs are used in a HOD model, the predicted correlation functions agree with observations. We show that the use of m_obs leads to less uncertain MRs than m_acc, and discuss implications of the obtained satellite-subhalo MR. For example, we show that the tension between abundance and dynamics of MW satellites in LCDM gives if the slope of the GSMF faint-end slope upturns to -1.6.

We use the Marcario Low Resolution Spectrograph (LRS) at the Hobby-Eberly-Telescope (HET) to study the kinematics of pseudobulges and classical bulges in the nearby universe. We present major-axis rotational velocities, velocity dispersions, and h3 and h4 moments derived from high-resolution (sigma ~ 39 km/s) spectra for 45 S0 to Sc galaxies; for 27 of the galaxies we also present minor axis data. We combine our kinematics with bulge-to-disk decompositions. We demonstrate for the first time that purely kinematic diagnostics of the bulge dichotomy agree systematically with those based on S\'ersic index. Low S\'ersic index bulges have both increased rotational support (higher v/sigma values) and on average lower central velocity dispersions. Furthermore, we confirm that the same correlation also holds when visual morphologies are used to diagnose bulge type. The previously noted trend of photometrically flattened bulges to have shallower velocity dispersion profiles turns to be significant and systematic if the S\'ersic index is used to distinguish between pseudobulges and classical bulges. The correlation between h3 and v/sigma observed in elliptical galaxies is also observed in intermediate type galaxies, irrespective of bulge type. Finally, we present evidence for formerly undetected counter rotation in the two systems NGC 3945 and NGC 4736. Based on observations obtained with the Hobby-Eberly Telescope, which is a joint project of the University of Texas at Austin, the Pennsylvania State University, Stanford University, Ludwig-Maximilians-Universit\"at M\"unchen, and Georg-August-Universit\"at G\"ottingen.

We present results of a statistical study of the cosmic evolution of the mass dependent major-merger rate since z=1. A stellar mass limited sample of close major-merger pairs (the CPAIR sample) was selected from the archive of the COSMOS survey. Pair fractions at different redshifts derived using the CPAIR sample and a local K-band selected pair sample show no significant variations with stellar mass. The pair fraction exhibits moderately strong cosmic evolution, with the best-fitting evolutionary index m=2.2+-0.2. The best-fitting function for the merger rate implies that galaxies with stellar mass between 1E+10 -- 3E+11 M_sun have undergone 0.5 -- 1.5 major-mergers since z=1. Our results show that, for massive galaxies at z<1, major mergers involving star forming galaxies (i.e. wet and mixed mergers) can account for the formation of both ellipticals and red quiescent galaxies (RQGs). On the other hand, major mergers cannot be responsible for the formation of most low mass ellipticals and RQGs. Our quantitative estimates indicate that major mergers have significant impact on the stellar mass assembly of the most massive galaxies, but for less massive galaxies the stellar mass assembly is dominated by the star formation. Comparison with the mass dependent (U)LIRG rates suggests that the frequency of major-merger events is comparable to or higher than that of (U)LIRGs.

We present an inventory of galaxy bulge types (elliptical galaxy, classical bulge, pseudobulge, and bulgeless galaxy) in a volume-limited sample within the local 11 Mpc volume using Spitzer 3.6 micron and HST data. We find that whether counting by number, star formation rate, or stellar mass, the dominant galaxy type in the local universe has pure disk characteristics (either hosting a pseudobulge or being bulgeless). Galaxies that contain either a pseudobulge or no bulge combine to account for over 80% of the number of galaxies above a stellar mass of 10^9 M_sun. Classical bulges and elliptical galaxies account for ~1/4, and disks for ~3/4 of the stellar mass in the local 11 Mpc. About 2/3 of all star formation in the local volume takes place in galaxies with pseudobulges. Looking at the fraction of galaxies with different bulge types as a function of stellar mass, we find that the frequency of classical bulges strongly increases with stellar mass, and comes to dominate above 10^{10.5} M_sun. Galaxies with pseudobulges dominate at 10^{9.5}-10^{10.5} M_sun. Yet lower-mass galaxies are most likely to be bulgeless. If pseudobulges are not a product of mergers, then the frequency of pseudobulges in the local universe poses a challenge for galaxy evolution models.

We study the escape of Ly-alpha photons from Ly-alpha emitting galaxies (LAEs) and the overall galaxy population using a sample of 99 LAEs at 1.9<z<3.8 detected through integral-field spectroscopy of blank fields by the HETDEX Pilot Survey. For 89 LAEs showing counterparts in deep broad-band images we measure the rest-frame UV luminosity and the UV slope, which we use to estimate E(B-V) under the assumption of a constant intrinsic UV slope for LAEs. These two quantities are used to measure the dust-corrected star formation rate (SFR). A comparison between the observed Ly-alpha luminosity and that predicted by the dust-corrected SFR yields the Ly-alpha escape fraction. We also measure the Ly-alpha luminosity function. Integration of the luminosity function provides a measurement of the Ly-alpha luminosity density across our redshift range. We combine our data with that from other surveys at 0.3<z<7.7 to trace the evolution of the Ly-alpha luminosity density. We then compare it to that expected from the star-formation history of the universe in order to characterize the evolution of the Ly-alpha escape fraction of the overall galaxy population [abriged]

We present Keck/NIRSPEC spectroscopic observations of three Lyman alpha emitting galaxies (LAEs) at z ~ 2.3 discovered with the HETDEX pilot survey. We detect Halpha, [OIII], and Hbeta emission from two galaxies at z = 2.29 and 2.49, designated HPS194 and HPS256, respectively, representing the first detection of multiple rest-frame optical emission lines in galaxies at high-redshift selected on the basis of their Lyman alpha emission. The redshifts of the Lyman alpha emission from these galaxies are offset redward of the systemic redshifts by Delta_v = 162 +/- 37 (photometric) +/- 42 (systematic) km/s for HPS194, and Delta_v = 36 +/- 35 +/- 18 km/s for HPS256. An interpretation for HPS194 is that a large-scale outflow may be occurring in its interstellar medium. The emission line ratios imply that neither LAE hosts an active galactic nucleus. Using the upper limits on the [NII] emission we place meaningful constraints on the gas-phase metallicities in these two LAEs of Z < 0.17 and < 0.28 Zsol (1 sigma). Measuring the stellar masses of these objects via spectral energy distribution (SED) fitting (~ 10^10 and 6 x 10^8 Msol, respectively), we study the nature of LAEs in a mass-metallicity plane. At least one of these two LAEs appears to be more metal poor than continuum-selected star-forming galaxies at the same redshift and stellar mass, implying that objects exhibiting Lyman alpha emission may be systematically less chemically enriched than the general galaxy population. We use the SEDs of these two galaxies to show that neglecting the emission lines when fitting stellar population models to the observed photometry can result in overestimates of the population age by orders of magnitude, and the stellar mass by a factor of ~ 2. This effect is particularly important at z > 7, where similarly strong emission lines may masquerade in the photometry as a 4000 A break (abridged).

We present a catalog of emission-line galaxies selected solely by their emission-line fluxes using a wide-field integral field spectrograph. This work is partially motivated as a pilot survey for the upcoming Hobby-Eberly Telescope Dark Energy Experiment (HETDEX). We describe the observations, reductions, detections, redshift classifications, line fluxes, and counterpart information for 397 emission-line galaxies detected over 169 sq.arcmin with a 3500-5800 Ang. bandpass under 5 Ang. full-width-half-maximum (FWHM) spectral resolution. The survey's best sensitivity for unresolved objects under photometric conditions is between 4-20 E-17 erg/s/sq.cm depending on the wavelength, and Ly-alpha luminosities between 3-6 E42 erg/s are detectable. This survey method complements narrowband and color-selection techniques in the search for high redshift galaxies with its different selection properties and large volume probed. The four survey fields within the COSMOS, GOODS-N, MUNICS, and XMM-LSS areas are rich with existing, complementary data. We find 104 galaxies via their high redshift Ly-alpha emission at 1.9<z<3.8, and the majority of the remainder objects are low redshift [OII]3727 emitters at z<0.56. The classification between low and high redshift objects depends on rest frame equivalent width, as well as other indicators, where available. Based on matches to X-ray catalogs, the active galactic nuclei (AGN) fraction amongst the Ly-alpha emitters (LAEs) is 6%. We also analyze the survey's completeness and contamination properties through simulations. We find five high-z, highly-significant, resolved objects with full-width-half-maximum sizes >44 sq.arcsec which appear to be extended Ly-alpha nebulae. We also find three high-z objects with rest frame Ly-alpha equivalent widths above the level believed to be achievable with normal star formation, EW(rest)>240 Ang.

We dissect giant Sc-Scd galaxies with Hubble Space Telescope photometry and Hobby-Eberly Telescope spectroscopy. We use HET's High Resolution Spectrograph (resolution = 15,000) to measure stellar velocity dispersions in the nuclear star clusters and pseudobulges of the pure-disk galaxies M33, M101, NGC 3338, NGC 3810, NGC 6503, and NGC 6946. We conclude: (1) Upper limits on the masses of any supermassive black holes are MBH <= (2.6+-0.5) * 10**6 M_Sun in M101 and MBH <= (2.0+-0.6) * 10**6 M_Sun in NGC 6503. (2) HST photometry shows that the above galaxies contain tiny pseudobulges that make up <~ 3 % of the stellar mass but no classical bulges. We inventory a sphere of radius 8 Mpc centered on our Galaxy to see whether giant, pure-disk galaxies are common or rare. In this volume, 11 of 19 galaxies with rotation velocity > 150 km/s show no evidence for a classical bulge. Four may contain small classical bulges that contribute 5-12% of the galaxy light. Only 4 of the 19 giant galaxies are ellipticals or have classical bulges that contribute 1/3 of the galaxy light. So pure-disk galaxies are far from rare. It is hard to understand how they could form as the quiescent tail of a distribution of merger histories. Recognition of pseudobulges makes the biggest problem with cold dark matter galaxy formation more acute: How can hierarchical clustering make so many giant, pure-disk galaxies with no evidence for merger-built bulges? This problem depends strongly on environment: the Virgo cluster is not a puzzle, because >2/3 of its stellar mass is in merger remnants.

The increasing abundance of passive "red-sequence" galaxies since z=1-2 is mirrored by a coincident rise in the number of galaxies with spheroidal morphologies. In this paper, however, we show that in detail the correspondence between galaxy morphology and color is not perfect, providing insight into the physical origin of this evolution. Using the COSMOS survey, we study a significant population of red sequence galaxies with disk-like morphologies. These passive disks typically have Sa-Sb morphological types with large bulges, but they are not confined to dense environments. They represent nearly one-half of all red-sequence galaxies and dominate at lower masses (log Mstar < 10) where they are increasingly disk-dominated. As a function of time, the abundance of passive disks with log Mstar < 11 increases, but not as fast as red-sequence spheroidals in the same mass range. At higher mass, the passive disk population has declined since z~1, likely because they transform into spheroidals. We estimate that as much as 60% of galaxies transitioning onto the red sequence evolve through a passive disk phase. The origin of passive disks therefore has broad implications for understanding how star formation shuts down. Because passive disks tend to be more bulge-dominated than their star-forming counterparts, a simple fading of blue disks does not fully explain their origin. We explore several more sophisticated explanations, including environmental effects, internal stabilization, and disk regrowth during gas-rich mergers. While previous work has sought to explain color and morphological transformations with a single process, these observations open the way to new insight by highlighting the fact that galaxy evolution may actually proceed through several separate stages.

We investigate scaling relations of bulges using bulge-disk decompositions at 3.6 micron and present bulge classifications for 173 E-Sd galaxies within 20 Mpc. Pseudobulges and classical bulges are identified using Sersic index, HST morphology, and star formation activity (traced by 8 micron emission). In the near-IR pseudobulges have n_b<2 and classical bulges have n_b>2, as found in the optical. Sersic index and morphology are essentially equivalent properties for bulge classification purposes. We confirm, using a much more robust sample, that the Sersic index of pseudobulges is uncorrelated with other bulge structural properties, unlike for classical bulges and elliptical galaxies. Also, the half-light radius of pseudobulges is not correlated with any other bulge property. We also find a new correlation between surface brightness and pseudobulge luminosity; pseudobulges become more luminous as they become more dense. Classical bulges follow the well known scaling relations between surface brightness, luminosity and half-light radius that are established by elliptical galaxies. We show that those pseudobulges (as indicated by Sersic index and nuclear morphology) that have low specific star formation rates are very similar to models of galaxies in which both a pseudobulge and classical bulge exist. Therefore, pseudobulge identification that relies only on structural indicators is incomplete. Our results, especially those on scaling relations, imply that pseudobulges are very different types of objects than elliptical galaxies.

The baryon content of high-density regions in the universe is relevant to two critical unanswered questions: the workings of nurture effects on galaxies and the whereabouts of the missing baryons. In this paper, we analyze the distribution of dark matter and semianalytical galaxies in the Millennium Simulation to investigate these problems. Applying the same density field reconstruction schemes as used for the overall matter distribution to the matter locked in halos we study the mass contribution of halos to the total mass budget at various background field densities, i.e., the conditional halo mass function. In this context, we present a simple fitting formula for the cumulative mass function accurate to ~ 5% for halo masses between 10^{10} and 10^{15}Msol/h. We find that in dense environments the halo mass function becomes top heavy and present corresponding fitting formulae for different redshifts. We demonstrate that the major fraction of matter in high-density fields is associated with galaxy groups. Since current X-ray surveys are able to nearly recover the universal baryon fraction within groups, our results indicate that the major part of the so-far undetected warm-hot intergalactic medium resides in low-density regions at low temperatures. Similarly, we show that the differences in galaxy mass functions with environment seen in observed and simulated data stem predominantly from differences in the mass distribution of halos. In particular, the hump in the galaxy mass function is associated with the central group galaxies, and the bimodality observed in the galaxy mass function is therefore interpreted as that of central galaxies versus satellites.

VENGA is a large-scale extragalactic IFU survey, which maps the bulges, bars and large parts of the outer disks of 32 nearby normal spiral galaxies. The targets are chosen to span a wide range in Hubble types, star formation activities, morphologies, and inclinations, at the same time of having vast available multi-wavelength coverage from the far-UV to the mid-IR, and available CO and 21cm mapping. The VENGA dataset will provide 2D maps of the SFR, stellar and gas kinematics, chemical abundances, ISM density and ionization states, dust extinction and stellar populations for these 32 galaxies. The uniqueness of the VIRUS-P large field of view permits these large-scale mappings to be performed. VENGA will allow us to correlate all these important quantities throughout the different environments present in galactic disks, allowing the conduction of a large number of studies in star formation, structure assembly, galactic feedback and ISM in galaxies.

We present a new analysis of stellar mass functions (MF) in the COSMOS field to fainter limits than has been previously probed to z~1. Neither the total nor the passive or star-forming MF can be well fit with a single Schechter function once one probes below 3e9 Msun. We observe a dip or plateau at masses ~1e10 Msun, and an upturn towards a steep faint-end slope of -1.7 at lower mass at any z<1. This bimodal nature of the MF is not solely a result of the blue/red dichotomy. The blue MF is by itself bimodal at z~1. This suggests a new dichotomy in galaxy formation that predates the appearance of the red sequence. We propose two interpretations for this bimodality. If the gas fraction increases towards lower mass, galaxies with M_baryon~1e10 Msun would shift to lower stellar masses, creating the observed dip. This would indicate a change in star formation efficiency, perhaps linked to supernovae feedback becoming much more efficient. Therefore, we investigate whether the dip is present in the baryonic (stars+gas) MF. Alternatively, the dip could be created by an enhancement of the galaxy assembly rate at ~1e11 Msun, a phenomenon that naturally arises if the baryon fraction peaks at M_halo ~1e12 Msun. In this scenario, galaxies occupying the bump around M* would be identified with central galaxies and the second fainter component having a steep faint-end slope with satellites. While the dip is apparent in the total MF at any z, it appears to shift from the blue to red population, likely as a result of transforming high-mass blue galaxies into red ones. At the same time, we detect a drastic upturn in the number of low-mass red galaxies. Their increase with time reflects a decrease in the number of blue systems and so we tentatively associate them with satellite dwarf galaxies that have undergone quenching.

The Cosmic Evolution Survey (COSMOS) allows for the first time a highly significant census of environments and structures up to redshift one, as well as a full morphological description of the galaxy population. In this paper we present a study aimed to constrain the evolution, in the redshift range 0.2 < z < 1.2, of the mass content of different morphological types and its dependence on the environmental density. We use a deep multicolor catalog, covering an area of ~0.7 square degrees inside the COSMOS field, with accurate photometric redshifts (i < 26.5 and dz/(z+1) ~ 0.035). We estimate galaxy stellar masses by fitting the multi-color photometry to a grid of composite stellar population models. We quantitatively describe the galaxy morphology by fitting PSF convolved Sersic profiles to the galaxy surface brightness distributions down to F814 = 24 mag for a sample of 41300 objects. We confirm an evolution of the morphological mix with redshift: the higher the redshift the more disk-dominated galaxies become important. We find that the morphological mix is a function of the local comoving density: the morphology density relation extends up to the highest redshift explored. The stellar mass function of disk-dominated galaxies is consistent with being constant with redshift. Conversely, the stellar mass function of bulge-dominated systems shows a decline in normalization with redshift. Such different behaviors of late-types and early-types stellar mass functions naturally set the redshift evolution of the transition mass. ABRIDGED

We study star formation rates (SFR) and stellar masses in bulges of nearby disk galaxies, using SFRs and stellar masses derived from Spitzer and GALEX data. At present day SFR the median pseudobulge could have grown the present day stellar mass in 8 Gyr. In almost all galaxies in our sample the specific SFR (SFR per unit stellar mass) of the bulge is higher than that of the outer disk, suggesting that almost all galaxies are increasing their B/T through internal star formation. In pseudobulges, SFR density correlates, positvely, with mass density, this is consistent with that stellar mass being formed by moderate, extended star fromation. As well, SFR density and stellar mass of pseudobulges are shown to be correlated with the stellar mass of the outer disk. Classical bulges have the lowest specific SFR implying a growth times that are longer than a Hubble time. We identify a class of bulges that have nuclear morphology similar to pseudobulges, significantly lower specific SFR than pseudobulges, and are closer to classical bulges in structural parameter correlations. Our results are consistent with a scenario in which bulge growth via internal star formation is a natural, and near ubiquitous phenomenon in disk galaxies. Some disk galaxies with out a large classical bulge, over long periods of extended star formation are able to growth a pseudobulge. In this sense, galaxies with pseudobulges may very well be bulgeless (or "quasi-bulgeless") galaxies, and galaxies with classical bulges are galaxies in which both internal evolution and hierarchical merging are responsible for the bulge mass by fractions that vary from galaxy-to-galaxy. [Abridged]

The Hobby-Eberly Telescope Dark Energy Experiment (HETDEX) will outfit the 10 m HET with a new wide field and an array of 150 integral-field spectrographs to survey a 420 sq. deg. area in the north Galactic cap. Each fiber-coupled unit spectrograph will cover 350-550 nm, simultaneously. This instrument, called VIRUS, will produce ~34,000 spectra per exposure, and will open up the emission-line universe to large surveys for the first time. The survey will detect 0.8 million Lyman-alpha emitting (LAE) galaxies with 1.9<z<3.5 and more than a million [OII] emitting galaxies with z<0.5. The 3-D map of LAE galaxies in 9 cubic Gpc volume will be used to measure the expansion history at this early epoch using baryonic acoustic oscillations and the shape of the power spectrum. The aim of HETDEX is to provide a direct detection of dark energy at z~3. The measurement will constrain the evolution of dark energy and will also provide 0.1%-level accuracy on the curvature of the Universe, ten times better than current. The prototype of the VIRUS unit spectrograph (VIRUS-P) is a powerful instrument in its own right. Used on the McDonald 2.7 m, it covers the largest area of any integral field spectrograph, and reaches wavelengths down to 340 nm. VIRUS-P is being used for a pilot survey to better measure the properties of LAE galaxies in support of HETDEX. We report initial results from this survey.

In this paper we study the properties of pseudobulges (bulges that appear similar to disk galaxies) and classical bulges (bulges which appear similar to E-type galaxies) in bulge-disk decompositions. We show that the distribution of bulge Sersic indices is bimodal, and this bimodality correlates with the morphology of the bulge. Pseudobulges have n<2 and classical bulges have n>2 with little overlap. Also, pseudobulges do not follow the correlations of Sersic index with structural parameters or the photometric projections of the fundamental plane in the same way as classical bulges and elliptical galaxies do. We find that pseudobulges are systematically flatter than classical bulges and thus more disk-like in both their morphology and shape. We do not find significant differences between different bulge morphologies that we are collectively calling pseudobulges (nuclear spirals, nuclear rings, nuclear bars, and nuclear patchiness) appear to behave similarly in all parameter correlations. In Sersic index, flattening, and bulge-to-total ratio, the distinction appears to be between classical bulges and pseudobulges, not between different pseudobulge morphologies. The Sersic index of the pseudobulges does not correlate with B/T, in contrast to classical bulges. Also, the half-light radius of the pseudobulge correlates with the scale length of the disk; this is not the case for classical bulges. The correlation of Sersic index and scale lengths with bulge morphology suggests that secular evolution is creating pseudobulges with low-Sersic index, and that other processes (e.g. major mergers) are responsible for the higher Sersic index in classical bulges and elliptical galaxies.

We derive the disk masses of 18 spiral galaxies of different luminosity and Hubble Type, both by mass modelling their rotation curves and by fitting their SED with spectro-photometric models. The good agreement of the estimates obtained from these two different methods allows us to quantify the reliability of their performance and to derive very accurate stellar mass-to-light ratio vs color (and stellar mass) relationships.

We present a formalism to infer the presence of merging by comparing the time derivative of the observed galaxy stellar mass function (MF) to the change of the MF expected from the star formation rate (SFR) in galaxies as a function of mass and time. We present the SFR in as a function of stellar mass and time spanning 9<logM<12 and 0<z<5. We show that at z>=3 the average SFR, is a power law of stellar mass (SFR~M^0.6). The average SFR in the most massive objects at this redshift is 100-500 Msun/yr. At z~3, the SFR starts to drop at the high mass end. As z decreases further, the SFR drops at progressively lower masses (downsizing), dropping most rapidly for high mass (logM>11) galaxies. The mass at which the SFR starts to deviate from the power-law form (break mass) progresses smoothly from logM~13 at z~5 to logM~10.9 at z~0.5. The break mass evolves with redshift as M(z)=2.7x10^10 (1+z)^2.1. We directly observe a relationship between SFH and mass. More massive galaxies have steeper and earlier onsets of SF, their SFR peaks earlier and higher, and the following decay has a shorter e-folding time. The SFR observed in high mass galaxies at z~4 is sufficient to explain their rapid increase in number density. Within large uncertainties, at most 0.8 major mergers per Gyr are consistent with the high-z data, yet enough to transform most high mass objects into ellipticals contemporaneously with their major star formation episode. In contrast, at z<1.5 and at logM>11, mergers contribute 0.1-0.2 Gyr^-1 to the relative increase in number density (~1 major merger per massive object at 1.5>z>0). At 10<logM<11, galaxies are being preferably destroyed in mergers at high z, while at later times the change in their numbers turns positive. This shows the top-down buildup of the red sequence suggested by other observations.

The global colors of galaxies have recently been shown to follow bimodal distributions. Galaxies separate into a ``red sequence'', populated prototypically by early-type galaxies, and a ``blue cloud'', whose typical objects are late-type disk galaxies. Intermediate-type (Sa-Sbc) galaxies populate both regions. It has been suggested that this bimodality reflects the two-component nature of disk-bulge galaxies. However, it has now been established that there are two types of bulges: ``classical bulges'' that are dynamically hot systems resembling (little) ellipticals, and ``pseudobulges'', dynamically cold, flattened, disk-like structures that could not have formed via violent relaxation. Therefore thee question is whether at types Sa-Sbc, where both bulge types are found, the red-blue dichotomy separates galaxies at some value of disk-to-bulge ratio, $B/T$, or, whether it separates galaxies of different bulge type, irrespective of their $B/T$. We identify classical bulges and pseudobulges morphologically with HST images in a sample of nearby galaxies. Detailed surface photometry reveals that: (1) The red -- blue dichotomy is a function of bulge type: at the same $B/T$, pseudobulges are in globally blue galaxies and classical bulges are in globally red galaxies. (2) Bulge type also predicts where the galaxy lies in other (bimodal) global structural parameters. (3) Hence, the red -- blue dichotomy is not due to decreasing bulge prominence alone, and the bulge type of a galaxy carries significance for the galaxy's evolutionary history ... (Abridged)

(Abridged) We present B, R, and I-band selected galaxy catalogues based on the Munich Near-Infrared Cluster Survey (MUNICS) which, together with the K-selected sample, serve as an important probe of galaxy evolution in the redshift range 0 < z < 2. Furthermore, used in comparison they are ideally suited to study selection effects. The construction of the B, R, and I-selected photometric catalogues, containing ~9000, ~9000, and ~6000 galaxies, respectively, is described in detail. The catalogues reach 50% completeness limits for point sources of B ~ 24.5mag, R ~ 23.5mag, and I ~ 22.5mag and cover an area of about 0.3 square degrees. Photometric redshifts are derived for all galaxies with an accuracy of dz/(1+z) ~ 0.057. We investigate the influence of selection band and environment on the specific star formation rate (SSFR). We find that K-band selection indeed comes close to selection in stellar mass, while B-band selection purely selects galaxies in star formation rate. We use a galaxy group catalogue constructed on the K-band selected MUNICS sample to study possible differences of the SSFR between the field and the group environment, finding a marginally lower average SSFR in groups as compared to the field, especially at lower redshifts. The field-galaxy luminosity function in the B and R band as derived from the R-selected sample evolves out to z ~ 2 in the sense that the characteristic luminosity increases but the number density decreases. This effect is smaller at longer rest-frame wavelengths and gets more pronounced at shorter wavelengths. Parametrising the redshift evolution of the Schechter parameters as M*(z) = M*(0) + a ln(1+z) and Phi*(z) = Phi*(0) (1+z)^b we find evolutionary parameters a ~ -2.1 and b ~ -2.5 for the B band, and a ~ -1.4 and b ~ -1.8 for the R band.

We present the morphological analysis based on HST-NIC2 (0.075 arcsec/pixel) images in the F160W filter of a sample of 9 massive field (> 10^{11} M_\odot) galaxies spectroscopically classified as early-types at 1.2<z<1.7. Our analysis shows that all of them are bulge dominated systems. In particular, 6 of them are well fitted by a de Vaucouleurs profile (n=4) suggesting that they can be considered pure elliptical galaxies. The remaining 3 galaxies are better fitted by a Sersic profile with index 1.9<n<2.3 suggesting that a disk-like component could contribute up to 30% to the total light of these galaxies. We derived the effective radius R_e and the mean surface brightness <mu_e> within R_e of our galaxies and we compared them with those of early-types at lower redshifts. We find that the surface brightness <mu_e> of our galaxies should get fainter by 2.5 mag from z~1.5 to z~0 to match the surface brightness of the local ellipticals with comparable R_e, i.e. the local Kormendy relation. Luminosity evolution without morphological changes can only explain half of this effect, as the maximum dimming expected for an elliptical galaxy is ~1.6 mag in this redshift range. Thus, other parameters, possibly structural, may undergo evolution and play an important role in reconciling models and observations. Hypothesizing an evolution of the effective radius of galaxies we find that R_e should increase by a factor 1.5 from z~1.5 to z~0.

Spectra of galaxies contain an enormous amount of information about the relative mixture of ages and metallicities of constituent stars. We present a comprehensive study designed to extract the maximum information from spectra of data quality typical in large galaxy surveys. We test fitting techniques using the Bruzual-Charlot 2003 high resolution simple stellar population synthesis models to simultaneously estimate the ages and metallicities of 101 globular clusters in M31 and the Magellanic Clouds. The clusters cover a wide range of ages and metallicities, 4 Myr to 20 Gyr and -1.6 < [Fe/H] < +0.3, estimated by other methods in the literature. We compare results from model fits to both the spectra and photometry and find that fits to continuum-normalized (CN) spectra over the entire range available, typically 350-1000 nm for this sample, provides the best results. For clusters older than 1 Gyr we agree with literature ages to 0.16 dex (35%) and [Fe/H] to 0.12 dex. For younger clusters we agree with literature ages to 0.3 dex (63%), but cannot constrain the metallicity. It is particularly important to use the entire CN spectrum to avoid problems with model continua for young objects and to break age-metallicity degeneracies of broadband photometry. Our required S/N = 15-30 for 20% age uncertainties and 30-55 for 10% uncertainties over spectral resolutions of 5-25 Angstroms. This technique should work well for the age-metallicity parameter space expected for early-type galaxies at z~1, although individual galaxy spectral S/N may require the coaddition of a few like objects.

The optical rest-frame sizes of 10 of the most massive (~5x10^{11}h_{70}^{-2}M_sun) galaxies found in the near-infrared MUNICS survey at 1.2<z<1.7 are analysed. Sizes were estimated both in the J and K' filters. These massive galaxies are at least a factor of 4_{-1.0}^{+1.9} (+-1 sigma) smaller in the rest-frame V-band than local counterparts of the same stellar mass. Consequently, the stellar mass density of these objects is (at least) 60 times larger than massive ellipticals today. Although the stellar populations of these objects are passively fading, their structural properties are rapidly changing since that redshift. This observational fact disagrees with a scenario where the more massive and passive galaxies are fully assembled at z~1.4 (i.e. a monolithic scenario) and points towards a dry merger scenario as the responsible mechanism for the subsequent evolution of these galaxies.

We study the evolution of the stellar mass density for the separate families of bulge-dominated and disk-dominated galaxies over the redshift range 0.25 < z < 1.15. We derive quantitative morphology for a statistically significant galaxy sample of 1645 objects selected from the FORS Deep and the GOODS-S Fields. We find that the morphological mix evolves monotonically with time: the higher the redshift, the more disk systems dominate the total mass content. At redshift about 1, massive objects (M_stellar > 7E10 M_solar) host about half of the mass contained in objects of similar mass in the local universe. The contribution from early and late type galaxies to the mass budget at z about 1 is nearly equal. We show that in situ star formation is not sufficient to explain the changing mass budget. Moreover we find that the star formation rate per unit stellar mass of massive galaxies increases with redshift only for the intermediate and early morphological types, while it stays nearly constant for late-type objects. This suggests that merging and/or frequent accretion of small mass objects has a key role in the shaping of the Hubble sequence as we observe it now, and also in decreasing the star formation activity of the bulge-dominated descendants of massive disk galaxies.

We present the redshift evolution of the restframe galaxy luminosity function (LF) in the red r', i', and z' bands as derived from the FORS Deep Field (FDF). Using the deep and homogeneous I-band selected dataset of the FDF we are able to follow the red LFs over the redshift range 0.5 < z < 3.5. The results are based on photometric redshifts for 5558 galaxies derived from the photometry in 9 filters achieving an accuracy of Delta z / (z_spec+1) ~ 0.03 with only ~ 1 % outliers. Because of the depth of the FDF we can give relatively tight constraints on the faint-end slope alpha of the LF: The faint-end of the red LFs does not show a large redshift evolution and is compatible within 1 sigma to 2 sigma with a constant slope over the redshift range 0.5 < z < 2.0. Moreover, the slopes in r', i', and z' are very similar with a best fitting value of alpha= -1.33 +- 0.03 for the combined bands. There is a clear trend of alpha to steepen with increasing wavelength: alpha_(UV & u')=-1.07 +- 0.04 -> alpha_(g' & B)=-1.25 +- 0.03 -> alpha_(r' & i' & z')=-1.33 +- 0.03. We show that the wavelength dependence of the LF slope can be explained by the relative contribution of different SED-type LFs to the overall LF, as different SED types dominate the LF in the blue and red bands. Furthermore we also derive and analyze the luminosity density evolution of the different SED types up to z ~ 2. Based on the FDF data, we find only a mild brightening of M_star and decrease of phi_star with increasing redshift. Therefore, from <z> ~ 0.5 to <z> \~ 3 the characteristic luminosity increases by ~0.8, ~0.4 and ~0.4 magnitudes in the r', i', and z' bands, respectively. Simultaneously the characteristic density decreases by about 40 % in all analyzed wavebands. [abridged]

We explore the build-up of stellar mass in galaxies over a wide redshift range 0.4 < z < 5.0 by studying the evolution of the specific star formation rate (SSFR), defined as the star formation rate per unit stellar mass, as a function of stellar mass and age. Our work is based on a combined sample of ~ 9000 galaxies from the FORS Deep Field and the GOODS-S field, providing high statistical accuracy and relative insensitivity against cosmic variance. As at lower redshifts, we find that lower-mass galaxies show higher SSFRs than higher mass galaxies, although highly obscured galaxies remain undetected in our sample. Furthermore, the highest mass galaxies contain the oldest stellar populations at all redshifts, in principle agreement with the existence of evolved, massive galaxies at 1 < z < 3. It is remarkable, however, that this trend continues to very high redshifts of z ~ 4. We also show that with increasing redshift the SSFR for massive galaxies increases by a factor of ~ 10, reaching the era of their formation at z ~ 2 and beyond. These findings can be interpreted as evidence for an early epoch of star formation in the most massive galaxies, and ongoing star-formation activity in lower mass galaxies.

We present results from a study to determine how star formation contributes to galaxy growth since redshift z=1.5. Using galaxies from the MUnich Near-Infrared Cluster Survey (MUNICS) and the FORS Deep Field (FDF), we investigate the specific star formation rate (SSFR, star formation rate [SFR] per unit galaxy stellar mass) as a function of galaxy stellar mass and redshift. We test the compatibility of our results with a sample drawn from a larger volume using the Sloan Digital Sky Survey. We find that the SSFR decreases as galaxy stellar mass increases, suggesting that star formation contributes more to the growth of low-mass galaxies than high-mass galaxies at all redshifts in this study. We also find a ridge in the SSFR that runs parallel to lines of constant SFR and decreases by a factor of 10 from z=1 to today, matching the results of the evolution in SFR density seen in the ``Lilly-Madau'' diagram. The ridge evolves independently of galaxy stellar mass to a particular turnover mass at the high mass end. Galaxies above the turnover mass show a sharp decrease in SSFR compared to the average at that epoch, and the turnover mass increases with redshift.

We present the analysis of 10 massive early-type galaxies at $z\sim1.5$. They have been identified by means of a near-IR low resolution spectroscopic follow-up of a complete sample of 36 bright (K' $<$ 18.5) Extremely Red Objects (EROs, R-K'$>$ 5) selected from the Munich Near-IR ClusterSurvey (MUNICS; Drory et al. 2001). The low resolution near-IR spectra constrain their redshift at $1.2<z<1.7$, implying absolute magnitudes M$_{K'}<-26.0$ and stellar masses well in excess of 10$^{11}$ M$_\odot$. Under the hypothesis of pure passive evolution from $z\sim1.5$ to $z=0$, in the local universe they would have luminosities L$_K\ge2.5$L$^*$. Thus, they are the high-z counterparts of the local old massive (10$^{11}-10^{12}$ M$_\odot$) early-type galaxies populating the bright end of the local luminosity function of galaxies. The comparison of their spectro-photometric properties with a grid of synthetic models suggests that the stellar populations in more than half of the sample are about $\sim$3-5 Gyr old and 1-2 Gyr old in the remaining part. These ages imply formation redshift $z_{f} > 2$ for all the galaxies and $z_{f} \geq 4$ for the oldest ones. The comparison of the 4000\AA break and of the overall spectral shape of the average spectrum of the 10 galaxies at $z\sim1.5$ with those of their local counterpartsconfirms that field massive early-type galaxies formed the bulk of their stellar mass at $2<z<4$, most likely over a short ($<$ 1 Gyr) star formation time scale, consistently with the results derived from the analysis of their individual spectro-photometric properties.

We present a study to determine how star formation contributes to galaxy growth since z=1.5 over five decades in galaxy stellar mass. We investigate the specific star formation rate (SSFR; star formation rate [SFR] per unit galaxy stellar mass) as a function of galaxy stellar mass and redshift. A sample of 175 K-band selected galaxies from the MUnich Near-Infrared Cluster Survey spectroscopic dataset provide intermediate to high mass galaxies (mostly M* > 10^10 Msun) to z=1. The FORS Deep Field provides 168 low mass galaxies (mostly M* < 10^10 Msun) to z=1.5. We use a Sloan Digital Sky Survey galaxy sample to test the compatibility of our results with data drawn from a larger volume. We find that at all redshifts, the SSFR decreases with increasing galaxy stellar mass suggesting that star formation contributes more to the growth of low mass galaxies than to the growth of high mass galaxies, and that high mass galaxies formed the bulk of their stellar content before z=1. At each epoch we find a ridge in SSFR versus stellar mass that is parallel to lines of constant SFR and evolves independently of galaxy stellar mass to a particular turnover mass. Galaxies above this turnover mass show a sharp decrease in the SFR compared to the average at each epoch and the turnover mass increases with redshift. The SFR along the SSFR ridge decreases by roughly a factor of 10, from 10 Msun/yr at z=1.5 to 1 Msun/yr at z=0. High mass galaxies could sustain the observed rates of star formation over the 10 Gyr observed, but low mass galaxies likely undergo episodic starbursts.

We present a measurement of the evolution of the stellar mass function (MF) of galaxies and the evolution of the total stellar mass density at 0<z<5. We use deep multicolor data in the Fors Deep Field (FDF; I-selected reaching I_AB=26.8) and the GOODS-S/CDFS region (K-selected reaching K_AB=25.4) to estimate stellar masses based on fits to composite stellar population models for 5557 and 3367 sources, respectively. The MF of objects from the GOODS-S sample is very similar to that of the FDF. Near-IR selected surveys hence detect the more massive objects of the same principal population as do I-selected surveys. We find that the most massive galaxies harbor the oldest stellar populations at all redshifts. At low z, our MF follows the local MF very well, extending the local MF down to 10^8 Msun. The faint end slope is consistent with the local value of alpha~1.1 at least up to z~1.5. Our MF also agrees very well with the MUNICS and K20 results at z<2. The MF seems to evolve in a regular way at least up to z~2 with the normalization decreasing by 50% to z=1 and by 70% to z=2. Objects having M>10^10 Msun which are the likely progenitors of todays L* galaxies are found in much smaller numbers above z=2. However, we note that massive galaxies with M>10^11 Msun are present even to the largest redshift we probe. Beyond z=2 the evolution of the mass function becomes more rapid. We find that the total stellar mass density at z=1 is 50% of the local value. At z=2, 25% of the local mass density is assembled, and at z=3 and z=5 we find that at least 15% and 5% of the mass in stars is in place, respectively. The number density of galaxies with M>10^11 Msun evolves very similarly to the evolution at lower masses. It decreases by 0.4 dex to z=1, by 0.6 dex to z=2, and by 1 dex to z=4.

We spectroscopically identified 7 massive evolved galaxies with magnitudes 17.8<K<18.4 at 1.3<z<1.7 over an area of ~160 arcmin^2 of the MUNICS survey. Their rest-frame K-band absolute magnitudes are -26.8<M$_K<-26.1 (5L*<L_K<10L*) and the resulting stellar masses are in the range 3-6.5x10^{11} M_sun. The analysis we performed unambiguously shows the early-type nature of their spectra. The 7 massive evolved galaxies account for a comoving density of \~5.5(+-2) x 10^{-5} Mpc^{-3} at <z>~1.5, a factor 1.5 lower than the density (8.4(+-1) x 10^{-5} Mpc^{-3}) of early-types with comparable masses at z=0. The incompleteness (30%) of our spectroscopic observations accounts for this discrepancy. Thus, our data do not support a decrease of the comoving density of early-type galaxies with masses comparable to the most massive ones in the local Universe up to z~1.7. This suggests that massive evolved galaxies do not play an important role in the evolution of the mass density outlined by recent surveys in this redshift range, evolution which instead has to be ascribed to the accretion of the stellar mass in late-type galaxies. Finally, the presence of such massive evolved galaxies at these redshifts suggests that the assembly of massive spheroids has taken place at z>2 supporting a high efficiency in the accretion of the stellar mass in massive halos in the early Universe.

We investigate the contribution of star formation to the growth of stellar mass in galaxies over the redshift range 0.5 < z < 1.1 by studying the redshift evolution of the specific star formation rate (SSFR), defined as the star formation rate per unit stellar mass. We use an I-band selected sample of 6180 field galaxies from the Munich Near-Infrared Cluster Survey (MUNICS) with spectroscopically calibrated photometric redshifts. The SSFR decreases with stellar mass at all redshifts. The low SSFRs of massive galaxies indicates that star formation does not significantly change their stellar mass over this redshift range: The majority of massive galaxies have assembled the bulk of their mass before redshift unity. Furthermore, these highest mass galaxies contain the oldest stellar populations at all redshifts. The line of maximum SSFR runs parallel to lines of constant star formation rate. With increasing redshift, the maximum SFR is generally increasing for all stellar masses, from SFR ~ 5 M_sun/yr at z = 0.5 to SFR ~ 10 M_sun/yr at z = 1.1. We also show that the large SSFRs of low-mass galaxies cannot be sustained over extended periods of time. Finally, our results do not require a substantial contribution of merging to the growth of stellar mass in massive galaxies over the redshift range probed. We note that highly obscured galaxies which remain undetected in our sample do not affect these findings for the bulk of the field galaxy population.

We present XMM-Newton data centered on one of the MUNICS Near Infrared Cluster Survey fields (S2F1) and we discuss the X-ray properties of the 6 X-ray emitting EROs found. For one of them we have already obtained the redshift using near-infrared spectroscopic data, while for the remaining 5 EROs the analysis is based on photometric redshifts. We find evidence for the presence of an X-ray obscured QSO in at least 5 out of the 6 X-ray emitting EROs. For these 5 objects we derive intrinsic (2-10 keV) luminosities in excess of 10^44 erg/s and intrinsic column densities higher than 10^22 cm^-2. These values have been obtained through a basic X-ray spectral analysis for the three brightest sources and through the analysis of the hardness ratios for the remaining two. All of these 5 X-ray emitting EROs appear extended in the optical/near-infrared bands indicating that the host galaxy emission dominates at these wavelengths. This suggests that the hosted AGNs are likely to be absorbed also in the optical/near-infrared bands: i.e. they are likely X-ray obscured possible type 2 QSOs. For the remaining ERO the presence of an AGN is suggested both by its high 0.5-2 keV luminosity (L(0.5-2 keV)~10^43 erg/s) and by its X-ray-to-optical flux ratio. In this case the quality of the present data prevents us from placing firm constraints on the AGN type hosted. Finally, the near-IR spectrum obtained for one of the 6 EROs classifies the host galaxy as an elliptical at z~1.7 with a stellar mass well in excess of 10^11 M_\odot. This result corroborates the possible link between the QSO activity and the formation of massive spheroids.

We measure the star formation rate (SFR) as a function of redshift z up to z \~4.5, based on B, I and (I+B) selected galaxy catalogues from the FORS Deep Field (FDF) and the K-selected catalogue from the GOODS-South field. Distances are computed from spectroscopically calibrated photometric redshifts accurate to (Delta_z / (z_spec+1)) ~0.03 for the FDF and ~0.056 for the GOODS-South field. The SFRs are derived from the luminosities at 1500 Angstroem. We find that the total SFR estimates derived from B, I and I+B catalogues agree very well ($\lsim 0.1$ dex) while the SFR from the K catalogue is lower by ~0.2 dex. We show that the latter is solely due to the lower star-forming activity of K-selected intermediate and low luminosity (L<L_*) galaxies. The SFR of bright (L>L_*) galaxies is independent of the selection band, i.e. the same for B, I, (I+B), and K-selected galaxy samples. At all redshifts, luminous galaxies (L>L_*) contribute only ~1/3 to the total SFR. There is no evidence for significant cosmic variance between the SFRs in the FDF and GOODs-South field, ~0.1 dex, consistent with theoretical expectations. The SFRs derived here are in excellent agreement with previous measurements provided we assume the same faint-end slope of the luminosity function as previous works (alpha ~ -1.6). However, our deep FDF data indicate a shallower slope of alpha=-1.07, implying a SFR lower by ~0.3 dex. We find the SFR to be roughly constant up to z ~4 and then to decline slowly beyond, if dust extinctions are assumed to be constant with redshift.

The purpose of this letter is to compare the quality of different methods for estimating stellar masses of galaxies. We compare the results of (a) fitting stellar population synthesis models to broad band colors from SDSS and 2MASS, (b) the analysis of spectroscopic features of SDSS galaxies (Kauffmann et al. 2003), and, (c) a simple dynamical mass estimate based on SDSS velocity dispersions and effective radii. Knowing that all three methods can have significant biases, a comparison can help to establish their (relative) reliability. In this way, one can also probe the quality of the observationally cheap broadband color mass estimators for galaxies at higher redshift. Generally, masses based on broad-band colors and spectroscopic features agree reasonably well, with a rms scatter of only ~ 0.25 dex over almost 4 decades in mass. However, as may be expected, systematic differences do exist and have an amplitude of ~ 0.15 dex, corrleting with Halpha emission strength. Interestingly, masses from broad-band color fitting are in better agreement with dynamical masses than masses based on the analysis of spectroscopic features. In addition, the differences between the latter and the dynamical masses correlate with Halpha equivalent width, while this much less the case for the broad-band masses. We conclude that broad band color mass estimators, provided they are based on a large enough wavelength coverage and use an appropriate range of ages, metallicities and dust extinctions, can yield fairly reliable stellar masses for galaxies. This is a very encouraging result as such mass estimates are very likely the only ones available at significant redshifts for some time to come.

We present a measurement of the evolution of the stellar mass function in four redshift bins at 0.4 < z < 1.2 using a sample of more than 5000 K-selected galaxies drawn from the MUNICS dataset. Our data cover the stellar mass range 10^10 < M/Msun < 10^12. We derive K-band mass-to-light ratios by fitting a grid of composite stellar population models of varying star formation history, age, and dust extinction to BVRIJK photometry. We discuss the evolution of the average mass-to-light ratio as a function of galaxy stellar mass in the K-band and in the B-band. We compare our stellar mass function at z > 0 to estimates obtained similarly at z=0. We find that the mass-to-light ratios in the K-band decline with redshift. This decline is similar for all stellar masses above $10^10 Msun. Lower mass galaxies have lower mass-to-light ratios at all redshifts. The stellar mass function evolves significantly to z = 1.2. The total normalization decreases by a factor of ~2, the characteristic mass (the knee) shifts towards lower masses and the bright end therefore steepens with redshift. The amount of number density evolution is a strong function of stellar mass, with more massive systems showing faster evolution than less massive systems. We discuss the total stellar mass density of the universe and compare our results to the values from the literature both at lower and higher redshift. We find that the stellar mass density at z~1 is roughly 50% of the local value. Our results imply that the mass assembly of galaxies continues well after $z \sim 1$. Our data favor a scenario in which the growth of the most massive galaxies is dominated by accretion and merging rather than star formation which plays a larger role in the growth of less massive systems.

We use the very deep and homogeneous I-band selected dataset of the FORS Deep Field (FDF) to trace the evolution of the luminosity function over the redshift range 0.5 < z < 5.0. We show that the FDF I-band selection down to I(AB)=26.8 misses of the order of 10 % of the galaxies that would be detected in a K-band selected survey with magnitude limit K(AB)=26.3 (like FIRES). Photometric redshifts for 5558 galaxies are estimated based on the photometry in 9 filters (U, B, Gunn g, R, I, SDSS z, J, K and a special filter centered at 834 nm). A comparison with 362 spectroscopic redshifts shows that the achieved accuracy of the photometric redshifts is (Delta z / (z_spec+1)) < 0.03 with only ~ 1 % outliers. This allows us to derive luminosity functions with a reliability similar to spectroscopic surveys. In addition, the luminosity functions can be traced to objects of lower luminosity which generally are not accessible to spectroscopy. We investigate the evolution of the luminosity functions evaluated in the restframe UV (1500 Angstroem and 2800 Angstroem), u', B, and g' bands. Comparison with results from the literature shows the reliability of the derived luminosity functions. Out to redshifts of z ~ 2.5 the data are consistent with a slope of the luminosity function approximately constant with redshift, at a value of -1.07 +- 0.04 in the UV (1500 Angstroem, 2800 Angstroem) as well as u', and -1.25 +- 0.03 in the blue (g', B). We do not see evidence for a very steep slope (alpha < -1.6) in the UV at z ~ 3.0 and z ~ 4.0 favoured by other authors. There may be a tendency for the faint-end slope to become shallower with increasing redshift but the effect is marginal. We find a brightening of M_star and a decrease of Phi_star with redshift for all analyzed wavelengths. [abridged]

We summarize the results obtained from the MUNICS K-band selected Galaxy survey thus far. MUNICS is a wide-area, medium-deep, photometric and spectroscopic survey selected in the K band, targeting randomly-selected high Galactic latitude fields. It covers an area of roughly one square degree in the K and J bands with optical imaging in the I, R, V, and B bands in 0.5 square degrees. The MUNICS photometric survey is complemented by spectroscopic follow-up observations down to limits of K<17.5 (wide area) and K<19.5 (smaller area). We have obtained 593 redshifts to this date. Here, we present results concerning the evolution of the K-band luminosity function to z ~ 1, both from the full photometric redshift sample and from the spectroscopic sample alone. We also report on new results concerning the evolution of the stellar mass function to z ~ 1. We fit stellar population synthesis models to our multicolor photometry to obtain M/L values for each source. We detect significant evolution in the stellar mass function and we find that more massive systems evolve faster in number density than less massive systems. We also measure the evolution of the total stellar mass density of the universe and find that about half the present day stellar mass formed since z ~ 1.

Here we present the results obtained from the analysis of 75 ksec of XMM-Newton observations of a sample of EROs selected from one MUNICS field (K'<19.5 mag). We find 6 EROs with a X-ray counterpart down to a 2--10 keV flux limit of ~10^{-15} cgs. For all of them the X-ray--to--optical flux ratios and the 2--10 keV luminosities suggest the presence of AGN. In particular, a complete X-ray spectral analysis shows that high luminosity, obscured AGNs (i.e. QSO2 candidates) are present in 3 of them.

We present the preliminary results of an on-going near-IR low resolution spectroscopic follow-up of a complete sample of 30 bright (K'<18.5) EROs selected over 360 arcmin^2 from the MUNICH survey. Among the 13 EROs observed so far, 7 turned out to be early-type galaxies at 1.2<z<1.45 with stellar masses in excess of 2x10^11 M\odot. These 7 early-types account for 40% of the local population of bright ellipticals providing a lower limit to the density of massive early-types at z>1.

We are carrying on a near-IR very low resolution spectroscopic follow-up in parallel with XMM-Newton observations of a complete sample of ~30 bright (K'<18.5) Extremely Red Objects (EROs) selected over an area of 360 arcmin^2 of the MUNICS survey. We here present the preliminary results of the spectroscopic and X-ray data analysis.

(Abriged) We present a measurement of the evolution of the rest-frame K-band luminosity function to z ~ 1.2 using a sample of more than 5000 K-selected galaxies drawn from the MUNICS dataset. Distances and absolute K-band magnitudes are derived using photometric redshifts from spectral energy distribution fits to BVRIJK photometry. These are calibrated using >500 spectroscopic redshifts. We obtain redshift estimates having a rms scatter of 0.055 and no mean bias. We use Monte-Carlo simulations to investigate the influence of the errors in distance associated with photometric redshifts on our ability to reconstruct the shape of the luminosity function. Finally, we construct the rest-frame K-band LF in four redshift bins spanning 0.4<z<1.2 and compare our results to the local luminosity function. We discuss and apply two different estimators to derive likely values for the evolution of the number density, Phi*, and characteristic luminosity, M*, with redshift. While the first estimator relies on the value of the luminosity function binned in magnitude and redshift, the second estimator uses the individually measured {M,z} pairs alone. In both cases we obtain a mild decrease in number density by \~ 25% to z=1 accompanied by brightening of the galaxy population by 0.5 to 0.7 mag. These results are fully consistent with an analogous analysis using only the spectroscopic MUNICS sample. The total K-band luminosity density is found to scale as dlog(rho_L)/dz = 0.24. We discuss possible sources of systematic errors and their influence on our parameter estimates.

We present spectroscopic follow-up observations of galaxies from the Munich Near-Infrared Cluster Survey (MUNICS). MUNICS is a wide-field medium-deep K'-band selected survey covering 1 square degree in the near-infrared K' and J pass-bands, and 0.35 square degrees in I, R, V, and - recently completed - B. The spectroscopic sample comprises observations of objects down to a limit of K' < 17.5 in five survey fields (0.17 square degrees in total), and a sparsely selected deeper sample (K' < 19.0) constructed in one of the survey patches (0.03 square degrees). Here we describe the selection procedure of objects for spectroscopic observations, the observations themselves, the data reduction, and the construction of the spectroscopic catalogue containing roughly 500 galaxies with secure redshifts. Furthermore we discuss global properties of the sample like its distribution in colour-redshift space, the accuracy of redshift determination, and the completeness function of the data. We derive the rest-frame K'-band luminosity function of galaxies at median redshifts of z = 0.2, z = 0.4, and z = 0.7. We find evidence for mild evolution of magnitudes (-0.70 mag) and number densities (35 per cent) to redshift one. Furthermore, we present the rest-frame J-band luminosity function of galaxies at these redshifts, the first determination of this quantity at higher redshifts, with a behaviour similar to the K-band luminosity function.

We present the results of TNG near-IR low resolution spectroscopy of two (S7F5-254 and S7F5-45) sources belonging to a complete sample of 15 EROs with K'<18 and R-K'>5 selected from the MUNICS Survey. Both the spectra show a sharp drop in the continuum which can be ascribed only to the Balmer break. This places them at 1.2<z<1.5. Their rest-frame z=1.2 K-band absolute magnitude is M_K ~ -26.6 (L ~ 7L*). The comparison of the spectra and the photometric data with a grid of synthetic template spectra provides a redshift z=1.22 for S7F5-254 and z=1.46 for S7F5-45. The resulting lower limits to their stellar mass are M_{stars}=6x10^{11} M_{sun} and M_{stars}=4x10^{11} M_{sun}. The minimum age of the last burst of star formation in S7F5-254 is 3.5 Gyr while it is 0.5 Gyr in S7F5-45 implying a minimum formation redshift z_f>3.5 and z_f>2 for the two EROs respectively.

The FORS Deep Field project is a multi-colour, multi-object spectroscopic investigation of an approx. 7 times 7 region near the south galactic pole based mostly on observations carried out with the FORS instruments attached to the VLT telescopes. It includes the QSO Q 0103-260 (z = 3.36). The goal of this study is to improve our understanding of the formation and evolution of galaxies in the young Universe. In this paper the field selection, the photometric observations, and the data reduction are described. The source detection and photometry of objects in the FORS Deep Field is discussed in detail. A combined B and I selected UBgRIJKs photometric catalog of 8753 objects in the FDF is presented and its properties are briefly discussed. The formal 50% completeness limits for point sources, derived from the co-added images, are 25.64, 27.69, 26.86, 26.68, 26.37, 23.60 and 21.57 in U, B, g, R, I, J and Ks (Vega-system), respectively. A comparison of the number counts in the FORS Deep Field to those derived in other deep field surveys shows very good agreement.

The Munich Near-IR Cluster Survey (MUNICS) is a wide-area, medium-deep, photometric survey selected in the K' band. The project's main scientific aims are the identification of galaxy clusters up to redshifts of unity and the selection of a large sample of field early-type galaxies up to z < 1.5 for evolutionary studies. We created a Large Scale Structure catalog, using a new structure finding technique specialized for photometric datasets, that we developed on the basis of a friends-of-friends algorithm. We tested the plausibility of the resulting galaxy group and cluster catalog with the help of Color-Magnitude Diagrams (CMD), as well as a likelihood- and Voronoi-approach.

Using a sample of 57 VLT FORS spectra in the redshift range 1.37<z<3.40 (selected mainly from the FORS Deep Field survey) and a comparison sample with 36 IUE spectra of local (z ~ 0) starburst galaxies we derive CIV and SiIV equivalent width values and estimate metallicities of starburst galaxies as a function of redshift. Assuming that a calibration of the CIV equivalent widths in terms of the metallicity based on the local sample of starburst galaxies is applicable to high-z objects, we find a significant increase of the average metallicities from about 0.16 Z_sun at the cosmic epoch corresponding to z ~ 3.2 to about 0.42 Z_sun at z ~ 2.3. A significant further increase in metallicity during later epochs cannot be detected in our data. Compared to the local starburst galaxies our high-redshift objects tend to be overluminous for a fixed metallicity. Our observational results are in good agreement with published observational data by other authors and with theoretical predictions of the cosmic chemical evolution.

We present the results of completeness simulations for the detection of point sources as well as redshifted elliptical and spiral galaxies in the K'-band images of the Munich Near-Infrared Cluster Survey (MUNICS). The main focus of this work is to quantify the selection effects introduced by threshold-based object detection algorithms used in deep imaging surveys. Therefore, we simulate objects obeying the well-known scaling relations between effective radius and central surface brightness, both for de Vaucouleurs and exponential profiles. The results of these simulations, while presented for the MUNICS project, are applicable in a much wider context to deep optical and near-infrared selected samples. We investigate the detection probability as well as the reliability for recovering the true total magnitude with Kron-like (adaptive) aperture photometry. The results are compared to the predictions of the visibility theory of Disney and Phillipps in terms of the detection rate and the lost-light fraction. Additionally, the effects attributable to seeing are explored. The results show a bias against detecting high-redshifted massive elliptical galaxies in comparison to disk galaxies with exponential profiles, and that the measurements of the total magnitudes for intrinsically bright elliptical galaxies are systematically too faint. Disk galaxies, in contrast, show no significant offset in the magnitude measurement of luminous objects. Finally we present an analytic formula to predict the completeness of point-sources using only basic image parameters.

A near-IR low-resolution spectroscopic follow-up of a complete sample of K<18.5 EROs has started at the 3.5 m Telescopio Nazionale Galileo (TNG, La Palma, Canary Islands) by using a prism disperser AMICI. The sample has been selected from the Munich Near-IR Cluster Survey (MUNICS). The project aims at revealing massive galaxies at z>1.2 and at estimating their spatial density to put constraints on the galaxy formation and evolution models. In this proceeding, the preliminary analysis of the first spectra obtained is presented.

We derive the number density evolution of massive field galaxies in the redshift range 0.4 < z < 1.2 using the K-band selected field galaxy sample from the Munich Near-IR Cluster Survey (MUNICS). We rely on spectroscopically calibrated photometric redshifts to determine distances and absolute magnitudes in the rest-frame K-band. To assign mass-to-light ratios, we use two different approaches. First, we use an approach which maximizes the stellar mass for any K-band luminosity at any redshift. We take the mass-to-light ratio of a Simple Stellar Population (SSP) which is as old as the universe at the galaxy's redshift as a likely upper limit. Second, we assign each galaxy a mass-to-light ratio by fitting the galaxy's colours against a grid of composite stellar population models and taking their M/L. We compute the number density of galaxies more massive than 2 x 10^10 h^-2 Msun, 5 x 10^10 h^-2 Msun, and 1 x 10^11 h^-2 Msun, finding that the integrated stellar mass function is roughly constant for the lowest mass limit and that it decreases with redshift by a factor of ~ 3 and by a factor of ~ 6 for the two higher mass limits, respectively. This finding is in qualitative agreement with models of hierarchical galaxy formation, which predict that the number density of ~ M* objects is fairly constant while it decreases faster for more massive systems over the redshift range our data probe.

We present the B-band Tully-Fisher relation (TFR) of 60 late-type galaxies with redshifts 0.1-1. The galaxies were selected from the FORS Deep Field with a limiting magnitude of R=23. Spatially resolved rotation curves were derived from spectra obtained with FORS2 at the VLT. High-mass galaxies with v_max>=150 km/s show little evolution, whereas the least massive systems in our sample are brighter by about 1-2 mag compared to their local counterparts. For the entire distant sample, the TFR slope is flatter than for local field galaxies (-5.77+-0.45 versus -7.92+-0.18). Thus, we find evidence for evolution of the slope of the TFR with redshift on the 3-sigma level. This is still true when we subdivide the sample into three redshift bins. We speculate that the flatter tilt of our sample is caused by the evolution of luminosities and an additional population of blue galaxies at z>=0.2. The mass dependence of the TFR evolution also leads to variations for different galaxy types in magnitude-limited samples, suggesting that selection effects can account for the discrepant results of previous TFR studies on the luminosity evolution of late-type galaxies.

We derive the number density evolution of massive field galaxies in the redshift range 0.4 < z < 1.2 using the K-band selected field galaxy sample from the Munich Near-IR Cluster Survey (MUNICS). We rely on spectroscopically calibrated photometric redshifts to determine distances and absolute magnitudes in the rest-frame K-band. To assign mass-to-light ratios, we use an approach which maximizes the stellar mass for any K-band luminosity at any redshift. We take the mass-to-light ratio, M/L_K, of a Simple Stellar Population (SSP) which is as old as the universe at the galaxy's redshift as a likely upper limit. This is the most extreme case of pure luminosity evolution and in a more realistic model M/L_K will probably decrease faster with redshift due to increased star formation. We compute the number density of galaxies more massive than 2 10^10 h^-2 solar masses, 5 10^10 h^-2 solar masses, and 1 10^11 h^-2 solar masses, finding that the integrated stellar mass function is roughly constant for the lowest mass limit and that it decreases with redshift by a factor of roughly 3 and by a factor of roughly 6 for the two higher mass limits, respectively. This finding is in qualitative agreement with models of hierarchical galaxy formation, which predict that the number density of ~ M* objects is fairly constant while it decreases faster for more massive systems over the redshift range our data probe.

We present the results of radio, optical and near-infrared observations of the field of TOC J0233.3+3021, a cluster of milliJansky radiosources from the TexOx Cluster survey. In an observation of this field with the Ryle Telescope (RT) at 15 GHz, we measure a decrement in the cosmic microwave background (CMB) of $-675 \pm 95 \mu$Jy on the RT's $\approx$ 0.65 k$\lambda$ baseline. Using optical and infrared imaging with the McDonald 2.7-m Smith Reflector, Calar Alto 3.5-m telescope and UKIRT, we identify the host galaxies of five of the radiosources and measure magnitudes of $R \approx 24$, $J \approx 20$, $K \approx 18$. The CMB decrement is consistent with the Sunyaev-Zel'dovich (SZ) effect of a massive cluster of galaxies, which if modelled as a spherical King profile of core radius $\theta_C = 20^{\prime\prime}$ has a central temperature decrement of $900 \mu$K. The magnitudes and colours of the galaxies are consistent with those of old ellipticals at $z \sim 1$. We therefore conclude that TOC J0233.3+3021 is a massive, high redshift cluster. These observations add to the growing evidence for a significant population of massive clusters at high redshift, and demonstrate the effectiveness of combining searches for AGN `signposts' to clusters with the redshift-independence of the SZ effect.

The Munich Near-IR Cluster Survey (MUNICS) is a wide-area, medium-deep, photometric survey selected in the K' band. It covers an area of roughly one square degree in the K' and J near-IR pass-bands. The survey area consists of 16 6' x 6' fields targeted at QSOs with redshifts 0.5 < z < 2 and 7 28' x 13' stripes targeted at `random' high Galactic latitude fields. Ten of the QSO fields were additionally imaged in R and I, and 0.6 square degrees of the randomly selected fields were also imaged in the V, R, and I bands. The resulting object catalogues were strictly selected in K', having a limiting magnitude (50 per cent completeness) of K' ~ 19.5 mag and J ~ 21 mag, sufficiently deep to detect passively evolving systems up to a redshift of z ~ 1.5 and luminosity of 0.5 L*. The optical data reach a depth of roughly R ~ 23.5 mag. The project's main scientific aims are the identification of galaxy clusters at redshifts around unity and the selection of a large sample of field early-type galaxies at 0 < z < 1.5 for evolutionary studies. In this paper - the first in a series - we describe the survey's concept, the selection of the survey fields, the near-IR and optical imaging and data reduction, object extraction, and the construction of photometric catalogues. Finally, we show the J-K' vs. K' colour-magnitude diagramme and the R-J vs. J-K', V-I vs. J-K', and V-I vs. V-R colour-colour diagrammes for MUNICS objects, together with stellar population-synthesis models for different star-formation histories, and conclude that the data set presented is suitable for extracting a catalogue of massive field galaxies in the redshift range 0.5 < z < 1.5 for evolutionary studies and follow-up observations.

Dedicating a major fraction of its guaranteed time, the FORS consortium established a FORS Deep Field which contains a known QSO at z = 3.36. It was imaged in UBgRIz with FORS at the VLT as well as in J and Ks with the NTT. Covering an area 6-8 times larger as the HDFs but with similar depth in the optical it is one of the largest deep fields up to date to investigate i) galaxy evolution in the field from present up to z $\sim$ 5, ii) the galaxy distribution in the line of sight to the QSO, iii) the high-z QSO environment and iv) the galaxy-galaxy lensing signal in such a large field. In this presentation a status report of the FORS Deep Field project is given. In particular, the field selection, the imaging results (number counts, photometric redshifts etc.) and the first spectroscopic results are presented.

The Munich Near-IR Cluster Survey (MUNICS) is a K' selected survey covering 1 square degree in the K' and J NIR bands with complementary optical photometry in the V, R, and I bands covering a subarea of 0.35 square degrees. The 3-sigma limiting magnitude is 19.5 in K'. The main goals of the project are the identification of clusters of galaxies at redshifts 0.6<z<1.0 and the study of the evolution of the early-type field population at similar redshifts. Here we present first results regarding color distributions and the surface densities of EROs as well as photometric redshifts and a first clustering analysis of the sample.

We present high-speed long term photometric observations from 1992 to 1996 (UBVRI), as well as recently obtained (1997) spectrophotometric and high resolution spectroscopic studies of the bright galactic supersoft X-ray source RXJ0019.8+2156. Our photometry reveals a highly variable object. The shortest observed quasi-periodic variations are humps with a period of approx 1.8h. The timings of the main minimum are not exactly regular and occur on average around phase phi=0.0. With our new data set we calculated an update of the orbital period. Our simultaneous spectroscopic and spectrophotometric studies from 1997 give an even more detailed insight: blue and red shifted satellite lines at HeII 4686 and at the strong Balmer lines, which are interpreted as high velocity outflows, can clearly be detected. We also present for the first time an analysis of this supersoft source by means of Doppler tomography. The lines of HeII (4542,4686), H_a and H_b show no emission from an accretion disk. These emission lines originate in regions within the binary system with very low velocities. The spatial origin of the emitting material is not quite clear and will be investigated in a follow-up paper.

We describe our preliminary results on a search for high-redshift ($z >$ 0.5) galaxy clusters using near-infrared photometry obtained with the Omega camera at the 3.5m telescope at Calar Alto.

The MUNICS project is an ongoing imaging survey designed to cover 3 sq. degrees in V,R,I,J,K'. We describe here partial results of the project concerning the clustering properties of K' < 19.5 galaxies in scales of 3.6'' to 63.0'' over an area of $\sim 800 arcmin^2$. We present K' data for a sample of 20 fields, five of which contain $z > 0.5$ radio-loud quasars with steep spectra, eight contain $z > 0.5$ radio-loud quasars with flat spectra and seven are high-galactic latitude fields with no quasars in them. The two-point angular correlation function for the total sample shows significant clustering at $\sim 5 \sigma$ level of K'=19.5 galaxies. The correlation angle of the galaxies is $\theta_0 = 1.7'' \pm 0.4''$ for K' < 19 mag and $\theta_0 = 1.0'' \pm 0.2''$ for K' < 19.5 mag. When the correlation functions for the subsamples are considered, the mean $\omega(\theta)$ amplitude of the fields which contain steep-spectra $z > 0.5$ radio-loud quasars is determined to be $\sim 2.0 - 2.5$ that of the high-galactic latitude fields.