HTTP_Request2_Exception Unable to connect to tcp://think-ws.ca.edps.org:85. Error: php_network_getaddresses: getaddrinfo failed: Name or service not known CO-driven activity constrains the origin of comets | Astronomy & Astrophysics (A&A)
Free Access
Issue
A&A
Volume 636, April 2020
Article Number L3
Number of page(s) 3
Section Letters to the Editor
DOI https://doi.org/10.1051/0004-6361/202037805
Published online 07 April 2020

© ESO 2020

1. Introduction

While the activity of comets inside the Jupiter orbit is mainly driven by water-ice sublimation, beyond Saturn, this activity is driven by the so-called supervolatiles, of which CO is the most abundant (Bockelee-Morvan et al. 2004). A detailed characterization of the CO-driven activity of comets is needed also in view of the ESA mission Comet Interceptor (Snodgrass & Jones 2019), which will target an Oort-cloud (or even interstellar) comet discovered as far as possible from the Sun. These models will support the choice of the probe trajectory to best fly-by its target.

Of the comets that are active beyond Saturn the incoming comet C/2017K2 PANSTARRS was discovered at 16 au from the Sun, with prediscovery images at 24 au, when it was already active, and a period of ≈3 Myr consistent with inner temperatures of a few K, typical of Oort comets before the first Sun approach (Jewitt et al. 2019). Models of the collected data suggest an activity onset at about 26 au from the Sun and a dust coma made of particles larger than 0.1 mm, released at a steady rate of ≈200 kg s−1 (Jewitt et al. 2019). Beyond Saturn, the gas temperature and density at the nucleus surface are too low to detach the observed dust from the nucleus; tensile strengths > 0.1 Pa bond each particle to the nucleus (Skorov & Blum 2012). Therefore, exotic processes (thermal stress or electrostatic charging) have been invoked to explain the dust ejection (Jewitt et al. 2019).

The dust cohesion bottleneck is a well-known paradox in cometary physics and has been explained only recently for a water-driven cometary activity (Fulle et al. 2020). Here, we apply this model to CO-driven activity, showing that the CO-gas pressure inside the pebbles of which cometary nuclei consist (Blum et al. 2017) explains the ejection of the observed dust and fits all the observed data. Here we make the following assumptions for CO: (i) its fraction trapped in other ices and clathrates is negligible, and (ii) it is present as pure ice, thus sublimating beyond Saturn. We find that the Oort comet C/2017K2 is similar to all Jupiter-family comets: it has high values of the refractory-to-all-ices mass ratio and of dust fallout (Fulle et al. 2019). Our results also constrain the evolution of cometary planetesimals.

2. Activity model

The first thermal model consistent with dust ejection from cometary nuclei (Fulle et al. 2020) is applied here to the specific case of a comet with CO-driven activity at heliocentric distances that are too large to allow any significant sublimation of H2O and CO2 ices. The model is consistent with a nucleus consisting of centimeter-sized pebbles (Blum et al. 2017) where ice sublimation occurs only inside the pebbles (Fulle et al. 2020). Pebbles are inhomogeneous clusters of dust particles, which are porous aggregates of ice and dust grains (Levasseur-Regourd et al. 2018; Güttler et al. 2019) and have a differential size distribution with a power index of −3 ± 1 (Blum et al. 2017). In this specific case, dust is composed not only of refractories, as usual, but also of H2O and CO2 ices. All the governing equations we apply here with parameters valid for CO have been obtained for a water-driven activity (Fulle et al. 2020). CO gas diffuses inside the pebbles, with a pressure inside the pebbles, P(s), and a CO-gas flux from the nucleus surface, Q, given by (Fulle et al. 2020)

P ( s ) = P 0 f ( s ) exp [ T 0 T s s T ] $$ \begin{aligned}&P(s) = P_0~f(s)~\exp \left[{- {T_0 \over {T_{\rm s} -\,s~\nabla T}}}\right] \end{aligned} $$(1)

Q = 14 r P ( R ) 3 R 2 m π k ( T s R T ) , $$ \begin{aligned}&Q = {{14~r~P(R)} \over {3~R}} \sqrt{{2~m} \over {\pi k~(T_{\rm s} - R~\nabla T)}}, \end{aligned} $$(2)

