Keck II Research Articles

The TESS–Keck Survey. XIII. An Eccentric Hot Neptune with a Similar-mass Outer Companion around TOI-1272

We report the discovery of an eccentric hot Neptune and a non-transiting outer planet around TOI-1272. We identified the eccentricity of the inner planet, with an orbital period of 3.3 days and R p,b = 4.1 ± 0.2 R ⊕, based on a mismatch between the observed transit duration and the expected duration for a circular orbit. Using ground-based radial velocity (RV) measurements from the HIRES instrument at the Keck Observatory, we measured the mass of TOI-1272b to be M p,b = 25 ± 2 M ⊕. We also confirmed a high eccentricity of e b = 0.34 ± 0.06, placing TOI-1272b among the most eccentric well-characterized sub-Jovians. We used these RV measurements to also identify a non-transiting outer companion on an 8.7 day orbit with a similar mass of M p,c sin i = 27 ± 3 M ⊕ and e c ≲ 0.35. Dynamically stable planet–planet interactions have likely allowed TOI-1272b to avoid tidal eccentricity decay despite the short circularization timescale expected for a close-in eccentric Neptune. TOI-1272b also maintains an envelope mass fraction of f env ≈ 11% despite its high equilibrium temperature, implying that it may currently be undergoing photoevaporation. This planet joins a small population of short-period Neptune-like planets within the “Hot Neptune Desert” with a poorly understood formation pathway.

An AGN with an Ionized Gas Outflow in a Massive Quiescent Galaxy in a Protocluster at z = 3.09

We report the detection of an ionized gas outflow from an X-ray active galactic nucleus hosted in a massive quiescent galaxy in a protocluster at z = 3.09 (J221737.29+001823.4). It is a type-2 QSO with broad (W 80 > 1000 km s−1) and strong (/erg s−1) ≈ 43.4) [O iii]λ λ 4959,5007 emission lines detected by slit spectroscopy in three-position angles using Multi-Object Infra-Red Camera and Spectrograph (MOIRCS) on the Subaru telescope and the Multi-Object Spectrometer For Infra-Red Exploration (MOSFIRE) on the Keck-I telescope. In the all slit directions, [O iii] emission is extended to ∼15 physical kpc and indicates a powerful outflow spreading over the host galaxy. The inferred ionized gas mass outflow rate is 22 ± 3 M ⊙ yr−1. Although it is a radio source, according to the line diagnostics using Hβ, [O ii], and [O iii], photoionization by the central QSO is likely the dominant ionization mechanism rather than shocks caused by radio jets. On the other hand, the spectral energy distribution of the host galaxy is well characterized as a quiescent galaxy that has shut down star formation several hundred Myr ago. Our results suggest a scenario that QSOs are powered after the shutdown of the star formation and help complete the quenching of massive quiescent galaxies at high redshift.

Regularized 3D spectroscopy with CubeFit: Method and application to the Galactic Center circumnuclear disk

Context. The Galactic Center black hole and the nuclear star cluster are surrounded by a clumpy ring of gas and dust, the circumnuclear disk (CND), that rotates about them at a standoff distance of ≃1.5 pc. The mass and density of individual clumps in the CND are disputed. Aims. We seek to use H2 to characterize the clump size distribution and to investigate the morphology and dynamics of the interface between the ionized interior layer of the CND and the molecular reservoir lying farther out (corresponding to the inner rim of the CND, illuminated in ultraviolet light by the central star cluster). Methods. We have observed two fields of approximately 20″ × 20″ in the CND at near-infrared wavelengths with the OSIRIS spectro-imager at the Keck Observatory. These two fields, located at the approaching and receding nodes of the CND, best display this interface. Our data cover two H2 lines as well as the Brγ line (tracing H II). We have developed the tool CubeFit, an original method for extracting maps of continuous physical parameters (such as the velocity field and velocity dispersion) from integral-field spectroscopy data, using regularization to largely preserve spatial resolution in regions of low signal-to-noise ratio. Results. This original method enables us to isolate compact, bright features in the interstellar medium of the CND. Several clumps in the southwestern field assume the appearance of filaments, many of which are parallel to one another. We conclude that these clumps cannot be self-gravitating.

