# THE GEOMETRIC ALBEDO OF (4179) TOUTATIS ESTIMATED FROM KMTNET DEEP-SOUTH OBSERVATIONS

^{1}; MASATERU ISHIGURO

^{1}; SUNHO JIN

^{1}; HONGU YANG

^{2}; HONG-KYU MOON

^{2}; YOUNG-JUN CHOI

^{2}; YOUNGMIN JEONGAHN

^{2}; MYUNG-JIN KIM

^{2}; SUNGWON KWAK

^{2}

2Korea Astronomy and Space Science Institute, 776, Daedeok-daero, Yuseong-gu, Daejeon 34055

Published under Creative Commons license CC BY-SA 4.0

## Abstract

We derive the geometric albedo of a near-Earth asteroid, (4179) Toutatis, to investigate its surface physical conditions. The asteroid has been studied rigorously not only via ground-based photometric, spectrometric, polarimetric, and radar observations but also via *in situ* observation by the Chinese Chang'e-2 space probe; however, its geometric albedo is not well understood. We conducted V-band photometric observations when the asteroid was at opposition in April 2018 using the three telescopes in the southern hemisphere that compose the Korea Microlensing Telescope Network (KMTNet). The observed time-variable cross section was corrected using the radar shape model. We find that Toutatis has a geometric albedo *p*_{V} = $$ {0.185}_{-0.039}^{+0.045}$$, which is typical of S-type asteroids. We compare the geometric albedo with archival polarimetric data and further find that the polarimetric slope−albedo law provides a reliable estimate for the albedo of this S-type asteroid. The thermal infrared observation also produced similar results if the size of the asteroid is updated to match the results from Chang'e-2. We conjecture that the surface of Toutatis is covered with grains smaller than that of the near-Sun asteroids including (1566) Icarus and (3200) Phaethon.

## Keywords:

minor planets, asteroids: general, minor planets, asteroids: individual: (4179) Toutatis## 1. INTRODUCTION

The geometric albedo is one of the most fundamental observed quantities for characterizing the physical properties and composition on asteroids. It provides information on the surface composition, degree of space weathering, and physical conditions such as the surface particle size if polarimetric data at large phase angles (Sun−asteroid−observer angle) are available (Shkuratov & Opanasenko 1992; Dollfus 1998). Because the absolute photometric magnitude (at the opposite direction from the Sun viewed at unit heliocentric and observer distances) is related to the product of the geometric albedo and cross section of the asteroid, the accuracy of the size determination is a critical factor for deriving the geometric albedo in most cases. The sizes of asteroids have been investigated thoroughly via various techniques, namely, direct imaging of asteroids with cameras onboard spacecraft, delay-Doppler imaging (Ostro 1993), thermal−infrared observations (Tedesco et al. 2002; Mainzer et al. 2011; Usui et al. 2011), and occultation (Shevchenko & Tedesco 2006); nevertheless, absolute magnitudes are rarely studied because of the scarcity of observational opportunities at opposition. In addition, magnitude variations due to an asteroid's rotation have not been considered in many cases for the derivation, making the albedo values less certain. For these reasons, the number of reliable albedo estimates for asteroids is limited.

Here, we examine the geometric albedo of an Apollo-type near-Earth asteroid (NEA), (4179) Toutatis. The asteroid has been thoroughly studied by various observational techniques (i.e., photometry, spectroscopy, polarimetry, radar imaging, and *in situ* observation). It is classified into S-type (Lazzarin et al. 1994; Howell et al. 1994) or, more strictly, S_{k}-type (Bus & Binzel 2002) and is possibly composed of undifferentiated L-chondrites (Reddy et al. 2012). Its physical dimensions are 4.6km × 2.4km × 1.9km according to radar observations (Hudson & Ostro 1995). Toutatis is in an extremely slow, non-principal-axis rotational state (Ostro et al. 1995; Takahashi et al. 2013) having rotation and precession periods of ~ 5.38 days and ~ 7.40 days, respectively (Zhao et al. 2015). An accurate shape model was established by radar delay-Doppler imaging (Hudson et al. 2003).

More recently, in December 2012, the Chinese spacecraft Chang'e-2 captured images of this asteroid with a maximum resolution better than 3m during the flyby (Huang et al. 2013). The study revealed, e.g., a previously unknown ~ 805m basin at the end of the body, lineaments over the surface especially near the rim of it, and a sharp connection silhouette in the neck, as well as providing valuable information on the asteroid's size (4.75 × 1.95km with nominal uncertainty of 10%). Moreover, more than 50 craters larger than 36m and 30 boulders larger than 10m were identified and the chronology of craters and size distribution of boulders were investigated. From the cumulative boulder size frequency distribution analyses compared with (25143) Itokawa, it has been suggested that Toutatis may have a different preservation state or a diverse formation history (Jiang et al. 2015). The shape is considered to have been produced by a low-speed impact between two components (Hu et al. 2018). The hemispherical albedo was found (Zhao et al. 2016) to be 0.2083, 0.1269, and 0.1346 in the R, G, and B bands of the Chang'e-2 CMOS censor, respectively. Moreover, the close-orbit dynamics was used to explore Toutatis' complex gravitational field. The asteroid is also a notable for its extremely chaotic orbit in a 3.1 mean motion resonance with Jupiter and a weak 1.4 resonance with Earth (Whipple & Shelus 1993). In addition, the asteroid has a unique potential field in the close orbit (Scheeres et al. 1998). Thus, Toutatis has been central to the developments in NEA research since the 1990s.

Although Toutatis has been studied intensively, its geometric albedo is uncertain by a factor of > 3, ranging from 0.13 (Lupishko et al. 1995) to 0.41 (Masiero et al. 2017). In this paper, we aim to establish a better estimate of the albedo value. Taking advantage of the continuous observation capability from KMTNet observatories in three different continents (Kim et al. 2016, see Section 2), we derived the geometric albedo of the slowly rotating asteroid. Together with the polarimetric properties given in Lupishko et al. (1995); Mukai et al. (1997); Ishiguro et al. (1997), we discuss the physical properties of the surface.

For future reference, Table 1 lists the some important symbols we used in Sections 2, 3, and Appendix A, their meanings, and their values and units. Most source codes of the software we used for the data reduction and analysis are publicly available (see Appendix B).

## 2. OBSERVATION AND DATA REDUCTION

Our observation journal is provided in Table 2. We used three 1.6-m telescopes in the southern hemisphere that together comprise the Korea Microlensing Telescope Network (or KMTNet, Kim et al. 2016) within the frame of the DEEP−South project (Moon et al. 2016). Each telescope provides a 2° × 2° field of view via identical mosaic CCD cameras with a pixel size of 0.4''. The telescopes are located on three continents: at the Cerro Tololo Inter−American Observatory (KMTNet-CTIO) in Chile, the South African Astronomical Observatory (KMTNet-SAAO) in South Africa, and the Siding Spring Observatory (KMTNet-SSO) in Australia. We employed the standard Johnson−Cousins *V* -band filter system. We obtained ~ 10 frames per night at each observatory with individual exposure times of 1 minute. The images were preprocessed by a reduction pipeline for KMTNet data at the KMTNet data center, which performs bias subtraction, flat-fielding, and cross-talk correction. A summary of the data reduction process is given below (see also Figure 1).

