Search

As exoplanet surveys reach ever-higher sensitivities and durations, planets analogous to the solar system giant planets are increasingly within reach. HD 28185 is a Sun-like star known to host a $m\sin i=6 M_J$ planet on an Earth-like orbit; more recently, a brown dwarf with a more distant orbit has been claimed. In this work we present a comprehensive reanalysis of the HD 28185 system, based on 22 years of radial velocity observations and precision Hipparcos-Gaia astrometry. We confirm the previous characterisation of HD 28185 b as a temperate giant planet, with its $385.92^{+0.06}_{-0.07}$ day orbital period giving it an Earth-like incident flux. In contrast, we substantially revise the parameters of HD 28185 c; with a new mass of $m=6.0\pm0.6 M_J$ we reclassify this companion as a super-jovian planet. HD 28185 c has an orbital period of $24.9^{+1.3}_{-1.1}$ years, a semi-major axis of $8.50^{+0.29}_{-0.26}$ AU, and a modest eccentricity of $0.15\pm0.04$, resulting in one of the most Saturn-like orbits of any known exoplanet. HD 28185 c lies at the current intersection of detection limits for RVs and direct imaging, and highlights how the discovery of giant planets at $\approx$10 AU separations is becoming increasingly routine., Comment: 17 pages, 8 figures, 7 tables. Published in MNRAS

We analyze near-infrared integral field spectropolarimetry of the AB Aurigae protoplanetary disk and protoplanet (AB Aur b), obtained with SCExAO/CHARIS in 22 wavelength channels covering the J, H, and K passbands ($\lambda_{\rm o}$ = 1.1--2.4 $\mu m$) over angular separations of $\rho$ $\approx$ 0.13" to 1.1" ($\sim$20--175 au). Our images resolve spiral structures in the disk in each CHARIS channel. At the longest wavelengths, the data may reveal an extension of the western spiral seen in previous polarimetric data at $\rho$ $<$ 0.3" out to larger distances clockwise from the protoplanet AB Aur b, coincident with the ALMA-detected $CO$ gas spiral. While AB Aur b is detectable in complementary total intensity data, it is a non-detection in polarized light at $\lambda$ $>$ 1.3 $\mu $m. While the observed disk color is extremely red across $JHK$, the disk has a blue intrinsic scattering color consistent with small dust grains. The disk's polarization spectrum is redder than AB Aur b's total intensity spectrum. The polarization fraction peaks at $\sim$ 0.6 along the major disk axis. Radiative transfer modeling of the CHARIS data shows that small, porous dust grains with a porosity of $p$ = 0.6--0.8 better reproduce the scattered-light appearance of the disk than more compact spheres ($p$ = 0.3), especially the polarization fraction. This work demonstrates the utility of integral field spectropolarimetry to characterize structures in protoplanetary disks and elucidate the properties of the disks' dust., Comment: 23 pages, 17 Figures, accepted to AAS Journals

We present Subaru/SCExAO + CHARIS broadband (JHK) integral field spectroscopy of HD 34700 A in polarized light. CHARIS has the unique ability to obtain polarized integral field images at 22 wavelength channels in broadband, as the incoming light is first split into different polarization states before passing though the lenslet array. We recover the transition disk around HD 34700 A in multiband polarized light in our data. We combine our polarized intensity data with previous total intensity data to examine the scattering profiles, scattering phase functions and polarized fraction of the disk at multiple wavelengths. We also carry out 3D Monte Carlo radiative transfer simulations of the disk using MCFOST, and make qualitative comparisons between our models and data to constrain dust grain properties. We find that in addition to micron-sized dust grains, a population of sub-micron grains is needed to match the surface brightness in polarized light and polarized fraction. This could indicate the existence of a population of small grains in the disk, or it could be caused by Mie theory simulations using additional small grains to compensate for sub-micron structures of real dust aggregates. We find models that match the polarized fraction of the data but the models do not apply strong constraints on the dust grain type or compositions. We find no models that can match all observed properties of the disk. More detailed modeling using realistic dust aggregates with irregular surfaces and complex structures is required to further constrain the dust properties.

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 Hipparcos and 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} M_{jup}$, an eccentricity of ${0.260}_{-0.059}^{+0.065}$, and a nearly face-on inclination of $174.81_{-0.50}^{+0.48}$ degrees. For HD 63754 B, we obtain an L' band absolute magnitude of $L' = 11.39\pm0.06$ mag, from which we infer a bolometric luminosity of $log(L_{bol}/L_{\odot})= -4.55 \pm0.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.

This tutorial covers the use of absolute astrometry, in particular from the combination of the Hipparcos and Gaia missions, to identify faint companions to nearby stars and to measure the masses and orbits of those companions. Absolute astrometry has been used with increasing success to discover new planets and brown dwarfs and to measure masses and orbits for systems with periods as long as centuries. This tutorial summarizes the nature of the underlying astrometric data, the approach typically used to fit orbits, and the assumptions about that data implicit throughout the process. It attempts to provide intuition for the sensitivity of astrometry as a function of stellar and companion properties and how the available constraints depend on the character and quantity of data available. This tutorial is written for someone with some background in astronomy but with no more than a minimal acquaintance with astrometry or orbit fitting., Comment: 17 pages, 5 figures, PASP accepted. Tutorial intended primarily as an introduction and reference for non-experts

With a dynamical mass of $3 \, M_\mathrm{Jup}$, the recently discovered giant planet AF Lep b is the lowest-mass imaged planet with a direct mass measurement. Its youth and spectral type near the L/T transition make it a promising target to study the impact of clouds and atmospheric chemistry at low surface gravities. In this work, we present JWST/NIRCam imaging of AF Lep b. Across two epochs, we detect AF Lep b in F444W ($4.4 \, \mathrm{\mu m}$) with S/N ratios of $9.6$ and $8.7$, respectively. At the planet's separation of $320 \, \mathrm{mas}$ during the observations, the coronagraphic throughput is ${\approx}7\%$, demonstrating that NIRCam's excellent sensitivity persists down to small separations. The F444W photometry of AF Lep b affirms the presence of disequilibrium carbon chemistry and enhanced atmospheric metallicity. These observations also place deep limits on wider-separation planets in the system, ruling out $1.1 \, M_\mathrm{Jup}$ planets beyond $15.6 \, \mathrm{au}$ ($0.58$ arcsec), $1.1 \, M_\mathrm{Sat}$ planets beyond $27 \, \mathrm{au}$ ($1$ arcsec), and $2.8 \, M_\mathrm{Nep}$ planets beyond $67 \, \mathrm{au}$ ($2.5$ arcsec). We also present new Keck/NIRC2 $L'$ imaging of AF Lep b; combining this with the two epochs of F444W photometry and previous Keck $L'$ photometry provides limits on the long-term $3{-}5 \, \mathrm{\mu m}$ variability of AF Lep b on months-to-years timescales. AF Lep b is the closest-separation planet imaged with JWST to date, demonstrating that planets can be recovered well inside the nominal (50\% throughput) NIRCam coronagraph inner working angle., Comment: 17 pages, 4 figures, accepted in ApJL

Measuring fundamental stellar parameters is key to fully comprehending the evolution of stars. However, current theoretical models over-predict effective temperatures, and under-predict radii, compared to observations of K and M dwarfs (radius inflation problem). In this work, we developed a model independent method to infer precise radii of single FGK and M dwarfs using Gaia DR3 parallaxes and photometry, and we used it to study the radius inflation problem. We calibrated nine surface brightness-color relations for the three Gaia magnitudes and colors using a sample of stars with angular diameter measurements. We achieved an accuracy of 4% in our angular diameter estimations, which Gaia's parallaxes allow us to convert to a physical radii. We validated our method by comparing our radius measurements with literature samples and the Gaia DR3 catalog, which confirmed the accuracy of our method and revealed systematic offsets in the Gaia measurements. Moreover, we used a sample with measured Halpha equivalent width (HaEW), a magnetic activity indicator, to study the radius inflation problem. We demonstrated that active stars have larger radii than inactive stars, showing that radius inflation is correlated with magnetic activity. We found a correlation between the radius inflation of active stars and HaEW for the mass bin 0.5

