EXTREME AO OBSERVATIONS OF TWO TRIPLE
ASTEROID SYSTEMS WITH SPHERE

Bin Yang, Z. Wahhaj, L. Beauvalet, F. Marchis, C. Dumas, M. Marsset, E.

Nielsen, F. Vachier

► To cite this version:

Bin Yang, Z. Wahhaj, L. Beauvalet, F. Marchis, C. Dumas, et al. EXTREME AO OBSERVATIONS OF TWO TRIPLE ASTEROID SYSTEMS WITH SPHERE. The Astrophysical Journal Letters, 2016, 820 (2), pp.L35. 10.3847/2041-8205/820/2/L35 . hal-02320061

HAL Id: hal-02320061

https://hal.science/hal-02320061v1

Submitted on 7 Nov 2025

**HAL** is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
EXTREME AO OBSERVATIONS OF TWO TRIPLE ASTEROID SYSTEMS WITH SPHERE

B. YANG<fnref n="1" />,<fnref n="2" />, Z. WAHHAJ<fnref n="1" />, L. BEAUVALET<fnref n="3" />, F. MARCHIS<fnref n="4" />, C. DUMAS<fnref n="1" />,<fnref n="5" />, M. MARSSET<fnref n="1" />, E. L. NIELSEN<fnref n="4" />,<fnref n="6" />, AND F. VACHIER<fnref n="7" />

Received 2015 November 25; accepted 2016 March 11; published 2016 March 29

ABSTRACT

We present the discovery of a new satellite of asteroid (130) Elektra—S/2014 (130) 1—in differential imaging and in integral field spectroscopy data over multiple epochs obtained with Spectro-Polarimetric High-contrast Exoplanet Research/Very Large Telescope. This new (second) moonlet of Elektra is about 2 km across, on an eccentric orbit, and about 500 km away from the primary. For a comparative study, we also observed another triple asteroid system, (93) Minerva. For both systems, component-resolved reflectance spectra of the satellites and primary were obtained simultaneously. No significant spectral difference was observed between the satellites and the primary for either triple system. We find that the moonlets in both systems are more likely to have been created by sub-disruptive impacts as opposed to having been captured.

Key words: infrared: planetary systems – minor planets, asteroids: individual (130 Elektra, 93 Minerva)

1. INTRODUCTION

Asteroids are the relics of the building blocks that formed the terrestrial planets in the early days of the solar system. Asteroids with satellites are of particular importance because their formation mechanisms, accretional and collisional processes, are critical in planet formation and evolution. In addition, multiple asteroids provide otherwise unattainable information about the intrinsic properties of the system: composition, interior structure, and evolutionary processes (Margot et al. 2015). Recently, a second-generation extreme- AO instrument, Spectro-Polarimetric High-contrast Exoplanet Research instrument (SPHERE; Beuzit et al. 2008), has been commissioned at one of European Southern Observatory’s 8 m Very Large Telescopes (VLT) at Cerro Paranal in Chile. SPHERE, with its unprecedented sensitivity and spatial resolution, offers a fresh opportunity for detecting and characterizing multiple asteroid systems.

More than 100 binary or triple asteroids have been detected. However, how and under what conditions these systems formed is poorly understood. Component-resolved spectroscopy of the primary and the moonlet obtained simultaneously is the key to understanding their formation. If the moonlets resulted from an impact, then they would be of the same composition as the primary. Whereas, if the multiple system originated from capture during a three-body encounter, then the composition of the satellites would differ from that of the primary. Most of the available component-resolved spectroscopic studies find no significant differences between the primary and secondary, for example (22) Kalliope (Laver et al. 2009), (90) Antiope (Marchis et al. 2011), (379) Huenna (DeMeo et al. 2011), and (809) Lundia (Birlan et al. 2014). However, a recent work on two triple asteroid systems, (45) Eugenia and (87) Sylvia, yields rather surprising results. Using SINFONI/VLT, Marchis et al. (2013a) found that the spectra of the moonlets are noticeably redder than their primaries. The difference in the spectral slopes of the individual components for Eugenia and Sylvia is not well understood. One possible

explanation is that the primary and the satellites may have experienced different surface alteration processes (such as space weathering effects; see Lantz et al. 2013) due to the difference in their dynamical ages. Alternatively, satellites may not result from impact but originate, instead, from a capture during a three-body encounter (Funato et al. 2004). To investigate these possibilities, we performed direct imaging and integral field spectroscopy of two multiple asteroids (i.e., (93) Minerva and (130) Elektra) using SPHERE.