where s is the depth from the nucleus surface and f ( s ) = 1 ( 1 s R ) 4 $ f(s) = 1 - (1 - {s \over R})^4 $ for s ≤ R, f(s) = 1 elsewhere. In the range of temperatures 40 ≤ Ts ≤ 50 K of the nucleus surface, the values P0 = 1.73 × 1010 Pa and T0 = 942 K fit the available CO-ice sublimation data best (Fray & Schmitt 2009). ∇T is the nucleus average gradient of the temperature at P ≥ 0.1 Pa, r ≈ 50 nm and R ≈ 5 mm are the radii of the grains of which cometary dust consists (Levasseur-Regourd et al. 2018; Mannel et al. 2019; Güttler et al. 2019) and of the pebbles of which cometary nuclei consist (Blum et al. 2017; Fulle et al. 2020), m is the mass of the CO molecule, and k is the Boltzmann constant (all symbols are listed in Table 1).

Table 1.

Symbols.

The value of Ts is given by the energy balance equation,

( 1 A ) I cos θ r h 2 = ϵ σ T s 4 + λ s T + Λ Q , $$ \begin{aligned} (1 - A)~I_\odot ~\cos \theta ~r_{\rm h}^{-2} = \epsilon \sigma T_{\rm s}^4 + \lambda _{\rm s}~\nabla T + \Lambda ~Q, \end{aligned} $$(3)

where A is the nucleus Bond albedo (e.g., A = 1.2% measured at 67P; Fornasier et al. 2015), I is the solar flux at the heliocentric distance of Earth, θ is the solar zenithal angle, rh is the heliocentric distance of the nucleus in astronomical units, ϵ ≈ 0.9 is the emissivity, σ is the Stefan–Boltzmann constant, λs is the heat conductivity inside the nucleus, and Λ = 2.27 × 105 J kg−1 is the latent heat of sublimation of CO ice (Huebner et al. 2006).

The values of λs and ∇T are provided by (Fulle et al. 2020)

T = Λ Q / σ 8 ( T s R T ) R $$ \begin{aligned}&\nabla T = {\sqrt{\Lambda ~Q~/~\sigma } \over {8~(T_{\rm s} - R ~\nabla T)~R}} \end{aligned} $$(4)

λ s = 32 3 σ R ( T s R T ) 3 $$ \begin{aligned}&\lambda _{\rm s} = {32 \over 3}~\sigma ~R~(T_{\rm s} - R~\nabla T)^3\end{aligned} $$(5)

because in a pebble-made nucleus the heat conductivity is dominated by the radiation among the pebbles (Blum et al. 2017). Equations (1)–(5) must be solved altogether, fixing Ts, ∇T, λs, and Q. The tensile strength, S, which bonds the superficial dust particles to the pebble, is S = 13 s−2/3 mPa (with s in meters; Skorov & Blum 2012, Fig. 1). Because the ejected dust is detached from the nucleus surface, we can identify the dust size with the depth s. P ≥ S defines the minimum, sm, and maximum, sM, ejected dust sizes. The nucleus at the depth sM maintains a constant temperature Ts −  sMT ≈ 38.5 K (Fig. 1), so that dust ejection occurs in thermal quasi-equilibrium. Therefore the rates at which the nucleus is eroded, E, and is CO-depleted, D, are (Fulle et al. 2020)

E = 32 ( 1 + δ ) σ R ( T s s M T ) 3 3 ( 1 2 + δ ) ρ d c p s M $$ \begin{aligned} E&= {32~(1 + \delta )~\sigma ~R~(T_{\rm s} -\,s_{\rm M}~\nabla T)^3 \over {3 ~({1 \over 2} + \delta )~\rho _{\rm d}~c_{\rm p}~s_{\rm M}}} \end{aligned} $$(6)

D = ( 1 + δ ) γ Q ρ n , $$ \begin{aligned} D&= {{(1 + \delta )~\gamma ~Q} \over \rho _{\rm n}}, \end{aligned} $$(7)

thumbnail Fig. 1.

