Papers
Topics
Authors
Recent
Search
2000 character limit reached

HD 143811(AB)b: Directly Imaged Circumbinary Companion

Updated 10 July 2026
  • HD 143811(AB)b is a directly imaged circumbinary substellar companion orbiting a young spectroscopic binary in the Scorpius–Centaurus association, characterized by a wide (~64 au) orbit and a mass of ~5.6–6.1 M_Jup.
  • Detection involved multi-epoch high-contrast imaging with Gemini Planet Imager, Keck/NIRC2, and SPHERE, confirming common proper motion and refining orbital parameters (a ≈ 64 au, e ≈ 0.23, P ≈ 326 yr).
  • Atmospheric modeling indicates an effective temperature of ~1040 K with a red, dusty spectrum, offering key insights into planet formation and dynamical evolution in close binary systems.

HD 143811(AB)b, also written HD 143811 AB b, is a directly imaged substellar companion orbiting the young spectroscopic binary HD 143811 AB in the Scorpius–Centaurus association. The system lies at a Gaia DR3 distance of 136.85±0.38136.85 \pm 0.38 pc from a parallax of 7.3065±0.02047.3065 \pm 0.0204 mas, and the age adopted for planet characterization is 13±413 \pm 4 Myr. The companion was confirmed through multi-epoch high-contrast imaging and common proper motion analysis, and current analyses place it on a wide circumbinary orbit with a semi-major axis near $64$ au and a mass of order $5.6$–6.1 MJup6.1~M_\mathrm{Jup}, depending on the adopted atmospheric and evolutionary framework (Jones et al., 8 Sep 2025).

1. Host system and stellar architecture

HD 143811 AB is a double-lined spectroscopic binary (SB2) identified with HIP 78663 and associated with Sco–Cen. Probabilistic membership from BANYAN Σ\Sigma favors Upper Scorpius at 72.1%72.1\%, with a 27.8%27.8\% chance of Upper Centaurus–Lupus and 0.2%0.2\% field. The space motion is 7.3065±0.02047.3065 \pm 0.02040 km s7.3065±0.02047.3065 \pm 0.02041. The adopted age used in the planet analyses is 7.3065±0.02047.3065 \pm 0.02042 Myr, synthesized from typical subgroup ages, although a stellar SED fit yielded a model-dependent age of 7.3065±0.02047.3065 \pm 0.02043 Myr (Peck et al., 8 Sep 2025).

The inner binary is unusually compact among hosts of directly imaged planets. Its spectroscopic orbit has period 7.3065±0.02047.3065 \pm 0.02044 days, eccentricity 7.3065±0.02047.3065 \pm 0.02045, systemic velocity 7.3065±0.02047.3065 \pm 0.02046 km s7.3065±0.02047.3065 \pm 0.02047, and mass ratio 7.3065±0.02047.3065 \pm 0.02048. The radial-velocity semi-amplitudes are 7.3065±0.02047.3065 \pm 0.02049 km s13±413 \pm 40 and 13±413 \pm 41 km s13±413 \pm 42, and the argument of periastron for the primary is 13±413 \pm 43 deg. Combining the spectroscopic constraints with stellar evolutionary modeling gives component masses 13±413 \pm 44 and 13±413 \pm 45, with total mass 13±413 \pm 46. Kepler’s law then implies a binary semi-major axis 13±413 \pm 47 au.

The stellar properties are consistent with a near-equal-mass F+F pair. PHOENIX-based spectral fitting yielded 13±413 \pm 48, 13±413 \pm 49 K, $64$0 K, $64$1 km s$64$2, and $64$3 km s$64$4. A separate MIST+ATLAS9 SED fit to blended Gaia and 2MASS photometry gave $64$5, $64$6 mag, and somewhat hotter temperatures, $64$7 K and $64$8 K. The authors note that the SED-based temperatures are systematically higher than the spectral-fitting values and attribute this to extinction–temperature covariance and model systematics; the masses are described as more robust than $64$9.

The planet therefore resides in a hierarchical circumbinary configuration: an $5.6$0-day eccentric SB2 at $5.6$1 au, orbited by a giant planetary-mass companion at $5.6$2 au. This architecture is rare among directly imaged systems and provides the dynamical context for all subsequent interpretation.

2. Detection history and confirmation as a bound companion

