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

HAL
open science

# The Trans-Neptunian Object (84922) 2003 VS 2 through Stellar Occultations

Gustavo Benedetti-Rossi, P. Santos-Sanz, J. Ortiz, M. Assafin, B. Sicardy, N. Morales, R. Vieira-Martins, R. Duffard, F. Braga-Ribas, F. Rommel, et al.

► To cite this version:

Gustavo Benedetti-Rossi, P. Santos-Sanz, J. Ortiz, M. Assafin, B. Sicardy, et al. The Trans-Neptunian Object (84922) 2003 VS 2 through Stellar Occultations. The Astronomical Journal, 2019, 158 (4), pp.159. 10.3847/1538-3881/ab3b05. hal-03105635

HAL Id: hal-03105635

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

Submitted on 5 Sep 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_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_3.png" /></figure>

# The Trans-Neptunian Object (84922) 2003 VS₂ through Stellar Occultations

Gustavo Benedetti-Rossi<sup>1,2,3</sup>, P. Santos-Sanz<sup>4</sup>, J. L. Ortiz<sup>4</sup>, M. Assafin<sup>5</sup>, B. Sicardy<sup>3,6</sup>, N. Morales<sup>4</sup>, R. Vieira-Martins<sup>1,2,7</sup>, R. Duffard<sup>4</sup>, F. Braga-Ribas<sup>1,2,3,8</sup>, F. L. Rommel<sup>1,2</sup>, J. I. B. Camargo<sup>1,2</sup>, J. Desmars<sup>3</sup>, A. F. Colas<sup>7</sup>, F. Vachier<sup>7</sup>, A. Alvarez-Candal<sup>1</sup>, E. Fernández-Valenzuela<sup>9</sup>, L. Almenares<sup>10</sup>, R. Artola<sup>11</sup>, T.-P. Baum<sup>12</sup>, R. Behrend<sup>13</sup>, D. Bérard<sup>3</sup>, F. B. Bianco<sup>14,15,16,17</sup>, N. Brosch<sup>18</sup>, A. Ceretta<sup>19</sup>, C. A. Colazo<sup>11</sup>, A. R. Gomes-Junior<sup>2,20</sup>, V. D. Ivanov<sup>21</sup>, E. Jehin<sup>22</sup>, S. Kaspi<sup>18</sup>, J. Lecacheux<sup>3</sup>, A. Maury<sup>23</sup>, R. Melia<sup>11</sup>, S. Moindrot<sup>24</sup>, B. Morgado<sup>1,2</sup>, C. Opitom<sup>25</sup>, A. Peyrot<sup>12</sup>, J. Pollock<sup>26</sup>, A. Pratt<sup>27</sup>, S. Roland<sup>28</sup>, J. Spagnotto<sup>29</sup>, G. Tancredi<sup>10</sup>, J.-P. Teng<sup>12</sup>, P. Cacella<sup>30</sup>, M. Emilio<sup>1,31</sup>, F. Feys<sup>27</sup>, R. Gil-Hutton<sup>32</sup>, C. Jacques<sup>33</sup>, D. I. Machado<sup>34,35</sup>, M. Malacarne<sup>36</sup>, I. Manulis<sup>37</sup>, A. C. Milone<sup>38</sup>, G. Rojas<sup>39</sup>, and R. Sfair<sup>20</sup>

<sup>1</sup>Observatório Nacional (ON/MCTIC), Rua Gal. José Cristino, 77—Bairro Imperial de São Cristóvão, Rio de Janeiro, Brazil; gugabrossi@gmail.com, gustavo.benedetti-rossi@obspm.fr, gustavorossi@on.br

<sup>2</sup>Laboratório Interinstitucional de e-Astronomia (LIneA) & INCT do e-Universo, Rio de Janeiro, Brazil

<sup>3</sup> LESIA, Observatoire de Paris—Section Meudon, 5 Place Jules Janssen—92195 Meudon Cedex, France

<sup>4</sup>Instituto de Astrofísica de Andalucía, IAA-CSIC, Glorieta de la Astronomía s/n, E-18008 Granada, Spain

<sup>5</sup>Observatório do Valongo/UFRJ, Ladeira Pedro Antonio 43, Rio de Janeiro, 20080-090, Brazil

<sup>6</sup> PSL Research University, CNRS, Sorbonne Université, UPMC Univ. Paris 06, Univ. Paris Diderot, Sorbonne Paris Cité, France

<sup>7</sup>IMCCE/Observatoire de Paris, CNRS UMR F-8028, Paris, France

<sup>8</sup>Federal University of Technology—Paraná (UTFPR/DAFIS), Av. Sete de Setembro, 3165, CEP 80230-901, Curitiba/PR, Brazil

<sup>9</sup>Florida Space Institute, University of Central Florida, 12354 Research Parkway, Partnership 1, Suite 211, Orlando, FL, USA

<sup>10</sup>Departamento de Astronomia, Facultad de Ciencias, Montevideo, Uruguay

<sup>11</sup>Estación Astrofísica de Bosque Alegre, Cordoba, Argentina

<sup>12</sup>Observatoire Astronomique des Makes, Les Makes, F-97421—La Rivière, Ile de la Réunion, France

<sup>13</sup>Observatoire de Genève, CH-1290 Sauvernay, Switzerland

<sup>14</sup>University of Delaware, Department of Physics and Astronomy, USA

<sup>15</sup>University of Delaware, Joseph R. Biden Jr. School for Public Policy and Administration, USA

<sup>16</sup>University of Delaware, Data Science Institute, USA

<sup>17</sup>New York University, Center for Urban Science and Progress, USA

<sup>18</sup>School of Physics & Astronomy and Wise Observatory, Tel Aviv University, Israel

<sup>19</sup>Observatorio Astronómico Los Molinos, MEC, Uruguay

<sup>20</sup>UNESP—São Paulo State University, Grupo de Dinâmica Orbital e Planetologia, CEP 12516-410, Guaratinguetá, SP, Brazil

<sup>21</sup>ESO, Karl-Schwarzschild-Str. 2, D-85748 Garching bei München, Germany

<sup>22</sup>STAR Institute—University of Liège, Alleé du 6 Août 17, B-4000 Liège, Belgium

<sup>23</sup>San Pedro de Ataçama Celestial Explorations, Casilla 21, San Pedro de Atacama 1410000, Chile

<sup>24</sup>Observatoire de Puimichel, F-04700—Puimichel, France

<sup>25</sup>European Southern Observatory—Alonso de Cordova 3107, Vitacura, Santiago, Chile

<sup>26</sup>Physics and Astronomy Department, Appalachian State University, Boone, NC 28608, USA

<sup>27</sup>IOTA/ES—International Occultation Timing Association/European Section, Germany

<sup>28</sup>CURE, Universidad de la República, Uruguay

<sup>29</sup>Observatorio El Catalejo, Santa Rosa, La Pampa, Argentina

<sup>30</sup>DogsHeaven Observatory X87, Brasília, Brazil

<sup>31</sup>Universidade Estadual de Ponta Grossa, O.A.—DEGEO, Avenida Carlos Cavalcanti 4748, Ponta Grossa 84030-900, PR, Brazil

<sup>32</sup>Grupo de Ciencias Planetarias, Departamento de Geofísica y Astronomía, San Juan National University and CONICET, Av. José I. de la Roza 590(0), J5402DCS, San Juan, Argentina

<sup>33</sup>Southern Observatory for Near Earth Asteroids Research (SONEAR), Brazil

<sup>34</sup>Universidade Estadual do Oeste do Paraná (Unioeste), Avenida Tarquínio Joslin dos Santos 1300, Foz do Iguaçu, PR 85870-650, Brazil

<sup>35</sup>Polo Astronômico Casimiro Montenegro Filho/FPTI-BR, Avenida Tanredo Neves 6731, Foz do Iguaçu, PR 85867-900, Brazil

<sup>36</sup>Universidade Federal do Espírito Santo (UFES), Brazil

<sup>37</sup>Weizmann Institute of Science, Rehovot, Israel

<sup>38</sup>Astrophysics Division, National Institute for Space Research (INPE), Av. dos Astronautas 1758, São José dos Campos, SP, 12227-010, Brazil

<sup>39</sup>Observatório Astronômico—Universidade Federal de São Carlos (UFSCar), São Carlos, Brazil

*Received 2019 May 27; revised 2019 August 8; accepted 2019 August 11; published 2019 September 27*

## Abstract

We present results from three world-wide campaigns that resulted in the detections of two single-chord and one multi-chord stellar occultations by the plutino object (84922) 2003 VS₂. From the single-chord occultations in 2013 and 2014 we obtained accurate astrometric positions for the object, while from the multi-chord occultation on 2014 November 7, we obtained the parameters of the best-fitting ellipse to the limb of the body at the time of occultation. We also obtained short-term photometry data for the body in order to derive its rotational phase during the occultation. The rotational light curve present a peak-to-peak amplitude of 0.141 ± 0.009 mag. This allows us to reconstruct the 3D shape of the body, with principal semi-axes of *a* = 313.8 ± 7.1 km, *b* = 265.5<sup>+8.8</sup><sub>-9.8</sub> km, and *c* = 247.3<sup>+26.6</sup><sub>-43.6</sub> km, which is not consistent with a Jacobi triaxial equilibrium figure. The derived spherical volume equivalent diameter of 548.3<sup>+29.5</sup><sub>-44.6</sub> km is about 5% larger than the radiometric diameter of 2003 VS₂ derived from Herschel data of 523 ± 35 km, but still compatible with it within error bars. From those results we can also derive the geometric albedo (0.123<sup>+0.015</sup><sub>-0.014</sub>) and, under the assumption that the object is a Maclaurin spheroid, the density *ρ* = 1400<sup>+1000</sup><sub>-300</sub> for the plutino. The disappearances and reappearances of the star
during the occultations do not show any compelling evidence for a global atmosphere considering a pressure upper limit of about 1 microbar for a pure nitrogen atmosphere, nor secondary features (e.g., rings or satellite) around the main body.

*Unified Astronomy Thesaurus concepts: Trans-Neptunian objects (1705); Kuiper belt (893)*

# 1. Introduction

Trans-Neptunian objects (TNOs) are bodies that orbit the Sun with an orbital semimajor axis larger than that of Neptune (Jewitt et al. 2008). Due to their large distance to the Sun and low spatial density, those objects do not experience extensive differentiation. Consequently, TNOs are the least evolved bodies in the solar system, at least from the composition point of view. Hence, knowing their physical parameters such as size, shape, albedo, density, presence of atmosphere, rings, and their evolution, yield important information on the nature of material, physical conditions, and history of the primitive solar nebula (Morbidelli et al. 2008), and about the formation and evolution of our solar system (Lykawka & Mukai 2008; Parker 2015; Santos-Sanz et al. 2016). Additionally, the Kuiper belt provides the natural connection with the study of protoplanetary disks observed around other stars (Anglada et al. 2017).