Pressure P inside the pebbles vs. depth s. The continuous line shows the tensile strength, S, bonding a homogeneous aggregate of dust grains. The dashed line represents P at Ts = 48 K and ∇T = 12 K cm−1, which crosses S at sm = 0.12 mm and sM = 8 mm. The dotted line shows P at Ts = 40 K and ∇T = 4 K cm−1 and is tangent to S at sm =  sM = 3 mm. For every value of Ts, the nucleus at the depth sM maintains a temperature Ts − ∇TsM = 38.5 K.

where ρd ≈ 800 kg m−3 is the average bulk density of the refractory particles (Fulle et al. 2017), cp ≈ 103 J kg−1 K−1 is the heat capacity of the pebbles (Blum et al. 2017), δ and γ are the average refractory-to-all-ices and all-ices-to-CO mass ratios in the nucleus, respectively, and ρn = 538 kg m−3 is the average bulk density of the nucleus (Pätzold et al. 2019). Here, with all-ices we mean H2O, CO2, and CO ices, where the first two behave as refractories at T <  50 K and have a bulk density close to half of that of the compacted refractories (Fulle et al. 2017). Measurements in comets Hale–Bopp and Hyakutake provided the gas mass fractions CO/H2O ≈40% and CO2/H2O ≈15% (Bockelee-Morvan et al. 2004), suggesting γ ≈ 4.

When D <  E, the CO-driven activity never stops because the CO-gas diffusion erodes the superficial pebbles before they become CO-depleted, exposing the CO-ice-rich ones underneath. In contrast, when D >  E, the CO-driven activity stops because the CO-depletion builds up an insulating crust before it is eroded by the dust ejection.

3. Discussion

The activity onset occurs when the P curve in Fig. 1 is tangent to the S line, that is when Eqs. (1)–(5) provide Ts = 40 K, ∇T = 4 K cm−1, Q = 10−7 kg m−2 s−1, λs = 1.7 10−4 W m−1 K−1, and sm =  sM = 3 mm. These values in Eq. (3) provide rh ≈ 85 au if cos θ ≈ 1. The observed activity onset at rh ≈ 26 au suggests cos θ ≈ 0.1, that is, a nucleus that is inactive in all areas with 0.1 ≪ cos θ ≤ 1. For instance, if the nucleus spin axis points to the Sun, then the active area would be a tight equatorial belt. Conversely, if the nucleus spin axis has a low obliquity, then the active area would be a small polar spot. The orbital period of ≈3 Myr makes the orbit sensitive to galactic tides and perturbations by Jupiter, with a nucleus spin obliquity probably different in the past, when subsolar activity may have covered all the areas that now have 0.1 ≪ cos θ ≤ 1 of inactive fallout. A fallout of all the dust ejected beyond 26 au from the Sun is improbable because the CO flux increases by a factor 10 only from the activity onset up to the observations. According to Eq. (6), the erosion rate remains about constant (within a factor 3) since the activity onset, explaining the steady brightness slope of the C/2017K2 dust coma (Jewitt et al. 2019).

With cos θ ≈ 0.1, Eqs. (1)–(5) at rh = 14 au provide Ts = 48 K, ∇T = 12 K cm−1, Q = 10−6 kg m−2 s−1, λs = 2.3 × 10−4 W m−1 K−1, sm = 0.12 mm, and sM = 8 mm (Fig. 1). The obtained value of sm is consistent with the available observations at rh = 14 au (Jewitt et al. 2019). The CO-driven activity never stops if δ ≤ 4.5 in Eqs. (6) and (7) (i.e., if the refractory-to-water-ice mass ratio is ≤7), in which case D ≤ E = 3 mm day−1. When no dust fallout is assumed, the observed dust-loss rate Qd ≈ 200 kg s−1 (Jewitt et al. 2019) implies a nucleus active area of ( 1 + δ ) Q d / [ ( 1 2 + δ ) ρ d E ] 7 $ (1 + \delta) Q_{\mathrm{d}} / [({1 \over 2} + \delta) \rho_{\mathrm{d}} E] \approx 7 $ km2, with a CO-gas loss rate of ≈7 kg s−1. However, the computed Q-values are too low to drag all the dust of sizes sm ≤  s ≤  sM up to the escape velocity (Zakharov et al. 2018). If the fallout amounts to ≥30% of the eroded dust mass (a fraction much lower than occurred on 67P/Churyumov–Gerasimenko; Fulle et al. 2019), then the nucleus active area becomes ≥10 km2, with a CO-gas loss rate ≥10 kg s−1, and a dust-to-gas mass ratio in the lost material ≤20, consistent with the average dust-to-CO-ice mass ratio (1 + δ)γ − 1 ≤ 21 in the nucleus (Fulle et al. 2019).

