This paper describes how the Keck Observatory Archive (KOA) is extending open source software components to develop new services. In August 2015, KOA deployed a program interface to discover public data from all instruments equipped with an imaging mode. The interface complies with version 2 of the Simple Imaging Access Protocol (SIAP), under development by the International Virtual Observatory Alliance (IVOA), which defines a standard mechanism for discovering images through spatial queries. The heart of the KOA service is an R-tree-based, database-indexing mechanism prototyped by the Virtual Astronomical Observatory (VAO) and further developed by the Montage Image Mosaic project, designed to provide fast access to large imaging data sets as a first step in creating wide-area image mosaics. The KOA service uses the results of the spatial R-tree search to create an SQLite data database for further relational filtering. The service uses a JSON configuration file to describe the association between instrument parameters and the service query parameters, and to make it applicable beyond the Keck instruments. The R-tree program was itself extended to support temporal (in addition to spatial) indexing, in response to requests from the planetary science community for a search engine to discover observations of Solar System objects. With this 3D-indexing scheme, the service performs very fast time and spatial matches between the target ephemerides, obtained from the JPL SPICE service. Our experiments indicate these matches can be more than 100 times faster than when separating temporal and spatial searches. Images of the tracks of the moving targets, overlaid with the image footprints, are computed with a new command-line visualization tool, mViewer, released with the Montage distribution. The service is currently in test and will be released in Fall 2016.

A collaboration between the W. M. Keck Observatory (WMKO) in Hawaii and the NASA Exoplanet Science Institute (NExScI) in California, the Keck Observatory Archive (KOA) was commissioned in 2004 to archive observing data from WMKO, which operates two classically scheduled 10 m ground-based telescopes. The observing data from Keck is not suitable for direct ingestion into the archive since the metadata contained in the original FITS headers lack the information necessary for proper archiving. Coupled with different standards among instrument builders and the heterogeneous nature of the data inherent in classical observing, in which observers have complete control of the instruments and their observations, the data pose a number of technical challenges for KOA. We describe the methodologies and tools that we have developed to successfully address these difficulties, adding content to the FITS headers and "retrofitting" the metadata in order to support archiving Keck data, especially those obtained before the archive was designed. With the expertise gained from having successfully archived observations taken with all eight currently active instruments at WMKO, we have developed lessons learned from handling this complex array of heterogeneous metadata that help ensure a smooth ingestion of data not only for current but also future instruments, as well as a better experience for the archive user.

The Infrared Processing and Analysis Center (IPAC) and the W. M. Keck Observatory (WMKO) operate an archive for the Keck Observatory. At the end of 2013, KOA completed the ingestion of data from all eight active observatory instruments. KOA will continue to ingest all newly obtained observations, at an anticipated volume of 4 TB per year. The data are transmitted electronically from WMKO to IPAC for storage and curation. Access to data is governed by a data use policy, and approximately two-thirds of the data in the archive are public.

We present the results of the 2MASS Redshift Survey (2MRS), a ten-year project to map the full three-dimensional distribution of galaxies in the nearby Universe. The 2 Micron All-Sky Survey (2MASS) was completed in 2003 and its final data products, including an extended source catalog (XSC), are available on-line. The 2MASS XSC contains nearly a million galaxies with Ks <= 13.5 mag and is essentially complete and mostly unaffected by interstellar extinction and stellar confusion down to a galactic latitude of |b|=5 deg for bright galaxies. Near-infrared wavelengths are sensitive to the old stellar populations that dominate galaxy masses, making 2MASS an excellent starting point to study the distribution of matter in the nearby Universe. We selected a sample of 44,599 2MASS galaxies with Ks <= 11.75 mag and |b|>=5 deg (>= 8 deg towards the Galactic bulge) as the input catalog for our survey. We obtained spectroscopic observations for 11,000 galaxies and used previously-obtained velocities for the remainder of the sample to generate a redshift catalog that is 97.6% complete to well-defined limits and covers 91% of the sky. This provides an unprecedented census of galaxy (baryonic mass) concentrations within 300 Mpc. Earlier versions of our survey have been used in a number of publications that have studied the bulk motion of the Local Group, mapped the density and peculiar velocity fields out to 50 Mpc, detected galaxy groups, and estimated the values of several cosmological parameters. Additionally, we present morphological types for a nearly-complete sub-sample of 20,860 galaxies with Ks < 11.25 mag and |b|>= 10 deg.

