# The Changing Rotational Light-curve Amplitude of Varuna and Evidence for a Close-in Satellite

Estela Fernández-Valenzuela, Jose Luis Ortiz, Nicolás Morales, Pablo Santos-Sanz, René Duffard, Amadeo Aznar, Vania Lorenzi, Noemí Pinilla-Alonso, Emmanuel Lellouch

► To cite this version:

Estela Fernández-Valenzuela, Jose Luis Ortiz, Nicolás Morales, Pablo Santos-Sanz, René Duffard, et al.. The Changing Rotational Light-curve Amplitude of Varuna and Evidence for a Close-in Satellite. The Astrophysical Journal Letters, 2019, 883 (1), pp.L21. 10.3847/2041-8213/ab40c2. hal-03131439

HAL Id: hal-03131439

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

Submitted on 25 Jun 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
<figure><img src="image_2.png" /></figure>

# The Changing Rotational Light-curve Amplitude of Varuna and Evidence for a Close-in Satellite

Estela Fernández-Valenzuela<fnref n="1" /> ID, Jose Luis Ortiz<fnref n="2" />, Nicolás Morales<fnref n="2" />, Pablo Santos-Sanz<fnref n="2" /> ID, René Duffard<fnref n="2" /> ID, Amadeo Aznar<fnref n="3" />, Vania Lorenzi<fnref n="4" />,<fnref n="5" />, Noemí Pinilla-Alonso<fnref n="1" />, and Emmanuel Lellouch<fnref n="6" /> ID

Received 2019 June 30; revised 2019 August 29; accepted 2019 August 29; published 2019 September 23

## Abstract

From CCD observations carried out with different telescopes, we present short-term photometric measurements of the large trans-Neptunian object Varuna in 10 epochs, spanning around 19 years. We observe that the amplitude of the rotational light curve has changed considerably during this period of time from 0.41 to 0.55 mag. In order to explain this variation, we constructed a model in which Varuna has a simple triaxial shape, assuming that the main effect comes from the change of the aspect angle as seen from Earth, due to Varuna's orbital motion in the 19 year time span. The best fits to the data correspond to a family of solutions with axial ratios *b*/*a* between 0.56 and 0.60. This constrains the pole orientation in two different ranges of solutions presented here as maps. Apart from the remarkable variation of the amplitude, we have detected changes in the overall shape of the rotational light curve over shorter timescales. After the analysis of the periodogram of the residuals to a 6.343572 hr double-peaked rotational light-curve fit, we find a clear additional periodicity. We propose that these changes in the rotational light-curve shape are due to a large and close-in satellite whose rotation induces the additional periodicity. The peak-to-valley amplitude of this oscillation is in the order of 0.04 mag. We estimate that the satellite orbits Varuna with a period of 11.9819 hr (or 23.9638 hr), assuming that the satellite is tidally locked, at a distance of ~1300 km (or ~2000 km) from Varuna, outside the Roche limit.

*Unified Astronomy Thesaurus concepts: Small solar system bodies (1469); Solar system formation (1530); CCD photometry (208); Trans-Neptunian objects (1705); Observational astronomy (1145); Optical observation (1169)*

*Supporting material: machine-readable table*

## 1. Introduction

Trans-Neptunian objects (TNOs) are solar system bodies that orbit the Sun with larger semimajor axes than that of Neptune, formed quite outside the so-called “snow line” where the temperature of the protoplanetary disk was low enough to allow the survival of molecules of chemical compounds with low sublimation points. Due to the ample distances that separate TNOs from the Sun, they have suffered less chemical processes than other solar system bodies and may preserve material from the nebula that generated the solar system. Hence, TNOs yield important information about the solar system formation and its evolution.

Currently, observational searches have found ~2500 TNOs (according to the data from MPC<fnref n="7" />), although dynamical models suggest that the TNO population located between 30 and 50 au from the Sun is around 10⁵ objects with diameters larger than 100 km (Trujillo et al. 2001; Petit et al. 2008, 2011).

20000 Varuna is one of the most interesting TNOs due to its peculiar physical properties. It rotates relatively fast with a period of 6.34358 ± 0.00005 hr (Belskaya et al. 2006; Santos-Sanz et al. 2013), producing a double-peaked rotational light