(93) Minerva is the fifth triple asteroid discovered in the Main Belt (Marchis et al. 2009). The primary, located in the mid-asteroid belt ($r = 2.75$ au), is a large and dark asteroid ($D = 147 \pm 2$ km, $p_v = 0.068 \pm 0.003$; Usui et al. 2011), and its spectrum is classified as a C-type (Lazzaro et al. 2004; DeMeo et al. 2009). (130) Elektra was one of the four binary systems with eccentric mutual orbits (Marchis et al. 2008). The primary is also a large and dark main belt asteroid ($D = 197 \pm 20$ km, $p_v = 0.064 \pm 0.013$; Marchis et al. 2012). Its optical spectrum is classified as a G- or Ch-type (Tholen & Barucci 1989; Bus & Binzel 2002). Both Minerva and Elektra are low-albedo primitive asteroids, as are Eugenia and Sylvia. If the space weathering effect is responsible for the spectral difference between the primary and the two satellites seen in the Eugenia and Sylvia systems, the same process may have also modified the surface properties of the Minerva and Elektra systems. If so, we would expect the spectra of the moonlets to be different from those of the primaries.

2. OBSERVATION AND DATA REDUCTION

Observations were carried out with the two infrared subsystems of SPHERE/VLT simultaneously, i.e., the infrared differential imager and spectrograph (IRDIS; Dohlen et al. 2008) and the integral field spectrograph (IFS; Claudi et al. 2008). Low-resolution spectra were obtained with the IFS from 0.95 to 1.65 µm (Y- to H-band, spectral resolution ~30), while simultaneous images were obtained in two narrow ($\Delta\lambda = 0.1$ µm) bands K1 and K2 ($\lambda = 2.11$ µm and 2.25 µm)

<footnote n="1">¹European Southern Observatory, Santiago, Chile; byang@eso.org</footnote>

<footnote n="2">²Yunnan Astronomical Observatory, Chinese Academy of Sciences, Kunming, China</footnote>

<footnote n="3">³National Observatory, Rio de Janeiro, Brazil</footnote>

<footnote n="4">⁴Carl Sagan Center at the SETI Institute, Mountain View, CA, USA</footnote>

<footnote n="5">⁵TMT International Observatory, Pasadena, CA, USA</footnote>

<footnote n="6">⁶Kavli Institute for Particle Astrophysics and Cosmology, Stanford University, Stanford, CA 94305, USA</footnote>

<footnote n="7">⁷Institut de Mécanique Céleste et de Calcul des Éphémérides, Paris, France</footnote>
with IRDIS. For (130) Elektra, we observed in the field-stabilized mode, where the sky remains fixed with respect to the detector. (93) Minerva was observed in the pupil-tracking mode (Liu 2004; Marois et al. 2006), where the sky rotates, but the pupil stays fixed with respect to the detector. The fields of view (FOVs) of IRDIS and IFS are 11" × 12"/5 and 1"/73 × 1"/73, respectively, while the pixel scale is 0"/0123 for both.

We used the SPHERE consortium's pipeline (Pavlov et al. 2008) for data reduction. See Mesa et al. (2015) and Vigan et al. (2015) for further details. IRDIS images were dark subtracted and flat fielded and bad pixels were removed. The two K-band channels were combined to improve the signal-to-noise ratio (S/N). The IFS data were dark subtracted, bad pixel treated, flat fielded, and wavelength calibrated. The data were then resampled into a cube of 39 images of 3.3% band width (Δλ/λ) over the spectral range and with a scale of 0"/0074 per spaxel. Isolated bright and dark pixels were found in the IFS data, which is due to imperfections in the wavelength calibration procedure. The IDL routine acre (by Marc Buie) was used to remove these via interpolation between neighboring pixels.

Both systems were observed between UT 2014 December 6 and 9 as part of SPHERE science verification. Additional observations on Elektra were obtained on UT 2014 December 30 and 31, using director's discretionary time. At each epoch, we also observed a nearby solar analog (G2V) star to enable effective removal of telluric absorption artifacts and the solar color gradient. This star also served as a point-spread function (PSF) estimate for each spectral channel.

### 3. ADDITIONAL IMAGE PROCESSING

The known moonlet of (130) Elektra is not far from the primary and is nearly buried in the primary's halo. If this halo is not removed carefully, it could contaminate the photometry of its moonlet. Below, we describe the technique used to minimize the contamination from the primary. A two-dimensional Gaussian was fitted to the primary to determine its optical center. The primary's halo was then removed beyond ~0"/15 from the primary using a process described in Wahhaj et al. (2013). In this method, local medians over arcs (centered on the primary) of length 15-40 pixels are removed from each pixel. The arc lengths used depended on the spatial scale of background structure, and were at least three times the FWHM of any field source.

