---
title: 'HD 86226 c: Short-Period Sub-Neptune'
url: https://www.emergentmind.com/topics/hd-86226-c
type: topic
---

# HD 86226 c: Short-Period Sub-Neptune

Searching arXiv for the specified papers to ground the article in current records.
HD 86226 c is a short-period sub-Neptune transiting the bright, solar-type star HD 86226. It was identified in TESS photometry and characterized through joint transit and radial-velocity analyses, which established a planet with radius \(2.16 \pm 0.08\,R_\oplus\), mass \(7.25^{+1.19}_{-1.12}\,M_\oplus\), orbital period \(3.98442 \pm 0.00018\ {\rm d}\), and equilibrium temperature \(1311 \pm 28\ {\rm K}\) assuming zero Bond albedo [2007.13927]. Subsequent HST transmission spectroscopy found that the planet’s near-infrared spectrum is essentially featureless, with a constant transit depth of \(418 \pm 14\ {\rm ppm}\) and a spectral modulation of only \(0.01^{+0.17}_{-0.01}\) scale heights for an H/He-dominated atmosphere, excluding a cloud-free solar-metallicity atmosphere at \(6.5\,\sigma\) confidence [2507.13439]. Within the broader architecture of the HD 86226 system, the planet is also notable as the inner transiting companion to a long-period giant planet, HD 86226 b, making it relevant to studies of inner light planets in systems with outer giant planets [2509.26232].

## 1. Discovery, confirmation, and system context

HD 86226 c was detected by the Transiting Exoplanet Survey Satellite mission around the bright (\(V=7.9\)) star HD 86226, also designated TOI-652 and TIC 22221375. TESS observed the star in short-cadence mode in Sector 9, and the SPOC pipeline flagged a periodic transit-like signal at about 3.98 days. An independent detrending of the TESS light curve followed by a box-fitting least squares periodogram recovered the same signal, with transit depth \(\sim 360\) ppm, signal-to-noise \(\sim 20\), and period near 3.98 days [2007.13927].

The validation sequence was extensive. The reported checks included TESS centroid tests, odd-even transit depth checks, ghost diagnostics, and bootstrap validation from the SPOC Data Validation Report. Ground-based follow-up photometry with LCOGT/SINISTRO ruled out neighboring eclipsing binaries, and speckle imaging with SOAR found no nearby companions within \(3''\) and constrained other companions to \(\Delta I \sim 6\) beyond \(0.5''\). Archival and new radial velocities confirmed that the signal is associated with the primary star, and TOI-652.01 was therefore identified as the planet HD 86226 c [2007.13927].

The system had already been known to host a long-period giant planet, HD 86226 b. Earlier Doppler work had reported the outer companion as a very eccentric giant, but later and more comprehensive analyses revised that picture substantially. The TESS-era joint fit yielded for HD 86226 b a period of \(1628^{+22}_{-21}\ {\rm d}\), RV semi-amplitude \(7.74^{+0.69}_{-0.70}\ {\rm m\,s^{-1}}\), eccentricity \(0.059^{+0.062}_{-0.039}\), minimum mass \(0.45^{+0.04}_{-0.05}\,M_{\rm Jup}\), and semi-major axis \(2.73 \pm 0.06\ {\rm au}\) [2007.13927]. A later homogeneous RV re-analysis using CORALIE, HARPS, and PFS treated HD 86226 c specifically as the inner planet in a system with a known outer giant and confirmed that the short-period planet is robustly present [2509.26232].

This architecture places HD 86226 c in a category of systems used to investigate whether outer giant planets correlate with or suppress the presence of inner light planets. The later survey paper identified HD 86226 as one of three systems in a 26-star Giant sample with an inner planet detection, while emphasizing that the sample size remained too small for a definitive correlation statement [2509.26232]. A plausible implication is that HD 86226 c is important not only as an atmospheric target but also as a benchmark for formation and system-architecture studies.

## 2. Stellar and planetary properties