The companion was initially detected in Gemini Planet Imager observations and later confirmed with Keck/NIRC2 and, in an independent analysis, with SPHERE. The primary discovery paper reports GPI integral-field spectroscopy in $5.6$3 band $5.6$4 on 2016 Apr 30 and 2019 Aug 11, together with Keck/NIRC2 $5.6$5-band imaging $5.6$6 on 2022 Jun 10. The GPI epochs comprised $5.6$7 s and $5.6$8 s integrations, while the NIRC2 data consisted of 175 frames of 60 coadds $5.6$9 0.5 s, for a total of 87.5 min (Jones et al., 8 Sep 2025).

Epoch Instrument / band Relative astrometry
2016 GPI 6.1 MJup6.1~M_\mathrm{Jup}0 6.1 MJup6.1~M_\mathrm{Jup}1 mas, 6.1 MJup6.1~M_\mathrm{Jup}2 deg
2019 GPI 6.1 MJup6.1~M_\mathrm{Jup}3 6.1 MJup6.1~M_\mathrm{Jup}4 mas, 6.1 MJup6.1~M_\mathrm{Jup}5 deg
2022 Keck/NIRC2 6.1 MJup6.1~M_\mathrm{Jup}6 6.1 MJup6.1~M_\mathrm{Jup}7 mas, 6.1 MJup6.1~M_\mathrm{Jup}8 deg

The 2016 GPI detection had signal-to-noise ratio 6.1 MJup6.1~M_\mathrm{Jup}9 using a forward model matched filter with an L-type template. Its measured Σ\Sigma0-band flux ratio was Σ\Sigma1, corresponding to apparent Σ\Sigma2, absolute Σ\Sigma3, and integrated Σ\Sigma4-band flux Σ\Sigma5 W mΣ\Sigma6 Σ\Sigma7mΣ\Sigma8. The 2019 GPI epoch had SNR Σ\Sigma9, flux ratio 72.1%72.1\%0, and integrated flux 72.1%72.1\%1 W m72.1%72.1\%2 72.1%72.1\%3m72.1%72.1\%4. The 2022 Keck/NIRC2 data gave an 72.1%72.1\%5-band flux ratio 72.1%72.1\%6, apparent 72.1%72.1\%7, absolute 72.1%72.1\%8, and flux 72.1%72.1\%9 W m27.8%27.8\%0 27.8%27.8\%1m27.8%27.8\%2.

Common proper motion was established by comparison with the Gaia DR3 proper motion and parallax of the host. The expected trajectory of a stationary background source is strongly inconsistent with the measured astrometry, and the observed separations disagree with the background track at 27.8%27.8\%3. Bound-companion orbital tracks generated by OFTI from the first epoch are consistent with the later epochs, confirming that HD 143811 AB b is gravitationally bound.

An independent SPHERE-based study reached the same conclusion with a longer temporal baseline. That analysis used 2016 and 2019 GPI data plus SPHERE/IRDIS and IFS observations obtained in 2025, reporting astrometry of 27.8%27.8\%4 mas, 27.8%27.8\%5 in 2016; 27.8%27.8\%6 mas, 27.8%27.8\%7 in 2019; and 27.8%27.8\%8 mas, 27.8%27.8\%9 in 2025. That work derived a relative proper motion 0.2%0.2\%0 mas yr0.2%0.2\%1, equivalent to 0.2%0.2\%2 km s0.2%0.2\%3 at projected separation 0.2%0.2\%4 au, below the escape speed 0.2%0.2\%5 km s0.2%0.2\%6, and reported a background false-alarm probability of 0.2%0.2\%7 (Squicciarini et al., 7 Sep 2025).

3. Orbital solution and circumbinary geometry

The planet’s orbit has been modeled with Orbits For The Impatient (OFTI). In the discovery analysis, the priors were uniform in 0.2%0.2\%8, 0.2%0.2\%9, 7.3065±0.02047.3065 \pm 0.020400, 7.3065±0.02047.3065 \pm 0.020401, 7.3065±0.02047.3065 \pm 0.020402, and 7.3065±0.02047.3065 \pm 0.020403, with parallax 7.3065±0.02047.3065 \pm 0.020404 mas and total mass 7.3065±0.02047.3065 \pm 0.020405. Using three relative astrometry epochs, the fitted posterior summary yielded semi-major axis 7.3065±0.02047.3065 \pm 0.020406 au, eccentricity 7.3065±0.02047.3065 \pm 0.020407, inclination 7.3065±0.02047.3065 \pm 0.020408, and period 7.3065±0.02047.3065 \pm 0.020409 yr. The posterior medians were 7.3065±0.02047.3065 \pm 0.020410 au, 7.3065±0.02047.3065 \pm 0.020411, 7.3065±0.02047.3065 \pm 0.020412, 7.3065±0.02047.3065 \pm 0.020413 yr, 7.3065±0.02047.3065 \pm 0.020414, 7.3065±0.02047.3065 \pm 0.020415, and 7.3065±0.02047.3065 \pm 0.020416, with broad credible intervals for the angular elements (Jones et al., 8 Sep 2025).

