skip to main content

Title: The Orbits and Dynamical Masses of the Castor System

Castor is a system of six stars in which the two brighter objects, Castor A and B, revolve around each other every ∼450 yr and are both short-period spectroscopic binaries. They are attended by the more distant Castor C, which is also a binary. Here we report interferometric observations with the Center for High Angular Resolution Astronomy (CHARA) array that spatially resolve the companions in Castor A and B for the first time. We complement these observations with new radial velocity measurements of A and B spanning 30 yr, with the Hipparcos intermediate data, and with existing astrometric observations of the visual AB pair obtained over the past three centuries. We perform a joint orbital solution to solve simultaneously for the three-dimensional orbits of Castor A and B as well as the AB orbit. We find that they are far from being coplanar: the orbit of A is nearly at right angles (92°) relative to the wide orbit, and that of B is inclined about 59° compared to AB. We determine the dynamical masses of the four stars in Castor A and B to a precision better than 1%. We also determine the radii of the primary stars of both subsystems from their angular diameters measured with the CHARA array, and use them together with stellar evolution models to infer an age for the system of 290 Myr. The new knowledge of the orbits enables us to measure the slow motion of Castor C as well, which may assist future studies of the dynamical evolution of this remarkable sextuple system.

more » « less
Award ID(s):
1909165 1506540 2034336
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ; ; ; ; ;
Publisher / Repository:
DOI PREFIX: 10.3847
Date Published:
Journal Name:
The Astrophysical Journal
Page Range / eLocation ID:
Article No. 8
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Some evolved binaries, namely post–asymptotic giant branch (AGB) binaries, are surrounded by stable and massive circumbinary disks similar to protoplanetary disks found around young stars. Around 10% of these disks are transition disks: they have a large inner cavity in the dust. Previous interferometric measurements and modeling have ruled out these cavities being formed by dust sublimation and suggested that they are due to massive circumbinary planets that trap dust in the disk and produce the observed depletion of refractory elements on the surfaces of the post-AGB stars. In this study, we test an alternative scenario in which the large cavities could be due to dynamical truncation from the inner binary. We performed near-infrared interferometric observations with the CHARA Array on the archetype of such a transition disk around a post-AGB binary: AC Her. We detect the companion at ten epochs over 4 yr and determine the three-dimensional orbit using these astrometric measurements in combination with a radial velocity time series. This is the first astrometric orbit constructed for a post-AGB binary system. We derive the best-fit orbit with a semimajor axis of 2.01 ± 0.01 mas (2.83 ± 0.08 au), inclination (142.9 ± 1.1)°, and longitude of the ascending node (155.1 ± 1.8)°. We find that the theoretical dynamical truncation and dust sublimation radii are at least ∼3× smaller than the observed inner disk radius (∼21.5 mas or 30 au). This strengthens the hypothesis that the origin of the cavity is due to the presence of a circumbinary planet.

    more » « less
  2. Context. Stellar evolution models are highly dependent on accurate mass estimates, especially for highly massive stars in the early stages of stellar evolution. The most direct method for obtaining model-independent stellar masses is derivation from the orbit of close binaries. Aims. Our aim was to derive the first astrometric plus radial velocity orbit solution for the single-lined spectroscopic binary star MWC 166 A, based on near-infrared interferometry over multiple epochs and ∼100 archival radial velocity measurements, and to derive fundamental stellar parameters from this orbit. A supplementary aim was to model the circumstellar activity in the system from K band spectral lines. Methods. The data used include interferometric observations from the VLTI instruments GRAVITY and PIONIER, as well as the MIRC-X instrument at the CHARA Array. We geometrically modelled the dust continuum to derive relative astrometry at 13 epochs, determine the orbital elements, and constrain individual stellar parameters at five different age estimates. We used the continuum models as a base to examine differential phases, visibilities, and closure phases over the Br γ and He  I emission lines in order to characterise the nature of the circumstellar emission. Results. Our orbit solution suggests a period of P  = 367.7 ± 0.1 d, approximately twice as long as found with previous radial velocity orbit fits. We derive a semi-major axis of 2.61 ± 0.04 au at d  = 990 ± 50 pc, an eccentricity of 0.498 ± 0.001, and an orbital inclination of 53.6 ± 0.3°. This allowed the component masses to be constrained to M 1  = 12.2 ± 2.2  M ⊙ and M 2  = 4.9 ± 0.5  M ⊙ . The line-emitting gas was found to be localised around the primary and is spatially resolved on scales of ∼11 stellar radii, where the spatial displacement between the line wings is consistent with a rotating disc. Conclusions. The large spatial extent and stable rotation axis orientation measured for the Br γ and He  I line emission are inconsistent with an origin in magnetospheric accretion or boundary-layer accretion, but indicate an ionised inner gas disc around this Herbig Be star. We observe line variability that could be explained either with generic line variability in a Herbig star disc or V/R variations in a decretion disc scenario. We have also constrained the age of the system, with relative flux ratios suggesting an age of ∼(7 ± 2)×10 5 yr, consistent with the system being composed of a main-sequence primary and a secondary still contracting towards the main-sequence stage. 
    more » « less
  3. null (Ed.)
    Context. Dynamical models of Solar System evolution have suggested that the so-called P- and D-type volatile-rich asteroids formed in the outer Solar System beyond Neptune’s orbit and may be genetically related to the Jupiter Trojans, comets, and small Kuiper belt objects (KBOs). Indeed, the spectral properties of P- and D-type asteroids resemble that of anhydrous cometary dust. Aims. We aim to gain insights into the above classes of bodies by characterizing the internal structure of a large P- and D-type asteroid. Methods. We report high-angular-resolution imaging observations of the P-type asteroid (87) Sylvia with the Very Large Telescope Spectro-Polarimetric High-contrast Exoplanet REsearch (SPHERE) instrument. These images were used to reconstruct the 3D shape of Sylvia. Our images together with those obtained in the past with large ground-based telescopes were used to study the dynamics of its two satellites. We also modeled Sylvia’s thermal evolution. Results. The shape of Sylvia appears flattened and elongated (a/b ~1.45; a/c ~1.84). We derive a volume-equivalent diameter of 271 ± 5 km and a low density of 1378 ± 45 kg m −3 . The two satellites orbit Sylvia on circular, equatorial orbits. The oblateness of Sylvia should imply a detectable nodal precession which contrasts with the fully-Keplerian dynamics of its two satellites. This reveals an inhomogeneous internal structure, suggesting that Sylvia is differentiated. Conclusions. Sylvia’s low density and differentiated interior can be explained by partial melting and mass redistribution through water percolation. The outer shell should be composed of material similar to interplanetary dust particles (IDPs) and the core should be similar to aqueously altered IDPs or carbonaceous chondrite meteorites such as the Tagish Lake meteorite. Numerical simulations of the thermal evolution of Sylvia show that for a body of such a size, partial melting was unavoidable due to the decay of long-lived radionuclides. In addition, we show that bodies as small as 130–150 km in diameter should have followed a similar thermal evolution, while smaller objects, such as comets and the KBO Arrokoth, must have remained pristine, which is in agreement with in situ observations of these bodies. NASA Lucy mission target (617) Patroclus (diameter ≈140 km) may, however, be differentiated. 
    more » « less
  4. Context . The Gl 486 system consists of a very nearby, relatively bright, weakly active M3.5 V star at just 8 pc with a warm transiting rocky planet of about 1.3 R ⊕ and 3.0 M ⊕ . It is ideal for both transmission and emission spectroscopy and for testing interior models of telluric planets. Aims . To prepare for future studies, we aim to thoroughly characterise the planetary system with new accurate and precise data collected with state-of-the-art photometers from space and spectrometers and interferometers from the ground. Methods . We collected light curves of seven new transits observed with the CHEOPS space mission and new radial velocities obtained with MAROON-X at the 8.1 m Gemini North telescope and CARMENES at the 3.5 m Calar Alto telescope, together with previously published spectroscopic and photometric data from the two spectrographs and TESS. We also performed near-infrared interferometric observations with the CHARA Array and new photometric monitoring with a suite of smaller telescopes (AstroLAB, LCOGT, OSN, TJO). This extraordinary and rich data set was the input for our comprehensive analysis. Results . From interferometry, we measure a limb-darkened disc angular size of the star Gl 486 at θ LDD = 0.390 ± 0.018 mas. Together with a corrected Gaia EDR3 parallax, we obtain a stellar radius R * = 0.339 ± 0.015 R ⊕ . We also measure a stellar rotation period at P rot = 49.9 ± 5.5 days, an upper limit to its XUV (5-920 A) flux informed by new Hubble /STIS data, and, for the first time, a variety of element abundances (Fe, Mg, Si, V, Sr, Zr, Rb) and C/O ratio. Moreover, we imposed restrictive constraints on the presence of additional components, either stellar or sub-stellar, in the system. With the input stellar parameters and the radial-velocity and transit data, we determine the radius and mass of the planet Gl 486 b at R p = 1.343 −0.062 +0.063 R ⊕ and M p = 3.00 −0.12 +0.13 M ⊕ , with relative uncertainties of the planet radius and mass of 4.7% and 4.2%, respectively. From the planet parameters and the stellar element abundances, we infer the most probable models of planet internal structure and composition, which are consistent with a relatively small metallic core with respect to the Earth, a deep silicate mantle, and a thin volatile upper layer. With all these ingredients, we outline prospects for Gl 486 b atmospheric studies, especially with forthcoming James Webb Space Telescope ( Webb ) observations. 
    more » « less
  5. Abstract We present the visual orbits of four spectroscopic binary stars, HD 61859, HD 89822, HD 109510, and HD 191692, using long baseline interferometry with the CHARA Array. We also obtained new radial velocities from echelle spectra using the APO 3.5 m, CTIO 1.5 m, and Fairborn Observatory 2.0 m telescopes. By combining the astrometric and spectroscopic observations, we solve for the full, three-dimensional orbits and determine the stellar masses to 1%–12% uncertainty and distances to 0.4%–6% uncertainty. We then estimate the effective temperature and radius of each component star through Doppler tomography and spectral energy distribution analyses. We found masses of 1.4–3.5 M ⊙ , radii of 1.5–4.7 R ⊙ , and temperatures of 6400–10,300 K. We then compare the observed stellar parameters to the predictions of the stellar evolution models, but found that only one of our systems fits well with the evolutionary models. 
    more » « less