One unavoidable side effect of this process is that the flux of a moonlet is partially subtracted along with the halo. To estimate this signal loss, we performed simulations on the science data set. We inserted Gaussian point sources into the reduced IFS images for wavelength. The surface of the moonlets of neither system were resolved by SPHERE, and their PSFs, which depended on the AO performance, had FWHMs close to the diffraction limit. The FWHM of the simulated sources at each spectral channel was determined by PSF fitting of the standard star that was observed right before or after the asteroid observations. We then defined annuli centered on the primary, extending 10 pixels radially inward and outward from each satellite. In each annulus, we randomly selected 200 positions and inserted two types of point sources at each position. A strong source with 10 times the estimated flux of the satellite was used to estimate the systematic flux loss. A weaker source with a flux similar to that of the satellite was used to estimate the random error of our photometric measurements. After removing the

primary's halo, we measured the brightness of the recovered sources using two methods: aperture photometry and PSF fitting. Our analysis showed that the PSF fitting method generally yields smaller photometric uncertainties. However, when the moonlets were too faint or exhibit an irregular shape due to high airmass or poor AO correction, the PSF fitting method failed to emulate the PSF shape. In these cases, aperture photometry was used instead.

### 4. DISCOVERY OF NEW SATELLITE S/2014 (130) 1

As a result of this process (see Figure 1), a new moonlet, S/2014 (130) 1, was detected both in the IFS and IRDIS data sets over several epochs. The strongest signal of the new moonlet was detected on UT December 9 when it was 0"/36 from the primary. On December 31, the new inner satellite was only detected in the IFS data. The weakened signal of the inner satellite was due to the less favorable geometry where the asteroid system was farther away from the Sun and the Earth. The relative positions of the two satellites from the primary are listed in Table 1. The diameters of the satellites, estimated by calculating the integrated flux ratio of the primary to each companion, are also listed in Table 1. The effective diameter of the primary of 197 ± 20 km is adopted from Marchis et al. (2012). Assuming that the satellites have the same albedo as the primary, the estimated diameter of the new moon, S/2014 (130) 1, is 2.0 ± 1.5 km, and the estimated diameter for S/2003 (130) 1 is 6.0 ± 1.5 km.

#### 4.1. Numerical Modeling and Fitting

We obtained the orbital parameters of the two satellites of (130) Elektra using the numerical code ODIN (Beauvalet et al. 2013), which simultaneously integrates the heliocentric motion of the system, and the motion of the individual components around their barycenter. ODIN has been used to study the Pluto system and the triple systems, such as (45) Eugenia and (87) Sylvia (Beauvalet & Marchis 2014). Using ODIN to fit Elektra's astrometric data and assuming that the masses of the satellites are negligible, we obtained satisfactory orbital solutions for the two moonlets; see Figure 2(a). For satellite S/2003(130) 1, we also included the previous observations from Marchis et al. (2008) and compared our results with those obtained by GENOID (Vachier et al. 2012). We found that the two models provided very similar results. The orbital periods for the inner and the outer moonlets are 1.256 ± 0.003 days and 5.287 ± 0.001 days, respectively. The residuals for the outer moonlet in $α * cos(δ)$ and $δ$ are: 20 mas and 25 mas, where $α$ is right ascension and $δ$ is declination. For the inner satellite, S/2014 (130) 1, the residuals in R.A. and decl. are 11 mas and 14 mas, respectively. We used the random hold-out method described in Marchis et al. (2010) to estimate the uncertainties of the orbital elements. We randomly removed three observations from the data set and fit the remaining observations to obtain an orbital solution. We repeated this process 100 times and adopted $3σ$ uncertainty for each orbital element. We attempted to fit the gravitational harmonic parameter $J_2$, which is the main indicator of the mass distribution of the primary. We adopted an initial value of 0.13, based on axial ratios from shape model estimates (Marchis et al. 2008) and assumed a homogeneous interior. However, we found reasonable orbital solutions for a wide range of $J_2$ values from 0 to 0.13. More astrometric
<figure><img src="image_2.png" /><figcaption>**Figure 1.** Reduced IRDIS and IFS images of (130) Elektra, the triple asteroid system. The IFS images are the median combination of 31 of the 39 spectral channels. The channels severely affected by the atmospheric absorption were not included. The IFS images are rotated by 100°7 clockwise to align with the corresponding IRIDS data with PA = 0. The pixel intensities within 0.′′25 of the primary have been divided by 1500 to make the faint satellites more easily viewable. A new satellite S/2014 (130) S1 was detected in both the IRDIS and IFS data sets in most epochs. This inner moon detection was weaker in the December 31 data and was only identified in the IFS data.</figcaption></figure>