The period estimate is consistent with Kepler’s third law. In astronomical units and solar masses, the relation is 7.3065±0.02047.3065 \pm 0.020417. Using 7.3065±0.02047.3065 \pm 0.020418 au and 7.3065±0.02047.3065 \pm 0.020419 gives 7.3065±0.02047.3065 \pm 0.020420 yr. The broad uncertainty in 7.3065±0.02047.3065 \pm 0.020421, rather than the small uncertainty in 7.3065±0.02047.3065 \pm 0.020422, dominates the reported period range.

Projected physical separations derived from 7.3065±0.02047.3065 \pm 0.020423 are 7.3065±0.02047.3065 \pm 0.020424 au in 2016, 7.3065±0.02047.3065 \pm 0.020425 au in 2019, and 7.3065±0.02047.3065 \pm 0.020426 au in 2022, fully consistent with a fitted semi-major axis near 7.3065±0.02047.3065 \pm 0.020427 au and modest eccentricity. A parallel orbit analysis using GPI plus SPHERE astrometry found closely similar values: 7.3065±0.02047.3065 \pm 0.020428 au, 7.3065±0.02047.3065 \pm 0.020429, 7.3065±0.02047.3065 \pm 0.020430 deg, and 7.3065±0.02047.3065 \pm 0.020431 yr. The agreement between the two studies indicates that, with present astrometric coverage, the outer orbit is already constrained to be wide, low-to-moderately eccentric, and mostly face-on, while 7.3065±0.02047.3065 \pm 0.020432, 7.3065±0.02047.3065 \pm 0.020433, and 7.3065±0.02047.3065 \pm 0.020434 remain weakly constrained (Squicciarini et al., 7 Sep 2025).

The relation between the inner and outer orbital planes remains unresolved. The host-binary characterization paper inferred a binary inclination 7.3065±0.02047.3065 \pm 0.020435 deg, with a spectroscopically degenerate retrograde alternative at 7.3065±0.02047.3065 \pm 0.020436 deg. Because spectroscopy does not constrain the inner node 7.3065±0.02047.3065 \pm 0.020437, and the planetary 7.3065±0.02047.3065 \pm 0.020438 posterior remains broad, the mutual inclination cannot yet be computed. The current estimates overlap within roughly 7.3065±0.02047.3065 \pm 0.020439, so coplanarity cannot be ruled out, but neither can it be established (Peck et al., 8 Sep 2025).

4. Photometry, spectroscopy, and atmospheric modeling

The planet’s spectro-photometric characterization currently rests on a GPI 7.3065±0.02047.3065 \pm 0.020440-band spectrum and thermal-infrared photometry. The 2016 GPI 7.3065±0.02047.3065 \pm 0.020441-band spectrum was extracted with KLIP-FM forward modeling and Bayesian PSF fitting, calibrated using satellite spots and the stellar SED, and was described as exhibiting strongly correlated spectral noise typical for GPI but robustly measured flux. The 2022 NIRC2 7.3065±0.02047.3065 \pm 0.020442 photometry was calibrated to stellar 7.3065±0.02047.3065 \pm 0.020443 through the system SED model and was reported as consistent with a cool, substellar companion (Jones et al., 8 Sep 2025).

In color–magnitude space, the combined measurements yield 7.3065±0.02047.3065 \pm 0.020444 mag and 7.3065±0.02047.3065 \pm 0.020445 mag. On a color–magnitude diagram, the object lies near late-L/early-T field objects and close to the directly imaged companions HD 206893 B and 2M 1207 B. The discovery paper states that this positioning suggests a very red, dusty atmosphere may be plausible. This is explicitly framed as an interpretation of the photometric placement rather than a direct measurement.

Atmospheric model comparison was carried out against 842 PHOENIX synthetic stellar spectra and 9575 Exo-REM exoplanet atmosphere models, including clouds and disequilibrium chemistry, across 7.3065±0.02047.3065 \pm 0.020446 K. The initial goodness-of-fit statistic was

7.3065±0.02047.3065 \pm 0.020447

