Relative Astrometry Research Articles

Orbits and masses in two triple systems

Abstract In an effort to determine accurate orbital and physical properties of a large number of bright stars, a method was developed to fit simultaneously stellar parameters (masses, luminosities, effective temperatures), distance, and orbits to the available data on multiple systems, namely the combined and differential photometry, positional measurements, radial velocities (RVs), accelerations, etc. The method is applied to a peculiar resolved triple system HIP 86286. The masses of its components estimated using observations and standard relations are 1.3, 0.9, and 0.9 M⊙; the main star is a G8IV subgiant, while its two companions are main-sequence dwarfs. The inner and outer orbital periods are 35 and 287 years, respectively, and the orbits are nearly coplanar. The second system, HIP 117258, is an accelerating star with a resolved companion; its 35.7-yr orbit based on relative astrometry and precise RVs yields the secondary mass of 0.95 M⊙, much larger than inferred from the photometry. The apparent paradox is explained by assuming that the secondary is a close pair of M-type dwarfs with yet unknown period.

Subaru/CHARIS High-resolution Mode Spectroscopy of the Brown Dwarf Companion HD 33632 Ab

Brown dwarfs (BDs) 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.

The Keck-HGCA Pilot Survey II: Direct imaging discovery of HD 63754 B, a ∼20 au massive companion near the hydrogen burning limit