The two planets of the HAT-P-11 system represent fascinating dynamical puzzles due to their significant eccentricities and orbital misalignments. In particular, HAT-P-11 b is on a close-in orbit that tides should have circularized well within the age of the system. Here we propose a two-step dynamical process that can reproduce all intriguing aspects of the system. We first invoke planet-planet scattering to generate significant eccentricities and mutual inclinations between the planets. We then propose that this misalignment initiated von-Zeipel-Lidov-Kozai cycles and high-eccentricity migration that ultimately brought HAT-P-11 b to its present-day orbit. We find that this scenario is fully consistent only when significant tidally-driven radius inflation is accounted for during the tidal migration. We present a suite of N-body simulations exploring each phase of evolution and show that this scenario is consistent with all observational posteriors and the reported age of the system., Comment: 17 pages, 8 figures, submitted to ApJ

Planet-star obliquity and planet-planet ]mutual inclination encode a planetary system's dynamical history, but both of their values are hard to measure for misaligned systems with close-in companions. HAT-P-11 is a K4 star with two known planets: a close-in, misaligned super-Neptune with a $\approx$5-day orbit, and an outer super-Jupiter with a $\approx$10-year orbit. In this work we present a joint orbit fit of the HAT-P-11 system with astrometry and S-index corrected RV data. By combining our results with previous constraints on the orientation of the star and the inner planet, we find that all three angular momenta -- those of the star, planet b, and planet c -- are significantly misaligned. We confirm the status of planet c as a super-Jupiter, with $2.68\pm0.41\, \mathrm{M_{\rm Jup}}$ at a semimajor axis of $4.10\pm0.06\,$AU, and planet b's mass of $\mathrm{M_b\sin{i_b}}=0.074\pm0.004\, \mathrm{M_{\rm Jup}}$. We present the posterior probability distribution of obliquity between star A and planet c, and mutual inclination between planet b and planet c., Comment: 14 pages, 10 figures, submitted to AJ

