OH-Suppressing Infrared Imaging Spectrograph Research Articles

Using Non-negative Matrix Factorization to Improve Calibration of the Keck OSIRIS Integral Field Spectrograph

Integral Field Spectrographs (IFS) often require non-trivial calibration techniques to process raw data. The OH Suppressing InfraRed Imaging Spectrograph at the W. M. Keck Observatory is a lenslet-based IFS that requires precise methods to associate the flux on the detector with both a wavelength and a position on the detector. During calibration scans, a single column lenslet mask is utilized to keep light from adjacent lenslet columns separate from the primary lenslet column, in order to uniquely determine spectral response of individual lenslets on the detector. Despite employing a single column lenslet mask, an issue associated with such calibration schemes may occur when light from adjacent masked lenslet columns leaks into the primary lenslet column. Incorrectly characterizing the flux due to additional light in the primary lenslet column results in one form of crosstalk between lenslet columns, which most clearly manifest as non-physical artifacts in the spectral dimension of the reduced data. We treat the problem of potentially blended calibration scans as a source separation problem and implement Non-negative Matrix Factorization (NMF) as a way to separate blended calibration scan spectra. After applying NMF to calibration scan data, extracted spectra from calibration scans show reduced crosstalk of up to 26.7% ± 0.5% while not adversely impacting the signal-to-noise ratio. Additionally, we determined the optimal number of calibration scans per lenslet column needed to create NMF factors, finding that greatest reduction crosstalk occurs when NMF factors are created using one calibration scan per lenslet column.

Observations of the global haze redistribution on Titan from 2006 to 2015 with OSIRIS at Keck

We observed Titan with the OH Suppressing InfraRed Imaging Spectrograph (OSIRIS) at the W. M. Keck observatory from 2006 through 2015 using adaptive optics. The sunlight scattered by atmospheric haze was spatially resolved in the 2.0 µm (K) band window, and the spectra were analyzed with a radiative transfer model to determine the vertical (altitude) and meridional (latitudinal) variation in the haze distribution over this time period. This study complements recent work by Karkoschka (2016) in the season of observations, in the time span and sampling interval, in wavelength coverage and spectral resolution, as well as in the radiative transfer methodology and analysis. We observe the largest meridional gradient in haze opacity above 20 km toward the northern hemisphere in January 2010. Individual observations can show significant deviations from a relatively smooth linear gradient in haze across the entire disk. The variation in haze below 20 km is rarely well-described by a simple model and there is a systematically smaller amount of haze opacity retrieved from the equator to 10° S when observing the disk with a sub-observer longitude near 150° W. This correlation with longitude suggests one of the following; a localized decrease in haze scattering, a localized increase in gas opacity, or a systematic over-estimate of the surface albedo in this region.

Reconstructing merger timelines using star cluster age distributions: the case of MCG+08-11-002

We present near infrared imaging and integral field spectroscopy of the centre of the dusty luminous infrared galaxy merger MCG+08-11-002, taken using the Near InfraRed Camera 2 (NIRC2) and the OH-Suppressing InfraRed Imaging Spectrograph (OSIRIS) on Keck II. We achieve a spatial resolution of ~25 pc in the K band, allowing us to resolve 41 star clusters in the NIRC2 images. We calculate the ages of 22/25 star clusters within the OSIRIS field using the equivalent widths of the CO 2.3$\mu$m absorption feature and the Br$\gamma$ nebular emission line. The star cluster age distribution has a clear peak at ages < 20 Myr, indicative of current starburst activity associated with the final coalescence of the progenitor galaxies. There is a possible second peak at ~65 Myr which may be a product of the previous close passage of the galaxy nuclei. We fit single and double starburst models to the star cluster age distribution and use Monte Carlo sampling combined with two-sided K-S tests to calculate the probability that the observed data are drawn from each of the best fit distributions. There is a >90 per cent chance that the data are drawn from either a single or double starburst star formation history, but stochastic sampling prevents us from distinguishing between the two scenarios. Our analysis of MCG+08-11-002 indicates that star cluster age distributions provide valuable insights into the timelines of galaxy interactions and may therefore play an important role in the future development of precise merger stage classification systems.

MAPPING THE CLUMPY STRUCTURES WITHIN SUBMILLIMETER GALAXIES USING LASER-GUIDE STAR ADAPTIVE OPTICS SPECTROSCOPY