Despite the fact that more than 25 yr have elapsed since this population was discovered (the first TNO discovered was (15760) Albion by Jewitt & Luu (1993), only Pluto and, more recently (in 2019 January 1), the small TNO 2014 MU<sub>69</sub> have been visited so far by a spacecraft, the NASA/New Horizons mission (Stern et al. 2019), making them unique bodies in the trans-Neptunian region with a more complete set of information. In summary, our knowledge of basic physical properties of the TNO population is still scarce and fragmentary, mainly due to the faintness and small angular sizes of these bodies as seen from Earth (Stansberry et al. 2008; Lellouch et al. 2013).

In this context, we have been conducting various kinds of Earth-based observations for nearly two decades to gather relevant physical information about TNOs (Lellouch et al. 2002; Ortiz et al. 2002, 2004, 2006, 2007, 2011, 2012, 2015, 2017; Belskaya et al. 2006; Duffard et al. 2008; Assafin et al. 2010, 2012; Thirouin et al. 2010; Sicardy et al. 2011; Braga-Ribas et al. 2013, 2014; Camargo et al. 2014, 2018; Benedetti-Rossi et al. 2016; Bérard et al. 2017; Dias-Oliveira et al. 2017; Leiva et al. 2017). In this framework, we predicted and detected three stellar occultations by the TNO (84922) 2003 VS2—VS2 hereafter—that we monitored through a collaboration with a large panel of amateurs. Among the three events, we obtained two single-chord detections in 2013 December and in 2014 March and one multi-chord detection in 2014 November. The latter observation was made from four well-separated sites, from which valuable physical information is derived, as described in this paper.

The object VS2 was discovered by the Near Earth Asteroid Tracking (NEAT)<fnref n="40" /> program on 2003 November 14. As a 2:3 resonant object with Neptune, it belongs to the plutino class (MPlanet Center Electronic Circular 2006; Gladman et al. 2008).<fnref n="41" /> Most of its physical properties were derived from thermal measurements using the space telescopes *Herschel* and *Spitzer* (Stansberry et al. 2008; Mommert et al. 2012; Santos-Sanz et al. 2017). Orbital parameters, absolute magnitude, *B*−*R* color, rotation period, photometric variation, taxonomy, diameter, and

geometric albedo, taken from previously published works, are listed in Table 1.

Here we present the results from those three stellar occultations by VS2, with the determination of size and shape of this body based on the multi-chord occultation of 2014 November 7. We also present a rotational light curve for VS2 derived from the 1.5 m telescope located at Sierra Nevada observatory in Granada (Spain), showing an amplitude that is smaller than reported previously in the literature.

We describe in Section 2 the circumstances of observation for the stellar occultations (predictions, updates, and the occultations themselves) and for the photometric runs to obtain the light curve. In Section 3 we describe our data reduction, analysis, and the results obtained from the stellar occultations. We discuss in Section 4 the results concerning the 3D shape of the body as well as other features by combining our occultation and photometric data. Finally, conclusions are presented in Section 5.

# 2. Observations

Observations performed in this work are separated into three groups: (2.1) observations performed to predict and refine the astrometric positions of the occulted stars and of the object for the occultations; (2.2) observations of the occultations; and (2.3) the observation runs that yielded VS2's rotational light curve.

## 2.1. Stellar Occultations Predictions and Updates

Stellar occultations require large efforts in order to determine accurate positions of the star and of the object's ephemeris. Despite the fact that the *Gaia* catalog—on its second data release, GDR2—presents positions and proper motions for stars with unprecedented accuracy (milliarcseconds levels; Gaia Collaboration et al. 2016a, 2016b, 2018; van Leeuwen et al. 2017), the object's ephemeris is usually determined at accuracies of a few hundred milliarcseconds, due to scarce observations. This is larger than the projected sizes of the objects in the sky plane, resulting in misses when attempting stellar occultation observations. Until data from large surveys, such as the Large Synoptic Survey Telescope (LSST),<fnref n="42" /> become available in a few years from now, regular observations of those small objects will still be needed for most of them in order to improve their orbital parameters. In fact, even with LSST data available, big efforts will still be needed to improve objects' ephemeris, as LSST will mostly cover the southern sky, leaving the objects in the northern hemisphere with big uncertainties. Note also that all those objects have long orbital periods (nearly two centuries for the closest TNOs) and that the time span of the observations are short (less than 20 yr) so the uncertainties in the orbital elements are very sensitive to new observations.

At the time of the three occultations reported here, we did not have the *Gaia* catalog in our hands so the candidate stars were identified in systematic surveys performed at the 2.2 m telescope of the European Southern Observatory (ESO) at La

<footnote n="40">⁴⁰ https://neat.jpl.nasa.gov/</footnote>

<footnote n="41">⁴¹ Another classification is given by M. Brown in http://web.gps.caltech.edu/~mbrown/dps.html.</footnote>

<footnote n="42">⁴² https://www.lsst.org/</footnote>
<table><caption>**Table 1**
Orbital and Physical Characteristics of the Plutino Object 2003 VS2</caption><thead><tr><th>a<br>(au)</th><th>q<br>(au)</th><th>i<br>(deg)</th><th>e</th><th>H<sub>v</sub><sup>*</sup><br>(mag)</th><th>B-R<br>(mag)</th><th>P<br>(hr)</th><th>Δm<br>(mag)</th><th>Taxon.</th><th>D<sup>**</sup><br>(km)</th><th>p<sub>V</sub><br>(%)</th><th>APmag</th></tr></thead><tbody><tr><td>39.25</td><td>36.43</td><td>14.80</td><td>0.07</td><td>4.130 ± 0.070</td><td>1.52 ± 0.03</td><td>7.4175285 ± 0.00001</td><td>0.224 ± 0.013</td><td>BB</td><td>523.0<sup>+35.1</sup><sub>-34.4</sub></td><td>14.7<sup>+6.3</sup><sub>-4.3</sub></td><td>19.95</td></tr></tbody></table>

**Note.** Orbital parameters: *a*: semimajor axis in astronomical units; *q*: perihelion distance in astronomical units; *i*: orbital inclination in degrees; *e*: eccentricity—from the JPL Small-Body Database Browser; *H<sub>v</sub>*: average absolute magnitude at the V band obtained from Alvarez-Candal et al. (2016) and private communication; *B-R* color from Sheppard (2007), Mommert et al. (2012), and references therein; *P*: preferred rotation period in hours from Santos-Sanz et al. (2017); Δ*m* [*mag*]: light curve peak-to-valley amplitude from Thirouin (2013); Taxon.: taxonomic color class from Perna et al. (2010), and references therein; *D*: area-equivalent diameter; *p<sub>V</sub>*: geometric albedo at the V band from Mommert et al. (2012) using *Herschel* and *Spitzer* data; and *APmag*: 2003 VS2's average apparent visual magnitude between 2013 and 2014 December from Horizons/JPL (calculated using *H<sub>v</sub>* = 4.2). * Sheppard (2007) and Mommert et al. (2012) obtained a value of *H<sub>v</sub>* = 4.110 ± 0.380; ** Stansberry et al. (2008) obtained a diameter of *D* = 725<sup>+199</sup><sub>-187</sub> km with *Spitzer* data only.

<table><caption>**Table 2**
J2000 Stars Coordinates, Proper Motions, Parallax, and Magnitudes</caption><thead><tr><th>Occultation Date</th><th>R.A. ICRS</th><th>errRA<br>(mas)</th><th>Decl. ICRS</th><th>errDEC<br>(mas)</th><th>pmRA<br>(mas yr<sup>-1</sup>)</th><th>pmDEC<br>(mas yr<sup>-1</sup>)</th><th>Plx<br>(mas)</th><th>G<br>(mag)</th><th>B<br>(mag)</th><th>V<br>(mag)</th><th>K<br>(mag)</th></tr></thead><tbody><tr><td>2013 Dec 12</td><td>04<sup>h</sup>36<sup>m</sup>19<sup>s</sup>830879</td><td>0.0502</td><td>+33° 55′ 55.45028</td><td>0.0434</td><td>0.878</td><td>-0.649</td><td>0.0733</td><td>16.0899</td><td>16.57</td><td>16.20</td><td>12.64</td></tr><tr><td>2014 Mar 4</td><td>04<sup>h</sup>32<sup>m</sup>04<sup>s</sup>463481</td><td>0.0548</td><td>+33° 24′ 41.95395</td><td>0.0422</td><td>2.925</td><td>-3.268</td><td>0.7682</td><td>15.9110</td><td>17.29</td><td>16.23</td><td>13.54</td></tr><tr><td>2014 Nov 7</td><td>04<sup>h</sup>48<sup>m</sup>32<sup>s</sup>138835</td><td>0.0406</td><td>+33° 58′ 36.40859</td><td>0.0303</td><td>7.404</td><td>-8.135</td><td>1.7609</td><td>15.1311</td><td>17.16</td><td>15.82</td><td>12.04</td></tr></tbody></table>

**Note.** R.A. and decl. ICRS: barycentric R.A. and decl. (ICRS) propagated to the occultation epoch; errRA and errDEC: standard error of R.A. (multiplied by cosDEC) and decl.; pmRA and pmDEC: proper motion in R.A. (multiplied by cosDEC) and in decl. direction; Plx: absolute stellar parallax; and G:magnitude. All of the values were obtained from GDR2 (Gaia Collaboration et al. 2016a, 2016b, 2018) and the B, V, and K:magnitudes were obtained from the NOMAD catalog (Zacharias et al. 2004).

Silla—IAU code 809, using the Wide Field Imager (WFI). The surveys yielded local astrometric catalogs for 5 Centaurs and 34 TNOs (plus Pluto and its moons) up to 2015, with stars with magnitudes as faint as R mag ~ 19, see details in Assafin et al. (2010, 2012) and Camargo et al. (2014). All three VS2 occultation candidate stars were then identified in the GDR2 catalog and their R.A., decl., proper motions, and G, B, V, and K magnitudes are presented in Table 2.

Astrometric updates of the candidate stars (and the TNO when possible) were performed using several telescopes close to the dates of events in Brazil at the Pico dos Dias Observatory—IAU code 874—(Perkin Elmer 1.6 m, Boller & Chivens 0.6 m, and ZEISS 0.6 m), in Spain at the Calar Alto Observatory—IAU code 493—(1.23 m and 2.2 m telescopes), at the Sierra Nevada Observatory—IAU code J86—(1.5 m telescope), and at the Observatorio de La Hitta—IAU code I95—(77 cm telescope), and in France at Pic du Midi—IAU code 586—(T100 cm telescope) so to reduce systematic errors caused by the reference catalogs. Astrometric positions obtained for VS2 with those observations are listed on Table 3. The final positions obtained for the star and the TNO before the occultations provided predictions with uncertainties as large as 50 mas, equivalent to about 1300 km when projected onto Earth's surface.

## 2.2. Stellar Occultations

For each of the three stellar occultations, an alert was triggered at several potential sites, resulting in data collected with a large diversity of instruments, see Tables 4–6. All sites used robust clock synchronizations, either by a global positioning system (GPS) or by setting up the clocks using Internet servers—network time Protocol (NTP), and acquisition times of each image were inserted on the image header or printed on individual video frames. Using any of those synchronization methods, the times

should not present absolute errors larger than 1 ms (Deeths & Brunette 2001). At all sites with favourable weather conditions, data were collected from about 10 minutes prior to the predicted occultation times until about 10 minutes after those times. No filters were used in any of the observations to maximize photon fluxes, and thus signal-to-noise-ratio (S/N).

Out of 23 stations distributed in nine countries (France, Israel, United Kingdom, Greece, Argentina, Uruguay, Chile, Brazil, and Bolivia), we obtained two single-chord events and one multi-chord occultation, for a total of seven positive detections. Note that the 2014 March occultation was detected by two telescopes at the same site in Israel. As such, they do not provide any constraint on the shape of the object, and we consider it as a single-chord detection. Figures 1–3 show the post-occultation, reconstructed maps for the three events. Note that in those figures we used a radius for VS2 of 564.8 km (determined from the multi-chord occultation), but we do not know the rotational phase for the two single-chord occultations, so the shadow path may be smaller than this value. Details of the data analysis for each occultation are given in Section 3.

## 2.3. Rotational Light Curve

In order to determine the rotational phase of VS2 at the time of the multi-chord occultation, we performed a photometric follow up of the object a few days after the event. We obtained a total of 97 images of VS2 over three nights—2014 November 15, 16, and 17—with the 2k × 2k CCD of the 1.5 m telescope at Sierra Nevada Observatory in Granada (Spain). The image scale of the instrument is 0.′′232 pixel⁻¹ with a field of view (FOV) of 7.92 × 7.92. All of the images were acquired without filter and in 2 × 2 binning mode to maximize S/N. The integration time was 300 s throughout the three nights with a Moon illumination of 36% during the first night and 18% during
<table><caption>**Table 3**

Astrometric Positions for VS2</caption><thead><tr><th colspan="3">Date</th><th colspan="3">R.A. ICRS</th><th colspan="3">Decl. ICRS</th><th rowspan="2">Site<br>[IAU Code]</th></tr><tr><th>Year</th><th>Month</th><th>Day</th><th>Deg</th><th>Min</th><th>Sec</th><th>Deg</th><th>Min</th><th>Sec</th></tr></thead><tbody><tr><td>2014</td><td>Feb</td><td>01.868243</td><td>68</td><td>03</td><td>09.4710</td><td>33</td><td>36</td><td>12.4600</td><td>J86</td></tr><tr><td>2014</td><td>Feb</td><td>01.882180</td><td>68</td><td>03</td><td>08.9115</td><td>33</td><td>36</td><td>12.1370</td><td>J86</td></tr><tr><td>2014</td><td>Feb</td><td>01.886832</td><td>68</td><td>03</td><td>08.8125</td><td>33</td><td>36</td><td>11.9780</td><td>J86</td></tr><tr><td>2014</td><td>Feb</td><td>01.891829</td><td>68</td><td>03</td><td>08.5650</td><td>33</td><td>36</td><td>11.8790</td><td>J86</td></tr><tr><td>2014</td><td>Feb</td><td>01.905770</td><td>68</td><td>03</td><td>08.1405</td><td>33</td><td>36</td><td>11.4950</td><td>J86</td></tr><tr><td>2014</td><td>Feb</td><td>03.000210</td><td>68</td><td>02</td><td>31.7055</td><td>33</td><td>35</td><td>44.4050</td><td>J86</td></tr><tr><td>2014</td><td>Feb</td><td>03.004858</td><td>68</td><td>02</td><td>31.5690</td><td>33</td><td>35</td><td>44.2810</td><td>J86</td></tr><tr><td>2014</td><td>Feb</td><td>03.009507</td><td>68</td><td>02</td><td>31.4220</td><td>33</td><td>35</td><td>44.1640</td><td>J86</td></tr><tr><td>2014</td><td>Feb</td><td>03.014154</td><td>68</td><td>02</td><td>31.2825</td><td>33</td><td>35</td><td>44.0510</td><td>J86</td></tr><tr><td>2014</td><td>Feb</td><td>03.018799</td><td>68</td><td>20</td><td>31.0935</td><td>33</td><td>35</td><td>43.9310</td><td>J86</td></tr><tr><td>2014</td><td>Feb</td><td>03.023446</td><td>68</td><td>02</td><td>30.9375</td><td>33</td><td>35</td><td>43.7860</td><td>J86</td></tr><tr><td>2014</td><td>Feb</td><td>03.028092</td><td>68</td><td>02</td><td>30.8565</td><td>33</td><td>35</td><td>43.7010</td><td>J86</td></tr><tr><td>2014</td><td>Feb</td><td>03.032739</td><td>68</td><td>02</td><td>30.6495</td><td>33</td><td>35</td><td>43.5740</td><td>J86</td></tr><tr><td>2014</td><td>Feb</td><td>03.037384</td><td>68</td><td>02</td><td>30.5475</td><td>33</td><td>35</td><td>43.4220</td><td>J86</td></tr><tr><td>2014</td><td>Feb</td><td>03.042029</td><td>68</td><td>02</td><td>30.4290</td><td>33</td><td>35</td><td>43.3400</td><td>J86</td></tr><tr><td>2014</td><td>Feb</td><td>19.920970</td><td>67</td><td>58</td><td>05.7000</td><td>33</td><td>29</td><td>06.1000</td><td>586</td></tr><tr><td>2014</td><td>Feb</td><td>19.920971</td><td>67</td><td>58</td><td>05.8260</td><td>33</td><td>29</td><td>06.2120</td><td>586</td></tr><tr><td>2014</td><td>Feb</td><td>19.924510</td><td>67</td><td>58</td><td>05.7000</td><td>33</td><td>29</td><td>06.0000</td><td>586</td></tr><tr><td>2014</td><td>Feb</td><td>19.924510</td><td>67</td><td>58</td><td>05.7705</td><td>33</td><td>29</td><td>06.1050</td><td>586</td></tr><tr><td>2014</td><td>Feb</td><td>19.928070</td><td>67</td><td>58</td><td>05.7000</td><td>33</td><td>29</td><td>06.0000</td><td>586</td></tr><tr><td>2014</td><td>Feb</td><td>19.928074</td><td>67</td><td>58</td><td>05.8995</td><td>33</td><td>29</td><td>06.0400</td><td>586</td></tr><tr><td>2014</td><td>Feb</td><td>19.931640</td><td>67</td><td>58</td><td>05.7000</td><td>33</td><td>29</td><td>05.8000</td><td>586</td></tr><tr><td>2014</td><td>Feb</td><td>19.931640</td><td>67</td><td>58</td><td>05.8170</td><td>33</td><td>29</td><td>05.9790</td><td>586</td></tr><tr><td>2014</td><td>Feb</td><td>19.935207</td><td>67</td><td>58</td><td>05.8155</td><td>33</td><td>29</td><td>05.9290</td><td>586</td></tr><tr><td>2014</td><td>Feb</td><td>19.935210</td><td>67</td><td>58</td><td>05.7000</td><td>33</td><td>29</td><td>05.8000</td><td>586</td></tr><tr><td>2014</td><td>Feb</td><td>19.938776</td><td>67</td><td>58</td><td>05.8560</td><td>33</td><td>29</td><td>05.8680</td><td>586</td></tr><tr><td>2014</td><td>Feb</td><td>19.938780</td><td>67</td><td>58</td><td>05.7000</td><td>33</td><td>29</td><td>05.7000</td><td>586</td></tr><tr><td>2014</td><td>Feb</td><td>19.942340</td><td>67</td><td>58</td><td>05.7000</td><td>33</td><td>29</td><td>05.7000</td><td>586</td></tr><tr><td>2014</td><td>Feb</td><td>19.942344</td><td>67</td><td>58</td><td>05.8230</td><td>33</td><td>29</td><td>05.7750</td><td>586</td></tr><tr><td>2014</td><td>Feb</td><td>19.945910</td><td>67</td><td>58</td><td>05.8500</td><td>33</td><td>29</td><td>05.6000</td><td>586</td></tr><tr><td>2014</td><td>Feb</td><td>19.945910</td><td>67</td><td>58</td><td>05.8620</td><td>33</td><td>29</td><td>05.7090</td><td>586</td></tr><tr><td>2014</td><td>Feb</td><td>19.949476</td><td>67</td><td>58</td><td>05.8080</td><td>33</td><td>29</td><td>05.5710</td><td>586</td></tr><tr><td>2014</td><td>Feb</td><td>19.949480</td><td>67</td><td>58</td><td>05.7000</td><td>33</td><td>29</td><td>05.4000</td><td>586</td></tr><tr><td>2014</td><td>Feb</td><td>19.953040</td><td>67</td><td>58</td><td>05.8500</td><td>33</td><td>29</td><td>05.4000</td><td>586</td></tr><tr><td>2014</td><td>Feb</td><td>19.953041</td><td>67</td><td>58</td><td>05.9040</td><td>33</td><td>29</td><td>05.4890</td><td>586</td></tr><tr><td>2014</td><td>Feb</td><td>19.956608</td><td>67</td><td>58</td><td>05.8545</td><td>33</td><td>29</td><td>05.3620</td><td>586</td></tr><tr><td>2014</td><td>Feb</td><td>19.956610</td><td>67</td><td>58</td><td>05.7000</td><td>33</td><td>29</td><td>05.2000</td><td>586</td></tr><tr><td>2014</td><td>Feb</td><td>19.963670</td><td>67</td><td>58</td><td>05.8500</td><td>33</td><td>29</td><td>05.1000</td><td>586</td></tr><tr><td>2014</td><td>Feb</td><td>19.967230</td><td>67</td><td>58</td><td>05.7000</td><td>33</td><td>29</td><td>05.0000</td><td>586</td></tr><tr><td>2014</td><td>Feb</td><td>19.970800</td><td>67</td><td>58</td><td>05.7000</td><td>33</td><td>29</td><td>05.1000</td><td>586</td></tr><tr><td>2014</td><td>Feb</td><td>23.832322</td><td>67</td><td>58</td><td>25.3935</td><td>33</td><td>27</td><td>42.0770</td><td>195</td></tr><tr><td>2014</td><td>Feb</td><td>23.833071</td><td>67</td><td>58</td><td>24.9570</td><td>33</td><td>27</td><td>41.8390</td><td>195</td></tr><tr><td>2014</td><td>Feb</td><td>23.835317</td><td>67</td><td>58</td><td>25.4850</td><td>33</td><td>27</td><td>41.7380</td><td>195</td></tr><tr><td>2014</td><td>Feb</td><td>23.836792</td><td>67</td><td>58</td><td>24.8685</td><td>33</td><td>27</td><td>41.7520</td><td>195</td></tr><tr><td>2014</td><td>Feb</td><td>23.840530</td><td>67</td><td>58</td><td>25.2300</td><td>33</td><td>27</td><td>41.8910</td><td>195</td></tr><tr><td>2014</td><td>Feb</td><td>23.841275</td><td>67</td><td>58</td><td>25.2330</td><td>33</td><td>27</td><td>42.0330</td><td>195</td></tr><tr><td>2014</td><td>Sep</td><td>24.148195</td><td>72</td><td>47</td><td>13.9545</td><td>33</td><td>53</td><td>41.9090</td><td>493</td></tr><tr><td>2014</td><td>Sep</td><td>24.157750</td><td>72</td><td>47</td><td>13.8450</td><td>33</td><td>53</td><td>42.0140</td><td>493</td></tr><tr><td>2014</td><td>Sep</td><td>24.167348</td><td>72</td><td>47</td><td>13.6620</td><td>33</td><td>53</td><td>42.1710</td><td>493</td></tr><tr><td>2014</td><td>Sep</td><td>24.172112</td><td>72</td><td>47</td><td>13.5675</td><td>33</td><td>53</td><td>42.2450</td><td>493</td></tr><tr><td>2014</td><td>Sep</td><td>24.181729</td><td>72</td><td>47</td><td>13.4475</td><td>33</td><td>53</td><td>42.3420</td><td>493</td></tr><tr><td>2014</td><td>Sep</td><td>24.186538</td><td>72</td><td>47</td><td>13.3410</td><td>33</td><td>53</td><td>42.4430</td><td>493</td></tr><tr><td>2014</td><td>Sep</td><td>24.196181</td><td>72</td><td>47</td><td>13.1865</td><td>33</td><td>53</td><td>42.5840</td><td>493</td></tr><tr><td>2014</td><td>Sep</td><td>24.200971</td><td>72</td><td>47</td><td>13.1190</td><td>33</td><td>53</td><td>42.6250</td><td>493</td></tr><tr><td>2014</td><td>Sep</td><td>24.205804</td><td>72</td><td>47</td><td>13.0065</td><td>33</td><td>53</td><td>42.7340</td><td>493</td></tr><tr><td>2014</td><td>Sep</td><td>24.210612</td><td>72</td><td>47</td><td>12.9735</td><td>33</td><td>53</td><td>42.7860</td><td>493</td></tr><tr><td>2014</td><td>Sep</td><td>24.215402</td><td>72</td><td>47</td><td>12.9045</td><td>33</td><td>53</td><td>42.8490</td><td>493</td></tr></tbody></table>

Note. Positions obtained with the ESO 2p2 telescope—IAU code 809—under mission 090.C-0118(A) are published in Camargo et al. (2014).
<table><caption>**Table 4**
Observation Details of the Single-chord Stellar Occultations</caption><thead><tr><th colspan="7">2013 Dec 12</th></tr><tr><th>Site Name<br>(Location)<br>[IAU Code]</th><th>Longitude (E)<br>Latitude (N)<br>Altitude (m)</th><th>Telescope Aperture<br>Detector/Instrument</th><th>Exposure Time<br>Cycle Time<br>(s)</th><th>Observers</th><th>Detection</th></tr></thead><tbody><tr><td>Les Makes Observatory<br>(La Réunion, France)<br>[181]</td><td>55° 24′ 36″0<br>-21° 11′ 58″4<br>992</td><td>60 cm<br>Raptor Merlin 246</td><td>0.75000<br>0.75031</td><td>J. Lecacheux<br>A. Peyrot, T.-P. Baum<br>J.-P. Teng</td><td>Positive</td></tr></tbody><thead><tr><th colspan="7">2014 Mar 4</th></tr><tr><th>Site Name<br>(Location)<br>[IAU Code]</th><th>Longitude (E)<br>Latitude (N)<br>Altitude (m)</th><th>Telescope Aperture<br>Detector/Instrument</th><th>Exposure Time<br>Cycle Time<br>(s)</th><th>Observers</th><th>Detection</th></tr></thead><tbody><tr><td>Wise Observatory<br>(Mitzpe Ramon, Israel)<br>[097]</td><td>34° 45′ 43″60<br>30° 35′ 50″38<br>884</td><td>100 cm<br>Roper PVCAM<br>PI 1300B/LN</td><td>4<br>7</td><td>S. Kaspi</td><td>Positive</td></tr><tr><td>Wise Observatory<br>(Mitzpe Ramon, Israel)<br>[097]</td><td>34° 45′ 43″92<br>30° 35′ 48″23<br>865</td><td>70 cm<br>FLI<br>PL16801</td><td>5<br>6.5</td><td>N. Brosch</td><td>Positive</td></tr><tr><td>West Park Observatory<br>(Leeds, UK)<br>[Z92]</td><td>358° 23′ 31″98<br>53° 50′ 15″42<br>114</td><td>20 cm<br>WATEC<br>WAT-910HX/RC (x256)</td><td>5.08<br>5.08<br>(video)</td><td>A. Pratt</td><td>Negative</td></tr><tr><td>Puimichel Observatory<br>(Puimichel, France)<br>...</td><td>06° 01′ 15″5<br>43° 58′ 48″7<br>724</td><td>104 cm<br>Raptor Merlin 246</td><td>0.25000<br>0.25031</td><td>J. Lecacheux<br>S. Moindrot</td><td>Negative</td></tr><tr><td>Weizmann<br>(Kraar, Israel)<br>...</td><td>34° 48′ 45″76<br>31° 54′ 29″1<br>107</td><td>41 cm<br>SBIG<br>ST-8XME CCD</td><td>...</td><td>I. Manulis</td><td>Overcast</td></tr><tr><td>Sasteria Observatory<br>(Crete, Greece)<br>...</td><td>25° 58′ 16″37<br>35° 04′ 11″12<br>395</td><td>25 cm<br>N/A</td><td>...</td><td>F. Feys</td><td>Overcast</td></tr></tbody></table>

Note. FLI stands for Finger Lakes instruments, PI for Princeton instruments, PL for pro line, and N/A for not available.

the third night. The average seeing during the three nights was 3"/3. We pointed the telescope at the same coordinates on the three nights in order to use the same reference stars and thus minimize systematic photometric errors.

Standard bias and flat-field corrections were applied on science images. Specific routines written in interactive data language (IDL) were developed to perform the aperture photometry of all the chosen reference stars and VS2. We tried different apertures for the target, calibration stars, and sky ring annulus in order to maximize the S/N on the object for each night and to minimize the dispersion of the photometry.

The flux of VS2 relative to the comparison stars is finally obtained versus the Julian date, accounting for light travel times. The procedures we used were identical to those detailed in Fernández-Valenzuela et al. (2016, 2017). We folded the final photometric data with the well-known rotational period for VS2 obtained in Santos-Sanz et al. (2017), $P = 7.4175285 \pm 0.00001$ hr. From this rotational light curve, we deduce that VS2 was near one of its absolute brightness maxima at the time of the 2014 November 7 stellar occultation (Figure 4), which implies that the object occulted the star when its apparent surface area was near its maximum. These folded data are also fitted with a second-order Fourier function in order to obtain the peak-to-valley amplitude of the rotational light curve, which turns out to

be $0.141 \pm 0.009$ mag (as listed in Table 10). The decrease in the peak-to-peak amplitude, compared to the value obtained by Thirouin (2013), is probably due to geometric effects (i.e., a change in the aspect angle with respect to the previous amplitude estimations). The epoch for zero phase was chosen to be near the time of occultation (2014 November 7 at 04:00:00 UTC).

### 3. Stellar Occultation Data Analysis

Some observations involved in the multi-chord event recorded only the integer part of the second in each image header. It was then necessary to retrieve the fractional part of the second for the mid-exposure time of each image by performing a linear fit to the set, see details in Sicardy et al. (2011) and Dias-Oliveira et al. (2017). From the linear fit we then retrieve the fraction of second for each image with an internal accuracy that depends on the square root of the number of images used, which is less than 0.1 s in practice.

Data obtained in video format from West Park Observatory was converted to flexible image transport system (fits) format using the same procedure as described in Benedetti-Rossi et al. (2016) and Buie & Keller (2016).

The two single-chord events were analyzed to obtain the stellar flux (plus the faint, background contribution due to the occulting body) ratioed to nearby calibration stars using
<table><caption>**Table 5**
Observation Details of the Multi-chord Stellar Occultation on 2014 November 7 with Data Acquisition</caption><thead><tr><th>Site Name<br>(Location)<br>(IAU Code)</th><th>Longitude (E)<br>Latitude (N)<br>Altitude (m)</th><th>Telescope Aperture<br>Detector/Instrument</th><th>Exposure Time<br>Cycle Time<br>(s)</th><th>Observers</th><th>Detection</th></tr></thead><tbody><tr><td>Bosque Alegre<br>(Cordoba, Argentina)<br>[821]</td><td>295° 27′ 01″5<br>-31° 35′ 54″0<br>1213.1</td><td>154 cm<br>CMOS ZWO ASI120MM</td><td>5<br>5.000032</td><td>C. A. Colazo<br>R. Artola<br>R. Melia</td><td>Positive</td></tr><tr><td>Rosário<br>(Rosário, Argentina)<br>...</td><td>299° 21′ 09″0<br>-32° 58′ 16″<br>31.0</td><td>49 cm<br>CCD KAF 8300</td><td>3<br>4</td><td>V. Buso</td><td>Positive</td></tr><tr><td>Montevideo<br>(Montevideo, Uruguay)<br>...</td><td>303° 48′ 37.1<br>-34° 45′ 20″0<br>130.0</td><td>46 cm<br>FLI<br>PL9000</td><td>3<br>4</td><td>G. Tancredi<br>S. Roland<br>L. Almenares<br>A. Ceretta</td><td>Positive</td></tr><tr><td>Santa Rosa<br>(Santa Rosa, Argentina)<br>[J48]</td><td>295° 40′ 32″5<br>-36° 38′ 16″0<br>182.0</td><td>20 cm<br>CCD Meade DSI-I</td><td>25<br>26</td><td>J. Spagnotto</td><td>Positive</td></tr><tr><td>ESO—VLT<br>(Cerro Paranal, Chile)<br>[309]</td><td>289° 35′ 49′9<br>-24° 37′ 31″5<br>2635.0</td><td>400 cm<br>HAWK-I</td><td>0.00664<br>0.00667</td><td>V. D. Ivanov</td><td>Negative</td></tr><tr><td>ESO - La Silla<br>(La Silla, Chile)<br>[809]</td><td>289° 15′ 58′5<br>-29° 15′ 32″3<br>2345.4.0</td><td>355 cm<br>NTT/SOFI</td><td>0.1<br>0.10003</td><td>B. Sicardy<br>D. Bérard</td><td>Negative</td></tr><tr><td>TRAPPIST-S<br>(La Silla, Chile)<br>[J40]</td><td>289° 15′ 38′2<br>-29° 15′ 16″6<br>2317.7</td><td>60 cm<br>TRAPPIST-S</td><td>3<br>4.5</td><td>C. Opitom<br>E. Jehin</td><td>Negative</td></tr><tr><td>PROMPT<br>(Cerro Tololo, Chile)<br>[807]</td><td>289° 11′ 05″5<br>-30° 09′ 56″3<br>2225.0</td><td>40 cm<br>PROMPT</td><td>4<br>6</td><td>J. Pollock</td><td>Negative</td></tr><tr><td>Las Cumbres*<br>(Las Cumbres, Chile)<br>[W85]</td><td>289° 11′ 42″7<br>-30° 10′ 02″6<br>2201.7</td><td>100 cm<br>FLI<br>MicroLine 4720</td><td>3/4<br>5/6</td><td>F. Bianco</td><td>Negative</td></tr><tr><td>San Pedro de Atacama<br>(S. P. de Atacama, Chile)<br>[J16]</td><td>291° 49′ 13″2<br>-22° 57′ 12″1<br>2398.50</td><td>Dall Kirkham 50 cm<br>Apogee U42<br>Ash2 Newtonian 40 cm<br>SBIG-STL11000</td><td>3<br>4<br>15<br>17.59</td><td>A. Maury<br>N. Morales</td><td>Negative</td></tr><tr><td>OPD/LNA<br>(Itajubá, Brazil)<br>[874]</td><td>314° 25′ 03″0<br>-22° 32′ 04″<br>1864.0</td><td>Perkin Elmer 160 cm<br>Andor IXON<br>DU-888E-C00-#BV</td><td>2<br>2.950</td><td>G. Benedetti-Rossi<br>B. Morgado<br>A. R. Gomes-Junior</td><td>Negative<br>(w/ clouds)</td></tr></tbody></table>

Note. FLI stands for Finger Lakes instruments, PI for Princeton instruments, and PL for pro line. * Three 1 m telescopes were used at Las Cumbres (IAU codes W85, W86, and W87). One of them used the *I* filter to minimize the moon contamination while the other two observed in Clear. Exposure time were 3 s for two telescopes (filters *I* and Clear) and 4 s for the third telescope. Readout time is not constant for this observation, but on average is 2 s. Computers are synchronized via NTP to the site GPS and the start of exposure in the three telescopes had a time shift in order to combine the three data sets and obtain one unique light curve without readout time.

procedures based on the standard Daophot IDL routines (Stetson 1987). The multi-chord occultation was harder to analyze. Since the S/N of the comparison and target stars were too low in all data sets, using the IDL routines the derived chords did not produce satisfactory results, as they presented sizes and positions that did not fit a reasonable and realistic ellipse in the sky plane. We then subjected the data sets to a careful treatment by the new photometry task version of the Program for Reduction of Astronomical Images Automatically (PRAIA; Assafin et al. 2011).

PRAIA uses differential aperture photometry and it differs from the IDL routines in many instances, presenting many

useful features for which description is beyond the scope of this work. Among them, the apertures, sky ring sizes, and widths are automatically optimized for better S/N (outside of the event) for each image, following a thorough aperture centering procedure, and the flux is precisely corrected for the variable apertures for all objects on each image. Calibration fluxes are smoothed and the best set of calibration stars is set automatically as a function of the best light curve standard deviation. The flux ratio versus time is further corrected by atmospheric systematic factors (using the observed flux before and after the event) by applying a polynomial fit (we used the third degree in all cases). Finally, the flux ratio is normalized
<table><caption>**Table 6**
Observation Details of the Multi-chord Stellar Occultation on 2014 November 7 with No Data Acquired</caption><thead><tr><th>Site Name<br>(Country)</th><th>Longitude (E)<br>Latitude (N)<br>Altitude (m)</th><th>Telescope Aperture<br>Detector/Instrument</th><th>Observers</th><th>Situation</th></tr></thead><tbody><tr><td>Brasília<br>(Brazil)</td><td>312° 07′ 04″02<br>-15° 47″ 40″37<br>1096</td><td>G11 35 cm<br>SBIG<br>ST8XMEI</td><td>P. Cacella</td><td>Overcast</td></tr><tr><td>Vitória (UFES)<br>(Brazil)</td><td>319° 41′ 33″18<br>-20° 16′ 41″43<br>6</td><td>GSO 35 cm<br>SBIG<br>ST8XME</td><td>M. Malacarne</td><td>Overcast</td></tr><tr><td>Oliveira (SONEAR)<br>(Brazil)</td><td>315° 10′ 27″88<br>-20° 42 54.27<br>984</td><td>45 cm<br>FLI<br>16803</td><td>C. Jacques</td><td>Overcast</td></tr><tr><td>U. Autónoma Tomás Frías<br>(Bolivia)</td><td>295° 22′ 33″00<br>-21° 35′ 45″90<br>1866</td><td>Cassegrain 60 cm<br>FLI<br>IMG1001E</td><td>R. Condori</td><td>Overcast</td></tr><tr><td>São Carlos (UFSCar)<br>(Brazil)</td><td>312° 06′ 50″11<br>-21° 58′ 49″01<br>847</td><td>Meade LX200 30 cm<br>SBIG<br>ST9</td><td>G. Rojas</td><td>Overcast</td></tr><tr><td>Guaratinguetá (UNESP)<br>(Brazil)</td><td>314° 48′ 29″50<br>-22° 48′ 10″02<br>573</td><td>Meade LX200 40 cm<br>SBIG<br>ST-7XME</td><td>R. Sfair</td><td>Overcast</td></tr><tr><td>Rio de Janeiro (ON)<br>(Brazil)</td><td>316° 46′ 32,83″<br>-22° 53′ 46,85″<br>29</td><td>Meade LX200 30 cm<br>Raptor Merlin 246</td><td>F. Braga-Ribas<br>J. I. B. Camargo</td><td>Overcast</td></tr><tr><td>São José dos Campos (INPE)<br>(Brazil)</td><td>314° 08′ 16″<br>-23° 12′ 33″<br>975</td><td>C11 28 cm<br>SBIG<br>ST-7XE</td><td>A. C. Milone</td><td>Overcast</td></tr><tr><td>Ponta Grossa (UEPG)<br>(Brazil)</td><td>309° 53′ 41″39<br>-25° 05′ 40″27<br>910</td><td>Meade RC40 40 cm<br>SBIG<br>ST-6303E</td><td>M. Emilio</td><td>Overcast</td></tr><tr><td>Foz do Iguaçu<br>(Brazil)</td><td>305° 24′ 22″55<br>-25° 26′ 05″45<br>185</td><td>C11 28 cm<br>SBIG<br>ST-7XME</td><td>D. I. Machado</td><td>Overcast</td></tr><tr><td>CASLEO<br>(Argentina)</td><td>290° 41′ 15″9<br>-31° 47′ 54″7<br>2552</td><td>215 cm<br>PI-2048B</td><td>R. Gil-Hutton</td><td>Overcast</td></tr><tr><td>Montevideo<br>(Uruguay)</td><td>303° 48′ 37.1<br>-34° 45′ 20″0<br>130.0</td><td>Meade LX200 30 cm<br>Raptor Merlin 246</td><td>G. Tancredi<br>S. Roland<br>L. Almenares<br>A. Ceretta</td><td>Star not detected</td></tr></tbody></table>

Note. FLI stands for Finger Lakes instruments and PI for Princeton instruments.

and the resulting normalized light curve stored with many kinds of information (errors, photometry parameters, etc.). Note that the three stars magnitudes are about 16, while it is almost 20 for the TNO. This means that VS2 is not visible in any of the equipment and exposure setup used in the observations and the stars disappear completely during the three occultations, with no residual light from VS2, i.e., we assume that the flux drops to zero in each of these observations. The resulting light curves with positive occultation detections are presented in Figure 5. The light curves flux standard deviation are 21.8% for the 2013 occultation, 12.9% and 6.6% for the two telescopes of 0.7 m and 1 m, respectively, from the 2014 March event, while it is 6.9% for Bosque Alegre, 14.7% for Rosário, 14.5% for

Santa Rosa, and 13.1% for Montevideo in the multi-chord occultation.

The start (star disappearance) and end (star reappearance) times of the occultation are determined by fitting each event by a sharp opaque edge model, after convolution by (i) Fresnel diffraction, (ii) finite CCD bandwidth, (iii) finite stellar diameter, and (iv) finite integration time (see details in Widemann et al. 2009; Braga-Ribas et al. 2013; Ortiz et al. 2017). Note that the Fresnel monochromatic diffraction operates over the Fresnel scale $F = \sqrt{\lambda D/2}$, where $\lambda$ is the wavelength of observation and $D$ is the geocentric distance of the object.

The stellar diameter projected at VS2's distance is estimated using the *B*, *V*, and *K* apparent magnitudes provided by the
<figure><img src="image_4.png" /><figcaption>**Figure 1.** Post-occultation map for the 2013 December 12 occultation. North is up and east is right. Blue lines represent the object equivalent diameter from the 2014 November 7 occultation of 564.8 km. Red dots represent the position of the center of the body spaced every 30 s, the bigger dot corresponding to 20:10:44 UTC. The direction of the shadow is shown by the arrow at the right corner. The green dot is the site position with positive detection (see Table 4).</figcaption></figure>

<figure><img src="image_5.png" /><figcaption>**Figure 2.** Post-occultation map for the 2014 March 4 occultation. North is up and east is right. Blue lines represent the object equivalent diameter from the 2014 November 7 occultation of 564.8 km. Red dots represent the position of the center of the body spaced every 30 s, the bigger dot corresponding to 19:55:54 UTC. The direction of the shadow is shown by the arrow at the right corner. The green dot is the site position with positive detection, while the blue dots are the negative detections and the white dots are the sites with weather overcast (see Table 4).</figcaption></figure>

NOMAD catalog (Zacharias et al. 2004) and the formulae of van Belle (1999)—see Table 2. Using the apparent TNO motion relative to the star, it is possible to express the integration times in actual distances traveled between adjacent data points in the sky plane. Each effect is evaluated in Table 7. We note that our data are dominated by the integration times, and not by Fresnel diffraction or stellar diameter. The occultation timings are then obtained by minimizing a classical $\chi^2$ function, using the same procedures as described in Sicardy et al. (2011).

<figure><img src="image_6.png" /><figcaption>**Figure 3.** Post-occultation map for the 2014 November 7 occultation. North is up and east is right. Blue lines represent the object equivalent diameter obtained from this occultation, of 564.8 km. Red dots represent the position of the center of the body spaced every 30 s with the bigger one corresponding to 04:22:00 UTC. The direction of the shadow is shown by the arrow at the right corner. The green dots are the site positions with positive detections, while the blue dots are the negative detections and the white dots are the sites with weather overcast (see Tables 5 and 6).</figcaption></figure>

<figure><img src="image_7.png" /><figcaption>**Figure 4.** Rotational light curve for VS2 obtained on 2014 November 15 (in red), 16 (in blue), and 17 (in green) at Sierra Nevada Observatory (see the text for details) just a few days after the occultation. Light blue dots on the bottom panel show the residuals between the observed points and the fit (blue line). The estimated rotational phase for the 2014 November stellar occultation is near one of the brightness maxima. We used the occultation date as the zero phase date.</figcaption></figure>

### 3.1. Single-chord Events

The single-chord occultation of 2013 December 12 was detected from La Réunion Island and had a shadow velocity of 24.91 km s⁻¹, while the 2014 March 4 event had a shadow velocity of 8.99 km s⁻¹ and it was detected by two telescopes at the same site in Israel (see Table 4 for details). The best fits to the ingress and egress profiles are presented in Figure 6. The occultation durations of 15.80±0.40 s (La Réunion) and 45.75 ± 1.85 s (Israel) then correspond to occultation chords with lengths of 393.6±10.0 km and 411.3 ± 16.6 km, respectively, see Table 8.
<figure><img src="image_8.png" /><figcaption>**Figure 5.** Normalized light curves of the positive detections of the stellar occultations—see Tables 4 and 5 for site details. Top: single-chord event in 2013 December 12; center: single-chord event in 2014 March 4, showing the detection with the two telescopes at the same site; bottom: multi-chord event on 2014 November. Chords are shifted in flux for better visualization. Vertical lines are the uncertainties in photometry for each point.</figcaption></figure>

The timings in Table 8 provide the extremities of the corresponding occultation chords projected in the sky plane. More precisely, they give the position (*f*, *g*) of the star projected in the sky plane and relative to the center of the object. The quantities *f* and *g* are expressed in kilometers, and counted positively toward local celestial east and celestial north, respectively. Note that the calculated position (*f*, *g*) depends on both the ephemeris used for the body and the adopted star position, which are both determined to within a finite accuracy. As such, (*f*, *g*) must be corrected by an offset that can be determined if several chords are available, thus pinning down the position of the object's center position relative to the star.

This is not possible for a single-chord event. In that case, a circle was adjusted to each chord using the area-equivalent

<figure><img src="image_9.png" /><figcaption>**Figure 6.** Best fit for the star disappearance (left) and reappearance (right) instants of the single-chord events—see Table 4 for sites details. Top: 2013 December 12 occultation; center: 2014 March 4 occultation observed with the 0.7 m telescope; bottom: 2014 March 4 occultation observed with the 1 m telescope. Data points are shown in black, the geometric square well models are in blue, the effect of diffraction is shown in green, the effect of finite integration time is added in the orange, and the red lines show all the effects accounted for and then used for comparison with the actual data points.</figcaption></figure>

<table><caption>**Table 7**
Highest Amplitude of Effects to Consider in Order to Determine the Start and End Instants of the Three Stellar Occultations</caption><thead><tr><th>Occultation Date</th><th>D<sup>a</sup><br>(au)</th><th>Fresnel<sup>b</sup><br>(km)</th><th>Star Diameter<sup>c</sup><br>(km)</th><th>Integration Time<sup>d</sup><br>(km)</th></tr></thead><tbody><tr><td>2013 Dec 12</td><td>35.58</td><td>1.73</td><td>0.41</td><td>18.7</td></tr><tr><td>2014 Mar 4</td><td>36.56</td><td>1.78</td><td>0.24</td><td>36.0</td></tr><tr><td>2014 Nov 7</td><td>35.72</td><td>1.74</td><td>0.56</td><td>63.5</td></tr></tbody></table>

**Notes.**

diameter of 564.8 km in order to determine the center of the TNO in the plane of the sky relative to the star at a prescribed time. Two solutions are then equally possible, depending on whether the center of the body went south or north of the star,

<footnote><sup>a</sup> Geocentric distance in astronomical units from Horizons/JPL (https://ssd.jpl.nasa.gov/horizons.cgi) using JPL#30 and DE431.</footnote>

<footnote><sup>b</sup> Fresnel diffraction effect ($F = \sqrt{\lambda D/2}$); see the text.</footnote>

<footnote><sup>c</sup> Stellar diameter projected at VS2's geocentric distance calculated using (van Belle 1999) and assuming supergiant stars.</footnote>

<footnote><sup>d</sup> The finite integration time is given for the smallest exposure time of the positive detections. It is expressed in kilometers after accounting for the velocity of the body projected in the sky plane; see the text.</footnote>
<table><caption>**Table 8**
Timings of Star Disappearances and Reappearances</caption><thead><tr><th>Occultation Event<sup>a</sup></th><th>Disappearance<sup>b</sup><br>(s)</th><th>Reappearance<sup>b</sup><br>(s)</th><th>Chord Size<br>(km)</th></tr></thead><tbody><tr><td>2013 Dec 12</td><td>72785.15<sup>+0.15</sup><sub>-0.20</sub></td><td>72800.95<sup>+0.15</sup><sub>-0.25</sub></td><td>393.6<sup>+8.7</sup><sub>-10.0</sub></td></tr><tr><td>2014 Mar 4 (1 m)</td><td>71865.15 ± 1.55</td><td>71910.9 ± 0.3</td><td>411.3 ± 16.6</td></tr><tr><td>2014 Mar 4 (70 cm)</td><td>71861.25 ± 1.15</td><td>71909.7<sup>+0.5</sup><sub>-0.6</sub></td><td>435.6<sup>+14.8</sup><sub>-15.7</sub></td></tr><tr><td>2014 Nov 7 (BA)</td><td>15903.354 ± 0.646</td><td>15927.400 ± 0.642</td><td>517.7 ± 27.8</td></tr><tr><td>2014 Nov 7 (RO)</td><td>15883.351 ± 0.887</td><td>15913.006 ± 0.670</td><td>638.5 ± 33.5</td></tr><tr><td>2014 Nov 7 (SR)</td><td>15901.0 ± 8.0</td><td>15928.0 ± 8.0</td><td>581 ± 345</td></tr><tr><td>2014 Nov 7 (MO)</td><td>15868.286 ± 0.800</td><td>15890.583 ± 0.800</td><td>480.1 ± 34.4</td></tr></tbody></table>

**Notes.** BA: Bosque Alegre; RO: Rosário; SR: Santa Rosa; MO: Montevideo.

<sup>a</sup> See Tables 4–6 for occultation details.

<sup>b</sup> Times are given in seconds after 00:00:00.000 UTC of the occultation day.

<table><caption>**Table 9**
Astrometric R.A. and Decl., and the Offsets (*f*, *g*) for VS2 Derived from the Three Stellar Occultations</caption><thead><tr><th>Occultation Date</th><th>R.A. ICRS<br>(JPL30 and DE431 + Offset)</th><th>Decl. ICRS</th><th>f<sup>a</sup><br>(mas)</th><th>g<sup>a</sup><br>(mas)</th></tr></thead><tbody><tr><td>2013 Dec 12 20:13:00.000</td><td>04<sup>b</sup>36<sup>m</sup>19<sup>s</sup>83798</td><td>+33° 55′ 55″2525</td><td>-15.7 ± 2.4</td><td>-17.5 ± 8.1</td></tr><tr><td>2014 Mar 4 19:55:00.000</td><td>04<sup>b</sup>32<sup>m</sup>04<sup>s</sup>47805</td><td>+33° 24′ 42″0523</td><td>-12.2 ± 5.4</td><td>-23.7 ± 6.1</td></tr><tr><td>2014 Nov 7 04:22:00.000</td><td>04<sup>b</sup>48<sup>m</sup>32<sup>s</sup>13788</td><td>+33° 58′ 36″1927</td><td>-60.1 ± 0.3</td><td>-24.5 ± 0.4</td></tr></tbody></table>

**Note.**

<sup>a</sup> The offsets depend also on the adopted star positions given on Table 2.

as seen from Earth. We consider here the average of the two solutions, still resulting in an accurate astrometric positions for VS2 at the time of the two occultations, with uncertainties smaller than 25 mas for the two single-chord occultations. Note that the diameter used here was obtained from the 2014 November 7 occultation, when the object was at one of the maxima in brightness of the light curve, but we do not know the rotational phase for the other two single-chord occultations. This means that probably the shadow path must be smaller than 564.8 km. Anyway, a difference of 50 km in radius will give an uncertainty of less than 1.7 mas in the object position, which is within the error bars. The object positions in R.A. and decl. corrected by the offsets obtained from the occultations are given in Table 9.

## 3.2. Multi-chord Event

The 2014 November 7 occultation had a typical shadow velocity of 21.53 km s⁻¹ and observers attempted to record the event from 21 different sites. Unfortunately, due to bad weather only eight sites were successful in gathering data, among which four detected the event (see Tables 5 and 6). Despite nondetec- tions, the chords obtained at Cerro Tololo (PROMPT), Las Cumbres, and La Silla (NTT/SOFI; Son of ISAAC Moorwood et al. 1998a, 1998b) provide important constraints on the size of the object, being close to the occultation path.

The resulting light curves with positive detection, shown in Figure 5, provide *N* = 8 chord extremities, from which we derived the disappearance and reappearance instants of the occultation (see Table 8). The best fits for the ingress and egress times are presented in Figure 7.

The shape of the body’s limb is assumed to be an ellipse, as adopted in previous works, see Braga-Ribas et al. (2013, 2014), Benedetti-Rossi et al. (2016), Dias-Oliveira et al. (2017), and Ortiz et al. (2017). The ellipse is defined by *M* = 5 adjustable

parameters: the position of the ellipse center, (*f*<sub>c</sub>, *g*<sub>c</sub>); the apparent semimajor axis, *a*′; the apparent oblateness, *ϵ*′ = (*a*′ − *b*′)/*a*′ (where *b*′ is the apparent semimajor axis); and the position angle of the pole *P*′ of *b*′, which is the apparent position angle of the pole measured eastward from celestial north. Note that the center (*f*<sub>c</sub>, *g*<sub>c</sub>) actually measures the offsets in R.A. and decl. to be applied to the adopted ephemeris, assuming that the star position is correct. Note also that the quantities *a*′, *b*′, *f*<sub>c</sub>, and *g*<sub>c</sub> are all expressed in kilometers.

The statistical significance of the fit is evaluated from the χ² per degree of freedom (pdf) defined as χ²<sub>pdf</sub> = χ²/(N − *M*), which should be close to unity for a correct fit. The individual 1σ error bar of each parameter is obtained by varying that parameter from its nominal solution value (keeping the other parameters free), so that χ² varies from its minimum value χ²<sub>min</sub> to χ²<sub>min</sub> + 1.

Using the timings for the star disappearance and reappearance obtained from the four positive detections (Table 8), we obtain a best-fitting ellipse with χ²<sub>pdf</sub> = 0.78 and its parameters *a*′ = 313.8 ± 7.1 km, *ϵ*′ = 0.190<sup>+0.052</sup><sub>-0.060</sub> (*b*′ = 254.8<sup>+25.0</sup><sub>-21.7</sub>), (*f*, *g*) = (1558.1 ± 8.1, 634.6 ± 11.0) km, and *P*′ = 5° ± 7° with an equivalent radius of *R*<sub>eq</sub> = 282.4<sup>+16.9</sup><sub>-15.1</sub> km, as shown in Figure 8.

# 4. Results

## 4.1. Size and Shape

Homogeneous rocky objects with diameters on the order of 1000 km or more are expected to have reached the hydrostatic equilibrium, assuming Maclaurin spheroidal or Jacobi ellipsoidal shapes (Chandrasekhar 1987). The critical diameter allowing to reach this equilibrium is discussed by Tancredi & Favre (2008). In particular, hydrostatic equilibrium could be reached for smaller objects, if made of ices or of a combination of ices and rocks. Nevertheless, the formation scenario and collisional history of
<figure><img src="image_10.png" /><figcaption>**Figure 7.** Best fit for the star disappearance (left) and reappearance (right) instants of the multi-chord event of 2014 November 7—see Table 5 for sites details. From top to bottom: Bosque Alegre, Rosario, Santa Rosa, and Montevideo. Line colors are the same as in Figure 6.</figcaption></figure>

individual objects are not known, allowing objects with complex internal structures (i.e., differentiated density layers, granular for apeeble accretion scenario, etc.). Thus, we cannot discard solutions that diverge somewhat from the expected equilibrium figure given by Chandrasekhar (1987). The case of the dwarf planet Haumea (Ortiz et al. 2017) is a good example of a body with sizes that does not match a hydrostatic equilibrium figure for a geologically homogeneous fluid body.

The rotational light curve of Figure 4 shows that VS2 is consistent with an object having a triaxial ellipsoidal shape. The short-term photometry of the body presented in Section 2.3 showed that the multi-chord occultation took place near one of the absolute brightness maxima of the rotational light curve, meaning that VS2 occulted the star when its apparent surface area was near its maximum. Assuming a triaxial shape, this occurs when its longest axis *a* was perpendicular to the line of sight, so that the semimajor axis of the projected ellipse *a*′ is

<figure><img src="image_11.png" /><figcaption>**Figure 8.** Occultation chords (blue) and their uncertainties (red) from the multi-chord event of 2014 November with the best elliptical fit (black). The offsets in both axis are with respect to ephemeride positions obtained using JPL#30 and DE431. The arrow shows the direction of the shadow motion and the green lines represent the two closest sites from where there were no detection of the occultation (see Table 5). The equivalent circle (that has the same area as the ellipse) is shown in orange with a diameter of 564.8 km.</figcaption></figure>

equal to the real 3D axes. With the known rotation period and phase, we verified if it was possible to use one or the two single-chord occultation results to improve the accuracy on the determination of size and shape for the TNO. The problem is that even if the rotational phase and period are well determined, the ephemeris for the object is not well constrained, with uncertainties of about 50 mas, in a way that the chord from the two occultations can fit anywhere in the multi-chord apparent ellipse. With this uncertainty we could not use any of the other chords to obtain further constraints on VS2's shape using the same method that Dias-Oliveira et al. (2017) did for the TNO 2003 AZ84.

Considering the limb fit to the apparent ellipse obtained from the multi-chord occultation and the light curve amplitude, we can use the same procedure as in Sicardy et al. (2011; also used in Braga-Ribas et al. 2013 and Ortiz et al. 2017) to derive the three axes of the body. Our best solution for a triaxial shape has semi-axis values of *a* = 313.8 ± 7.1 km, *b* = 265.5<sup>+8.8</sup><sub>-9.8</sub> km (β = 0.846<sup>+0.031</sup><sub>-0.028</sub>), and *c* = 247.3<sup>+26.6</sup><sub>-43.6</sub> km (γ = 0.788<sup>+0.085</sup><sub>-0.139</sub>), with the *c*-axis inclined θ = 65<sup>+15°</sup><sub>-10°</sub> with respect to the observer, which is not consistent with the Jacobi equilibrium figure—see the Appendix for details. This solution gives a spherical volume equivalent diameter of 548.3<sup>+29.5</sup><sub>-44.6</sub>, from which we can derive (using *H*<sub>v</sub> = 4.130 ± 0.070 mag) the geometric albedo *p*<sub>V</sub> = 0.131<sup>+0.024</sup><sub>-0.013</sub>.

Since we cannot make any constraint about VS2’s density using the triaxial approach without making any other assumption of its mass, we can make an assumption that VS2 has an oblate shape, as a Maclaurin equilibrium figure with a semi-axis of *a* = 313.8 km and 254.8 ≤ *c* ≤ 313.8 km. Note that with this assumption, we are also assuming that the light curve amplitude variation is due to albedo features. Making this assumption and combining with the rotation period (given in
<table><caption>**Table 10**
Solution for 2003 VS₂ Derived from the Multi-chord Stellar Occultation and the Rotational Light Curve</caption><thead><tr><th>Parameter</th><th>Value</th></tr></thead><tbody><tr><td>a' × b' (km)</td><td>313.8 ± 7.1 × 254.8<sup>+25.0</sup><sub>-21.7</sub></td></tr><tr><td>(f, g) (km)<sup>1</sup></td><td>(-1558.1 ± 8.1; -634.6 ± 11.0)</td></tr><tr><td>Position angle of projected ellipse (deg)</td><td>5 ± 7</td></tr><tr><td>Light curve amplitude<sup>2</sup>—Δm—(mag)</td><td>0.141 ± 0.009</td></tr><tr><td>Geometric albedo—p<sub>V</sub></td><td>0.131<sup>+0.024</sup><sub>-0.013</sub></td></tr><tr><td>a × b × c (km)</td><td>313.8 ± 7.1 × 265.5<sup>+8.8</sup><sub>-9.8</sub> × 247.3<sup>+26.6</sup><sub>-43.6</sub></td></tr><tr><td>Aspect angle (deg)</td><td>65<sup>+15</sup><sub>-10</sub></td></tr><tr><td>Equivalent diameter (km)</td><td>548.3<sup>+29.5</sup><sub>-44.6</sub></td></tr><tr><td>Density—Maclaurin (kg m<sup>-3</sup>)</td><td>1400<sup>+1000</sup><sub>-300</sub></td></tr></tbody></table>

Note. ¹ Offsets obtained with respect to JPL30 + DE431. It also depends on the star position given in Table 2. ² Period used to fit was 7.4175285 ± 0.00001 hr from Santos-Sanz et al. (2017).

Table 1), we can derive its density of $1400_{-300}^{+1000}$ kg m$^{-3}$. All those values are summarized in Table 10.

## 4.2. Atmosphere and Secondary Detections

On the largest TNOs the most common volatiles found through spectroscopic studies are water ice ($H_2O$), methane ($CH_4$), ammonia ($NH_3$), molecular nitrogen ($N_2$), and even methanol ($CH_3OH$) and ethane ($C_2H_6$) for a few objects (Guilbert et al. 2009; Barucci et al. 2011). In the case of VS2, near-infrared spectra of this body shows the presence of exposed water ice (Barkume et al. 2008), but no other volatile is reported to be detected. Data from the three stellar occultations showed no compelling evidence for a global atmosphere.

Using the same method for Quaoar (Braga-Ribas et al. 2013), we have modeled synthetic occultation light curves using a ray tracing code as described in Sicardy et al. (2011) and Widemann et al. (2009). Several synthetic light curves are then compared to our best light curve obtained at the Wise Observatory on the single-chord occultation on 2014 March 4 providing a range of $\chi^2$ with the detection threshold. We test two possible atmospheric models: one with an isothermal $N_2$ atmosphere, as this gas is the most volatile one among those listed above. We consider a typical temperature of 40 K for this model, but the result is weakly dependent on this particular value. The other model is a $N_2$ atmosphere with $CH_4$ as a minor species, starting near 40 K at the surface and ramping up to 100 K near 25 km altitude due to methane infrared absorption, as is observed in Pluto's atmosphere (Hinson et al. 2017).

Our $\chi^2$-tests show that for both models, the Wise data are consistent with no atmosphere, with upper limits of 0.2 microbar and 1 microbar at the $1\sigma$ and $3\sigma$ levels, respectively (see Figure 9).

We also searched for secondary events that may be related to a satellite, ring, or some material orbiting VS2. Using our highest time-resolution light curve—obtained with the 3.5 m European Southern Observatory (ESO) New Technology Telescope (NTT) telescope with 0.1 s of cycle time and roughly positioned 400 km from VS2's center—we can estimate an upper limit of 0.2 km for an opaque object (or 18.5 km for any material with optical depth of 0.1). Considering the best chord with positive detection—from Bosque Alegre with 5 s cycle time—those limits are 107.7 km for opaque and 732 km for material with optical depth of 0.1.

One secondary feature was also detected in the NTT light curve, as shown in Figure 10. The drop has a relative depth of less than 10% (compared to the light curve standard deviation of

<figure><img src="image_12.png" /><figcaption>**Figure 9.** Light curve from the single-chord stellar occultation of 2014 March 4 obtained at the Wise Observatory (black) and the two models for the isothermic pure nitrogen (red) and Pluto-like (blue) atmospheres within a $3\sigma$ limit.</figcaption></figure>

<figure><img src="image_13.png" /><figcaption>**Figure 10.** Light curve from the multi-chord stellar occultation of 2014 November 7 obtained at La Silla/NTT. All the flux drops in the curve are due to the star being too close to the edge of the CCD or being influenced by the bad row of pixels in the center of the CCD. Note the relative flux drop of 0.09 that happens between instants 15885 and 15915 s, fitted using a square well model, represented by the red line.</figcaption></figure>

8.6%) and a duration of about 29 s, corresponding to a chord length of 625 km. However, this drop remains marginally significant and cannot be confirmed or rejected using the nearby chords from Cerro Tololo (insufficient S/N) and Pico dos Dias (overcast). This drop could be due to an occultation of a companion star by the object, causing a secondary occultation as happened during a stellar occultation by the Centaur object Chariklo in 2014—see Leiva et al. (2017) and Bérard et al. (2017). Note that with the size of VS2 and the geometry of the occultation, the Bosque Alegre positive light curve (data set with the best S/N) should not present any sign of this possible
secondary occultation. Also, considering the main star magnitude ($V = 15.82$) and the drop of 10%, the companion would have a magnitude of $V = 18.21$. If this star is closer than $0.8$ it is likely not to be present in the GDR2 catalog (Brandeker & Cataldi 2019). We estimate that the companion star should be closer than $0.5$, so the existence of a stellar companion remains consistent with its absence in the GDR2.

Another possibility is that this drop is due to some diffuse material around VS2, but since none of the other light curves show such a secondary event to within the $5\sigma$ level, this hypothesis is less likely. Also, a careful look on the NTT images presents some issues. The images have 320 versus 32 pixels with a dead row in the middle of the CCD (row 69). There was only one reference star in the field for photometric calibration, which is fainter than the occulted star, and the telescope tracking was not perfect, in a way that sometimes one or both stars were close to the edge of the CCD. So we cannot completely discard instrumental effects or reduction artifacts on this data set. In conclusion, more occultations are needed in order to confirm or discard the presence of any secondary feature.

## 5. Conclusions

We observed two single-chord and one multi-chord stellar occultations by the plutino (84922) 2003 VS2. Observations were performed on 2013 December 12 in La Réunion, on 2014 March 4 in Israel, and from 11 sites in South America on 2014 November 7, with four positive detections. The two single-chord events are consistent with the multi-chord solutions but do not bring further constraints to the object shape. However, they provide good astrometric constraints, thus improving the ephemeris of the body.

Using the timings of the 2014 November 7 event, we find the apparent ellipse with area-equivalent diameter of $564.8_{-30.2}^{+33.8}$ km. Assuming that VS2 has a triaxial shape and combining the occultation data with the light curve amplitude ($\Delta m = 0.141 \pm 0.009$) obtained from photometric observations of the TNO, we derive the object 3D shape, with semi-axes of $a = 313.8 \pm 7.1$ km, $b = 265.5_{-9.8}^{+8.8}$ km, and $c = 247.3_{-43.6}^{+26.6}$ km, with the c-axis inclined $\theta = 65_{-10}^{+15}^\circ$ with respect to the observer. Those values are not consistent with a Jacobi equilibrium figure, but are in accordance, at the $1\sigma$ level, with previous values published by Mommert et al. (2012), Stansberry et al. (2008), and Lellouch et al. (2013) (area-equivalent diameter of $523_{-34.4}^{+35.1}$ km).

We also considered the possibility that VS2 is an oblate Maclaurin equilibrium figure (with $a = b > c$). With this assumption, the light curve amplitude variation is due to albedo features and VS2 would have a density of $\rho = 1400_{-300}^{+1000}$ kg m$^{-3}$. The light curve amplitude can be explained in the extreme cases to be due only to the shape (with the triaxial shape model) or to a darker or brighter spot on the surface of about 100 km—16% of object diameter—if VS2 is a Maclaurin.

A possible secondary event was detected in the NTT data, that may be caused by the presence of some feature in VS2's surroundings or due to a secondary occultation occurred by the presence of a companion star, but the data are insufficient to make any deep conclusions. We also derived an upper limit for a global atmosphere (either isothermic pure nitrogen or Pluto-like) of 1 microbar ($3\sigma$).

Although the rotation period is well defined, and values for the geometric albedo and size of VS2 are well determined with the occultations, its 3D shape (and as a consequence, its

density) still need more observations in order to constrain the values obtained in this work.

G.B.-R. is thankful for the support of the CAPES and FAPERJ/PAPDRJ (E26/203.173/2016) grant. Part of the research leading to these results has received funding from the European Research Council under the European Community's H2020 (2014-2020/ERC grant agreement No. 669416 "LUCKY STAR"). The research leading to these results has received funding from the European Union's Horizon 2020 Research and Innovation Programme, under grant agreement No. 687378 (SBNAF). P.S.-S. and J.L.O. acknowledge the financial support by the Spanish grant AYA-2017-84637-R and the Proyecto de Excelencia de la Junta de Andalucía J.A. 2012-FQM1776. P.S.-S., J.L.O., and R.D. acknowledge financial support from the State Agency for Research of the Spanish MCIU through the "Center of Excellence Severo Ochoa" award for the Instituto de Astrofísica de Andalucía (SEV-2017-0709). Based on observations made with ESO Telescopes at the La Silla Paranal Observatory under program ID 094.C-0352. M.A. thanks CNPq (grants 427700/2018-3, 310683/2017-3, and 473002/2013-2) and FAPERJ (grant E-26/111.488/2013). J.I.B.C. acknowledges CNPq grant 308150/2016-3. R.V.-M. thanks grants: CNPq-304544/2017-5, 401903/2016-8, Faperj: PAPDRJ-45/2013, and E-26/203.026/2015. F.B.-R. acknowledges CNPq grant 309578/2017-5. E.F.-V. acknowledges UFC 2017 Preeminent Postdoctoral Program (P³). TRAPPIST-South is a project funded by the Belgian Fonds (National) de la Recherche Scientifique (F.R.S.-FNRS) under grant FRFC 2.5.594.09. F. E.J. is a FNRS Senior Research Associate. A.A.C. acknowledges support from FAPERJ (grant E-26/203.186/2016) and CNPq grants (304971/2016-2 and 401669/2016-5). B.M. thanks the CAPES/Cofecub-394/2016-05 grant. A.R.G.J. and R.S. thank the financial support of FAPESP (proc. 2018/11239-8, proc. 2011/08171-3, proc. 2016/24561-0). A.M. thanks Caisey Harlingten for the use of his 50 cm telescope. We thank V. Buso and R. Condori for the observation efforts.

## Appendix
Deriving VS2's Shape and Density

We describe here the procedure used to determine the values of the semi-axis of a triaxial ellipsoid ($a > b > c$) and the angle between the c-axis and the observer ($\theta$).

First we have the relations from Sicardy et al. (2011; also used in Braga-Ribas et al. 2013 and Ortiz et al. 2017) between the real c and b semiminor axis and the observed $b'$:

$$b' = \sqrt{c^2 \sin^2 \theta + b^2 \cos^2 \theta}, \qquad (1)$$

where in our case we have $b' = 254.8_{-21.7}^{+25.0}$ km as the semiminor axis of the ellipse observed in occultation. Multiplying both sides by $1/a$ and replacing $b/a$ by $\beta$ and $c/a$ by $\gamma$, we can rewrite Equation (1) and rearrange the terms to have a direct relation of $\gamma$ as a function of $\beta$:

$$\gamma^2 = \frac{\beta'^2 - \beta^2 \cos^2 \theta}{\sin^2 \theta}. \qquad (2)$$

Note that since the multi-chord occultation took place near the maximum brightness, we have $a' = a$, and so $\beta' = b'/a$. Second, using Equation (8) from Sicardy et al. (2011; SI)
<figure><img src="image_14.png" /><figcaption>**Figure 11.** Relation between $\beta = b/a$ and $\gamma = c/a$ obtained using Equation (2). Black lines correspond to a value of $\theta$ between 5° and 90° (from left to right), every 5°. The orange line is the value for $\beta = \gamma$. The red dotted vertical line is the value of $\beta' = 0.811$, observed in occultation.</figcaption></figure>

we have

$$ \Delta m = -1.25 \cdot \log_{10} \left[ \frac{1 + \gamma^2 \tan^2 \theta}{1 + (\gamma/\beta)^2 \tan^2 \theta} \right], \quad (3) $$

where $\Delta m$ is the amplitude of the rotational light curve ($0.141 \pm 0.009$ mag; see Table 1). We can rearrange the terms in the equation and obtain another relation of $\gamma$ as function of $\beta$:

$$ \gamma^2 = \frac{\beta^2 (\xi - 1)}{\operatorname{tg}^2 \theta (\beta^2 - \xi)}, \quad (4) $$

where $\xi = 10^{-\Delta m/1.25}$. Those two equations combined (Equations (2) and (4)) give us constraints on the values of $\beta$ for every $\theta$, but without any further assumptions, $\gamma$ can assume any value from zero (or more specifically, undefined) to $\beta$ and $\theta$ can have values from 0° to 90°.

Now, assuming that VS2 is large enough to achieve hydrostatic equilibrium, we have the limits for $\beta$ and $\gamma$ as a Jabobi-shape object. Tanredi & Favre (2008) present the relations $\Gamma$—associated with the angular momentum $L$—and $\Omega$—associated with the angular velocity $\omega$—given respectively by:

$$ \Gamma = \frac{L}{\sqrt{GM^3 R}}, \quad (5) $$

and

$$ \Omega = \frac{\omega^2}{\pi G \rho}, \quad (6) $$

where $G$ is the gravitational constant, $M$, $R$, and $\rho$ are the mass, equivalent radius, and density of the body, respectively. For a Jacobi object, there is a lower and an upper limit for the two quantities: $0.303 \le \Gamma \le 0.390$, and $0.284 \le \Omega \le 0.374$; which also limit the values of $\beta$ between 0.432 and 1, and $\gamma$ between 0.345 and 0.583. The shapes of Jacobi ellipsoids in terms of the semi-axes ($a$, $b$, and $c$) can only be obtained by solving Chandrasekhar (1969)—also in Lacerda & Jewitt (2007) and Sicardy et al. (2011) supplementary

information—integrals:

$$ \begin{aligned} & \beta^2 \int_0^\infty \frac{du}{(1+u)(\beta^2+u)\Delta(\beta, \gamma, u)} \\ &= \gamma^2 \int_0^\infty \frac{du}{(\gamma^2+u)\Delta(\beta, \gamma, u)} \end{aligned} \quad (7) $$

and

$$ \frac{\omega^2 a^3}{GM} = \frac{3}{2} \int_0^\infty \frac{u \, du}{(1+u)(\beta^2+u)\Delta(\beta, \gamma, u)}, \quad (8) $$

where $\Delta(\beta, \gamma, u) = [(1+u)(\beta^2+u)(\gamma^2+u)]^{1/2}$. Once $\beta = b/a$ is given, Equation (7) yields $\gamma$, which in turn allows for the calculation of $\omega$ (or $\rho$), using Equation (8).

We can plot the two relations between $\gamma$ and $\beta$ for a triaxial body (Equation (2) is represented in Figure 11 and Equation (4) in Figure 12) and verify if there exists a set of possible common solutions for every value of $\theta$. Each of the intersections between the two curves for a same $\theta$ are presented as gray lines in Figure 13. If we also plot the Jacobi relation between $\gamma$ and $\beta$ (blue line in Figure 13) we should see a solution that intersects the three curves, which is not the case for VS2. This means that there is a solution for a triaxial shape but this solution is not a Jacobi shape.

In the case of VS2 the best solution for a triaxial shape have values for $\theta = 65^{+15}_{-10}$ degrees, with $\beta = 0.846^{+0.031}_{-0.028}$ ($b = 265.5^{+8.8}_{-9.8}$ km), $\gamma = 0.788^{+0.085}_{-0.139}$ ($c = 247.3^{+26.6}_{-43.6}$ km). Those values are also given on Table 10.

In order to have a Jacobi solution, we can explore different values for $\Delta m$, from zero to the nominal value of 0.141 mag. In fact we are assuming that some of the light curve contribution is due to VS2's shape and some due to albedo variation in the surface. When we try values for $\Delta m$ smaller than 0.141 the lines in Figure 12 will move to the right and so more intersections with the lines from Figure 11 will be available, i.e., there will be an intersection between the gray and the blue lines in Figure 13. For $\Delta m = 0.015$ mag we find a Jacobi solution with $\beta = 0.908$ ($b = 284.9$ km), $\gamma = 0.553$ ($c = 173.5$ km), and $\theta = 75^\circ$.

We can also explore the other extreme and assume that VS2 has an oblate Maclaurin shape (with semi-axis
<figure><img src="image_15.png" /><figcaption>**Figure 12.** Relation between $\beta = b/a$ and $\gamma = c/a$ obtained using Equation (4). Black lines correspond to a value of $\theta$ between 0° and 90° (from left to right), every 5°. The orange line is the value for $\beta = \gamma$. The red dotted vertical line is the value of $\beta' = 0.811$, observed in occultation.</figcaption></figure>

<figure><img src="image_16.png" /><figcaption>**Figure 13.** Relation between $\beta = b/a$ and $\gamma = c/a$ combining Equation (2) (black dotted lines—as in Figure 11) and Equation (4) (red dotted lines—as in Figure 12). The blue curve is the intersection between the black and red lines for each value of $\theta$ and the uncertainties (dotted blue lines—which depend on the determination of $\beta'$ and $\Delta m$). The orange line is the value for $\beta = \gamma$ while the gray line represents the relation for $\beta$ and $\gamma$ for the Jacobi shape. For $\theta$ between 55° and 80°, defined by the black and red full lines, there is no intersection between the blue and gray lines, meaning that there is no Jacobi solution for VS2.</figcaption></figure>

$a = b = 313.8 \pm 7.1$ km and $254.8 \le c \le 313.8$ km). This assumption automatically imposes that the light curve amplitude variation is only due to albedo features on VS2's surface. Considering the light curve amplitude of $0.141 \pm 0.009$ mag—see Table 1—the presence of surface irregularities (lumps) and some albedo variation (spots) on the object with a size approximately

of 100 km (~16% of the object equivalent area) is needed. Note that the *Sputnik Planitia* in Pluto is nearly 1000 × 800 km across (~15% of its equivalent area; Hamilton et al. 2016). Considering the Maclaurin shape and using the rotational period ($P = 7.4175285 \pm 0.00001$ hr; see Table 1), we can derive VS2's density of $\rho = 1.4_{-0.3}^{+1.0}$ g cm$^{-3}$, as shown in Figure 14.
<figure><img src="image_17.png" /><figcaption>**Figure 14.** Relation between $\epsilon$ and density for a Maclaurin object with rotation period of 7.4175285 hr (black line). The oblateness obtained for VS2 (blue horizontal full and dotted lines) gives limits for the minimum value of the density of $1.4 \text{ g cm}^{-3}$. Blue dotted lines represent the values considering the uncertainties in *a* and *b*. The red line is the interval between the limits.</figcaption></figure>

## ORCID iDs

Gustavo Benedetti-Rossi https://orcid.org/0000-0002-4106-476X

P. Santos-Sanz https://orcid.org/0000-0002-1123-983X

M. Assafin https://orcid.org/0000-0002-8211-0777

B. Sicardy https://orcid.org/0000-0003-1995-0842

R. Duffard https://orcid.org/0000-0001-5963-5850

F. Braga-Ribas https://orcid.org/0000-0003-2311-2438

F. Vachier https://orcid.org/0000-0002-4289-4466

E. Fernández-Valenzuela https://orcid.org/0000-0003-2132-7769

D. Bérard https://orcid.org/0000-0001-7654-6809

F. B. Bianco https://orcid.org/0000-0003-1953-8727

V. D. Ivanov https://orcid.org/0000-0002-5963-1283

M. Emilio https://orcid.org/0000-0001-5589-9015

## References

Alvarez-Candal, A., Pinilla-Alonso, N., Ortiz, J. L., et al. 2016, A&A, 586, A155

Anglada, G., Amado, P. J., Ortiz, J. L., et al. 2017, ApJL, 850, L6

Assafin, M., Camargo, J. I. B., Vieira Martins, R., et al. 2010, A&A, 515, A32

Assafin, M., Camargo, J. I. B., Vieira Martins, R., et al. 2012, A&A, 541, A142

Assafin, M., Vieira-Martins, R., Camargo, J. I. B., et al. 2011, in GAIA FUN-SSO Workshop Proc., Gaia Follow-up Network for the Solar System Objects, ed. P. Tanga & W. Thuillot (Paris: Paris Observatory), 85

Barkume, K. M., Brown, M. E., & Schaller, E. L. 2008, AJ, 135, 55

Barucci, M. A., Dotto, E., & Levasseur-Regourd, A. C. 2011, A&ARv, 19, 48

Belskaya, I. N., Ortiz, J. L., Rousselot, P., et al. 2006, Icar, 184, 277

Benedetti-Rossi, G., Sicardy, B., Buie, M. W., et al. 2016, AJ, 152, 156

Bérard, D., Sicardy, B., Camargo, J. I. B., et al. 2017, AJ, 154, 144

Binzel, R. P., Farinella, P., Zappala, V., & Cellino, A. 1989, in Asteroids II, ed. R. P., Binzel et al. (Tucson, AZ: Univ. Arizona Press), 416

Braga-Ribas, F., Sicardy, B., Ortiz, J. L., et al. 2013, ApJ, 773, 26

Braga-Ribas, F., Sicardy, B., Ortiz, J. L., et al. 2014, Natur, 508, 72

Brandeker, A., & Cataldi, G. 2019, A&A, 621, A86

Buie, M. W., & Keller, J. M. 2016, AJ, 151, 73

Camargo, J. I. B., Desmars, J., Braga-Ribas, F., et al. 2018, P&SS, 154, 59

Camargo, J. I. B., Vieira-Martins, R., Assafin, M., et al. 2014, A&A, 561, A37

Chandrasekhar, S. 1969, Ellipsoidal Figures of Equilibrium (New Haven, CT: Yale Univ. Press)

Chandrasekhar, S. 1987, Ellipsoidal Figures of Equilibrium (New York: Dover)

Deeths, D., & Brunette, G. 2001, Using NTP to Control and Synchronize System Clocks Part I: Introduction to NTP (Palo Alto, CA: Sun

Microsystems, Inc.), http://www-it.desy.de/common/documentation/cd-docs/sun/blueprints/0701/NTP.pdf

Dias-Oliveira, A., Sicardy, B., Ortiz, J. L., et al. 2017, AJ, 154, 22

Duffard, R., Ortiz, J. L., Santos Sanz, P., et al. 2008, A&A, 479, 877

Fernández-Valenzuela, E., Ortiz, J. L., Duffard, R., et al. 2016, MNRAS, 456, 2354

Fernández-Valenzuela, E., Ortiz, J. L., Duffard, R., et al. 2017, MNRAS, 466, 4147

Gaia Collaboration, Brown, A. G. A., Vallenari, A., et al. 2016a, A&A, 595, A2

Gaia Collaboration, Brown, A. G. A., & Vallenari, A. 2018, A&A, 616, A1

Gaia Collaboration, Prusti, T., de Bruijne, J. H. J., et al. 2016b, A&A, 595, A1

Gladman, B., Marsden, B. G., & Vanlaerhoven, C. 2008, in The Solar System Beyond Neptune, ed. M. A. Barucci et al. (Tucson, AZ: Univ. Arizona Press), 43

Guilbert, A., Alvarez-Candal, A., Merlin, F., et al. 2009, Icar, 201, 272

Hamilton, Douglas P., Stern, S. A., Moore, J. M., et al. 2016, Natur, 540, 97

Hinson, D. P., Linscott, I. R., Young, L. A., et al. 2017, Icar, 290, 96

Jewitt, D., & Luu, J. 1993, Natur, 362, 6422

Jewitt, D., Morbidelli, A., & Rauer, H. 2008, Trans-Neptunian Objects and Comets (Berlin: Springer)

Lacerda, P., & Jewitt, D. C. 2007, AJ, 133, 1393

Leiva, R., Sicardy, B., Camargo, J. I. B., et al. 2017, AJ, 154, 159

Lellouch, E., Hubbard, W. B., Sicardy, B., et al. 1986, Natur, 324, 227

Lellouch, E., Moreno, R., Ortiz, J. L., et al. 2002, A&A, 391, 1133

Lellouch, E., Santos-Sanz, P., Lacerda, P., et al. 2013, A&A, 557, A60

Lykawka, P. S., & Mukai, T. 2008, AJ, 135, 1161

Minor Planet Center Electronic Circular 2006, MPEC 2006-X45: DISTANT MINOR PLANETS (2006 DEC. 21.0 TT), https://www.minorplanetcenter.net/mpec/K06/K06X45.html

Mommert, M., Harris, A. W., Kiss, C., et al. 2012, A&A, 541, A93

Moorwood, A., Cuby, J.-G., Biereichel, P., et al. 1998a, Msngr, 94, 7

Moorwood, A., Cuby, J.-G., & Lidman, C. 1998b, Msngr, 91, 9

Morbidelli, A., Levison, H. F., & Gomes, R. 2008, in The Solar System Beyond Neptune, ed. M. A. Barucci et al. (Tucson: Univ. Arizona Press), 275

Ortiz, J. L., Baumont, S., Gutiérrez, P. J., et al. 2002, A&A, 388, 661

Ortiz, J. L., Cikota, A., Cikota, S., et al. 2011, A&A, 525, A31

Ortiz, J. L., Duffard, R., Pinilla-Alonso, N., et al. 2015, A&A, 576, A18

Ortiz, J. L., Gutiérrez, P. J., Santos-Sanz, P., et al. 2006, A&A, 447, 1131

Ortiz, J. L., Santos-Sanz, P., Gutiérrez, P. J., et al. 2007, A&A, 468, L13

Ortiz, J. L., Santos-Sanz, P., Sicardy, B., et al. 2017, Natur, 550, 219

Ortiz, J. L., Sicardy, B., Braga-Ribas, F., et al. 2012, Natur, 491, 566

Ortiz, J. L., Sota, A., Moreno, R., et al. 2004, A&A, 420, 383

Parker, A. H. 2015, Icar, 247, 112

Perna, D., Barucci, M. A., Fornasier, S., et al. 2010, A&A, 510, A53

Santos-Sanz, P., French, R. G., Pinilla-Alonso, N., et al. 2016, PASP, 128, 018011

Santos-Sanz, P., Lellouch, E., Groussin, O., et al. 2017, A&A, 604, A95
Sheppard, S. S. 2007, AJ, 134, 787

Sicardy, B., Ortiz, J. L., Assafin, M., et al. 2011, Natur, 478, 493

Stansberry, J., Grundy, W., Brown, M., et al. 2008, in The Solar System
Beyond Neptune, ed. M. A. Barucci et al. (Tucson, AZ: Univ. Arizona
Press), 161

Stern, A., Bagenal, F., Ennico, K., et al. 2015, Sci, 350, aad1815

Stern, S. A., Weaver, H. A., Spencer, J. R., et al. 2019, Sci, 364, eaaw9771

Stetson, Peter B. 1987, PASP, 99, 191

Tancredi, G., & Favre, S. 2008, Icar, 195, 851

Thirouin, A. 2013, Study of Trans-Neptunian Objects using Photometric
Techniques and Numerical Simulations (Granada: Universidad de Granada)

Thirouin, A., Ortiz, J. L., Duffard, R., et al. 2010, A&A, 522, A93

van Belle, G. T. 1999, PASP, 111, 1515

van Leeuwen, F., de Bruijne, J. H. J., & Arenou, F. 2017, Gaia DR1 Documentation, European Space Agency; Gaia Data Processing and Analysis Consortium, https://gaia.esac.esa.int/documentation/GDR1/index.html

Widemann, T., Sicardy, B., Dusser, R., et al. 2009, Icar, 199, 458

Zacharias, N., Monet, D. G., Levine, S. E., et al. 2004, AAS Meeting, 205, 48.15