The CO-driven erosion of 3 mm day−1 implies that cometary planetesimals, with radius ≤50 km, a pristine season-independent activity, and buildup at heliocentric distances < 85 au by the gravitational collapse of centimeter-sized pebbles (Blum et al. 2017), were completely eroded within 50 kyr after the protoplanetary disk became optically thin to solar radiation. This erosion lifetime is 103 times shorter than the timescale of catastrophic collisions that occurred in the first 0.5 Gyr after planet formation (Morbidelli & Rickman 2015), which in turn excludes that comets were destroyed by collisions before a complete erosion. This suggests a natural source of the observed CO and debris disks (Hughes et al. 2017), that is, CO-driven cometary activity, which eroded most cometary planetesimals into millimeter-sized dust, with a Poynting–Robertson timescale of the orbital collapse of ≈1 Gyr (Burns et al. 1979). Because cometary planetesimals probably formed after the giant planets (Fulle & Blum 2017), the only ones surviving into comets were those that closely encountered migrating planets before a complete erosion, with a storage into orbits of perihelion > 85 au, where the timescale of catastrophic collisions is > 5 Gyr (Morbidelli & Rickman 2015). It follows that comets observed today cannot be products of catastrophic collisions, consistent with the constraints provided by the Rosetta mission (Fulle & Blum 2017).

4. Conclusions

The activity model based on gas diffusion inside the pebbles of which cometary nuclei consist (Fulle et al. 2020) was applied here to comet C/2017K2 PANSTARRS, whose activity is driven by CO-ice sublimation assuming CO as pure ice and neglecting the fraction that is trapped in other ices and clathrates. We list our results below.

(i) CO-gas pressure inside the pebbles reaches 20 Pa and erodes the nucleus surface into ejected dust (composed of refractories, H2O ice, and CO2 ice), without invoking other processes to explain the observed dust coma (Jewitt et al. 2019).

(ii) CO-driven activity onset occurs up to heliocentric distances of 85 au, ejecting dust of ≈3 mm size.

(iii) The activity onset of C/2017K2, observed at ≈26 au (Jewitt et al. 2019), suggests a low obliquity of the nucleus spin axis if the activity mainly occurs in a polar summer.

(iv) At 14 au, the smallest size of the ejected dust is ≈0.1 mm, consistent with observations (Jewitt et al. 2019).

(v) The observed dust-loss rate of ≈200 kg s−1 (Jewitt et al. 2019) implies a fallout ≥30%, a nucleus surface active area ≥10 km2, and a CO-gas loss rate ≥10 kg s−1.

(vi) CO-driven activity never stops if the average refractory-to-all-ices mass ratio in the nucleus is δ ≤ 4.5 for a nucleus all-ices-to-CO mass ratio γ ≈ 4, as observed in comets Hale–Bopp and Hyakutake. This δ-value is consistent with the value found for comet 67P (Fulle et al. 2020) and with the constraints given by pebble accretion models in protoplanetary disks (Lorek et al. 2016), namely δ ≥ 3.

(vii) The CO-driven erosion lifetime of cometary planetesimals with δ ≤ 4.5 and a pristine season-independent activity is a factor 103 shorter than the timescale of catastrophic collisions (Morbidelli & Rickman 2015). This means that comets observed today cannot be products of catastrophic collisions (Fulle & Blum 2017).

Acknowledgments

We thank an anonymous referee for having improved the first version of this letter. Part of this research was supported by the ESA Express Procurement (EXPRO) RFP for IPL-PSS/JD/190.2016 and by the Italian Space Agency (ASI) within the ASI-INAF agreements I/032/05/0 and I/024/12/0. J.B. thanks the Deutsche Forschungsgemeinschaft for support under grant BL298/26-1 as part of the international CoPhyLab collaboration. CoPhyLab is jointly funded through DFG (Germany), FWF (Austria) and SNF (Switzerland).