observations, especially of the inner moonlet, are needed to constrain the $J_2$ of the Elektra system.

## 4.2. Stability

We checked the stability of the orbits of the both satellites of Elektra using ODIN, the accuracy of which has been tested using other well-observed systems by simulating them for >100 years. Our main goal is to verify whether the orbits of the two moonlets are stable between 2006 and 2014. We performed two different simulations over 20 years. Neither of these simulations included the Yarkovsky, YORP, or tidal effects, which are only needed for very long-term integrations. In one simulation, the masses of the satellites are negligible. In the other simulation, the two satellites have masses of 7 × 10¹² kg and 2 × 10¹⁴ kg, respectively, assuming the same density as the primary of 1.7 g cm⁻³ (Marchis et al. 2008). Neither simulation showed significant instability over time. The semimajor axis for S/2014 (130) 1 has a small oscillation spanning less than 100 m. For S/2003 (130) 1, semimajor axis oscillates less than 1.4 km.

# 5. COMPONENT-RESOLVED SPECTROSCOPY

## 5.1. The (130) Elektra System

The Elektra system was observed on four nights. However, we have reliable spectra for the inner moonlet from December 9 to 30 only. On UT December 9, the primary's halo varied significantly

in intensity over small spatial scales near the inner moonlet. Hence, we used an arc length of 15 pixels for halo removal. The solar analog star, HD 20926, had an average PSF FWHM over 39 spectral channels of 5.3 pixels. The outer moonlet is roughly 0.′′6 away from the primary where the primary's halo is more uniform and so an arc of length 30 pixels was used to remove the background flux. On UT December 30, given a different AO performance, an arc of length 40 pixels was used for both moonlets for halo removal. The solar analog star, HD 19315, had an average PSF FWHM of 4.9 pixels. The spectra of the Elektra system are presented in Figure 3, which have been corrected for the systematic flux loss estimated from our simulations. We calculated the spectral slope, $s'$, the fractional intensity change per 1000 Å, using a weighted linear regression. We performed students' t-test and calculated *P* to compare the spectra of the satellites and that of the primary, where *P* is used as an indicator to determine whether the spectra of the individual component are significantly different. We found *P* = 0.56 and 0.82 on UT December 9 and 30, respectively, for the inner moonlet and *P* = 1.0 and 0.56 for the outer moonlet. Since the *p*-values are much larger than the significance level of 0.05, we conclude that there is no significant spectral difference between the moonlets and the primary.

## 5.2. The (93) Minerva System

As shown in Figure 4(a), the two moons are well separated from the primary in the IFS data and the halo of the primary has
<table><caption>**Table 1**