The Dwarf Galaxy Population at z ∼ 0.7: A Catalog of Emission Lines and Redshifts from Deep Keck Observations

We present a catalog of spectroscopically measured redshifts over 0 < z < 2 and emission-line fluxes for 1440 galaxies. The majority (∼65%) of the galaxies come from the HALO7D survey, with the remainder from the DEEPwinds program. This catalog includes redshifts for 646 dwarf galaxies with . Eight-hundred and ten catalog galaxies did not have previously published spectroscopic redshifts, including 454 dwarf galaxies. HALO7D used the DEIMOS spectrograph on the Keck II telescope to take very deep (up to 32 hr exposure, with a median of ∼7 hr) optical spectroscopy in the COSMOS, EGS, GOODS-North, and GOODS-South CANDELS fields, and in some areas outside CANDELS. We compare our redshift results to existing spectroscopic and photometric redshifts in these fields, finding only a 1% rate of discrepancy with other spectroscopic redshifts. We measure a small increase in median photometric redshift error (from 1.0% to 1.3%) and catastrophic outlier rate (from 3.5% to 8%) with decreasing stellar mass. We obtained successful redshift fits for 75% of massive galaxies, and demonstrate a similar 70%–75% successful redshift measurement rate in galaxies, suggesting similar survey sensitivity in this low-mass range. We describe the redshift, mass, and color–magnitude distributions of the catalog galaxies, finding HALO7D galaxies representative of CANDELS galaxies up to i-band magnitudes of 25. The catalogs presented will enable studies of star formation, the mass–metallicity relation, star formation–morphology relations, and other properties of the z ∼ 0.7 dwarf galaxy population.

AIROPA III: testing simulated and on-sky data

Adaptive optics (AO) images from the W. M. Keck Observatory have delivered numerous influential scientific results, including detection of multi-system asteroids, the supermassive black hole at the center of the Milky Way, and directly imaged exoplanets. Specifically, the precise and accurate astrometry these images yield was used to measure the mass of the supermassive black hole using orbits of the surrounding star cluster. Despite these successes, one of the major obstacles to improved astrometric measurements is the spatial and temporal variability of the point-spread function delivered by the instruments. Anisoplanatic and Instrumental Reconstruction of Off-axis PSFs for AO (AIROPA) is a software package for the astrometric and photometric analysis of AO images using point-spread function fitting together with the technique of point-spread function reconstruction. In AO point-spread function reconstruction, the knowledge of the instrument performance and of the atmospheric turbulence is used to predict the long-exposure point-spread function of an observation. We present the results of our tests using AIROPA on both simulated and on-sky images of the Galactic Center. We find that our method is very reliable in accounting for the static aberrations internal to the instrument, but it does not improve significantly the accuracy on sky, possibly due to uncalibrated telescope aberrations.

Kinematics and Metallicity of Red Giant Branch Stars in the Northeast Shelf of M31* *The data presented herein were obtained at the W. M. Keck Observatory, which is operated as a scientific partnership among the California Institute of Technology, the University of California, and the National Aeronautics and Space Administration. The Observatory was made possible by the generous financial support

