Source author record

Matthew A. Bershady

Matthew A. Bershady appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

49works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

49 published item(s)

preprint2022arXiv

Optimum telescope focal ratios for microlens-to-fiber coupled integral field spectrographs

We describe the optimum telescope focal ratio for a two-element, three surface, telecentric image-transfer microlens-to-fiber coupled integral field unit within the constraints imposed by micro-optics fabrication and optical aberrations. We create a generalized analytical description of the micro-optics optical parameters from first principals. We find that the optical performance, including all aberrations, of a design constrained by an analytic model considering only spherical aberration and diffraction matches within +/-4% of a design optimized by ray-tracing software such as Zemax. The analytical model does not require any compromise on the available clear aperture; about 90% mechanical aperture of a hexagonal microlens is available for light collection. The optimum telescope f-ratio for a 200μm core fiber fed at f/3.5 is between f/7 and f/12. We find the optimum telescope focal ratio changes as a function of fiber core diameter and fiber input beam speed. A telescope focal ratio of f/8 would support the largest range of fiber diameters (100 to 500μm) and fiber injection speeds (between f/3 and f/5). The optimization of telescope and lenslet-coupled fibers is relevant for the design of high-efficiency dedicated survey telescopes, and for retro-fitting existing facilities via introducing focal macro-optics to match the instrument input requirements.

preprint2022arXiv

PetroFit: A Python Package for Computing Petrosian Radii and Fitting Galaxy Light Profiles

PetroFit is an open-source Python package, based on Astropy and Photutils, that can calculate Petrosian profiles and fit galaxy images. It offers end-to-end tools for making accurate photometric measurements, estimating morphological properties, and fitting 2D models to galaxy images. Petrosian metric radii can be used for model parameter estimation and aperture photometry to provide accurate total fluxes. Correction tools are provided for improving Petrosian radii estimates affected by galaxy morphology. PetroFit also provides tools for sampling Astropy-based models (including custom profiles and multi-component models) onto image grids and enables PSF convolution to account for the effects of seeing. These capabilities provide a robust means of modeling and fitting galaxy light profiles. We have made the PetroFit package publicly available on GitHub (PetroFit/petrofit) and PyPi (pip install petrofit).

preprint2022arXiv

SDSS-IV MaNGA: Exploring the local scaling relations for N/O

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

preprint2022arXiv

SDSS-IV MaNGA: Understanding Ionized Gas Turbulence using Integral Field Spectroscopy of 4500 Star-Forming Disk Galaxies

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

preprint2022arXiv

Slit Mask Integral Field Units for the Southern African Large Telescope

Two fibre integral field units (IFU) are being built in the SAAO fibre-lab for the Robert Stobie Spectrograph's visible arm and the future red arm. Each IFU sits in its own slit-mask cassette and is referred to as a slit-mask IFU (SMI). They will be available some time in 2022. The smaller, 200 micron fibre IFU has 309 X 0.9 arcsec diameter spatial elements covering an elongated hexagonal footprint of 18 X 23 arcsec. The larger, 400 micron fibre IFU has 178 X 1.8 arcsec diameter spatial elements covering an on-sky area of 21 X 44 arcsec. In both cases there are two groups of 13 fibres offset by roughly 50 arcsec on either side of the primary array to sample sky. The 1.8 and 0.9 arcsec spatial resolution SMIs provide median spectral resolution of 1200 and 2400 respectively at H alpha wavelengths in the low resolution mode covering 320 to 740 nm bandpass. At a higher grating angle the SMI will deliver spectral resolution up to 5000 and 10000 with 400 and 200 micron core fibre respectively. A future red-arm will extend the simultaneous wavelength coverage up to 900 nm at a median resolution of 3000/6000 for the same flavors of IFUs. SMIs are inserted in the same fashion as the existing longslit cassettes at the SALT focal plane. Prismatic fold mirrors direct the focal plane into the fibre IFU and then back into the RSS collimator after the fibres are routed 180 deg within the cassette and formatted into a pseudo-slit. Fold prisms ensure that the spectrograph collimator continues to see the same focal plane. In this paper we describe the design, fabrication, assembly and characterization of Slit Mask IFUs.

preprint2022arXiv

The MaNGA FIREFLY Value-Added-Catalogue: resolved stellar populations of 10,010 nearby galaxies

We present the MaNGA FIREFLY Value-Added-Catalogue (VAC) - a catalogue of ~3.7 million spatially resolved stellar population properties across 10,010 nearby galaxies from the final data release of the MaNGA survey. The full spectral fitting code firefly is employed to derive parameters such as stellar ages, metallicities, stellar and remnant masses, star formation histories, star formation rates and dust attenuation. In addition to Voronoi-binned measurements, our VAC also provides global properties, such as central values and radial gradients. Two variants of the VAC are available: presenting the results from fits using the M11-MILES and the novel MaStar stellar population models. MaStar allows to constrain the fit over the whole MaNGA wavelength range, extends the age-metallicity parameter space, and uses empirical spectra from the same instrument as MaNGA. The fits employing MaStar models find on average slightly younger ages, higher mass-weighted metallicities and smaller colour excesses. These differences are reduced when matching wavelength range and converging template grids. We further report that FIREFLY stellar masses are systematically lower by ~0.3 dex than masses from the MaNGA PCA and Pipe3D VACs, but match masses from the NSA best with only ~0.1 dex difference. Finally, we show that FIREFLY stellar ages correlate with spectral index age indicators H$δ_A$ and $D_n$(4000), though with a clear additional metallicity dependence.

preprint2020arXiv

Fiber positioning in microlens-fiber coupled integral field unit

A generic fiber positioning strategy and a fabrication path are presented for microlens-fiber-coupled integral field units. It is assumed that microlens-produced micro-images are carried to the spectrograph input through step-index,multi-mode fiber, but our results apply to micro-pupil reimaging applications as well. Considered are the performance trades between the filling percentage of the fiber core with the micro-image versus throughput and observing efficiency.A merit function is defined as the product of the transmission efficiency and the etendue loss. For a hexagonal packing of spatial elements, the merit function has been found to be maximized to 94% of an ideal fiber IFU merit value (which has zero transmission loss and does not increase the etendue) with a microlens-fiber alignment (centering) tolerance of 1 um RMS. The maximum acceptable relative tilt between the fiber and the microlens face has been analyzed through optical modeling and found to be 0.3 degree RMS for input f-ratio slower than f/3.5 but it is much more relaxed for faster beams. Several options of fabricating fiber holders have been compared to identify cost-effective solutions that deliver the desired fiber positioning accuracy. Femto-second laser-drilling methods deliver holes arrayed on plates with a position and diameter accuracy of 1.5 um RMS, and with an aspect ratio of 1:10. A commercial vendor produces plates with thickness of 5 mm, but with similar (1 um RMS) positioning accuracy. Both of these techniques are found to be moderately expensive. A purely photo-lithographic technique performed at WCAM (a facility at the University of Wisconsin, Madison), in tandem with deep reactive ion etching, has been used to produce a repeatable recipe with 100% yield. Photo-lithography is more precise (0.5 um RMS) in terms of hole positioning and similar diameter accuracy (1 um RMS).