The locations of Toutatis were specified using the Moving Object Detection Program (MODP) developed by KASI. MODP updates the World Coordinate System (WCS) information in the header of each image first and specifies the position of the asteroid in each image. For convenience, we analyzed images clipped to 2000×2000 pixels (800''×800'') centered on the Toutatis locations calculated by the MODP. Cosmic-ray signals were removed with **astroscrappy**^{1} which uses the **L. A. Cosmic** algorithm (van Dokkum 2001). We coadded all images taken in the same night at the same observatory while shifting the individual images to increase the signal-to-noise (S/N) ratio. Eventually, two composite images per night per observatory were generated: one co-added image with sidereal offset and one by centering the individual images on Toutatis.

The signals from sources detected in each composite image are compared with those of field stars in the Pan-STARRS Data Release 1 (PS1, Flewelling et al. 2016). We first cone-searched PS1 objects within 0.06° centered at Toutatis and with more than 5 observations in the *g* and *r* bands. We then used the criterion **iMeanPSFMag**–**iMeanKronMag** > 0.05 to reject non-stellar objects,^{2} where **iMeanPSFMag** and **iMeanKronMag** are the mean i-band magnitudes obtained from using point spread function and Kron radius photometry, respectively. If saturated stars resulted in blooming and/or contaminated nearby pixels, all affected pixels were manually masked, and stars located closer than ~ 6*d*, where *d* is the full-width at half-maximum of the seeing disk, were rejected. We fixed *d* to be 4 pixels (~ 1.6''). In addition, nearly colocated stars were identified and rejected by using the DAOGROUP algorithm (Stetson 1987 implemented via **photutils** by Bradley et al. 2017) with a minimum separation distance of 7.5*d*.

Next, we used the transformation formula

$$$V=g-0.59\left(g-r\right)-0.01$$$ | (1) |

from Lupton (2005)^{3} to obtain the Johnson−Cousins *V* -band magnitudes of the stars, where *g* and *r* are the mean aperture photometry magnitudes (**MeanApMag**) of the *g*-and *r*-band, respectively. We fit a line to *V* and the instrumental magnitude from circular aperture photometry of the stars (aperture radius 1.5*d*, inner and outer sky radii 4*d* and 6*d*, respectively) in the sidereally coadded images; only stars with 15 ≤ *V* ≤ 20 are used for the analysis. The regression line provides the amount of atmospheric extinction in relative photometry, and we used it to obtain the V-band magnitude of Toutatis in the Toutatis-centered images via extrapolation. The minimum number of stars used in the regression (in the April 7 data from CTIO) was nine. Due to the extrapolation, the uncertainty in *V*, ∆*V* (obtained from the 1-*σ* confidence interval of the regression line), must be understood as a lower limit of the true uncertainty.

After data reduction, we compared the results with the raw images and found an artificial pattern contaminating one of the datasets (2018-04-11 SAAO). For this special case, we did manual photometry (selecting appropriate sky regions by visual inspection); all other data sets were analyzed as described above. We summarize the observational quantities in Table 2 and plotted them as a function of time in Figure 2.

## 3. DERIVATION OF THE GEOMETRIC ALBEDO

### 3.1. Correction for Rotation and Distances

The observed *V* varies not only because of rotation and variation of the solar phase angle (*α*, the Sun−asteroid−observer angle) but also because of changes in the heliocentric and observer distances (*r*_{h} and *r*_{g}, respectively). We first corrected the effects of the heliocentric and observer distances by deriving the reduced magnitude *H*_{V}(*α*), which is defined as a magnitude when the asteroid is viewed at unit heliocentric and observers distance (i.e., *r*_{h} = *r*_{g} = 1 au) but at arbitrary phase angles. This parameter is given by

$$${H}_{\mathrm{V}}\left(\alpha \right)=V\left({r}_{\mathrm{h}},{r}_{\mathrm{g}},\alpha \right)-5{\mathrm{log}}_{10}\left(\frac{{r}_{\mathrm{h}}{r}_{\mathrm{g}}}{1\mathrm{a}\mathrm{u}\cdot 1\mathrm{a}\mathrm{u}}\right).$$$ | (2) |

The IAU *H*, *G* magnitude system (Bowell et al. 1989) describes the *α*-dependence of *H*_{V}(*α*) using two parameters, the absolute magnitude *H*_{V} = *H*_{V}(*α* = 0), and the slope parameter *G*:

$$${H}_{\mathrm{V}}\left(\alpha \right)={H}_{\mathrm{V}}+2.5{\mathrm{log}}_{10}\left[\left(1-G\right){\mathrm{\Phi}}_{1}\left(\alpha \right)+G{\mathrm{\Phi}}_{2}\left(\alpha \right)\right],$$$ | (3) |

where

$$$\begin{array}{c}{\mathrm{\Phi}}_{1}\left(\alpha \right)=\mathrm{exp}\left[-3.33{\mathrm{tan}}^{0.63}\left(\alpha /2\right)\right],\\ {\mathrm{\Phi}}_{2}\left(\alpha \right)=\mathrm{exp}\left[-1.87{\mathrm{tan}}^{1.22}\left(\alpha /2\right)\right].\end{array}$$$ | (4) |

Our observations were made at very small *α* values. For a spherical object, the fraction of dark area (outside the terminator) is *f*_{dark} = (1 – cos *α*)/2. Since our observations were made at *α* < 3°, *f*_{dark} < 0.07%, and a non-spherical shape, such as that of Toutatis, does not significantly increase this value. Therefore, we assume that the total projected area of the target viewed from a given direction will be the same as that of the sunlit surface.

Furthermore, the reflected intensity per projected area (which is closely related to the radiance factor; see Section 3.2) is nearly independent of the incidence and emission angles and depends only on the phase angle *α*, when the object is viewed near opposition. Thus, the fluxux reflected from a surface patch on Toutatis is proportional to its projected area, regardless of its orientation. This was first conjectured observationally by Markov & Barabashev (1926) from lunar observations as mentioned in the appendix of Dollfus & Bowell (1971). It was confirmed by Lee & Ishiguro (2018) using a modern theoretical model: their Figure 4 shows that the radiance factor (defined in Section 3.2) is reduced only by ~ 1% or ~ 3% when the incidence angle is chagned from 0° to 70° or 90° (the azimuthal angle and phase angle are fixed to 0). This means, if we correct the phase effect using Equation (3), the obtained *H*_{V} will be linearly related to the logarithm of the geometrical cross-sectional area (Equation 5 below).

Figure 2 shows *H*_{V}(*α*) and *H*_{V} as function of time, using *G* = 0.10 following Spencer et al. (1995). The best rotational model, based on the high-resolution radar shape model of Hudson et al. (2004), is shown as red solid line, and other models within the 1-*σ* confidence interval are shown as faint black dotted lines. The process of selecting rotational models is described in Appendix A. We confirmed that any *G* with 0.10 ≤ *G* ≤ 0.50 does not change our results. Our reduced magnitudes and the amplitude of the light curve agree with previously reported results (Spencer et al. 1995).

### 3.2. Radiance Factor and Geometric Albedo