NGC 3147, NGC 4698 and 1ES 1927+654 are active galaxies that are classified as Seyfert 2s, based on the line ratios of strong narrow emission lines in their optical spectra. However, they exhibit rapid X-ray spectral variability and/or little indication of obscuration in X-ray spectral fitting, contrary to expectation from the active galactic nucleus (AGN) unification model. Using optical spectropolarimetry with LRIS and near-infrared spectroscopy with NIRSPEC at the W. M. Keck Observatory, we conducted a deep search for hidden polarized broad H-alpha and direct broad Pa-beta or Br-gamma emission lines in these objects. We found no evidence for any broad emission lines from the active nucleus of these galaxies, suggesting that they are unobscured, completely "naked" AGNs that intrinsically lack broad-line regions.

We present 432 low-dispersion optical spectra of 32 Type Ia supernovae (SNe Ia) that also have well-calibrated light curves. The coverage ranges from 6 epochs to 36 epochs of spectroscopy. Most of the data were obtained with the 1.5m Tillinghast telescope at the F. L. Whipple Observatory with typical wavelength coverage of 3700-7400A and a resolution of ~7A. The earliest spectra are thirteen days before B-band maximum; two-thirds of the SNe were observed before maximum brightness. Coverage for some SNe continues almost to the nebular phase. The consistency of the method of observation and the technique of reduction makes this an ideal data set for studying the spectroscopic diversity of SNe Ia.

This paper presents a summary of the first 10 years of operating the Hobby-Eberly Telescope (HET) in queue mode. The scheduling can be quite complex but has worked effectively for obtaining the most science possible with this uniquely designed telescope. The queue must handle dozens of separate scientific programs, the involvement of a number of institutions with individual Telescope Allocation Committees as well as engineering and instrument commissioning. We have continuously revised our queue operations as we have learned from experience. The flexibility of the queue and the simultaneous availability of three instruments, along with a staff trained for all aspects of telescope and instrumentation operation, have allowed optimum use to be made of variable weather conditions and have proven to be especially effective at accommodating targets of opportunity and engineering tasks. In this paper we review the methodology of the HET queue along with its strengths and weaknesses.

We present the reconstructed real-space density and the predicted velocity fields from the Two Mass Redshift Survey (2MRS). The 2MRS is the densest all-sky redshift survey to date and includes about 23,200 galaxies with extinction corrected magnitudes brighter than K = 11.25. Our method is based on the expansion of these fields in Fourier-Bessel functions. Within this framework, the linear redshift distortions only affect the density field in the radial direction and can easily be deconvolved using a distortion matrix. Moreover, in this coordinate system, the velocity field is related to the density field by a simple linear transformation. The shot noise errors in the reconstructions are suppressed by means of a Wiener filter which yields a minimum variance estimate of the density and velocity fields. Using the reconstructed real-space density fields, we identify all major superclusters and voids. At 50 Mpc/h, our reconstructed velocity field indicates a back-side infall to the Great Attractor region of vi = (491 +/- 200)(beta/0.5) km/sec in the Local Group frame and v = (64 +/- 205)(beta/0.5) km/sec in the cosmic microwave background (CMB) frame and beta is the redshift distortion parameter. The direction of the reconstructed dipole agrees well with the dipole derived by Erdogdu et al. (2006). The misalignment between the reconstructed 2MRS and the CMB dipoles drops to 13 degrees at around 5000 km/sec but then increases at larger distances. A version of this paper with high resolution figures can be obtained from http://www.nottingham.ac.uk/~ppzzpe

We estimate the flux weighted acceleration on the Local Group (LG) from the near-infrared Two Micron All Sky Redshift Survey (2MRS). The near-infrared flux weighted dipoles are very robust because they closely approximate a mass weighted dipole, bypassing the effects of redshift distortions and require no preferred reference frame. We use this method with the redshift information to determine the change in dipole with distance. The LG dipole seemingly converges by 60 Mpc/h. Assuming convergence, the comparison of the 2MRS flux dipole and the CMB dipole provides a value for the combination of the mass density and luminosity bias parameters Omega_m^0.6/b_L= 0.40+/-0.09.

We estimate the acceleration on the Local Group (LG) from the Two Micron All Sky Redshift Survey (2MRS). The sample used includes about 23,200 galaxies with extinction corrected magnitudes brighter than K_s=11.25 and it allows us to calculate the flux weighted dipole. The near-infrared flux weighted dipoles are very robust because they closely approximate a mass weighted dipole, bypassing the effects of redshift distortions and require no preferred reference frame. This is combined with the redshift information to determine the change in dipole with distance. The misalignment angle between the LG and the CMB dipole drops to 12 degrees at around 50 Mpc/h, but then increases at larger distances, reaching 21 degrees at around 130 Mpc/h. Exclusion of the galaxies Maffei 1, Maffei 2, Dwingeloo 1, IC342 and M87 brings the resultant flux dipole to 14 degrees away from the CMB velocity dipole In both cases, the dipole seemingly converges by 60 Mpc/h. Assuming convergence, the comparison of the 2MRS flux dipole and the CMB dipole provides a value for the combination of the mass density and luminosity bias parameters Omega_m^0.6/b_L=0.40+/-0.09.