preprint2020arXiv

SDSS-IV MaNGA: Modeling the Spectral Line Spread Function to Sub-Percent Accuracy

The SDSS-IV Mapping Nearby Galaxies at APO (MaNGA) program has been operating from 2014-2020, and has now observed a sample of 9,269 galaxies in the low redshift universe (z ~ 0.05) with integral-field spectroscopy. With rest-optical (λλ0.36 - 1.0 um) spectral resolution R ~ 2000 the instrumental spectral line-spread function (LSF) typically has 1sigma width of about 70 km/s, which poses a challenge for the study of the typically 20-30 km/s velocity dispersion of the ionized gas in present-day disk galaxies. In this contribution, we present a major revision of the MaNGA data pipeline architecture, focusing particularly on a variety of factors impacting the effective LSF (e.g., undersampling, spectral rectification, and data cube construction). Through comparison with external assessments of the MaNGA data provided by substantially higher-resolution R ~ 10,000 instruments we demonstrate that the revised MPL-10 pipeline measures the instrumental line spread function sufficiently accurately (<= 0.6% systematic, 2% random around the wavelength of Halpha) that it enables reliable measurements of astrophysical velocity dispersions sigma_Halpha ~ 20 km/s for spaxels with emission lines detected at SNR > 50. Velocity dispersions derived from [O II], Hbeta, [O III], [N II], and [S II] are consistent with those derived from Halpha to within about 2% at sigma_Halpha > 30 km/s. Although the impact of these changes to the estimated LSF will be minimal at velocity dispersions greater than about 100 km/s, scientific results from previous data releases that are based on dispersions far below the instrumental resolution should be reevaulated.

preprint2019arXiv

SDSS-IV MaNGA: Variations in the N/O -- O/H relation bias metallicity gradient measurements

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.

preprint2016arXiv

Does the Milky Way Obey Spiral Galaxy Scaling Relations?

It is crucial to understand how the Milky Way, the galaxy we can study in the most intimate detail, fits in amongst other galaxies. Key examples include the Tully-Fisher relation (TFR) --- i.e., the tight correlation between luminosity ($L$) and rotational velocity ($V_\textrm{rot}$) --- and the 3-dimensional luminosity-velocity-radius ($LVR$) scaling relation. Several past studies have characterized the MW as a 1--1.5$σ$ outlier to the TFR. This study reexamines such comparisons using new estimates of MW properties that are robust to many of the systematic uncertainties that have been a problem in the past and are based on assumptions consistent with those used for other spiral galaxies. Comparing to scaling relations derived from modern extragalactic data, we find that our Galaxy's properties are in excellent agreement with TFRs defined using any SDSS-filter absolute magnitude, stellar mass, or baryonic mass as the $L$ proxy. We next utilize disk scale length ($R_\textrm{d}$) measurements to extend this investigation to the $LVR$ relation. Here we find that our Galaxy lies farther from the relation than $\sim$90\% of other spiral galaxies, yielding $\sim$9.5$σ$ evidence that it is unusually compact for its $L$ and $V_\textrm{rot}$ (based on MW errors alone), a result that holds for all of the $L$ proxies considered. The expected $R_\textrm{d}$ for the MW from the $LVR$ relation is $\sim$5 kpc, nearly twice as large as the observed value, with error estimates placing the two in tension at the $\sim$1.4$σ$ level. The compact scale length of the Galactic disk could be related to other ways in which the MW has been found to be anomalous.

preprint2016arXiv

SDSS IV MaNGA - Spatially resolved diagnostic diagrams: A proof that many galaxies are LIERs

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$α$ equivalent width of line emission ($<$ 3 Å) 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$α$ surface brightness radial profiles are always shallower than $\rm 1/r^{2}$ and the line ratio [OIII]$λ$5007/[OII]$λ$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.

preprint2016arXiv

SDSS-IV MaNGA IFS Galaxy Survey --- Survey Design, Execution, and Initial Data Quality

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.

preprint2016arXiv

SDSS-IV MaNGA: Bulge-Disc Decomposition of IFU Datacubes (BUDDI)

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.

preprint2016arXiv