For the GPI spectrum alone, PHOENIX and Exo-REM both matched moderately well, with PHOENIX reaching a maximum 7.3065±0.02047.3065 \pm 0.020448, and Exo-REM slightly preferred. For the joint GPI+NIRC2 fit, Exo-REM was strongly preferred, with best-fit 7.3065±0.02047.3065 \pm 0.020449, while the PHOENIX templates were rejected. A representative best-fit Exo-REM model had 7.3065±0.02047.3065 \pm 0.020450 K, 7.3065±0.02047.3065 \pm 0.020451, 7.3065±0.02047.3065 \pm 0.020452 dex, and 7.3065±0.02047.3065 \pm 0.020453, but the paper explicitly states that 7.3065±0.02047.3065 \pm 0.020454, metallicity, and C/O are not constrained by the data.

A more complete inference used nested sampling with dynesty and a Gaussian-process noise model for the correlated GPI spectral noise. That analysis yielded 7.3065±0.02047.3065 \pm 0.020455 K, radius 7.3065±0.02047.3065 \pm 0.020456, and luminosity 7.3065±0.02047.3065 \pm 0.020457. Surface gravity, metallicity, and C/O ratio remained unconstrained across the Exo-REM grid. Under hot-start evolutionary models of Baraffe et al. (2003), and adopting age 7.3065±0.02047.3065 \pm 0.020458 Myr, the derived luminosity corresponds to a mass 7.3065±0.02047.3065 \pm 0.020459, with an evolutionary radius 7.3065±0.02047.3065 \pm 0.020460, consistent with the atmospheric-fit radius within uncertainties.

A separate GPI+SPHERE analysis used the HADES coupled atmosphere–interior grid with an affine-invariant MCMC, including SPHERE/IRDIS 7.3065±0.02047.3065 \pm 0.020461 photometry and SPHERE/IFS 7.3065±0.02047.3065 \pm 0.020462 upper limits. That work reported 7.3065±0.02047.3065 \pm 0.020463 mag, 7.3065±0.02047.3065 \pm 0.020464 mag, 7.3065±0.02047.3065 \pm 0.020465 mag, and a 7.3065±0.02047.3065 \pm 0.020466 non-detection limit 7.3065±0.02047.3065 \pm 0.020467 mag, implying 7.3065±0.02047.3065 \pm 0.020468 mag. It derived intrinsic temperature 7.3065±0.02047.3065 \pm 0.020469 K, mass 7.3065±0.02047.3065 \pm 0.020470, radius 7.3065±0.02047.3065 \pm 0.020471, atmospheric metallicity 7.3065±0.02047.3065 \pm 0.020472 dex, cloud sedimentation efficiency 7.3065±0.02047.3065 \pm 0.020473, and vigorous vertical mixing 7.3065±0.02047.3065 \pm 0.020474. It interpreted the H-band shape and red 7.3065±0.02047.3065 \pm 0.020475 color as consistent with strong 7.3065±0.02047.3065 \pm 0.020476 absorption and muted 7.3065±0.02047.3065 \pm 0.020477 near the L/T transition in a young, cloudy, chemically disequilibrium atmosphere (Squicciarini et al., 7 Sep 2025).

5. Dynamical regime and formation significance

HD 143811(AB)b belongs to a small class of directly imaged planets orbiting binaries rather than single stars. The discovery paper explicitly places it among few such systems and compares it with HD 106906 b, b Cen b, and WISPIT 1bc, while also noting borderline planetary-mass companions around binaries such as Ross 458C, SR 12 C, 2MASS J01033563–5515561 (AB)b, and ROXs 42B b. Within that group, HD 143811 AB b is distinguished by a moderate separation of about 7.3065±0.02047.3065 \pm 0.020478 au, modest eccentricity near 7.3065±0.02047.3065 \pm 0.020479, and youth at 7.3065±0.02047.3065 \pm 0.020480 Myr (Jones et al., 8 Sep 2025).

The inner binary is sufficiently tight that the planet resides deep in the stable circumbinary regime. The host-binary study invokes the Holman & Wiegert (1999) critical semimajor axis framework, writing 7.3065±0.02047.3065 \pm 0.020481 with 7.3065±0.02047.3065 \pm 0.020482. Using 7.3065±0.02047.3065 \pm 0.020483 au, 7.3065±0.02047.3065 \pm 0.020484, and 7.3065±0.02047.3065 \pm 0.020485, it concludes that 7.3065±0.02047.3065 \pm 0.020486 remains orders of magnitude smaller than the planet’s 7.3065±0.02047.3065 \pm 0.020487 au orbit, so the system is comfortably within the stable circumbinary regime. A separate analysis, adopting a compact binary with 7.3065±0.02047.3065 \pm 0.020488 au, similarly states that the planet’s orbit at 7.3065±0.02047.3065 \pm 0.020489 au comfortably exceeds any plausible critical stability radius (Peck et al., 8 Sep 2025).

