Search

The Keck Planet Imager and Characterizer (KPIC) is a series of upgrades for the Keck II Adaptive Optics (AO) system and the NIRSPEC spectrograph to enable diffraction limited, high resolution (R>30000) spectroscopy of exoplanets and low mass companions in the K and L bands. Phase I consisted of single mode fiber injection/extraction units (FIU/FEU) used in conjunction with a H band pyramid wavefront sensor. The use of single mode fibers provides a gain in stellar rejection, a substantial reduction in sky background, and an extremely stable line spread function in the spectrograph. Phase II, deployed and commissioned in 2022, brought a 1000 actuator deformable mirror, beam shaping optics, a vortex mask, and other upgrades to the FIU/FEU. An additional service mission in 2024 extended operations down to y band, delivered an atmospheric dispersion corrector, and provided access to two laser frequency combs. KPIC phase II brings higher planet throughput, lower stellar leakage and many new observing modes which extend its ability to characterize exoplanets at high spectral resolution, building on the success of phase I. In this paper we present a description of the final phase II version of KPIC, along with results of system level laboratory testing and characterization showing the instrument's phase II throughput, stability, repeatability, and other key performance metrics prior to delivery and during installation at Keck. We outlined the capabilities of the various observing modes enabled by the new modules as well as efforts to compensate for static aberrations and non common path errors at Keck, which were issues that plagued phase I. Finally, we show results from commissioning., Comment: 36 pages, 16 figures

We present an atmospheric characterization and orbital analysis of HD 206893 B, an exceptionally red, L/T-transition substellar companion in a multiplanetary system, via Keck Planet Imager and Characterizer (KPIC) high-resolution (R $\sim$ 35,000) K-band spectroscopy. Using PHOENIX atmospheric models in a forward-model framework that fits the spectrum of the companion and diffracted starlight simultaneously, we detect HD 206893 B at $>8\sigma$ significance via cross-correlation in two epochs. We find an effective temperature for the companion of $1634^{+72}_{-38}$ K and a log(g) of $4.55^{+0.17}_{-0.22}$. Only accounting for statistical uncertainties, we measure the carbon-oxygen ratio (C/O) of this companion to be $0.57 \pm 0.02$, or near-solar while assuming solar metallicity. The C/O ratio we measure fits the tentative trend of $>4 M_{Jup}$ companions having near-solar C/O ratios while less massive companions have greater-than-solar C/O ratios. Using substellar evolution models, we find an age of $112^{+36}_{-22}$ Myr, a mass of $22.7^{+2.5}_{-1.7} M_{Jup}$, and a radius of $1.11 \pm 0.03 R_{Jup}$ for this companion. We also use KPIC radial velocity data to fit the orbit of HD 206893 B and analyze the orbital stability of this system. We find that the orbital stability is relatively independent of the mass of HD 206893 B, and favors an orbital configuration where B and its interior planetary companion, HD 206893 c, are co-planar. The measured C/O ratio coupled with the current architecture of the system cannot rule out a core accretion scenario, nor a disk fragmentation scenario regarding the formation pathway of HD 206893 B., Comment: 37 pages, 23 figures

We present Keck/KPIC phase II $K$-band observations of the non-transiting hot Jupiter HD 143105 b. Using a cross-correlation approach, we make the first detection of the planetary atmosphere at $K_p = 185^{+11}_{-13}\rm km\ s^{-1}$ and an inferior conjunction time 2.5 hours before the previously-published ephemeris. The retrieved $K_p$ value, in combination with orbital period, mass of the host star, and lack of transit detection, gives an orbital inclination of $78^{\circ+2}_{-12}$ and a true planet mass of 1.23$\pm0.10\rm\ M_J$. While the equilibrium temperature of HD 143105 b is in the transition regime between non-inverted and inverted atmospheres, our analysis strongly prefers a non-inverted atmosphere. Retrieval analysis indicates the atmosphere of HD 143105 b is cloud-free to approximately 1 bar and dominated by H$_2$O absorption ($\log \rm H_2O_{MMR} = -3.9^{+0.8}_{-0.5}$), placing only an upper limit on the CO abundance ($\log \rm CO_{MMR} < -3.7$ at 95% confidence). We place no constraints on the abundances of Fe, Mg, or $^{13}$CO. From these abundances, we place an upper limit on the carbon-to-oxygen ratio for HD 143105 b, $\rm C/O < 0.2$ at 95% confidence, and find the atmospheric metallicity is approximately $0.1\times$ solar. The low metallicity may be responsible for the lack of a thermal inversion, which at the temperature of HD 143105 b would likely require significant opacity from TiO and/or VO. With these results, HD 143105 b joins the small number of non-transiting hot Jupiters with detected atmospheres., Comment: 19 pages, 7 figures, 2 tables. Accepted in AJ

