Abstract
A stellar occultation of Gaia DR3 2646598228351156352 by the Centaur (2060) Chiron was observed from the South African Astronomical Observatory on 2018 November 28 UT. Here we present a positive detection of material surrounding Chiron from the 74-inch telescope for this event. Additionally, a global atmosphere is ruled out at the tens of microbars level for several possible atmospheric compositions. There are multiple 3σ drops in the 74-inch light curve: three during immersion and two during emersion. Occulting material is located between 242 and 270 km from the center of the nucleus in the sky plane. Assuming the ring-plane orientation proposed for Chiron from the 2011 occultation, the flux drops are located at 352, 344, and 316 km (immersion) and 357 and 364 km (emersion) from the center, with normal optical depths of 0.26, 0.36, and 0.22 (immersion) and 0.26 and 0.18 (emersion) and equivalent widths between 0.7 and 1.3 km. This detection is similar to the previously proposed two-ring system and is located within the error bars of that ring-pole plane; however, the normal optical depths are less than half of the previous values, and three features are detected on immersion. These results suggest that the properties of the surrounding material have evolved between the 2011, 2018, and 2022 observations.
1. Introduction
The Centaur (2060) Chiron was discovered in 1977 (Gehrels et al. 1977; Kowal 1979). Its orbital semimajor axis ranges from approximately 8.5 to 19 au, with an eccentricity of 0.38 and an inclination of 7°. Chiron's rotational period has been reported as 5.92 hr (low error; Marcialis & Buratti 1993), with a slightly shorter period of approximately 5.5 hr reported from more recent data (Fornasier et al. 2013; Galiazzo et al. 2016). Chiron has exhibited outbursting behavior and is codesignated as comet 95P/Chiron (e.g., Hartmann et al. 1990; Luu & Jewitt 1990; Meech & Belton 1990; Bauer et al. 2004; Jewitt 2009). Uncorrelated with heliocentric distance, Chiron's brightness varies over both hourly and decadal timescales (e.g., Tholen et al. 1988; Belskaya et al. 2010). For rotational variations, the typical peak-to-valley magnitude changes have not exceeded 0.09 mag (Belskaya et al. 2010). Chiron is one of the largest Centaurs; the Herschel Space Telescope determined a diameter of 218 ± 20 km (Fornasier et al. 2013), while data from the Atacama Large Millimeter/submillimeter Array returned a spherical equivalent diameter of km and a preferred elliptical size with semiaxes of 114 × 98 × 62 km (Lellouch et al. 2017).
Stellar occultations provide a method that can directly and accurately determine the sizes and shapes of small bodies in the outer solar system, as well as detecting and measuring atmospheres or rings (e.g., Elliot et al. 2010; Sicardy 2022). Successful observations of stellar occultations by Chiron in 1993 and 1994 were consistent with a nucleus size of >180 km in diameter (Elliot et al. 1995; Bus et al. 1996). Occultation data from 2018 and 2019 were more recently used to derive a Jacobi ellipsoid assuming a hydrostatic equilibrium shape with semiaxes a = 126 ± 22, b = 109 ± 19, and c = 68 ± 12 km, implying a volume equivalent radius of Rvol = 98 ± 17 km (Braga-Ribas et al. 2023).
Notably, based on stellar occultation data from 2011, Chiron is one of only a few small bodies in the solar system known to possibly host a ring system (Ortiz et al. 2015; Sickafoose et al. 2020). A multichord stellar occultation in 2013 by the Centaur (10199) Chariklo revealed two thin rings (Braga-Ribas et al. 2014) and sparked a new field of study. Chariklo is the largest Centaur at ∼250 km in diameter, and the rings are roughly 385–390 and 400–405 km from the center of the nucleus, with widths of roughly 7 and 3 km and normal optical depths of 0.31–0.46 and 0.04–0.14, respectively (e.g., Braga-Ribas et al. 2014; Morgado et al. 2021). A single 70-km wide ring was later detected around the trans-Neptunian object (TNO) Haumea (Ortiz et al. 2017), and distant, inhomogeneous material was recently reported around the TNO (50000) Quaoar (Morgado et al. 2023; Pereira et al. 2023). Haumea and Quaoar are an order of magnitude larger in size and located farther away from the Sun than Chariklo and Chiron. Nonetheless, these bodies share the traits of being vastly different in scale and composition when compared to the giants previously known as the sole ring-bearing planets in our solar system. Study of each of these small bodies is thus merited.
In addition to providing an upper limit on Chiron's size, the stellar occultations in the 1990s detected what seemed to be jetlike features and an asymmetric dust coma around Chiron's nucleus (Elliot et al. 1995; Bus et al. 1996). The 2011 stellar occultation showed discrete, symmetric features on either side of the nucleus, which were initially interpreted as a near-circular arc or dust shell based on the limited, two-chord data set (Ruprecht et al. 2015). Ortiz et al. (2015) combined the 2011 occultation results with rotational light curves and long-term photometric and spectroscopic variations to propose that Chiron has a two-ring system similar to Chariklo, with a mean radius of 324 ± 10 km and a preferred pole orientation in ecliptic coordinates of λ = 144° ± 10°, β = 24° ± 10°. Sickafoose et al. (2020) further analyzed the 2011 occultation data, characterizing the proposed ring material to be at 300 and 309 km from the nucleus center in the ring plane and 2.5–4.5 km in width, given an assumed central chord. The ring widths varied azimuthally by 1.5 (inner ring) and 0.5 (outer ring) km, and the normal optical depths ranged between 0.6 and 0.85 (inner) and 0.5 and 0.71 (outer; Sickafoose et al. 2020). Recently, Braga-Ribas et al. (2023) reported stellar occultations from 2018 and 2019 that did not show any features other than the nucleus: they concluded that a permanent ring similar in optical depth and extension to that at Chariklo was ruled out.
Here we report on data from the same 2018 stellar occultation by Chiron and at the same successful site of the South African Astronomical Observatory (SAAO) as discussed in Braga-Ribas et al. (2023) but from a larger telescope. Surrounding material is detected in these data. The observations and data analyses are described in Sections 2 and 3. Results are presented in Section 4, and a discussion is provided in Section 5.
2. Observations
Characteristics of the 2018 November 28 UT Chiron occultation star are listed in Table 1. The prediction was based on the Gaia DR2 catalog, which was the most current star information at the time (Gaia Collaboration et al. 2016), and Chiron's position was from the JPL Horizons Ephemeris System JPL#121 with offsets applied from our Ephemeris Correction Model (Elliot et al. 2009; Zuluaga et al. 2015; Person et al. 2021). The predicted shadow path is shown in Figure 1. This occultation was also predicted by the Lucky-Star project, using their Numerical Integration of the Motion of an Asteroid (NIMA), as described in Desmars et al. (2015) and Braga-Ribas et al. (2023).
Figure 1. The predicted shadow path for the 2018 November 28 UT Chiron occultation. The solid black lines indicate the northern, center, and southern extent of the shadow path with 3σ error bars as dashed lines. Chiron's shadow was assumed to be 108 km in radius for this prediction.
Download figure:
Standard image High-resolution imageTable 1. Occultation Star Information
Parameter | Value |
---|---|
Designation (GDR3) a | 2646598228351156352 |
α (ICRS, epoch date; hh:mm:ss.ss) a | 23:46:04.34 |
δ (ICRS, epoch date; deg:arcmin:arcsec) a | 02:13:05.5 |
μα (mas yr–1) a | −13.1 ± 0.2 |
μδ (mas yr–1) a | −2.7 ± 0.1 |
Parallax (mas) a | 1.691 ± 0.115 |
Magnitudes a , b | G c = 17.3, GBP = 18.1, GRP = 16.3, J = 15.2, H = 14.8, K = 14.3 |
Predicted midtime (UT) d | 2018 November 28, 20:50:59 ± 00:01:07 |
Predicted closest approach (km) d | 265 ± 395 N |
Relative velocity e | 5.96 km s−1 |
Position angle | 133![]() |
Notes.
a From Gaia Data Release 3 (Gaia Collaboration et al. 2021). b Near-infrared magnitudes from 2MASS (Cutri et al. 2003; Skrutskie et al. 2006). c For comparison, the stars for the 2011 and 2022 Chiron occultations were G = 14.9 and 12.7, respectively. d Based on Gaia Data Release 2 (Gaia Collaboration et al. 2016), with 3σ error bars. e Between the star and Chiron, from the SAAO.Download table as: ASCIITypeset image
This occultation was predicted separately by two different groups. To allow independent campaigns and analyses, observations were planned on two telescopes at the same site: the 40 and 74 inch at the SAAO in Sutherland, South Africa (N latitude −32 22 54, E longitude 20 48 36, altitude 1760 m, WGS84, MPC site code K94). Two identical instruments were used on these telescopes: the Sutherland High-speed Optical Cameras (SHOC; Coppejans et al. 2013). The foundation of each system is an Andor iXon 888 camera system based around an e2v frame-transfer electron-multiplying charge-coupled device. Every occultation observation frame was triggered by a temperature-stabilized and frequency-tuned GPS. The timing accuracy is on the order of microseconds.
Data from the 40-inch telescope taken with a cycle time of 1.5 s were reported in Braga-Ribas et al. (2023). This work is focused on the 74-inch data set, although the 40-inch data are also included here in order to show the light curve and investigate light-curve stability. For the 74-inch observations, the camera was cooled to −50°C, no filters were used, and images were taken using the conventional amplifier, read out at 1 MHz, and digitized with the 16-bit analog-to-digital converter. Data were taken with a cycle time of 0.5 s (including dead time of 3.4 ms) from 20:20:00 to 21:20:00 UT. The field of view for SHOC on the 74-inch telescope is roughly . Seeing varied throughout the observations from 2'' to 3''. The images were binned 16 × 16 pixels for a plate scale of 122.
3. Data Analysis
Raw data were analyzed in order to extract light curves from both the 40 and 74 inch data using circular-aperture photometry. Example images are shown in Figure 2. Four hours of time was requested on each telescope, centered on the occultation midtime; therefore, twilight time was not available to take flat fields. Reduced data using bias subtraction alone or dithered sky flats added noise in terms of increasing the standard deviations of the light curves. Dark images were not required, as the camera cooling has undetectable dark current at these exposure times (Coppejans et al. 2013). For the 40-inch analysis, the comparison star was centroided on each frame, and an offset was applied to determine the location of the target star. Photometry was carried out on both sources. For the 74-inch analysis, the field of view was too small to contain a comparison star of adequate magnitude (see Section 5 for additional discussion). Photometry was carried out on the target star alone, which was centroided on each image. For both data sets, four hand-selected boxes well outside of the apertures and without background stars were averaged to determine the sky value, which was subtracted out. A sequence of aperture sizes were tested, and the highest signal-to-noise ratios (S/Ns) were obtained with apertures of 8 (98) and 6 (8
0) superpixels in diameter for the 74 and 40 inch, respectively. For the 40 inch, the target star signal was divided by the comparison star for differential photometry. The light curves were normalized to 1 by dividing by the mean value of the baseline (flux = 1) and scaled such that the occulted portion had a mean value of zero. Data outside of the occultation were fit by a third-order polynomial to ensure a flat baseline. The resulting light curves are shown in Figure 3. There is an obvious detection of Chiron's nucleus in both data sets, and there are significant dips in the 74-inch light curve that are not present in the 40-inch data.
Figure 2. Example images from the 74- and 40-inch telescope data cubes. (Left) The 74-inch data, taken at 2-Hz cadence. The occultation star is the bright source, and a single, faint comparison (G = 18.7) is detectable just above the noise to the NE. (Right) The 40-inch data, taken at 0.67-Hz cadence. The occultation star and a bright comparison star are labeled.
Download figure:
Standard image High-resolution imageFigure 3. Light curves from the (top) 74-inch and (bottom) 40-inch SAAO telescopes. The portion near the occultation by Chiron's nucleus is shown in more detail in Figures 5 and 6, with significant flux drops highlighted in Figure 7.(The data used to create this figure are available.)
Download figure:
Standard image High-resolution imageThe important size scales for diffraction effects are the Fresnel length, , the stellar diameter, and the integration length of the observations (integration time multiplied by the velocity in the sky plane). λ is the wavelength of the observation, here assumed to be 700 nm, and D is the distance from Earth, 18.355 au at the midtime of the occultation. For this event, the Fresnel length is 0.98 km. The maximum stellar diameter is 0.1 km; the same value is obtained using the TESS stellar diameter (Stassun et al. 2019) projected to Chiron's distance or using NOMAD magnitudes (Zacharias et al. 2004) and the equations for the angular size of a giant star from van Belle (1999). The lengths of the integrations correspond to approximately 3 and 9 km for the 74- and 40-inch data, respectively. The diffraction effects are thus dominated by the integration times, and the stellar profile can be neglected.
We fit the immersion and emersion of the star behind the nucleus following French & Gierasch (1976) for a thin atmosphere with diffraction on a solid surface. This model is based on a radially symmetric atmosphere with a constant scale height much less than the body's radius. The model is applicable when the scale height is much greater than the Fresnel scale, which is the case for this event. Immersion and emersion are fit separately. The fits are carried out at three representative scale heights, 10, 50, and 100 km. The fit parameters are the zero and full light-curve levels, the geometric edge, and the bending parameter, b, which represents the amount of differential bending of light causing a change in flux at the surface. The maximum 3σ bending parameter over all of the tested scale heights is used as the upper limit on the flux change at the surface of Chiron. We then consider the maximum bending parameter along with different constituent molecules to place corresponding upper limits on a possible atmosphere. For the atmospheric calculations, we assume a temperature of 125 K (Campins et al. 1994), as well as a mass of 4.8 ± 2.3 × 1018 kg and volume equivalent radius of 98 ± 17 km from Braga-Ribas et al. (2023). We use this volume equivalent radius in order to have the best estimate of Chiron's gravitational pull. The errors on the mass and size of Chiron's nucleus are carried through when calculating the atmospheric upper limits and their errors.
The closest approach of the chord to the center of the nucleus was determined to be 50 km by assuming a sphere of 210 km in diameter (Lellouch et al. 2017) and using the measured ingress and egress times from the diffraction fitting (results shown in Table 2). This geometry is consistent with that reported in Braga-Ribas et al. (2023, their Figure 9), in which a sphere is shown in the sky-plane figure. Because there is only a single chord for this event, a unique solution for the closest approach by a body of unknown shape cannot be determined. Chiron is likely to be a triaxial ellipsoid, the values for which were constrained with large error bars by Braga-Ribas et al. (2023). Assuming a spherical nucleus is a good first assumption, for example, following Ortiz et al. (2023) and including error bars on the closest approach to allow for a different nucleus shape. See Section 5 for additional discussion. The sky-plane view for the nominal occultation geometry is shown in Figure 4. As noted in Braga-Ribas et al. (2023), the null detection from Boyden Observatory for this event dictates that the chord passed south of the center.
Figure 4. Sky-plane projection of Chiron at the time of the occultation, along with the nominal location of the observed chord and the proposed rings. Chiron's orientation is that of the preferred ring-plane pole from Ortiz et al. (2015), and the chord has a closest approach of 50 km calculated from assuming a spherical nucleus with a radius of 105 km (following Lellouch et al. 2017) and using the ingress and egress times in Table 2. The thick black line, which appears continuous at this resolution, is a plot of individual black dots that each represent one image taken using the 74-inch telescope. Large red dots indicate locations where the light-curve data deviate by 3σ or more below the baseline, with the chord going behind the solid nucleus in this plot. Solid gray ellipses represent the two most prominent dips from Sickafoose et al. (2020), while the dotted gray ellipse represents the ring location and width from Ortiz et al. (2015).
Download figure:
Standard image High-resolution imageTable 2. Results from Diffraction Fitting to the Nucleus
Parameter | Value |
---|---|
Immersion, geometric edge (UT) a | 20: 50: 12.58 ± 0.04 s |
Emersion, geometric edge (UT) a | 20: 50: 43.81 ± 0.02 s |
Occultation duration (s) | 31.23 ± 0.04 |
Chord length (km) | 186.13 ± 0.27 |
Immersion b (10, 50, 100 km; 10−2) | −4.3 ± 5.6, −3.1 ± 4.3, −4.7 ± 5.9 |
Emersion b (10, 50, 100 km; 10−2) | −12.6 ± 6.3, −3.8 ± 5.1, −4.0 ± 6.9 |
Note.
a Following French & Gierasch (1976).Download table as: ASCIITypeset image
To search for symmetric features around Chiron in the light curve, we plot the immersion and emersion portions versus distance from the body center and overlay them in Figure 5. The standard deviations of the baselines are 0.08 and 0.09 normalized flux for the 74- and 40- inch data, respectively. Features that are greater than or equal to 3σ drops from the baseline are flagged. There are three such features during immersion and two during emersion (as highlighted in red in Figure 4). The top and middle panels of Figure 5 respectively demonstrate the data in the sky plane and with the nominal ring-plane geometry. As discussed in Section 5, the bottom panel explores a wider range of possible ring-pole positions and closest-approach values.
Figure 5. Test for symmetric features. The immersion (gray) and emersion (black) portions of the light curves are overlaid to check for features that are symmetric on either side of the nucleus. The overlay is centered on the image closest to the fitted midtime. Dashed lines are provided at the locations of all data points that are 3σ or more below the baseline. (top) Data with respect to distance in the sky plane. (Middle) Data with respect to distance in the nominal ring plane, with the preferred ring-orbital pole from Ortiz et al. (2015). In both panels, data from the 74-inch telescope are shown on the top, and data from the 40-inch telescope are shown on the bottom. No significant features were detected in the 40-inch data. (Bottom) Example alternative geometries for the 74-inch data set; data with respect to the ring plane for the extremes of the ring-orbital pole errors in Ortiz et al. (2015). In order for the nucleus immersion and emersion to align in this panel, the closest-approach (CA) value was adjusted from the nominal value to the values given in the plot.
Download figure:
Standard image High-resolution imageWe fit each of the 3σ deviations in the sky plane using a square-well diffraction model for occultation profiles that treats the blocking material as a uniformly transmitting gray screen with sharp edges, following Elliot et al. (1984). Fit parameters for the model are the apparent optical depth, τ, and the width in the sky plane. Distances from the nucleus and widths in the nominal ring plane are determined using the preferred solution of the proposed Chiron ring system from Ortiz et al. (2015) and our best-fit closest approach of 50 km. We convert to normal optical depths, τN, and equivalent widths by using the opening angle from the preferred ring solution, B = 473.
4. Results
The results from the diffraction fitting of the nucleus are plotted in Figure 6 and listed in Table 2. The nucleus occultation lasted 31.23 s, or slightly over 186 km at Chiron. The maximum 3σ bending parameter is 0.13 on immersion and 0.17 on emersion (the maximum values of b in Table 2 considering +3σ). The considered atmospheric constituents and limits on an atmosphere at immersion and emersion are listed in Table 3. These 3σ upper limits are roughly between 10 and 30 μbar, with constraints on emersion being approximately 20% larger than those on immersion.
Figure 6. Diffraction model fit to the 74-inch data for immersion (left) and emersion (right). Data are shown as black dots. The best-fit models are shown as dashed green, solid red, and dotted blue lines for 10, 50, and 100 km scale heights, respectively. The models are very similar, with a slightly larger diffraction spike on emersion for the 10 km scale height (consistent with the lowest bending parameter, b, for this case in Table 2).
Download figure:
Standard image High-resolution imageTable 3. 3σ Upper Limits on a Possible Chiron Atmosphere
Molecule | Immersion (μbar) | Emersion (μbar) |
---|---|---|
CH4 | 12.7 | 16.5 |
CO | 22.1 | 28.7 |
CO2 | 20.7 | 26.9 |
H2 | 14.3 | 18.6 |
H2O | 23.6 | 30.7 |
N2 | 24.9 | 32.4 |
Download table as: ASCIITypeset image
Figure 7 shows the portions of the light curve containing the 3σ dips, along with plots of the best-fit diffraction models. Table 4 contains the measured and fitted characteristics for each of these significant features. The double dips, two overlapping features on both immersion and emersion, are located between 344 and 365 km from the center of the nucleus in the nominal ring plane. There is an approximately 8km gap between the centers of the two features on each side, although the immersion dips are roughly 13 km closer to the nucleus than the emersion. The immersion dips are wider, at 4.5 and 4.0 km for the exterior and interior dips, respectively, versus 2.2 km ring-plane widths for both on emersion. The normal optical depths are also slightly higher on immersion, at 0.26 and 0.36 versus 0.19 and 0.26. The third, stand-alone dip on immersion is closer to the nucleus, at 317 km from the center. It is the widest feature, at 5.1 km, and has a similar normal optical depth to the other dips of 0.22.
Figure 7. Diffraction model fits to the 3σ drops in the 74 inch light curve. Dots represent data points, and lines represent the best-fit models. The black lines are the diffraction models at high time resolution, and the red lines represent the model fits at the sampling rate of the data set. Overlapping features were fit simultaneously.
Download figure:
Standard image High-resolution imageTable 4. Characteristics of Significant Light-curve Features
Sky-plane | Line-of-sight | Ring-plane | Ring-plane | Normal | Equivalent | ||
---|---|---|---|---|---|---|---|
Midtime | Distance | Width | Optical | Distance | Width | Optical Depth | Width |
(s after 20:20 UT) | (km) a | (km) b | depth τb | (km) c | (km) c | τN d | (km) e |
1783.25 | 269.16 | 5.19 ± 0.26 | 0.36 ± 0.11 | 351.50 | 4.54 ± 0.23 | 0.26 ± 0.08 | 1.16 ± 0.36 |
1784.25 | 263.19 | 4.59 ± 0.30 | 0.46 ± 0.11 | 343.67 | 4.03 ± 0.26 | 0.36 ± 0.08 | 1.24 ± 0.31 |
1787.75 | 242.33 | 5.76 ± 0.57 | 0.30 ± 0.10 | 316.30 | 5.10 ± 0.51 | 0.22 ± 0.07 | 1.11 ± 0.42 |
1871.75 | 258.31 | 3.37 ± 0.75 | 0.35 ± 0.10 | 356.64 | 2.23 ± 0.50 | 0.26 ± 0.07 | 0.76 ± 0.33 |
1872.75 | 264.27 | 3.36 ± 0.78 | 0.25 ± 0.14 | 364.48 | 2.23 ± 0.52 | 0.19 ± 0.07 | 0.67 ± 0.32 |
Notes.
a Distance from the center of the nucleus, as measured in the sky plane. The error is ±half a time step: 0.25 s or 1.49 km. b Based on least-squares fits to the data using a square-well diffraction model for occultation profiles from Elliot et al. (1984). Error bars were determined by fitting the light-curve data at ±1σ flux. c The conversion from the sky plane to the ring plane assumes the preferred pole position from Ortiz et al. (2015) with a closest approach of 50 km. d Normal optical depths assume a ring opening angle of B = 47
Download table as: ASCIITypeset image
5. Discussion
Observations of this 2018 November 28 UT stellar occultation by Chiron were successful from two telescopes at the SAAO. An occultation by the nucleus was observed from each telescope. No other significant flux drops were detected in the 40 inch data set (as reported in Braga-Ribas et al. 2023 and shown in Figure 5), while five significant flux drops are reported here in the 74-inch data set. This event was the first Chiron occultation observation since 2011 in which surrounding material has been detected.
From the 74-inch data, the occultation by the nucleus has a chord length of 186.1 ± 0.3 km. This size is a few kilometers longer than the value obtained in Braga-Ribas et al. (2023); however, (i) we use a diffraction fit in which the geometric edge is located near 0.25 flux in the shadow (derived from the diffraction model for a knife edge), and (ii) based on Table 5 in Braga-Ribas et al. (2023), a different relative velocity was assumed than the 5.96 km s−1 that we have calculated at the SAAO (the egress–ingress times and chord length in their Table 5 imply a velocity of 5.77 km s−1). Our chord length is consistent with previous size estimates for Chiron.
The 74-inch field of view was too small to provide an adequate comparison star for differential photometry. Therefore, care must be taken to prevent misinterpreting variations in flux that are due to changes in sky conditions as opposed to occultations at the targeted body. The detected flux drops are located near the proposed ring locations, and we find no other 3σ drops in the light curve out to more than 1000 km from the nucleus. In addition, the 40-inch data set did have a bright comparison star, and there were no 3σ flux drops in that star's signal within a few thousand kilometers on either side of the occultation midtime (see Figure 8). We thus determine that the star signal was sufficiently stable during the occultation that the detected drops were caused by material near Chiron.
Figure 8. Light curve for the comparison star in the 40-inch data set. The dotted–dashed line indicates the level of 3σ drops from the baseline. The solid line shows the spatial extent of the occultation features discussed in this work.
Download figure:
Standard image High-resolution imageFigure 4 shows the nominal occultation geometry for this event. In the nominal ring plane, the 2018 double-dip features are outside of the previously proposed distances from the nucleus, at 20–40 km exterior to the 324 ± 10 km ring location from Ortiz et al. (2015). However, the geometry varies depending on the exact location of the ring-pole position and the closest-approach distance of the chord. Ortiz et al. (2023) recently presented results from a stellar occultation by Chiron in 2022 December. They found that the structure of the surrounding material evolved since 2011. The 2022 data indicated a broad, tenuous disk of 580 km with ecliptic coordinates of the pole of λ = 151° ± 8° and β = 18° ± 11°, as well as concentrations of material at 325 ± 16 and 423 ± 11 km (Ortiz et al. 2023). This pole position is within the 10° errors of the previous estimate and has similarly sized error bars. To compare occultation geometries between 2011 and 2022, Ortiz et al. (2023) assumed a spherical nucleus of the same size assumed in this work with a radial uncertainty of 10 km to account for any displacement due to a triaxial shape. Sky-plane projections are given in Figure 9 to explore the impact of the different possible ring geometries for the 2018 data set. The left panel of Figure 9 shows the range of ring and ring-pole positions, including the error bars from the preferred solution in Ortiz et al. (2015) along with the 2018 chord closest approach shifted by ±10 km. The right panel in Figure 9 shows the triaxial nucleus shape from Braga-Ribas et al. (2023) with the newest ring-pole position and the locations of ring concentrations from Ortiz et al. (2023). Changing the closest approach (the equivalent of considering a nonspherical nucleus) moves the locations of the significant features by a maximum of 1 km on ingress and 7 km on egress for the nominal ring plane. The possible locations in the sky plane for ring material are much broader when considering the ring-pole position errors (left panel of Figure 9).
Figure 9. Sky-plane projections of Chiron at the time of the 2018 occultation, along with the location of the observed chord for different ring-pole positions and closest approaches. The black and gray lines, which appear continous at this resolution, are plots of individual dots that each represents one image taken using the 74-inch telescope. Large red and pink dots indicate locations where the light-curve data deviate by 3σ or more below the baseline, with the chord going behind the solid nucleus in these plots. (Left) Range of ring-pole positions based on the errors on ecliptic latitude and longitude for the preferred solution in Ortiz et al. (2015). The solid gray ellipse represents the ring location for the preferred pole position (the pole indicated by a black arrow, as shown in Figure 4). Dotted gray ellipses indicate the different ring locations for the extremes of the ring-pole error bars (the poles indicated by gray arrows). The chords are plotted with closest-approach distances of ±10 km from the nominal 50 km solution, and each data point is plotted in gray with the significant light-curve deviations represented in pink. (Right) The best-fit triaxial shape for Chiron's nucleus from Braga-Ribas et al. (2023) with the updated nominal ring-pole position (not including errors) and ring concentrations (including errors) from Ortiz et al. (2023). The chord is plotted in black and red for the nominal closest approach, with gray and pink for closest approaches of ±10 km.
Download figure:
Standard image High-resolution imageLikewise, the bottom panel of Figure 5 provides an example of how the locations of significant features can vary as a function of the ring-pole position; two extreme cases at the edges of the error bars of the ring-pole position ecliptic coordinates are plotted in the bottom panel. They each require lower closest-approach distances than the nominal value in order to align the occultation by the nucleus. The locations of the significant features in the ring planes change by more than 100 km, depending on the ring-pole location. For extreme case 1, the 2018 double-dip features appear well aligned at approximately 300 and 306 km from the center. This orientation falls within the error bars of the newer Ortiz et al. (2023) ring-pole position. For extreme case 2, the nucleus would be more than 300 km in diameter in the ring plane, which is not consistent with any reported size measurements for Chiron. While we can test for locations of surrounding material given a nominal ring-pole position, more accurate ring-pole coordinates are needed to narrow down the range of possible locations and more accurately measure positional changes over time.
The occultation star in 2018 was significantly fainter than in 2011 or 2022 (see Table 1). S/Ns of occultation light curves are typically reported as the mean over the standard deviation of the normalized flux; to compare data quality between data sets, the spatial scale of the sampling must also be considered. For example, for the 2018 occultation data, the 0.5 s cadence at the 74-inch telescope had a spatial sampling of 2.98 km compared to the 1.5 s or 8.9 km sampling in the 40 inch data set. We select 3 km as a comparative spatial scale as the highest-resolution ring widths at Chiron from Sickafoose et al. (2020). The 2018 data sets have an S/N per 3 km at Chiron of 12.5 and 7.2 for the 74 and 40 inch, respectively. The S/N was insufficient to detect material in the 2018 data from the 40-inch telescope, while the material was apparent in the 74-inch data set. In 2011, the velocity was 9.87 km s−1, and the S/Ns per 3 km were 25.9 and 18.7 for the 0.2 s FTN and 0.5 s MORIS data, respectively (Sickafoose et al. 2020). In 2022, the velocity was 4.25 km s−1, and the highest-quality data set was from Kottamia Observatory with a sampling of 4.5 s and a standard deviation of 0.01 estimated from Figure 1 of Ortiz et al. (2023), corresponding to a robust S/N per 3 km of 40.
From the 2011 data, two distinct drops in flux with ring-plane widths of 2.6–4.4 km and normal optical depths of 0.5–0.85 were reported (Ruprecht et al. 2015; Sickafoose et al. 2020). The highest-quality 2011 light curve had a spatial resolution of 2 km at Chiron. The best 2018 light curve has a spatial resolution of 3 km at Chiron, and the significant drops have similar widths but at lower optical depths than those detected in 2011. The 2022 data have a resolution of only 19 km; optical depths were not reported, but the depth and breadth of flux attenuation (Ortiz et al. 2023) suggests higher optical depths in 2022 for features that would have similar widths to those observed in 2018. Indeed, the detection of surrounding material in 2022 led to the conclusion from Ortiz et al. (2023) that there was more material present than in 2011. In terms of location, the 2018 double-dip features are nominally outside of the previously proposed distance from the nucleus (see Figure 4); however, they do fall within the range of possible locations given the error bars on the ring-pole position (see the left panel of Figure 9). The 2018 features are slightly interior to the newly determined ring concentrations from Ortiz et al. (2023; see right panel of Figure 9), noting that the errors in the newer ring-pole position are similar to those from 2011 and thus likely encompass the 2018 detections. The gap between the two features remained roughly the same at 8 km in 2018 and 9–10 km in 2011 (Sickafoose et al. 2020); the resolution was too low to observe such a gap in the 2022 data (Ortiz et al. 2023). An additional single drop in flux was detected in 2018 within the proposed ring location at 316 km from the center in the nominal ring plane. Sickafoose et al. (2020) also reported one drop in flux on immersion that was interior to the proposed two-ring system, but it was broader (17.6 km ring-plane width) and more optically thick (τN = 0.29) than the third feature on immersion observed here.
The orbital distances and widths of the flux drops at Chiron are roughly similar to Chariklo's established rings, with the latter system at ∼400 km and widths of ∼3 and 7 km (Braga-Ribas et al. 2014). Note that observed values of optical depths for small-body rings have not been consistently calculated in the literature. Published values for optical depths for material surrounding Chiron (Sickafoose et al. 2020) and in this work are based on the standard relationship of I/I0 = e−τ , where I and I0 are the transmitted and incident fluxes. The normal optical depth is calculated using the opening angle, B, by for a polylayer ring (e.g., Elliot et al. 1984). Reported optical depths for Chariklo's rings, the proposed ring at Quaoar, and the recent recalculation of the characteristics of the features at Chiron have included a factor of 2, (Braga-Ribas et al. 2014; Bérard et al. 2017; Braga-Ribas et al. 2023; Morgado et al. 2023). According to Bérard et al. (2017), the factor of 2 stems from Cuzzi (1985), who found that the optical thickness measured from stellar occultations for the Uranian rings was larger by an extinction efficiency factor of 2 than the fractional area physically filled by particles (i.e., the Mie coefficient was Qe = 2; Cuzzi 1985). Typically, the Mie coefficient comes into play when inferring the physical characteristics of ring particles, which we have not attempted here. This factor of 2 was noted as the primary difference between the optical depths reported in Sickafoose et al. (2020) and the recalculated versions in Braga-Ribas et al. (2023). This difference also means that the τN values reported for Chariklo need to be increased by a factor of 2 for direct comparison with the reported values for Chiron from this work. Here, we find normal optical depths of roughly 0.2–0.4, and the comparable values for Chariklo's rings would be such that C1R is much more optically thick at τN = 0.8, and C2R is more optically thin at τN = 0.12.
Ring-type material around small bodies in the outer solar system has been discovered well beyond the Roche limit, a location beyond which basic theory would suggest that particles should disperse or accrete (e.g., Melita et al. 2017; Morgado et al. 2023; Pereira et al. 2023). The Roche radius can be defined as , where M is the body's mass, is the density of the orbiting material, and the factor describing particle shape γ = 0.85 for classical calculations, while γ = 1.6 has been preferred to represent the tidal destruction limit (e.g., Sicardy et al. 2019; Morgado et al. 2023). For Chiron's mass M = 4.8 ± 2.3 × 1018 kg (Braga-Ribas et al. 2023), water-ice particles with , and either value of γ, Chiron's Roche limit is <260 km. Considering the extreme case of the largest mass including errors and , the density of Saturnian satellites (Thomas & Helfenstein 2020), the Roche limit can reach approximately 380 km. Like Chariklo, the surrounding material at Chiron is thus very near or beyond the classical Roche limit.
The locations of the significant material surrounding Chiron in 2018 are tens of kilometers different from the proposed rings in 2011; however, they do fall within the errors of that ring-pole position. The optical depths in 2018 are lower than in 2011. Material around Chiron's nucleus has been previously observed while not in a ringlike configuration (e.g., Elliot et al. 1995; Meech et al. 1997). In addition, multiple thin and broad features were reported in 2011 by Sickafoose et al. (2020). Ortiz et al. (2023) most recently found that the structure of the surrounding material changed between 2011 and 2022. The 2018 results presented here differ enough from the proposed two-ring system from 2011 of Ortiz et al. (2015) to support the conclusion that the surrounding material at Chiron is evolving on relatively short timescales. More occultation observations are needed to continue studying Chiron and other intriguing and evolving small-body ring systems.
Acknowledgments
We thank two anonymous reviewers for comments that improved the manuscript. This work was partially supported by National Science Foundation (NSF) Astronomy and Astrophysics Research Grants award No. 2206306 and NASA grant 80NSSC21K0432. This paper uses observations made at the South African Astronomical Observatory (SAAO). A portion this work was supported by the MIT-Germany University of Stuttgart Seed Fund (https://misti.mit.edu).
Facility: SAAO - .