The radiance factor is given by *I*/*F* where *I* is the measured irradiance of the object, whereas *F* is that of a perfectly diffuse (Bond albedo of unity) Lambertian plate^{4} of the same projected area at the identical geometric configuration to that of the object, but oriented such that the incidence angle of sunlight is 0. Since the projected area is the same for both *I* and *F*, the ratio *I*/*F* will be identical to the ratio of the bidirectional reflectance of the object (which is unknown) and the Lambertian plate (which is cos(0)/*π* = 1/*π*). Thus, at a perfect opposition (*α* = 0°), *I*/*F* coincides with *p*_{V} by definition.

As described in Section 3.1, the *I*/*F* value remains almost constant over the surface when the object is observed near opposition. Thus, the observed fluxux *F*_{obs} is proportional to *I*/*F*(*α*)×*S*_{proj}(*t*), where *I*/*F*(*α*) ≈ const over the surface of Toutatis. The obtained reduced magnitude, *H*_{V}(*α*), can then be converted into a logarithm of *I*/*F*:

$$$-2.5{\mathrm{log}}_{10}\left(\frac{I}{F}\right)={H}_{\mathrm{V}}\left(\alpha \right)-{V}_{\odot}-2.5{\mathrm{log}}_{10}\left(\frac{\pi}{{S}_{\mathrm{p}\mathrm{r}\mathrm{o}\mathrm{j}}\left(t\right)}\right)+{m}_{c},$$$ | (5) |

where *V*_{⊙} = –26.762 is the V-band magnitude of the Sun at 1 au (Mann & von Braun 2015), *S*_{proj}(*t*) is the geometrical cross section in m^{2} as a function of time, and *m _{c}* = –5 log

_{10}(1 au/1m) = –55.87 is a constant to adjust the length unit.

The calculated *I*/*F* values are plotted against the phase angle *α* in the upper panel of Figure 3. We used a single power-law (*I*/*F*(*α*) = *p*_{V} × 10* ^{bα}*) to fit the points approximately. From this plot, we considered four parameters of interest: the geometric albedo (

*p*

_{V}), the albedo at phase angle of 5° (

*A*

_{5}), and two factors

$$$\begin{array}{ccc}{f}_{1}=\frac{I\left({0.3}^{\circ}\right)}{I\left({5}^{\circ}\right)}& ;& {f}_{2}=\frac{I/F\left({0.7}^{\circ}\right)-I/F\left({2.5}^{\circ}\right)}{{2.5}^{\circ}-{0.7}^{\circ}}\end{array},$$$ | (6) |

which were used to check the validity of our results by comparing with previous works (explained below).

We used a brute-force grid search algorithm (using a fixed grid) to find the best-fit power-law parameters *p*_{V} and *b*. The best-fit parameters were first found by simple *χ*^{2} minimization where the reduced *χ*^{2} is defined as

$$${\chi}_{\mathrm{r}\mathrm{e}\mathrm{d}}^{2}\u2254\frac{1}{{N}_{\mathrm{o}\mathrm{b}\mathrm{s}}-2}\sum _{i=1}^{{N}_{\mathrm{o}\mathrm{b}\mathrm{s}}}{\left(\frac{I/F\left({\alpha}_{i}\right)-{p}_{\mathrm{V}}\times {10}^{b{\alpha}_{i}}}{{\sigma}_{i}}\right)}^{2}.$$$ | (7) |

Here, *N*_{obs} = 9 is the total number of data points, *N*_{obs} – 2 = 7 is the number of degrees of freedom, and *σ _{i}* is the uncertainty in

*I*/

*F*for the

*i*-th observation. The $$ {\chi}_{\mathrm{r}\mathrm{e}\mathrm{d}}^{2}$$ contours in the 2-D parameter space (bottom panel of Figure 3) are carefully investigated to assure that we searched a large enough parameter space with a suffciently fine grid. The uncertainties of

*f*

_{1},

*f*

_{2}, and

*A*

_{5}are estimated based on the minimum and maximum values of each of them for the models within the 1-

*σ*confidence interval (the gray region in the top panel of Figure 3).