Astrometric and Photometric Measurements for the (93) Minerva and (130) Elektra Triple Systems</caption><thead><tr><th>Object</th><th>Date (UT)</th><th>Time (UT)</th><th>δX<sup>a</sup> (arcsec)</th><th>δY<sup>a</sup> (arcsec)</th><th>δD<sup>b</sup> (arcsec)</th><th>Δm<sup>c</sup> (mag)</th><th>D<sub>s</sub><sup>d</sup> (km)</th></tr></thead><tbody><tr><td>S/2003 (130) 1</td><td>2014 Dec 06</td><td>03:25:22</td><td>-0.351 ± 0.002</td><td>-0.625 ± 0.004</td><td>0.717</td><td>8.25</td><td>4.4</td></tr><tr><td>S/2003 (130) 1</td><td>2014 Dec 06</td><td>03:36:20</td><td>-0.339 ± 0.002</td><td>-0.624 ± 0.002</td><td>0.710</td><td>8.13</td><td>4.7</td></tr><tr><td>S/2003 (130) 1</td><td>2014 Dec 06</td><td>03:47:43</td><td>-0.327 ± 0.003</td><td>-0.623 ± 0.004</td><td>0.703</td><td>8.09</td><td>4.8</td></tr><tr><td>S/2003 (130) 1</td><td>2014 Dec 09</td><td>01:33:19</td><td>0.317 ± 0.001</td><td>0.555 ± 0.003</td><td>0.639</td><td>7.65</td><td>5.8</td></tr><tr><td>S/2003 (130) 1</td><td>2014 Dec 09</td><td>01:36:13</td><td>0.315 ± 0.001</td><td>0.557 ± 0.003</td><td>0.640</td><td>7.67</td><td>5.8</td></tr><tr><td>S/2003 (130) 1</td><td>2014 Dec 09</td><td>01:39:11</td><td>0.313 ± 0.001</td><td>0.557 ± 0.003</td><td>0.639</td><td>7.70</td><td>5.7</td></tr><tr><td>S/2003 (130) 1</td><td>2014 Dec 09</td><td>01:42:05</td><td>0.309 ± 0.001</td><td>0.554 ± 0.005</td><td>0.635</td><td>7.65</td><td>5.8</td></tr><tr><td>S/2003 (130) 1</td><td>2014 Dec 09</td><td>01:45:00</td><td>0.307 ± 0.001</td><td>0.557 ± 0.002</td><td>0.637</td><td>7.62</td><td>5.9</td></tr><tr><td>S/2003 (130) 1</td><td>2014 Dec 09</td><td>01:47:56</td><td>0.305 ± 0.001</td><td>0.558 ± 0.001</td><td>0.636</td><td>7.58</td><td>6.0</td></tr><tr><td>S/2003 (130) 1</td><td>2014 Dec 09</td><td>01:50:50</td><td>0.302 ± 0.001</td><td>0.554 ± 0.004</td><td>0.631</td><td>7.61</td><td>5.9</td></tr><tr><td>S/2003 (130) 1</td><td>2014 Dec 09</td><td>01:53:45</td><td>0.301 ± 0.001</td><td>0.558 ± 0.001</td><td>0.634</td><td>7.55</td><td>6.1</td></tr><tr><td>S/2003 (130) 1</td><td>2014 Dec 09</td><td>01:56:40</td><td>0.298 ± 0.001</td><td>0.556 ± 0.003</td><td>0.631</td><td>7.56</td><td>6.0</td></tr><tr><td>S/2003 (130) 1</td><td>2014 Dec 09</td><td>01:59:36</td><td>0.299 ± 0.003</td><td>0.565 ± 0.006</td><td>0.639</td><td>7.58</td><td>6.0</td></tr><tr><td>S/2003 (130) 1</td><td>2014 Dec 30</td><td>01:05:42</td><td>0.491 ± 0.002</td><td>0.456 ± 0.003</td><td>0.670</td><td>7.53</td><td>6.1</td></tr><tr><td>S/2003 (130) 1</td><td>2014 Dec 30</td><td>01:29:57</td><td>0.475 ± 0.002</td><td>0.458 ± 0.004</td><td>0.660</td><td>7.56</td><td>6.1</td></tr><tr><td>S/2003 (130) 1</td><td>2014 Dec 30</td><td>01:39:00</td><td>0.467 ± 0.002</td><td>0.458 ± 0.003</td><td>0.654</td><td>7.59</td><td>6.0</td></tr><tr><td>S/2003 (130) 1</td><td>2014 Dec 31</td><td>03:45:33</td><td>-0.683 ± 0.002</td><td>0.120 ± 0.003</td><td>0.693</td><td>7.83</td><td>5.4</td></tr><tr><td>S/2003 (130) 1</td><td>2014 Dec 31</td><td>03:54:35</td><td>-0.688 ± 0.001</td><td>0.117 ± 0.005</td><td>0.698</td><td>7.72</td><td>5.6</td></tr><tr><td>S/2003 (130) 1</td><td>2014 Dec 31</td><td>04:03:35</td><td>-0.691 ± 0.002</td><td>0.112 ± 0.003</td><td>0.700</td><td>7.65</td><td>5.8</td></tr><tr><td>S/2014 (130) 1</td><td>2014 Dec 06</td><td>03:25:22</td><td>-0.374 ± 0.005</td><td>-0.062 ± 0.006</td><td>0.379</td><td>10.26</td><td>1.7</td></tr><tr><td>S/2014 (130) 1</td><td>2014 Dec 06</td><td>03:36:20</td><td>-0.369 ± 0.005</td><td>-0.073 ± 0.005</td><td>0.376</td><td>10.21</td><td>1.8</td></tr><tr><td>S/2014 (130) 1</td><td>2014 Dec 06</td><td>03:47:43</td><td>-0.359 ± 0.007</td><td>-0.080 ± 0.008</td><td>0.367</td><td>10.98</td><td>1.3</td></tr><tr><td>S/2014 (130) 1</td><td>2014 Dec 09</td><td>01:33:19</td><td>0.361 ± 0.005</td><td>-0.002 ± 0.005</td><td>0.361</td><td>9.74</td><td>2.2</td></tr><tr><td>S/2014 (130) 1</td><td>2014 Dec 09</td><td>01:36:13</td><td>0.361 ± 0.004</td><td>0.001 ± 0.005</td><td>0.361</td><td>9.74</td><td>2.2</td></tr><tr><td>S/2014 (130) 1</td><td>2014 Dec 09</td><td>01:39:11</td><td>0.363 ± 0.003</td><td>0.004 ± 0.004</td><td>0.363</td><td>9.77</td><td>2.2</td></tr><tr><td>S/2014 (130) 1</td><td>2014 Dec 09</td><td>01:42:05</td><td>0.362 ± 0.003</td><td>0.003 ± 0.006</td><td>0.362</td><td>9.66</td><td>2.3</td></tr><tr><td>S/2014 (130) 1</td><td>2014 Dec 09</td><td>01:45:00</td><td>0.365 ± 0.003</td><td>0.010 ± 0.003</td><td>0.365</td><td>9.66</td><td>2.3</td></tr><tr><td>S/2014 (130) 1</td><td>2014 Dec 09</td><td>01:47:56</td><td>0.366 ± 0.003</td><td>0.011 ± 0.003</td><td>0.366</td><td>9.76</td><td>2.2</td></tr><tr><td>S/2014 (130) 1</td><td>2014 Dec 09</td><td>01:50:50</td><td>0.365 ± 0.003</td><td>0.013 ± 0.004</td><td>0.365</td><td>9.75</td><td>2.2</td></tr><tr><td>S/2014 (130) 1</td><td>2014 Dec 09</td><td>01:53:45</td><td>0.367 ± 0.003</td><td>0.018 ± 0.003</td><td>0.368</td><td>9.77</td><td>2.2</td></tr><tr><td>S/2014 (130) 1</td><td>2014 Dec 09</td><td>01:56:40</td><td>0.370 ± 0.003</td><td>0.020 ± 0.004</td><td>0.371</td><td>9.78</td><td>2.2</td></tr><tr><td>S/2014 (130) 1</td><td>2014 Dec 09</td><td>01:59:36</td><td>0.375 ± 0.004</td><td>0.032 ± 0.007</td><td>0.377</td><td>9.90</td><td>2.1</td></tr><tr><td>S/2014 (130) 1</td><td>2014 Dec 30</td><td>01:05:42</td><td>-0.322 ± 0.004</td><td>-0.103 ± 0.005</td><td>0.338</td><td>9.70</td><td>2.3</td></tr><tr><td>S/2014 (130) 1</td><td>2014 Dec 30</td><td>01:29:57</td><td>-0.320 ± 0.006</td><td>-0.130 ± 0.007</td><td>0.337</td><td>9.59</td><td>2.4</td></tr><tr><td>S/2014 (130) 1</td><td>2014 Dec 30</td><td>01:39:00</td><td>-0.307 ± 0.003</td><td>-0.141 ± 0.004</td><td>0.338</td><td>9.25</td><td>2.8</td></tr><tr><td>S/2014 (130) 1</td><td>2014 Dec 31</td><td>03:45:33</td><td>-0.315 ± 0.011</td><td>-0.017 ± 0.011</td><td>0.316</td><td>10.46</td><td>1.6</td></tr><tr><td>S/2014 (130) 1</td><td>2014 Dec 31</td><td>03:54:35</td><td>-0.316 ± 0.009</td><td>-0.018 ± 0.010</td><td>0.317</td><td>10.76</td><td>1.4</td></tr><tr><td>S/2014 (130) 1</td><td>2014 Dec 31</td><td>04:03:35</td><td>-0.316 ± 0.004</td><td>-0.034 ± 0.005</td><td>0.318</td><td>10.17</td><td>1.8</td></tr><tr><td>Aegis</td><td>2014 Dec 08</td><td>01:48:17</td><td>-0.323 ± 0.001</td><td>0.085 ± 0.002</td><td>0.333</td><td>7.25</td><td>4.3</td></tr><tr><td>Gorgoneion</td><td>2014 Dec 08</td><td>01:48:17</td><td>-0.037 ± 0.002</td><td>-0.231 ± 0.002</td><td>0.234</td><td>8.65</td><td>3.7</td></tr></tbody></table>