We obtained Keck/DEIMOS spectra of 556 individual red giant branch stars in four spectroscopic fields spanning 13−31 projected kpc along the northeast (NE) shelf of M31. We present the first detection of a complete wedge pattern in the space of projected M31-centric radial distance versus line-of-sight velocity for this feature, which includes the returning stream component of the shelf. This wedge pattern agrees with expectations of a tidal shell formed in a radial merger and provides strong evidence in favor of predictions of Giant Stellar Stream (GSS) formation models in which the NE shelf originates from the second orbital wrap of the tidal debris. The observed concentric wedge patterns of the NE, west (W), and southeast (SE) shelves corroborate this interpretation independently of the models. We do not detect a kinematical signature in the NE shelf region corresponding to an intact progenitor core, favoring GSS formation models in which the progenitor is completely disrupted. The shelf’s photometric metallicity ([Fe/H]phot) distribution implies that it is dominated by tidal material, as opposed to the phase-mixed stellar halo or the disk. The metallicity distribution ([Fe/H]phot = −0.42 ± 0.01) also matches the GSS, and consequently the W and SE shelves, further supporting a direct physical association between the tidal features.

An investigation of the magnetic activity of HD 134319 based on TESS photometry and ground-based spectroscopy

ABSTRACT We present an analysis of the starspot properties and chromospheric activity on HD 134319 using high-precision photometry by Transiting Exoplanet Survey Satellite (TESS) in sectors 14–16 (T1) and 21–23 (T2) and high-resolution spectroscopy by the ELODIE spectrograph on 1.93-m telescope of Observatoire de Haute-Provence (OHP/ELODIE) and the High Resolution Echelle Spectrometer on Keck 10-m telescope of W. M. Keck Observatory (Keck/HIRES) during the years 1995–2013. We applied a two-spot model with generalized Lomb–Scargle (GLS) periodogram determined period of P = 4.436391 ± 0.00137 d to model chunks sliding over TESS light curve. We also measured the relative equivalent widths of Ca ii H and K, Hβ and Hα emissions by subtracting the overall spectrum from individual spectra. It was found that a two-spot configuration, namely a primary, slowly evolving and long-lasting spot (P) plus a secondary, rapidly evolving spot (S), was capable of explaining the data, although the actual starspot distribution cannot be derived from the collected data. Despite the spot radius–latitude degeneracy revealed in the best-fitting solutions, a sudden alternation between P and S radii followed by a gradual decrease of S in T1 and a decrease of both P and S from T1 to T2 were significant, corresponding to the evolution of magnetic activity. In addition, S revealed rotation and an oscillatory longitude migration synchronized to P in T1, but had a much larger migration than P in T2. This might indicate the evolution of the internal magnetic configuration. Chromospheric activity indicators were found to be tightly correlated with each other and revealed rotational modulation as well as a long-term decrease of emissions, implying the existence and evolution of magnetic activity on HD 134319.

Discovery of three new near-pristine absorption clouds at z = 2.6–4.4

ABSTRACT We report the discovery of three new ‘near-pristine’ Lyman limit systems (LLSs), with metallicities ≈1/1000 solar, at redshifts 2.6, 3.8, and 4.0, with a targeted survey at the Keck Observatory. High-resolution echelle spectra of eight candidates yielded precise column densities of hydrogen and weak, but clearly detected, metal lines in seven LLSs; we previously reported the one remaining, apparently metal-free LLS, to have metallicity &amp;lt;1/10 000 solar. Robust photoionization modelling provides metallicities [Si/H] = −3.05 to −2.94, with 0.26 dex uncertainties (95 per cent confidence) for three LLSs, and $\textrm {[Si/H]} \gtrsim -2.5$ for the remaining four. Previous simulations suggest that near-pristine LLSs could be the remnants of PopIII supernovae, so comparing their detailed metal abundances with nucleosynthetic yields from supernovae models is an important goal. Unfortunately, at most two different metals were detected in each new system, despite their neutral hydrogen column densities ($10^{19.2}\textrm {--}10^{19.4}\, \textrm {cm}^{-2}$) being two orders of magnitude larger than the two previous, serendipitously discovered near-pristine LLSs. Nevertheless, the success of this first targeted survey for near-pristine systems demonstrates the prospect that a much larger, future survey could identify clear observational signatures of PopIII stars. With a well-understood selection function, such a survey would also yield the number density of near-pristine absorbers that, via comparison to future simulations, could reveal the origin(s) of these rare systems.

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.