The formation interpretation remains open. The discovery paper states that the planet’s moderate separation, modest eccentricity, and youthful age provide a valuable data point for understanding planet formation and dynamical evolution in circumbinary environments where core accretion is challenged but not excluded. The host-binary study adds that the tight inner binary likely truncated and structured the primordial circumbinary disk, setting initial conditions for planet formation and migration, and that the large binary–planet separation ratio of 7.3065±0.02047.3065 \pm 0.020490 is consistent with long-term stability. It further states that the observed configuration suggests formation in a circumbinary disk at tens of au, or post-formation outward migration or scattering, but does not distinguish among these scenarios.

The system may also include circum-system debris. Spitzer detections indicate a cold debris belt with fractional infrared luminosity between 7.3065±0.02047.3065 \pm 0.020491 and 7.3065±0.02047.3065 \pm 0.020492, inferred temperatures of 7.3065±0.02047.3065 \pm 0.020493 K, and blackbody radii of 7.3065±0.02047.3065 \pm 0.020494 au, although the host-binary paper emphasizes that grain properties and the then-unknown binary nature can shift radii outward and that a 70 7.3065±0.02047.3065 \pm 0.020495m non-detection was reported. This suggests that the architecture may eventually permit joint analysis of binary, planet, and disk dynamics, but current disk parameters remain too uncertain for strong conclusions.

6. Open questions and future characterization

Several central parameters remain weakly constrained. For the planet, the current three-epoch orbit fits already limit 7.3065±0.02047.3065 \pm 0.020496, 7.3065±0.02047.3065 \pm 0.020497, 7.3065±0.02047.3065 \pm 0.020498, and 7.3065±0.02047.3065 \pm 0.020499, but 13±413 \pm 400, 13±413 \pm 401, and 13±413 \pm 402 remain broad. For the inner binary, spectroscopy fixes the period, eccentricity, and mass ratio with high precision but cannot determine the longitude of the ascending node, and radial velocities alone leave a mirror degeneracy between prograde and retrograde orbital inclinations. Consequently, the full three-dimensional architecture and mutual inclination are not yet known (Peck et al., 8 Sep 2025).

Continued high-precision relative astrometry is expected to refine the planetary orbital parameters 13±413 \pm 403, 13±413 \pm 404, 13±413 \pm 405, 13±413 \pm 406, 13±413 \pm 407, and 13±413 \pm 408, and may eventually permit dynamical mass constraints when combined with absolute astrometry and a well-determined binary orbit. The discovery paper specifically notes that Hipparcos–Gaia accelerations could contribute once the host-binary solution is improved. The host-binary paper recommends VLTI/GRAVITY, since the binary apastron corresponds to approximately 13±413 \pm 409 mas at 13±413 \pm 410 pc, potentially allowing recovery of the inner 13±413 \pm 411 and true 13±413 \pm 412, and cites Gaia DR4 individual-scan astrometry as another route to constraining the 13±413 \pm 413 mas photocenter wobble.

Additional spectroscopy across broader wavelength coverage is a major priority. The discovery analysis recommends further observations in 13±413 \pm 414, 13±413 \pm 415, and 13±413 \pm 416 bands and medium-resolution spectroscopy to constrain 13±413 \pm 417, 13±413 \pm 418, and CO, as well as cloud properties, metallicity, and C/O, while reducing the impact of correlated-noise systematics. It explicitly recommends follow-up with upgraded high-contrast imagers such as Gemini Planet Imager 2.0, along with complementary Keck observations. The GPI+SPHERE study further identifies JWST NIRSpec/G395M and MIRI/LRS–MRS as promising for resolving molecular and cloud signatures across 13±413 \pm 419m, and also points to variability, polarimetry, and ELT-class high-resolution spectroscopy as routes to probing cloud patchiness, particle size, winds, rotation, and disequilibrium chemistry (Squicciarini et al., 7 Sep 2025).

The scientific importance of these follow-up programs is twofold. First, they can test whether the apparent consistency of the binary and planetary inclinations reflects genuine coplanarity expected from formation in a common circumbinary disk. Second, they can determine whether HD 143811(AB)b is best interpreted as a planet formed in situ in a circumbinary disk, a product of migration, or an object whose present orbit encodes later dynamical evolution. Present data suggest several of these possibilities, but do not yet isolate a unique formation pathway.

Topic to Video (Beta)

No one has generated a video about this topic yet.

Whiteboard

No one has generated a whiteboard explanation for this topic yet.

Follow Topic

Get notified by email when new papers are published related to HD 143811(AB)b.