SDSS-IV MaNGA: Faint quenched galaxies I- Sample selection and evidence for environmental quenching

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α} <2$ Å, and all have red colours $(u-r)>1.9$. They lie on the size-magnitude and $σ$-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 $ρ_{proj} = 8.2\pm2.0$ Mpc$^{-2}$, compared to $ρ_{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.

preprint2016arXiv

SDSS-IV MaNGA: stellar population gradients as a function of galaxy environment

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.

preprint2016arXiv

The Data Reduction Pipeline for the SDSS-IV MaNGA IFU Galaxy Survey

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.

preprint2016arXiv

The SDSS-IV extended Baryon Oscillation Spectroscopic Survey: Overview and Early Data

The Extended Baryon Oscillation Spectroscopic Survey (eBOSS) will conduct novel cosmological observations using the BOSS spectrograph at Apache Point Observatory. Observations will be simultaneous with the Time Domain Spectroscopic Survey (TDSS) designed for variability studies and the Spectroscopic Identification of eROSITA Sources (SPIDERS) program designed for studies of X-ray sources. eBOSS will use four different tracers to measure the distance-redshift relation with baryon acoustic oscillations (BAO). Using more than 250,000 new, spectroscopically confirmed luminous red galaxies at a median redshift z=0.72, we project that eBOSS will yield measurements of $d_A(z)$ to an accuracy of 1.2% and measurements of H(z) to 2.1% when combined with the z>0.6 sample of BOSS galaxies. With ~195,000 new emission line galaxy redshifts, we expect BAO measurements of $d_A(z)$ to an accuracy of 3.1% and H(z) to 4.7% at an effective redshift of z= 0.87. A sample of more than 500,000 spectroscopically-confirmed quasars will provide the first BAO distance measurements over the redshift range 0.9<z<2.2, with expected precision of 2.8% and 4.2% on $d_A(z)$ and H(z), respectively. Finally, with 60,000 new quasars and re-observation of 60,000 quasars known from BOSS, we will obtain new Lyman-alpha forest measurements at redshifts z>2.1; these new data will enhance the precision of $d_A(z)$ and H(z) by a factor of 1.44 relative to BOSS. Furthermore, eBOSS will provide improved tests of General Relativity on cosmological scales through redshift-space distortion measurements, improved tests for non-Gaussianity in the primordial density field, and new constraints on the summed mass of all neutrino species. Here, we provide an overview of the cosmological goals, spectroscopic target sample, demonstration of spectral quality from early data, and projected cosmological constraints from eBOSS.

preprint2016arXiv

The Time-Domain Spectroscopic Survey: Understanding the Optically Variable Sky with SEQUELS in SDSS-III

The Time-Domain Spectroscopic Survey (TDSS) is an SDSS-IV eBOSS subproject primarily aimed at obtaining identification spectra of ~220,000 optically-variable objects systematically selected from SDSS/Pan-STARRS1 multi-epoch imaging. We present a preview of the science enabled by TDSS, based on TDSS spectra taken over ~320 deg^2 of sky as part of the SEQUELS survey in SDSS-III, which is in part a pilot survey for eBOSS in SDSS-IV. Using the 15,746 TDSS-selected single-epoch spectra of photometrically variable objects in SEQUELS, we determine the demographics of our variability-selected sample, and investigate the unique spectral characteristics inherent in samples selected by variability. We show that variability-based selection of quasars complements color-based selection by selecting additional redder quasars, and mitigates redshift biases to produce a smooth quasar redshift distribution over a wide range of redshifts. The resulting quasar sample contains systematically higher fractions of blazars and broad absorption line quasars than from color-selected samples. Similarly, we show that M-dwarfs in the TDSS-selected stellar sample have systematically higher chromospheric active fractions than the underlying M-dwarf population, based on their H-alpha emission. TDSS also contains a large number of RR Lyrae and eclipsing binary stars with main-sequence colors, including a few composite-spectrum binaries. Finally, our visual inspection of TDSS spectra uncovers a significant number of peculiar spectra, and we highlight a few cases of these interesting objects. With a factor of ~15 more spectra, the main TDSS survey in SDSS-IV will leverage the lessons learned from these early results for a variety of time-domain science applications.

preprint2015arXiv

Observing Strategy for the SDSS-IV/MaNGA IFU Galaxy Survey

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

preprint2015arXiv

P-MaNGA: Gradients in Recent Star Formation Histories as Diagnostics for Galaxy Growth and Death

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.

preprint2015arXiv

The Time Domain Spectroscopic Survey: Variable Object Selection and Anticipated Results

We present the selection algorithm and anticipated results for the Time Domain Spectroscopic Survey (TDSS). TDSS is an SDSS-IV eBOSS subproject that will provide initial identification spectra of approximately 220,000 luminosity-variable objects (variable stars and AGN) across 7,500 square degrees selected from a combination of SDSS and multi-epoch Pan-STARRS1 photometry. TDSS will be the largest spectroscopic survey to explicitly target variable objects, avoiding pre-selection on the basis of colors or detailed modeling of specific variability characteristics. Kernel Density Estimate (KDE) analysis of our target population performed on SDSS Stripe 82 data suggests our target sample will be 95% pure (meaning 95% of objects we select have genuine luminosity variability of a few magnitudes or more). Our final spectroscopic sample will contain roughly 135,000 quasars and 85,000 stellar variables, approximately 4,000 of which will be RR Lyrae stars which may be used as outer Milky Way probes. The variability-selected quasar population has a smoother redshift distribution than a color-selected sample, and variability measurements similar to those we develop here may be used to make more uniform quasar samples in large surveys. The stellar variable targets are distributed fairly uniformly across color space, indicating that TDSS will obtain spectra for a wide variety of stellar variables including pulsating variables, stars with significant chromospheric activity, cataclysmic variables and eclipsing binaries. TDSS will serve as a pathfinder mission to identify and characterize the multitude of variable objects that will be detected photometrically in even larger variability surveys such as LSST.

preprint2014arXiv

Near-Infrared Structure of Fast and Slow Rotating Disk Galaxies

We investigate the stellar disk structure of six nearby edge-on spiral galaxies using high-resolution JHKs-band images and 3D radiative transfer models. To explore how mass and environment shape spiral disks, we selected galaxies with rotational velocities between 69 < Vrot < 245 km/sec, and two with unusual morphologies. We find a wide diversity of disk structure. Of the fast-rotating (Vrot > 150 km/sec) galaxies, only NGC 4013 has the super-thin+thin+thick nested disk structure seen in NGC 891 and the Milky Way, albeit with decreased oblateness, while NGC 1055, a disturbed massive spiral galaxy, contains disks with hz $\lesssim$ 200 pc. NGC 4565, another fast-rotator, contains a prominent ring at a radius ~5 kpc but no super-thin disk. Despite these differences, all fast-rotating galaxies in our sample have inner truncations in at least one of their disks. These truncations lead to Freeman Type II profiles when projected face-on. Slow-rotating galaxies are less complex, lacking inner disk truncations and requiring fewer disk components to reproduce their light distributions. Super-thin disk components in undisturbed disks contribute ~25% of the total Ks-band light, up to that of the thin-disk contribution. The presence of super-thin disks correlates with infrared flux ratios; galaxies with super-thin disks have f(Ks)/f(60 $μ$m) $\leq$ 0.12 for integrated light, consistent with super-thin disks being regions of on-going star-formation. Attenuation-corrected vertical color gradients in (J-Ks) correlate with the observed disk structure and are consistent with population gradients with young-to-intermediate ages closer to the midplane, indicating that disk heating-or cooling-is a ubiquitous phenomenon.

preprint2014arXiv

Overview of the SDSS-IV MaNGA Survey: Mapping Nearby Galaxies at Apache Point Observatory

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.

preprint2014arXiv

Spatial and Kinematic Distributions of Transition Populations in Intermediate Redshift Galaxy Clusters

We analyze the spatial and velocity distributions of confirmed members in five massive clusters of galaxies at intermediate redshift (0.5 < z < 0.9) to investigate the physical processes driving galaxy evolution. Based on spectral classifications derived from broad- and narrow-band photometry, we define four distinct galaxy populations representing different evolutionary stages: red sequence (RS) galaxies, blue cloud (BC) galaxies, green valley (GV) galaxies, and luminous compact blue galaxies (LCBGs). For each galaxy class, we derive the projected spatial and velocity distribution and characterize the degree of subclustering. We find that RS, BC, and GV galaxies in these clusters have similar velocity distributions, but that BC and GV galaxies tend to avoid the core of the two z~0.55 clusters. GV galaxies exhibit subclustering properties similar to RS galaxies, but their radial velocity distribution is significantly platykurtic compared to the RS galaxies. The absence of GV galaxies in the cluster cores may explain their somewhat prolonged star-formation history. The LCBGs appear to have recently fallen into the cluster based on their larger velocity dispersion, absence from the cores of the clusters, and different radial velocity distribution than the RS galaxies. Both LCBG and BC galaxies show a high degree of subclustering on the smallest scales, leading us to conclude that star formation is likely triggered by galaxy-galaxy interactions during infall into the cluster.

preprint2014arXiv

The DiskMass Survey. VIII. On the Relationship Between Disk Stability and Star Formation

We study the relationship between the stability level of late-type galaxy disks and their star-formation activity using integral-field gaseous and stellar kinematic data. Specifically, we compare the two-component (gas+stars) stability parameter from Romeo & Wiegert (Q_RW), incorporating stellar kinematic data for the first time, and the star-formation rate estimated from 21cm continuum emission. We determine the stability level of each disk probabilistically using a Bayesian analysis of our data and a simple dynamical model. Our method incorporates the shape of the stellar velocity ellipsoid (SVE) and yields robust SVE measurements for over 90% of our sample. Averaging over this subsample, we find a meridional shape of sigma_z/sigma_R = 0.51^{+0.36}_{-0.25} for the SVE and, at 1.5 disk scale lengths, a stability parameter of Q_RW = 2.0 +/- 0.9. We also find that the disk-averaged star-formation-rate surface density (Sigma-dot_e,*) is correlated with the disk-averaged gas and stellar mass surface densities (Sigma_e,g and Sigma_e,*) and anti-correlated with Q_RW. We show that an anti-correlation between Sigma-dot_e,* and Q_RW can be predicted using empirical scaling relations, such that this outcome is consistent with well-established statistical properties of star-forming galaxies. Interestingly, Sigma-dot_e,* is not correlated with the gas-only or star-only Toomre parameters, demonstrating the merit of calculating a multi-component stability parameter when comparing to star-formation activity. Finally, our results are consistent with the Ostriker et al. model of self-regulated star-formation, which predicts Sigma-dot_e,*/Sigma_e,g/sqrt(Sigma_e,*). Based on this and other theoretical expectations, we discuss the possibility of a physical link between disk stability level and star-formation rate in light of our empirical results.

preprint2014arXiv

The Link Between Light and Mass in Late-type Spiral Galaxy Disks

We present the correlation between the extrapolated central disk surface brightness (mu) and extrapolated central surface mass density (Sigma) for galaxies in the DiskMass sample. This mu-Sigma-relation has a small scatter of 30% at the high-surface-brightness (HSB) end. At the low surface brightness (LSB) end, galaxies fall above the mu-Sigma-relation, which we attribute to their higher dark matter content. After correcting for the dark matter, as well as for the contribution of gas and the effects of radial gradients in the disk, the LSB end falls back on the linear mu-Sigma-relation. The resulting scatter about the corrected mu-Sigma-relation is 25% at the HSB end, and about 50% at the LSB end. The intrinsic scatter in the mu-Sigma-relation is estimated to be 10% to 20%. Thus, if the surface brightness is known, the stellar surface mass density is known to within 10-20% (random error). Assuming disks have an exponential vertical distribution of mass, the average (M_L)_K is 0.24 Msun/Lsun, with an intrinsic scatter around the mean of at most 0.05 Msun/Lsun. This value for (M/L)_K is 20% smaller than we found in Martinsson et al., mainly due to the correction for dark matter applied here. This small scatter means that among the galaxies in our sample variations in scale height, vertical density profile shape, and/or the ratio of vertical over radial velocity dispersion must be small.

preprint2013arXiv

A Pilot for a VLA HI Deep Field

High-resolution 21-cm HI deep fields provide spatially and kinematically resolved neutral gas maps at different redshifts, which are key to understanding galaxy evolution across cosmic time and testing predictions of cosmological simulations. Here we present results from a pilot for the COSMOS HI Large Extragalactic Survey (CHILES) done with the Karl G. Jansky Very Large Array (VLA). We take advantage of the newly expanded capabilities of the telescope to probe the redshift interval 0<z<0.193 in one observation. We observe the COSMOS field for 50 hours, which contains 413 galaxies with optical spectroscopic redshifts in the imaged field of view of 34' x 34' and the observed redshift interval. We have detected neutral hydrogen gas in 33 galaxies in different environments spanning the probed redshift range, including three without a previously known spectroscopic redshift. The detections have a range of HI and stellar masses, indicating the diversity of galaxies we are probing. We discuss the observations, data reduction, results and highlight interesting detections. We find that the VLA's B-array is the ideal configuration for HI deep fields since its long spacings mitigate RFI. This pilot shows that the VLA is ready to carry out such a survey, and serves as a test for future HI deep fields planned with other SKA pathfinders.

preprint2013arXiv

Discovery of Super-Thin Disks in Nearby Edge-on Spiral Galaxies

We report the identification of a super-thin disk (h_z~60 pc) in the edge-on spiral galaxy NGC 891. This component is only apparent after we perform a physically motivated attenuation correction, based on detailed radiation transfer models, to our sub-arcsecond resolution near-infrared imaging. In addition to the super-thin disk, we are also find several structural features near the center of NGC 891, including an inner disk truncation at ~3 kpc. Inner disk truncations may be commonplace among massive spiral galaxies, possibly due to the effects of instabilities, such as bars. Having successfully demonstrated our methods, we are poised to apply them to a small sample of nearby edge-on galaxies, consisting both of massive and low-mass spirals.

preprint2013arXiv

Near-Infrared Detection of a Super-Thin Disk in NGC 891

We probe the disk structure of the nearby, massive, edge-on spiral galaxy NGC 891 with sub-arcsecond resolution JHKs-band images covering ~+/-10 kpc in radius and +/-5 kpc in height. We measure intrinsic surface brightness profiles using realistic attenuation corrections constrained from near- and mid-infrared (Spitzer) color maps and three-dimensional Monte-Carlo radiative-transfer models. In addition to the well-known thin and thick disks, a super-thin disk with 60-80 pc scale-height - comparable to the star-forming disk of the Milky Way - is visibly evident and required to fit the attenuation-corrected light distribution. Asymmetries in the super-thin disk light profile are indicative of young, hot stars producing regions of excess luminosity and bluer (attenuation-corrected) near-infrared color. To fit the inner regions of NGC 891, these disks must be truncated within ~3 kpc, with almost all their luminosity redistributed in a bar-like structure 50% thicker than the thin disk. There appears to be no classical bulge but rather a nuclear continuation of the super-thin disk. The super-thin, thin, thick, and bar components contribute roughly 30%, 42%, 13%, and 15% (respectively) to the total Ks-band luminosity. Disk axial ratios (length/height) decrease from 30 to 3 from super-thin to thick components. Both exponential and sech^2 vertical SB profiles fit the data equally well. We find that the super-thin disk is significantly brighter in the Ks-band than typically assumed in integrated SED models of NGC 891: it appears that in these models the excess flux, likely produced by young stars in the super-thin disk, has been mistakenly attributed to the thin disk.

preprint2013arXiv

The DiskMass Survey. VI. Gas and stellar kinematics in spiral galaxies from PPak integral-field spectroscopy

We present ionized-gas (OIII) and stellar kinematics (velocities and velocity dispersions) for 30 nearly face-on spiral galaxies out to as much as three disk scale lengths (h_R). These data have been derived from PPak IFU spectroscopy (4980-5370A), observed at a mean resolution of R=7700 (sigma_inst=17km/s). These data are a fundamental product of our survey and will be used in companion papers to, e.g., derive the detailed (baryonic+dark) mass budget of each galaxy in our sample. Our presentation provides a comprehensive description of the observing strategy, data reduction, and analysis. Along with a clear presentation of the data, we demonstrate: (1) The OIII and stellar rotation curves exhibit a clear signature of asymmetric drift with a rotation difference that is 11% of the maximum rotation speed of the galaxy disk, comparable to measurements in the solar neighborhood in the Milky Way. (2) The e-folding length of the stellar velocity dispersion is two times h_R on average, as expected for a disk with a constant scale height and mass-to-light ratio, with a scatter that is notably smaller for massive, high-surface-brightness disks in the most luminous galaxies. (3) At radii larger than 1.5 h_R, the stellar velocity dispersion tends to decline slower than the best-fitting exponential function, which may be due to an increase in the disk mass-to-light ratio, disk flaring, or disk heating by the dark-matter halo. (4) A strong correlation exists between the central vertical stellar velocity dispersion of the disks and their circular rotational speed at 2.2 h_R, with a zero point indicating that galaxy disks are submaximal. Moreover, weak but consistent correlations exist such that disks with a fainter central surface brightness in bluer and less luminous galaxies of later morphological types are kinematically colder with respect to their rotational velocities.

preprint2013arXiv

The DiskMass Survey. VII. The distribution of luminous and dark matter in spiral galaxies

We present dynamically-determined rotation-curve mass decompositions of 30 spiral galaxies, which were carried out to test the maximum-disk hypothesis and to quantify properties of their dark-matter (DM) halos. We used measured vertical velocity dispersions of the disk stars to calculate dynamical mass surface densities. Together with our atomic and molecular gas mass surface densities, we derived the stellar mass surface densities, and thus have absolute measurements of all dominant baryonic components. Using K-band surface brightness profiles, we calculated the K-band mass-to-light ratio of the stellar disks (M/L). Our result is consistent with all galaxies in the sample having equal M/L, with a sample average and scatter of <M/L>=0.31+/-0.07. Rotation-curves of the baryonic components were calculated from their mass surface densities, and used with circular-speed measurements to derive the structural parameters of the DM halos, modeled as either a pseudo-isothermal sphere (pISO) or an NFW halo. All galaxies in our sample are submaximal, such that at 2.2 disk scale lengths (hR) the ratios between the baryonic and total rotation-curves (Fb^{2.2hR}) are less than 0.75. We find this ratio to be nearly constant between 1-6 hR within individual galaxies. We find a sample average and scatter of <Fb^{2.2hR}>=0.57+/-0.07, with trends of larger Fb^{2.2hR} for more luminous and higher-surface-brightness galaxies. To enforce these being maximal, we need to scale M/L by a factor 3.6 on average. The DM rotation curves are marginally better fit by a pISO than by an NFW halo. For the nominal-M/L (submaximal) case, the derived NFW-halo parameters have values consistent with LCDM N-body simulations, suggesting that the baryonic matter has only had a minor effect on the DM distribution. In contrast, maximum-M/L decompositions yield halo concentrations that are too low compared to the LCDM simulations.

preprint2013arXiv

The Stability of Galaxy Disks

We calculate the stellar surface mass density (Sigma_*) and two-component (gas+stars) disk stability (Q_RW) for 25 late-type galaxies from the DiskMass Survey. These calculations are based on fits of a dynamical model to our ionized-gas and stellar kinematic data performed using a Markov Chain Monte Carlo sampling of the Bayesian posterior. Marginalizing over all galaxies, we find a median value of Q_RW=2.0+/-0.9 at 1.5 scale lengths. We also find that Q_RW is anti-correlated with the star-formation rate surface density (Sigma_SFR), which can be predicted using a closed set of empirical scaling relations. Finally, we find that the star-formation efficiency (Sigma_SFR/Sigma_g) is correlated with Sigma_* and weakly anti-correlated with Q_RW. The former is consistent with an equilibrium prediction of Sigma_SFR/Sigma_g propto Sigma_*^{1/2}. Despite its order-of-magnitude range, we find no correlation of Sigma_SFR/Sigma_g/Sigma_*^{1/2} with any other physical quantity derived by our study.

preprint2012arXiv

HexPak and GradPak: variable-pitch dual-head IFUs for the WIYN 3.5m Telescope Bench Spectrograph

We describe the design, construction, and expected performance of two new fiber integral field units (IFUs) --- HexPak and GradPak --- for the WIYN 3.5m Telescope Nasmyth focus and Bench Spectrograph. These are the first IFUs to provide formatted fiber integral field spectroscopy with simultaneous sampling of varying angular scales. HexPak and GradPak are in a single cable with a dual-head design, permitting easy switching between the two different IFU heads on the telescope without changing the spectrograph feed: the two heads feed a variable-width double-slit. Each IFU head is comprised of a fixed arrangement of fibers with a range of fiber diameters. The layout and diameters of the fibers within each array are scientifically-driven for observations of galaxies: HexPak is designed to observe face-on spiral or spheroidal galaxies while GradPak is optimized for edge-on studies of galaxy disks. HexPak is a hexagonal array of 2.9 arcsec fibers subtending a 40.9 arcsec diameter, with a high-resolution circular core of 0.94 arcsec fibers subtending 6 arcsec diameter. GradPak is a 39 by 55 arcsec rectangular array with rows of fibers of increasing diameter from angular scales of 1.9 arcsec to 5.6 arcsec across the array. The variable pitch of these IFU heads allows for adequate sampling of light profile gradients while maintaining the photon limit at different scales.

preprint2012arXiv

Impact of surface-polish on the angular and wavelength dependence of fiber focal ratio degradation

We present measurements of how multimode fiber focal-ratio degradation (FRD) and throughput vary with levels of fiber surface polish from 60 to 0.5 micron grit. Measurements used full-beam and laser injection methods at wavelengths between 0.4 and 0.8 microns on 17 meter lengths of Polymicro FBP 300 and 400 micron core fiber. Full-beam injection probed input focal-ratios between f/3 and f/13.5, while laser injection allowed us to isolate FRD at discrete injection angles up to 17 degrees (f/1.6 marginal ray). We find (1) FRD effects decrease as grit size decreases, with the largest gains in beam quality occurring at grit sizes above 5 microns; (2) total throughput increases as grit size decreases, reaching 90% at 790 nm with the finest polishing levels; (3) total throughput is higher at redder wavelengths for coarser polishing grit, indicating surface-scattering as the primary source of loss. We also quantify the angular dependence of FRD as a function of polishing level. Our results indicate that a commonly adopted micro-bending model for FRD is a poor descriptor of the observed phenomenon.

preprint2012arXiv

The 3-Dimensional Structure of NGC 891 and M51

We investigate the three-dimensional structure of the nearby edge-on spiral galaxy NGC 891 using 3D Monte Carlo radiative transfer models, with realistic spiral structure and fractally clumped dust. Using the spiral and clumpiness parameters found from recently completed scattered light models we produce lower resolution SED models which reproduce the global UV-to-FIR SED of NGC 891. Our models contain a color gradient across the major axis of the galaxy - similar to what is seen in images of the NGC 891. With minor adjustment our SED models are able to match the majority of M51's SED, a similar galaxy at a near face-on different inclination.

preprint2011arXiv

Galaxy Disks are Submaximal

We measure the contribution of galaxy disks to the overall gravitational potential of 30 nearly face-on intermediate-to-late-type spirals from the DiskMass Survey. The central vertical velocity dispersion of the disk stars, sigma(z,R=0), is related to the maximum rotation speed (Vmax) as sigma(z,R=0) ~ 0.26 Vmax, consistent with previous measurements for edge-on disk galaxies and a mean stellar velocity ellipsoid axial ratio sigma(z) / sigma(R) = 0.6. For reasonable values of disk oblateness, this relation implies these galaxy disks are submaximal. We find disks in our sample contribute only 15% to 30% of the dynamical mass within 2.2 disk scale-lengths (hR), with percentages increasing systematically with luminosity, rotation speed and redder color. These trends indicate the mass ratio of disk-to-total matter remains at or below 50% at 2.2 hR even for the most extreme, fast-rotating disks (Vmax > 300 km/s), of the reddest rest-frame, face-on color (B-K ~ 4 mag), and highest luminosity (M(K)<-26.5 mag). Therefore, spiral disks in general should be submaximal. Our results imply that the stellar mass-to-light ratio and hence the accounting of baryons in stars should be lowered by at least a factor of 3.

preprint2011arXiv

The 3-Dimensional Distribution of Dust in NGC 891

We produce three-dimensional Monte-Carlo radiative transfer models of the edge-on spiral galaxy NGC 891, a fast-rotating galaxy thought to be an analogue to the Milky Way. The models contain realistic spiral arms and a fractal distribution of clumpy dust. We fit our models to Hubble Space Telescope images corresponding to the B and I bands, using shapelet analysis and a genetic algorithm to generate 30 statistically best-fitting models. These models have a strong preference for spirality and clumpiness, with average face-on attenuation decreasing from 0.24(0.16) to 0.03(0.03) mag in the B(I) band between 0.5 and 2 radial scale-lengths. Most of the attenuation comes from small high-density clumps with low (<10%) filling factors. The fraction of dust in clumps is broadly consistent with results from fitting NGC 891's spectral energy distribution. Because of scattering effects and the intermixed nature of the dust and starlight, attenuation is smaller and less wavelength-dependent than the integrated dust column-density. Our clumpy models typically have higher attenuation at low inclinations than previous radiative transfer models using smooth distributions of stars and dust, but similar attenuation at inclinations above 70 degrees. At all inclinations most clumpy models have less attenuation than expected from previous estimates based on minimizing scatter in the Tully-Fisher relation. Mass-to-light ratios are higher and the intrinsic scatter in the Tully-Fisher relation is larger than previously expected for galaxies similar to NGC 891. The attenuation curve changes as a function of inclination, with R_(B,B-I)=A_(B)/E(B-I) increasing by ~0.75 from face-on to near-edge-on orientations.

preprint2011arXiv

The DiskMass Survey. IV. The Dark-Matter-Dominated Galaxy UGC 463

We present a detailed and unique mass budget for the high-surface-brightness galaxy UGC 463, showing it is dominated by dark matter (DM) at radii beyond one scale length (h_R) and has a baryonic-to-DM mass ratio of approximately 1:3 within 4.2 h_R. Assuming a constant scale height (h_z, calculated via an empirical oblateness relation), we calculate dynamical disk mass surface densities from stellar kinematics, which provide vertical velocity dispersions after correcting for the shape of the stellar velocity ellipsoid (measured to have sigma_theta/sigma_R=1.04 +/- 0.22 and sigma_z/sigma_R=0.48 +/- 0.09). We isolate the stellar mass surface density by accounting for all gas mass components and find an average K-band mass-to-light ratio of 0.22 +/- 0.09 (ran) ^{+0.16}_{-0.15} (sys) M_{sun}/L_{sun}^{K}; Zibetti et al. and Bell et al. predict, respectively, 0.56 and 3.6 times our dynamical value based on stellar-population-synthesis modeling. The baryonic matter is submaximal by a factor of ~3 in mass and the baryonic-to-total circular-speed ratio is 0.61^{+0.07}_{-0.09} (ran) ^{+0.12}_{-0.18} (sys) at 2.2 h_R; however, the disk is globally stable with a multi-component stability that decreases asymptotically with radius to Q~2. We directly calculate the circular speed of the DM halo by subtracting the baryonic contribution to the total circular speed; the result is equally well described by either a Navarro-Frenk-White halo or a pseudo-isothermal sphere. The volume density is dominated by DM at heights of |z|>1.6 h_z for radii of R > h_R. As is shown in follow-up papers, UGC 463 is just one example among nearly all galaxies we have observed that contradict the hypothesis that high-surface-brightness spiral galaxies have maximal disks.

preprint2010arXiv

The DiskMass Survey. I. Overview

We present a survey of the mass surface-density of spiral disks, motivated by outstanding uncertainties in rotation-curve decompositions. Our method exploits integral-field spectroscopy to measure stellar and gas kinematics in nearly face-on galaxies sampled at 515, 660, and 860 nm, using the custom-built SparsePak and PPak instruments. A two-tiered sample, selected from the UGC, includes 146 nearly face-on galaxies, with B<14.7 and disk scale-lengths between 10 and 20 arcsec, for which we have obtained H-alpha velocity-fields; and a representative 46-galaxy subset for which we have obtained stellar velocities and velocity dispersions. Based on re-calibration of extant photometric and spectroscopic data, we show these galaxies span factors of 100 in L(K) (0.03 < L/L(K)* < 3), 8 in L(B)/L(K), 10 in R-band disk central surface-brightness, with distances between 15 and 200 Mpc. The survey is augmented by 4-70 micron Spitzer IRAC and MIPS photometry, ground-based UBVRIJHK photometry, and HI aperture-synthesis imaging. We outline the spectroscopic analysis protocol for deriving precise and accurate line-of-sight stellar velocity dispersions. Our key measurement is the dynamical disk-mass surface-density. Star-formation rates and kinematic and photometric regularity of galaxy disks are also central products of the study. The survey is designed to yield random and systematic errors small enough (i) to confirm or disprove the maximum-disk hypothesis for intermediate-type disk galaxies, (ii) to provide an absolute calibration of the stellar mass-to-light ratio well below uncertainties in present-day stellar-population synthesis models, and (iii) to make significant progress in defining the shape of dark halos in the inner regions of disk galaxies.

preprint2010arXiv

The DiskMass Survey. II. Error Budget

We present a performance analysis of the DiskMass Survey. The survey uses collisionless tracers in the form of disk stars to measure the surface-density of spiral disks, to provide an absolute calibration of the stellar mass-to-light ratio, and to yield robust estimates of the dark-matter halo density profile in the inner regions of galaxies. We find a disk inclination range of 25-35 degrees is optimal for our measurements, consistent with our survey design to select nearly face-on galaxies. Uncertainties in disk scale-heights are significant, but can be estimated from radial scale-lengths to 25% now, and more precisely in the future. We detail the spectroscopic analysis used to derive line-of-sight velocity dispersions, precise at low surface-brightness, and accurate in the presence of composite stellar populations. Our methods take full advantage of large-grasp integral-field spectroscopy and an extensive library of observed stars. We show that the baryon-to-total mass fraction (F_b) is not a well-defined observational quantity because it is coupled to the halo mass model. This remains true even when the disk mass is known and spatially-extended rotation curves are available. In contrast, the fraction of the rotation speed supplied by the disk at 2.2 scale lengths (disk maximality) is a robust observational indicator of the baryonic disk contribution to the potential. We construct the error-budget for the key quantities: dynamical disk mass surface-density, disk stellar mass-to-light ratio, and disk maximality (V_disk / V_circular). Random and systematic errors in these quantities for individual galaxies will be ~25%, while survey precision for sample quartiles are reduced to 10%, largely devoid of systematic errors outside of distance uncertainties.

preprint2007arXiv

Integral Field Unit Observations of NGC 4302: Kinematics of the Diffuse Ionized Gas Halo

We present moderate resolution spectroscopy of extraplanar diffuse ionized gas (EDIG) emission in the edge-on spiral galaxy NGC 4302. The spectra were obtained with the SparsePak integral field unit (IFU) at the WIYN Observatory. The spectra are used to construct position-velocity (PV) diagrams at several ranges of heights above the midplane. Azimuthal velocities are directly extracted from the PV diagrams using the envelope tracing method, and indicate an extremely steep dropoff in rotational velocity with increasing height, with magnitude ~30 km/s/kpc. We find evidence for a radial variation in the velocity gradient on the receding side. We have also performed artificial observations of galaxy models in an attempt to match the PV diagrams. The results of a statistical analysis also favor a gradient of ~30 km/s/kpc. We compare these results with an entirely ballistic model of disk-halo flow, and find a strong dichotomy between the observed kinematics and those predicted by the model. The disagreement is worse than we have found for other galaxies in previous studies. The conclusions of this paper are compared to results from two other galaxies, NGC 5775 and NGC 891. We find that the vertical gradient in rotation speed, per unit EDIG scale height, for all three galaxies is consistent with a constant magnitude (within the errors) of approximately 15-25 km/s/scaleheight, independent of radius. This relationship is also true within the galaxy NGC 4302. We also discuss how the gradient depends on the distribution and morphology of the EDIG and the star formation rates of the galaxies, and consequences for the origin of the gas.

preprint2006arXiv

A Search for Low Surface Brightness Structure Around Compact Narrow Emission Line Galaxies

As the most extreme members of the rapidly evolving faint blue galaxy population at intermediate redshift, the compact narrow emission line galaxies (CNELGs) are intrinsically luminous (-22 < M_B < -18) with narrow emission linewidths (30 < σ< 125 km/s). Their nature is heavily debated: they may be low-mass starbursting galaxies that will fade to present-day dwarf galaxies or bursts of star formation temporarily dominating the flux of more massive galaxies, possibly related to in situ bulge formation or the formation of cores of galaxies. We present deep, high-quality (~0.6 - 0.8 arcsec) images with CFHT of 27 CNELGs. One galaxy shows clear evidence for a tidal tail; the others are not unambiguously embedded in galactic disks. Approximately 55% of the CNELGS have sizes consistent with local dwarfs of small-to-intermediate sizes, while 45% have sizes consistent with large dwarfs or disks galaxies. At least 4 CNELGs cannot harbor substantial underlying disk material; they are low-luminosity galaxies at the present epoch (M_B > -18). Conversely, 15 are not blue enough to fade to low-luminosity dwarfs (M_B > -15.2). The majority of the CNELGs are consistent with progenitors of intermediate-luminosity dwarfs and low-luminosity spiral galaxies with small disks. CNELGs are a heterogeneous progenitor population with significant fractions (up to 44%) capable of fading into today's faint dwarfs (M_B > -15.2), while 15 to 85% may only experience an apparently extremely compact CNELG phase at intermediate redshift but remain more luminous galaxies at the present epoch.

preprint2006arXiv

Integral Field Unit Observations of NGC 891: Kinematics of the Diffuse Ionized Gas Halo

We present high and moderate spectral resolution spectroscopy of diffuse ionized gas (DIG) emission in the halo of NGC 891. The data were obtained with the SparsePak integral field unit at the WIYN Observatory. The wavelength coverage includes the [NII]6548,6583, Halpha, and [SII]6716,6731 emission lines. Position-velocity (PV) diagrams, constructed using spectra extracted from four SparsePak pointings in the halo, are used to examine the kinematics of the DIG. Using two independent methods, a vertical gradient in azimuthal velocity is found to be present in the northeast quadrant of the halo, with magnitude approximately 15-18 km/s/kpc, in agreement with results from HI observations. The kinematics of the DIG suggest that this gradient begins at approximately 1 kpc above the midplane. In another part of the halo, the southeast quadrant, the kinematics are markedly different, and suggest rotation at about 175 km/s, much slower than the disk but with no vertical gradient. We utilize an entirely ballistic model of disk-halo flow in an attempt to reproduce the kinematics observed in the northeast quadrant. Analysis shows that the velocity gradient predicted by the ballistic model is far too shallow. Based on intensity cuts made parallel to the major axis in the ballistic model and an Halpha image of NGC 891 from the literature, we conclude that the DIG halo is much more centrally concentrated than the model, suggesting that hydrodynamics dominate over ballistic motion in shaping the density structure of the halo. Velocity dispersion measurements along the minor axis of NGC 891 seem to indicate a lack of radial motions in the halo, but the uncertainties do not allow us to set firm limits.

preprint2006arXiv

The Photometric and Kinematic Structure of Face-On Disk Galaxies. I. Sample Definition, H-alpha Integral Field Spectroscopy, and HI Line-Widths

We present a survey of the photometric and kinematic properties of 39 nearby, nearly face-on disk galaxies. Our approach exploits echelle-resolution integral-field spectroscopy of the H-alpha regions, obtained with DensePak on the WIYN 3.5m telescope Bench Spectrograph. This data is complemented by HI line-profiles observed with the Nancay radio telescope for 25 of these sample galaxies. Twelve additional line-widths are available for sample galaxies from the literature. In this paper, we introduce the goals of this survey, define the sample selection algorithm, and amass the integral field spectroscopic data and HI line-widths. We establish spatially-integrated H-alpha line-widths for the sample. We test the veracity of these spatially-integrated line profiles by convolving narrow-band imaging data with velocity field information for one of the sample galaxies, PGC 38268, and also by comparing to HI line profiles. We find HI and H-alpha line profiles to be similar in width but different in shape, indicating we are observing different spatial distributions of ionized and neutral gas in largely axisymmetric systems with flat outer rotation-curves. We also find vertical velocity dispersions of the ionized disk gas within several disk scale-lengths have a median value of 18 km/s and an 80% range of 12-26 km/s. This is only a factor of ~2 larger than what is observed for neutral atomic and molecular gas. With standard assumptions for intrinsic and thermal broadening for H-alpha, this translates into a factor of three range in turbulent velocities, between 8 and 25 km/s.

preprint2005arXiv

Are There Local Analogs of Lyman Break Galaxies?

To make direct comparisons in the rest-far-ultraviolet between LBGs at z~3 and more local star-forming galaxies, we use HST/STIS to image a set of 12 nearby (z<0.05) HII galaxies in the FUV and a set of 14 luminous compact blue galaxies (LCBGs) at moderate redshift (z~0.5) in the NUV, corresponding to the rest-FUV. We then subject both sets of galaxy images and those of LBGs at z~3 to the same morphological and structural analysis. We find many qualitative and quantitative similarities between the rest-FUV characteristics of distant LBGs and of the more nearby starburst samples, including general morphologies, sizes, asymmetries, and concentrations. Along with some kinematic similarities, this implies that nearby HII galaxies and LCBGs may be reasonable local analogs of distant Lyman break galaxies.

preprint2005arXiv

Asymmetric Drift and the Stellar Velocity Ellipsoid

We present the decomposition of the stellar velocity ellipsoid using stellar velocity dispersions within a 40 deg wedge about the major-axis (sigma_maj), the epicycle approximation, and the asymmetric drift equation. Thus, we employ no fitted forms for sigma_maj and escape interpolation errors resulting from comparisons of the major and minor axes. We apply the theoretical construction of the method to integral field data taken for NGC 3949 and NGC 3982. We derive the vertical-to-radial velocity dispersion ratio (sigma_z / sigma_R) and find (1) our decomposition method is accurate and reasonable, (2) NGC 3982 appears to be rather typical of an Sb type galaxy with sigma_z / sigma_R = 0.73 (+0.13/-0.11) despite its high surface brightness and small size, and (3) NGC 3949 has a hot disk with sigma_z / sigma_R = 1.18 (+0.36/-0.28).

preprint2005arXiv

Total Galaxy Magnitudes and Effective Radii from Petrosian Magnitudes and Radii

Petrosian magnitudes were designed to help with the difficult task of determining a galaxy's total light. Although these magnitudes (taken here as the flux within 2R_P, with the inverted Petrosian index 1/eta(R_P)=0.2) can represent most of an object's flux, they do of course miss the light outside of the Petrosian aperture (2R_P). The size of this flux deficit varies monotonically with the shape of a galaxy's light-profile, i.e., its concentration. In the case of a de Vaucouleurs R^{1/4} profile, the deficit is 0.20 mag; for an R^{1/8} profile this figure rises to 0.50 mag. Here we provide a simple method for recovering total (Sersic) magnitudes from Petrosian magnitudes using only the galaxy concentration (R_90/R_50 or R_80/R_20) within the Petrosian aperture. The corrections hold to the extent that Sersic's model provides a good description of a galaxy's luminosity profile. We show how the concentration can also be used to convert Petrosian radii into effective half-light radii, enabling a robust measure of the mean effective surface brightness. Our technique is applied to the SDSS DR2 Petrosian parameters, yielding good agreement with the total magnitudes, effective radii, and mean effective surface brightnesses obtained from the NYU--VAGC Sersic R^{1/n} fits by Blanton et al. (2005). Although the corrective procedure described here is specifically applicable to the SDSS DR2 and DR3, it is generally applicable to all imaging data where any Petrosian index and concentration can be constructed.

preprint2004arXiv

Internal Kinematics of Luminous Compact Blue Galaxies

We describe the dynamical properties which may be inferred from HST/STIS spectroscopic observations of luminous compact blue galaxies (LCBGs) between 0.1<z<0.7. While the sample is homogeneous in blue rest-frame color, small size and line-width, and high surface-brightness, their detailed morphology is eclectic. Here we determine the amplitude of rotation versus random, or disturbed motions of the ionized gas. This information affirms the accuracy of dynamical mass and M/L estimates from Keck integrated line-widths, and hence also the predictions of the photometric fading of these unusual galaxies. The resolved kinematics indicates this small subset of LCBGs are dynamically hot, and unlikely to be embedded in disk systems.

preprint2000arXiv

The Asymmetry of Galaxies: Physical Morphology for Nearby and High Redshift Galaxies

We present a detailed study of rotational asymmetry in galaxies for both morphological and physical diagnostic purposes. An unambiguous method for computing asymmetry is developed, robust for both distant and nearby galaxies. By degrading real galaxy images, we test the reliability of this asymmetry measure over a range of observational conditions, e.g. spatial resolution and signal-to-noise (S/N). Compared to previous methods, this new algorithm avoids the ambiguity associated with choosing a center by using a minimization method, and successfully corrects for variations in S/N. There is, however, a strong relationship between the rotational asymmetry and physical resolution (distance at fixed spatial resolution); objects become more symmetric when less well-resolved. We further investigate asymmetry as a function of galactic radius an rotation. We find the asymmetry index has a strong radial dependence that differs vastly between Hubble types. As a result, a meaningful asymmetry index must be specified within a well-defined radius representative of the physical galaxy scale. We enumerate several viable alternatives, which excludes the use of isophotes. Asymmetry as a function of angle (A[phi]) is also a useful indicator of ellipticity and higher-order azimuthal structure. In general, we show the power of asymmetry as a morphological parameter lies in the strong correlation with (B-V) color for galaxies undergoing normal star formation, spanning all Hubble types from ellipticals to irregular galaxies. Interacting galaxies do not fall on this asymmetry-color ``fiducial sequence,'' as these galaxies are too asymmetric for their color. We propose to use this fact to distinguish between `normal' galaxies and galaxies undergoing an interaction or merger at high redshift.