On-sky Reconstruction of Keck Primary Mirror Piston Offsets Using a Zernike Wavefront Sensor

The next generation of large ground- and space-based optical telescopes will have segmented primary mirrors. Co-phasing the segments requires a sensitive wavefront sensor capable of measuring phase discontinuities. The Zernike wavefront sensor (ZWFS) is a passive wavefront sensor that has been demonstrated to sense segmented-mirror piston, tip, and tilt with picometer precision in laboratory settings. We present the first on-sky results of an adaptive optics fed ZWFS on a segmented aperture telescope, W.M. Keck Observatory's Keck II. Within the Keck Planet Imager and Characterizer light path, the ZWFS mask operates in the H band using an InGaAs detector (CRED2). We piston segments of the primary mirror by a known amount and measure the mirror's shape using both the ZWFS and a phase retrieval method on data acquired with the facility infrared imager, NIRC2. In the latter case, we employ slightly defocused NIRC2 images and a modified Gerchberg–Saxton phase retrieval algorithm to estimate the applied wavefront error. We find good agreement when comparing the phase retrieval and ZWFS reconstructions, with average measurements of 408 ± 23 and 394 ± 46 nm, respectively, for three segments pistoned by 400 nm of optical path difference (OPD). Applying various OPDs, we find that we are limited to ∼100 nm OPD of applied piston, due to insufficient averaging of the adaptive optics residuals of our observations. We also present simulations of the ZWFS that help to explain the systematic offset observed in the ZWFS reconstructed data.

Analyzing long-term performance of the Keck-II adaptive optics system

We present an analysis of the long-term performance of the W. M. Keck observatory laser guide star adaptive optics (LGS-AO) system and explore factors that influence the overall AO performance most strongly. Astronomical surveys can take years or decades to finish, so it is worthwhile to characterize the AO performance on such timescales in order to better understand future results. The Keck telescopes have two of the longest-running LGS-AO systems in use today, and as such they represent an excellent test-bed for processing large amounts of AO data. We use a Keck-II near infrared camera 2 (NIRC2) LGSAO surve of the Galactic Center (GC) from 2005 to 2019 for our analysis, combining image metrics with AO telemetry files, multiaperture scintillation sense/differential imaging motion monitor turbulence profiles, seeing information, weather data, and temperature readings in a compiled dataset to highlight areas of potential performance improvement. We find that image quality trends downward over time, despite multiple improvements made to Keck-II and its AO system, resulting in a 9 mas increase in the average full width at half maximum (FWHM) and a 3% decrease in the average Strehl ratio over the course of the survey. Image quality also trends upward with ambient temperature, possibly indicating the presence of uncorrected turbulence in the beam path. Using nine basic features from our dataset, we train a simple machine learning (ML) algorithm to predict the delivered image quality of NIRC2 given current atmospheric conditions, which could eventually be used for real-time observation planning and exposure time adjustments. A random forest algorithm trained on this data can predict the Strehl ratio of an image to within 18% and the FWHM to within 7%, which is a solid baseline for future applications involving more advanced ML techniques. The assembled dataset and coding tools are released to the public as a resource for testing new predictive control and point spread function-reconstruction algorithms.

SRG/ART-XC discovery of SRGA J204318.2+443815: Towards the complete population of faint X-ray pulsars