We present the first integral-field spectroscopic observations of high-redshift submillimeter-selected galaxies (SMGs) using Laser Guide Star Adaptive Optics (LGS-AO). We target H-alpha emission of three SMGs at redshifts z~1.4-2.4 with the OH-Suppressing Infrared Imaging Spectrograph (OSIRIS) on Keck. The spatially-resolved spectroscopy of these galaxies reveals unresolved broad H-alpha line regions (FWHM>1000 km/s) likely associated with an AGN and regions of diffuse star formation traced by narrow-line H-alpha emission (FWHM<500 km/s) dominated by multiple Halpha-bright stellar clumps, each contributing 1-30% of the total clump-integrated H-alpha emission. We find that these SMGs host high star-formation rate surface densities, similar to local extreme sources, such as circumnuclear starbursts and luminous infrared galaxies. However, in contrast to these local environments, SMGs appear to be undergoing such intense activity on significantly larger spatial scales as revealed by extended H-alpha emission over 4-16 kpc. H-alpha kinematics show no evidence of ordered global motion as would be found in a disk, but rather large velocity offsets (~few x 100 km/s) between the distinct stellar clumps. Together with the asymmetric distribution of the stellar clumps around the AGN in these objects, it is unlikely that we are unveiling a clumpy disk structure as has been suggested in other high-redshift populations of star-forming galaxies. The SMG clumps in this sample may correspond to remnants of originally independent gas-rich systems that are in the process of merging, hence triggering the ultraluminous SMG phase.

A STELLAR DYNAMICAL MASS MEASUREMENT OF THE BLACK HOLE IN NGC 3998 FROM KECK ADAPTIVE OPTICS OBSERVATIONS

We present a new stellar dynamical mass measurement of the black hole in the nearby, S0 galaxy NGC 3998. By combining laser guide star adaptive optics observations obtained with the OH-Suppressing Infrared Imaging Spectrograph on the Keck II telescope with long-slit spectroscopy from the Hubble Space Telescope and the Keck I telescope, we map out the stellar kinematics on both small spatial scales, well within the black hole sphere of influence, and on large scales. We find that the galaxy is rapidly rotating and exhibits a sharp central peak in the velocity dispersion. Using the kinematics and the stellar luminosity density derived from imaging observations, we construct three-integral, orbit-based, triaxial stellar dynamical models. We find the black hole has a mass of M_BH = (8.1_{-1.9}^{+2.0}) x 10^8 M_sun, with an I-band stellar mass-to-light ratio of M/L = 5.0_{-0.4}^{+0.3} M_sun/L_sun (3-sigma uncertainties), and that the intrinsic shape of the galaxy is very round, but oblate. With the work presented here, NGC 3998 is now one of a very small number of galaxies for which both stellar and gas dynamical modeling have been used to measure the mass of the black hole. The stellar dynamical mass is nearly a factor of four larger than the previous gas dynamical black hole mass measurement. Given that this cross-check has so far only been attempted on a few galaxies with mixed results, carrying out similar studies in other objects is essential for quantifying the magnitude and distribution of the cosmic scatter in the black hole mass - host galaxy relations.

Scaling relations of star-forming regions: from kpc-sized clumps to H ii regions

We present the properties of eight star-forming regions, or ‘clumps,’ in three galaxies at z∼ 1.3 from the WiggleZ Dark Energy Survey, which are resolved with the OH Suppressing InfraRed Imaging Spectrograph (OSIRIS) integral field spectrograph. Within turbulent discs, σ∼ 90 km s−1, clumps are measured with average sizes of 1.5 kpc and average Jeans masses of 4.2 × 109 M⊙, in total accounting for 40–60 per cent of the stellar mass of the discs. These findings lend observational support to models that predict larger clumps will form as a result of higher disc velocity dispersions driven-up by cosmological gas accretion. As a consequence of the changes in global environment, it may be predicted that star-forming regions at high redshift should not resemble star-forming regions locally. Yet despite the increased sizes and dispersions, clumps and H ii regions are found to follow tight scaling relations over the range z= 0–2 for Hα size, velocity dispersion, luminosity and mass when comparing >2000 H ii regions locally and 30 clumps at z > 1 (σ∝r0.42 ± 0.03, LHα∝r2.72 ± 0.04, LHα∝σ4.18 ± 0.21 and ⁠). We discuss these results in the context of the existing simulations of clump formation and evolution, with an emphasis on the processes that drive-up the turbulent motions in the interstellar medium. Our results indicate that while the turbulence of discs may have important implications for the size and luminosity of regions which form within them, the same processes govern their formation from high redshift to the current epoch.