**Notes.**

<sup>a</sup> The relative Cartesian coordinates (δX and δY) of the satellite with respect to the optical center of the primary.

<sup>b</sup> Angular separation of the satellite from the primary.

<sup>c</sup> Integrated flux ratio between the primary and the satellite.

<sup>d</sup> Estimated diameter of the satellite using the integrated flux ratio of the primary and the satellite, assuming the same albedo.

been removed effectively using the same method described above. We adopted the effective diameter of the primary of 154 ± 6 km from Marchis et al. (2013b) and estimated the sizes of the two satellites using the ratios of the integrated fluxes between the satellites and the primary. Assuming the satellites share the same albedos as the primary, the estimated diameter for Aegis is 4.3 ± 1.5 km and the diameter of Gorgoneion is 3.7 ± 1.5 km, which are consistent with the previous diameter estimates of 3.6 ± 1.0 km and 3.2 ± 0.9 km in Marchis et al. (2013b). The large differences between the primary spectrum and the satellites' spectra, at wavelengths

between 1.3 and 1.5 µm, are a result of strong atmospheric absorption. We also performed a t-test for the Minerva system. We found no discernible difference between satellite and primary, both for Aegis (P = 0.64) and for Gorgoneion (P = 0.79).

## 6. DISCUSSION

The major challenge in obtaining reliable component-resolved spectroscopy of tight asteroid systems is the bright halo of the primary and the variable sky background. The outer satellites, i.e., S/2003 (130) 1 and I Aegis that are orbiting at
<figure><img src="image_3.png" /><figcaption>**Figure 2.** Left: the gray dashed lines show the best-fit orbit for S/2003 S1, as it evolved over many epochs, and the pink line is the best-fit orbit for S/2014 S1. The blue and green solid circles are the astrometric measurements of the two moonlets. The red and pink pluses are the simulated positions using the dynamical model ODIN. The uncertainties of the orbital elements are $3\sigma$ errors derived using the random hold-out method. Right: the eccentricities of the satellites that belong to the five known triple systems; $a_s$ is the semimajor axis of the satellite's orbit and $R_{\text{Hill}}$ is the Hill radius of the primary. References: Marchis et al. (2010; M10), Descamps et al. (2011; D11), Fang et al. (2012; F12), Marchis et al. (2013b; M13b).</figcaption></figure>