This paper implements likelihood-based jump detection for detectors read out up-the-ramp, using the entire set of reads to compute likelihoods. The approach compares the $\chi^2$ value of a fit with and without a jump for every possible jump location. I show that this approach can be substantially more sensitive than one that only uses the difference between sequential groups of reads, especially for long ramps and for jumps that occur in the middle of a group of reads. It can also be implemented for a computational cost that is linear in the number of resultants. I provide and describe a pure Python implementation that can process a 10-resultant ramp on a $4096 \times 4096$ detector in $\approx$20 seconds, including iterative cosmic ray detection and removal, on a single core of a 2020 Macbook Air. This Python implementation, together with tests and a tutorial notebook, are available at \url{https://github.com/t-brandt/fitramp}. I also provide tests and demonstrations of the full ramp fitting and cosmic ray rejection approach on data from JWST., Comment: 17 pages, 10 figures, PASP accepted. Code available at https://github.com/t-brandt/fitramp. This submission supersedes Sections 6 and 7 and Appendix C of arXiv:2309.08753v1. arXiv admin note: substantial text overlap with arXiv:2309.08753

We present the direct-imaging discovery of a substellar companion (a massive planet or low-mass brown dwarf) to the young, $\gamma$ Doradus-type variable star, HIP 39017 (HD 65526). The companion's SCExAO/CHARIS JHK ($1.1-2.4\mu$m) spectrum and Keck/NIRC2 L$^{\prime}$ photometry indicate that it is an L/T transition object. A comparison of the JHK+L$^{\prime}$ 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 semi-major axis of $23.8^{+8.7}_{-6.1}$ au, a dynamical companion mass of $30^{+31}_{-12}$~M$_J$, and a mass ratio of $\sim$1.9\%, properties most consistent with low-mass brown dwarfs. However, its mass estimated from luminosity models is a lower $\sim$13.8 $M_{\rm J}$ due to an estimated young age ($\lesssim$ 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 $\sim$1.4\%. Analysis of the host star's multi-frequency $\gamma$ Dor-type pulsations, astrometric monitoring of HIP 39017b, 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., Comment: Published in AJ: April 9, 2024

Multi-planet systems exhibit a diversity of architectures that diverge from the solar system and contribute to the topic of exoplanet demographics. Radial velocity (RV) surveys form a crucial component of exoplanet surveys, as their long observational baselines allow searches for more distant planetary orbits. This work provides a significantly revised architecture for the multi-planet system HD 134606 using both HARPS and UCLES RVs. We confirm the presence of previously reported planets b, c, and d with periods $12.0897^{+0.0019}_{-0.0018}$, $58.947^{+0.056}_{-0.054}$, and $958.7^{+6.3}_{-5.9}$ days, and masses $9.14^{+0.65}_{-0.63}$, $11.0\pm1$, and $44.5\pm2.9$ Earth masses respectively, with the planet d orbit significantly revised to over double that originally reported. We report two newly detected super-Earths, e and f, with periods $4.31943^{+0.00075}_{-0.00068}$ and $26.9^{+0.019}_{-0.017}$ days, and masses $2.31^{+0.36}_{-0.35}$ and $5.52^{+0.74}_{-0.73}$ Earth masses, respectively. In addition, we identify a linear trend in the RV time series, and the cause of this acceleration is deemed to be a newly detected sub-stellar companion at large separation. HD 134606 now displays four low mass planets in a compact region near the star, one gas giant further out in the Habitable Zone, an additional massive companion in the outer regime, and a low mass M dwarf stellar companion at large separation, making it an intriguing target for system formation/evolution studies. The location of planet d in the Habitable Zone proves to be an exciting candidate for future space-based direct imaging missions, whereas continued RV observations of this system are recommended for understanding the nature of the massive, long period companion., Comment: 23 pages, 10 figures, 4 tables, accepted for publication in the Astronomical Journal

We present the discoveries of two giant planets orbiting the red giant branch (RGB) star HD 112570 and the red clump (RC) star HD 154391, based on the radial velocity (RV) measurements from Xinglong station and Okayama Astrophysical Observatory (OAO). Spectroscopic and asteroseismic analyses suggest that HD 112570 has a mass of $1.15\pm0.12\,M_{\odot}$, a radius of $9.85\pm0.23\,R_{\odot}$, a metallicity [Fe/H] of $-0.46\pm0.1$ and a ${\rm log}\,g$ of $2.47\pm0.1$. With the joint analysis of RV and Hipparcos-Gaia astrometry, we obtain a dynamical mass of $M_{\rm p}={3.42}_{-0.84}^{+1.4}\ M_{\rm Jup}$, a period of $P={2615}_{-77}^{+85}$ days and a moderate eccentricity of $e={0.20}_{-0.14}^{+0.16}$ for the Jovian planet HD 112570 b. For HD 154391, it has a mass of $2.07\pm0.03\,M_{\odot}$, a radius of $8.56\pm0.05\,R_{\odot}$, a metallicity [Fe/H] of $0.07\pm0.1$ and a ${\rm log}\,g$ of $2.86\pm0.1$. The super-Jupiter HD 154391 b has a mass of $M_{\rm p}={9.1}_{-1.9}^{+2.8}\ M_{\rm Jup}$, a period of $P={5163}_{-57}^{+60}$ days and an eccentricity of $e={0.20}_{-0.04}^{+0.04}$. We found HD 154391 b has one of the longest orbital period among those ever discovered orbiting evolved stars, which may provide a valuable case in our understanding of planetary formation at wider orbits. Moreover, while a mass gap at $4\,M_{\rm Jup}$ seems to be present in the population of giant stars, there appears to be no significant differences in the distribution of metallicity among giant planets with masses above or below this threshold. Finally, The origin of the abnormal accumulation near 2 au for planets around large evolved stars ($R_{\star}>21\,R_{\odot}$), remains unclear., Comment: 26 pages, 18 figures, 8 tables, Accepted for publication in The Astronomical Journal

The effect of stellar multiplicity on planetary architecture and orbital dynamics provides an important context for exoplanet demographics. We present a volume-limited catalog 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 indicate 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 motion anomalies show nearly four times more false positives due to Eclipsing Binaries compared to TOIs with marginal proper motion anomalies. In addition, we find tentative evidence that orbital periods of planets orbiting TOIs with significant proper motion anomalies are shorter than those orbiting TOIs without significant proper motion anomalies, consistent with the scenario that stellar companions can truncate planet-forming disks. Furthermore, TOIs with significant proper motion anomalies exhibit lower Gaia differential velocities in comparison to field stars with significant proper motion 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 deg between the orbit of BC around A and that of C around B. The coplanarity may explain why multiple planets around LTT 1445 A survive in the dynamically hostile environment of this system., Comment: 28 pages, 9 figures, published by AJ

This paper derives the optimal fit to a pixel's count rate in the case of an ideal detector read out nondestructively in the presence of both read and photon noise. The approach is general for any readout scheme, provides closed-form expressions for all quantities, and has a computational cost that is linear in the number of resultants (groups of reads). I also derive the bias of the fit from estimating the covariance matrix and show how to remove it to first order. The ramp-fitting algorithm I describe provides the $\chi^2$ value of the fit of a line to the accumulated counts, which can be interpreted as a goodness-of-fit metric. I provide and describe a pure Python implementation of these algorithms that can process a 10-resultant ramp on a $4096 \times 4096$ detector in $\approx$8 seconds with bias removal on a single core of a 2020 Macbook Air. This Python implementation, together with tests and a tutorial notebook, are available at \url{https://github.com/t-brandt/fitramp}. A companion paper describes a jump detection algorithm based on hypothesis testing of ramp fits and demonstrates all algorithms on data from JWST., Comment: 27 pages, 10 figures, PASP accepted. Code available at https://github.com/t-brandt/fitramp. Sections 6, 7 and Appendix C of v1 moved to arxiv:2404.01326

We present the pyKLIP-CHARIS post-processing pipeline, a Python library that reduces high contrast imaging data for the CHARIS integral field spectrograph used with the SCExAO project on the Subaru Telescope. The pipeline is a part of the pyKLIP package, a Python library dedicated to the reduction of direct imaging data of exoplanets, brown dwarfs, and discs. For PSF subtraction, the pyKLIP-CHARIS post-processing pipeline relies on the core algorithms implemented in pyKLIP but uses image registration and calibrations that are unique to CHARIS. We describe the pipeline procedures, calibration results, and capabilities in processing imaging data acquired via the angular differential imaging and spectral differential imaging observing techniques. We showcase its performance on extracting spectra of injected synthetic point sources as well as compare the extracted spectra from real data sets on HD 33632 and HR 8799 to results in the literature. The pipeline is a python-based complement to the SCExAO project supported, widely used (and currently IDL-based) CHARIS data post-processing pipeline (CHARIS DPP) and provides an additional approach to reducing CHARIS data and extracting calibrated planet spectra., Comment: 17 pages, 13 figures

In this study, we revisit 32 planetary systems around evolved stars observed within the framework of the Okayama Planet Search Program and its collaborative framework of the EAPS-Net to search for additional companions and investigate the properties of stars and giant planets in multiple-planet systems. With our latest radial velocities obtained from Okayama Astrophysical Observatory (OAO), we confirm an additional giant planet in the wide orbit of 75 Cet system ($P_{\rm{c}} = 2051.62_{-40.47}^{+45.98}\ \rm{d}$, $M_{\rm{c}}\sin i=0.912_{-0.090}^{+0.088}\ M_{\rm{J}}$, and $a_{\rm{c}}=3.929_{-0.058}^{+0.052}\ \rm{au}$), along with five stars exhibiting long-term radial velocity accelerations, which indicates massive companions in the wide orbits. We have also found that the radial velocity variations of several planet-harboring stars may indicate additional planet candidates, stellar activities, or other understudied sources. These stars include $\epsilon$ Tau, 11 Com, 24 Boo, 41 Lyn, 14 And, HD 32518, and $\omega$ Ser. We further constrain the orbital configuration of the HD 5608, HD 14067, HD 120084, and HD 175679 systems by combining radial velocities with astrometry, as their host central stars exhibit significant astrometric accelerations. For other systems, we simply refine their orbital parameters. Moreover, our study indicates that the OPSP planet-harboring stars are more metal-poor compared to the currently known planet-harboring stars, and this is likely due to the $B-V$ color upper limit at 1.0 for star selection in the beginning of the survey. Finally, by investigating the less-massive giant planets ($< 5 M_{\rm{J}}$) around currently known planet-harboring evolved stars, we have found that metallicity positively correlates with the multiplicity and total planet mass of the system, which can be evidence for the core-accretion planet formation model., Comment: 49 figures, 4 tables, accepted by PASJ, RV data will be available online as supplementary after the publication

Being one of the most fundamental physical parameter of astronomical objects, mass plays a vital role in the study of exoplanets, including their temperature structure, chemical composition, formation, and evolution. However, nearly a quarter of the known confirmed exoplanets lack measurements of their masses. This is particularly severe for those discovered via the radial-velocity (RV) technique, which alone could only yield the minimum mass of planets. In this study, we use published RV data combined with astrometric data from a cross-calibrated Hipparcos-Gaia Catalog of Accelerations (HGCA) to jointly constrain the masses of 115 RV-detected substellar companions, by conducting full orbital fits using the public tool \texttt{orvara}. Among them, 9 exoplanets with $M_{\rm p}\,{\rm sin}\,i<13.5\ M_{\rm Jup}$ are reclassified to the brown dwarf (BD) regime, and 16 BD candidates ($13.5\leqslant M_{\rm p}\,{\rm sin}\,i<80\,M_{\rm Jup}$) turn out to be low-mass M dwarfs. We point out the presence of a transition in the BD regime as seen in the distributions of host star metallicity and orbital eccentricity with respect to planet masses. We confirm the previous findings that companions with masses below $42.5\ M_{\rm Jup}$ might primarily form in the protoplanetary disc through core accretion or disc gravitational instability, while those with masses above $42.5\ M_{\rm Jup}$ formed through the gravitational instability of molecular cloud like stars. Selection effects and detection biases which may affect our analysis to some extent, are discussed., Comment: 37 pages, 14 figures, accepted by Research in Astronomy and Astrophysics

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 HST data, and obtain the first dynamical mass of ${0.798}_{-0.041}^{+0.10}$ $M_{\odot}$, 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 \pm 0.011$ $M_{\odot}$) and HD 169889 B (${0.526}_{-0.037}^{+0.039}$ $M_{\odot}$), but their age precisions are limited by their uncertain temperatures. For HD 27786 B, the unusually small mass of $0.443 \pm 0.012$ $M_{\odot}$ suggests a history of rapid mass loss, possibly due to binary interaction in its progenitor's AGB phase. The orbits of HD 118475 and HD 136138 from our RV 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_{\odot}$ 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., Comment: 21 pages, 7 figures

We present the direct imaging discovery of a giant planet orbiting the young star AF Lep, a 1.2 $M_{\odot}$ member of the 24 $\pm$ 3 Myr $\beta$ Pic moving group. AF Lep was observed as part of our ongoing high-contrast imaging program targeting stars with astrometric accelerations between Hipparcos and Gaia that indicate the presence of substellar companions. Keck/NIRC2 observations in $L'$ with the Vector Vortex Coronagraph reveal a point source, AF Lep b, at ${\approx}340$ mas which exhibits orbital motion at the 6-$\sigma$ level over the course of 13 months. A joint orbit fit yields precise constraints on the planet's dynamical mass of 3.2$^{+0.7}_{-0.6}$ $M_\mathrm{Jup}$, semi-major axis of $8.4^{+1.1}_{-1.3}$ au, and eccentricity of $0.24^{+0.27}_{-0.15}$. AF Lep hosts a debris disk located at $\sim$50 au, but it is unlikely to be sculpted by AF Lep b, implying there may be additional planets in the system at wider separations. The stellar inclination ($i_* = 54^{+11}_{-9} {}^\circ$) and orbital inclination ($i_o = 50^{+9}_{-12} {}^\circ$) are in good agreement, which is consistent with the system having spin-orbit alignment. AF Lep b is the lowest-mass imaged planet with a dynamical mass measurement and highlights the promise of using astrometric accelerations as a tool to find and characterize long-period planets., Comment: 14 pages, 3 figures, accepted to ApJL

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}$ Gyrs at a distance of $36.99 \pm 0.03$ pc. In advance of our high-contrast imaging observations, we combined precision HARPS 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'_p = 9.20\pm0.06$ mag at a projected separation of $\approx$0.$\!\!''35$ ($\approx$13 AU) from the host star. We revise our orbital fit by incorporating our dual-epoch relative astrometry using the open-source MCMC orbit fitting code $\tt orvara$. 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(L_{bol}/L_{\odot}) = -5.26\pm0.07$ and a model-dependent effective temperature of $980 \pm 35$ K for HD 176535 B. 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/KPIC, and further spectroscopic characterization with instruments like the CHARIS/SCExAO/Subaru integral field spectrograph.