The WiggleZ Dark Energy Survey: high-resolution kinematics of luminous star-forming galaxies

We report evidence of ordered orbital motion in luminous star-forming galaxies at z~1.3. We present integral field spectroscopy (IFS) observations, performed with the OH Suppressing InfraRed Imaging Spectrograph (OSIRIS) system, assisted by laser guide star adaptive optics on the Keck telescope, of 13 star-forming galaxies selected from the WiggleZ Dark Energy Survey. Selected via ultraviolet and [OII] emission, the large volume of the WiggleZ survey allows the selection of sources which have comparable intrinsic luminosity and stellar mass to IFS samples at z>2. Multiple 1-2 kpc size sub-components of emission, or 'clumps', are detected within the Halpha spatial emission which extends over 6-10 kpc in 4 galaxies, resolved compact emission (r<3 kpc) is detected in 5 galaxies, and extended regions of Halpha emission are observed in the remaining 4 galaxies. We discuss these data in the context of different snapshots in a merger sequence and/or the evolutionary stages of coalescence of star-forming regions in an unstable disk. We find evidence of ordered orbital motion in galaxies as expected from disk models and the highest values of velocity dispersion (\sigma>100 km/s) in the most compact sources. This unique data set reveals that the most luminous star-forming galaxies at z>1 are gaseous unstable disks indicating that a different mode of star formation could be feeding gas to galaxies at z>1, and lending support to theories of cold dense gas flows from the intergalactic medium.

Integral Field Spectroscopy of a Candidate Disk Galaxy atz∼ 1.5 Using Laser Guide Star Adaptive Optics

We present 0.1'' resolution near-infrared integral field spectroscopy of Hα in a z = 1.4781 star-forming galaxy, Q2343-BM133. These observations were obtained with the OH Suppressing Infra-Red Imaging Spectrograph (OSIRIS) using the W. M. Keck Observatory laser guide star adaptive optics (LGS AO) system. Hα emission is resolved over a 0.8''(6.8 kpc) × 0.5''(4.3 kpc) region with a 0.1'' spatial resolution. We find a global flux of 4.2 ± 0.6 × 10-16 ergs s-1 cm-2, and detect a spatially resolved velocity gradient of ~134 km s-1 across the galaxy and a global velocity dispersion of 73 ± 9 km s-1. An upper limit of [N ]/Hα ≲ 0.12 is inferred, which implies that this galaxy is not dominated by an active galactic nucleus and has a metallicity at or below 1/2 solar metallicity. We derive a star formation rate (SFR) of 47 ± 6 M☉ yr-1 and a dereddened SFR of 66 ± 9 M☉ yr-1. Two-dimensional kinematics for Q2343-BM133 fit well with an inclined disk model, with which we estimate an enclosed mass of 4.3 × 109 M☉ within 5.5 kpc. A possible merger scenario is also presented, and cannot be fully ruled out. We derive a virial mass of 1.1 × 1010 M☉ for a disk geometry, using the observed velocity dispersion. We propose that Q2343-BM133 is currently at an early stage of disk formation at a look-back time of 9.3 Gyr.

Predictions and Strategies for Integral-Field Spectroscopy of High-Redshift Galaxies

We investigate the ability of infrared integral-field spectrographs to map the velocity fields of high redshift galaxies, presenting a formalism which may be applied to any telescope and imaging spectrograph system. We discuss the 5-sigma limiting line fluxes which current integral-field spectrographs will reach, and extend this discussion to consider future large aperture telescopes with cryogenically cooled adaptive reimaging optics. In particular, we simulate observations of spectral line emission from star-forming regions at redshifts z = 0.5 to 2.5 using a variety of spatial sampling scales and give predictions for the signal-to-noise ratio expected as a function of redshift. Using values characteristic of the W.M. Keck II telescope and the new OH-Suppressing InfraRed Imaging Spectrograph (OSIRIS) we calculate integral-field signal-to-noise ratio maps for a sample of U_nGR color-selected star-forming galaxies at redshift z ~ 2 - 2.6 and demonstrate that OSIRIS will be able to reconstruct the two-dimensional projected velocity fields of these galaxies on scales of 100 mas (~ 1 kpc at redshift z ~ 2). With signal-to-noise ratios per spatial sample up to ~ 20, OSIRIS will in some cases be able to distinguish between merger activity and ordered disk rotation. Structures on scales smaller than 1 kpc may be detected by OSIRIS for particularly bright sources, and will be easy targets for future 30m class telescopes.