(Abridged) We present spectroscopic and photometric observations of the eclipsing system V1061 Cyg (P = 2.35 days). We show that it is a hierarchical triple in which the third star is visible in the spectra. We combine the radial velocities for the three stars, times of eclipse, and intermediate astrometric data from the HIPPARCOS mission (abscissa residuals) to establish the elements of the outer orbit (P = 15.8 yr) and accurate values for the masses and radii (1-2% errors) and the effective temperatures of the binary components. Both stars are rotating rapidly and have their rotation synchronized with the orbital motion. There are signs of activity including strong X-ray emission and possibly spots. Current stellar evolution models agree well with the properties of the primary, but show a very large discrepancy in the radius of the secondary: the predicted values are about 10% smaller than observed (a 5-sigma effect). Also, the temperature is cooler than predicted by some 200 K. These discrepancies are quite remarkable given that the star is only 7% less massive than the Sun, the calibration point of all stellar models. We identify the chromospheric activity as the likely cause of the effect. Inactive stars agree very well with the models, while active ones such as V1061 Cyg Ab appear systematically too large and too cool. Theory provides an understanding of this in terms of the strong magnetic fields commonly associated with stellar activity, which tend to inhibit convective heat transport.

We report the results of an extensive imaging and spectroscopic survey in the GOODS-North field completed using DEIMOS on the Keck II telescope. Observations of 2018 targets in a magnitude-limited sample of 2911 objects to R=24.4 yield secure redshifts for a sample of 1440 galaxies and AGN plus 96 stars. In addition to redshifts and associated quality assessments, our catalog also includes photometric and astrometric measurements for all targets detected in our R-band imaging survey of the GOODS-North region. We investigate various sources of incompleteness and find the redshift catalog to be 53% complete at its limiting magnitude. The median redshift of z=0.65 is lower than in similar deep surveys because we did not select against low-redshift targets. Comparison with other redshift surveys in the same field, including a complementary Hawaii-led DEIMOS survey, establishes that our velocity uncertainties are as low as 40 km/s for red galaxies and that our redshift confidence assessments are accurate. The distributions of rest-frame magnitudes and colors among the sample agree well with model predictions out to and beyond z=1. We will release all survey data, including extracted 1-D and sky-subtracted 2-D spectra, thus providing a sizable and homogeneous database for the GOODS-North field which will enable studies of large scale structure, spectral indices, internal galaxy kinematics, and the predictive capabilities of photometric redshifts.

We present new photometric observations in the BVRI bands of the double-lined eclipsing binary BD+05 706 conducted over three observing seasons, as well as new X-ray observations obtained with the ROSAT satellite covering a full orbital cycle (P = 18.9 days). A detailed light-curve analysis of the optical data shows the system to be semidetached, confirming indications from an earlier analysis by Torres et al. (1998), with the less massive and cooler star filling its Roche lobe. The system is a member of the rare class of cool Algol systems, which are different from the "classical" Algol systems in that the mass-gaining component is also a late-type star rather than a B- or A-type star. By combining the new photometry with a reanalysis of the spectroscopic observations reported by Torres et al. (1998) we derive accurate absolute masses for the components of M1 = 2.633 +/- 0.028 Msun and M2 = 0.5412 +/- 0.0093 Msun, radii of R1 = 7.55 +/- 0.20 Rsun and R2 = 11.02 +/- 0.21 Rsun, as well as effective temperatures of 5000 +/- 100 K and 4640 +/- 150 K for the primary and secondary, respectively. There are obvious signs of activity (spottedness) in the optical light curve of the binary. Our X-ray light curve clearly shows the primary eclipse but not the secondary eclipse, suggesting that the primary star is the dominant source of the activity in the system. The depth and duration of the eclipse allow us to infer some of the properties of the X-ray emitting region around that star.