The $\sim$5 Myr PDS 70 is the only known system with protoplanets residing in the cavity of the circumstellar disk from which they formed, ideal for studying exoplanet formation and evolution within its natal environment. Here we report the first spin constraint and C/O measurement of PDS 70b from Keck/KPIC high-resolution spectroscopy. We detected CO (3.8 $\sigma$) and H$_2$O (3.5 $\sigma$) molecules in the PDS 70b atmosphere via cross-correlation, with a combined CO and H$_2$O template detection significance of 4.2 $\sigma$. Our forward model fits, using BT-Settl model grids, provide an upper limit for the spin-rate of PDS 70b ($<$29 km s$^{-1}$). The atmospheric retrievals constrain the PDS 70b C/O ratio to ${0.28}^{+0.20}_{-0.12}$ ($<$0.63 under 95$\%$ confidence level) and a metallicity [C/H] of ${-0.2}^{+0.8}_{-0.5}$ dex, consistent with that of its host star. The following scenarios can explain our measured C/O of PDS 70b in contrast with that of the gas-rich outer disk (for which C/O $\gtrsim$ 1). First, the bulk composition of PDS 70b might be dominated by dust+ice aggregates rather than disk gas. Another possible explanation is that the disk became carbon-enriched $\textit{after}$ PDS 70b was formed, as predicted in models of disk chemical evolution and as observed in both very low mass star and older disk systems with $\textit{JWST}$/MIRI. Because PDS 70b continues to accrete and its chemical evolution is not yet complete, more sophisticated modeling of the planet and the disk, and higher quality observations of PDS 70b (and possibly PDS 70c), are necessary to validate these scenarios., Comment: Accepted to ApJ Letters; 13 pages, 3 figures

We present several nonlinear wavefront sensing techniques for few-mode sensors, all of which are empirically calibrated and agnostic to the choice of wavefront sensor. The first class of techniques involves a straightforward extension of the linear phase retrieval scheme to higher order; the resulting Taylor polynomial can then be solved using the method of successive approximations, though we discuss alternate methods such as homotopy continuation. In the second class of techniques, a model of the WFS intensity response is created using radial basis function interpolation. We consider both forward models, which map phase to intensity and can be solved with nonlinear least-squares methods such as the Levenberg-Marquardt algorithm, as well as backwards models which directly map intensity to phase and do not require a solver. We provide demonstrations for both types of techniques in simulation using a quad-cell sensor and a photonic lantern wavefront sensor as examples. Next, we demonstrate how the nonlinearity of an arbitrary sensor may studied using the method of numerical continuation, and apply this technique both to the quad-cell sensor and a photonic lantern sensor. Finally, we briefly consider the extension of nonlinear techniques to polychromatic sensors., Comment: Accepted to Applied Optics

Brown dwarfs with measured dynamical masses and spectra from direct imaging are benchmarks that anchor substellar atmosphere cooling and evolution models. We present Subaru SCExAO/CHARIS infrared spectroscopy of HIP 93398 B, a brown dwarf companion recently discovered by Li et al. 2023 as part of an informed survey using the Hipparcos-Gaia Catalog of Accelerations. This object was previously classified as a T6 dwarf based on its luminosity, with its independently-derived age and dynamical mass in tension with existing models of brown dwarf evolution. Spectral typing via empirical standard spectra, temperatures derived by fitting substellar atmosphere models, and J-H, J-K and H-L' colors all suggest that this object has a substantially higher temperature and luminosity, consistent with classification as a late-L dwarf near the L/T transition (T = 1200$^{+140}_{-119}$ K) with moderate to thick clouds possibly present in its atmosphere. When compared with the latest generation of evolution models that account for clouds with our revised luminosity and temperature for the object, the tension between the model-independent mass/age and model predictions is resolved., Comment: Accepted to AJ, 19 pages

Spectroastrometry, which measures wavelength-dependent shifts in the center of light, is well-suited for studying objects whose morphology changes with wavelength at very high angular resolutions. Photonic lantern (PL)-fed spectrometers have potential to enable measurement of spectroastrometric signals because the relative intensities between the PL output SMFs contain spatial information on the input scene. In order to use PL output spectra for spectroastrometric measurements, it is important to understand the wavelength-dependent behaviors of PL outputs and develop methods to calibrate the effects of time-varying wavefront errors in ground-based observations. We present experimental characterizations of the 3-port PL on the SCExAO testbed at the Subaru Telescope. We develop spectral response models of the PL and verify the behaviors with lab experiments. We find sinusoidal behavior of astrometric sensitivity of the 3-port PL as a function of wavelength, as expected from numerical simulations. Furthermore, we compare experimental and numerically simulated coupling maps and discuss their potential use for offsetting pointing errors. We then present a method of building PL spectral response models (solving for the transfer matrices as a function of wavelength) using coupling maps, which can be used for further calibration strategies., Comment: Accepted for publication in the Journal of Astronomical Telescopes, Instruments, and Systems (JATIS)

Liger is an adaptive optics (AO) fed imager and integral field spectrograph (IFS) designed to take advantage of the Keck All-sky Precision Adaptive-optics (KAPA) upgrade to the Keck I telescope. Liger adapts the design of the InfraRed Imaging Spectrograph (IRIS) for the Thirty Meter Telescope (TMT) to Keck by implementing a new imager and re-imaging optics. The performance of the imager is critical as it sequentially feeds the spectrograph and contains essential components such as the pupil wheel, filter wheel, and pupil viewing camera. We present the design and structural analysis of the Liger imager optical assembly including static, modal, and thermal simulations. We present the fabrication as well as the full assembly and characterization plan. The imager will be assembled bench-top in a clean room utilizing a coordinate-measuring machine (CMM) for warm alignment. To ensure optimal performance, the imager will be characterized in a test cryostat before integration with the full Liger instrument. This comprehensive approach to characterization ensures the precision and reliability of the imager, enhancing the observational capabilities of Liger and W.M. Keck Observatory.