We report the discovery of the new long-period X-ray pulsar SRGA J204318.2+443815/SRGe J204319.0+443820 in a Be binary system. The source was found in the second all-sky survey by the Mikhail Pavlinsky ART-XC telescope on board the SRG mission. The follow-up observations with XMM-Newton, NICER, and NuSTAR allowed us to discover a strong coherent signal in the source light curve with a period of ~742 s. The pulsed fraction was found to depend on an increase in energy from ~20% in soft X-rays to &gt;50% at high energies, as is typical for X-ray pulsars. The source has a quite hard spectrum with an exponential cutoff at high energies and a bolometric luminosity of LX ≃ 4 × 1035 ergs-1. The X-ray position of the source is found to be consistent with the optical transient ZTF18abjpmzf, located at a distance of ~8.0 kpc. Dedicated optical and infrared observations with the RTT-150, NOT, Keck, and Palomar telescopes revealed a number of emission lines (Ha, He I, and the Paschen and Braket series) with a strongly absorbed continuum. According to the SRG scans and archival XMM-Newton data, the source flux is moderately variable (by a factor of 4-10) on timescales of several months and years. All this suggests that SRGA J204318.2+443815/SRGe J204319.0+443820 is a new quasi-persistent low-luminosity X-ray pulsar in a distant binary system with a Be-star of the B0-B2e class. Thus the SRG observatory allowed us to unveil a hidden population of faint objects, including a population of slowly rotating X-ray pulsars in Be systems.

Searching for Islands of Reionization: A Potential Ionized Bubble Powered by a Spectroscopic Overdensity at z = 8.7

We present the results from a spectroscopic survey using the MOSFIRE near-infrared spectrograph on the 10 m Keck telescope to search for Lyα emission from candidate galaxies at z ∼ 9–10 in four of the CANDELS fields (GOODS-N, EGS, UDS, and COSMOS). We observed 11 target galaxies, detecting Lyα from one object in ∼8.1 hr of integration, at z = 8.665 ± 0.001 with an integrated signal-to-noise ratio > 7. This galaxy is in the CANDELS Extended Groth Strip (EGS) field and lies physically close (3.5 physical Mpc [pMpc]) to another confirmed galaxy in this field with Lyα detected at z = 8.683. The detection of Lyα suggests the existence of large (∼1 pMpc) ionized bubbles fairly early in the reionization process. We explore the ionizing output needed to create bubbles of this size at this epoch and find that such a bubble requires more than the ionizing power provided by the full expected population of galaxies (by integrating the UV luminosity function down to M UV = −13). The Lyα we detect would be able to escape the predominantly neutral intergalactic medium at this epoch if our detected galaxy is inhabiting an overdensity, which would be consistent with the photometric overdensity previously identified in this region by Finkelstein et al. This implies that the CANDELS EGS field is hosting an overdensity at z = 8.7 that is powering one or more ionized bubbles, a hypothesis that will be imminently testable with forthcoming James Webb Space Telescope observations in this field.

Scientific vs. Social Media Impact in Astronomy: Study of Gemini, Keck, and PRL Publications

This study was carried out to assess whether a correlation exists between the social media visibility (altmetrics) and scientific impact (citations) of published refereed articles in astronomy. We have compiled citations and altmetrics of refereed articles between 2016-2019 at Gemini Observatory/NSF’s NOIRLab (Gemini), W.M. Keck Observatory (Keck) and the Physical Research Laboratory (PRL). The study was aimed to identify the impact of social media, if any, on academic recognition by peers. This collaborative study identifies the existence of a positive correlation between the citations and altmetrics and also explores the characteristics of the correlation for the publications from the three institutions.

SHARP – VIII. J0924+0219 lens mass distribution and time-delay prediction through adaptive-optics imaging