Abstract 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 Hipparcosand 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} \rm {M_{\rm Jup}}$, an eccentricity of ${0.260}_{-0.059}^{+0.065}$, and a nearly face-on inclination of $174.\!\!^\circ 81_{-0.50}^{+0.48}$. For HD 63754 B, we obtain an L′ band absolute magnitude of L′ = 11.39 ± 0.06 mag, from which we infer a bolometric luminosity of $\mathrm{log(L_{\rm bol}/\rm {\mathrm{L}_{\odot }})= -4.55 \pm 0.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.

PHANGS-JWST: Data-processing Pipeline and First Full Public Data Release

The exquisite angular resolution and sensitivity of JWST are opening a new window for our understanding of the Universe. In nearby galaxies, JWST observations are revolutionizing our understanding of the first phases of star formation and the dusty interstellar medium. Nineteen local galaxies spanning a range of properties and morphologies across the star-forming main sequence have been observed as part of the PHANGS-JWST Cycle 1 Treasury program at spatial scales of ∼5–50 pc. Here, we describe pjpipe, an image-processing pipeline developed for the PHANGS-JWST program that wraps around and extends the official JWST pipeline. We release this pipeline to the community as it contains a number of tools generally useful for JWST NIRCam and MIRI observations. Particularly for extended sources, pjpipe products provide significant improvements over mosaics from the MAST archive in terms of removing instrumental noise in NIRCam data, background flux matching, and calibration of relative and absolute astrometry. We show that slightly smoothing F2100W MIRI data to 0.″9 (degrading the resolution by about 30%) reduces the noise by a factor of ≈3. We also present the first public release (DR1.1.0) of the pjpipe processed eight-band 2–21 μm imaging for all 19 galaxies in the PHANGS-JWST Cycle 1 Treasury program. An additional 55 galaxies will soon follow from a new PHANGS-JWST Cycle 2 Treasury program.

HIP 65426 is a High-frequency Delta Scuti Pulsator in Plausible Spin–Orbit Alignment with its Directly Imaged Exoplanet

HIP 65426 hosts a young giant planet that has become the first exoplanet directly imaged with JWST. Using time-series photometry from the Transiting Exoplanet Survey Satellite (TESS), we classify HIP 65426 as a high-frequency δ Scuti pulsator with a possible large-frequency separation of Δν = 7.23 ± 0.02 cycles day−1. We check the TESS data for pulsation-timing variations and use the nondetection to estimate a 95% dynamical mass upper limit of 12.8 M Jup for HIP 65426 b. We also identify a low-frequency region of signal that we interpret as stellar latitudinal differential rotation with two rapid periods of 7.85 ± 0.08 hr and 6.67 ± 0.04 hr. We use our TESS rotation periods together with published values of radius and vsini to jointly measure the inclination of HIP 65426 to i⋆=107−11+12 °. Our stellar inclination is consistent with the orbital inclination of HIP 65426 b ( 108−3+6 °) at the 68% percent level based on our orbit fit using published relative astrometry. The lack of significant evidence for spin–orbit misalignment in the HIP 65426 system supports an emerging trend consistent with preferential alignment between imaged long-period giant planets and their host stars.

Astrometric detection of a Neptune-mass candidate planet in the nearest M-dwarf binary system GJ65 with VLTI/GRAVITY

The detection of low-mass planets orbiting the nearest stars is a central stake of exoplanetary science, as they can be directly characterized much more easily than their distant counterparts. Here, we present the results of our long-term astrometric observations of the nearest binary M-dwarf Gliese 65 AB (GJ65), located at a distance of only 2.67 pc. We monitored the relative astrometry of the two components from 2016 to 2023 with the VLTI/GRAVITY interferometric instrument. We derived highly accurate orbital parameters for the stellar system, along with the dynamical masses of the two red dwarfs. The GRAVITY measurements exhibit a mean accuracy per epoch of 50−60 ms in 1.5 h of observing time using the 1.8 m Auxiliary Telescopes. The residuals of the two-body orbital fit enable us to search for the presence of companions orbiting one of the two stars (S-type orbit) through the reflex motion they imprint on the differential A–B astrometry. We detected a Neptune-mass candidate companion with an orbital period of p = 156 ± 1 d and a mass of mp = 36 ± 7 M⊕. The best-fit orbit is within the dynamical stability region of the stellar pair. It has a low eccentricity, e = 0.1 − 0.3, and the planetary orbit plane has a moderate-to-high inclination of i > 30° with respect to the stellar pair, with further observations required to confirm these values. These observations demonstrate the capability of interferometric astrometry to reach microarcsecond accuracy in the narrow-angle regime for planet detection by reflex motion from the ground. This capability offers new perspectives and potential synergies with Gaia in the pursuit of low-mass exoplanets in the solar neighborhood.

Direct-imaging Discovery of a Substellar Companion Orbiting the Accelerating Variable Star HIP 39017

We present the direct-imaging discovery of a substellar companion (a massive planet or low-mass brown dwarf) to the young, γ Doradus (γ Dor)-type variable star HIP 39017 (HD 65526). The companion’s SCExAO/CHARIS JHK (1.1–2.4 μm) spectrum and Keck/NIRC2 L′ photometry indicate that it is an L/T transition object. A comparison of the JHK+L ′ spectrum to several atmospheric model grids finds a significantly better fit to cloudy models than cloudless models. Orbit modeling with relative astrometry and precision stellar astrometry from Hipparcos and Gaia yields a semimajor axis of 23.8−6.1+8.7 au, a dynamical companion mass of 30−12+31 M J, and a mass ratio of ∼1.9%, properties most consistent with low-mass brown dwarfs. However, its mass estimated from luminosity models is a lower ∼13.8 M J due to an estimated young age (≲115 Myr); using a weighted posterior distribution informed by conservative mass constraints from luminosity evolutionary models yields a lower dynamical mass of 23.6−7.4+9.1 M J and a mass ratio of ∼1.4%. Analysis of the host star’s multifrequency γ Dor-type pulsations, astrometric monitoring of HIP 39017 b, and Gaia Data Release 4 astrometry of the star will clarify the system age and better constrain the mass and orbit of the companion. This discovery further reinforces the improved efficiency of targeted direct-imaging campaigns informed by long-baseline, precision stellar astrometry.

Beyond Point Masses. II. Non-Keplerian Shape Effects Are Detectable in Several TNO Binaries

About 40 trans-Neptunian binaries (TNBs) have fully determined orbits with about 10 others being solved except for breaking the mirror ambiguity. Despite decades of study, almost all TNBs have only ever been analyzed with a model that assumes perfect Keplerian motion (e.g., two point masses). In reality, all TNB systems are non-Keplerian due to nonspherical shapes, possible presence of undetected system components, and/or solar perturbations. In this work, we focus on identifying candidates for detectable non-Keplerian motion based on sample of 45 well-characterized binaries. We use MultiMoon, a non-Keplerian Bayesian inference tool, to analyze published relative astrometry allowing for nonspherical shapes of each TNB system’s primary. We first reproduce the results of previous Keplerian fitting efforts with MultiMoon, which serves as a comparison for the non-Keplerian fits and confirms that these fits are not biased by the assumption of a Keplerian orbit. We unambiguously detect non-Keplerian motion in eight TNB systems across a range of primary radii, mutual orbit separations, and system masses. As a proof of concept for non-Keplerian fitting, we perform detailed fits for (66652) Borasisi-Pabu, possibly revealing a J 2 ≈ 0.44, implying Borasisi (and/or Pabu) may be a contact binary or an unresolved compact binary. However, full confirmation of this result will require new observations. This work begins the next generation of TNB analyses that go beyond the point mass assumption to provide unique and valuable information on the physical properties of TNBs with implications for their formation and evolution.

The orbit of HD 142527 B is too compact to explain many of the disc features

Context. HD 142527 A is a young and massive Herbig Ae/Be star surrounded by a highly structured disc. The disc shows numerous morphological structures, such as spiral arms, a horseshoe region of dust emission, a set of shadows cast by an inner disc on the outer disc, and a large cavity extending from ≃30 au to ≃130 au. HD 142527 A also has a lower mass companion, HD 142527 B (M = 0.13 ± 0.03 M⊙), which is thought to be responsible for most of the structures observed in the surrounding disc. Aims. We aim to fully constrain the orbit of HD 142527 B and determine whether the binary alone is truly responsible for the observed morphology of the HD 142527 disc. Methods. We gathered VLTI/GRAVITY observations of HD 142527, either from our own programmes or from the ESO archive. We used this inhomogeneous set of data to extract a total of seven high-precision measurements of the relative astrometry between HD 142527 A and B, spread from mid-2017 to early 2021. Combined with what is available in the literature, these new measurements offer a total of 9 yr of astrometric monitoring on HD 142527. We used orbit fitting tools to determine the orbital parameters of HD 142527 B, and used them as inputs for a 3D hydrodynamical model of the disc to determine whether or not the binary is able to create the structures observed in the disc. Results. Our VLTI/GRAVITY astrometry gives excellent constraints on the orbit of HD 142527 B. We show that the secondary is following an orbit of semi-major axis a = 10.80 ± 0.22 au, with moderate eccentricity (e = 0.47 ± 0.01), and has recently passed its periapsis (τ = 2020.42). With such a compact orbit, we show that HD 142527 B can only generate a gap and spiral arms of ~30 au in the disc, which is much smaller than what is revealed by observations. Conclusions. Even from a theoretical standpoint, the observed cavity size of ~ 100 au far exceeds even the most generous predictions for a companion like HD 142527 B on such a compact orbit. Thus, we conclude that the low-mass companion cannot be solely responsible for the observed morphology of the disc surrounding the system.

Orbit of the Patroclus–Menoetius Binary System and Predictions for the 2024/2025 Mutual Events Season

We report on the ephemeris development for Menoetius, the satellite of Patroclus. Our data set consisted of ground–based and Hubble Space Telescope relative astrometry, as well as 42 lightcurves from the mutual events seasons in 2007, 2012, and 2017/2018. Our dynamical model included the effects of oblate, nonspherical shapes of the components, and we assumed that Menoetius contained ∼22% of the system’s mass. We numerically integrated the equations of motion and obtained a set of dynamical parameters that fit the data. We report the fit results in terms of residuals, state vectors, orbital elements and their 1σ uncertainties. The mean osculating semimajor axis is a = 692.5 ± 4.0 km, the mean eccentricity is e = 0.004 ± 0.004, and the International Celestial Reference Frame pole direction in R.A. and decl. is R. A. = 178.0 ± 0.5 deg, δ = −74.1 ± 0.2 deg. We determined the siderial orbital period of P = 4.282753 ± 0.000023 days. The fit yielded the system GM = 0.0950 ± 0.0012 km3 s−2, which, in combination with the system volume determined from the stellar occultation and the assumed volume uncertainty of 20%, suggests a system bulk density of 1.05 ± 0.21 g cm−3(1σ). The next season of mutual events starts in February of 2024 and lasts until January of 2025. The Patroclus system is in opposition for the observers on Earth in late September and is suitable for observations of the mutual events with an edge-on geometry in October.

Dynamical Architectures of S-type Transiting Planets in Binaries. I. Target Selection Using Hipparcos and Gaia Proper Motion Anomalies* * Also referred to as Hipparcos and Gaia astrometric accelerations.

The effect of stellar multiplicity on planetary architecture and orbital dynamics provides an important context for exoplanet demographics. We present a volume-limited catalog of up to 300 pc of 66 stars hosting planets and planet candidates from Kepler, K2, and TESS with significant Hipparcos-Gaia proper motion anomalies, which indicates the presence of companions. We assess the reliability of each transiting planet candidate using ground-based follow-up observations, and find that the TESS Objects of Interest (TOIs) with significant proper anomalies show nearly four times more false positives due to eclipsing binaries compared to TOIs with marginal proper anomalies. In addition, we find tentative evidence that orbital periods of planets orbiting TOIs with significant proper anomalies are shorter than those orbiting TOIs without significant proper anomalies, consistent with the scenario that stellar companions can truncate planet-forming disks. Furthermore, TOIs with significant proper anomalies exhibit lower Gaia differential velocities in comparison to field stars with significant proper anomalies, suggesting that planets are more likely to form in binary systems with low-mass substellar companions or stellar companions at wider separation. Finally, we characterize the three-dimensional architecture of LTT 1445 ABC using radial velocities, absolute astrometry from Gaia and Hipparcos, and relative astrometry from imaging. Our analysis reveals that LTT 1445 is a nearly flat system, with a mutual inclination of ∼2.°88 between the orbit of BC around A and that of C around B. This coplanarity may explain why multiple planets around LTT 1445 A survive in the dynamically hostile environments of this system.

CD –27°11535: Evidence for a Triple System in the β Pictoris Moving Group

We present new spatially resolved astrometry and photometry of the CD –27°11535 system, a member of the β Pictoris moving group consisting of two resolved K-type stars on a ∼20 yr orbit. We fit an orbit to relative astrometry measured from NIRC2, GPI, and archival NaCo images, in addition to literature measurements. However, the total mass inferred from this orbit is significantly discrepant from that inferred from stellar evolutionary models using the luminosity of the two stars. We explore two hypotheses that could explain this discrepant mass sum: a discrepant parallax measurement from Gaia due to variability, and the presence of an additional unresolved companion to one of the two components. We find that the ∼20 yr orbit could not bias the parallax measurement, but that variability of the components could produce a large-amplitude astrometric motion, an effect that cannot be quantified exactly without the individual Gaia measurements. The discrepancy could also be explained by an additional star in the system. We jointly fit the astrometric and photometric measurements of the system to test different binary and triple architectures for the system. Depending on the set of evolutionary models used, we find an improved goodness of fit for a triple system architecture that includes a low-mass (M = 0.177 ± 0.055 M ⊙) companion to the primary star. Further studies of this system will be required in order to resolve this discrepancy, either by refining the parallax measurement with a more complex treatment of variability-induced astrometric motion or by detecting a third companion.

ELemental abundances of Planets and brown dwarfs Imaged around Stars (ELPIS). I. Potential Metal Enrichment of the Exoplanet AF Lep b and a Novel Retrieval Approach for Cloudy Self-luminous Atmospheres

AF Lep A+b is a remarkable planetary system hosting a gas-giant planet that has the lowest dynamical mass among directly imaged exoplanets. We present an in-depth analysis of the atmospheric composition of the star and planet to probe the planet’s formation pathway. Based on new high-resolution spectroscopy of AF Lep A, we measure a uniform set of stellar parameters and elemental abundances (e.g., [Fe/H] = −0.27 ± 0.31 dex). The planet’s dynamical mass ( M Jup) and orbit are also refined using published radial velocities, relative astrometry, and absolute astrometry. We use petitRADTRANS to perform chemically consistent atmospheric retrievals for AF Lep b. The radiative–convective equilibrium temperature profiles are incorporated as parameterized priors on the planet’s thermal structure, leading to a robust characterization for cloudy self-luminous atmospheres. This novel approach is enabled by constraining the temperature–pressure profiles via the temperature gradient , a departure from previous studies that solely modeled the temperature. Through multiple retrievals performed on different portions of the 0.9–4.2 μm spectrophotometry, along with different priors on the planet’s mass and radius, we infer that AF Lep b likely possesses a metal-enriched atmosphere ([Fe/H] > 1.0 dex). AF Lep b’s potential metal enrichment may be due to planetesimal accretion, giant impacts, and/or core erosion. The first process coincides with the debris disk in the system, which could be dynamically excited by AF Lep b and lead to planetesimal bombardment. Our analysis also determines T eff ≈ 800 K, dex, and the presence of silicate clouds and disequilibrium chemistry in the atmosphere. Straddling the L/T transition, AF Lep b is thus far the coldest exoplanet with suggested evidence of silicate clouds.

Direct Measurements of Stellar Masses with the Habitable World Observatory

Background stars within ∼1″ of a nearby massive object will have their centroids shift by about 1 mas due to microlensing. This provides mass measurements of stellar mass objects in the solar neighborhood. Habitable Worlds Observatory (HWO) targets toward the Galactic Bulge will have a sufficient density of background stars that HWO will be able to measure their masses independently of stellar evolutionary models, which will allow more accurate habitable histories of their exoplanets to be evaluated. This opportunity is maximized for designs of HWO with very high contrasts (Δmag ≳25) and the capability for precise relative astrometry with respect to the target star. HWO can apply this technique to any nearby stellar mass object in the direction of the Galactic Bulge, enhancing its science case and that of precursor missions like Roman.

VLTI/GRAVITY Observations and Characterization of the Brown Dwarf Companion HD 72946 B

Tension remains between the observed and modeled properties of substellar objects, but objects in binary orbits, with known dynamical masses, can provide a way forward. HD 72946 B is a recently imaged brown dwarf companion to a nearby, solar-type star. We achieve ∼100 μas relative astrometry of HD 72946 B in the K band using VLTI/GRAVITY, unprecedented for a benchmark brown dwarf. We fit an ensemble of measurements of the orbit using orbitize! and derive a strong dynamical mass constraint M B = 69.5 ± 0.5 M Jup assuming a strong prior on the host star mass M A = 0.97 ± 0.01 M ⊙ from an updated stellar analysis. We fit the spectrum of the companion to a grid of self-consistent BT-Settl-CIFIST model atmospheres, and perform atmospheric retrievals using petitRADTRANS. A dynamical mass prior only marginally influences the sampled distribution of effective temperature, but has a large influence on the surface gravity and radius, as expected. The dynamical mass alone does not strongly influence retrieved pressure–temperature or cloud parameters within our current retrieval setup. Independently of the cloud prescription and prior assumptions, we find agreement within ±2σ between the C/O of the host (0.52 ± 0.05) and brown dwarf (0.43–0.63), as expected from a molecular cloud collapse formation scenario, but our retrieved metallicities are implausibly high (0.6–0.8) in light of the excellent agreement of the data with the solar-abundance model grid. Future work on our retrieval framework will seek to resolve this tension. Additional study of low surface gravity objects is necessary to assess the influence of a dynamical mass prior on atmospheric analysis.

Visual Orbits and Alignments of Planet-hosting Binary Systems

Roughly half of Solar-type planet hosts have stellar companions, so understanding how these binary companions affect the formation and evolution of planets is an important component to understanding planetary systems overall. Measuring the dynamical properties of planet host binaries enables a valuable test of planet formation in multistar systems and requires knowledge of the binary orbital parameters. Using high-resolution imaging, we have measured the relative astrometry and visual orbits of 13 binary systems where one of the stars is known to host a transiting exoplanet. Our results indicate that the mutual inclination between the orbits of the binary hosts and the transiting planets are well aligned. Our results for close binary systems (a < 100 au) complement past work for wide planet host binaries from Gaia.

Micro-arcsecond Astrometry Technology: Detector and Field Distortion Calibration

Microarcsecond (μas) astrometry is an indispensable technique to detect earth-like exoplanets, fully characterize exoplanetary orbits, and measure their masses—information critical for assessing their habitability. Highly accurate astrometric measurements can also probe the nature of dark matter, the early universe, black holes, and neutron stars, thus providing unique data for new astrophysics. This paper presents technologies of calibrating detectors and field distortions for achieving narrow field μas relative astrometry with a focal plane array detector on a 6 m telescope.

Dynamical masses and ages of Sirius-like systems

ABSTRACT We measure precise orbits and dynamical masses and derive age constraints for six confirmed and one candidate Sirius-like systems, including the Hyades member HD 27483. Our orbital analysis incorporates radial velocities, relative astrometry, and Hipparcos–Gaia astrometric accelerations. We constrain the main-sequence lifetime of a white dwarf’s progenitor from the remnant’s dynamical mass and semi-empirical initial–final mass relations and infer the cooling age from mass and effective temperature. We present new relative astrometry of HD 27483 B from Keck/NIRC2 observations and archival Hubble Space Telescope data, and obtain the first dynamical mass of ${0.798}_{-0.041}^{+0.10}$ M⊙, and an age of ${450}_{-180}^{+570}$ Myr, consistent with previous age estimates of Hyades. We also measure precise dynamical masses for HD 114174 B (0.591 ± 0.011 M⊙) and HD 169889 B (${0.526}_{-0.037}^{+0.039}$ M⊙), but their age precisions are limited by their uncertain temperatures. For HD 27786 B, the unusually small mass of 0.443 ± 0.012 M⊙ suggests a history of rapid mass-loss, possibly due to binary interaction in its progenitor’s asymtotic giant branch phase. The orbits of HD 118475 and HD 136138 from our radial velocity fitting are overall in good agreement with Gaia DR3 astrometric two-body solutions, despite moderate differences in the eccentricity and period of HD 136138. The mass of ${0.580}_{-0.039}^{+0.052}$ M⊙ for HD 118475 B and a speckle imaging non-detection confirms that the companion is a white dwarf. Our analysis shows examples of a rich number of precise WD dynamical mass measurements enabled by Gaia DR3 and later releases, which will improve empirical calibrations of the white dwarf initial–final mass relation.

Dynamical Mass of the Young Brown Dwarf Companion PZ Tel B

Dynamical masses of giant planets and brown dwarfs are critical tools for empirically validating substellar evolutionary models and their underlying assumptions. We present a measurement of the dynamical mass and an updated orbit of PZ Tel B, a young brown dwarf companion orbiting a late-G member of the β Pic moving group. PZ Tel A exhibits an astrometric acceleration between Hipparcos and Gaia EDR3, which enables the direct determination of the companion’s mass. We have also acquired new Keck/NIRC2 adaptive optics imaging of the system, which increases the total baseline of relative astrometry to 15 yr. Our joint orbit fit yields a dynamical mass of , semimajor axis of , eccentricity of , and inclination of . The companion’s mass is consistent within 1.1σ of predictions from four grids of hot-start evolutionary models. The joint orbit fit also indicates a more modest eccentricity of PZ Tel B than previous results. PZ Tel joins a small number of young (<200 Myr) systems with benchmark substellar companions that have dynamical masses and precise ages from moving group membership.

Surveying nearby brown dwarfs with HGCA: direct imaging discovery of a faint, high-mass brown dwarf orbiting HD 176535 A

ABSTRACT Brown dwarfs with well-measured masses, ages, and luminosities provide direct benchmark tests of substellar formation and evolutionary models. We report the first results from a direct imaging survey aiming to find and characterize substellar companions to nearby accelerating stars with the assistance of the Hipparcos–Gaia Catalog of Accelerations (HGCA). In this paper, we present a joint high-contrast imaging and astrometric discovery of a substellar companion to HD 176535 A, a K3.5V main-sequence star aged approximately $3.59_{-1.15}^{+0.87}$ Gyr at a distance of 36.99 ± 0.03 pc. In advance of our high-contrast imaging observations, we combined precision High Accuracy Radial velocity Planet Searcher (HARPS) Radial Velocities (RVs) and HGCA astrometry to predict the potential companion’s location and mass. We thereafter acquired two nights of KeckAO/NIRC2 direct imaging observations in the L′ band, which revealed a companion with a contrast of $\Delta L^{\prime }_p = 9.20\pm 0.06$ mag at a projected separation of ≈0.35 arcsec (≈13 au) from the host star. We revise our orbital fit by incorporating our dual-epoch relative astrometry using the open-source Markov chain Monte Carlo orbit fitting code orvara. We obtain a dynamical mass of $65.9_{-1.7}^{+2.0} M_{\rm Jup}$ that places HD 176535 B firmly in the brown dwarf regime. HD 176535 B is a new benchmark dwarf useful for constraining the evolutionary and atmospheric models of high-mass brown dwarfs. We found a luminosity of $\rm log(\mathit{ L}_{bol}/L_{\odot }) = -5.26\pm 0.07$ and a model-dependent effective temperature of 980 ± 35 K for HD 176535 B. We infer HD 176535 B to be a T dwarf from its mass, age, and luminosity. Our dynamical mass suggests that some substellar evolutionary models may be underestimating luminosity for high-mass T dwarfs. Given its angular separation and luminosity, HD 176535 B would make a promising candidate for Aperture Masking Interferometry with JWST and GRAVITY/Keck Planet Imager and Characterizer, and further spectroscopic characterization with instruments like the CHARIS/SCExAO/Subaru integral field spectrograph.