The host star is consistently described as a near-solar G-type star. In the atmospheric characterization study, the adopted stellar properties were \(T_{\rm eff} = 5980 \pm 70\) K, \(\log g = 4.55^{+0.40}_{-0.48}\), \([\mathrm{Fe/H}] = -0.07^{+0.15}_{-0.14}\), and \(R_\star = 1.05 \pm 0.02\,R_\odot\), and the star was noted to be relatively quiet, consistent with a 25-day rotation period and no strong photometric variability in the ground-based monitoring [2507.13439]. Other analyses reported closely similar stellar parameters. The TESS discovery paper derived \(M_* = 1.019^{+0.061}_{-0.066}\,M_\odot\), \(R_* = 1.053^{+0.026}_{-0.026}\,R_\odot\), \(T_{\rm eff} = 5863^{+88}_{-88}\,{\rm K}\), \([{\rm Fe/H}] = 0.018^{+0.057}_{-0.043}\), age \(= 4.6^{+3.7}_{-2.7}\) Gyr, and distance \(d = 45.57 \pm 0.12\) pc [2007.13927]. The homogeneous RV study gave spectral type G2 V, distance \(45.68^{+0.09}_{-0.08}\) pc, effective temperature \(6007^{+75}_{-71}\) K, metallicity \([\mathrm{Fe/H}] = -0.015 \pm 0.041\), radius \(1.03^{+0.02}_{-0.02}\,R_\odot\), mass \(1.08^{+0.05}_{-0.05}\,M_\odot\), chromospheric activity \(\log R'_{\rm HK} = -4.946 \pm 0.012\), and rotation period \(22.8 \pm 3.0\) d [2509.26232].

For HD 86226 c itself, the joint transit and RV fit reported orbital period \(P = 3.98442 \pm 0.00018\ {\rm d}\), mid-transit time \(T_0 = 8543.2539 \pm 0.0007\) \(({\rm BJD}_{\rm TDB} - 2450000)\), scaled semi-major axis \(a/R_* = 10.11^{+0.07}_{-0.08}\), impact parameter \(b = 0.63^{+0.06}_{-0.08}\), inclination \(i_p = 86.45^{+0.26}_{-0.16}\ {\rm deg}\), RV semi-amplitude \(K = 2.89^{+0.46}_{-0.43}\ {\rm m\,s^{-1}}\), eccentricity \(e = 0.075^{+0.065}_{-0.048}\), argument of periastron \(\omega = 196^{+60}_{-90}\ {\rm deg}\), radius \(R_p = 2.16 \pm 0.08\,R_\oplus\), mass \(M_p = 7.25^{+1.19}_{-1.12}\,M_\oplus\), semi-major axis \(a = 0.049 \pm 0.001\ {\rm au}\), density \(\rho_p = 3.97^{+0.78}_{-0.73}\ {\rm g\,cm^{-3}}\), equilibrium temperature \(T_{\rm eq} = 1311 \pm 28\ {\rm K}\), and transit duration \(T_D = 3.12^{+0.12}_{-0.05}\) hr [2007.13927].

The HST atmospheric study adopted a broadband-fit planet radius of \(R_p = 2.313 \pm 0.051\,R_\oplus\), with the mass taken as \(7.25\,M_\oplus\) with \(\sim 1.19\,M_\oplus\) uncertainty [2507.13439]. This suggests that modest differences in reported radius reflect differences in analysis context rather than a disagreement over the existence or basic character of the planet.

The planet’s location near both the radius gap and the hot Neptune desert was already emphasized at discovery [2007.13927]. Because the orbital period is about 4 days and the equilibrium temperature is about \(1310\)–\(1311\) K, HD 86226 c occupies a strongly irradiated regime in which envelope retention, atmospheric loss, and compositional diversity are all astrophysically consequential.

## 3. Observational basis for atmospheric characterization

The principal atmospheric dataset came from the Sub-neptune Planetary Atmosphere Characterization Experiment (SPACE) Program. The observations combined HST/WFC3 G141 transmission spectroscopy over \(1.1\)–\(1.7 \,\mu\)m with HST/STIS ultraviolet characterization of the host star, including reconstruction of the stellar UV/EUV spectrum and Ly\(\alpha\) [2507.13439].

Nine transits were observed with WFC3. Because the target is bright and the scan speed was high, the spectrum shifted on the detector by up to \(\sim 2\) pixels. The analysis therefore had to model a row-position-dependent systematic. Two independent reductions were performed, using PACMAN and Eureka!, and both recovered a flat spectrum. PACMAN was treated as the more reliable result because it explicitly handled the detector-position systematics better [2507.13439].

The systematics modeling is an important part of the interpretation because the reported atmospheric result is not the detection of a spectral feature but the robust recovery of a nearly constant transit depth in the presence of instrumental structure. PACMAN modeled orbit-long systematics with
\[
\exp(-r_1 \times t_\mathrm{orbit} - r_2),
\]
where \(t_\mathrm{orbit}\) is time since orbit start. Eureka!’s ramp model used
\[
h_0\times \exp(-h_1 \times t_\mathrm{batch}), \quad t_\mathrm{batch} = (t_\mathrm{local}-h_5)\%h_4 .
\]
These expressions are instrumental rather than atmospheric, but they delimit the reduction framework under which the transmission spectrum was extracted [2507.13439].

The ultraviolet stellar characterization also matters physically. The paper notes that the star is relatively quiet and that the STIS data were used to reconstruct the stellar UV/EUV spectrum. Within the paper’s logic, this constrains the host-star environment in which the planet’s atmosphere is irradiated. The combination of near-infrared transmission spectroscopy and ultraviolet stellar characterization is therefore not merely observationally convenient; it is part of the broader attempt to place the atmosphere in the context of host-star forcing [2507.13439].

## 4. Featureless transmission spectrum and atmospheric retrievals

The main empirical atmospheric result is that HD 86226 c has a featureless near-infrared transmission spectrum. For the PACMAN spectrum, a constant-depth fit gives
\[
D = 418 \pm 14\ \mathrm{ppm},
\]
and the data are consistent with a constant transit depth at only \(0.4\,\sigma\) from flat. The broadband WFC3 depth is
\[
410 \pm 11\ \mathrm{ppm},
\]
and the spectroscopic bins mostly scatter around \(\sim 400\)–450 ppm with no compelling molecular structure [2507.13439].

The spectral modulation is reported as
\[
0.01^{+0.17}_{-0.01}\ \text{scale heights}
\]
for an H/He-dominated atmosphere. In the paper’s interpretation, this is extraordinarily small and is far flatter than expected for a clear hydrogen-rich atmosphere [2507.13439]. That conclusion is strengthened by retrieval calculations using petitRADTRANS, which show that a cloud-free solar-metallicity H/He atmosphere is ruled out at
\[
6.5\,\sigma .
\]

The paper notes an important caveat. If the planet mass is allowed to float freely, the exclusion of the cloud-free solar-metallicity H/He case weakens to \(2.9\sigma\), but that solution requires
\[
M_p \approx 11.1 \pm 0.8\,M_\oplus ,
\]
which conflicts with the dynamical mass estimate [2507.13439]. The robust conclusion is therefore not that every H/He-rich interpretation is impossible, but that the specific case of a clear, solar-metallicity H/He atmosphere is not viable under the observed mass constraint.

The retrieval setup assumed an isothermal atmosphere at \(T_\mathrm{eq}=1310\) K, hydrostatic equilibrium, a reference radius fixed by the broadband transit depth, a free reference pressure \(P_\mathrm{ref}\), and optional gray clouds at \(P_\mathrm{cloud}\). The gray-cloud retrieval has a near one-to-one degeneracy between the cloud-top pressure \(P_\mathrm{cloud}\) and the reference pressure \(P_\mathrm{ref}\), since both shift the effective transit radius [2507.13439]. This clarifies why the data strongly disfavor some classes of models while leaving a residual degeneracy between high mean molecular weight and cloud opacity.

## 5. Atmospheric interpretations: metal enrichment and refractory clouds

The atmospheric analysis explored two principal explanatory classes. The first is a cloud-free but metal-rich atmosphere. Using scaled-solar equilibrium chemistry and no gray cloud deck, the flat spectrum can be matched if the atmosphere has
\[
[\mathrm{M}/\mathrm{H}] = 2.9^{+0.1}_{-0.2},
\]
with a \(3\sigma\) lower limit
\[
[\mathrm{M}/\mathrm{H}] > 2.3 .
\]
This corresponds to a minimum mean molecular weight of about \(6\) u [2507.13439].

The paper defines metallicity as
\[
\left[\frac{\mathrm{M}}{\mathrm{H}}\right] = \log_{10}\left(\frac{N_\mathrm{M}}{N_\mathrm{H}}\right)_\mathrm{planet} - \log_{10}\left(\frac{N_\mathrm{M}}{N_\mathrm{H}}\right)_\odot ,
\]
where \(N_\mathrm{M}\) is the abundance of elements heavier than He and \(N_\mathrm{H}\) is hydrogen [2507.13439]. In this framework, the flatness of the spectrum is explained by the reduction of atmospheric scale height through a substantially increased mean molecular weight.

The second class of explanation is cloud opacity. The featureless spectrum can be reproduced by clouds of MnS, MgSiO\(_3\), or Fe. These species were selected because their condensation curves intersect the relevant \(T\)-\(P\) regime of HD 86226 c [2507.13439]. The analysis used radiative-convective temperature-pressure profiles and compared them with vapor-pressure curves for candidate condensates. The condensation levels of MnS, MgSiO\(_3\), and Fe cross the observable atmospheric region, roughly
\[
10^{-2} \text{ to } 10^{-4}\ \mathrm{bar},
\]
while Na\(_2\)S, KCl, and ZnS do not [2507.13439].

This distinction is central because HD 86226 c is hot enough that methane-based haze formation is suppressed. The paper explicitly states that the equilibrium temperature of about \(1310\) K suppresses CH\(_4\)-based haze formation and favors CO as the dominant carbon reservoir [2507.13439]. One might therefore have expected a clearer atmosphere than in cooler sub-Neptunes, yet the observed spectrum remains nearly flat. The authors conclude that the planet likely falls into one of two categories: a very metal-rich atmosphere with large mean molecular weight, or an atmosphere obscured by high-altitude refractory clouds more typical of hot-Jupiter-like condensation physics than of cooler sub-Neptunes [2507.13439].

A common misconception would be to interpret the flat spectrum as proof that the atmosphere is absent. The papers do not support that conclusion. The discovery paper states that the density is low enough to imply a small volatile envelope [2007.13927], while the atmospheric paper shows that the spectrum is compatible with either strong metal enrichment or specific cloud species rather than with an atmosphere-free rocky body [2507.13439].

## 6. Interior structure, population context, and future tests

From mass and radius alone, HD 86226 c is not consistent with a purely rocky interpretation. The discovery paper states that the planet lies slightly above the silicate-only mass-radius curve and interprets the density as low enough to imply at least a small volatile envelope [2007.13927]. The interior model adopted a layered structure with iron core, silicate mantle, water layer, and H/He envelope, with the structure equations solved using a generalized Bayesian inference approach. A notable modeling detail is that the transit radius was defined by the chord optical depth condition
\[
\tau_{\rm ch} = \frac{2}{3}.
\]

The inferred interior composition is highly degenerate. Under priors on stellar Fe/Si and Mg/Si, the Bayesian interior modeling gave core mass fraction
\[
0.35^{+0.22}_{-0.16},
\]
mantle mass fraction
\[
0.33^{+0.24}_{-0.20},
\]
water mass fraction
\[
0.32^{+0.21}_{-0.17},
\]
and atmospheric mass fraction
\[
6.2\times 10^{-5}\,M_{\rm total},
\]
corresponding to an H/He mass of about \(4.6\times 10^{-4}\,M_\oplus\) and a thickness of about \(0.39\,R_\oplus\) [2007.13927]. These values show that the presence of volatiles is favored, while the relative contributions of rock, iron, and water are not tightly determined.

In population context, HD 86226 c is especially noteworthy because the atmospheric study concludes that it does not follow the aerosol trend of sub-Neptunes found by previous studies [2507.13439]. At \(T_\mathrm{eq}\approx 1310\) K, a clearer atmosphere might have been expected than on cooler sub-Neptunes, but the observed transmission spectrum remains nearly featureless. This suggests that temperature alone is not a sufficient predictor of observable spectral amplitude in the sub-Neptune regime.

The system is also used in the literature on planetary architecture. The homogeneous RV survey describes HD 86226 c as one of only three inner planet detections in a 26-star Giant sample, with all three detected inner planets having periods below 6 days [2509.26232]. The authors state that their relatively low number of detections seems to contradict previous studies that found a strong outer giant planet–inner light planet correlation, although they do not yet conclude a definitive correlation or anticorrelation because the sample is still small and completeness corrections are not yet fully done. HD 86226 c therefore functions as an informative but non-dispositive data point in that debate.

The atmospheric paper identifies specific future observations that could break the present degeneracy. JWST/NIRSpec G395H can look for CO\(_2\) above \(\sim 3\,\mu\)m, and a strong CO\(_2\) feature would support the metal-rich scenario, whereas a weak or absent feature would point more toward cloud opacity or an even more unusual composition. MIRI could test cloud composition because MgSiO\(_3\) has a prominent feature around \(8\)–\(12 \,\mu\)m, while Fe and MnS are less distinctive in that region [2507.13439]. The paper also notes that optical or visible observations and phase curves could probe Rayleigh slopes, albedo, and cloud coverage.

Taken together, the available measurements define HD 86226 c as a hot, short-period sub-Neptune around a quiet G-type star, with a flat WFC3 transmission spectrum, a constant transit depth of \(418\pm14\) ppm, a spectral amplitude of only \(\sim 0.01\) scale heights, and a \(6.5\sigma\) exclusion of a cloud-free solar-metallicity H/He atmosphere [2507.13439]. Its significance lies in the conjunction of precise bulk properties, anomalously muted transmission features, and a system architecture that links inner sub-Neptune evolution to the presence of an outer giant companion.

Source: https://www.emergentmind.com/topics/hd-86226-c