arXiv:1911.06271v2 [astro-ph.EP] 18 Nov 2019 h ubeSaeTlsoeoedyatrMKye al. et McKay after board day on one camera Telescope the Space a in by Hubble resolved taken the 2I/Borisov, image of disseminated coma widely inner the that noted nrdcn abnmnxd stedie fatvt in activity of in of observations. driver need is pre-discovery the the that the as without al. water, monoxide et of carbon Bolin abundance scenario introducing or high a a present al. with I et following, line McKay the by In the considered taken. not of well area is al., surface po- McKay et comet’s by total the independently of the expressed suggestion per- hyperactivity, the tential of 140 that than limit appears It more upper nucleus. of the area of active cent an finds one distance, ihBlne l’ siae pe ii f14km of et rate 1.4 McKay sublimation of with water averaged and limit an diameter 0.37 upper of nuclear estimate estimated comet’s al.’s al.’s the et on Bolin with 2019). al. et (Fitzsimmons rate production ihBlne l’ 21)wtrpouto aeestimate s rate kg production water 100 (2019) of al.’s et Bolin with of certainty yia CN/H typical coe –5 rvse ta.(09 siae nOH an estimated on 3.3 (2019) MHz of of 1667 al. time rate at et integration line production Crovisier hour OH 2–25, 15 the October of a observations from radio hand, their other the On utdi hi eemnto fapouto aeof rate re- production water a in- of of and production Sun) determination the the their 0.63 of from in AU terms sulted 11 2.38 October in was 2019 terpretation on comet 2I/Borisov the in (when oxygen atomic of line lcrncades Zdenek.Sekanina@jpl.nasa.gov. address: Electronic 2. narcn hr omncto Sknn 09)I 2019a) (Sekanina communication short recent a In obnn rvse ta.spouto aeresult rate production al.’s et Crovisier Combining ca ta.s(09 eeto fterdforbidden red the of detection (2019) al.’s et McKay eso oebr1,21;Adnu oebr1,2019 18, L November using Addendum typeset 2019; Preprint 14, November Version × × AOO AG C-UTGAN N CHUNKS AND GRAINS ICY-DUST LARGE OF HALO 10 10 ULMTO FWTRIEFO OUAINO AG,LONG-LA LARGE, OF POPULATION A FROM ICE WATER OF SUBLIMATION ae c rmtenceshsbe nraigsnetetm fdisco of o time acceleration the nongravitational since headings: measurable Subject increasing a T exert been suspec to has strongly disintegration. a enough nucleus comet’s rapid high cm the the from to to 2-3 ice contributing than subjected water halo smaller some icy-dust low initial grains and exceedingly the All was 2019 rate excessively mid-October sublimation recently. appearing by the more comet when Sun, condensation the the the for nuclear of from responsible AU th population been 5–6 that The have at suggests effect. to comparison likely sublimation This integrated is AU. the 2 wit mitigates 2I/Borisov of 2I comparing distance by perihelion investigated equal are at ice, activity water of amorphous product a grains. 27 27 e rplinLbrtr,Clfri nttt fTechno of Institute California Laboratory, Propulsion Jet oeta ffcso ulmto fwtriefo eysol oigm moving slowly very from ice water of sublimation of effects Potential oeue km molecules oeue s molecules − ± 1 7pret hsrsl svrulyidentical virtually is result This percent. 27 2 bnac ai n naalbeCN available an and ratio abundance O ae nterapoiaesaigo a of scaling approximate their on based 1. A INTRODUCTION T − × E tl mltajv 08/22/09 v. emulateapj style X 1 − 10 iha ro of error an with 2 oes niiul(IBrsv otCodcmt)—mtos da methods: — comets) Cloud Oort (2I/Borisov, individual comets: 27 s − oeue s molecules 1 ttegvnheliocentric given the at eso oebr1,21;Adnu oebr1,2019 18, November Addendum 2019; 14, November Version ERTENCESO 2I/BORISOV? OF NUCLEUS THE NEAR ∼ − 1 0A rfrhrfo h u rvnpeual yanaigof annealing by presumably driven Sun the from farther or AU 10 ± iha un- an with 4percent. 24 dnkSekanina Zdenek ABSTRACT oy 80OkGoeDie aaea A919 U.S.A. 91109, CA Pasadena, Drive, Grove Oak 4800 logy, odabd f00 n h msiiyo ..Because for assumed 0.9. one of the from emissivity differs the the regime wa- and assume sublimation 0.04 and I this of reradiation arbitrarily albedo thermal Somewhat Bond the sublimation. by ice the losses ter from the input and radiation by the Sun distance between heliocentric balance uniform given energy fairly the a a at having dictated equilibrium, mate- consisting temperature thermal refractory grains, in by contaminated are dark (i.e., rial), that ice water assume dirty I of 2I/Borisov, to sis paper earlier my of in cloud 1982). a proposed (Sekanina to grains water of sizable rates pre-existing attributed production they high 1980b); observed E1 = the C/1980 I comet by 1980 Cloud designation considered Oort old an was (Bowell; ac- for scenario (1984) is al. similar nucleus et very rel- A’Hearn small A a relatively active. a only tually of if fraction even contribute small hyperactivity thereby atively could comet’s production, water atmosphere the reported explaining its the in to per- substantially chunks to larger way and and their centimeter-sized grains on sublimating comets that and Cloud ihelion Oort Oort do the as evolution from arriving 2009). just al. amor- comets et of (Meech in annealing Cloud by ice driven water be releases to phous that believed milli- activity is beyond The a debris or characteristic this 1975). near from (Sekanina the perihelia range line with snow in a comets the Cloud commonplace (in Oort from is debris of more) tails up) of or AU type across 10 the This meter of when order times the Sun. at (on nearly- the far nucleus of still the grains was halo from comet icy-dust a released by larger pebbles occupied or and be centimeter-sized could stationary, of coma side inner projected this vector. the that orbital-velocity suggest negative a considerations in the Orbital condensation of nuclear direction the from general arcseconds extending few feature a a over shows observation, their made (2019) oeaietepasblt fapyn hshypothe- this applying of plausibility the examine To similar undergone has 2I/Borisov that propose now I h ria oino 2I/Borisov. of motion orbital the n oia otCodcmtof comet Cloud Oort nominal a h e yeatvt.Sbiainof Sublimation hyperactivity. ted togyhproi oinof motion hyperbolic strongly e egan ertenceso 2I of nucleus the near grains se ey u h aehsntbeen not has rate the but very, rgti r-icvr images pre-discovery in bright i etol agrcuk of chunks larger only left his swl sfrteprominent the for as well as , rs a endevolatilized been had cross liee-ie n larger and illimeter-sized TN GRAINS STING aanalysis ta 2 Sekanina the nucleus by McKay et al. (2019), the sublimation area Λ˙ of the putative source of detected water is also different, (cm/day) COLUMNAR SUBLIMATION RATE 2 2 equaling 5.8 km instead of 1.7 km . OF OORT CLOUD COMET ······················································· Let the sublimation rate of water ice implied by the en- −1 AND 2I/BORISOV ······································· 10 ·························· ········· ergy balance at time t, when the comet is at heliocentric ············ ······ − − ······· ···· distance r(t), be Z˙ (t), expressed in molecules cm 2 s 1. ······· ···· ········· ····· Next, I introduce a columnar sublimation rate, dΛ(t)/dt ········ ···· −2 ····· ··· or Λ(˙ t), at which the thickness of a layer or the length of 10 ······· ···· −1 ····· ··· a column of water ice contracts (in cm s , for example) ······· ···· OORT CLOUD ···· ·· due to its sublimation, by ···· ··· 2I/BORISOV COMET ····· ··· −3 ···· ·· µZ˙ (t) 10 ······ ··· Λ(˙ t)= , (1) ····· ··· δ ······ ··· ······ ··· where µ is the mass of a water molecule and δ is the bulk ······ ··· −4 ··· ·· density of water ice in the grain. Integrating (1) from an 10 ····· ··· ······ ··· early time before the sublimation began, the length of ······ ··· a column of water ice that has sublimated away by a ······ ··· ······· ···· reference time, tref , is −5 ···· ·· 10 ···· ··· tref rref · µ −1 Λ(tref)= Λ(˙ t) dt = Z˙ (r)r ˙ dr, (2) −∞ δ ∞ −300 −200 −100 0 Z Z TIME FROM PERIHELION AT 2.006 AU (days) where tref is reckoned from perihelion, tπ, and rref is the Figure 1. Comparison of preperihelion columnar rates of water heliocentric distance at time tref. In general, there are no ice sublimation from grains in the atmospheres of comet 2I/Borisov constraints on tref , but in the following I only consider the and a nominal Oort Cloud comet of equal perihelion distance of contraction of a column of water ice by its sublimation 2.006 AU. Note that because of the higher orbital velocity, 2I’s rate is lower by one order of magnitude 90 days before perihelion, along the preperihelion leg of the orbit, in which case but by five orders of magnitude 140 days before perihelion. tref tπ 0 andr ˙ 0. To− facilitate≤ the≤ integration, I introduce a new, dimen- sionless variable, 0 z 1, by substituting ≤ ≤ for 2I/Borisov and by τ = 300 days for the nominal Oort τ 0 z = exp , (3) Cloud comet in the parabolic orbit. τ0 An obvious choice for approximating a function de- fined by a sequence of (x, y) data pairs is a polynomial where τ = t tπ and τ0 > 0 is a constant. The length of a column of− water ice sublimated by time t is then y(x), whose coefficients are determined by least squares. ref The advantage of this approach is twofold: first, one can exp(τref /τ0) optimize the polynomial’s degree by searching for the Λ(tref )= τ0 (t) dz, (4) lowest possible degree that secures the requested quality F Z0 of fit; and second, the integration is straightforward. To where τref = tref tπ and follow this procedure, I write − n Λ(˙ t) τ k (t)= = Λ(˙ t) exp . (5) j(t)= ak,j z , (6) F z −τ0 F Xk=0 The columnar rates Λ(˙ t) of water ice sublimation for where j = 1 for 2I/Borisov and j = 2 for the Oort Cloud 2I/Borisov and an Oort Cloud comet of the perihelion comet. In accordance with this notation, I refer from distance of 2.006 AU, are compared in Figure 1 as a func- now on to Λ and τ0 for 2I/Borisov as Λ1 and τ0,1, for tion of time reckoned from perihelion. The lower rates the Oort Cloud comet as Λ2 and τ0.2, respectively. Ex- for 2I/Borisov stem from the higher orbital velocity in perimentation with a number of reference points (z, ) its strongly hyperbolic path, i.e., at an equal time from has shown that a cubic polynomial fits both 2I/BorisovF perihelion 2I is farther from the Sun than the Oort Cloud and the Oort Cloud comet quite adequately, if a typi- −1 comet of equal perihelion distance. cal error of approximately 0.001 cm day for j is The length of a column of water ice sublimated by the acceptable, but only over the± range of heliocentricF dis- reference time, defined by Equation (2), is obviously also tances not exceeding r =3.4 AU, equivalent to about ∗ ∗ shorter for 2I than for the nominal Oort Cloud comet. z 0.43 and t tπ = 125 days for 2I/Borisov and to 1∗ ≃ ∗1 − − To determine this quantitatively, the integration of the z2 0.45 and t∗2 tπ = 240 days for the nominal Oort expression on the right-hand side of Equation (2) or (4) Cloud≃ comet. The− failure− of any fit of this kind for all requires decisions on two issues: (i) to devise an appro- values of z smaller (and times from perihelion larger) priate method of integration; and (ii) to choose the con- than these limits has two major implications. One is the stant τ0 separately for 2I and the nominal comet in the need to devise an integration approach other than via parabolic orbit once Equation (4) is used. Equation (4) for all reference times preceding t∗ (i.e., The choice of τ0 is not at all critical, the only concern for heliocentric distances larger than r∗). And two, all being that, for the purpose of smooth integration, is reference times between t∗ and the perihelion time tπ approximately flat near perihelion, where Λ˙ reachesF its call obviously for a revision of the integration limits in maximum. This condition is satisfied by τ0 = 150 days Equation (4). Sublimation of Water from Grains in 2I/Borisov? 3
′ The remedy of the first problem is greatly facilitated by polynomials, Λj , I find (dropping the subscript ref ) the fact that at heliocentric distances exceeding 3.4 AU the incident solar radiation is spent overwhelmingly on exp(τj /τ0) Λj (t)= τ0,j j dz +Λj(t ) reradiation, so that the columnar rate of water ice subli- ∗ ∗j exp(τ /τ0) F mation follows closely a relation Z j (14) ′ ′ ˙ = τ0,j Λ (t) Λ (t )] + o[Λj (t) , j =1, 2, Λ= A exp B√r , (7) j − j ∗j − 1 13 −1 − 2 where A =6.7 10 cm day andB = 21.587 AU . where τj∗ = t∗j tπ, o[x] means a negligible contribution One bypasses Equation× (4) by integrating directly Equa- to x, and again− j = 1 for 2I/Borisov and j = 2 for the tion (2) to obtain for the length of the sublimated column nominal Oort Cloud comet. of water ice by time t∗ an expression I already remarked that cubic polynomials represent ∗ adequate approximations for both objects; the columns t ′ Λj sublimated by time t (t∗j