ABSTRACT Strongly lensed quasars can provide measurements of the Hubble constant (H0) independent of any other methods. One of the key ingredients is exquisite high-resolution imaging data, such as Hubble Space Telescope (HST) imaging and adaptive-optics (AO) imaging from ground-based telescopes, which provide strong constraints on the mass distribution of the lensing galaxy. In this work, we expand on the previous analysis of three time-delay lenses with AO imaging (RX J1131−1231, HE 0435−1223, and PG 1115+080), and perform a joint analysis of J0924+0219 by using AO imaging from the Keck telescope, obtained as part of the Strong lensing at High Angular Resolution Program (SHARP) AO effort, with HST imaging to constrain the mass distribution of the lensing galaxy. Under the assumption of a flat Λ cold dark matter (ΛCDM) model with fixed Ωm = 0.3, we show that by marginalizing over two different kinds of mass models (power-law and composite models) and their transformed mass profiles via a mass-sheet transformation, we obtain $\Delta t_{\rm BA}=6.89\substack{+0.8\\-0.7}\, h^{-1}\hat{\sigma }_{v}^{2}$ d, $\Delta t_{\rm CA}=10.7\substack{+1.6\\-1.2}\, h^{-1}\hat{\sigma }_{v}^{2}$ d, and $\Delta t_{\rm DA}=7.70\substack{+1.0\\-0.9}\, h^{-1}\hat{\sigma }_{v}^{2}$ d, where $h=H_{0}/100\,\rm km\, s^{-1}\, Mpc^{-1}$ is the dimensionless Hubble constant and $\hat{\sigma }_{v}=\sigma ^{\rm ob}_{v}/(280\,\rm km\, s^{-1})$ is the scaled dimensionless velocity dispersion. Future measurements of time delays with 10 per cent uncertainty and velocity dispersion with 5 per cent uncertainty would yield a H0 constraint of ∼15 per cent precision.

Kinematics and Feedback in H ii Regions in the Dwarf Starburst Galaxy IC 10

We present a survey of the central region of the nearest starburst galaxy, IC 10, using the W. M. Keck Observatory Keck Cosmic Web Imager (KCWI) at high spectral and spatial resolution. We map the central starburst of IC 10 to sample the kinematic and ionization properties of the individual star-forming regions. Using the low spectral resolution mode of KCWI, we map the oxygen abundance, and with the high spectral resolution mode, we identify 46 individual H II regions. These H II regions have an average radius of 4.0 pc, star formation rate ∼1.3 × 10−4 M ⊙ yr−1, and velocity dispersion ∼16 km s−1. None of the H II regions appear to be virialized (α vir ≫ 1), and on average, they show evidence of ongoing expansion. IC 10's H II regions are offset from the star-forming-region size–luminosity scaling relationships, as well as Larson’s Law that relates size and velocity dispersion. We investigate the balance of inward and outward pressure, P in and P out, finding P out > P in in 89% of H II regions, indicating feedback-driven expansion even in these low-mass H II regions. We find warm gas pressure (P gas) provides the dominant contribution to the outward pressure (P out). This counteracts the inward pressure, which is dominated by turbulence in the surrounding gas rather than self-gravity. Five H II regions show evidence of outflows that are most likely supported by either stellar winds (two regions) or champagne flows (three regions). These observations provide new insights into the state of the star-forming regions in IC 10 and negative feedback from low-mass clusters.

Checking the possibility of determining the relative orbits of stars rotating around the center body of the Galaxy

Abstract The results of observations of selected stars orbiting the central body of our Galaxy, which, according to the latest data, has a mass of about 4 million of solar mass, are considered. Using the methods of observation and processing of binary stars used at the Pulkovo Observatory and the published exact relative positions of the star-central body obtained with the Keck and VLT telescopes, we checked the possibility of determining the preliminary orbits of the above objects in some cases from observations covering sufficient short arc in relation to the turnover period. Examples are the stars S02 and S102 with published positions and the graphical data for S27 and S4711. We attempted to calculate the assumed relativistic effects in accordance with the orbital elements obtained by Pulkovo methods. In this case, the estimates of the O–C errors for these orbits turned out to be comparable with the O–C errors for the orbits published by the Keck and VLT observers, in some cases about several milliarcseconds. The theoretical relativistic parameters also turned out to be of the same order of value with the results of these authors.

Saturn's Weather‐Driven Aurorae Modulate Oscillations in the Magnetic Field and Radio Emissions

