Showing 1–50 of 148 results for author: Brandt, D

Astrometry as a Tool for Discovering and Weighing Faint Companions to Nearby Stars

Authors: Timothy D. Brandt

Abstract: 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 centu… ▽ More 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. △ Less

Submitted 17 June, 2024; originally announced June 2024.

Comments: 17 pages, 5 figures, PASP accepted. Tutorial intended primarily as an introduction and reference for non-experts

JWST/NIRCam 4-5 $μ$m Imaging of the Giant Planet AF Lep b

Authors: Kyle Franson, William O. Balmer, Brendan P. Bowler, Laurent Pueyo, Yifan Zhou, Emily Rickman, Zhoujian Zhang, Sagnick Mukherjee, Tim D. Pearce, Daniella C. Bardalez Gagliuffi, Lauren I. Biddle, Timothy D. Brandt, Rachel Bowens-Rubin, Justin R. Crepp, James W. Davidson, Jr., Jacqueline Faherty, Christian Ginski, Elliott P. Horch, Marvin Morgan, Caroline V. Morley, Marshall D. Perrin, Aniket Sanghi, Maissa Salama, Christopher A. Theissen, Quang H. Tran , et al. (1 additional authors not shown)

Abstract: 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… ▽ More 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{μ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{μ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. △ Less

Submitted 13 June, 2024; originally announced June 2024.

Comments: 17 pages, 4 figures, submitted to ApJL

Accurate and Model Independent Radius Determination of Single FGK and M Dwarfs Using Gaia DR3 Data

Authors: Rocio Kiman, Timothy D. Brandt, Jacqueline K. Faherty, Mark Popinchalk

Abstract: 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 ph… ▽ More 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<M[Msun]<= 0.6, but we found no correlation for lower masses. This could be due to lack of precision in our radius estimation or a physical reason. Radius measurements with smaller uncertainties are necessary to distinguish between the two scenarios. △ Less

Submitted 31 May, 2024; originally announced June 2024.

Comments: 21 pages, 10 figures, 1 table, Accepted for publication in AJ

Planet-Planet Scattering and ZLK Migration -- The Dynamical History of HAT-P-11

Authors: Tiger Lu, Qier An, Gongjie Li, Sarah C. Millholland, G. Mirek Brandt, Timothy D. Brandt

Abstract: 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 sc… ▽ More 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. △ Less

Submitted 29 May, 2024; originally announced May 2024.

Comments: 17 pages, 8 figures, submitted to ApJ

Significant mutual inclinations between the stellar spin and the orbits of both planets in the HAT-P-11 system

Authors: Qier An, Tiger Lu, G. Mirek Brandt, Timothy D Brandt, Gongjie Li

Abstract: Planet-star and planet-planet obliquity encode a planetary system's dynamical history, but both obliquities 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 HAT-P-11 sys… ▽ More Planet-star and planet-planet obliquity encode a planetary system's dynamical history, but both obliquities 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 HAT-P-11 system with astrometry and 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 3.06 pm 0.42 Jupiter mass, at a semimajor axis of 4.192 pm 0.07 AU, and planet b's minimum mass of 0.073 pm 0.0053 Jupiter mass. We present the posterior probability distribution of obliquity between star A and planet c, and between planet b and planet c. △ Less

Submitted 29 May, 2024; originally announced May 2024.

Comments: 12 pages, 8 figures, submitted to AJ

Gliese 12 b: A temperate Earth-sized planet at 12 pc ideal for atmospheric transmission spectroscopy

Authors: M. Kuzuhara, A. Fukui, J. H. Livingston, J. A. Caballero, J. P. de Leon, T. Hirano, Y. Kasagi, F. Murgas, N. Narita, M. Omiya, Jaume Orell-Miquel, E. Palle, Q. Changeat, E. Esparza-Borges, H. Harakawa, C. Hellier, Yasunori Hori, Kai Ikuta, H. T. Ishikawa, T. Kodama, T. Kotani, T. Kudo, J. C. Morales, M. Mori, E. Nagel , et al. (81 additional authors not shown)

Abstract: Recent discoveries of Earth-sized planets transiting nearby M dwarfs have made it possible to characterize the atmospheres of terrestrial planets via follow-up spectroscopic observations. However, the number of such planets receiving low insolation is still small, limiting our ability to understand the diversity of the atmospheric composition and climates of temperate terrestrial planets. We repor… ▽ More Recent discoveries of Earth-sized planets transiting nearby M dwarfs have made it possible to characterize the atmospheres of terrestrial planets via follow-up spectroscopic observations. However, the number of such planets receiving low insolation is still small, limiting our ability to understand the diversity of the atmospheric composition and climates of temperate terrestrial planets. We report the discovery of an Earth-sized planet transiting the nearby (12 pc) inactive M3.0 dwarf Gliese 12 (TOI-6251) with an orbital period ($P_{\rm{orb}}$) of 12.76 days. The planet, Gliese 12b, was initially identified as a candidate with an ambiguous $P_{\rm{orb}}$ from TESS data. We confirmed the transit signal and $P_{\rm{orb}}$ using ground-based photometry with MuSCAT2 and MuSCAT3, and validated the planetary nature of the signal using high-resolution images from Gemini/NIRI and Keck/NIRC2 as well as radial velocity (RV) measurements from the InfraRed Doppler instrument on the Subaru 8.2 m telescope and from CARMENES on the CAHA 3.5 m telescope. X-ray observations with XMM-Newton showed the host star is inactive, with an X-ray-to-bolometric luminosity ratio of $\log L_{\rm X}/L_{\rm bol} \approx -5.7$. Joint analysis of the light curves and RV measurements revealed that Gliese 12b has a radius of 0.96 $\pm$ 0.05 $R_\oplus$, a 3$σ$ mass upper limit of 3.9 $M_\oplus$, and an equilibrium temperature of 315 $\pm$ 6 K assuming zero albedo. The transmission spectroscopy metric (TSM) value of Gliese 12b is close to the TSM values of the TRAPPIST-1 planets, adding Gliese 12b to the small list of potentially terrestrial, temperate planets amenable to atmospheric characterization with JWST. △ Less

Submitted 23 May, 2024; originally announced May 2024.

Comments: 29 pages (20 pages in main body), 13 figures (10 figures in main body). Equal contributions from M. K. and A. F.. Accepted for Publication in ApJL at 2024 March 21

Journal ref: Published on 2024 May 23 by Astrophysical Journal Letters (ApJL) 967 L21

Likelihood-Based Jump Detection and Cosmic Ray Rejection for Detectors Read Out Up-the-Ramp

Authors: Timothy D. Brandt

Abstract: 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 $χ^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… ▽ More 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 $χ^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. △ Less

Submitted 28 March, 2024; originally announced April 2024.

Comments: 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

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

Authors: Taylor L. Tobin, Thayne Currie, Yiting Li, Jeffrey Chilcote, Timothy D. Brandt, Brianna Lacy, Masayuki Kuzuhara, Maria Vincent, Mona El Morsy, Vincent Deo, Jonathan P. Williams, Olivier Guyon, Julien Lozi, Sebastien Vievard, Nour Skaf, Kyohoon Ahn, Tyler Groff, N. Jeremy Kasdin, Taichi Uyama, Motohide Tamura, Aidan Gibbs, Briley L. Lewis, Rachel Bowens-Rubin, Maïssa Salama, Qier An , et al. (1 additional authors not shown)

Abstract: We present the direct-imaging discovery of a substellar companion (a massive planet or low-mass brown dwarf) to the young, $γ$ Doradus-type variable star, HIP 39017 (HD 65526). The companion's SCExAO/CHARIS JHK ($1.1-2.4μ$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… ▽ More We present the direct-imaging discovery of a substellar companion (a massive planet or low-mass brown dwarf) to the young, $γ$ Doradus-type variable star, HIP 39017 (HD 65526). The companion's SCExAO/CHARIS JHK ($1.1-2.4μ$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 $γ$ 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. △ Less

Submitted 15 May, 2024; v1 submitted 6 March, 2024; originally announced March 2024.

Comments: Published in AJ: April 9, 2024

Journal ref: 2024, AJ, 167, 205

Revised Architecture and Two New Super-Earths in the HD 134606 Planetary System

Authors: Zhexing Li, Stephen R. Kane, Timothy D. Brandt, Tara Fetherolf, Paul Robertson, Jinglin Zhao, Paul A. Dalba, Robert A. Wittenmyer, R. Paul Butler, Matias R. Diaz, Steve B. Howell, Jeremy Bailey, Brad Carter, Elise Furlan, Crystal L. Gnilka, Hugh R. A. Jones, Simon O'Toole, Chris Tinney

Abstract: 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 1… ▽ More 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. △ Less

Submitted 1 February, 2024; v1 submitted 30 January, 2024; originally announced January 2024.

Comments: 23 pages, 10 figures, 4 tables, accepted for publication in the Astronomical Journal

Two long-period giant planets around two giant stars: HD 112570 and HD 154391

Authors: Guang-Yao Xiao, Huan-Yu Teng, Jianzhao Zhou, Bun'ei Sato, Yu-Juan Liu, Shaolan Bi, Takuya Takarada, Masayuki Kuzuhara, Marc Hon, Liang Wang, Masashi Omiya, Hiroki Harakawa, Fei Zhao, Gang Zhao, Eiji Kambe, Hideyuki Izumiura, Hiroyasu Ando, Kunio Noguchi, Wei Wang, Meng Zhai, Nan Song, Chengqun Yang, Tanda Li, Timothy D. Brandt, Michitoshi Yoshida , et al. (2 additional authors not shown)

Abstract: 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 meta… ▽ More 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. △ Less

Submitted 3 December, 2023; originally announced December 2023.

Comments: 26 pages, 18 figures, 8 tables, Accepted for publication in The Astronomical Journal

Dynamical Architectures of S-type Transiting Planets in Binaries I: Target Selection using Hipparcos and Gaia proper motion anomalies

Authors: Jingwen Zhang, Lauren M. Weiss, Daniel Huber, Eric L. N. Jensen, Timothy D. Brandt, Karen Collins, Dennis M. Conti, Howard Isaacson, Pablo Lewin, Giuseppe Marino, Bob Massey, Felipe Murgas, Enric Palle, Don J. Radford, Howard M. Relles, Gregor Srdoc, Chris Stockdale, Thiam-Guan Tan, Gavin Wang

Abstract: 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 tr… ▽ More 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. △ Less

Submitted 5 March, 2024; v1 submitted 4 October, 2023; originally announced October 2023.

Comments: 28 pages, 9 figures, published by AJ

Journal ref: AJ, 167, 8 (2024)

Optimal Fitting and Debiasing for Detectors Read Out Up-the-Ramp

Authors: Timothy D. Brandt

Abstract: 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 esti… ▽ More 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 $χ^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. △ Less

Submitted 3 April, 2024; v1 submitted 15 September, 2023; originally announced September 2023.

Comments: 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

Post-processing CHARIS integral field spectrograph data with PyKLIP

Authors: Minghan Chen, Jason J. Wang, Timothy D. Brandt, Thayne Currie, Julien Lozi, Jeffrey Chilcote, Maria Vincent

Abstract: 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… ▽ More 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. △ Less

Submitted 7 September, 2023; originally announced September 2023.

Comments: 17 pages, 13 figures

Journal ref: RAS Techniques and Instruments, rzad039 (2023)

Revisiting Planetary Systems in Okayama Planet Search Program: A new long-period planet, RV astrometry joint analysis, and multiplicity-metallicity trend around evolved stars

Authors: Huan-Yu Teng, Bun'ei Sato, Masayuki Kuzuhara, Takuya Takarada, Masashi Omiya, Hiroki Harakawa, Hideyuki Izumiura, Eiji Kambe, Mesut Yilmaz, Ilfan Bikmaev, Selim O. Selam, Timothy D. Brandt, Guang-Yao Xiao, Michitoshi Yoshida, Yoichi Itoh, Hiroyasu Ando, Eiichiro Kokubo, Shigeru Ida

Abstract: 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 con… ▽ More 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 $ε$ Tau, 11 Com, 24 Boo, 41 Lyn, 14 And, HD 32518, and $ω$ 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. △ Less

Submitted 10 August, 2023; originally announced August 2023.

Comments: 49 figures, 4 tables, accepted by PASJ, RV data will be available online as supplementary after the publication

The Masses of a Sample of Radial-Velocity Exoplanets with Astrometric Measurements

Authors: Guang-Yao Xiao, Yu-Juan Liu, Huan-Yu Teng, Wei Wang, Timothy D. Brandt, Gang Zhao, Fei Zhao, Meng Zhai, Qi Gao

Abstract: 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,… ▽ More 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. △ Less

Submitted 22 March, 2023; originally announced March 2023.

Comments: 37 pages, 14 figures, accepted by Research in Astronomy and Astrophysics

Dynamical Masses and Ages of Sirius-like Systems

Authors: Hengyue Zhang, Timothy D. Brandt, Rocio Kiman, Alexander Venner, Qier An, Minghan Chen, Yiting Li

Abstract: We measure precise orbits and dynamical masses and derive age constraints for six confirmed and one candidate Sirius-like systems, including the Hyades member HD 27483. Our orbital analysis incorporates radial velocities, relative astrometry, and Hipparcos-Gaia astrometric accelerations. We constrain the main-sequence lifetime of a white dwarf's progenitor from the remnant's dynamical mass and sem… ▽ More 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. △ Less

Submitted 9 July, 2023; v1 submitted 14 March, 2023; originally announced March 2023.

Comments: 21 pages, 7 figures

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 524, Issue 1, September 2023, Pages 695-715

Astrometric Accelerations as Dynamical Beacons: A Giant Planet Imaged Inside the Debris Disk of the Young Star AF Lep

Authors: Kyle Franson, Brendan P. Bowler, Yifan Zhou, Tim D. Pearce, Daniella C. Bardalez Gagliuffi, Lauren Biddle, Timothy D. Brandt, Justin R. Crepp, Trent J. Dupuy, Jacqueline Faherty, Rebecca Jensen-Clem, Marvin Morgan, Aniket Sanghi, Christopher A. Theissen, Quang H. Tran, Trevor A. Wolf

Abstract: 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 $β$ 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… ▽ More 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 $β$ 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-$σ$ 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. △ Less

Submitted 25 May, 2023; v1 submitted 10 February, 2023; originally announced February 2023.

Comments: 14 pages, 3 figures, accepted to ApJL

Surveying Nearby Brown Dwarfs with HGCA: Direct Imaging Discovery of a Faint, High-Mass Brown Dwarf Orbiting HD 176535 A

Authors: Yiting Li, Timothy D. Brandt, G. Mirek Brandt, Qier An, Kyle Franson, Trent J. Dupuy, Minghan Chen, Rachel Bowens-Rubin, Briley L. Lewis, Brendan P. Bowler, Aidan Gibbs, Rocio Kiman, Jacqueline Faherty, Thayne Currie, Rebecca Jensen-Clem, Hengyue Zhang Ezequiel Contreras-Martinez, Michael P. Fitzgerald, Benjamin A. Mazin, Maxwell Millar-Blanchaer

Abstract: Brown dwarfs with well-measured masses, ages and luminosities provide direct benchmark tests of substellar formation and evolutionary models. We report the first results from a direct imaging survey aiming to find and characterize substellar companions to nearby accelerating stars with the assistance of the Hipparcos-Gaia Catalog of Accelerations (HGCA). In this paper, we present a joint high-cont… ▽ More 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 $Δ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. △ Less

Submitted 16 May, 2023; v1 submitted 25 January, 2023; originally announced January 2023.

TESS-Gaia Light Curve: a PSF-based TESS FFI light curve product

Authors: Te Han, Timothy D. Brandt

Abstract: 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 astro… ▽ More 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. △ Less

Submitted 9 January, 2023; originally announced January 2023.

Comments: 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

Journal ref: AJ 165 71 (2023)

Direct Imaging and Astrometric Detection of a Gas Giant Planet Orbiting an Accelerating Star

Authors: Thayne Currie, G. Mirek Brandt, Timothy D. Brandt, Brianna Lacy, Adam Burrows, Olivier Guyon, Motohide Tamura, Ranger Y. Liu, Sabina Sagynbayeva, Taylor Tobin, Jeffrey Chilcote, Tyler Groff, Christian Marois, William Thompson, Simon Murphy, Masayuki Kuzuhara, Kellen Lawson, Julien Lozi, Vincent Deo, Sebastien Vievard, Nour Skaf, Taichi Uyama, Nemanja Jovanovic, Frantz Martinache, N. Jeremy Kasdin , et al. (9 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 13 April, 2023; v1 submitted 30 November, 2022; originally announced December 2022.

Comments: 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

Astrometric Accelerations as Dynamical Beacons: Discovery and Characterization of HIP 21152 B, the First T-Dwarf Companion in the Hyades

Authors: Kyle Franson, Brendan P. Bowler, Mariangela Bonavita, Timothy D. Brandt, Minghan Chen, Matthias Samland, Zhoujian Zhang, Anna Lueber, Kevin Heng, Daniel Kitzmann, Trevor Wolf, Brandon A. Jones, Quang H. Tran, Daniella C. Bardalez Gagliuffi, Beth Biller, Jeffrey Chilcote, Justin R. Crepp, Trent J. Dupuy, Jacqueline Faherty, Clemence Fontanive, Tyler D. Groff, Raffaele Gratton, Olivier Guyon, Rebecca Jensen-Clem, Nemanja Jovanovic , et al. (6 additional authors not shown)

Abstract: 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.… ▽ More 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σ$ 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. △ Less

Submitted 17 November, 2022; originally announced November 2022.

Comments: 40 pages, 19 figures, accepted to AJ

The McDonald Accelerating Stars Survey (MASS): Architecture of the Ancient Five-Planet Host System Kepler-444

Authors: Zhoujian Zhang, Brendan P. Bowler, Trent J. Dupuy, Timothy D. Brandt, G. Mirek Brandt, William D. Cochran, Michael Endl, Phillip J. MacQueen, Kaitlin M. Kratter, Howard T. Isaacson, Kyle Franson, Adam L. Kraus, Caroline V. Morley, Yifan Zhou

Abstract: 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… ▽ More 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. △ Less

Submitted 13 October, 2022; originally announced October 2022.

Comments: AJ in press

A non-interacting Galactic black hole candidate in a binary system with a main-sequence star

Authors: Sukanya Chakrabarti, Joshua D. Simon, Peter A. Craig, Henrique Reggiani, Timothy D. Brandt, Puragra Guhathakurta, Paul A. Dalba, Evan N. Kirby, Philip Chang, Daniel R. Hey, Alessandro Savino, Marla Geha, Ian B. Thompson

Abstract: 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… ▽ More 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. △ Less

Submitted 5 June, 2023; v1 submitted 10 October, 2022; originally announced October 2022.

Comments: 30 pages, 17 figures, accepted to AAS journals (various changes in response to referee comments)

Effects of Rotation on the Spectra of Brown Dwarfs

Authors: Mikhail Lipatov, Timothy D. Brandt, Natasha E. Batalha

Abstract: 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… ▽ More 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. △ Less

Submitted 29 September, 2022; originally announced September 2022.

Comments: 11 pages, 4 figures, MNRAS accepted

Precise dynamical masses of new directly imaged companions from combining relative astrometry, radial velocities, and Hipparcos-Gaia eDR3 accelerations

Authors: E. L. Rickman, E. Matthews, W. Ceva, D. Ségransan, G. M. Brandt, H. Zhang, T. D. Brandt, T. Forveille, J. Hagelberg, S. Udry

Abstract: Aims. With an observing time span of more than 20 years, the CORALIE radial-velocity survey is able to detect long-term trends in data corresponding to companions with masses and separations accessible to direct imaging. Combining exoplanet detection techniques such as radial velocities from the CORALIE survey, astrometric accelerations from Hipparcos and Gaia eDR3, and relative astrometry from di… ▽ More Aims. With an observing time span of more than 20 years, the CORALIE radial-velocity survey is able to detect long-term trends in data corresponding to companions with masses and separations accessible to direct imaging. Combining exoplanet detection techniques such as radial velocities from the CORALIE survey, astrometric accelerations from Hipparcos and Gaia eDR3, and relative astrometry from direct imaging, removes the degeneracy of unknown orbital parameters. This allows precise model-independent masses of detected companions to be derived, which provides a powerful tool to test models of stellar and substellar mass-luminosity relations. Methods. Long-term precise Doppler measurements with the CORALIE spectrograph reveal radial-velocity signatures of companions on long-period orbits. The long baseline of radial-velocity data allows the detectability of such companion candidates to be assessed with direct imaging. We combine long-period radial-velocity data with absolute astrometry from Hipparcos and Gaia eDR3 and relative astrometry derived from new direct imaging detections with VLT/SPHERE to fit orbital parameters and derive precise dynamical masses of these companions. Results. In this paper we report the discovery of new companions orbiting HD~142234, HD~143616, and HIP~22059, as well as the first direct detection of HD~92987~B, and update the dynamical masses of two previously directly imaged companions; HD~157338~B and HD~195010~B. The companions span a period range of 32 to 279 years and are all very low mass stellar companions, ranging from 218 to 487~$M_{\rm{Jup}}$. We compare the derived dynamical masses to mass-luminosity relations of very low mass stars (<0.5~$M_{\odot}$), and discuss the importance of using precursor radial-velocity and astrometric information to inform the future of high-contrast imaging of exoplanets and brown dwarfs △ Less

Submitted 3 November, 2022; v1 submitted 26 September, 2022; originally announced September 2022.

Comments: Accepted for publication in A&A, updated version including language edits

Journal ref: A&A 668, A140 (2022)

Probing Photon Statistics in Adaptive Optics Images with SCExAO/MEC

Authors: Sarah Steiger, Timothy D. Brandt, Olivier Guyon, Noah Swimmer, Alexander B. Walter, Clinton Bockstiegel, Julien Lozi, Vincent Deo, Sebastien Vievard, Nour Skaf, Kyohoon Ahn, Nemanja Jovanovic, Frantz Martinache, Benjamin A. Mazin

Abstract: 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-… ▽ More 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. △ Less

Submitted 13 September, 2022; originally announced September 2022.

Comments: 11 pages, 9 figures, accepted for publication in AJ

SCExAO and Keck Direct Imaging Discovery of a Low-Mass Companion Around the Accelerating F5 Star HIP 5319

Authors: Noah Swimmer, Thayne Currie, Sarah Steiger, Gregory Mirek Brandt, Timothy D. Brandt, Olivier Guyon, Masayuki Kuzuhara, Jeffrey Chilcote, Taylor Tobin, Tyler D. Groff, Julien Lozi, John I. Bailey III, Alexander B. Walter, Neelay Fruitwala, Nicholas Zobrist, Jennifer Pearl Smith, Gregoire Coiffard, Rupert Dodkins, Kristina K. Davis, Miguel Daal, Bruce Bumble, Sebastien Vievard, Nour Skaf, Vincent Deo, Nemanja Jovanovic , et al. (4 additional authors not shown)

Abstract: 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 $μ$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… ▽ More 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 $μ$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. △ Less

Submitted 30 July, 2022; originally announced August 2022.

Comments: 19 pages, 9 figuresm 5 tables

Rotational Variation Allows for Narrow Age Spread in the Extended Main Sequence Turnoff of Massive Cluster NGC 1846

Authors: Mikhail Lipatov, Timothy D. Brandt, Seth Gossage

Abstract: 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 stel… ▽ More 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. △ Less

Submitted 16 June, 2022; originally announced June 2022.

Comments: 34 pages, 13 figures, ApJ accepted. Python code is available at https://github.com/mlipatov/calc_cluster

wdwarfdate: A Python Package to Derive Bayesian Ages of White Dwarfs

Authors: Rocio Kiman, Siyi Xu, Jacqueline K. Faherty, Jonathan Gagne, Ruth Angus, Timothy D. Brandt, Sarah L. Casewell, Kelle L. Cruz

Abstract: 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… ▽ More 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<Teff<90,000 K and 7.9<logg<9.3. Finally, the code assumes a uniform mixture of C/O in the core and single star evolution, which is reliable in the range of white dwarf masses 0.45-1.1 Msun (7.73<logg<8.8). △ Less

Submitted 10 June, 2022; originally announced June 2022.

Comments: 19 pages, 12 figures, 2 tables, Accepted for publication in AJ

A Super-Earth Orbiting Near the Inner Edge of the Habitable Zone around the M4.5-dwarf Ross 508

Authors: Hiroki Harakawa, Takuya Takarada, Yui Kasagi, Teruyuki Hirano, Takayuki Kotani, Masayuki Kuzuhara, Masashi Omiya, Hajime Kawahara, Akihiko Fukui, Yasunori Hori, Hiroyuki Tako Ishikawa, Masahiro Ogihara, John Livingston, Timothy D. Brandt, Thayne Currie, Wako Aoki, Charles A. Beichman, Thomas Henning, Klaus Hodapp, Masato Ishizuka, Hideyuki Izumiura, Shane Jacobson, Markus Janson, Eiji Kambe, Takanori Kodama , et al. (24 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 24 May, 2022; originally announced May 2022.

Comments: Accepted for publication in PASJ (May 23, 2022)

Precise Dynamical Masses of Epsilon Indi Ba and Bb: Evidence of Slowed Cooling at the L/T Transition

Authors: Minghan Chen, Yiting Li, Timothy D. Brandt, Trent J. Dupuy, Cátia V. Cardoso, Mark J. McCaughrean

Abstract: 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 FOR… ▽ More 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. △ Less

Submitted 17 May, 2022; originally announced May 2022.

Comments: 21 pages, 14 figures, to be published in the Astronomical Journal

Direct Imaging Discovery and Dynamical Mass of a Substellar Companion Orbiting an Accelerating Hyades Sun-like Star with SCExAO/CHARIS

Authors: Masayuki Kuzuhara, Thayne Currie, Takuya Takarada, Timothy D. Brandt, Bun'ei Sato, Taichi Uyama, Markus Janson, Jeffrey Chilcote, Taylor Tobin, Kellen Lawson, Yasunori Hori, Olivier Guyon, Tyler D. Groff, Julien Lozi, Sebastien Vievard, Ananya Sahoo, Vincent Deo, Nemanja Jovanovic, Kyohoon Ahn, Frantz Martinache, Nour Skaf, Eiji Akiyama, Barnaby R. Norris, Mickael Bonnefoy, Krzysztof G. Hełminiak , et al. (11 additional authors not shown)

Abstract: 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 c… ▽ More 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. △ Less

Submitted 12 June, 2022; v1 submitted 5 May, 2022; originally announced May 2022.

Comments: 21 pages (11 pages in main body), 8 figures (4 figures in main body). Accepted for Publication in ApJL at July 9, 2022 (UT)

An Optical Spectrum of the Diffuse Galactic Light from BOSS and IRIS

Authors: Blake Chellew, Timothy D. Brandt, Brandon S. Hensley, Bruce T. Draine, Eve Matthaey

Abstract: 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… ▽ More 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. △ Less

Submitted 4 July, 2022; v1 submitted 4 January, 2022; originally announced January 2022.

Comments: 16 pages, 8 figures. Revised version after reviewer comments

Journal ref: The Astrophysical Journal, Volume 932, Number 2. June 22, 2022

The Gliese 86 Binary System: A Warm Jupiter Formed in a Disk Truncated at $\approx$2 AU

Authors: Yunlin Zeng, Timothy D. Brandt, Gongjie Li, Trent J. Dupuy, Yiting Li, G. Mirek Brandt, Jay Farihi, Jonathan Horner, Robert A. Wittenmyer, R. Paul. Butler, Christopher G. Tinney, Bradley D. Carter, Duncan J. Wright, Hugh R. A. Jones, Simon J. O'Toole

Abstract: 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… ▽ More 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. △ Less

Submitted 12 December, 2021; originally announced December 2021.

Comments: 17 pages, 10 figures

14 Her: a likely case of planet-planet scattering

Authors: Daniella C. Bardalez Gagliuffi, Jacqueline K. Faherty, Yiting Li, Timothy D. Brandt, Lauryn Williams, G. Mirek Brandt, Christopher R. Gelino

Abstract: 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-Hipparc… ▽ More 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. △ Less

Submitted 12 November, 2021; v1 submitted 10 November, 2021; originally announced November 2021.

Dynamical Mass of the Young Substellar Companion HD 984 B

Authors: Kyle Franson, Brendan P. Bowler, Timothy D. Brandt, Trent J. Dupuy, Quang H. Tran, G. Mirek Brandt, Yiting Li, Adam L. Kraus

Abstract: 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… ▽ More 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. △ Less

Submitted 2 November, 2021; originally announced November 2021.

Comments: 23 pages, accepted to AJ

Limits on the Mass and Initial Entropy of 51 Eri b from Gaia EDR3 Astrometry

Authors: Trent J. Dupuy, G. Mirek Brandt, Timothy D. Brandt

Abstract: 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 ti… ▽ More 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$σ$) 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 $β$ 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. △ Less

Submitted 25 October, 2021; originally announced October 2021.

Comments: Accepted for publication in MNRAS (9 pages, 6 figures)

Orbital Parameters and Binary Properties of 37 FGK stars in the Cores of Open Clusters NGC 2516 and NGC 2422

Authors: Isabel Lipartito, John I. Bailey, Timothy D. Brandt, Benjamin A. Mazin, Mario Mateo, Meghin E. Spencer, Ian U. Roederer

Abstract: 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. W… ▽ More 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. △ Less

Submitted 5 October, 2021; originally announced October 2021.

Comments: 38 pages, 47 figures (incl. 37 Appendix figures), 4 tables. AJ accepted

SCExAO/CHARIS Direct Imaging of A Low-Mass Companion At A Saturn-Like Separation from an Accelerating Young A7 Star

Authors: Jeffrey Chilcote, Taylor Tobin, Thayne Currie, Timothy D. Brandt, Tyler D. Groff, Masayuki Kuzuhara, Olivier Guyon, Julien Lozi, Nemanja Jovanovic, Ananya Sahoo, Vincent Deo, Eiji Akiyama, Markus Janson, Jill Knapp, Jungmi Kwon, Michael W. McElwain, Jun Nishikawa, Kevin Wagner, Krzysztof Hełminiak, Motohide Tamura

Abstract: 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 $μ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-$σ$ astrometric acceleratio… ▽ More 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 $μ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-$σ$ 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. △ Less

Submitted 24 September, 2021; originally announced September 2021.

Comments: 14 pages, 13 figures, Accepted for publication in AJ

Precise Masses and Orbits for Nine Radial Velocity Exoplanets

Authors: Yiting Li, Timothy D. Brandt, G. Mirek Brandt, Trent J. Dupuy, Daniel Michalik, Rebecca Jensen-Clem, Yunlin Zeng, Jacqueline Faherty, Elena L. Mitra

Abstract: 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 Mo… ▽ More 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. △ Less

Submitted 4 October, 2021; v1 submitted 21 September, 2021; originally announced September 2021.

Multiband imaging of the HD 36546 debris disk: a refined view from SCExAO/CHARIS

Authors: Kellen Lawson, Thayne Currie, John P. Wisniewski, Motohide Tamura, Jean-Charles Augereau, Timothy D. Brandt, Olivier Guyon, N. Jeremy Kasdin, Tyler D. Groff, Julien Lozi, Vincent Deo, Sebastien Vievard, Jeffrey Chilcote, Nemanja Jovanovic, Frantz Martinache, Nour Skaf, Thomas Henning, Gillian Knapp, Jungmi Kwon, Michael W. McElwain, Tae-Soo Pyo, Michael L. Sitko, Taichi Uyama, Kevin Wagner

Abstract: We present the first multi-wavelength (near-infrared; $1.1 - 2.4$ $μ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… ▽ More We present the first multi-wavelength (near-infrared; $1.1 - 2.4$ $μ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 $ρ\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 ($α_{in} \approx 1.0$ and $α_{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}$. △ Less

Submitted 2 November, 2021; v1 submitted 18 September, 2021; originally announced September 2021.

Comments: 23 pages, 19 figures; corrected an error in the abstract, added a reference in Section 2

Improved Dynamical Masses for Six Brown Dwarf Companions Using Hipparcos and Gaia EDR3

Authors: G. Mirek Brandt, Trent J. Dupuy, Yiting Li, Minghan Chen, Timothy D. Brandt, Tin Long Sunny Wong, Thayne Currie, Brendan P. Bowler, Michael C. Liu, William M. J. Best, Mark W. Phillips

Abstract: 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 comp… ▽ More 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$σ$ 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. △ Less

Submitted 30 September, 2021; v1 submitted 15 September, 2021; originally announced September 2021.

Comments: 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

htof: A new open-source tool for analyzing Hipparcos, Gaia, and future astrometric missions

Authors: G. Mirek Brandt, Daniel Michalik, Timothy D. Brandt, Yiting Li, Trent J. Dupuy, Yunlin Zeng

Abstract: 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… ▽ More 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 $β$~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 . △ Less

Submitted 30 September, 2021; v1 submitted 14 September, 2021; originally announced September 2021.

Comments: 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]..."

High-contrast integral field spectropolarimetry of planet-forming disks with SCExAO/CHARIS

Authors: Kellen Lawson, Thayne Currie, John P. Wisniewski, Jun Hashimoto, Olivier Guyon, N. Jeremy Kasdin, Tyler D. Groff, Julien Lozi, Timothy D. Brandt, Jeffrey Chilcote, Vincent Deo, Taichi Uyama, Sebastien Vievard

Abstract: 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… ▽ More 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$ $μ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. △ Less

Submitted 19 August, 2021; originally announced August 2021.

Comments: 20 pages, 8 figures

Asteroseismology of iota Draconis and Discovery of an Additional Long-Period Companion

Authors: Michelle L. Hill, Stephen R. Kane, Tiago L. Campante, Zhexing Li, Paul A. Dalba, Timothy D. Brandt, Timothy R. White, Benjamin J. S. Pope, Keivan G. Stassun, Benjamin J. Fulton, Enrico Corsaro, Tanda Li, J. M. Joel Ong, Timothy R. Bedding, Diego Bossini, Derek L. Buzasi, William J. Chaplin, Margarida S. Cunha, Rafael A. Garcia, Sylvain N. Breton, Marc Hon, Daniel Huber, Chen Jiang, Cenk Kayhan, James S. Kuszlewicz , et al. (3 additional authors not shown)

Abstract: 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 $ι$ Draconis ($ι$ Dra), kno… ▽ More 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 $ι$ Draconis ($ι$ 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. △ Less

Submitted 28 July, 2021; originally announced July 2021.

Comments: 20 pages, 11 Figures, Accepted for publication in the Astronomical Journal

The First Dynamical Mass Measurement in the HR 8799 System

Authors: G. Mirek Brandt, Timothy D. Brandt, Trent J. Dupuy, Daniel Michalik, Gabriel-Dominique Marleau

Abstract: 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 (… ▽ More 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$σ$) 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 $β$~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. △ Less

Submitted 26 May, 2021; originally announced May 2021.

Comments: Accepted to ApJ Letters

orvara: An Efficient Code to Fit Orbits using Radial Velocity, Absolute, and/or Relative Astrometry

Authors: Timothy D. Brandt, Trent J. Dupuy, Yiting Li, G. Mirek Brandt, Yunlin Zeng, Daniel Michalik, Daniella C. Bardalez Gagliuffi, Virginia Raposo-Pulido

Abstract: 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 obse… ▽ More 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. △ Less

Submitted 15 September, 2021; v1 submitted 25 May, 2021; originally announced May 2021.

Comments: 24 pages, 5 figures, 5 tables. AJ accepted with minor changes. orvara is available at https://github.com/t-brandt/orvara

The Hipparcos-Gaia Catalog of Accelerations: Gaia EDR3 Edition

Authors: Timothy D. Brandt

Abstract: 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… ▽ More 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 $μ$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. △ Less

Submitted 25 May, 2021; originally announced May 2021.

Comments: 18 pages, 8 figures, 4 tables. ApJS accepted. Catalog available pre-publication at http://physics.ucsb.edu/~tbrandt/HGCA_vEDR3.fits

exoplanet: Gradient-based probabilistic inference for exoplanet data & other astronomical time series

Authors: Daniel Foreman-Mackey, Rodrigo Luger, Eric Agol, Thomas Barclay, Luke G. Bouma, Timothy D. Brandt, Ian Czekala, Trevor J. David, Jiayin Dong, Emily A. Gilbert, Tyler A. Gordon, Christina Hedges, Daniel R. Hey, Brett M. Morris, Adrian M. Price-Whelan, Arjun B. Savel

Abstract: "exoplanet" is a toolkit for probabilistic modeling of astronomical time series data, with a focus on observations of exoplanets, using PyMC3 (Salvatier et al., 2016). PyMC3 is a flexible and high-performance model-building language and inference engine that scales well to problems with a large number of parameters. "exoplanet" extends PyMC3's modeling language to support many of the custom functi… ▽ More "exoplanet" is a toolkit for probabilistic modeling of astronomical time series data, with a focus on observations of exoplanets, using PyMC3 (Salvatier et al., 2016). PyMC3 is a flexible and high-performance model-building language and inference engine that scales well to problems with a large number of parameters. "exoplanet" extends PyMC3's modeling language to support many of the custom functions and probability distributions required when fitting exoplanet datasets or other astronomical time series. While it has been used for other applications, such as the study of stellar variability, the primary purpose of "exoplanet" is the characterization of exoplanets or multiple star systems using time-series photometry, astrometry, and/or radial velocity. In particular, the typical use case would be to use one or more of these datasets to place constraints on the physical and orbital parameters of the system, such as planet mass or orbital period, while simultaneously taking into account the effects of stellar variability. △ Less

Submitted 24 June, 2021; v1 submitted 5 May, 2021; originally announced May 2021.

Comments: Published in the Journal of Open Source Software. Comments (still) welcome. Software available at https://docs.exoplanet.codes

The McDonald Accelerating Stars Survey (MASS): Discovery of a Long-Period Substellar Companion Orbiting the Old Solar Analog HD 47127

Authors: Brendan P. Bowler, Michael Endl, William D. Cochran, Phillip J. MacQueen, Justin R. Crepp, Greg W. Doppmann, Shannon Dulz, Timothy D. Brandt, G. Mirek Brandt, Yiting Li, Trent J. Dupuy, Kyle Franson, Kaitlin M. Kratter, Caroline V. Morley, Yifan Zhou

Abstract: Brown dwarfs with well-determined ages, luminosities, and masses provide rare but valuable tests of low-temperature atmospheric and evolutionary models. We present the discovery and dynamical mass measurement of a substellar companion to HD 47127, an old ($\approx$7-10 Gyr) G5 main sequence star with a mass similar to the Sun. Radial velocities of the host star with the Harlan J. Smith Telescope u… ▽ More Brown dwarfs with well-determined ages, luminosities, and masses provide rare but valuable tests of low-temperature atmospheric and evolutionary models. We present the discovery and dynamical mass measurement of a substellar companion to HD 47127, an old ($\approx$7-10 Gyr) G5 main sequence star with a mass similar to the Sun. Radial velocities of the host star with the Harlan J. Smith Telescope uncovered a low-amplitude acceleration of 1.93 $\pm$ 0.08 m s$^{-1}$ yr$^{-1}$ based on 20 years of monitoring. We subsequently recovered a faint ($ΔH$=13.14 $\pm$ 0.15 mag) co-moving companion at 1.95$''$ (52 AU) with follow-up Keck/NIRC2 adaptive optics imaging. The radial acceleration of HD 47127 together with its tangential acceleration from Hipparcos and Gaia EDR3 astrometry provide a direct measurement of the three-dimensional acceleration vector of the host star, enabling a dynamical mass constraint for HD 47127 B (67.5-177 $M_\mathrm{Jup}$ at 95% confidence) despite the small fractional orbital coverage of the observations. The absolute $H$-band magnitude of HD 47127 B is fainter than the benchmark T dwarfs HD 19467 B and Gl 229 B but brighter than Gl 758 B and HD 4113 C, suggesting a late-T spectral type. Altogether the mass limits for HD 47127 B from its dynamical mass and the substellar boundary imply a range of 67-78 $M_\mathrm{Jup}$ assuming it is single, although a preference for high masses of $\approx$100 $M_\mathrm{Jup}$ from dynamical constraints hints at the possibility that HD 47127 B could itself be a binary pair of brown dwarfs or that another massive companion resides closer in. Regardless, HD 47127 B will be an excellent target for more refined orbital and atmospheric characterization in the future. △ Less

Submitted 3 May, 2021; originally announced May 2021.

Comments: Accepted to ApJ Letters