The Transiting Exoplanet Survey Satellite (TESS) is continuing its second extended mission after 55 sectors of observations. TESS publishes full-frame images (FFI) at a cadence of 1800, 600, or 200 seconds, allowing light curves to be extracted for stars beyond a limited number of pre-selected stars. Simulations show that thousands of exoplanets, eclipsing binaries, variable stars, and other astrophysical transients can be found in these FFI light curves. To obtain high-precision light curves, we forward model the FFI with the effective point spread function to remove contamination from nearby stars. We adopt star positions and magnitudes from Gaia DR3 as priors. The resulting light curves, called TESS-Gaia Light Curves (TGLC), show a photometric precision closely tracking the pre-launch prediction of the noise level. TGLC's photometric precision reaches <~2% at 16th TESS magnitude even in crowded fields. We publish TGLC Aperture and PSF light curves for stars down to 16th TESS magnitude through the Mikulski Archive for Space Telescopes (MAST) for all available sectors and will continue to deliver future light curves via DOI: 10.17909/610m-9474. The open-source package tglc is publicly available to enable any user to produce customized light curves., Comment: 23 pages, 13 figures, 4 tables, AJ accepted. Light curves are available at https://archive.stsci.edu/doi/resolve/resolve.html?doi=10.17909/610m-9474 tglc package is pip-installable and available at https://github.com/TeHanHunter/TESS_Gaia_Light_Curve