We present spectra of the optical transient of GRB021004 obtained with the Hobby-Eberly telescope starting 15.48, 20.31 hours, and 4.84 days after the burst and a spectrum obtained with the H. J. Smith 2.7 m Telescope starting 14.31 hours after the burst. GRB021004 is the first afterglow whose spectrum is dominated by absorption lines from high ionization species with multiple velocity components separated by up to 3000 km/s. We argue that these lines are likely to come from shells around a massive progenitor star. The high velocities and high ionizations arise from a combination of acceleration and flash-ionization by the burst photons and the wind velocity and steady ionization by the progenitor. We also analyze the broad-band spectrum and the light curve. We distinguish six components along the line of sight: (1) The z~2.293 absorption lines arise from the wind of a massive star. For a mass loss rate of ~6 x 10^{-5} solar masses per year, this component also provides the external medium to create the afterglow light. (2) A second shell produces absorption lines with a relative velocity of 560 km/s, and this is associated with the shell created by the fast massive star wind blowing a bubble in the preceding slow wind at a radial distance of order 10 pc. (3) More distant clouds within the host galaxy lie between 30-2500 pc, where they have been ionized by the burst. (4-6) The massive star wind has clumps with radii and over-densities of 0.022, 0.063, and 0.12 parsecs and 50%, 10%, and 10% respectively. The immediate progenitor of the burster could either be a WC-type Wolf-Rayet star or a highly evolved star whose original mass was just too small for it to become a WN-type Wolf-Rayet star.

We use a matched filter algorithm to find and study clusters in both N-body simulations artificially populated with galaxies and the 2MASS survey. In addition to numerous checks of the matched filter algorithm, we present results on the halo multiplicity function and the cluster number function. For a subset of our identified clusters we have information on X-ray temperatures and luminosities which we cross-correlate with optical richness and galaxy velocity dispersions. With all quantities normalized by the spherical radius corresponding to a mass overdensity of Delta_M=200 or the equivalent galaxy number overdensity of Delta_N=200Omega^{-1}~666, we find that the number of L>L_* galaxies in a cluster of mass M_{200} is log N_{*666} = (1.44 +/- 0.17)+(1.10 +/- 0.09)log(M_{200}h/10^{15}Msun) where the uncertainties are dominated by the scatter created by three choices for relating the observed quantities to the cluster mass. The region inside the virial radius has a K-band cluster mass-to-light ratio of (M/L)_K=(116 +/- 46)h which is essentially independent of cluster mass. Integrating over all clusters more massive than M_{200}=10^{14}Msun/h, the virialized regions of clusters contain ~7% of the local stellar luminosity, quite comparable to the mass fraction in such objects in currently popular LambdaCDM models.

We measured the K-band luminosity function using a complete sample of 4192 morphologically-typed 2MASS galaxies with 7 < K < 11.25 mag spread over 2.12 str. Early-type (T < -0.5) and late-type (T > -0.5) galaxies have similarly shaped luminosity functions, alpha_e=-0.92+/-0.10 and alpha_l=-0.87+/-0.09. The early-type galaxies are brighter, M_*e=-23.53+/-0.06 mag compared to M_*l=-22.98\pm0.06 mag, but less numerous, n_*e=(0.0045+/-0.0006)h^3/Mpc^3 compared to n_*l=(0.0101+/-0.0013)h^3/Mpc^3 for H_0=100h km/s Mpc, such that the late-type galaxies slightly dominate the K-band luminosity density, j_late/j_early=1.17+/-0.12. Our morphological classifications are internally consistent, consistent with previous classifications and lead to luminosity functions unaffected by the estimated uncertainties in the classifications. These luminosity functions accurately predict the K-band number counts and redshift distributions for K < 18 mag, beyond which the results depend on galaxy evolution and merger histories.

The Two Micron All Sky Survey (2MASS) is now underway and will provide a complete census of galaxies as faint as 13.5 mag (3 mJy) at 2.2 microns for most of the sky, and ~12.1 mag (10 mJy) for regions veiled by the Milky Way. This census has already discovered nearby galaxies previously hidden behind our Galaxy, and will allow delineation of large-scale structures in the distribution of galaxies across the whole sky. Here we report the detection and discovery of new extended sources from this survey for fields incorporating the Galactic plane at longitudes between 40 and 70 deg. The area-normalized detection rate is ~1-2 galaxies per deg^2 brighter than 12.1 mag (10 mJy), roughly constant with Galactic latitude throughout the "Zone of Avoidance", of which 85-95% are newly discovered sources. In conjunction with the deep HI surveys, 2MASS will greatly increase the current census of galaxies hidden behind the Milky Way.

In the construction of an X-ray selected sample of galaxy clusters for cosmological studies, we have assembled a sample of 495 X-ray sources found to show extended X-ray emission in the first processing of the ROSAT All-Sky Survey. The sample covers the celestial region with declination $\delta \ge 0\deg $ and galactic latitude $|b_{II}| \ge 20\deg $ and comprises sources with a count rate $\ge 0.06$ counts s$^{-1}$ and a source extent likelihood of 7. In an optical follow-up identification program we find 378 (76%) of these sources to be clusters of galaxies. ...