We investigate the potential of photonic lantern (PL) fiber fed spectrometers for two-dimensional spectroastrometry. Spectroastrometry, a technique for studying small angular scales by measuring centroid shifts as a function of wavelength, is typically conducted using long-slit spectrographs. However, slit-based spectroastrometry requires observations with multiple position angles to measure two-dimensional spectroastrometric signals. In a typical configuration of PL-fed spectrometers, light from the focal plane is coupled into the few-moded PL, which is then split into several single-mode outputs, with the relative intensities containing astrometric information. The single-moded beams can be fed into a high-resolution spectrometer to measure wavelength-dependent centroid shifts. We perform numerical simulations of a standard 6-port PL and demonstrate its capability of measuring spectroastrometric signals. The effects of photon noise, wavefront errors, and chromaticity are investigated. When the PL is designed to have large linear responses to tip-tilts at the wavelengths of interest, the centroid shifts can be efficiently measured. Furthermore, we provide mock observations of detecting accreting protoplanets. PL spectroastrometry is potentially a simple and efficient technique for detecting spectroastrometric signals., Comment: Accepted for publication in the Journal of Astronomical Telescopes, Instruments, and Systems (JATIS)

Photonic lanterns are waveguide devices enabling high throughput single mode spectroscopy and high angular resolution. We aim to present the first on-sky demonstration of a photonic lantern (PL) operating in visible light, to measure its throughput and assess its potential for high-resolution spectroscopy of compact objects. We used the SCExAO instrument (a double stage extreme AO system installed at the Subaru telescope) and FIRST mid-resolution spectrograph (R 3000) to test the visible capabilities of the PL on internal source and on-sky observations. The best averaged coupling efficiency over the PL field of view was measured at 51% +/- 10% with a peak at 80%. We also investigate the relationship between coupling efficiency and the Strehl ratio for a PL, comparing them with those of a single-mode fiber (SMF). Findings show that in the AO regime, a PL offers better coupling efficiency performance than a SMF, especially in the presence of low spatial frequency aberrations. We observed Ikiiki (alpha Leo - mR = 1.37) and `Aua (alpha Ori - mR = -1.17) at a frame rate of 200 Hz. Under median seeing conditions (about 1 arcsec measured in H band) and large tip/tilt residuals (over 20 mas), we estimated an average light coupling efficiency of 14.5% +/- 7.4%, with a maximum of 42.8% at 680 nm. We were able to reconstruct both star's spectra, containing various absorption lines. The successful demonstration of this device opens new possibilities in terms of high throughput single-mode fiber-fed spectroscopy in the Visible. The demonstrated on-sky coupling efficiency performance would not have been achievable with a single SMF injection setup under similar conditions, partly because the residual tip/tilt alone exceeded the field of view of a visible SMF (18 mas at 700 nm). Thus emphasizing the enhanced resilience of PL technology to such atmospheric disturbances. The additional, Comment: Accepted in Astronomy and Astrophysics journal on 9/11/2024

GQ Lup B is one of the few substellar companions with a detected cicumplanetary disk, or CPD. Observations of the CPD suggest the presence of a cavity, possibly formed by an exo-satellite. Using the Keck Planet Imager and Characterizer (KPIC), a high contrast imaging suite that feeds a high resolution spectrograph (1.9-2.5 microns, R$\sim$35,000), we present the first dedicated radial velocity (RV) observations around a high-contrast, directly imaged substellar companion, GQ Lup B, to search for exo-satellites. Over 11 epochs, we find a best and median RV error of 400-1000 m/s, most likely limited by systematic fringing in the spectra due to transmissive optics within KPIC. With this RV precision, KPIC is sensitive to exomoons 0.6-2.8% the mass of GQ Lup B ($\sim 30 M_{\text{Jup}}$) at separations between the Roche limit and $65 R_{\text{Jup}}$, or the extent of the cavity inferred within the CPD detected around GQ Lup B. Using simulations of HISPEC, a high resolution infrared spectrograph planned to debut at W.M. Keck Observatory in 2026, we estimate future exomoon sensitivity to increase by over an order of magnitude, providing sensitivity to less massive satellites potentially formed within the CPD itself. Additionally, we run simulations to estimate the amount of material that different masses of satellites could clear in a CPD to create the observed cavity. We find satellite-to-planet mass ratios of $q > 2 \times 10^{-4}$ can create observable cavities and report a maximum cavity size of $\sim 51 \, R_{\text{Jup}}$ carved from a satellite., Comment: 15 pages, 5 figures

The Keck Planet Imager and Characterizer (KPIC) combines high contrast imaging with high resolution spectroscopy (R$\sim$35,000 in K band) to study directly imaged exoplanets and brown dwarfs in unprecedented detail. KPIC aims to spectrally characterize substellar companions through measurements of planetary radial velocities, spins, and atmospheric composition. Currently, the dominant source of systematic noise for KPIC is fringing, or oscillations in the spectrum as a function of wavelength. The fringing signal can dominate residuals by up to 10% of the continuum for high S/N exposures, preventing accurate wavelength calibration, retrieval of atmospheric parameters, and detection of planets with flux ratios less than 1% of the host star. To combat contamination from fringing, we first identify its three unique sources and adopt a physically informed model of Fabry-P\'{e}rot cavities to apply to post-processed data. We find this strategy can effectively model the fringing in observations of A0V/F0V stars, reducing the residual systematics caused by fringing by a factor of 2. Next, we wedge two of the transmissive optics internal to KPIC to eliminate two sources of fringing and confirm the third source as the entrance window to the spectrograph. Finally, we apply our previous model of the Fabry-P\'{e}rot cavity to new data taken with the wedged optics to reduce the amplitude of the residuals by a factor of 10., Comment: 13 pages, 6 figures

We present the joint astrometric and direct imaging discovery, mass measurement, and orbital analysis of HD 63754 B (HIP 38216 B), a companion near the stellar-substellar boundary orbiting ~20 AU from its Sun-like host. HD 63754 was observed in our ongoing high-contrast imaging survey targeting stars with significant proper-motion accelerations between Hipparcos and Gaia consistent with wide-separation substellar companions. We utilized archival HIRES and HARPS radial velocity (RV) data, together with the host star's astrometric acceleration extracted from the Hipparcos-Gaia Catalog of Accelerations (HGCA), to predict the location of the candidate companion around HD 63754 A. We subsequently imaged HD 63754 B at its predicted location using the Near Infrared Camera 2 (NIRC2) in the $L'$ band at the W. M. Keck Observatory. We then jointly modeled the orbit of HD 63754 B with RVs, Hipparcos-Gaia accelerations, and our new relative astrometry, measuring a dynamical mass of ${81.9}_{-5.8}^{+6.4} M_{jup}$, an eccentricity of ${0.260}_{-0.059}^{+0.065}$, and a nearly face-on inclination of $174.81_{-0.50}^{+0.48}$ degrees. For HD 63754 B, we obtain an L' band absolute magnitude of $L' = 11.39\pm0.06$ mag, from which we infer a bolometric luminosity of $log(L_{bol}/L_{\odot})= -4.55 \pm0.08$ dex using a comparison sample of L and T dwarfs with measured luminosities. Although uncertainties linger in age and dynamical mass estimates, our analysis points toward HD 63754 B's identity as a brown dwarf on the L/T transition rather than a low-mass star, indicated by its inferred bolometric luminosity and model-estimated effective temperature. Future RV, spectroscopic, and astrometric data such as those from JWST and Gaia DR4 will clarify HD 63754 B's mass, and enable spectral typing and atmospheric characterization.

Young, self-luminous super-Jovian companions discovered by direct imaging provide a challenging test of planet formation and evolution theories. By spectroscopically characterizing the atmospheric compositions of these super-Jupiters, we can constrain their formation histories. Here we present studies of the recently discovered HIP 99770 b, a 16 MJup high-contrast companion on a 17 au orbit, using the fiber-fed high-resolution spectrograph KPIC (R~35,000) on the Keck II telescope. Our K-band observations led to detections of H2O and CO in the atmosphere of HIP 99770 b. We carried out free retrieval analyses using petitRADTRANS to measure its chemical abundances, including the metallicity and C/O ratio, projected rotation velocity (vsini), and radial velocity (RV). We found that the companion's atmosphere has C/O=0.55(-0.04/+0.06) and [M/H]=0.26(-0.23/+0.24) (1{\sigma} confidence intervals), values consistent with those of the Sun and with a companion formation via gravitational instability or core accretion. The projected rotation velocity < 7.8 km/s is small relative to other directly imaged companions with similar masses and ages. This may imply a near pole-on orientation or effective magnetic braking by a circumplanetary disk. In addition, we added the companion-to-primary relative RV measurement to the orbital fitting and obtained updated constraints on orbital parameters. Detailed characterization of super-Jovian companions within 20 au like HIP 99770 b is critical for understanding the formation histories of this population., Comment: 18 pages, 5 figures, accepted to AJ

Brown dwarfs (BD) are model degenerate in age and mass. High-contrast imaging and spectroscopy of BD companions to host stars where the mass and age can be independently constrained by dynamics and stellar age indicators respectively provide valuable tests of BD evolution models. In this paper, we present a new epoch of Subaru/CHARIS H- and K-band observations of one such previously discovered system, HD 33632 Ab. We reanalyze the mass and orbit using our new epoch of extracted relative astrometry, and fit extracted spectra to the newest generation of equilibrium, disequilibrium, and cloudy spectral and evolution models for BDs. No spectral model perfectly agrees with evolutionary tracks and the derived mass and age, instead favoring a somewhat younger BD than the host star's inferred age. This tension can potentially be resolved using atmosphere and evolution models that consider both clouds and disequilibrium chemistry simultaneously, or by additional future spectra at higher resolution or in other band passes. Photometric measurements alone remain consistent with the luminosity predicted by evolutionary tracks. Our work highlights the importance of considering complexities like clouds, disequilibrium chemistry, and composition when comparing spectral models to evolutionary tracks., Comment: 15 pages, 7 figures, accepted for publication in The Astrophysical Journal

A Photonic Lantern (PL) is a novel device that efficiently converts a multi-mode fiber into several single-mode fibers. When coupled with an extreme adaptive optics (ExAO) system and a spectrograph, PLs enable high throughput spectroscopy at high angular resolution. The Subaru Coronagraphic Extreme Adaptive Optics (SCExAO) system of the Subaru Telescope recently acquired a PL that converts its multi-mode input into 19 single-mode outputs. The single mode outputs feed a R~4,000 spectrograph optimized for the 600 to 760 nm wavelength range. We present here the integration of the PL on SCExAO, and study the device performance in terms of throughput, field of view, and spectral reconstruction. We also present the first on-sky demonstration of a Visible PL coupled with an ExAO system, showing a significant improvement of x12 in throughput compared to the use of a sole single-mode fiber. This work paves the way towards future high throughput photonics instrumentation at small angular resolution., Comment: Proceeding published in SPIE Astronomical Telescopes + Instrumentation (2024) - Paper number 13096-25

The 16 Myr-old A0V star HD 156623 in the Scorpius--Centaurus association hosts a high-fractional-luminosity debris disk, recently resolved in scattered light for the first time by the Gemini Planet Imager (GPI) in polarized intensity. We present new analysis of the GPI H-band polarimetric detection of the HD 156623 debris disk, with particular interest in its unique morphology. This debris disk lacks a visible inner clearing, unlike the majority of low-inclination disks in the GPI sample and in Sco-Cen, and it is known to contain CO gas, positioning it as a candidate ``hybrid'' or ``shielded'' disk. We use radiative transfer models to constrain the geometric parameters of the disk based on scattered light data and thermal models to constrain the unresolved inner radius based on the system's spectral energy distribution (SED). We also compute a measurement of the polarized scattering phase function, adding to the existing sample of empirical phase function measurements. We find that HD 156623's debris disk inner radius is constrained to less than 26.6 AU from scattered light imagery and less than 13.4 AU from SED modeling at a 99.7% confidence interval, and suggest that gas drag may play a role in retaining sub-blowout size dust grains so close to the star., Comment: 18 pages, 10 figures. Accepted for publication in AJ

The Keck Planet Imager and Characterizer (KPIC), a series of upgrades to the Keck II Adaptive Optics System and Instrument Suite, aims to demonstrate high-resolution spectroscopy of faint exoplanets that are spatially resolved from their host stars. In this paper, we measure KPIC's sensitivity to companions as a function of separation (i.e., the contrast curve) using on-sky data collected over four years of operation. We show that KPIC is able to reach contrasts of $1.3 \times 10^{-4}$ at 90 mas and $9.2 \times 10^{-6}$ at 420 mas separation from the star, and that KPIC can reach planet-level sensitivities at angular separations within the inner working angle of coronagraphic instruments such as GPI and SPHERE. KPIC is also able to achieve more extreme contrasts than other medium-/high-resolution spectrographs that are not as optimized for high-contrast performance. We decompose the KPIC performance budget into individual noise terms and discuss limiting factors. The fringing that results from combining a high-contrast imaging system with a high-resolution spectrograph is identified as an important source of systematic noise. After mitigation and correction, KPIC is able to reach within a factor of 2 of the photon noise limit at separations < 200 mas. At large separations, KPIC is limited by the background noise performance of NIRSPEC., Comment: 16 pages, 6 figures, submitted to the proceedings of SPIE Astronomical Telescopes + Instrumentation 2024, 13096-69

Using Keck Planet Imager and Characterizer (KPIC) high-resolution ($R$~35000) spectroscopy from 2.29-2.49 $\mu$m, we present uniform atmospheric retrievals for eight young substellar companions with masses of ~10-30 $M_\textrm{Jup}$, orbital separations spanning ~50-360 au, and $T_\textrm{eff}$ between ~1500-2600 K. We find that all companions have solar C/O ratios, and metallicities, to within the 1-2$\sigma$ level, with the measurements clustered around solar composition. Stars in the same stellar associations as our systems have near-solar abundances, so these results indicate that this population of companions is consistent with formation via direct gravitational collapse. Alternatively, core accretion outside the CO snowline would be compatible with our measurements, though the high mass ratios of most systems would require rapid core assembly and gas accretion in massive disks. On a population level, our findings can be contrasted with abundance measurements for directly imaged planets with m<10 $M_\textrm{Jup}$, which show tentative atmospheric metal enrichment. In addition, the atmospheric compositions of our sample of companions are distinct from those of hot Jupiters, which most likely form via core accretion. For two companions with $T_\textrm{eff}$~1700-2000 K (kap And b and GSC 6214-210 b), our best-fit models prefer a non-gray cloud model with >3$\sigma$ significance. The cloudy models yield 2-3$\sigma$ lower $T_\textrm{eff}$ for these companions, though the C/O and [C/H] still agree between cloudy and clear models at the $1\sigma$ level. Finally, we constrain 12CO/13CO for three companions with the highest S/N data (GQ Lup b, HIP 79098 b, and DH Tau b), and report $v$sin($i$) and radial velocities for all companions., Comment: Accepted to ApJ, 52 pages, 20 figures

We used the Keck Planet Imager and Characterizer (KPIC) to obtain high-resolution (R$\sim$35,000) K-band spectra of kappa Andromedae b, a planetary-mass companion orbiting the B9V star, kappa Andromedae A. We characterized its spin, radial velocity, and bulk atmospheric parameters through use of a forward modeling framework to jointly fit planetary spectra and residual starlight speckles, obtaining likelihood-based posterior probabilities. We also detected H$_{2}$O and CO in its atmosphere via cross correlation. We measured a $v\sin(i)$ value for kappa And b of $38.42\pm{0.05}$ km/s, allowing us to extend our understanding of the population of close in bound companions at higher rotation rates. This rotation rate is one of the highest spins relative to breakup velocity measured to date, at close to $50\%$ of breakup velocity. We identify a radial velocity $-17.35_{-0.09}^{+0.05}$ km/s, which we use with existing astrometry and RV measurements to update the orbital fit. We also measure an effective temperature of $1700\pm{100}$ K and a $\log(g)$ of $4.7\pm{0.5}$ cgs dex., Comment: 14 pages, 7 figures, accepted for publication in AJ

We present the projected rotational velocity and molecular abundances for HD 33632 Ab obtained via Keck Planet Imager and Characterizer high-resolution spectroscopy. HD 33632 Ab is a nearby benchmark brown dwarf companion at a separation of $\sim$20 au that straddles the L/T transition. Using a forward-modeling framework with on-axis host star spectra, self-consistent substellar atmospheric and retrieval models for HD 33632 Ab, we derive a projected rotational velocity of 53 $\pm$ 3 km/s and carbon/water mass fractions of log CO = $-$2.3 $\pm$ 0.3 and log H$_2$O = $-$2.7 $\pm$ 0.2. The inferred carbon-to-oxygen ratio (C/O = 0.58 $\pm$ 0.14), molecular abundances, and metallicity ([C/H] = 0.0 $\pm$ 0.2 dex) of HD 33632 Ab are consistent with its host star. Although detectable methane opacities are expected in L/T transition objects, we did not recover methane in our KPIC spectra, partly due to the high $v\sin{i}$ and to disequilibrium chemistry at the pressures we are sensitive to. We parameterize the spin as the ratio of rotation over break-up velocity, and compare HD 33632 Ab to a compilation of >200 very low-mass objects (M$\lesssim$0.1 M$_{\odot}$) that have spin measurements in the literature. There appears to be no clear trend for the isolated field low-mass objects versus mass, but a tentative trend is identified for low-mass companions and directly imaged exoplanets, similar to previous findings. A larger sample of close-in gas giant exoplanets and brown dwarfs will critically examine our understanding of their formation and evolution through rotation and chemical abundance measurements., Comment: Accepted for publication in the Astrophysical Journal. 36 pages, 15 figures, 5 tables

The 1RXS J034231.8+121622 system consists of an M dwarf primary and a directly imaged low-mass stellar companion. We use high resolution spectroscopic data from Keck/KPIC to estimate the objects' atmospheric parameters and radial velocities (RVs). Using PHOENIX stellar models, we find that the primary has a temperature of 3460 $\pm$ 50 K a metallicity of 0.16 $\pm$ 0.04, while the secondary has a temperature of 2510 $\pm$ 50 K and a metallicity of $0.13\substack{+0.12 \\ -0.11}$. Recent work suggests this system is associated with the Hyades, placing it an older age than previous estimates. Both metallicities agree with current $[Fe/H]$ Hyades measurements (0.11 -- 0.21). Using stellar evolutionary models, we obtain significantly higher masses for the objects, of 0.30 $\pm$ 0.15 $M_\odot$ and 0.08 $\pm$ 0.01 $M_\odot$ (84 $\pm$ 11 $M_{Jup}$) respectively. Using the RVs and a new astrometry point from Keck/NIRC2, we find that the system is likely an edge-on, moderately eccentric ($0.41\substack{+0.27 \\ -0.08}$) configuration. We also estimate the C/O ratio of both objects using custom grid models, obtaining 0.42 $\pm$ 0.10 (primary) and 0.55 $\pm$ 0.10 (companion). From these results, we confirm that this system most likely went through a binary star formation process in the Hyades. The significant changes in this system's parameters since its discovery highlight the importance of high resolution spectroscopy for both orbital and atmospheric characterization of directly imaged companions., Comment: 30 pages, 18 figures, accepted for publication in AJ

Photonic lantern nulling (PLN) is a method for enabling the detection and characterization of close-in exoplanets by exploiting the symmetries of the ports of a mode-selective photonic lantern (MSPL) to cancel out starlight. A six-port MSPL provides four ports where on-axis starlight is suppressed, while off-axis planet light is coupled with efficiencies that vary as a function of the planet's spatial position. We characterize the properties of a six-port MSPL in the laboratory and perform the first testbed demonstration of the PLN in monochromatic light (1569 nm) and in broadband light (1450 nm to 1625 nm), each using two orthogonal polarizations. We compare the measured spatial throughput maps with those predicted by simulations using the lantern's modes. We find that the morphologies of the measured throughput maps are reproduced by the simulations, though the real lantern is lossy and has lower throughputs overall. The measured ratios of on-axis stellar leakage to peak off-axis throughput are around 10^(-2), likely limited by testbed wavefront errors. These null-depths are already sufficient for observing young gas giants at the diffraction limit using ground-based observatories. Future work includes using wavefront control to further improve the nulls, as well as testing and validating the PLN on-sky., Comment: 30 pages, 12 figures

Vortex fiber nulling (VFN) is a technique for detecting and characterizing faint companions at small separations from their host star. A near-infrared ($\sim2.3 \mu$m) VFN demonstrator mode was deployed on the Keck Planet Imager and Characterizer (KPIC) instrument at the Keck Observatory and presented earlier. In this paper, we present the first VFN companion detections. Three targets, HIP 21543 Ab, HIP 94666 Ab, and HIP 50319 B, were detected with host-companion flux ratios between 70 and 430 at and within one diffraction beamwidth ($\lambda/D$). We complement the spectra from KPIC VFN with flux ratio and position measurements from the CHARA Array to validate the VFN results and provide a more complete characterization of the targets. This paper reports the first direct detection of these three M dwarf companions, yielding their first spectra and flux ratios. Our observations provide measurements of bulk properties such as effective temperatures, radial velocities, and v$\sin{i}$, and verify the accuracy of the published orbits. These detections corroborate earlier predictions of the KPIC VFN performance, demonstrating that the instrument mode is ready for science observations., Comment: 13 pages, 2 figures; Accepted to ApJ Letters

Photonic Lanterns (PLs) are tapered waveguides that gradually transition from a multi-mode fiber geometry to a bundle of single-mode fibers (SMFs). They can efficiently couple multi-mode telescope light into a multi-mode fiber entrance at the focal plane and convert it into multiple single-mode beams. Thus, each SMF samples its unique mode (lantern principal mode) of the telescope light in the pupil, analogous to subapertures in aperture masking interferometry (AMI). Coherent imaging with PLs can be enabled by interfering SMF outputs and applying phase modulation, which can be achieved using a photonic chip beam combiner at the backend (e.g., the ABCD beam combiner). In this study, we investigate the potential of coherent imaging by interfering SMF outputs of a PL with a single telescope. We demonstrate that the visibilities that can be measured from a PL are mutual intensities incident on the pupil weighted by the cross-correlation of a pair of lantern modes. From numerically simulated lantern principal modes of a 6-port PL, we find that interferometric observables using a PL behave similarly to separated-aperture visibilities for simple models on small angular scales ($<\lambda/D$) but with greater sensitivity to symmetries and capability to break phase angle degeneracies. Furthermore, we present simulated observations with wavefront errors and compare them to AMI. Despite the redundancy caused by extended lantern principal modes, spatial filtering offers stability to wavefront errors. Our simulated observations suggest that PLs may offer significant benefits in the photon noise-limited regime and in resolving small angular scales at low contrast regime., Comment: Accepted for publication in ApJ

Spatially-resolved images of debris disks are necessary to determine disk morphological properties and the scattering phase function (SPF) which quantifies the brightness of scattered light as a function of phase angle. Current high-contrast imaging instruments have successfully resolved several dozens of debris disks around other stars, but few studies have investigated trends in the scattered-light, resolved population of debris disks in a uniform and consistent manner. We have combined Karhunen-Loeve Image Projection (KLIP) with radiative-transfer disk forward modeling in order to obtain the highest quality image reductions and constrain disk morphological properties of eight debris disks imaged by the Gemini Planet Imager at H-band with a consistent and uniformly-applied approach. In describing the scattering properties of our models, we assume a common SPF informed from solar system dust scattering measurements and apply it to all systems. We identify a diverse range of dust density properties among the sample, including critical radius, radial width, and vertical width. We also identify radially narrow and vertically extended disks that may have resulted from substellar companion perturbations, along with a tentative positive trend in disk eccentricity with relative disk width. We also find that using a common SPF can achieve reasonable model fits for disks that are axisymmetric and asymmetric when fitting models to each side of the disk independently, suggesting that scattering behavior from debris disks may be similar to Solar System dust., Comment: 23+5 pages, 12+6 figures, 15 pages of Online Supplemental Material included; Accepted for publication in MNRAS

The direct imaging of an Earth-like exoplanet will require sub-nanometric wavefront control across large light-collecting apertures, to reject host starlight and detect the faint planetary signal. Current adaptive optics (AO) systems, which use wavefront sensors that reimage the telescope pupil, face two challenges that prevent this level of control: non-common-path aberrations (NCPAs), caused by differences between the sensing and science arms of the instrument; and petaling modes: discontinuous phase aberrations caused by pupil fragmentation, especially relevant for the upcoming 30-m class telescopes. Such aberrations drastically impact the capabilities of high-contrast instruments. To address these issues, we can add a second-stage wavefront sensor to the science focal plane. One promising architecture uses the photonic lantern (PL): a waveguide that efficiently couples aberrated light into single-mode fibers (SMFs). In turn, SMF-confined light can be stably injected into high-resolution spectrographs, enabling direct exoplanet characterization and precision radial velocity measurements; simultaneously, the PL can be used for focal-plane wavefront sensing. We present a real-time experimental demonstration of the PL wavefront sensor on the Subaru/SCExAO testbed. Our system is stable out to around ~400 nm of low-order Zernike wavefront error, and can correct petaling modes. When injecting ~30 nm RMS of low order time-varying error, we achieve ~10x rejection at 1 s timescales; further refinements to the control law and lantern fabrication process should make sub-nanometric wavefront control possible. In the future, novel sensors like the PLWFS may prove to be critical in resolving the wavefront control challenges posed by exoplanet direct imaging., Comment: Accepted to ApJL

The unprecedented medium-resolution (R~1500-3500) near- and mid-infrared (1-18um) spectrum provided by JWST for the young (140+/-20Myr) low-mass (12-20MJup) L-T transition (L7) companion VHS1256b gives access to a catalogue of molecular absorptions. In this study, we present a comprehensive analysis of this dataset utilizing a forward modelling approach, applying our Bayesian framework, ForMoSA. We explore five distinct atmospheric models to assess their performance in estimating key atmospheric parameters: Teff, log(g), [M/H], C/O, gamma, fsed, and R. Our findings reveal that each parameter's estimate is significantly influenced by factors such as the wavelength range considered and the model chosen for the fit. This is attributed to systematic errors in the models and their challenges in accurately replicating the complex atmospheric structure of VHS1256b, notably the complexity of its clouds and dust distribution. To propagate the impact of these systematic uncertainties on our atmospheric property estimates, we introduce innovative fitting methodologies based on independent fits performed on different spectral windows. We finally derived a Teff consistent with the spectral type of the target, considering its young age, which is confirmed by our estimate of log(g). Despite the exceptional data quality, attaining robust estimates for chemical abundances [M/H] and C/O, often employed as indicators of formation history, remains challenging. Nevertheless, the pioneering case of JWST's data for VHS1256b has paved the way for future acquisitions of substellar spectra that will be systematically analyzed to directly compare the properties of these objects and correct the systematics in the models., Comment: 32 pages, 16 figures, 6 tables, 2 appendices

M dwarfs are common host stars to exoplanets but often lack atmospheric abundance measurements. Late-M dwarfs are also good analogs to the youngest substellar companions, which share similar $T_{\rm eff}\sim2300-2800~K$. We present atmospheric analyses for the M7.5 companion HIP 55507 B and its K6V primary star with Keck/KPIC high-resolution ($R\sim35,000$) $K$ band spectroscopy. First, by including KPIC relative radial velocities between the primary and secondary in the orbit fit, we improve the dynamical mass precision by 60% and find $M_B=88.0_{-3.2}^{+3.4}$ $M_{\rm Jup}$, putting HIP 55507 B above the stellar-substellar boundary. We also find that HIP 55507 B orbits its K6V primary star with $a=38^{+4}_{-3}$ AU and $e=0.40\pm0.04$. From atmospheric retrievals of HIP 55507 B, we measure $\rm [C/H]=0.24\pm0.13$, $\rm [O/H]=0.15\pm0.13$, and $\rm C/O=0.67\pm0.04$. Moreover, we strongly detect $\rm ^{13}CO$ ($7.8\sigma$ significance) and tentatively detect $\rm H_2^{18}O$ ($3.7\sigma$ significance) in companion's atmosphere, and measure $\rm ^{12}CO/^{13}CO=98_{-22}^{+28}$ and $\rm H_2^{16}O/H_2^{18}O=240_{-80}^{+145}$ after accounting for systematic errors. From a simplified retrieval analysis of HIP 55507 A, we measure $\rm ^{12}CO/^{13}CO=79_{-16}^{+21}$ and $\rm C^{16}O/C^{18}O=288_{-70}^{+125}$ for the primary star. These results demonstrate that HIP 55507 A and B have consistent $\rm ^{12} C/^{13}C$ and $\rm ^{16}O/^{18}O$ to the $<1\sigma$ level, as expected for a chemically homogeneous binary system. Given the similar flux ratios and separations between HIP 55507 AB and systems with young, substellar companions, our results open the door to systematically measuring $\rm ^{13}CO$ and $\rm H_2^{18}O$ abundances in the atmospheres of substellar or even planetary-mass companions with similar spectral types., Comment: Accepted to ApJ, 28 pages, 14 figures

Photonic lanterns (PLs) are tapered waveguides that gradually transition from a multi-mode fiber geometry to a bundle of single-mode fibers. In astronomical applications, PLs can efficiently couple multi-mode telescope light into a multi-mode fiber entrance and convert it into multiple single-mode beams. The output beams are highly stable and suitable for feeding into high-resolution spectrographs or photonic chip beam combiners. For instance, by using relative intensities in the output cores as a function of wavelength, PLs can enable spectroastrometry. In addition, by interfering beams in the output cores with a beam combiner in the backend, PLs can be used for high-throughput interferometric imaging. When used on an Extremely Large Telescope (ELT), with its increased sensitivity and angular resolution, the imaging and spectroastrometric capabilities of PLs will be extended to higher contrast and smaller angular scales. We study the potential spectroastrometry and imaging science cases of PLs on ELTs, including study of exomoons, broad-line regions of quasars, and inner circumstellar disks., Comment: AO4ELT7 conference proceedings 2023

We present high-resolution $K$-band emission spectra of the quintessential hot Jupiter HD 189733 b from the Keck Planet Imager and Characterizer (KPIC). Using a Bayesian retrieval framework, we fit the dayside pressure-temperature profile, orbital kinematics, mass-mixing ratios of H$_2$O, CO, CH$_4$, NH$_3$, HCN, and H$_2$S, and the $\rm ^{13}CO/^{12}CO$ ratio. We measure mass fractions of $\rm \log H_2O = -2.0^{+0.4}_{-0.4}$ and $\rm \log CO = -2.2^{+0.5}_{-0.5}$, and place upper limits on the remaining species. Notably, we find $\rm \log CH_4 < -4.5$ at 99\% confidence, despite its anticipated presence at the equilibrium temperature of HD 189733 b assuming local thermal equilibrium. We make a tentative ($\sim3\sigma$) detection of $\rm ^{13}CO$, and the retrieved posteriors suggest a $\rm ^{12}C/^{13}C$ ratio similar to or substantially less than the local interstellar value. The possible $\rm ^{13}C$ enrichment would be consistent with accretion of fractionated material in ices or in the protoplanetary disk midplane. The retrieved abundances correspond to a substantially sub-stellar atmospheric $\rm C/O = 0.3\pm0.1$, while the carbon and oxygen abundances are stellar to slightly super-stellar, consistent with core-accretion models which predict an inverse correlation between C/O and metallicity. The specific combination of low C/O and high metallicity suggests significant accretion of solid material may have occurred late in the formation process of HD 189733 b., Comment: 17 pages, 7 figures, 2 tables, accepted in AJ