Direct imaging of gas giant exoplanets provides key information on planetary atmospheres and the architectures of planetary systems. However, few planets have been detected in blind surveys used to achieve imaging detections. Using Gaia and Hipparcos astrometry we identified dynamical evidence for a gas giant planet around the nearby star HIP 99770 and then confirmed this planet by direct imaging with the Subaru Coronagraphic Extreme Adaptive Optics Project. HIP 99770 b orbits 17 astronomical units from its host star, with an insolation comparable to Jupiter's and a dynamical mass of 13.9--16.1 Jupiter masses. Its planet-to-star mass ratio (7--8$\times$10$^{-3}$) is comparable to that other directly-imaged planets. The planet's atmosphere resembles an older, less-cloudy analogue of the atmospheres of previously-imaged exoplanets around HR 8799., Comment: Authors' version of revised paper, published in Science on April 14, 2023 (passed independent peer review/was recommended for publication in Science by external referees on July 5, 2022). First joint direct imaging + astrometric discovery of an exoplanet. 49 pages, 18 figures, 6 tables

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^{\prime\prime}$ (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^{+6}_{-4}\,\mathrm{M_{Jup}}$, which is $1{-}2{\sigma}$ 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^{\prime\prime}$ or $7.9 \times 10^4 \, \mathrm{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_{\mathrm{sed}}=2$, indicating the presence of clouds, $T_{\mathrm{eff}}=1400 \, \mathrm{K}$, and $\log{g}=4.5 \, \mathrm{dex}$. These results are consistent with the object's spectral type of $\mathrm{T0\pm1}$. 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., Comment: 40 pages, 19 figures, accepted to AJ

We present the latest and most precise characterization of the architecture for the ancient ($\approx 11$ Gyr) Kepler-444 system, which is composed of a K0 primary star (Kepler-444 A) hosting five transiting planets, and a tight M-type spectroscopic binary (Kepler-444 BC) with an A-BC projected separation of 66 au. We have measured the system's relative astrometry using the adaptive optics imaging from Keck/NIRC2 and Kepler-444 A's radial velocities from the Hobby Eberly Telescope, and re-analyzed relative radial velocities between BC and A from Keck/HIRES. We also include the Hipparcos-Gaia astrometric acceleration and all published astrometry and radial velocities into an updated orbit analysis of BC's barycenter. These data greatly extend the time baseline of the monitoring and lead to significant updates to BC's barycentric orbit compared to previous work, including a larger semi-major axis ($a = 52.2^{+3.3}_{-2.7}$ au), a smaller eccentricity ($e = 0.55 \pm 0.05$), and a more precise inclination ($i =85.4^{+0.3}_{-0.4}$ degrees). We have also derived the first dynamical masses of B and C components. Our results suggest Kepler-444~A's protoplanetary disk was likely truncated by BC to a radius of $\approx 8$ au, which resolves the previously noticed tension between Kepler-444 A's disk mass and planet masses. Kepler-444 BC's barycentric orbit is likely aligned with those of A's five planets, which might be primordial or a consequence of dynamical evolution. The Kepler-444 system demonstrates that compact multi-planet systems residing in hierarchical stellar triples can form at early epochs of the Universe and survive their secular evolution throughout cosmic time., Comment: AJ in press

We describe the discovery of a solar neighborhood (d=468 pc) binary system with a main-sequence sunlike star and a massive non-interacting black hole candidate. The spectral energy distribution (SED) of the visible star is described by a single stellar model. We derive stellar parameters from a high signal-to-noise Magellan/MIKE spectrum, classifying the star as a main-sequence star with $T_{\rm eff} = 5972 \rm K$, $\log{g} = 4.54$, and $M = 0.91$ \msun. The spectrum shows no indication of a second luminous component. To determine the spectroscopic orbit of the binary, we measured radial velocities of this system with the Automated Planet Finder, Magellan, and Keck over four months. We show that the velocity data are consistent with the \textit{Gaia} astrometric orbit and provide independent evidence for a massive dark companion. From a combined fit of our spectroscopic data and the astrometry, we derive a companion mass of $11.39^{+1.51}_{-1.31}$\msun. We conclude that this binary system harbors a massive black hole on an eccentric $(e =0.46 \pm 0.02)$, $185.4 \pm 0.1$ d orbit. These conclusions are independent of \cite{ElBadry2022Disc}, who recently reported the discovery of the same system. A joint fit to all available data (including \cite{ElBadry2022Disc}'s) yields a comparable period solution, but a lower companion mass of $9.32^{+0.22}_{-0.21} M_{\odot}$. Radial velocity fits to all available data produce a unimodal solution for the period that is not possible with either data set alone. The combination of both data sets yields the most accurate orbit currently available., Comment: 30 pages, 17 figures, accepted to AAS journals (various changes in response to referee comments)

Measured rotational speeds of giant planets and brown dwarfs frequently constitute appreciable fractions of the breakup limit, resulting in centrifugal expansion of these objects at the equator. According to models of internal energy transport, this expansion ought to make the poles of a rotator significantly hotter than the equator, so that inclination of the rotational axis greatly affects both spectral shape and total flux. In this article, we explore the dependence of a substellar object's observables on its rotational speed and axis inclination. To do so, we combine PICASO (a Planetary Intensity Code for Atmospheric Spectroscopy Observations) with software PARS (Paint the Atmospheres of Rotating Stars). The former computer program models radiative transfer within plane-parallel planetary atmospheres, while the latter computes disk-integrated spectra of centrifugally deformed gaseous masses. We find that the specific flux of a typical fast-rotating brown dwarf can increase by as much as a factor of 1.5 with movement from an equator-on to a pole-on view. On the other hand, the distinctive effect of rotation on spectral shape increases toward the equator-on view. The latter effect also increases with lower effective temperature. The bolometric luminosity estimate for a typical fast rotator at extreme inclinations has to be adjusted by as much as about 20% due to the anisotropy of the object's observed flux. We provide a general formula for the calculation of the corresponding adjustment factor in terms of rotational speed and inclination., Comment: 11 pages, 4 figures, MNRAS accepted

We present an experimental study of photon statistics for high-contrast imaging with the Microwave Kinetic Inductance Detector (MKID) Exoplanet Camera (MEC) located behind the Subaru Coronagraphic Extreme Adaptive Optics System (SCExAO) at the Subaru Telescope. We show that MEC measures the expected distributions for both on-axis companion intensity and off-axis intensity which manifests as quasi-static speckles in the image plane and currently limits high-contrast imaging performance. These statistics can be probed by any MEC observation due to the photon-counting capabilities of MKID detectors. Photon arrival time statistics can also be used to directly distinguish companions from speckles using a post-processing technique called Stochastic Speckle Discrimination (SSD). Here, we we give an overview of the SSD technique and highlight the first demonstration of SSD on an extended source -- the protoplanetary disk AB Aurigae. We then present simulations that provide an in-depth exploration as to the current limitations of an extension of the SSD technique called Photon-Counting SSD (PCSSD) to provide a path forward for transitioning PCSSD from simulations to on-sky results. We end with a discussion of how to further improve the efficacy of such arrival time based post-processing techniques applicable to both MKIDs, as well as other high speed astronomical cameras., Comment: 11 pages, 9 figures, accepted for publication in AJ

We present the direct imaging discovery of a low-mass companion to the nearby accelerating F star, HIP 5319, using SCExAO coupled with the CHARIS, VAMPIRES, and MEC instruments in addition to Keck/NIRC2 imaging. CHARIS $JHK$ (1.1-2.4 $\mu$m) spectroscopic data combined with VAMPIRES 750 nm, MEC $Y$, and NIRC2 $L_{\rm p}$ photometry is best matched by an M3--M7 object with an effective temperature of T=3200 K and surface gravity log($g$)=5.5. Using the relative astrometry for HIP 5319 B from CHARIS and NIRC2 and absolute astrometry for the primary from $Gaia$ and $Hipparcos$ and adopting a log-normal prior assumption for the companion mass, we measure a dynamical mass for HIP 5319 B of $31^{+35}_{-11}M_{\rm J}$, a semimajor axis of $18.6^{+10}_{-4.1}$ au, an inclination of $69.4^{+5.6}_{-15}$ degrees, and an eccentricity of $0.42^{+0.39}_{-0.29}$. However, using an alternate prior for our dynamical model yields a much higher mass of 128$^{+127}_{-88}M_{\rm J}$. Using data taken with the LCOGT NRES instrument we also show that the primary HIP 5319 A is a single star in contrast to previous characterizations of the system as a spectroscopic binary. This work underscores the importance of assumed priors in dynamical models for companions detected with imaging and astrometry and the need to have an updated inventory of system measurements., Comment: 19 pages, 9 figuresm 5 tables

The color-magnitude diagrams (CMDs) of intermediate-age star clusters (less than ~ 2 Gyr) are much more complex than those predicted by coeval, nonrotating stellar evolution models. Their observed extended main sequence turnoffs (eMSTOs) could result from variations in stellar age, stellar rotation, or both. The physical interpretation of eMSTOs is largely based on the complex mapping between stellar models -- themselves functions of mass, rotation, orientation, and binarity -- and the CMD. In this paper, we compute continuous probability densities in three-dimensional color, magnitude, and vsini space for individual stars in a cluster's eMSTO, based on a rotating stellar evolution model. These densities enable the rigorous inference of cluster properties from a stellar model, or, alternatively, constraints on the stellar model from the cluster's CMD. We use the MIST stellar evolution models to jointly infer the age dispersion, the rotational distribution, and the binary fraction of the Large Magellanic Cloud cluster NGC 1846. We derive an age dispersion of ~ 70-80 Myr, approximately half the earlier estimates due to nonrotating models. This finding agrees with the conjecture that rotational variation is largely responsible for eMSTOs. However, the MIST models do not provide a satisfactory fit to all stars in the cluster and achieve their best agreement at an unrealistically high binary fraction. The lack of agreement near the main-sequence turnoff suggests specific physical changes to the stellar evolution models, including a lower mass for the Kraft break and potentially enhanced main sequence lifespans for rapidly rotating stars., Comment: 34 pages, 13 figures, ApJ accepted. Python code is available at https://github.com/mlipatov/calc_cluster

White dwarfs have been successfully used as cosmochronometers in the literature, however their reach has been limited in comparison to their potential. We present wdwarfdate, a publicly available Python package to derive the Bayesian age of a white dwarf, based on its effective temperature (Teff) and surface gravity (logg). We make this software easy to use with the goal of transforming the usage of white dwarfs as cosmochronometers into an accessible tool. The code estimates the mass and cooling age of the white dwarf, as well as the mass and main-sequence age of the progenitor star, allowing for a determination of the total age of the object. We test the reliability of the method by estimating the parameters of white dwarfs from previous studies, and find agreement with the literature within measurement errors. By analyzing the limitation of the code we find a typical uncertainty of 10% on the total age when both input parameters have uncertainties of 1%, and an uncertainty of 25% on the total age when Teff has an uncertainty of 10% and logg of 1%. Furthermore, wdwarfdate assumes single star evolution, and can be applied to calculate the total age of a white dwarf with parameters in the range 1,500

We report the near-infrared radial-velocity (RV) discovery of a super-Earth planet on a 10.77-day orbit around the M4.5 dwarf Ross 508 ($J_\mathrm{mag}=9.1$). Using precision RVs from the Subaru Telescope IRD (InfraRed Doppler) instrument, we derive a semi-amplitude of $3.92^{+0.60}_{-0.58}$ ${\rm m\,s}^{-1}$, corresponding to a planet with a minimum mass $m \sin i = 4.00^{+0.53}_{-0.55}\ M_{\oplus}$. We find no evidence of significant signals at the detected period in spectroscopic stellar activity indicators or MEarth photometry. The planet, Ross 508 b, has a semimajor-axis of $0.05366^{+0.00056}_{-0.00049}$ au. This gives an orbit-averaged insolation of $\approx$1.4 times the Earth's value, placing Ross 508 b near the inner edge of its star's habitable zone. We have explored the possibility that the planet has a high eccentricity and its host is accompanied by an additional unconfirmed companion on a wide orbit. Our discovery demonstrates that the near-infrared RV search can play a crucial role to find a low-mass planet around cool M dwarfs like Ross 508., Comment: Accepted for publication in PASJ (May 23, 2022)

We report individual dynamical masses of $66.92 \pm 0.36 \; M_{Jup}$ and $53.25 \pm 0.29 \; M_{Jup}$ for the binary brown dwarfs $\varepsilon$ Indi Ba and Bb, measured from long term ($\approx 10$ yr) relative orbit monitoring and absolute astrometry monitoring data on the VLT. Relative astrometry with NACO fully constrains the Keplerian orbit of the binary pair, while absolute astrometry with FORS2 measures the system's parallax and mass ratio. We find a parallax consistent with the Hipparcos and Gaia values for $\varepsilon$ Indi A, and a mass ratio for $\varepsilon$ Indi Ba to Bb precise to better than $0.2\%$. $\varepsilon$ Indi Ba and Bb have spectral types T1-1.5 and T6, respectively. With an age of $3.5^{+0.8}_{-1.0}$ Gyr from $\varepsilon$ Indi A's activity, these brown dwarfs provide some of the most precise benchmarks for substellar cooling models. Assuming coevality, the very different luminosities of the two brown dwarfs and our moderate mass ratio imply a steep mass-luminosity relationship $L \propto M^{5.37 \pm 0.08}$ that can be explained by a slowed cooling rate in the L/T transition, as previously observed for other L/T binaries. Finally, we present a periodogram analysis of the near-infrared photometric data, but find no definitive evidence of periodic signals with a coherent phase., Comment: 21 pages, 14 figures, to be published in the Astronomical Journal

We present the direct-imaging discovery of a substellar companion in orbit around a Sun-like star member of the Hyades open cluster. So far, no other substellar companions have been unambiguously confirmed via direct imaging around main-sequence stars in Hyades. The star HIP 21152 is an accelerating star as identified by the astrometry from the Gaia and Hipparcos satellites. We have detected the companion, HIP 21152 B, in multi-epoch using the high-contrast imaging from SCExAO/CHARIS and Keck/NIRC2. We have also obtained the stellar radial-velocity data from the Okayama 188cm telescope. The CHARIS spectroscopy reveals that HIP 21152 B's spectrum is consistent with the L/T transition, best fit by an early T dwarf. Our orbit modeling determines the semi-major axis and the dynamical mass of HIP 21152 B to be 17.5$^{+7.2}_{-3.8}$ au and 27.8$^{+8.4}_{-5.4}$ $M_{\rm{Jup}}$, respectively. The mass ratio of HIP 21152 B relative to its host is $\approx$2\%, near the planet/brown dwarf boundary suggested from recent surveys. Mass estimates inferred from luminosity evolution models are slightly higher (33--42 $M_{\rm{Jup}}$). With a dynamical mass and a well-constrained age due to the system's Hyades membership, HIP 21152 B will become a critical benchmark in understanding the formation, evolution, and atmosphere of a substellar object as a function of mass and age. Our discovery is yet another key proof-of-concept for using precision astrometry to select direct imaging targets., Comment: 21 pages (11 pages in main body), 8 figures (4 figures in main body). Accepted for Publication in ApJL at July 9, 2022 (UT)

We present a spectrum of the diffuse Galactic light (DGL) between 3700 and 10,000 A, obtained by correlating optical sky intensity with far-infrared dust emission. We use nearly 250,000 blank-sky spectra from BOSS/SDSS-III together with IRIS-reprocessed maps from the IRAS satellite. The larger sample size compared to SDSS-II results in a factor-of-two increase in signal to noise. We combine these data sets with a model for the optical/far-infrared correlation that accounts for self-absorption by dust. The spectral features of the DGL agree remarkably well with features present in stellar spectra. There is evidence for a difference in the DGL continuum between the regions covered by BOSS in the northern and southern Galactic hemisphere. We interpret the difference at red wavelengths as the result of a difference in stellar populations, with mainly old stars in both regions but a higher fraction of young stars in the south. There is also a broad excess in the southern DGL spectrum over the prediction of a simple radiative transfer model, without a clear counterpart in the north. We interpret this excess, centered at ~6500 A, as evidence for luminescence in the form of extended red emission (ERE). The observed strength of the 4000 A break indicates that at most ~7% of the dust-correlated light at 4000 A can be due to blue luminescence. Our DGL spectrum provides constraints on dust scattering and luminescence independent of measurements of extinction., Comment: 16 pages, 8 figures. Revised version after reviewer comments

Gliese 86 is a nearby K dwarf hosting a giant planet on a $\approx$16-day orbit and an outer white dwarf companion on a $\approx$century-long orbit. In this study we combine radial velocity data (including new measurements spanning more than a decade) with high angular resolution imaging and absolute astrometry from Hipparcos and Gaia to measure the current orbits and masses of both companions. We then simulate the evolution of the Gl 86 system to constrain its primordial orbit when both stars were on the main sequence; the closest approach between the two stars was then about $9\,$AU. Such a close separation limited the size of the protoplanetary disk of Gl 86 A and dynamically hindered the formation of the giant planet around it. Our measurements of Gl 86 B and Gl 86 Ab's orbits reveal Gl 86 as a system in which giant planet formation took place in a disk truncated at $\approx$2$\,$AU. Such a disk would be just big enough to harbor the dust mass and total mass needed to assemble Gl 86 Ab's core and envelope, assuming a high disk accretion rate and a low viscosity. Inefficient accretion of the disk onto Gl 86 Ab, however, would require a disk massive enough to approach the Toomre stability limit at its outer truncation radius. The orbital architecture of the Gl 86 system shows that giant planets can form even in severely truncated disks and provides an important benchmark for planet formation theory., Comment: 17 pages, 10 figures

In this Letter, we measure the full orbital architecture of the two-planet system around the nearby K0 dwarf 14 Herculis. 14 Her (HD 145675, HIP 79248) is a middle-aged ($4.6^{+3.8}_{-1.3}$ Gyr) K0 star with two eccentric giant planets identified in the literature from radial velocity (RV) variability and long-term trends. Using archival RV data from Keck/HIRES in concert with \textit{Gaia-Hipparcos} acceleration in the proper motion vector for the star, we have disentangled the mass and inclination of the b planet to ${9.1}_{-1.1}^{+1.0}$ $M_\mathrm{Jup}$ and ${32.7}_{-3.2}^{+5.3}$ degrees. Despite only partial phase coverage for the c planet's orbit, we are able to constrain its mass and orbital parameters as well to ${6.9}_{-1.0}^{+1.7}$ $M_\mathrm{Jup}$ and ${101}_{-33}^{+31}$ degrees. We find that coplanarity of the b and c orbits is strongly disfavored. Combined with the age of the system and the comparable masses of its planets, this suggests that planet-planet scattering may be responsible for the current configuration of the system.

Model-independent masses of substellar companions are critical tools to validate models of planet and brown dwarf cooling, test their input physics, and determine the formation and evolution of these objects. In this work, we measure the dynamical mass and orbit of the young substellar companion HD 984 B. We obtained new high-contrast imaging of the HD 984 system with Keck/NIRC2 which expands the baseline of relative astrometry from 3 to 8 years. We also present new radial velocities of the host star with the Habitable-Zone Planet Finder spectrograph at the Hobby-Eberly Telescope. Furthermore, HD 984 exhibits a significant proper motion difference between Hipparcos and Gaia EDR3. Our joint orbit fit of the relative astrometry, proper motions, and radial velocities yields a dynamical mass of $61 \pm 4$ $\mathrm{M_{Jup}}$ for HD 984 B, placing the companion firmly in the brown dwarf regime. The new fit also reveals a higher eccentricity for the companion ($e = 0.76 \pm 0.05$) compared to previous orbit fits. Given the broad age constraint for HD 984, this mass is consistent with predictions from evolutionary models. HD 984 B's dynamical mass places it among a small but growing list of giant planet and brown dwarf companions with direct mass measurements., Comment: 23 pages, accepted to AJ

51 Eri b is one of the only young planets consistent with a wide range of possible initial entropy states, including the cold-start scenario associated with some models of planet formation by core accretion. The most direct way to constrain the initial entropy of a planet is by measuring its luminosity and mass at a sufficiently young age that the initial conditions still matter. We present the tightest upper limit on 51 Eri b's mass yet (M < 11 Mjup at 2$\sigma$) using a cross-calibration of Hipparcos and Gaia EDR3 astrometry and the orbit-fitting code orvara. We also reassess its luminosity using a direct, photometric approach, finding log(Lbol/Lsun) = -5.5$\pm$0.2 dex. Combining this luminosity with the 24$\pm$3 Myr age of the $\beta$ Pic moving group, of which 51 Eri is a member, we derive mass distributions from a grid of evolutionary models that spans a wide range of initial entropies. We find that 51 Eri b is inconsistent with the coldest-start scenarios, requiring an initial entropy of >8 $k_B$/baryon at 97% confidence. This result represents the first observational constraint on the initial entropy of a potentially cold-start planet, and it continues the trend of dynamical masses for directly imaged planets pointing to warm- or hot-start formation scenarios., Comment: Accepted for publication in MNRAS (9 pages, 6 figures)

We present orbits for 24 binaries in the field of open cluster NGC 2516 (~150 Myr) and 13 binaries in the field of open cluster NGC 2422 (~130 Myr) using results from a multi-year radial velocity survey of the cluster cores. Six of these systems are double-lined spectroscopic binaries (SB2s). We fit these RV variable systems with orvara, a MCMC-based fitting program that models Keplerian orbits. We use precise stellar parallaxes and proper motions from Gaia EDR3 to determine cluster membership. We impose a barycentric radial velocity prior on all cluster members; this significantly improves our orbital constraints. Two of our systems have periods between 5 and 15 days, the critical window in which tides efficiently damp orbital eccentricity. These binaries should be included in future analyses of circularization across similarly-aged clusters. We also find a relatively flat distribution of binary mass ratios, consistent with previous work. With the inclusion of TESS lightcurves for all available targets, we identity target 378-036252 as a new eclipsing binary. We also identify a field star whose secondary has a mass in the brown dwarf range, as well as two cluster members whose RVs suggest the presence of an additional companion. Our orbital fits will help constrain the binary fraction and binary properties across stellar age and across stellar environment., Comment: 38 pages, 47 figures (incl. 37 Appendix figures), 4 tables. AJ accepted

We present the SCExAO direct imaging discovery and characterization of a low-mass companion to the nearby young A7IV star, HD 91312. SCExAO/CHARIS $JHK$ (1.1-2.4 $\mu m$) spectra and SCExAO/HiCIAO $H$ band imaging identify the companion over a two year baseline in a highly inclined orbit with a maximum projected separation of 8 au. The companion, HD 91312 B, induces an 8.8-$\sigma$ astrometric acceleration on the star as seen with the Gaia & Hipparcos satellites and a long-term radial velocity trend as previously identified by Borgniet et al. (2019). HD 91312 B's spectrum is consistent with that of an early-to-mid M dwarf. Hipparcos and Gaia absolute astrometry, radial-velocity data, and SCExAO/CHARIS astrometry constrain its dynamical mass to be $0.337^{+0.042}_{-0.044}$M$_\odot$, consistent with - but far more precise than - masses derived from spectroscopy, and favors a nearly edge-on orbit with a semi-major axis of $\sim$9.7 au. This work is an example of precisely characterizing properties of low-mass companions at solar system-like scales from a combination of direct imaging, astrometry, and radial-velocity methods., Comment: 14 pages, 13 figures, Accepted for publication in AJ

Radial velocity (RV) surveys have discovered hundreds of exoplanetary systems but suffer from a fundamental degeneracy between planet mass $M_p$ and orbital inclination $i$. In this paper we break this degeneracy by combining RVs with complementary absolute astrometry taken from the Gaia EDR3 version of the cross-calibrated Hipparcos-Gaia Catalog of Accelerations (HGCA). We use the Markov Chain Monte Carlo orbit code $\tt orvara$ to simultaneously fit literature RVs and absolute astrometry from the HGCA. We constrain the orbits, masses, and inclinations of nine single and massive RV companions orbiting nearby G and K stars. We confirm the planetary nature of six companions: HD 29021 b ($4.47_{-0.65}^{+0.67}\,M_{\rm Jup}$), HD 81040 b ($7.24_{-0.37}^{+1.0}\,M_{\rm Jup}$), HD 87883 b ($6.31_{-0.32}^{+0.31}\,M_{\rm Jup}$), HD 98649 b ($9.7_{-1.9}^{+2.3}\,M_{\rm Jup}$), HD 106252 b ($10.00_{-0.73}^{+0.78}\,M_{\rm Jup}$), and HD 171238 b ($8.8_{-1.3}^{+3.6}\,M_{\rm Jup}$). We place one companion, HD 196067 b ($12.5_{-1.8}^{+2.5}\,M_{\rm Jup}$) on the planet-brown dwarf boundary, and two companions in the low mass brown dwarf regime: HD 106515 Ab ($18.9_{-1.4}^{+1.5}\,M_{\rm Jup}$), and HD 221420 b (${20.6}_{-1.6}^{+2.0}\,M_{\rm Jup}$). The brown dwarf HD 221420 b, with a semi-major axis of ${9.99}_{-0.70}^{+0.74}$ AU, a period of ${27.7}_{-2.5}^{+3.0}$ years, and an eccentricity of $0.162_{-0.030}^{+0.035}$ represents a promising target for high-contrast imaging. The RV orbits of HD 87883 b, HD 98649 b, HD 171238 b, and HD 196067 b are not fully constrained yet because of insufficient RV data. We find two possible inclinations for each of these orbits due to difficulty in separating prograde from retrograde orbits, but we expect this will change decisively with future Gaia data releases.

We present the first multi-wavelength (near-infrared; $1.1 - 2.4$ $\mu m$) imaging of HD 36546's debris disk, using the Subaru Coronagraphic Extreme Adaptive Optics (SCExAO) system coupled with the Coronagraphic High Angular Resolution Imaging Spectrograph (CHARIS). As a 3-10 Myr old star, HD 36546 presents a rare opportunity to study a debris disk at very early stages. SCExAO/CHARIS imagery resolves the disk over angular separations of $\rho \sim 0.25" - 1.0"$ (projected separations of $\rm{r_{proj}} \sim 25 - 101$ $\rm{au}$) and enables the first spectrophotometric analysis of the disk. The disk's brightness appears symmetric between its eastern and western extents and it exhibits slightly blue near-infrared colors on average (e.g. $J-K =-0.4\pm0.1$) $-$ suggesting copious sub-micron sized or highly porous grains. Through detailed modeling adopting a Hong scattering phase function (SPF), instead of the more common Henyey-Greenstein function, and using the differential evolution optimization algorithm, we provide an updated schematic of HD 36546's disk. The disk has a shallow radial dust density profile ($\alpha_{in} \approx 1.0$ and $\alpha_{out} \approx -1.5$), a fiducial radius of $r_0 \approx 82.7$ au, an inclination of $i \approx 79.1^\circ$, and a position angle of $\rm PA \approx 80.1^\circ$. Through spine tracing, we find a spine that is consistent with our modeling, but also with a "swept-back wing" geometry. Finally, we provide constraints on companions, including limiting a companion responsible for a marginal Hipparcos-Gaia acceleration to a projected separation of $\lesssim 0.2''$ and to a minimum mass of $\lesssim 11$ $\rm M_{Jup}$., Comment: 23 pages, 19 figures; corrected an error in the abstract, added a reference in Section 2

We present comprehensive orbital analyses and dynamical masses for the substellar companions Gl~229~B, Gl~758~B, HD~13724~B, HD~19467~B, HD~33632~Ab, and HD~72946~B. Our dynamical fits incorporate radial velocities, relative astrometry, and most importantly calibrated Hipparcos-Gaia EDR3 accelerations. For HD~33632~A and HD~72946 we perform three-body fits that account for their outer stellar companions. We present new relative astrometry of Gl~229~B with Keck/NIRC2, extending its observed baseline to 25 years. We obtain a $<$1\% mass measurement of $71.4 \pm 0.6\,M_{\rm Jup}$ for the first T dwarf Gl~229~B and a 1.2\% mass measurement of its host star ($0.579 \pm 0.007\,M_{\odot}$) that agrees with the high-mass-end of the M dwarf mass-luminosity relation. We perform a homogeneous analysis of the host stars' ages and use them, along with the companions' measured masses and luminosities, to test substellar evolutionary models. Gl~229~B is the most discrepant, as models predict that an object this massive cannot cool to such a low luminosity within a Hubble time, implying that it may be an unresolved binary. The other companions are generally consistent with models, except for HD~13724~B that has a host-star activity age 3.8$\sigma$ older than its substellar cooling age. Examining our results in context with other mass-age-luminosity benchmarks, we find no trend with spectral type but instead note that younger or lower-mass brown dwarfs are over-luminous compared to models, while older or higher-mass brown dwarfs are under-luminous. The presented mass measurements for some companions are so precise that the stellar host ages, not the masses, limit the analysis., Comment: Accepted for publication in AJ. References updated in version 2. See the journal version for the full quality figures. Figure sets and the MCMC chains (reduced to just 1000 samples however) are included with the journal version of the article, and pre-publication at https://drive.google.com/drive/folders/1_A8QYn9NyPgmGqJaY5sMHyT_wAS3uRRK?usp=sharing

We present htof, an open-source tool for interpreting and fitting the intermediate astrometric data (IAD) from both the 1997 and 2007 reductions of Hipparcos, the scanning-law of Gaia, and future missions such as the Nancy Grace Roman Space Telescope (NGRST). htof solves for the astrometric parameters of any system for any arbitrary combination of absolute astrometric missions. In preparation for later Gaia data releases, htof supports arbitrarily high-order astrometric solutions (e.g. five-, seven-, nine-parameter fits). Using htof, we find that the IAD of 6617 sources in Hipparcos 2007 might have been affected by a data corruption issue. htof integrates an ad-hoc correction that reconciles the IAD of these sources with their published catalog solutions. We developed htof to study masses and orbital parameters of sub-stellar companions, and we outline its implementation in one orbit fitting code (orvara, https://github.com/t-brandt/orvara). We use htof to predict a range of hypothetical additional planets in the $\beta$~Pic system, which could be detected by coupling NGRST astrometry with Gaia and Hipparcos. htof is pip installable and available at https://github.com/gmbrandt/htof ., Comment: Accepted to AJ. References updated in version 2. The Hipparcos 2007 Re-reduction Java Tool Intermediate Astrometric Data are available at , via the "zip file" link at https://www.cosmos.esa.int/web/hipparcos/hipparcos-2 : "...human readable version of the IAD of the Java tool in a zip file [warning: ~350 MB]..."

We describe a new high-contrast imaging capability well suited for studying planet-forming disks: near-infrared (NIR) high-contrast spectropolarimetric imaging with the Subaru Coronagraphic Extreme Adaptive Optics (SCExAO) system coupled with the Coronagraphic High Angular Resolution Imaging Spectrograph (CHARIS) integral field spectrograph (IFS). The advent of extreme adaptive optics (AO) systems, like SCExAO, has enabled recovery of planet-mass companions at the expected locations of gas-giant formation in young disks alongside disk structures (such as gaps or spirals) that may indicate protoplanet formation. In combination with SCExAO, the CHARIS IFS in polarimetry mode allows characterization of these systems at wavelengths spanning the NIR J, H, and K bands ($1.1-2.4$ $\mu m$, $R\sim20$) and at angular separations as small as 0.04". By comparing the resulting images with forward-modeled scattered light or 3D radiative-transfer models, the likely origins of any observed features can be assessed. Utilization of swift optimization algorithms, such as differential evolution (DE), to identify model parameters that best reproduce the observations allows plausible disk geometries to be explored efficiently. The recent addition of CHARIS's unique integral field spectropolarimetry mode has further facilitated the study of planet-forming disks -- aiding in the confirmation of candidate protoplanets, the diagnosis of disk structures, and the characterization of dust grain populations. We summarize preliminary results for two young planet-forming disk systems based on observations with the novel integral field spectropolarimetry mode for SCExAO/CHARIS., Comment: 20 pages, 8 figures

Giant stars as known exoplanet hosts are relatively rare due to the potential challenges in acquiring precision radial velocities and the small predicted transit depths. However, these giant host stars are also some of the brightest in the sky and so enable high signal-to-noise follow-up measurements. Here we report on new observations of the bright (V ~ 3.3) giant star $\iota$ Draconis ($\iota$ Dra), known to host a planet in a highly eccentric ~511 day period orbit. TESS observations of the star over 137 days reveal asteroseismic signatures, allowing us to constrain the stellar radius, mass, and age to ~2%, ~6%, and ~28%, respectively. We present the results of continued radial velocity monitoring of the star using the Automated Planet Finder over several orbits of the planet. We provide more precise planet parameters of the known planet and, through the combination of our radial velocity measurements with Hipparcos and Gaia astrometry, we discover an additional long-period companion with an orbital period of ~$68^{+60}_{-36}$ years. Mass predictions from our analysis place this sub-stellar companion on the border of the planet and brown dwarf regimes. The bright nature of the star combined with the revised orbital architecture of the system provides an opportunity to study planetary orbital dynamics that evolve as the star moves into the giant phase of its evolution., Comment: 20 pages, 11 Figures, Accepted for publication in the Astronomical Journal

HR 8799 hosts four directly imaged giant planets, but none has a mass measured from first principles. We present the first dynamical mass measurement in this planetary system, finding that the innermost planet HR~8799~e has a mass of $9.6^{+1.9}_{-1.8} \, M_{\rm Jup}$. This mass results from combining the well-characterized orbits of all four planets with a new astrometric acceleration detection (5$\sigma$) from the Gaia EDR3 version of the Hipparcos-Gaia Catalog of Accelerations. We find with 95\% confidence that HR~8799~e is below $13\, M_{\rm Jup}$, the deuterium-fusing mass limit. We derive a hot-start cooling age of $42^{+24}_{-16}$\,Myr for HR~8799~e that agrees well with its hypothesized membership in the Columba association but is also consistent with an alternative suggested membership in the $\beta$~Pictoris moving group. We exclude the presence of any additional $\gtrsim$5-$M_{\rm Jup}$ planets interior to HR~8799~e with semi-major axes between $\approx$3-16\,au. We provide proper motion anomalies and a matrix equation to solve for the mass of any of the planets of HR~8799 using only mass ratios between the planets., Comment: Accepted to ApJ Letters

We present a cross-calibration of Hipparcos and Gaia EDR3 intended to identify astrometrically accelerating stars and to fit orbits to stars with faint, massive companions. The resulting catalog, the EDR3 edition of the Hipparcos-Gaia Catalog of Accelerations (HGCA), provides three proper motions with calibrated uncertainties on the EDR3 reference frame: the Hipparcos proper motion, the Gaia EDR3 proper motion, and the long-term proper motion given by the difference in position between Hipparcos and Gaia EDR3. Our approach is similar to that for the Gaia DR2 edition of the HGCA, but offers a factor of ~3 improvement in precision thanks to the longer time baseline and improved data processing of Gaia EDR3. We again find that a 60/40 mixture of the two Hipparcos reductions outperforms either reduction individually, and we find strong evidence for locally variable frame rotations between all pairs of proper motion measurements. The substantial global frame rotation seen in DR2 proper motions has been removed in EDR3. We also correct for color- and magnitude-dependent frame rotations at a level of up to ~50 $\mu$as/yr in Gaia EDR3. We calibrate the Gaia EDR3 uncertainties using a sample of radial velocity standard stars without binary companions; we find an error inflation factor (a ratio of total to formal uncertainty) of 1.37. This is substantially lower than the position dependent factor of ~1.7 found for Gaia DR2 and reflects the improved data processing in EDR3. While the catalog should be used with caution, its proper motion residuals provide a powerful tool to measure the masses and orbits of faint, massive companions to nearby stars., Comment: 18 pages, 8 figures, 4 tables. ApJS accepted. Catalog available pre-publication at http://physics.ucsb.edu/~tbrandt/HGCA_vEDR3.fits

We present an open-source Python package, Orbits from Radial Velocity, Absolute, and/or Relative Astrometry (orvara), to fit Keplerian orbits to any combination of radial velocity, relative astrometry, and absolute astrometry data from the Hipparcos-Gaia Catalog of Accelerations. By combining these three data types, one can measure precise masses and sometimes orbital parameters even when the observations cover a small fraction of an orbit. orvara achieves its computational performance with an eccentric anomaly solver five to ten times faster than commonly used approaches, low-level memory management to avoid python overheads, and by analytically marginalizing out parallax, barycenter proper motion, and the instrument-specific radial velocity zero points. Through its integration with the Hipparcos and Gaia intermediate astrometry package htof, orvara can properly account for the epoch astrometry measurements of Hipparcos and the measurement times and scan angles of individual Gaia epochs. We configure orvara with modifiable .ini configuration files tailored to any specific stellar or planetary system. We demonstrate orvara with a case study application to a recently discovered white dwarf/main sequence (WD/MS) system, HD 159062. By adding absolute astrometry to literature RV and relative astrometry data, our comprehensive MCMC analysis improves the precision of HD 159062B's mass by more than an order of magnitude to $0.6083^{+0.0083}_{-0.0073}\,M_\odot$. We also derive a low eccentricity and large semimajor axis, establishing HD 159062AB as a system that did not experience Roche lobe overflow., Comment: 24 pages, 5 figures, 5 tables. AJ accepted with minor changes. orvara is available at https://github.com/t-brandt/orvara