skip to main content


Title: An advanced multipole model for (216) Kleopatra triple system
Aims. To interpret adaptive-optics observations of (216) Kleopatra, we need to describe an evolution of multiple moons orbiting an extremely irregular body and include their mutual interactions. Such orbits are generally non-Keplerian and orbital elements are not constants. Methods. Consequently, we used a modified N -body integrator, which was significantly extended to include the multipole expansion of the gravitational field up to the order ℓ = 10. Its convergence was verified against the ‘brute-force’ algorithm. We computed the coefficients C ℓm , S ℓm for Kleopatra’s shape, assuming a constant bulk density. For Solar System applications, it was also necessary to implement a variable distance and geometry of observations. Our χ 2 metric then accounts for the absolute astrometry, the relative astrometry (second moon with respect to the first), angular velocities, and silhouettes, constraining the pole orientation. This allowed us to derive the orbital elements of Kleopatra’s two moons. Results. Using both archival astrometric data and new VLT/SPHERE observations (ESO LP 199.C-0074), we were able to identify the true periods of the moons, P 1 = (1.822359 ± 0.004156) d, P 2 = (2.745820 ± 0.004820) d. They orbit very close to the 3:2 mean-motion resonance, but their osculating eccentricities are too small compared to other perturbations (multipole, mutual), meaning that regular librations of the critical argument are not present. The resulting mass of Kleopatra, m 1 = (1.49 ± 0.16) × 10 −12 M ⊙ or 2.97 × 10 18 kg, is significantly lower than previously thought. An implication explained in the accompanying paper is that (216) Kleopatra is a critically rotating body.  more » « less
Award ID(s):
1743015
NSF-PAR ID:
10295804
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; more » ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; « less
Date Published:
Journal Name:
Astronomy & Astrophysics
Volume:
653
ISSN:
0004-6361
Page Range / eLocation ID:
A56
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Aims. The orbit of the outer satellite Alexhelios of (216) Kleopatra is already constrained by adaptive-optics astrometry obtained with the VLT/SPHERE instrument. However, there is also a preceding occultation event in 1980 attributed to this satellite. Here, we try to link all observations, spanning 1980–2018, because the nominal orbit exhibits an unexplained shift by + 60° in the true longitude. Methods. Using both a periodogram analysis and an ℓ = 10 multipole model suitable for the motion of mutually interacting moons about the irregular body, we confirmed that it is not possible to adjust the respective osculating period P 2 . Instead, we were forced to use a model with tidal dissipation (and increasing orbital periods) to explain the shift. We also analysed light curves spanning 1977–2021, and searched for the expected spin deceleration of Kleopatra. Results. According to our best-fit model, the observed period rate is Ṗ 2 = (1.8 ± 0.1) × 10 −8 d d −1 and the corresponding time-lag Δ t 2 = 42 s of tides, for the assumed value of the Love number k 2 = 0.3. This is the first detection of tidal evolution for moons orbiting 100 km asteroids. The corresponding dissipation factor Q is comparable with that of other terrestrial bodies, albeit at a higher loading frequency 2| ω − n |. We also predict a secular evolution of the inner moon, Ṗ 1 = 5.0 × 10 −8 , as well as a spin deceleration of Kleopatra, Ṗ 0 = 1.9 × 10 −12 . In alternative models, with moons captured in the 3:2 mean-motion resonance or more massive moons, the respective values of Δ t 2 are a factor of between two and three lower. Future astrometric observations using direct imaging or occultations should allow us to distinguish between these models, which is important for our understanding of the internal structure and mechanical properties of (216) Kleopatra. 
    more » « less
  2. null (Ed.)
    Context. The recent estimates of the 3D shape of the M/Xe-type triple asteroid system (216) Kleopatra indicated a density of ~5 g cm −3 , which is by far the highest for a small Solar System body. Such a high density implies a high metal content as well as a low porosity which is not easy to reconcile with its peculiar “dumbbell” shape. Aims. Given the unprecedented angular resolution of the VLT/SPHERE/ZIMPOL camera, here, we aim to constrain the mass (via the characterization of the orbits of the moons) and the shape of (216) Kleopatra with high accuracy, hence its density. Methods. We combined our new VLT/SPHERE observations of (216) Kleopatra recorded during two apparitions in 2017 and 2018 with archival data from the W. M. Keck Observatory, as well as lightcurve, occultation, and delay-Doppler images, to derive a model of its 3D shape using two different algorithms (ADAM, MPCD). Furthermore, an N -body dynamical model allowed us to retrieve the orbital elements of the two moons as explained in the accompanying paper. Results. The shape of (216) Kleopatra is very close to an equilibrium dumbbell figure with two lobes and a thick neck. Its volume equivalent diameter (118.75 ± 1.40) km and mass (2.97 ± 0.32) × 10 18 kg (i.e., 56% lower than previously reported) imply a bulk density of (3.38 ± 0.50) g cm −3 . Such a low density for a supposedly metal-rich body indicates a substantial porosity within the primary. This porous structure along with its near equilibrium shape is compatible with a formation scenario including a giant impact followed by reaccumulation. (216) Kleopatra’s current rotation period and dumbbell shape imply that it is in a critically rotating state. The low effective gravity along the equator of the body, together with the equatorial orbits of the moons and possibly rubble-pile structure, opens the possibility that the moons formed via mass shedding. Conclusions. (216) Kleopatra is a puzzling multiple system due to the unique characteristics of the primary. This system certainly deserves particular attention in the future, with the Extremely Large Telescopes and possibly a dedicated space mission, to decipher its entire formation history. 
    more » « less
  3. 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
  4. Abstract

    HR 8799 is a young A5/F0 star hosting four directly imaged giant planets at wide separations (∼16–78 au), which are undergoing orbital motion and have been continuously monitored with adaptive optics imaging since their discovery over a decade ago. We present a dynamical mass of HR 8799 using 130 epochs of relative astrometry of its planets, which include both published measurements and new medium-band 3.1μm observations that we acquired with NIRC2 at Keck Observatory. For the purpose of measuring the host-star mass, each orbiting planet is treated as a massless particle and is fit with a Keplerian orbit using Markov chain Monte Carlo. We then use a Bayesian framework to combine each independent total mass measurement into a cumulative dynamical mass using all four planets. The dynamical mass of HR 8799 is1.470.17+0.12Massuming a uniform stellar mass prior, or1.460.15+0.11Mwith a weakly informative prior based on spectroscopy. There is a strong covariance between the planets’ eccentricities and the total system mass; when the constraint is limited to low-eccentricity solutions ofe< 0.1, which are motivated by dynamical stability, our mass measurement improves to1.430.07+0.06M. Our dynamical mass and other fundamental measured parameters of HR 8799 together with Modules for Experiments in Stellar Astrophysics Isochrones and Stellar Tracks grids yields a bulk metallicity most consistent with [Fe/H] ∼ −0.25–0.00 dex and an age of 10–23 Myr for the system. This implies hot-start masses of 2.7–4.9MJupfor HR 8799 b and 4.1–7.0MJupfor HR 8799 c, d, and e, assuming they formed at the same time as the host star.

     
    more » « less
  5. Abstract

    Exoplanet systems with multiple transiting planets are natural laboratories for testing planetary astrophysics. One such system is HD 191939 (TOI 1339), a bright (V= 9) and Sun-like (G9V) star, which TESS found to host three transiting planets (b, c, and d). The planets have periods of 9, 29, and 38 days each with similar sizes from 3 to 3.4R. To further characterize the system, we measured the radial velocity (RV) of HD 191939 over 415 days with Keck/HIRES and APF/Levy. We find thatMb= 10.4 ± 0.9MandMc= 7.2 ± 1.4M, which are low compared to most known planets of comparable radii. The RVs yield only an upper limit onMd(<5.8Mat 2σ). The RVs further reveal a fourth planet (e) with a minimum mass of 0.34 ± 0.01MJupand an orbital period of 101.4 ± 0.4 days. Despite its nontransiting geometry, secular interactions between planet e and the inner transiting planets indicate that planet e is coplanar with the transiting planets (Δi< 10°). We identify a second high-mass planet (f) with 95% confidence intervals on mass between 2 and 11MJupand period between 1700 and 7200 days, based on a joint analysis of RVs and astrometry from Gaia and Hipparcos. As a bright star hosting multiple planets with well-measured masses, HD 191939 presents many options for comparative planetary astronomy, including characterization with JWST.

     
    more » « less