References

  1. Blum, J., Gundlach, B., Krause, M., et al. 2017, MNRAS, 469, S755 [CrossRef] [Google Scholar]
  2. Burns, J. A., Lamy, P. L., & Soter, S. 1979, Icarus, 40, 1 [NASA ADS] [CrossRef] [Google Scholar]
  3. Bockelee-Morvan, D., Crovisier, J., Mumma, M. J., & Weaver, H. A. 2004, in Comets II, eds. M. C. Festou, H. U. Keller, & H. A. Weaver (Tucson: Univ. of Arizona Press, Tucson), 391 [Google Scholar]
  4. Fornasier, S., Hasselmann, P. H., Barucci, M. A., et al. 2015, A&A, 583, A30 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  5. Fray, N., & Schmitt, B. 2009, Plan. Space Sci., 57, 2053 [NASA ADS] [CrossRef] [Google Scholar]
  6. Fulle, M., & Blum, J. 2017, MNRAS, 469, S39 [CrossRef] [Google Scholar]
  7. Fulle, M., Della Corte, V., Rotundi, A., et al. 2017, MNRAS, 469, S45 [Google Scholar]
  8. Fulle, M., Blum, J., Green, S., et al. 2019, MNRAS, 482, 3326 [NASA ADS] [CrossRef] [Google Scholar]
  9. Fulle, M., Blum, J., Rotundi, A., et al. 2020, MNRAS, 493, 4039 [NASA ADS] [CrossRef] [Google Scholar]
  10. Güttler, C., Mannel, T., Rotundi, A., et al. 2019, A&A, 630, A24 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  11. Huebner, W. F., Benkhoff, J., Capria, M. T., et al. 2006, The International Space Science Institute, Bern, Switzerland (Noordwijk, The Netherlands: ESA Publications Division) [Google Scholar]
  12. Hughes, A. M., Lieman-Sifry, J., Flaherty, K. M., et al. 2017, ApJ, 839, 86 [NASA ADS] [CrossRef] [Google Scholar]
  13. Jewitt, D., Agarwal, J., Hui, M.-T., et al. 2019, AJ, 157, 65 [NASA ADS] [CrossRef] [Google Scholar]
  14. Levasseur-Regourd, A.-C., Agarwal, J., Cottin, H., et al. 2018, Space Sci. Rev., 214, 64 [NASA ADS] [CrossRef] [Google Scholar]
  15. Lorek, S., Gundlach, B., Lacerda, P., & Blum, J. 2016, A&A, 587, A128 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  16. Mannel, T., Bentley, M. S., Boakes, P. D., et al. 2019, A&A, 630, A26 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  17. Morbidelli, A., & Rickman, H. 2015, A&A, 583, A43 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  18. Pätzold, M., Andert, T. P., Hahn, M., et al. 2019, MNRAS, 483, 2337 [NASA ADS] [CrossRef] [Google Scholar]
  19. Snodgrass, C., & Jones, G. H. 2019, Nat. Comm., 10, 5418 [Google Scholar]
  20. Skorov, Y. V., & Blum, J. 2012, Icarus, 221, 1 [NASA ADS] [CrossRef] [Google Scholar]
  21. Zakharov, V. V., Ivanovski, S. L., Crifo, J. F., et al. 2018, Icarus, 312, 121 [NASA ADS] [CrossRef] [Google Scholar]

All Tables

All Figures

thumbnail Fig. 1.

Pressure P inside the pebbles vs. depth s. The continuous line shows the tensile strength, S, bonding a homogeneous aggregate of dust grains. The dashed line represents P at Ts = 48 K and ∇T = 12 K cm−1, which crosses S at sm = 0.12 mm and sM = 8 mm. The dotted line shows P at Ts = 40 K and ∇T = 4 K cm−1 and is tangent to S at sm =  sM = 3 mm. For every value of Ts, the nucleus at the depth sM maintains a temperature Ts − ∇TsM = 38.5 K.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.