<figure><img src="image_4.png" /><figcaption>**Figure 3.** Relative reflectance spectra of the (130) Elektra system, normalized at 1.0 µm. Left: spectra of the two satellites and the primary obtained on UT December 9, using PSF fitting. Right: spectra of the two satellites and the primary obtained on UT December 30. The spectrum of the inner moon, obtained via aperture photometry. The spectra of the primary and the outer moon were obtained via PSF fitting. For comparison, we also plot the near-infrared spectrum of the primary obtained with the 3.0 m IRTF telescope by Takir & Emery (2012). The wavelength range affected by telluric absorption is shaded in the plot.</figcaption></figure>

13 and 8 primary radii and have detection $S/N > 10$, yield spectra that are not very sensitive to our image processing methods. In contrast, the spectral slopes of the inner moons (S/2014 (130) 1 and II Gorgoneion) are heavily influenced by the halo removal process because these two are closer to the primary and are faint in comparison to the background noise. As such, it is necessary to understand the impact of the halo removal methods at different wavelengths. Inserting simulated sources with known PSFs and fluxes into the real data set and processing the data in the same manner as the science is perhaps the most reliable way to correct photometric biases and estimate uncertainties.

Our simulations show large uncertainties in the spectral measurements. Therefore, we are not able to further investigate space weathering effects on the individual components of the triple systems. In terms of the satellites, we detect no absorption band near 1.0 µm, and the spectral slopes of these moonlets are flat (within a few percent). Our IFS observations

suggest that the components of both the Elektra and the Minerva systems are similar in composition and the four satellites are best classified as C-types. The rotational periods for (130) Elektra and (93) Minerva are 5.225 and 5.982 hr$^8$, respectively. Therefore, break-up from a high spin rate is unlikely to have led to formation these triple systems. Given the small separations of the satellites from their primaries and the large mass ratios between the moonlets and primaries, an erosive impact is the most likely explanation for both the Elektra and the Minerva triple systems.

As of now, the eccentricity of the satellite orbits for most multiple asteroid systems have not been estimated because of the difficulty in obtaining sufficient position measurements of the satellites. For the few systems for which we have enough observations to characterize the satellites' orbits, a large range

<footnote n="8">⁸ Taken from minor planet light curve parameters, A. W. Harris and B. D. Warner, http://cfa-www.harvard.edu/iau/lists/LightcurveDat.html.</footnote>
<figure><img src="image_5.png" /><figcaption>**Figure 4.** Left: reduced IFS image of the Minerva system, which is the median combination of 31 of the 39 spectral channels. The pixel intensities within 0.1" of the primary have been divided by 300 to show the two satellites clearly. Right: the relative reflectance spectra of Aegis and Gorgoneion were obtained via aperture photometry, normalized at 1.0 µm. The reflectance spectra of the individual components of (93) Minerva are presented in Figure 4(b). A comparison IRTF spectrum of the primary obtained by Marchis et al. (2013b) is shown in green.</figcaption></figure>

of eccentricities have been observed; see Figure 2(b). Our observations of (130) Elektra reveal a unique case: both moonlets are on eccentric orbits. We note that asteroid Elektra is elongated with *a/b* ~ 1.5. However, it is known that the elongation of the primary has no long-term effect on the eccentricity, the semimajor axis, or inclination (Murray & Dermott 2000). More detailed dynamical simulations of the orbital evolution of the asteroid triple systems are needed to understand the difference in eccentricity between the inner and the outer moon.