AbstractThe Cassini spacecraft revealed that Saturn's magnetic field displayed oscillations at a period originally thought to match the planetary rotation rate but later found not to. One of many proposed theories predicts that a polar twin‐cell neutral weather system drives this variation, producing observable differences in flows within Saturn's ionosphere. Here, using spectral observations of auroral emission lines taken by the Keck Observatory's Near Infrared Echelle Spectrograph (Keck‐NIRSPEC) in 2017, we derive ion line‐of‐sight velocity maps after grouping spectra into rotational quadrants matching phases of the planetary magnetic field. We measure 0.5 km s−1 wind systems in the ionosphere consistent with predicted neutral twin‐vortex flow patterns. These findings demonstrate that neutral winds in Saturn's polar regions cause the rotational period, as determined via the magnetic field, to exhibit differences and time variabilities relative to the planet's true period of rotation in a process never before seen within planetary atmospheres.

Dynamical Mass of the Exoplanet Host Star HR 8799

Abstract HR 8799 is a young A5/F0 star hosting four directly imaged giant planets at wide separations (∼16–78 au), which are undergoing orbital motion and have been continuously monitored with adaptive optics imaging since their discovery over a decade ago. We present a dynamical mass of HR 8799 using 130 epochs of relative astrometry of its planets, which include both published measurements and new medium-band 3.1 μm observations that we acquired with NIRC2 at Keck Observatory. For the purpose of measuring the host-star mass, each orbiting planet is treated as a massless particle and is fit with a Keplerian orbit using Markov chain Monte Carlo. We then use a Bayesian framework to combine each independent total mass measurement into a cumulative dynamical mass using all four planets. The dynamical mass of HR 8799 is 1.47 − 0.17 + 0.12 M ⊙ assuming a uniform stellar mass prior, or 1.46 − 0.15 + 0.11 M ⊙ with a weakly informative prior based on spectroscopy. There is a strong covariance between the planets’ eccentricities and the total system mass; when the constraint is limited to low-eccentricity solutions of e &lt; 0.1, which are motivated by dynamical stability, our mass measurement improves to 1.43 − 0.07 + 0.06 M ⊙. Our dynamical mass and other fundamental measured parameters of HR 8799 together with Modules for Experiments in Stellar Astrophysics Isochrones and Stellar Tracks grids yields a bulk metallicity most consistent with [Fe/H] ∼ −0.25–0.00 dex and an age of 10–23 Myr for the system. This implies hot-start masses of 2.7–4.9 M Jup for HR 8799 b and 4.1–7.0 M Jup for HR 8799 c, d, and e, assuming they formed at the same time as the host star.

Volatile Abundances, Extended Coma Sources, and Nucleus Ice Associations in Comet C/2014 Q2 (Lovejoy)

Abstract High-resolution infrared spectra of comet C/2014 Q2 Lovejoy were acquired with NIRSPEC at the W. M. Keck Observatory on two post-perihelion dates (UT 2015 February 2 and 3). H2O was measured simultaneously with CO, CH3OH, H2CO, CH4, C2H6, C2H4, C2H2, HCN, and NH3 on both dates, and rotational temperatures, production rates, relative abundances, H2O ortho-to-para ratios, and spatial distributions in the coma were determined. The first detection of C2H4 in a comet from ground-based observations is reported. Abundances relative to H2O for all species were found to be in the typical range compared with values for other comets in the overall population to date. There is evidence of variability in rotational temperatures and production rates on timescales that are small compared with the rotational period of the comet. Spatial distributions of volatiles in the coma suggest complex outgassing behavior. CH3OH, HCN, C2H6, and CH4 spatial distributions in the coma are consistent with direct release from associated ices in the nucleus and are peaked in a more sunward direction compared with co-measured dust. H2O spatial profiles are clearly distinct from these other four species, likely due to a sizable coma contribution from icy grain sublimation. Spatial distributions for C2H2, H2CO, and NH3 suggest substantial contributions from extended coma sources, providing further evidence for distinct origins and associations for these species in comets. CO shows a different spatial distribution compared with other volatiles, consistent with jet activity from discrete nucleus ice sources.