From the fit, the amplitude parameter is the geometric albedo by definition (*p*_{V} ≡ *I*/*F*(0)). We obtained the geometric albedo in the V-band *p*_{V} = $$ {0.185}_{-0.039}^{+0.045}$$. This is a typical albedo value for an S-type asteroid.

The albedo at *α* = 5°, *A*_{5}, is important since this is used in experimental studies to compare with celestial bodies (see Section 4.4). It is defined as the ratio of the measured fluxux of the obeject of interest at *α* = 5° to that of a perfectly diffuse Lambertian plate at the same geometrical configuration with the same projected area.^{5} Thus,

$$${A}_{5}=\frac{I}{F}\left(\alpha ={5}^{\circ}\right)={0.104}_{-0.034}^{+0.050},$$$ | (8) |

as the black marker shown in Figure 3. The two factors, *f*_{1} and *f*_{2} are

$$$\begin{array}{ccc}{f}_{1}={1.708}_{-0.740}^{+1.285}& ;& {f}_{2}={0.018}_{-0.019}^{+0.017}{\u2215}^{\circ}.\end{array}$$$ | (9) |

*f*_{1} was investigated in Belskaya & Shevchenko (2000) and found to be *f*_{1} = 1.34 ± 0.05 for C-type and *f*_{1} = 1.44 ± 0.04 for S-type asteroids. From the V-filter result in Figure 7 of Tatsumi et al. (2018), *f*_{1} ~ 1.3 for the asteroid (25143) Itokawa, an S-type asteroid. *f*_{2} was determined for asteroid Itokawa in Lee & Ishiguro (2018). For Itokawa, they found *f*_{2} ~ 0.01 ± 0.01 for *α* > 1.4° and *f*_{2} ~ 0.04±0.01 for *α* < 1.4° (note that the sign convention is opposite from the one for their *s*_{OE}). We employed 0.7° and 2.5° because the phase angle difference is larger than 0.3° as Lee & Ishiguro (2018) and this is the interpolated region. Both *f*_{1} and *f*_{2} coincide with the expected values from previous works, but it is diffcult to arrive at further conclusions due to the uncertainties.

Our parameter values are summarized in Table 3. In the last column of the table, we summarized the values of the same parameters but derived by excluding one of the data points (and updating the number of degrees of freedom accordingly), 2018-04-07 CTIO with *α* = 2.79°, since this is the observation with the largest deviation from the shape model prediction (Figure 2). Comparing the two results, it can be seen that the derived parameters are robust against the exclusion of the deviant data point (parameters changed much less than the 1-*σ* confidence range).

## 4. DISCUSSION

The geometric albedos of asteroids have been derived with several techniques, including (1) infrared observations of thermal radiation, (2) time-delay and Doppler measurements of radar signals emitted from ground-based stations, (3) measurements of the time duration during occultation of background stars, (4) direct imaging with adaptive optics techniques, (5) *in situ* measurements of the bidirectional reflectance using spacecraft cameras, and (6) an empirical method using polarimetric slopes. Since techniques (1)−(4) are contrived to derive the sizes (and shapes) rather than the albedo values of asteroids, it is essential to combine the visible magnitudes to derive the albedos using these four techniques. It is, however, true that optical observations have not been conducted strenuously despite recent comprehensive surveys, especially by means of (1) and (2). In general, the visible magnitudes have been obtained using uncalibrated observation data on the presumption that asteroids have a phase function (typically *G* = 0.15 for an IAU *H*, *G* magnitude system). In fact, it is quite diffcult to make optical observations of asteroids at the opposition point (or at least within a hypothetical circle having the equal apparent diameter of the Sun) because the orbital planes of almost all asteroids do not coincide with Earth's orbital plane (i.e., *i* ≠ 0). Fast time variations of the cross-sectional areas of asteroids also make it diffcult to derive the albedo values even when the phase function and dimension of the asteroids are known. This effect is especially significant for elongated asteroids. In contrast, direct observations of asteroidal surfaces via space missions provide more reliable information about the albedos, eliminating the above uncertainties regarding the cross-sectional areas. However, note that the calibration processes of onboard data have some uncertainties, such as the standardization of magnitude systems and unavailability of flat-fielding data in space. Some *in situ* data were acquired at large phase angles (e.g., (433) Eros by the NEAR-Shoemaker mission). Thus, it is not straightforward to derive the accurate geometric albedo values.

### 4.1. Sources of Uncertainties

The parameters we obatained in this study, especially the albedos, may have suffered from errors that are difficult to quantify correctly. For one of our observations (2018-04-07 CTIO with *α* = 2.79°), the magnitude is ~ 1-*σ* fainter than the model prediction (Figure 2). Although a 1-*σ* of deviation is not necessarily significant, there is a possibility that the value actually was affected by an unknown systematic shift, such as those caused by imperfect phase function modeling (Equation 3) and weather effects (Table 2). Although it was almost invisible, some artifact as in 2018-04-11 SAAO (see the last part of Section 2) may have affected this image in a non-trivial way. This data, however, does not seem to critically affect the quantities derived in this work, as shown in Table 3.

In addition, as the radiance factor calculation requires the absolute physical value of the projected area (Equation 5), the uncertainty in the physical size affects the *I*/*F* value. From Equation (5), *I*/*F* ∝ 1/*S*_{proj} (~ *D*^{–2} to first order), if all other quantities including the observed values are kept constant. The radar models we used give the maximum sizes along each principal axes as (4.60, 2.28, 1.92) ± 0.10km (Hudson et al. 2003), which is slightly smaller but within the uncertainty ranges obtained by Chang'e-2, (4.75 ± 0.48, 1.95 ± 0.19)km (Huang et al. 2013). Ignoring the uncertainty ranges, *S*_{proj} may have been underestimated systematically by ~ 5% in our estimation. Hence, our albedo estimations, *p*_{V} ≡ *I*/*F*(*α* = 0) and *A*_{5} ≡ *I*/*F*(*α* = 5°), could have been overestimated by 0.009 and 0.005, respectively. These uncertainties are negligible compared to the statistical uncertainties (only ~ 10–25% of the error-bars; see Table 3).

### 4.2. Comparison with Flyby Observations

The Chinese Chang'e-2 spacecraft succesfully obtained optical images of Toutatis during its December 2012 flyby at a distance of ~ 1.3km, and the hemispherical albedos were derived (Zhao et al. 2016). These results allow testing of the albedo estimates gained from different approaches including this work. It is, however, diffcult to directly compare their results to ours since they estimated the hemispherical albedo (assumed to be equal to the Bond albedo, *A*_{B}), different from bidirectional albedos derived in this work (e.g., *p*_{V} and *A*_{5}).

Using the IAU H, G system (Bowell et al. 1989), the Bond albedo can be obtained from *p*_{V} using a numerical approximation of the phase integral:

$$${A}_{\mathrm{B},\mathrm{V}}=\left(0.286+0.656G\right){p}_{\mathrm{V}},$$$ | (10) |

where *A*_{B,V} is the Bond albedo in V-band.^{6} Using the uncertainty in *p*_{V} and *G* = 0.0–0.2, we obtain extrema of *A*_{B,V} = 0.04–0.10. This is smaller than the values derived in Zhao et al. (2016): 0.2083, 0.1269, and 0.1346 in the R, G, and B bands of Chang'e-2, respectively. Using Equation (10), these can be transformed into the geometric albedos of 0.50–0.73, 0.30–0.44, and 0.32–0.47 assuming *G* = 0.0–0.2 and identical phase functions, respectively.

The discrepancy may have been caused by the assumptions used to compare the two physically different albedos (bidirectional and Bond). To measure the Bond albedo directly, one needs to observe the object at all the possible 4*π* sr directions, which is impractical. Thus, we assumed a phase function (Equation 10) to convert the *p*_{V} to *A*_{B,V}, while Zhao et al. (2016) assumed illumination conditions (incidence angles) to convert the radiance (Wm^{–2} sr^{–1} *μ*m^{–1}) into the albedo, for example. The phase angle coverages of the observations in both works are also different.

Therefore, this discrepancy is a result of imperfect assumptions used to indirectly infer the Bond albedo in the two works, while these systematic uncertainties are not taken into account in the error bars. We emphasize that our *p*_{V} and *A*_{5} should be accurate within the derived uncertainty range, because these were obtained from observations made at *α* = 0.6°–2.8°, close to where these quantities are defined.

### 4.3. Comparison with Thermal Infrared Studies

The sizes and albedos provided by infrared catalogs such as IRAS (Tedesco et al. 2002), AKARI (Usui et al. 2011), and NEOWISE (Masiero et al. 2017) provide an invaluable overview on asteroids of different taxonomic types distributed across the solar system; nonetheless, it is also important to note that some asteroids with a small number of detections and elongated shapes have large uncertainties due to rotation.

Masiero et al. (2017)^{7} obtained log_{10} *p*_{V} = –0.392± 0.166 with a diameter *D* = 1.788±0.376km for Toutatis from 6 NEOWISE infrared observations in both W1 and W2 bands at phase angle *α* ~ 48°, in combination with *H* = 15.30, *G* = 0.10, and a beaming parameter *η* = 1.007 ± 0.223. The *p*_{V} value can be transformed to *p*_{V} = $$ {0.406}_{-0.129}^{+0.189}$$, if we assume the posterior distribution of their log_{10} *p*_{V} estimation is truly, or at least similar to, Gaussian. It should be noted that the observations were made at a single phase angle, so the Near Earth Asteroid Thermal Model (NEATM) used in Masiero et al. (2017) may not give an accurate estimation, but other advanced modelings may be more adequate, such as the advanced thermophysical model (ATPM; Rozitis & Green 2011, which is used for space mission related targets, e.g., Yu et al. 2014, 2015).

Meanwhile, Usui et al. (2014) noted a 22% deviation for the 1-*σ* confidence range among three radiometric survey catalogs by IRAS, AKARI, and WISE (Tedesco et al. 2002; Usui et al. 2011; Mainzer et al. 2011). The deviation is caused not only by the different types of thermal models but also by the observed rotational phases of the asteroids. It is worth noting that the size and albedo determination can be more uncertain for individual targets than the estimated 22% from the *average* trend.

Note also the currently ongoing debate about the accuracy of the NEOWISE results (see, e.g., Myhrvold 2018a; Wright et al. 2018; Myhrvold 2018c). According to Wright et al. (2018), however, the results we used for Toutatis (Masiero et al. 2017) are not affected by any problem. Meanwhile, Myhrvold (2018b) discusses a new approach to the NEATM by introducing a new parametrization and a physically valid yet simple generalization for the case that the reflected sunlight is not negligible in the thermal spectra. This new approach may improve the NEOWISE results with respect to accuracy and consistency.

In case of Toutatis, one needs to pay close attention to the rotational phase because of the very elongated shape of the asteroid. In fact, the size given in the NEOWISE catalog is only 73 % of the mean diameter derived from radar observations (i.e., *D* = 2.45km). Using the relationship *D* = 1329 [km]*p*_{V}^{–0.5} 10^{–HV/5} (derived in, e.g., the appendix of Pravec & Harris 2007), the *p*_{V} value appears to be overestimated by a factor of ~ 1=0.73^{2} = 1.88. Despite this not being an accurate statement since the change in *p*_{V} changes the shape of the thermal spectrum of the asteroid, we may crudely estimate the updated albedo to be *p*_{V} ~ 0.406=1.88 = 0.22. This coincides with the albedo expected for an S-type asteroid.

### 4.4. Comparison with Polarimetric Studies

Lupishko et al. (1995) studied the polarimetric properties of Toutatis and derived *p*_{V} = 0.13 using the empirical polarimetric slope−albedo law stated in Zellner et al. (1977). Subsequently, Mukai et al. (1997) derived the same value, *p*_{V} = 0.13, using the polarimetric slope−albedo law in Dollfus & Zellner (1979). The empirical relationship, first noted in Widorn (1967), is

$$${\mathrm{log}}_{10}{p}_{\mathrm{V}}={C}_{1}{\mathrm{log}}_{10}h+{C}_{2},$$$ | (11) |

where *h* denotes the polarimetric slope around the inversion angle in %/° and where *C*_{1} and *C*_{2} are constants. Note that the albedo in Zellner et al. (1977) and Dollfus & Zellner (1979) is different from the one that we currently apply. Thus, these authors derived a set of constants (*C*_{1}, *C*_{2}) = (–0.93, –1.78) based on the laboratory polarimetric measurements, where the measurements defined the albedo at the normal (0°) incidence and 5° angle of emergence (i.e., *α* = 5°, not 0°). It is challenging to measure the geometric albedo, which is defined at *α* = 0°, because the light source is in the line of sight of the detector, thus hiding the sample. Although we acknowledge these early efforts, the polarimetric albedo in Lupishko et al. (1995) and Mukai et al. (1997) needs to be updated to match the definition of the geometric albedo.

Multiple studies made efforts to derive the set of constants for obtaining the geometric albedo from the polarimetric slope (Masiero et al. 2012; Cellino et al. 2015; Lupishko 2018), involving better ways to calibrate the data from different sources. We estimated the geometric albedos using the polarimetric slope−albedo law using these sets of parameters, and the results are summarized in Table 4. To first order, the Gaussian 1-*σ* confidence range of *p*_{V} is estimated from Equation (11) and

$$$\u2206{p}_{V}~\sqrt{{\left(\frac{\partial {p}_{V}}{\partial {C}_{1}}\right)}^{2}\u2206{C}_{1}^{2}+{\left(\frac{\partial {p}_{V}}{\partial {C}_{2}}\right)}^{2}\u2206{C}_{2}^{2}+{\left(\frac{\partial {p}_{V}}{\partial h}\right)}^{2}\u2206{h}^{2}},$$$ | (12) |

where ∆*X* is the Gaussian 1-*σ* confidence range for the parameter *X*. We summarize the geometric albedo derived using different observations and different sets of parameters in Table 4. These polarimetric results are consistent with the geometric albedo that we derived through the direct photometric observation at opposition.

Finally, we derived the grain size on the surface. The empirical relationships regarding the grain sizes (Geake & Dollfus 1986) include *A*_{5}, not *p*_{V}. In Section 3.2, we derived *A*_{5} = $$ {0.104}_{-0.034}^{+0.050}$$ in the V-band. From previous polarimetric data, it has been found that *A*_{5} = 0.13 ± 0.02 via the slope−albedo law (Lupishko et al. 1995; Mukai et al. 1997). By comparing with terrestrial rock powders (Figure 4), the grain size is estimated as *d* ≲ 50 *μ*m.

The grain size is smaller than that of (3200) Phaethon (*d* ≳ 360 *μ*m, Ito et al. 2018) and (1566) Icarus (*d* ~ 100 – 130 *μ*m, Ishiguro et al. 2017) that were calculated by using the same relationship. In addition, we compare the *P*_{max}−*A*_{5} relationship of these three asteroids and laboratory sample data (Figure 4). It is likely that (3200) Phaethon and (1566) Icarus are covered with grains larger than that of Toutatis, although other factors such as the porosity of grains can affect the size estimate. These results may suggest that the near-Sun environment at less than approximately 0.2 au affects the paucity of small grains on the surface.

Note that the observations or calculations of *A*_{5} and *P*_{max} for these asteroids were sometimes made at wavelengths different from those used in the laboratory experiments in Geake & Dollfus (1986). Nevertheless, the deviations that can be expected are within the measurement uncertainties.

## 5. CONCLUSIONS

In this work, we analyzed imaging data of the S_{k}-type asteroid (4179) Toutatis at near-opposition configuration obtained by KMTNet. We obtained a robust geometric albedo value of *p*_{V} = $$ {0.185}_{-0.039}^{+0.045}$$ from the observed radiance factors and the shape model. Our value is typical for an S-complex asteroid, and thus shows the future potential of using KMTNet for the opposition monitoring of asteroids. Other measures of albedo using the radiance factors (*f*_{1} and *f*_{2}), which were introduced in previous works, agree with our results within the uncertainties albeit it is substantial. Thermal modeling and polarimetry give consistent results if we use appropriate diameter values or empirical relationships, respectively. The albedo at *α* = 5° is also determined from the radiance factor and is consistent with Toutatis being covered with particles smaller than the asteroids in the near-Sun environment found from empirical relationships involving the grain size.

## Acknowledgments

We thank the two anonymous reviewers for their valuable comments which led to important improvements of the manuscript. This research has made use of the KMTNet system operated by the Korea Astronomy and Space Science Institute (KASI) and the data were obtained at three host sites of CTIO in Chile, SAAO in South Africa, and SSO in Australia. The observations around the opposition were made using telescopes in the frame of contract research between Seoul National University (SNU) and the Korea Astronomy and Space Science Institute (KASI) for the DEEP−South project. This research was supported by the Korea Astronomy and Space Science Institute under the R&D program supervised by the Ministry of Science, ICT and Future Planning.

^{1}https://github.com/astropy/astroscrappy version 1.0.5 using a separable median filter.

^{4}Using a Lambertian

*plate*rather than a Lambertian

*sphere*has historical rather than mathematical or physical reasons (see, e.g., Russell 1916). If we used a Lambertian

*sphere*, the

*F*value would be 1/2

*π*rather than 1/

*π*.

^{5}The angle of

*α*= 5° is chosen by experimentalists, since the albedo measurements were found to give the most reproducible results with

*α*= 5° (see Geake & Dollfus 1986, and references therein). A smoked MgO screen illuminated at an incidence angle of 0° is used as an analog of Lambertian plate. The

*A*

_{5}is the ratio of the fluxux of the target of interest and the MgO screen, while both fluxuxes are measured at

*α*= 5° with zero incidence angle.

^{6}Conventionally,

*A*

_{B,V}= (0.290 + 0.684G)

*p*

_{V}suggested by Bowell et al. (1989) is widely being used. However, Myhrvold (2016) found that the best fitting function is

*A*

_{B,V}= (0.286 + 0.656G)

*p*

_{V}. We confirmed the latter result and used it here.

^{7}The machine-readable table is available at http://iopscience.iop.org/1538-3881/154/4/168/suppdata/ajaa89ect1_mrt.txt

## References

- Astropy Collaboration: Price-Whelan, A. M., Sipőcz, B. M., et al. , (2018), The Astropy Project: Building an Open-science Project and Status of the v2.0 Core Package, AJ, 156, p123.
- Astropy Collaboration: Robitaille, T. P., Tollerud, E. J., et al. , (2013), Astropy: A community Python Package for Astronomy, A&A, 558, pA33.
- Belskaya, I. N., & Shevchenko, V. G., (2000), Opposition Effect of Asteroids, Icarus, 147, p94. [https://doi.org/10.1006/icar.2000.6410]
- Bowell, E., Hapke, B., Domingue, D., et al. , (1989), Application of Photometric Models to Asteroids, in Asteroids II, Binzel, R. P., Gehrels, T., & Matthews, M. S., ed.524.
- Bradley, L., Sipocz, B., Robitaille, T., et al. , (2017), astropy/photutils: v0.4.
- Bus, S. J., & Binzel, R. P., (2002), Phase II of the Small Main-Belt Asteroid Spectroscopic Survey. A Feature-Based Taxonomy, Icarus, 158, p146. [https://doi.org/10.1006/icar.2002.6856]
- Cellino, A., Bagnulo, S., Gil-Hutton, R., et al. , (2015), On the Calibration of the Relation between Geometric Albedo and Polarimetric Properties for the Asteroids, MNRAS, 451, p3473. [https://doi.org/10.1093/mnras/stv1188]
- Craig, M., Crawford, S., Seifert, M., et al. , (2017), astropy/ccdproc: v1.3.0.post1.
- DeMeo, F. E., Binzel, R. P., Slivan, S. M., et al. , (2009), An Extension of the Bus Asteroid Taxonomy into the Near-Infrared, Icarus, 202, p160. [https://doi.org/10.1016/j.icarus.2009.02.005]
- Dollfus, A., (1998), Lunar Surface Imaging Polarimetry: I. Roughness and Grain Size, Icarus, 136, p69. [https://doi.org/10.1006/icar.1998.6008]
- Dollfus, A., & Bowell, E., (1971), Polarimetric Properties of the Lunar Surface and its Interpretation Part I. Telescopic Observations, A&A, 10, p29.
- Dollfus, A., & Zellner, B., (1979), Optical Polarimetry of Asteroids and Laboratory Samples, T. Gehrels, & M. S. Matthews, ed.p170.
- Flewelling, H. A., Magnier, E. A., Chambers, K. C., et al. , (2016), The Pan-STARRS1 Database and Data Products, arXiv:1612.05243.
- Geake, J. E., & Dollfus, A., (1986), Planetary Surface Texture and Albedo from Parameter Plots of Optical Polarization Data, MNRAS, 218, p75. [https://doi.org/10.1093/mnras/218.1.75]
- Ginsburg, A., Sipocz, B., Parikh, M., et al. , (2018), astropy/astroquery: v0.3.8 release.
- Harris, A. W., (1998), A Thermal Model for Near-Earth Asteroids, Icarus, 131, p291. [https://doi.org/10.1006/icar.1997.5865]
- Howell, E. S., Britt, D. T., Bell, J. F., et al. , (1994), Visible and Near-infrared Spectral Observations of 4179 Toutatis, Icarus, 111, p468. [https://doi.org/10.1006/icar.1994.1157]
- Hu, S., Ji, J., Richardson, D. C., et al. , (2018), The Formation Mechanism of 4179 Toutatis' Elongated Bilobed Structure in a Close Earth Encounter Scenario, MNRAS, 478, p501. [https://doi.org/10.1093/mnras/sty1073]
- Huang, J., Ji, J., Ye, P., et al. , (2013), The Ginger-shaped Asteroid 4179 Toutatis: New Observations from a Successful Flyby of Chang'e-2, Sci. Rep., 3, p3411. [https://doi.org/10.1038/srep03411]
- Hudson, R. S., & Ostro, S. J., (1995), Shape and Non-Principal Axis Spin State of Asteroid 4179 Toutatis, Science, 270, p84. [https://doi.org/10.1126/science.270.5233.84]
- Hudson, R. S., Ostro, S. J., & Scheeres, D. J., (2003), High-resolution Model of Asteroid 4179 Toutatis, Icarus, 161, p346. [https://doi.org/10.1016/s0019-1035(02)00042-8]
- Hudson, R. S., Ostro, S. J., & Scheeres, D. J., (2004), Asteroid Radar Shape Models, 4179 Toutatis, HIRES, NASA Planetary Data System, 16, EAR.
- Ishiguro, M., Kuroda, D., Watanabe, M., et al. , (2017), Polarimetric Study of Near-Earth Asteroid (1566) Icarus, AJ, 154, p180. [https://doi.org/10.3847/1538-3881/aa8b1a]
- Ishiguro, M., Nakayama, H., Kogachi, M., et al. , (1997), Maximum Visible Polarization of 4179 Toutatis in the Apparition of 1996, PASJ, 49, pL31. [https://doi.org/10.1093/pasj/49.5.l31]
- Ito, T., Ishiguro, M., Arai, T., et al. , (2018), Extremely Strong Polarization of an Active Asteroid (3200) Phaethon, Nat. Comm., 9, p2486. [https://doi.org/10.1038/s41467-018-04727-2]
- Ji, J., Jiang, Y., Zhao, Y., et al. , (2016), Chang'e-2 Spacecraft Observations of Asteroid 4179 Toutatis, IAUS, 318, p144.
- Jiang, Y., Ji, J., Huang, J., et al. , (2015), Boulders on Asteroid Toutatis as Observed by Chang'e-2, NatSR, 5, 16029. [https://doi.org/10.1038/srep16029]
- Kim, S.-L., Lee, C.-U., Park, B.-G., et al. , (2016), KMTNET: A Network of 1.6 m Wide-Field Optical Telescopes Installed at Three Southern Observatories, JKAS, 49, p37. [https://doi.org/10.5303/jkas.2016.49.1.37]
- Lazzarin, M., Barbieri, C., & Marzari, F., (1994), CCD Reflectance Spectra of Apollo Asteroid 4179 Toutatis, Planet. Space Sci., 42, p327. [https://doi.org/10.1016/0032-0633(94)90105-8]
- Lee, M., & Ishiguro, M., (2018), Opposition Effect on S-type Asteroid (25143) Itokawa, A&A, 616, A178. [https://doi.org/10.1051/0004-6361/201832721]
- Lupishko, D. F., (2018), Generalized Calibration of the Polarimetric Albedo Scale of Asteroids, Sol. Syst. Res., 52, p98. [https://doi.org/10.1134/s0038094618010069]
- Lupishko, D. F., Vasilyev, S. V., Efimov, J. S., et al. , (1995), UBVRI-polarimetry of Asteroid 4179 Toutatis, Icarus, 113, p200. [https://doi.org/10.1006/icar.1995.1016]
- Mainzer, A., Bauer, J., Grav, T., et al. , (2011), Preliminary Results from NEOWISE: An Enhancement to the Wide-field Infrared Survey Explorer for Solar System Science, ApJ, 731, p53. [https://doi.org/10.1088/0004-637x/731/1/53]
- Mann, A. W., & von Braun, K., (2015), Revised Filter Profiles and Zero Points for Broadband Photometry, PASP, 127, p102. [https://doi.org/10.1086/680012]
- Markov, A., & Barabashev, N. P., (1926), On the Reflection of Light from the Lunar Surface, Astron. Zh., 3, p55.
- Masiero, J. R., Mainzer, A. K., Grav, T., et al. , (2012), A Revised Asteroid Polarization-Albedo Relationship Using WISE/NEOWISE Data, ApJ, 749, p104. [https://doi.org/10.1088/0004-637x/749/2/104]
- Masiero, J. R., Nugent, C., Mainzer, A. K., et al. , (2017), NEOWISE Reactivation Mission Year Three: Asteroid Diameters and Albedos, AJ, 154, p168. [https://doi.org/10.3847/1538-3881/aa89ec]
- Moon, H.-K., Kim, M.-J., Yim, H.-S., et al. , (2016), in IAU Symposium, Vol. 318, Asteroids: New Observations, New Models, ed.Chesley, S. R., Morbidelli, A., Jedicke, R., & Farnocchia, D., p306-310.
- Mueller, B. E. A., Samarasinha, N. H., & Belton, M. J. S., (2002), The Diagnosis of Complex Rotation in the Lightcurve of 4179 Toutatis and Potential Applications to Other Asteroids and Bare Cometary Nuclei, Icarus, 158, p305. [https://doi.org/10.1006/icar.2002.6892]
- Mukai, T., Iwata, T., Kikuchi, S., et al. , (1997), Polarimetric Observations of 4179 Toutatis in 1992/1993, Icarus, 127, p452. [https://doi.org/10.1006/icar.1997.5700]
- Myhrvold, N., (2016), Comparing NEO Search Telescopes, PASP, 128, p045004. [https://doi.org/10.1088/1538-3873/128/962/045004]
- Myhrvold, N., (2018), An Empirical Examination of WISE/NEOWISE Asteroid Analysis and Results, Icarus, 314, p64. [https://doi.org/10.1016/j.icarus.2018.05.004]
- Myhrvold, N., (2018), Asteroid Thermal Modeling in the Presence of Reflected Sunlight, Icarus, 303, p91. [https://doi.org/10.1016/j.icarus.2017.12.024]
- Myhrvold, N. P., (2018), Response to Wright et al. 2018: Even More Serious Problems with NEOWISE, arXiv:1812.06516.
- Nakamura, T., & Fuse, T., (1998), Asteroid Lightcurve Observations near Stationary Points: (4179) Toutatis in 1997 February, PASJ, 50, p531. [https://doi.org/10.1093/pasj/50.5.531]
- Ostro, S. J., (1993), Planetary Radar Astronomy, Rev. Mod. Phys., 65, p1235. [https://doi.org/10.1103/revmodphys.65.1235]
- Ostro, S. J., Hudson, R. S., Jurgens, R. F., et al. , (1995), Radar Images of Asteroid 4179 Toutatis, Science, 270, p80. [https://doi.org/10.1126/science.270.5233.80]
- Pravec, P., & Harris, A. W., (2007), Binary Asteroid Population. 1. Angular Momentum Content, Icarus, 190, p250.
- Reddy, V., Sanchez, J. A., Gaffey, M. J., et al. , (2012), Composition of Near-Earth Asteroid (4179) Toutatis, Icarus, 221, p1177. [https://doi.org/10.1016/j.icarus.2012.10.005]
- Robitaille, T., (2018), reproject: Astronomical Image Reprojection in Python.
- Rozitis, B., & Green, S. F., (2011), Directional Characteristics of Thermal-infrared Beaming from Atmosphereless Planetary Surfaces - a New Thermophysical Model, MNRAS, 415, p2042. [https://doi.org/10.1111/j.1365-2966.2011.18718.x]
- Russell, H. N., (1916), On the Albedo of the Planets and Their Satellites, ApJ, 43, p173. [https://doi.org/10.1086/142244]
- Scheeres, D. J., Ostro, S. J., Hudson, R. S., et al. , (1998), Dynamics of Orbits Close to Asteroid 4179 Toutatis, Icarus, 132, p53. [https://doi.org/10.1006/icar.1997.5870]
- Shevchenko, V. G., & Tedesco, E. F., (2006), Asteroid Albedos Deduced from Stellar Occultations, Icarus, 184, p211. [https://doi.org/10.1016/j.icarus.2006.04.006]
- Shkuratov, I. G., & Opanasenko, N. V., (1992), Polarimetric and Photometric Properties of the Moon: Telescope Observation and Laboratory Simulation. II - The Positive Polarization, Icarus, 99, p468.
- Spencer, J. R., Akimov, L. A., Angeli, C., et al. , (1995), The Lightcurve of 4179 Toutatis: Evidence for Complex Rotation., Icarus, 117, p71. [https://doi.org/10.1006/icar.1995.1143]
- Stetson, P. B., (1987), DAOPHOT - A Computer Program for Crowded-field Stellar Photometry, PASP, 99, p191. [https://doi.org/10.1086/131977]
- Takahashi, Y., Busch, M. W., & Scheeres, D. J., (2013), Spin State and Moment of Inertia Characterization of 4179 Toutatis, AJ, 146, p95. [https://doi.org/10.1088/0004-6256/146/4/95]
- Tatsumi, E., Domingue, D., Hirata, N., et al. , (2018), Vis-NIR Disk-integrated Photometry of Asteroid 25143 Itokawa around Opposition by AMICA/Hayabusa, Icarus, 311, p175. [https://doi.org/10.1016/j.icarus.2018.04.001]
- Tedesco, E. F., Noah, P. V., Noah, M., et al. , (2002), The Supplemental IRAS Minor Planet Survey, AJ, 123, p1056.
- Usui, F., Hasegawa, S., Ishiguro, M., et al. , (2014), A Comparative Study of Infrared Asteroid Surveys: IRAS, AKARI, and WISE, PASJ, 66, p56. [https://doi.org/10.1093/pasj/psu037]
- Usui, F., Kuroda, D., Muller, T. G., et al. , (2011), Asteroid Catalog Using Akari: AKARI/IRC Mid-Infrared Asteroid Survey, PASJ, 63, p1117. [https://doi.org/10.1093/pasj/63.5.1117]
- van Dokkum, P. G., (2001), Cosmic-Ray Rejection by Laplacian Edge Detection, PASP, 113, p1420.
- Whipple, A. L., & Shelus, P. J., (1993), Long-term Dynamical Evolution of the Minor Planet (4179) Toutatis, Icarus, 105, p408. [https://doi.org/10.1006/icar.1993.1137]
- Widorn, T., (1967), Zur photometrischen Bestimmung der Durchmesser der Kleinen Planeten, Annalen der Universitäts-Sternwarte Wien, 27, p109.
- Wright, E., Mainzer, A., Masiero, J., et al. , (2018), Response to "An Empirical Examination of WISE/NEOWISE Asteroid Analysis and Results", arXiv:1811.01454.
- Yu, L. L., Ji, J., & Wang, S., (2014), Shape, Thermal and Surface Properties Determination of a Candidate Spacecraft Target Asteroid (175706) 1996 FG3, MNRAS, 439, p3357. [https://doi.org/10.1093/mnras/stu164]
- Yu, L. L., & Ji, J., (2015), Surface Thermophysical Properties Determination of OSIRIS-REx Target Asteroid (101955) Bennu, MNRAS, 452, p368. [https://doi.org/10.1093/mnras/stv1270]
- Zellner, B., Leake, M., Lebertre, T., et al. , (1977), Lunar and Planetary Science Conference Proceedings, 8, Lunar and Planetary Science Conference Proceedings, ed.Merril, R. B., p1091.
- Zhao, Y., Ji, J., Huang, J., et al. , (2015), Orientation and Rotational Parameters of Asteroid 4179 Toutatis: New Insights from Chang'e-2's Close Flyby, MNRAS, 450, p3620. [https://doi.org/10.1093/mnras/stv792]
- Zhao, D.-F., Liu, P., Zhao, W., et al. , (2016), Reflectance of Asteroid 4179 Toutatis Based on Its Space Optical Image, AcASn, 40, p555.

## Appendix

## APPENDIX A. ROTATIONAL MODEL

We calculated the projected geometrical cross-sectional area, *S*_{proj}(*t*; *φ*_{1}, *φ*_{2}) of the asteroid as seen from Earth at the epoch of the observations (*t*), assuming arbitrary rotational and precessional offsets, *φ*_{1} and *φ*_{2}. In the calculation, we employed the high-resolution shape model^{8} (Hudson et al. 2003, 2004). We considered the long-axis mode rotation and precession along (*λ*, *β*) = (180°, –52°) in the ecliptic coordinate (Hudson & Ostro 1995). Although these rotation and precession periods have been reported to be *P*_{1} = 5.38 days and *P*_{2} = 7.40 days, respectively (Zhao et al. 2015), it is necessary to perform fine tuning by adding the offsets *φ*_{1} and *φ*_{2} to fit the observed magnitudes to the modeled ones because the periods in the model are not determined well enough to be extrapolated to the epoch of our observations (uncertainty in periods multiplied by the time difference is much larger than the periods themselves). We thus defined the rotational angle along the long-axis and the precession axis (*θ*_{1} and *θ*_{2}, respectively) as follows:

$$${\theta}_{j}\left(t,{\phi}_{j}\right)=\frac{2\pi \left(t-{t}_{0}\right)}{{P}_{j}}+{\phi}_{j},$$$ | (13) |

for *j* = 1, 2, where *t*_{0} denotes the time of the reference epoch.

We used IDL POLYSHADE to calculate the total projected area *S*_{proj}(*t*; *φ*_{1}, *φ*_{2}) for a given pair of (*φ*_{1}, *φ*_{2}) with grid size of ∆*φ*_{1} = 5° and ∆*φ*_{2} = 2° and time step size of 1=3 days. The projected area as a function of the direction to the observer is shown in Figure 5. The color and contour show the projected area viewed by an observer at an infinite distance to the direction specified the body-fixed frame, (*θ*_{bf}, *ϕ*_{bf}). As can be seen, an uncertainty of ~ 5° in the viewing direction causes a fractional uncertainty (normally ≪ 10%) which is much smaller than the photometric uncertainty of *H*_{V}(*α*) in Figure 2 and Table 2. Thus, the choices for the grid size of the offsets (∆*φ*_{1} and ∆*φ*_{2}) are reasonable.

Each of the 12,960 models for (*φ*_{1}, *φ*_{2}) pair is then spline interpolated when necessary. As mentioned in Section 3.1, the total reflected fluxux is proportional to the projected area when *α* is small. Hence, for each rotational model, we fit a such that *I*_{0}(*t*) = *aS*_{proj}(*t*), where *I*_{0}(*t*) = 10^{–HV(t)/2.5} and its uncertainty is ∆*I*_{0}(*t*) = ln(10)*I*_{0}(*t*)∆*V* (*t*)/2.5 (using Equation (3) with fixed *G*). Here, *H*_{V}(*t*) is the phase-corrected absolute magnitude, *H*_{V}, calculated using Equation (3) at epoch *t*. We used intensity units rather than magnitudes since uncertainties in magnitude are not Gaussian, precluding the use of *χ*^{2} minimization. The reduced *χ*^{2} statistic is calculated for each case by using

$$${\chi}_{\mathrm{r}\mathrm{e}\mathrm{d}}^{2}\left({\phi}_{1},{\phi}_{2}\right)=\frac{1}{N-1}\sum _{i=1}^{N}{\left(\frac{{I}_{0}\left({t}_{i}\right)-a{S}_{\mathrm{p}\mathrm{r}\mathrm{o}\mathrm{j}}\left({t}_{i};{\phi}_{1},{\phi}_{2}\right)}{\u2206{I}_{0}\left({t}_{i}\right)}\right)}^{2},$$$ | (14) |

where the denominator *N* – 1 is the number of degrees of freedom, *N* = 9 is the number of data points, and *t _{i}* is the

*i*-th observational time (Table 2). By varying

*φ*

_{1}and

*φ*

_{2}, we find an

*a*

_{0}for which

*a*

_{0}= arg min

_{a}$$ {\chi}_{\mathrm{r}\mathrm{e}\mathrm{d}}^{2}$$(

*φ*

_{1},

*φ*

_{2}). The minimum $$ {\chi}_{\mathrm{r}\mathrm{e}\mathrm{d}}^{2}$$ is then saved for each model.

The "good" models are selected such that the reduced *χ*^{2}-statistics is smaller than 1 + $$ {\chi}_{\mathrm{r}\mathrm{e}\mathrm{d},\mathrm{m}\mathrm{i}\mathrm{n}}^{2}$$, where $$ {\chi}_{\mathrm{r}\mathrm{e}\mathrm{d},\mathrm{m}\mathrm{i}\mathrm{n}}^{2}$$ is the minimum of such reduced *χ*^{2} statistics of all the models. Because we have only one free parameter (a), this gives the models with 1-*σ* confidence interval. The models are then converted into magnitude system using –2.5 log_{10}(*aI*_{0}). The minimum *χ*^{2} model is the red solid line and other "good" models are black faint dotted lines in Figure 2. We confirmed that changes of *G* in the range G ∈ [0.10, 0.50] do not significantly distort the results, except for a nearly parallel shift caused by the nearly constant change (since all *α*'s are similar) in all *H*_{V} values (Equation 3).

## APPENDIX B. ON-LINE MATERIAL

This work was carried using mostly open-source software, and we here disclose the codes we used to obtain our results. The original codes are available in a dedicated GitHub directory^{9} which contains five notebooks:

- •
**01_data_reduction:**the data reduction process described in Section 2 and the generation of Table 2. - •
**02_photometry_comparison:**the comparison of our reduced data with Spencer et al. (1995) and internally developed Moving Object Detection Pipeline (MODP) to cross-check our photometry results. - •
**03_radiance_factor:**the calculation of*I*/*F*,*f*_{1},*f*_{2}, and*A*_{5}values from the photometry results and projected area calculation (see Appendix A), and the generation of Figures 2 and 3. - •
**04_cross_section:**the calculation of*S*_{proj}(*θ*_{bf},*ϕ*_{bf}) in Figure 5 using**ParaView**software to check whether the*S*_{proj}(*t*) values in Appendix A (used**IDL POLYSHADE**) is reasonable. - •
**05_polarimetry_plot:**the generation of Figure 4 and a short discussion of the results in Geake & Dollfus (1986). - •
**06_KMTNet_Seeing:**a short description about weather and seeing condition of our observations in Table 2.