In summary, our SPHERE observations yield the following findings: (1) a new moonlet of (130) Elektra was detected and the estimated diameter of this moon is ~2 km, (2) both the inner and the outer moonlet of (130) Elektra are on orbits with non-negligible eccentricities, and (3) the reflectance spectra of the satellites of the (130) Elektra and the (93) Minerva systems are similar to those of their primaries, suggesting that both triple systems were created by disruptive impacts.

We thank Julien Milli and Arthur Vigan for helping with the SPHERE data reduction. Special thanks to David Jewitt for commenting on the paper. F.M. was supported in part by NASA cooperative agreement NNX14AJ80G.

## REFERENCES

* Beauvalet, L., & Marchis, F. 2014, Icar, 241, 13

* Beauvalet, L., Robert, V., Lainey, V., Arlot, J. E., & Colas, F. 2013, A&A, 553, A14

* Beuzit, J.-L., Feldt, M., Dohlen, K., et al. 2008, Proc. SPIE, 7014, 701418

* Birlan, M., Nedelcu, D. A., Popescu, M., et al. 2014, MNRAS, 437, 176

* Bus, S. J., & Binzel, R. P. 2002, Icar, 158, 146

* Claudi, R. U., Turatto, M., Gratton, R. G., et al. 2008, Proc. SPIE, 7014, 70143E

* DeMeo, F. E., Binzel, R. P., Slivan, S. M., & Bus, S. J. 2009, Icar, 202, 160

* DeMeo, F. E., Carry, B., Marchis, F., et al. 2011, Icar, 212, 677

* Descamps, P., Marchis, F., Berthier, J., et al. 2011, Icar, 211, 1022

* Dohlen, K., Langlois, M., Saisse, M., et al. 2008, Proc. SPIE, 7014, 70143L

* Fang, J., Margot, J.-L., & Rojo, P. 2012, AJ, 144, 70

* Funato, Y., Makino, J., Hut, P., Kokubo, E., & Kinoshita, D. 2004, Natur, 427, 518

* Lantz, C., Clark, B. E., Barucci, M. A., & Lauretta, D. S. 2013, A&A, 554, A138

* Laver, C., de Pater, I., Marchis, F., Ádámkovics, M., & Wong, M. H. 2009, Icar, 204, 574

* Lazzaro, D., Angeli, C. A., Carvano, J. M., et al. 2004, Icar, 172, 179

* Liu, M. C. 2004, Sci, 305, 1442

* Marchis, F., Descamps, P., Berthier, J., et al. 2008, Icar, 195, 295

* Marchis, F., Enriquez, J. E., Emery, J. P., et al. 2011, Icar, 213, 252

* Marchis, F., Enriquez, J. E., Emery, J. P., et al. 2012, Icar, 221, 1130

* Marchis, F., Lainey, V., Descamps, P., et al. 2010, Icar, 210, 635

* Marchis, F., Macomber, B., Berthier, J., Vachier, F., & Emery, J. P. 2009, IAUC, 9069, 1

* Marchis, F., Ruffio, J., Vachier, F., & Berthier, J. 2013a, AGUFM, 23, 1760

* Marchis, F., Vachier, F., Ñurech, J., et al. 2013b, Icar, 224, 178

* Margot, J.-L., Pravec, P., Taylor, P., Carry, B., & Jacobson, S. 2015, Asteroid Systems: Binaries, Triples, and Pairs, chapter in the book ASTEROIDS IV, in press

* Marois, C., Lafrenière, D., Doyon, R., Macintosh, B., & Nadeau, D. 2006, ApJ, 641, 556

* Mesa, D., Gratton, R., Zurlo, A., et al. 2015, A&A, 576, A121

* Murray, C. D., & Dermott, S. F. (ed.) 2000, Solar System Dynamics (Cambridge: Cambridge Univ. Press)

* Pavlov, A., Möller-Nilsson, O., Feldt, M., et al. 2008, Proc. SPIE, 7019, 701939

* Takir, D., & Emery, J. P. 2012, Icar, 219, 641

* Tholen, D. J., & Barucci, M. A. 1989, in Asteroids II (Tucson, AZ: Univ. Arizona Press), 298

* Usui, F., Kuroda, D., Müller, T. G., et al. 2011, PASJ, 63, 1117

* Vachier, F., Berthier, J., & Marchis, F. 2012, A&A, 543, A68

* Vigan, A., Gry, C., Salter, G., et al. 2015, MNRAS, 454, 129

* Wahhaj, Z., Liu, M. C., Biller, B., et al. 2013, ApJ, 779, 80