Search

Vortex fiber nulling (VFN) is a technique for detecting and characterizing faint companions at small separations from their host star. A near-infrared (∼2.3 μ m) VFN demonstrator mode was deployed on the Keck Planet Imager and Characterizer (KPIC) instrument at the Keck Observatory and presented earlier. In this Letter, 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 ( λ / 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 Letter 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.

GQ Lup B is one of the few substellar companions with a detected cicumplanetary disk (CPD). Observations of the CPD suggest the presence of a cavity, possibly formed by an exosatellite. Using the Keck Planet Imager and Characterizer (KPIC), a high-contrast imaging suite that feeds a high-resolution spectrograph (1.9–2.5 µ m, R ∼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 exosatellites. Over 11 epochs, we find a best and median RV error of 400–1000 m s ^−1 , 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 (∼30 M _Jup ) at separations between the Roche limit and 65 R _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 × 10 ^−4 can create observable cavities and report a maximum cavity size of ∼51 R _Jup carved from a satellite.

We used the Keck Planet Imager and Characterizer to obtain high-resolution ( R ∼ 35,000) K -band spectra of κ Andromedae b, a planetary-mass companion orbiting the B9V star, κ 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 κ Andromedae b of 38.42 ± 0.05 km s ^−1 , 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 ^−1 , which we use with existing astrometry and radial velocity measurements to update the orbital fit. We also measure an effective temperature of 1700 ± 100 K and a $\mathrm{log}(g)$ of 4.7 ± 0.5 cgs dex.

Young, self-luminous super-Jovian companions discovered by direct imaging provide a challenging test for 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 M _Jup high-contrast companion on a 17 au orbit, using the fiber-fed high-resolution spectrograph KPIC ( ${ \mathcal R }$ ∼ 35,000) on the Keck II telescope. Our K -band observations led to detections of H _2 O 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 ( $v\sin i$ ), 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 σ confidence intervals), values consistent with those of the Sun and with a companion formation via gravitational instability or core accretion. The projected rotation velocity $v\sin (i)\lt 7.8$ km s ^−1 is small relative to other directly imaged companions with similar masses and ages. This may imply a nearly 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.

We present the projected rotational velocity and molecular abundances for HD 33632 Ab obtained via Keck Planet Imager and Characterizer (KPIC) high-resolution spectroscopy. HD 33632 Ab is a nearby benchmark brown dwarf companion at a separation of ∼20 au that straddles the L–T transition. Using a forward-modeling framework with on-axis host star spectra, which provides self-consistent substellar atmospheric and retrieval models for HD 33632 Ab, we derive a projected rotational velocity of 53 ± 3 km s ^−1 and carbon monoxide and water mass fractions of logCO = −2.3 ± 0.3 and logH _2 O = −2.7 ± 0.2, respectively. The inferred carbon-to-oxygen ratio (C/O = 0.58 ± 0.14), molecular abundances, and metallicity ([C/H] = 0.0 ± 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 to which we are sensitive. We parameterize the spin as the ratio of rotation to the breakup velocity, and compare HD 33632 Ab to a compilation of >200 very low-mass objects ( M ≲ 0.1 M _⊙ ) that have spin measurements in the literature. There appears to be no clear trend for the isolated low-mass field 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.

Using Keck Planet Imager and Characterizer high-resolution ( R ∼ 35,000) spectroscopy from 2.29 to 2.49 μ m, we present uniform atmospheric retrievals for eight young substellar companions with masses of ∼10–30 M _Jup , orbital separations spanning ∼50–360 au, and T _eff between ∼1500 and 2600 K. We find that all companions have solar C/O ratios and metallicities to within the 1 σ –2 σ 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 _Jup , which show tentative atmospheric metal enrichment compared to their host stars. 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 _eff ∼ 1700–2000 K ( κ And b and GSC 6214–210 b), our best-fit models prefer a nongray cloud model with >3 σ significance. The cloudy models yield 2 σ −3 σ lower T _eff for these companions, though the C/O and [C/H] still agree between cloudy and clear models at the 1 σ level. Finally, we constrain ^12 CO/ ^13 CO for three companions with the highest signal-to-noise ratio data (GQ Lup b, HIP 79098b, and DH Tau b) and report $v\sin i$ and radial velocities for all companions.

We present the first on-sky segmented primary mirror closed-loop piston control using a Zernike wavefront sensor (ZWFS) installed on the Keck II telescope. Segment cophasing errors are a primary contributor to contrast limits on Keck and will be necessary to correct for the next generation of space missions and ground-based extremely large telescopes, which will all have segmented primary mirrors. The goal of the ZWFS installed on Keck is to monitor and correct primary mirror cophasing errors in parallel with science observations. The ZWFS is ideal for measuring phase discontinuities such as segment cophasing errors and is one of the most sensitive WFSs, but has limited dynamic range. The vector-ZWFS at Keck works on the adaptive-optics-corrected wavefront and consists of a metasurface focal plane mask that imposes two different phase shifts on the core of the point-spread function to two orthogonal light polarizations, producing two pupil images. This design extends the dynamic range compared with the scalar ZWFS. The primary mirror segment pistons were controlled in closed loop using the ZWFS, improving the Strehl ratio on the NIRC2 science camera by up to 10 percentage points. We analyze the performance of the closed-loop tests, the impact on NIRC2 science data, and discuss the ZWFS measurements.

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 ± 50 K and a metallicity of 0.16 ± 0.04, while the secondary has a temperature of 2510 ± 50 K and a metallicity of ${0.13}_{-0.11}^{+0.12}$ . Recent work suggests this system is associated with the Hyades, giving it an older age than previous estimates. Both metallicities agree with current Hyades [Fe/H] measurements (0.11–0.21). Using stellar evolutionary models, we obtain significantly higher masses for the objects, 0.30 ± 0.15 M _⊙ and 0.08 ± 0.01 M _⊙ (84 ± 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}_{-0.08}^{+0.27}$ ) configuration. We also estimate the C/O ratio of both objects using custom grid models, obtaining 0.42 ± 0.10 (primary) and 0.55 ± 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.

Photonic lanterns (PLs) are tapered waveguides that gradually transition from a multimode fiber geometry to a bundle of single-mode fibers (SMFs). They can efficiently couple multimode telescope light into a multimode 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 the interference of 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 the interference of 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 (< λ / D ) but with greater sensitivity to symmetries and capability to break phase angle degeneracies. Furthermore, we present simulated observations with wave front errors (WFEs) and compare them to AMI. Despite the redundancy caused by extended lantern principal modes, spatial filtering offers stability to WFEs. Our simulated observations suggest that PLs may offer significant benefits in the photon-noise-limited regime and in resolving small angular scales at the low-contrast regime.

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 _eff ∼ 2300–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 ∼ 35,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}_{\mathrm{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}_{-3}^{+4}$ au and e = 0.40 ± 0.04. From atmospheric retrievals of HIP 55507 B, we measure [C/H] = 0.24 ± 0.13, [O/H] = 0.15 ± 0.13, and C/O = 0.67 ± 0.04. Moreover, we strongly detect ^13 CO (7.8 σ significance) and tentatively detect ${{\rm{H}}}_{2}^{18}{\rm{O}}$ (3.7 σ significance) in the companion’s atmosphere and measure ${}^{12}\mathrm{CO}{/}^{13}\mathrm{CO}={98}_{-22}^{+28}$ and ${{\rm{H}}}_{2}^{16}{\rm{O}}/{{\rm{H}}}_{2}^{18}{\rm{O}}={240}_{-80}^{+145}$ after accounting for systematic errors. From a simplified retrieval analysis of HIP 55507 A, we measure ${}^{12}\mathrm{CO}{/}^{13}\mathrm{CO}={79}_{-16}^{+21}$ and ${{\rm{C}}}^{16}{\rm{O}}/{{\rm{C}}}^{18}{\rm{O}}={288}_{-70}^{+125}$ for the primary star. These results demonstrate that HIP 55507 A and B have consistent ^12 C/ ^13 C and ^16 O/ ^18 O to the

We present high-resolution K -band emission spectra of the quintessential hot Jupiter HD 189733 b from the Keck Planet Imager and Characterizer. 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 ^13 CO/ ^12 CO ratio. We measure mass fractions of ${\mathrm{logH}}_{2}{\rm{O}}=-{2.0}_{-0.4}^{+0.4}$ and $\mathrm{logCO}=-{2.2}_{-0.5}^{+0.5}$ , and place upper limits on the remaining species. Notably, we find logCH _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 (∼3 σ ) detection of ^13 CO, and the retrieved posteriors suggest a ^12 C/ ^13 C ratio similar to or substantially less than the local interstellar value. The possible ^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 substellar atmospheric C/O = 0.3 ± 0.1, while the carbon and oxygen abundances are stellar to slightly superstellar, 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.

The direct imaging of an Earth-like exoplanet will require sub-nanometric wave-front control across large light-collecting apertures to reject host starlight and detect the faint planetary signal. Current adaptive optics systems, which use wave-front sensors that reimage the telescope pupil, face two challenges that prevent this level of control: non-common-path aberrations, 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 wave-front 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 wave-front sensing. We present a real-time experimental demonstration of the PL wave-front sensor on the Subaru/SCExAO testbed. Our system is stable out to around ±400 nm of low-order Zernike wave-front error and can correct petaling modes. When injecting ∼30 nm rms of low-order time-varying error, we achieve ∼10× rejection at 1 s timescales; further refinements to the control law and lantern fabrication process should make sub-nanometric wave-front control possible. In the future, novel sensors like the PL wave-front sensor may prove to be critical in resolving the wave-front control challenges posed by exoplanet direct imaging.

Nulling interferometry is a technique combining lights from different telescopes or apertures to observe weak sources nearby bright ones. The authors report the first nulling interferometer implemented in a photonic chip doing spectrally dispersed nulling on several baselines, simultaneously.

Photonic technologies offer numerous functionalities that can be used to realize astrophotonic instruments. The most spectacular example to date is the ESO Gravity instrument at the Very Large Telescope in Chile that combines the light-gathering power of four 8 m telescopes through a complex photonic interferometer. Fully integrated astrophotonic devices stand to offer critical advantages for instrument development, including extreme miniaturization when operating at the diffraction-limit, as well as integration, superior thermal and mechanical stabilization owing to the small footprint, and high replicability offering significant cost savings. Numerous astrophotonic technologies have been developed to address shortcomings of conventional instruments to date, including for example the development of photonic lanterns to convert from multimode inputs to single mode outputs, complex aperiodic fiber Bragg gratings to filter OH emission from the atmosphere, complex beam combiners to enable long baseline interferometry with for example, ESO Gravity, and laser frequency combs for high precision spectral calibration of spectrometers. Despite these successes, the facility implementation of photonic solutions in astronomical instrumentation is currently limited because of (1) low throughputs from coupling to fibers, coupling fibers to chips, propagation and bend losses, device losses, etc, (2) difficulties with scaling to large channel count devices needed for large bandwidths and high resolutions, and (3) efficient integration of photonics with detectors, to name a few. In this roadmap, we identify 24 key areas that need further development. We outline the challenges and advances needed across those areas covering design tools, simulation capabilities, fabrication processes, the need for entirely new components, integration and hybridization and the characterization of devices. To realize these advances the astrophotonics community will have to work cooperatively with industrial partners who have more advanced manufacturing capabilities. With the advances described herein, multi-functional integrated instruments will be realized leading to novel observing capabilities for both ground and space based platforms, enabling new scientific studies and discoveries.

Wolf-Rayet (WR) 140 is the archetypal periodic dust-forming colliding-wind binary that hosts a carbon-rich WR (WC) star and an O-star companion with an orbital period of 7.93 yr and an orbital eccentricity of 0.9. Throughout the past few decades, multiple dust-formation episodes from WR 140 have been observed that are linked to the binary orbit and occur near the time of periastron passage. Given its predictable dust-formation episodes, WR 140 presents an ideal astrophysical laboratory to investigate the formation and evolution of dust in the hostile environment around a massive binary system. In this paper, we present near- and mid-infrared (IR) spectroscopic and imaging observations of WR 140 with Subaru/SCExAO+CHARIS, Keck/NIRC2+PyWFS, and Subaru/Cooled Mid-Infrared Camera and Spectrograph taken between 2020 June and September that resolve the circumstellar dust emission linked to its most recent dust-formation episode in 2016 December. Our spectral energy distribution analysis of WR 140's resolved circumstellar dust emission reveals the presence of a hot ( T _d ∼ 1000 K) near-IR dust component that is co-spatial with the previously known and cooler ( T _d ∼ 500 K) mid-IR dust component composed of 300–500 Å sized dust grains. We attribute the hot near-IR dust emission to the presence of nano-sized (nanodust) grains and suggest they were formed from grain–grain collisions or the rotational disruption of the larger grain size population by radiative torques in the strong radiation field from the central binary. Lastly, we speculate on the astrophysical implications of nanodust formation around colliding-wind WC binaries, which may present an early source of carbonaceous nanodust in the interstellar medium.

We present Keck Planet Imager and Characterizer (KPIC) high-resolution ( R ∼35,000) K -band thermal emission spectroscopy of the ultrahot Jupiter WASP-33b. The use of KPIC’s single-mode fibers greatly improves both blaze and line-spread stabilities relative to slit spectrographs, enhancing the cross-correlation detection strength. We retrieve the dayside emission spectrum with a nested-sampling pipeline, which fits for orbital parameters, the atmospheric pressure–temperature profile, and the molecular abundances. We strongly detect the thermally inverted dayside and measure mass-mixing ratios for CO ( ${\mathrm{logCO}}_{\mathrm{MMR}}=-{1.1}_{-0.6}^{+0.4}$ ), H _2 O ( ${\mathrm{logH}}_{2}{{\rm{O}}}_{\mathrm{MMR}}\,=-{4.1}_{-0.9}^{+0.7}$ ), and OH ( ${\mathrm{logOH}}_{\mathrm{MMR}}=-{2.1}_{-1.1}^{+0.5}$ ), suggesting near-complete dayside photodissociation of H _2 O. The retrieved abundances suggest a carbon- and possibly metal-enriched atmosphere, with a gas-phase C/O ratio of ${0.8}_{-0.2}^{+0.1}$ , consistent with the accretion of high-metallicity gas near the CO _2 snow line and post-disk migration or with accretion between the soot and H _2 O snow lines. We also find tentative evidence for ^12 CO/ ^13 CO ∼ 50, consistent with values expected in protoplanetary disks, as well as tentative evidence for a metal-enriched atmosphere (2–15 × solar). These observations demonstrate KPIC’s ability to characterize close-in planets and the utility of KPIC’s improved instrumental stability for cross-correlation techniques.

Benchmark brown dwarf companions with well-determined ages and model-independent masses are powerful tools to test substellar evolutionary models and probe the formation of giant planets and brown dwarfs. Here, we report the independent discovery of HIP 21152 B, the first imaged brown dwarf companion in the Hyades, and conduct a comprehensive orbital and atmospheric characterization of the system. HIP 21152 was targeted in an ongoing high-contrast imaging campaign of stars exhibiting proper-motion changes between Hipparcos and Gaia, and was also recently identified by Bonavita et al. (2022) and Kuzuhara et al. (2022). Our Keck/NIRC2 and SCExAO/CHARIS imaging of HIP 21152 revealed a comoving companion at a separation of 0.″37 (16 au). We perform a joint orbit fit of all available relative astrometry and radial velocities together with the Hipparcos-Gaia proper motions, yielding a dynamical mass of ${24}_{-4}^{+6}\,{M}_{\mathrm{Jup}}$ , which is 1–2 σ lower than evolutionary model predictions. Hybrid grids that include the evolution of cloud properties best reproduce the dynamical mass. We also identify a comoving wide-separation (1837″ or 7.9 × 10 ^4 au) early-L dwarf with an inferred mass near the hydrogen-burning limit. Finally, we analyze the spectra and photometry of HIP 21152 B using the Saumon & Marley (2008) atmospheric models and a suite of retrievals. The best-fit grid-based models have f _sed = 2, indicating the presence of clouds, T _eff = 1400 K, and $\mathrm{log}g=4.5\,\mathrm{dex}$ . These results are consistent with the object’s spectral type of T0 ± 1. As the first benchmark brown dwarf companion in the Hyades, HIP 21152 B joins the small but growing number of substellar companions with well-determined ages and dynamical masses.

The detection of satellites around extrasolar planets, so called exomoons, remains a largely unexplored territory. In this work, we study the potential of detecting these elusive objects from radial velocity monitoring of self-luminous, directly imaged planets. This technique is now possible thanks to the development of dedicated instruments combining the power of high-resolution spectroscopy and high-contrast imaging. First, we demonstrate a sensitivity to satellites with a mass ratio of 1%–4% at separations similar to the Galilean moons from observations of a brown-dwarf companion (HR 7672 B; K _mag = 13; 0.″7 separation) with the Keck Planet Imager and Characterizer ( R ∼ 35,000 in the K band) at the W. M. Keck Observatory. Current instrumentation is therefore already sensitive to large unresolved satellites that could be forming from gravitational instability akin to binary star formation. Using end-to-end simulations, we then estimate that future instruments such as the Multi-Object Diffraction-limited High-resolution Infrared Spectrograph, planned for the Thirty Meter Telescope, should be sensitive to satellites with mass ratios of ∼10 ^−4 . Such small moons would likely form in a circumplanetary disk similar to the Jovian satellites in the solar system. Looking for the Rossiter–McLaughlin effect could also be an interesting pathway to detecting the smallest moons on short orbital periods. Future exomoon discoveries will allow precise mass measurements of the substellar companions that they orbit and provide key insight into the formation of exoplanets. They would also help constrain the population of habitable Earth-sized moons orbiting gas giants in the habitable zone of their stars.

The formation and evolution pathway for the directly imaged multiplanetary system HR 8799 remains mysterious. Accurate constraints on the chemical composition of the planetary atmosphere(s) are key to solving the mystery. We perform a detailed atmospheric retrieval on HR 8799 c to infer the chemical abundances and abundance ratios using a combination of photometric data along with low- and high-resolution spectroscopic data ( R ∼ 20–35,000). We specifically retrieve [C/H], [O/H], and C/O and find them to be ${0.55}_{-0.39}^{+0.36}$ , ${0.47}_{-0.32}^{+0.31}$ , and ${0.67}_{-0.15}^{+0.12}$ at 68% confidence. The superstellar C and O abundances, yet a stellar C/O ratio, reveal a potential formation pathway for HR 8799 c. Planet c, and likely the other gas giant planets in the system, formed early on (likely within ∼1 Myr), followed by further atmospheric enrichment in C and O through the accretion of solids beyond the CO ice line. The enrichment either preceded or took place during the early phase of the inward migration to the current planet locations.