curve dominated by its shape (e.g., Jewitt & Sheppard 2002; Lellouch et al. 2002; Hicks et al. 2005; Belskaya et al. 2006). The rotational light-curve amplitude reported in these studies is quite large, ~0.45 mag, indicating that the body is highly elongated (i.e., *b*/*a* < 0.66, being *a* and *b* the larger semiaxes of a triaxial body). The fast period and the elongated shape require a density of ~1000 kg m⁻<fnref n="3" /> (assuming hydrostatic equilibrium; Chandrasekhar 1987). This is somewhat high compared to other TNOs of similar sizes, considering that Varuna's equivalent diameter is ~700 km, given by thermo-physical models using Herschel Space Telescope measurements (Lellouch et al. 2013); see the supplementary material in Ortiz et al. (2012) to compare the density of objects with similar sizes. However, a stellar occultation by Varuna, detected in 2010, results in a long chord of 1004 km (Sicardy et al. 2010). This value is somewhat in tension with the one given by Lellouch et al. (2013), but fits better with Varuna's density.

Motivated by all of this, we have been monitoring Varuna for nearly two decades now, resulting in a collection of data with several rotational light curves (which we will simply call light curves in the following). In this work, we present new observations of Varuna at different epochs (Section 2). The results from these observations, together with other observations from the literature, are presented in Section 3. Section 4 contains the analysis of the set of light curves in order to obtain Varuna's pole orientation and shape. Section 5 displays the

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

<footnote n="1">¹ Florida Space Institute, 12354 Research Parkway, Partnership 1, Room 211, Orlando, FL 32826-0650, USA; estela@ucf.edu, ortiz@iaa.es</footnote>

<footnote n="2">² Instituto de Astrofísica de Andalucía (IAA-CSIC), Glorieta de la Astronomía s/n, Granada, E-18008, Spain</footnote>

<footnote n="3">³ Observatorio Isaac Aznar, Grupo de Observatorios APT, C/La Plana, 44, 13, E-46530 Puçol, Valencia, Spain</footnote>

<footnote n="4">⁴ Fundación Galileo Galilei—Istituto Nazionale di Astrofisica, Rambla José Ana Fernández Pérez, 7, E-38712 Breña Baja, TF, Spain</footnote>

<footnote n="5">⁵ Instituto de Astrofísica de Canarias, C/Va Láctea s/n, E-38205 La Laguna, Spain</footnote>

<footnote n="6">⁶ LESIA, Observatoire de Paris, Université PSL, CNRS, Univ. Paris Diderot, Sorbonne Paris, 1Cité, Sorbonne Université, 5 Place J. Janssen, F-92195 Meudon
Pricipal Cedex, France</footnote>

<footnote n="7">7 https://minorplanetcenter.net</footnote>

<footnote>Original content from this work may be used under the terms of the Creative Commons Attribution 3.0 licence. Any further distribution of this work must maintain attribution to the author(s) and the title of the work, journal citation and DOI.</footnote>
<table><caption>**Table 1**

Varunas Light-curve Amplitudes from Different Epochs and Log of Time Series Observations</caption><thead><tr><th>Year</th><th>Δm<br>(mag)</th><th>σ<br>(mag)</th><th>Tel.</th><th>Seeing<br>(′′)</th><th>N<sub>im</sub></th><th>Reference</th></tr></thead><tbody><tr><td>2001</td><td>0.41 ± 0.01</td><td>0.037</td><td></td><td></td><td></td><td>Jewitt & Sheppard (2002)</td></tr><tr><td>2002A</td><td>0.45 ± 0.01</td><td>0.057</td><td></td><td></td><td></td><td>Lellouch et al. (2002)</td></tr><tr><td>2002B</td><td>0.47 ± 0.02</td><td>0.064</td><td></td><td></td><td></td><td>Hicks et al. (2005)</td></tr><tr><td>2005</td><td>0.46 ± 0.02</td><td>0.036</td><td>1.5 m OSN</td><td>2.01</td><td>75</td><td>This work</td></tr><tr><td>2009</td><td>0.50 ± 0.01</td><td>0.024</td><td>1.5 m OSN</td><td>1.85</td><td>48</td><td>This work</td></tr><tr><td></td><td></td><td></td><td>2.2 m CAHA</td><td>1.40</td><td>28</td><td></td></tr><tr><td>2011</td><td>0.50 ± 0.01</td><td>0.042</td><td>1.5 m OSN</td><td>2.8</td><td>120</td><td>This work</td></tr><tr><td></td><td></td><td></td><td>2.2 m CAHA</td><td>1.8</td><td>38</td><td></td></tr><tr><td>2012</td><td>0.49 ± 0.02</td><td>0.040</td><td>1.5 m OSN</td><td>2.6</td><td>64</td><td>This work</td></tr><tr><td>2014A</td><td>0.543 ± 0.009</td><td>0.030</td><td>TNG</td><td>1.6</td><td>119</td><td>This work</td></tr><tr><td>2014B</td><td>0.57 ± 0.02</td><td>0.047</td><td>1.23 m CAHA</td><td>1.7</td><td>116</td><td>This work</td></tr><tr><td>2018</td><td>0.545 ± 0.007</td><td>0.052</td><td>1.23 m CAHA</td><td>1.8</td><td>135</td><td>This work</td></tr><tr><td></td><td></td><td></td><td>1.5 m OSN</td><td>2.1</td><td>255</td><td></td></tr><tr><td>2019</td><td>0.550 ± 0.008</td><td>0.047</td><td>2.2 m CAHA</td><td>1.5</td><td>119</td><td>This work</td></tr><tr><td></td><td></td><td></td><td>1.5 m OSN</td><td>3.9</td><td>64</td><td></td></tr><tr><td></td><td></td><td></td><td>1.23 m CAHA</td><td>1.7</td><td>75</td><td></td></tr></tbody></table>

Note. Abbreviations are defined as follows: light-curve amplitude (Δm); dispersion of the residual values of the second-order Fourier function fit to the observational data (σ), telescope used (Tel.), and number of images used (N<sub>im</sub>).

evidence for a possible satellite. A brief summary is presented in Section 6.

## 2. Observations and Data Reduction

We carried out eight observing campaigns from 2005 to 2019 using four different telescopes. A log of observations is shown in Table 1.

The 1.5 m telescope at the Sierra Nevada Observatory in Granada (Spain) has a CCD VersArray with 2048 × 2048 pixels, while the field of view (FoV) is 7.92 × 7.92 and the image scale is 0.′′232 pixel⁻¹. Observations were obtained in 2 × 2 binning mode. The typical seeing was ~1.′′5, being the point-spread function well sampled for the photometric goals.

The 1.23 m and the 2.2 m telescopes at Calar Alto Observatory are located in Almería (Spain). The 1.23 m telescope uses a DLR-MKIII camera with 4096 × 4096 pixels, while the FoV and the image scale are 21.5 × 21.5 and 0.′′315 pixel⁻¹, respectively. At the 2.2 m telescope we used two different instruments: the Bonn University Simultaneous CAmera (BUSCA) and the Calar Alto Faint Object Spectrograph (CAFOS). BUSCA allows the simultaneous direct imaging of the same sky area in four colors; each CCD has 4096 × 4096 pixels, with an FoV of 12′ × 12′ and an image scale of 0.′′176 pixel⁻¹. We used 2 × 2 binning mode. CAFOS is equipped with a 2048 × 2048 pixel CCD; we used the SITe#1d chip, which produces a 0.′′53 pixel⁻¹ image scale and a circular FoV of 16′.

The Telescopio Nazionale Galileo (TNG), with 3.58 m diameter, is located at the Roque de los Muchachos Observatory in Canary Islands (Spain). We used the Device Optimized for the LOw RESolution (DOLORES) with a detector of 2048 × 2048 pixels. The image scale is 0.′′252 pixel⁻¹, which yields an FoV of 8.6 × 8.6.

Bias and twilight flat-field calibration images were taken at the beginning of each observation night. The science images were processed by subtracting a median bias and dividing by a median flat-field. Using the Interactive Data Language, we developed our own code to perform these reduction steps, as

well as the synthetic aperture photometry for the time series data.

We performed the synthetic aperture photometry using different aperture radii, as described in Fernández-Valenzuela et al. (2016), and choosing the one in which the dispersion of the photometry is minimal. To construct the light curves, Varuna was compared to a set of reference stars located as close as possible to the object. Whenever possible during each campaign, we employed the same set of comparison stars each night in order to minimize systematic photometric errors.

## 3. Results from Observations

From the time series observations obtained between 2005 and 2019, we have built eight different light curves. Additionally, we used data published in the literature by Jewitt & Sheppard (2002), Lellouch et al. (2002), and Hicks et al. (2005), thus incorporating three additional light curves to our study from previous years (Figure 1). All the light curves were corrected from light travel time. Because Varuna's body is assumed to have an ellipsoidal shape (e.g., Jewitt & Sheppard 2002; Lellouch et al. 2002), data from each light curve were fitted to a Fourier series $m = \sum_i [a_i \sin(2i\pi\phi) + b_i \cos(2i\pi\phi)]$, where m is the theoretical value of the relative magnitude obtained from the fit, φ is the rotational phase (calculated as the fractional part of $(JD - JD_0)/P$, where JD is the Julian Date, $JD_0 = 2451957.0$ is the initial Julian Date, and P is the rotation period in days), and $(a_i, b_i)$ are the coefficients of the Fourier function (with $i = 0, 1, 2, ...$). In our specific case, we used up to second-order ($i = 2$) or up to fourth-order ($i = 4$) Fourier functions. The second order is the minimum order that allows a double-peaked fit; however, higher orders take into account small deviations on inhomogeneous objects and can be used to fit light curves that are highly sampled. The 2001, 2002A, 2011, 2018, and 2019 light curves were fitted to a fourth-order function, while the remaining light curves were fitted to a second-order Fourier function, because the number of data points in those runs was not large enough to use a higher order. Data were folded using
<figure><img src="image_4.png" /><figcaption>**Figure 1.** Varuna's light curves from 2001 to 2019. Each color corresponds to a different epoch. 2001, 2002A, and 2002B data were taken from Jewitt & Sheppard (2002), Lellouch et al. (2002), and Hicks et al. (2005), respectively. The lines overplotted to the data points represent the fit of the observational data to the Fourier series. Error bars are not shown to avoid a cluttered plot; the data with the uncertainties are given in Table 2.</figcaption></figure>

Varuna's rotation period of 6.343572 ± 0.000006 hr, which is obtained using the Lomb periodogram analysis of all our data, in agreement with the previous one reported in Belskaya et al. (2006). The peak-to-valley amplitudes Δ*m* (amplitudes in the following) of each light curve are given by the absolute maximum and minimum produced by the fits.

Table 1 contains the results from the fit to each light curve, i.e., the amplitude and the dispersion of the residuals of the Fourier function fit to the observational data. One of the evident results is that the amplitude has changed considerably along these 19 years, with an increase of ~0.13 mag. Table 2 presents all the relative photometry observations from 2005 to 2019.

## 4. Varuna's Pole Orientation and Shape

The aspect angle of a solar system body, formed by the line of sight and the rotational axis of the body, changes due to the orbital motion of the body around the Sun. This change produces a variation in the geometric cross section of the object's rotational modulation that is translated into a variation in the amplitude which, for a triaxial ellipsoid, is described as

$$ \Delta m = -2.5 \log \left[ \frac{b}{a} \left( \frac{a^2 \cos^2(\delta) + c^2 \sin^2(\delta)}{b^2 \cos^2(\delta) + c^2 \sin^2(\delta)} \right)^{1/2} \right], \quad (1) $$

where δ is the aspect angle and a, b, c are the semiaxes of a triaxial body (with a > b > c; Tegler et al. 2005). The aspect angle of the body can be written as

$$ \delta = \frac{\pi}{2} - \arcsin[\sin(\beta_e)\sin(\beta_p) + \cos(\beta_e)\cos(\beta_p)\cos(\lambda_e - \lambda_p)], \quad (2) $$

where λ_e and β_e are the ecliptic longitude and latitude of Varuna's-centered reference frame (given by the ephemeris),

<table><caption>**Table 2**

Results from the Relative Photometry and Observational Circumstances</caption><thead><tr><th>JD<br>(d)</th><th>Rel. mag.<br>(mag)</th><th>σ<br>(mag)</th><th>r<sub>H</sub><br>(au)</th><th>Δ<br>(au)</th><th>α<br>(deg)</th></tr></thead><tbody><tr><td>2453376.23</td><td>0.1137</td><td>0.023</td><td>43.2515</td><td>42.3639</td><td>0.5659</td></tr><tr><td>2453376.236</td><td>0.061</td><td>0.0225</td><td>43.2515</td><td>42.3639</td><td>0.5659</td></tr><tr><td>2453376.242</td><td>0.0272</td><td>0.021</td><td>43.2515</td><td>42.3639</td><td>0.5659</td></tr><tr><td>2453376.248</td><td>-0.0367</td><td>0.0197</td><td>43.2515</td><td>42.3639</td><td>0.5659</td></tr><tr><td>2453376.254</td><td>-0.1125</td><td>0.0205</td><td>43.2515</td><td>42.3639</td><td>0.5659</td></tr></tbody></table>

**Note.** Abbreviations are defined as follows: Julian date (JD), relative magnitude (Rel. mag.), error of the relative magnitude (σ), heliocentric and topocentric distances (r<sub>H</sub> and Δ, respectively), and phase angle (α).
(This table is available in its entirety in machine-readable form.)

and λ<sub>p</sub> and β<sub>p</sub> are the ecliptic longitude and latitude of the pole orientation (Schroll et al. 1976). These coordinates change with the observing epochs. In order to fit the observational data to Equation (1), we carried out a grid search for the free parameters in the equations: b/a, c/b, λ<sub>p</sub>, and β<sub>p</sub>, which gave theoretical values for Δ*m* (as done in Fernández-Valenzuela et al. 2017). In a first step, we assumed hydrostatic equilibrium, thus reducing the number of free parameters to three, as c/b is then related to b/a (Chandrasekhar 1987). We explored b/a in the range [0.44, 0.60] at intervals of 0.02. The upper limit of 0.60 is imposed by the measurement of the largest observed amplitude up to date (see Table 1). The lower limit of 0.44 is imposed by theoretical models of rotational equilibrium configurations of strengthless bodies (Leone et al. 1984). We explored the possible values of λ<sub>p</sub> and β<sub>p</sub> on the entire sky at intervals of 0.5°.
<figure><img src="image_5.png" /><figcaption>**Figure 2.** Left panel: model of the amplitude given by the fit to the observational data using $b/a = 0.60$ (pink line, best fit), 0.58 (orange line), and 0.56 (yellow line). The gray line represents the Varuna's aspect angle given by the best fit. Right panel: $\chi^2_{pdft}$ map of possible values for $b/a = 0.60$. Similar configurations occur for values of the axis ratios of 0.58 and 0.56.</figcaption></figure>

The goodness of the fit is evaluated using the $\chi^2_{pdft}$ test. The possible solutions are those that have values of $\chi^2_{pdft}$ within the range $\chi^2_{pdft, min} + 0.584$ (from the $\chi^2$ distribution with a 0.9 level of confidence for 3 degrees of freedom; Feller 1971). As a result, the subsequent range of values for the axis ratio $b/a$ is [0.56, 0.60], where the lower limit is given by the goodness of the fit. Considering the results from the different axis ratios, we constrained the pole orientation with $\lambda_p \in [43.0, 83.5]^{\circ}$ and $\beta_p \in [-70.0, -58.0]^{\circ}$ and $\lambda_p \in [209, 220]^{\circ}$ and $\beta_p \in [-49.0, -35.5]^{\circ}$ (with their complementary directions also possible).

Specifically, we obtained the $\chi^2_{pdft, min} = 1.904$ for the following parameters: semiaxis ratio $b/a = 0.60$ (with $c/b = 0.72$ imposed by the hydrostatic equilibrium), and pole orientation $\lambda_p = 53^{\circ} \pm 10^{\circ}$ and $\beta_p = -64^{\circ} \pm 6^{\circ}$ (the complementary direction $\lambda_p = 233^{\circ}$ and $\beta_p = 64^{\circ}$ is also possible for the same $\chi^2_{pdft}$ value). The errors of $\lambda_p$ and $\beta_p$ have been obtained considering the values that fit within the interval of $\chi^2_{pdft, min} + 0.584$ for this specific axis ratio. Figure 2 shows the best fit given by the above values (left panel) and a $\chi^2_{pdft}$ map as a function of $\lambda_p$ and $\beta_p$ (right panel) for the specific value of $b/a = 0.60$. The solution in the range of $\lambda_p \in [209, 211]^{\circ}$ and $\beta_p \in [-42.0, -35.5]^{\circ}$ can be discarded making use of results from the stellar occultation that occurred in 2010 (Sicardy et al. 2010). Even if a detailed analysis will be published elsewhere, we can advance that this yields restricted intervals for the values of $\lambda_p$ and $\beta_p$. It turns out that the first possibility $\lambda_p = 53^{\circ} \pm 10^{\circ}$, $\beta_p = -64^{\circ} \pm 6^{\circ}$ (as well as its complementary direction) lies inside these intervals, for which the resulting area of Varuna does not produce a geometric albedo being unnaturally small,<fnref n="8" /> but this is not the case for $\lambda_p \in [209, 211]^{\circ}$, $\beta_p \in [-42.0, -35.5]^{\circ}$, for which the size of Varuna would be too large, corresponding to an albedo well below the lower limit of 0.04 for TNOs.

In a second step, we relaxed the assumption of hydrostatic equilibrium, adding the parameter $c/b$ to the grid search in the range [0.5, 1.0] using steps with a length of 0.1. However, this does not lead to appreciable differences. The smallest value of $\chi^2_{pdft}$ found is 1.905, which results in values of $\lambda_p = 54^{\circ}$ and $\beta_p = -65^{\circ}$, with the axis ratios $b/a = 0.60$ and $c/b = 0.70$, which is the equivalent result to the case of hydrostatic equilibrium. Currently, from the best fit, Varuna is close to an “edge-on” configuration (see the left panel in Figure 2), where the aspect angle is close to 90°. Therefore, the c-axis, which is the one that would determine whether Varuna is in hydrostatic equilibrium or not, does not produce a strong effect on the light-curve shape, being difficult to distinguish between both scenarios.

There also exists the possibility that Varuna could be a contact binary, as happened with 2014 MU19 (Stern et al. 2019). However, this scenario was explored by Jewitt & Sheppard (2002) and was considered very unlikely for Varuna. Note that contact binaries have specific V-shaped light curves (Thirouin & Sheppard 2017; Thirouin et al. 2017), which is not the case of Varuna.

## 5. Indications for Varuna's Close-in Satellite

Even though the signal-to-noise ratio (S/N) of the data in most of the runs was good (well above 50), the residuals of the Fourier fits were considerably scattered, indicating that the model may have been missing some aspect. We focused on the data obtained during 2019, with the best S/N (~80) for Varuna. In Figure 3, we plot the relative magnitude of each day versus Julian date. While the four light curves present a good fit, at localized points there are deviations of the model with respect to the data. Moreover, when merging the four days, the quality of the assembled light curve decreases considerably, as the scatter of residuals increases with respect to single night fits and the parameters of the fits change slightly from night to night. A satellite might be responsible for changes in the shape of the main Varuna's light curve. In order to test the potential presence of a rotating satellite that could be generating an additional periodicity in the photometry, we used the residuals from the fit of the 2011, 2014A, 2014B, 2018, and 2019 light

<footnote n="8">⁸ The relation between the equivalent diameter and the geometric albedo is $D = Cp^{-1/2}10^{-H/5}$, where C = 1329 km is a constant, p is the geometric albedo and H is the absolute magnitude.</footnote>
<figure><img src="image_6.png" /><figcaption>**Figure 3.** Light curves observed during 2019; each panel corresponds to a different day in which the corresponding data are represented by circles. The solid lines overplotted to the circles represent the fourth-order Fourier fit to the observational data, with each day fitted independently. The x-axis represents the fraction of JD, which is given within the corresponding panel. The residuals to the fit can be seen in the bottom panels.</figcaption></figure>

<figure><img src="image_7.png" /><figcaption>**Figure 4.** Left panel: periodogram of the residuals of the fit to Varuna's observational data from years 2011, 2014, 2018, and 2019. A maximum spectral power of 45 is obtained at 2.003 cycles/day (11.9819 hr). Right panel: residuals of Varuna's observational data from years 2011, 2014, 2018, and 2019 (blue circles), folded to the 11.9819 hr period. The blue line represents a one-order Fourier function fit to the points.</figcaption></figure>

curves (which are the ones with higher quality that are better sampled) and searched for a period using the Lomb periodogram technique (Lomb 1976) as implemented in Press et al. (1992). This technique has been largely applied to detect satellites in asteroids (e.g., Pravec & Hahn 1997; Pravec et al. 2002; Margot et al. 2015). In Figure 4 we plotted the Lomb periodogram of the residuals (left panel), which gives the maximum spectral power of 45 for the frequency of 2.003 cycles/day (11.9819 hr), and the residuals folded to that period (right panel), where the blue line shows a one-order Fourier function fit. The amplitude obtained from this fit, 0.04 mag, is related to the rotational modulation of the satellite. Assuming that the satellite's light curve is due to an albedo spot and that its spin and orbit are synchronized (we will check below that this is a reasonable assumption), the distance *d* at which the

satellite orbits from Varuna is given by the equation:

$$d^3 = \frac{T_S^2 GM}{4\pi}, \qquad (3)$$

where $T_S$ is the satellite's orbital period (in seconds), $G$ is the gravitational constant and $M = M_V + M_S$ is the mass of the system. Varuna's mass, ~10²¹ kg, is estimated taking into account a density of 1100 kg m⁻³, assuming hydrostatic equilibrium in a Jacobi ellipsoid with axis ratio $b/a \in [0.60 - 0.56]$, and volume-equivalent diameter of 700 km. This results in an estimation of $d \sim 1300$ km, which places
the satellite outside the Roche limit,

$$d_R = a_V \left( 2 \frac{\rho_V}{\rho_S} \right)^{1/3}, \qquad (4)$$

where $a_V \simeq 550$ km and $\rho_V$ are Varuna's largest semiaxis and density, respectively, and $\rho_S$ is the satellite's density. Assuming that the satellite has a density of $\sim 300$ kg m$^{-3}$ (similar to 2014 MU$_{69}$'s density; Stern et al. 2019), this results in a Roche limit of 1000 km.

We should note that the satellite could have an elongated shape with a double-peaked light curve. This implies that the rotation period of the satellite would be double ($\sim 24$ hr), and so would its orbital period. This would locate the satellite $\sim 2000$ km from Varuna, with no strong implications for the following calculations. The angular distance between Varuna and the satellite is then smaller than the resolution of the Wide Field Camera 3 installed at the Hubble Space Telescope, thus making it impossible to resolve the system.

The time required by the satellite to be tidally locked, $t_{\text{lock}}$, can be obtained using the equation given by Hubbard (1984),

$$t_{\text{lock}} = \frac{2\pi M_V d^6}{3k_V GM_S^2 a_V^3 T_{V0} \delta}, \qquad (5)$$

where $M_V$ is Varuna's mass, $k_V$ is the secular Lomb number, which has a value of $3/2$ for homogeneous bodies, $M_S$ is the mass of the satellite, $T_{V0}$ is the initial rotation rate of Varuna, and $\delta$ is expressed as arctan$(1/Q)$ with $Q$ is the dissipation. Assuming a dissipation of $Q = 100$, and that $T_{V0}$ was greater than or equal to $T_V$ currently measured, it is straightforward to calculate the dimensionless number $t_{\text{lock}}/t_{\text{ss}}$ with $t_{\text{ss}} \simeq 1.5 \times 10^{17}$ s the age of the solar system, obtaining $t_{\text{lock}}/t_{\text{ss}} \lesssim 5 \times 10^9$ kg $\times$ ($M_V/M_S^2$). Using Varuna's largest semiaxis and its density, it follows that, as long as $M_S/M_V \gg 10^{-6}$, then the dimensionless number $t_{\text{lock}}/t_{\text{ss}}$, is small. This validates the assumption that the satellite is tidally locked, as the time required by the satellite to enter this configuration is much smaller than the age of the solar system.

The value obtained for the pole orientation (Section 4) places the body almost in an “edge-on” configuration, implying that mutual events produced by the satellite may occur in the coming years, assuming that the satellite’s orbit is in Varuna’s equatorial plane.

## 6. Conclusions

Time-series photometry of Varuna at different epochs have revealed a remarkable change in the amplitude of the light curve of $\sim 0.13$ mag. A geometric model that explains this change as the rotational variability due to Varuna's shape provides a range of solutions for its pole orientation of $\lambda_p \in [43.0, 83.5]^{\circ}$ and $\beta_p \in [-70.0, -58.0]^{\circ}$ and its axes ratio $b/a \in [0.56, 0.60]$. Nowadays, only around four TNOs have well determined pole orientations (Pluto, Charon, Haumea, and 2014 MU$_{69}$), with three of them visited by spacecraft. Additionally, we noticed that Varuna's light curves have deviations from the expected behavior of a simple ellipsoidal body, and we analyzed the periodogram of the residuals of the photometric fits to the data, finding a peak of high spectral power for a period of 11.9819 hr. This appears consistent with the idea of the existence of a close-in satellite, separated from

Varuna a distance of $\sim 1300$ km (depending on Varuna's mass, which is not accurately known). Until now, satellite discoveries in the trans-Neptunian region have been only accomplished using direct imaging (Noll et al. 2008). With this work, we have illustrated that photometric techniques can be used to reveal the existence of a satellite for an object in the trans-Neptunian region. The light-curve technique has the potential to reveal objects much closer-in than the direct imaging technique, helping to remove biases in our current knowledge of the population of binaries or bodies with satellites.

E.F.V. acknowledges UCF 2017 Preeminent Postdoctoral Program (P$^3$). Part of the research leading to these results has received funding from the European Unions Horizon 2020 Research and Innovation Programme, under Grant Agreement No. 687378 (SBNAF). We would like to acknowledge financial support by the Spanish grant AYA-2017-84637-R and the 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). This work is partially based on observations made with the Italian Telescopio Nazionale Galileo (TNG) operated on the island of La Palma by the Fundación Galileo Galilei of the INAF (Istituto Nazionale di Astrofisica) at the Spanish Observatorio del Roque de los Muchachos of the Instituto de Astrofísica de Canarias. We thank the director of the Telescopio Nazionale Galileo for allocation of Director's Discretionary Time. We acknowledge Raúl Carballo-Rubio for providing useful comments that helped improve this manuscript and Pedro Gutiérrez for helpful discussions.

## ORCID iDs

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

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

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

Emmanuel Lellouch https://orcid.org/0000-0001-7168-1577

## References

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

* Chandrasekhar, S. 1987, (New York: Dover), 1987

* Feller, W. 1971, Wiley Series in Probability and Mathematical Statistics (New York: Wiley)

* Fernández-Valenzuela, E., Ortiz, J. L., Duffard, R., Morales, N., & Santos-Sanz, P. 2017, MNRAS, 466, 4147

* Fernández-Valenzuela, E., Ortiz, J. L., Duffard, R., Santos-Sanz, P., & Morales, N. 2016, MNRAS, 456, 2354

* Hicks, M. D., Simonelli, D. P., & Buratti, B. J. 2005, Icar, 176, 492

* Hubbard, W. B. 1984, Planetary Interiors (New York: Van Nostrand Reinhold)

* Jewitt, D. C., & Sheppard, S. S. 2002, AJ, 123, 2110

* 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

* Leone, G., Paolicchi, P., Farinella, P., & Zappala, V. 1984, A&A, 140, 265

* Lomb, N. R. 1976, Ap&SS, 39, 447

* Margot, J.-L., Pravec, P., Taylor, P., Carry, B., & Jacobson, S. 2015, in Asteroids IV, ed. P. Michel et al. (Tucson, AZ: Univ. Arizona Press), 355

* Noll, K. S., Grundy, W. M., Chiang, E. I., Margot, J.-L., & Kern, S. D. 2008, in The Solar System Beyond Neptune, ed. M. A. Barucci et al. (Tucson, AZ: Univ. Arizona Press), 345

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

* Petit, J.-M., Kavelaars, J. J., Gladman, B., & Loredo, T. 2008, in The Solar System Beyond Neptune, ed. M. A. Barucci et al. (Tucson, AZ: Univ. Arizona Press), 71

* Petit, J.-M., Kavelaars, J. J., Gladman, B. J., et al. 2011, AJ, 142, 131

* Pravec, P., & Hahn, G. 1997, Icar, 127, 431
Pravec, P., Šarounová, L., Hicks, M. D., et al. 2002, Icar, 158, 276

Press, W. H., Teukolsky, S. A., Vetterling, W. T., & Flannery, B. P. 1992,
Numerical recipes in FORTRAN. The Art of Scientific Computing (2nd
ed.; Cambridge: Cambridge University Press)

Santos-Sanz, P., Ortiz, J. L., Lellouch, E., et al. 2013, EPSC, 8, 665

Schroll, A., Haupt, H. F., & Maitzen, H. M. 1976, Icar, 27, 147

Sicardy, B., Colas, F., Maquet, L., et al. 2010, AAS/DPS Meeting, 42, 23.11

Stern, S. A., Spencer, J. R., Weaver, H. A., et al. 2019, arXiv:1901.02578

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

Tegler, S. C., Romanishin, W., Consolmagno, G. J., et al. 2005, Icar, 175, 390

Thirouin, A., & Sheppard, S. S. 2017, AJ, 154, 241

Thirouin, A., Sheppard, S. S., & Noll, K. S. 2017, ApJ, 844, 135

Trujillo, C. A., Jewitt, D. C., & Luu, J. X. 2001, AJ, 122, 457