with the newest extreme adaptive optics systems can reveal substellar or even planet-mass companions on solar system-like scales at an increased frequency compared to blind surveys.
Facility and now extreme adaptive optics systems have provided the first direct detections of exoplanets around nearby, young stars (Marois et al. 2008 (Marois et al. , 2010 Kuzuhara et al. 2013; Currie et al. 2014a Currie et al. , 2015 Macintosh et al. 2015; Keppler et al. 2018) . Blind direct imaging surveys have provided the first constraints on the frequency of jovian planets and low-mass brown dwarfs at orbital distances of 10-500 au (Brandt et al. 2014; Nielsen et al. 2019) .
Unfortunately, the low yield of blind direct imaging surveys shows that exoplanets detectable with current instruments are rare (e.g. Nielsen et al. 2019) , especially around solar and subsolar-mass stars. The few companions found by these surveys typically exceed five jovian masses and orbit well beyond separations where the jovian planet frequency peaks (Fernandes et al. 2019) . Improving blind survey yields will only be possible with substantial contrast gains at small separations and/or sensitivity gains at wider separations (Crepp & Johnson 2011) . Limited sample sizes and sparse coverage of ages, temperatures, and surface gravities impede our understanding of the atmospheric evolution of gas giant planets.
Targeted searches focused on stars showing evidence of gravitational pulls from massive planets could improve survey yields. While radial-velocity data have been used for sample selection with some success for finding (sub-)stellar companions (Crepp et al. 2014) , astrometry can also select direct imaging targets, even around stars whose activity and spectral type preclude precise radial-velocity measurements. The Hipparcos-Gaia Catalog of Accelerations (HGCA; Brandt 2018) provides absolute astrometry for 115,000 nearby stars, including those with clear dynamical evidence for unseen massive companions. HGCA-derived accelerations can provide dynamical masses of imaged exoplanets and low-mass brown dwarfs independent of luminosity evolution models and irrespective of stellar age uncertainties Dupuy et al. 2019) .
In this Letter, we report the direct imaging discovery of a low-mass substellar companion 20 au from the nearby Sun-like star HD 33632 Aa using the Subaru Coronagraphic Extreme Adaptive Optics Project (SCExAO, Jovanovic et al. 2015b) coupled to the CHARIS integral field spectrograph (Groff et al. 2016) and using the NIRC2 camera on Keck. The HGCA identifies an astrometric acceleration for the primary induced by HD 33632 Ab, providing a dynamical mass of 46 ± 8 M J for the companion. CHARIS spectra and NIRC2 photometry reveal HD 33632 Ab to be a highergravity, higher-mass counterpart to young L/T transition planet-mass objects such as HR 8799 cde with a low eccentricity more typical of directly-imaged planets than brown dwarfs.
The HD 33632 primary is a 26.56 pc distant F8V 1.1 M main sequence star (Takeda et al. 2007; Gaia Collaboration et al. 2018 Brandt et al. (2014) favor the low end of this range (∼ 1-2.5 Gyr). HD 33632 Aa's fractional X-ray luminosity (log(L X /L bol = −5.2) overlaps with those for 800 Myr-old Hyades stars (Brandt & Huang 2015; Freund et al. 2020) . Li abundances are consistent with the Hyades but marginally consistent with 4 Gyrold M67 (Ramírez et al. 2012; Castro et al. 2016) . Ages estimated using neutron capture elements abundances individually (e.g. [Zr/Fe]) or in aggregate (e.g. the average abundances of "heavy" elements Ce, Ba, and La; [Y/Mg]) cluster around 1.5-2.5 Gyr (Spina et al. 2018) . Given activity, gyrochronology, and neutron capture abundance-derived estimates, we adopt a broad age range of 1.5 +3.0 −0.7 Gyr with a preference for 1.0-2.5 Gyr. HGCA reveals the HD 33632 primary to have a statistically significant acceleration (χ 2 ≈ 116) for a model of constant proper motion, equivalent to 10.5 σ with two degrees of freedom: evidence for the dynamical pull of a companion. The acceleration is almost 1000 times larger than that expected from the M-dwarf companion. The Lick radial-velocity survey monitored HD 33632 Aa for eleven years (18 Jan 1998 to 1 Feb 2009) without finding evidence of a planet or brown dwarf (Fischer et al. 2014) : the companion must be widely separated. Following es- Note-a) For CHARIS data, this column refers to the wavelength range. For NIRC2 imaging data, it refers to the central wavelength. b) Twenty six frames (∼ 20 minutes) of the 31 August 2020 data were also obtained with the star dithered by 1. 1 to the north and east. c) From the Canada France Hawaii Telescope seeing monitor. d) The SNR ratio of HD 33632 Ab in the wavelength-collapsed image reduced using "conservative" ADI/A-LOCI algorithm settings or SDI/A-LOCI for the case of 2018 October SCExAO/CHARIS data. Figure 1 . HD 33632 Ab detections in 2020 (top panels) and 2018 (bottom panels) reduced using A-LOCI with different combinations of ADI, SDI, and ADI+SDI. Our adopted rotation gap/radial scaling criteria precludes access to masked regions.
timates in Brandt et al. (2019) , its minimum mass at angles accessible by SCExAO/CHARIS (ρ ∼ 0. 1-1. 05) ranges from a few jovian masses to over 50 M J . Table 1 summarizes our observations with SCExAO/CHARIS and Keck/NIRC2 in 2018 and 2020. The seeing ranged between θ V = 0. 4 and 1. 0; conditions were photometric each night. SCExAO achieved a highfidelity AO correction; we utilized the CHARIS integral field spectrograph (Groff et al. 2016 ) in low spectral resolution (broadband) mode covering the JHK passbands simultaneously (1.16-2.37 µm at R ∼ 18). NIRC2 data were taken in the L p broadband filter (λ o = 3.78 µm) using Keck's facility AO system.
All CHARIS data utilized satellite spots for precise astrometric and spectrophotometric calibration (Jovanovic et al. 2015a ). The CHARIS data were taken with a Lyot coronagraph (0. 139 radius occulting mask); the NIRC2 data were taken behind a Lyot coronagraph with a 0. 2 radius occulting mask. All observations were conducted in pupil tracking mode, enabling angular differential imaging (ADI; Marois et al. 2006) ; the CHARIS data also enable spectral differential imaging (SDI; Marois et al. 2000) .
We extracted CHARIS data cubes using the standard CHARIS pipeline (Brandt et al. 2017 ) and perfomed basic reduction steps including sky subtraction, image registration, and spectrophotometric calibration as in Currie et al. (2018) . For spectrophotometric calibration, we adopted a Kurucz stellar atmosphere model appropriate for an F8V star. We reduced the NIRC2 data using the ADI-based pipeline from Currie et al. (2011) .
For all data, we used A-LOCI for point-spread function (PSF) subtraction (Currie et al. 2018 ). Due to a lack of field rotation, we subtracted the PSF for the 2018 CHARIS data in SDI-mode only. The 2020 CHARIS data enabled PSF subtraction with ADI, SDI, and ADI+SDI (i.e. performing SDI on the A-LOCI PSF subtraction residuals produced from ADI). In all cases, we used a large rotation gap/radial scaling criterion of δ > 1 λ/D to limit self-subtraction. For the CHARIS SDI reduction (component), we employed a pixel mask over the subtraction zone (Marois et al. 2010; Currie et al. 2015) to significantly reduce biasing of a companion spectrum. Figure 1 shows the detection of a faint point source, hereafter HD 33632 Ab, ρ ≈ 0. 75 west from the primary in each epoch. The detection's signal-to-noise ratios (SNR) in conservative ADI/A-LOCI or SDI/A-LOCI reductions range between 11-22 in the 2018 discovery epoch to 41 for the 31 August 2020 CHARIS broadband follow-up data. CHARIS data do not identify any other companion within ρ ∼ 1. 05. In addition to resolving HD 33632 Ab, NIRC2 data identify a wider separation object at [E,N] = [−1. 22, −1. 85] and L p ∼ 16.3 (not shown). However, 2020 CHARIS data reveal this object to be a background star based on its astrometry and its near-infrared brightness. To derive throughput-corrected spectrophotometry, we forward-modeled point sources at HD 33632 Ab's locations using stored LOCI coefficients (Currie et al. 2018) . We focus on ADI-only reductions due to the more straightforward applicability of ADI-reduced data to forward-modeling (Pueyo 2016) . We adopt our highestquality data set (31 August 2020) for our analyses. Our results agree with those from other data sets and with SDI reductions.
We extract spectra from the best-fit centroid position of the wavelength-collapsed image. Throughputs typically exceed 75%; because we use ADI only, they are in-dependent of HD 33632 Ab's spectrum. HD 33632 Ab's broadband photometry in standard Mauna Kea Observatory filters is J = 16.91 ± 0.11, H = 16.00 ± 0.09, K s = 15.37 ± 0.09, and L p = 13.67 ± 0.15 (∆J, H, K s = 11.52, 10.83, 10.21). Photometry obtained at different epochs or with different reductions agrees with these results to within errors.
HD 33632 Ab lies at the transition from cloudy L dwarfs to (nearly) cloud free T dwarfs (Burrows et al. 2006) . Figure 2 places it on an infrared color-magnitude diagram. The companion's J/J-K s position is consistent with field L/T objects, while its L p /H-L p position agrees with both field objects and the low gravity, cloudy exoplanets HR 8799 cde (Currie et al. 2011; Barman et al. 2015) .
The CHARIS spectrum of HD 33632 Ab (Figure 2 , bottom-left panel; Table 2 ) is highly peaked, characteristic of a substellar atmosphere shaped by water opacity. The 2018 CHARIS spectrum is significantly noisier and offset by ∼ 10% compared to the 2020 epoch data. Otherwise, spectra extracted using different reduction techniques and taken on different epochs agree to within errors. The detection significance of HD 33632 Ab reaches SNR = 20-30 in some channels for the 31 August 2020/ADI-reduced spectrum. However, contributions from uncertainties in the absolute flux calibration (i.e. the satellite spot flux uncertainty) limit the spectrophotometric precision to about 10%.
We compared HD 33632 Ab's CHARIS spectrum to objects in the Montreal Spectral Library 2 (e.g. Gagné et al. 2014) , considering the impact of spatially and spectrally correlated noise (Greco & Brandt 2016 ). The spectral covariance at HD 33632 Ab's angular separation in the 31 Aug 2020 ADI/A-LOCI reduced data is nearly zero for off-diagonal elements indicating largely uncorrelated noise.
The L9.5 field brown dwarfs SIMPJ0956-1447 and SIMPJ0150+3827 provide the best fit to HD 33632 Ab's spectrum (Figure 2, bottom-right panel) . While the Montreal library contains few objects later than L5, the χ 2 distribution for library objects exhibits a local mininum between L6.5 and T0.5. Thus, conservatively we assign HD 33632 Ab a spectral type of L9.5 +1.0 −3.0 . Adopting the mapping between spectral type and temperature from Stephens et al. (2009) , the object's effective temperature is 1300 +100 −100 K: slighty higher than that derived for the inner three HR 8799 planets (Currie et al. 2011; Barman et al. 2015) . Adopting the polynomial fit from Note-Throughput-corrected HD 33632 Ab spectrum extracted from 31 August 2020 data, reduced using ADI/A-LOCI. Golimowski et al. (2004) and a distance of 26.56 pc, HD 33632 Ab's luminosity is log 10 (L/L ) = −4.62 +0.04 −0.08 .
As our forward-modeling revealed negligible astrometric bias and different reductions yielded consistent astrometry, we used gaussian centroiding to determine the position of HD 33632 Ab. HD 33632 Ab is lo- 
We use the open-source code orvara (Brandt et al., submitted) to fit HD 33632 Ab's mass and orbit using a combination of HGCA measurements, Lick Observatory radial-velocity data (Fischer et al. 2014) , and relative astrometry from CHARIS and NIRC2. No other bound companion is detected with CHARIS, Lick data do not identify an inner companion, and the only other object seen with NIRC2 is a background star. Given its separation, the M dwarf companion contributes <1% of the acceleration seen by the HGCA. Therefore, we assume that HD 33632 Ab is solely responsible for the HGCA acceleration. The astrometric acceleration of ∼0.1 mas yr −2 is more than 1000 times smaller than the apparent acceleration due to parallactic motion: it negligibly affects the Gaia parallax. We fit for a radial velocity jitter, the six Keplerian orbital elements, and the masses of both components. We adopt a prior on HD 33632 Aa's mass of 1.1 ± 0.1 M and a log-uniform prior on the mass of HD 33632 Ab. We analytically marginalize out the nuisance parameters of barycenter proper motion, parallax, and radial velocity zero point. Figure 3 shows our posterior distributions for selected orbital parameters, the primary mass (almost identical to our prior of 1.1 ± 0.1 M ), and HD 33632 Ab's mass. Table 3 lists all priors and posteriors from the fit. The companion has a best-fit semimajor axis of 21.2 +4.1 −4.5 au with an orbit inclined by i = 39. • 4 +8.0 −20 . The semimajor axis is double-peaked because our relative astrometry spans only a 2 year baseline, while 25 years elapse between absolute astrometric measurements: Hipparcos mainly constrains HD 33632 Ab's 1991 position. Continued orbital monitoring will remove the double peak, improving the companion's mass measurement. HD 33632 Ab's median posterior eccentricity is e = 0.19 but the distribution peaks at zero (i.e. a circular orbit): the 68% and 95% upper limits are 0.29 and 0.46, respectively.
HD 33632 Ab has a dynamical mass of 46.4 +8.1 −7.6 M J and a mass ratio of q ∼ 4.0% ± 0.7%. We may compare this to masses derived by matching its observed luminosity and age with substellar evolutionary models. For its full (preferred) age range of 1.5 +3.0 −0.7 Gyr (1-2.5 Gyr), the Baraffe et al. (2003) and Burrows et al. (1997) luminosity evolution models imply masses of 48 +17 The orbit fits used Hipparcos and Gaia absolute astrometry, Lick Observatory precision radial-velocity measurements for HD 33632 Aa, and relative astrometry of HD 33632 Ab. The posterior on HD 33632 Aa's mass is dominated by our adopted prior of 1.1 ± 0.1 M . The insets display the best-fit orbit along with 100 random orbits drawn from our MCMC in both relative astrometry and radial velocity, color-coded by HD 33632 Ab's mass. Zero radial velocity corresponds to the system barycenter of the best-fit orbit (black curve); all radial velocity data points and curves are shifted by the best-fit radial velocity offsets to align with the black curve. model-dependent masses agree well with the dynamical mass if we adopt our preferred age of 1-2.5 Gyr, but show moderate tension for ages 2.5 Gyr.
Our SCExAO/CHARIS and complementary Keck/NIRC2 imaging identify a low-mass ratio brown dwarf companion at r proj ∼20 au around the main sequence Sun-like star, HD 33632 Aa. HD 33632 Ab's infrared colors and spectrum are best matched by a field L/T transition object. Combining measurements from Hipparcos, Gaia, and precision radial-velocity with our CHARIS/NIRC2 relative astrometry yields a dynamical mass of 46.4 +8.1 −7.6 M J . This is more precise than estimates from luminosity evolution models alone, given age uncertainties, and it will improve with continued astrometric monitoring and future Gaia data releases. Provided that the system age and companion mass are better constrained in future work, HD 33632 Ab adds to the still-small sample of objects enabling us to test substellar evolutionary models across masses and ages.
HD 33632 Ab lies at the field L/T transition tracking the dissipation of clouds/dust in substellar atmospheres. In L p /H-L p , its colors overlap with and its temperature is just slightly exceeds those of the young exoplanets HR 8799 cde. HD 33632 Ab represents an older, higher mass, higher gravity, and (likely) less cloudy/dusty counterpart to some of the first imaged exoplanets. Thermal infrared follow-up imaging and spectroscopy on the ground and with the James Webb Space Telescope could further clarify how its atmospheric chemistry compares with HR 8799 cde and other companions. Theoretical modeling incorporating these higher-resolution data will quantify constraints on HD 33632 Ab's atmosphere.
HD 33632 Ab's mass, mass ratio, and orbital separation are nearly identical to those of the older, colder brown dwarf companion GJ 758 B (Thalmann et al. 2009; Brandt et al. 2019) , which likewise orbits a Sunlike star. Its mass and mass ratio bracket those of other companions to mid F to early G stars orbiting at 20-50 au: e.g. the more massive, older T dwarf companion HD 19467 B (Crepp et al. 2014 ) and (near-) planet mass companions HR 2562 b (Konopacky et al. 2016) and GJ 504 b (Kuzuhara et al. 2013 ) with comparable or younger ages.
HD 33632 Ab's eccentricity posterior probability distribution has a 68% upper limit of 0.29. This is smaller than the eccentricities of most brown dwarfs studied in Bowler et al. (2020) , who argue that systematically higher eccentricities of brown dwarf-mass companions relative to exoplanets point to a different formation mechanism. Future CHARIS/NIRC2 astrometry for HD 33632 Ab together with more precise HD 33632 Aa astrometry in future Gaia data releases will better constrain the companion's mass and eccentricity.
In addition to HD 33632 Aa, hundreds of nearby stars show evidence for a statistically significant acceleration plausibly caused by an unidentified substellar companion. Upcoming Gaia data releases will identify thousands more accelerating stars. Some will be young enough that extreme AO systems like SCExAO can image self-luminous jovian exoplanets responsible for the astrometric trends. Our work is a proof-in-concept that direct imaging searches targeting HGCA-selected stars may identify new substellar and potentially even planetary companions, increasing the yield of discoveries.
We thank Eric Mamajek and William Cochran for helpful comments regarding the HD 33632 system properties.
The authors acknowledge the very significant cultural role and reverence that the summit of Mauna Kea holds within the Hawaiian community. We are most fortunate to have the opportunity to conduct observations from this mountain.
We acknowledge the critical importance of the current and recent Subaru and Keck Observatory daycrew, technicians, support astronomers, telescope operators, computer support, and office staff employees, especially during the challenging times presented by the COVID-19 pandemic. Their expertise, ingenuity, and dedication is indispensable to the continued successful operation of these observatories.
We thank the Subaru and NASA Keck Time Allocation Committees for their generous support of this program. TC was supported by a NASA Senior Postdoctoral Fellowship and NASA/Keck grant LK-2663-948181. TB gratefully acknowledges support from the Heising-Simons foundation and from NASA under grant #80NSSC18K0439. MT is supported by JSPS KAK-ENHI Grant # 18H05442.
The development of SCExAO was supported by JSPS (Grant-in-Aid for Research #23340051, #26220704 & #23103002), Astrobiology Center of NINS, Japan, the Mt Cuba Foundation, and the director's contingency fund at Subaru Telescope. CHARIS was developed under the support by the Grant-in-Aid for Scientific Research on Innovative Areas #2302. Some of the data presented herein were obtained at the W. M. Keck Observatory, which is operated as a scientific partnership among the California Institute of Technology, the University of California and the National Aeronautics and Space Administration. The Observatory was made possible by the generous financial support of the W. M.
Keck Foundation.