Exocomets size distribution in the $$beta$$ β Pictoris planetary system | Scientific Reports – Nature.com

https://ift.tt/tlcmqBa

Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

Advertisement

*Scientific Reports* **volume 12**, Article number: 5855 (2022) Cite this article

1527 Accesses

311 Altmetric

Metrics details

The star (beta) Pictoris harbors a young planetary system of about 20 million years old, which is characterized by the presence of a gaseous and dusty debris disk, at least two massive planets and many minor bodies. For more than thirty years, exocomets transiting the star have been detected using spectroscopy, probing the gaseous part of the cometary comas and tails. The detection of the dusty component of the tails can be performed through photometric observations of the transits. Since 2018, the Transiting Exoplanet Survey Satellite has observed (beta) Pic for a total of 156 days. Here we report an analysis of the TESS photometric data set with the identification of a total of 30 transits of exocomets. Our statistical analysis shows that the number of transiting exocomet events (*N*) as a function of the absorption depth (*AD*) in the light curve follows a power law in the form (dN(AD) propto AD^{-alpha }), where (alpha =2.3pm 0.4). This distribution of absorption depth leads to a differential comet size distribution proportional to (R^{-gamma }), where (gamma =3.6 pm 0.8), showing a striking similarity to the size distribution of comets in the Solar system and the distribution of a collisionally relaxed population ((gamma _}= 3.5)).

Since the mid-1980’s, spectroscopic observations of the bright southern star (beta) Pictoris, or (beta) Pic, have revealed variations in the calcium H and K lines, which have been interpreted as the transit of the gaseous tails of exocomets^{1,2,3,4,5}. Several decades of exocomets transit observations, mostly on (beta) Pic but also on a few other stars, provided original insight on the activity of minor bodies in planetary system within the first tens of million years^{6}. For that, the (beta) Pic system is exceptional in presenting simultaneously several favorable properties : it is only a few tens million years old, very nearby at 19.3 parsecs and seen exactly edge-on from the Earth^{7,8,9,10,11,12,13,14,15,16}. Nonetheless, even for (beta) Pic, the observations of exocomets were limited to the gaseous part of the tails probed by transit spectroscopy; observations of the dusty part of the tails indeed require high accuracy photometric measurements, which was until recently not available^{17,18,19,20,21}.

The situation changed when, using data obtained by the Transiting Exoplanet Survey Satellite (TESS^{21}) between October 2018 and February 2019, three photometric events were discovered^{20} and interpreted as due to the transit of the dust component of exocomets transiting (beta) Pic. In support of this interpretation, the observed light curves are almost identical to the predictions made twenty years earlier^{17,18}.

Since 2019, (beta) Pic has been re-observed by TESS. Here we present the analysis of the complete data set gathered up to February 2021 in order to perform a deep search for exocomet transits and determine the size distribution of the (beta) Pic comets to be compared with the distribution observed in the Solar system.

Considering all TESS observations of (beta) Pic until February 2021, after the cleaning process of the (delta) Scuti variations and other slower variations (see Methods) we end up with a 156 days light curve clearly showing some dips, which are similar to what is expected for exocomet transits (Extended Data Fig. 4). To make sure that the observed dips in the (beta) Pic light curve are real and not noise residuals nor artifacts due to the reduction process, we used the TESS observations of the nearby star (alpha) Pic. We downloaded the PDC flux time series of (alpha) Pic and applied the same procedure to remove the (delta) Scuti and slower variations. The observations of (alpha) Pic provide an excellent data-set to test our procedure because (alpha) Pic is in the same region of the sky as (beta) Pic (and has hence overlapping TESS observations epochs), it has the same spectral type (A8V versus A5V for (beta) Pic), similar magnitude (3.3 versus 3.85 for (beta) Pic) and a similarly high (vsin i) projected rotational velocity (205 km/s for (alpha) Pic versus 139 km/s for (beta) Pic). Thus, (alpha) Pic is almost a nearby stellar twin of (beta) Pic except for the presence of the young planetary system. (alpha) Pic has already been successfully used as a reference star for analysis of (beta) Pic observations^{22}.

The (alpha) Pic light curve is used to check for systematics that could mimic transit of exocomets. The light curve of (alpha) Pic shows noisy excursions from the mean value both in the positive and negative deviations with the same pattern. The light curve of (beta) Pic also shows noisy excursions in the positive direction that are similar to the ones observed in the light curve of (alpha) Pic, but noisy excursions that are more pronounced and more frequent in the negative direction. The latest are typical signatures of small, transiting objects with extended dust tails. We do not detect similar variations in (alpha) Pic that the ones seen in (beta) Pic that can be due to the passage of exocomets in front of the star.

To characterize this excess of transit-like features in the light curve of (beta) Pic and identify the corresponding individual events, we calculated the correlation of the light curves with a simple model of an exocomet transit photometric event, assuming a 1D transit of a translucent dust cloud with an exponential decrease of the optical thickness from the head of the comet. This model has only four parameters : *K*, the cloud optical thickness at the leading head, (Delta t), the transit duration corresponding to the time needed to cover the chord length of the stellar disk at the transit velocity, (beta), the speed of the transit of one scale length of the cometary tail (the inverse of the scale length of the cometary tail divided by the transit velocity), and finally (t_0), the time of the beginning of the transit. With this model of the transit of the exocometary tail, the corresponding decrease in relative flux at the time *t* is given by (Delta F/F, (t) = K times left( exp (-Delta ) – exp (-Delta ^prime )right)), where (Delta =beta (t-t_0)) if (tge t_0), (Delta =0) if (tle t_0), and (Delta ^prime =beta (t-t_0-Delta t)) if (tge t_0 +Delta t), and (Delta ^prime =0) if (t le t_0 +Delta t).

We calculated the correlation of this 1D-model with the observed light curves of (beta) Pic and (alpha) Pic by varying the value of the transit time (t_0) and using various plausible values for the (Delta t) and (beta) parameters characterizing the shapes of the exocomet transit light curves. We used (beta) from 2 to 20 (hbox {days}^{-1}) and (Delta t) from 0.15 to 0.5 days, corresponding to periastron distances ranging from about 0.08 to 0.85 au. The value of *K* in the model can be arbitrarily chosen because it only changes the amplitude of the light variations, hence it has no consequence in the position of the correlation maximum nor on the identification of exocomet transits events.

In the case of (beta) Pic the correlation reaches large positive values for some of the values of (t_0). This behaviour is not observed with the (alpha) Pic light curve, showing that the photometric variations with cometary transit shapes are specific to (beta) Pic. We check the negative values of the correlation of the model with the (beta) Pic light curve itself, and found that they are much less numerous and significantly smaller than the positive values. This confirms that the light curve of (beta) Pic shows photometric variations with decrease of the star brightness that are typical of the transits of exocomets and the absence of variations with increase of the star brightness with similar shape and amplitude. We interpret the correlation peaks as the signatures of potential exocomet transits. We keep only the peaks with correlation values that are higher than the maximum value obtained with the (alpha) Pic data and the maximum negative value with the (beta) Pic data. This is our conservative criterion to consider the observed variations as a detection of an exocomet transit.

We identified a total of 30 significant detections of exocomet transit events. Although this is not surprising, acknowledging the well known ubiquity of comets in that young planetary system, this is the first time such a large number of exocomets are detected in photometry. This allows statistical analysis of their properties. The light curves of these detected exocomet transits are plotted in Extended Data Fig. 4. Their characteristics and the parameters of the best fits with the 1-D model are given in Extended Data Table 2. The last column of the table gives the square root of the (chi ^2) improvements of the fits to the light curve between a model with no transit and the 1-D transit model using the best fit parameters. This shows that the identified exocomet transits are all detected at least at 4-(sigma) level.

Our detection of 30 photometric transits of exocomets allows a statistical analysis of their properties. Here we call “absorption depth”, noted *AD*, the decrease in relative flux at the minimum of a transit light curve. The numerical value of the absorption depth is estimated from the best fit with the 1-D model with (AD=K(1-exp (-beta Delta t))). A plot of the events frequency as a function of the absorption depth shows that there is a steep decrease of the number of events toward the larger absorption depths (Fig. 1). The differential number of transiting exocomet events (*dN*) detected with an observation of duration (delta t) as a function of the absorption depth can be fitted by a power law in the form (dN(AD)=N_{0}cdot (AD/10^{-4})^{-alpha }cdot (delta t/100, })cdot (dAD/10^{-4})). Considering 29 events detected in 156 days of observations, that is all the 30 events except the deepest one (see below), we find that (alpha =2.3 pm 0.4) and (N_{0}=33^{+16}_{-11}), where the uncertainties have been evaluated using a Poisson distribution for the number of events in each bin of width (d AD=1.5cdot 10^{-4}).

The deepest transit event^{20} of Julian Day JD = 2457000 + 1486 looks exceptional with an absorption depth of about (20times 10^{-4}). This event could be produced by a member of another family of exocomets than the one which produces the 29 other shallower events, as we know from spectroscopic transit observations the presence of several families of (beta) Pic comets^{5}. Nonetheless, with the distribution derived above, the expected number of events with absorption depths in the range [10–20](times 10^{-4}) in a 156 days observation is (1.2^{+0.6}_{-0.4}). Even in the range [15–20](times 10^{-4}) the expected number of detections in 156 days is (0.38^{+0.19}_{-0.12}), corresponding to a probability of 26% to have one single event in this range as observed. Therefore the event of JD = 2457000 + 1486 can simply be a normal event, only the deepest, within the same distribution of the other 29 events. If this event is taken into account, the estimates of the distribution power law indexes (alpha) and (gamma) change by about 0.1 and 0.2, towards slightly shallower distributions. Nonetheless, to remain conservative, this deepest event is not taken into account in the derivation of the distribution of absorption depths and exocomets sizes considered here.

The modeling of the exocomet transit light curves shows that the transit absorption depth, *AD*, is directly proportional to ({dot{M}}), the dust evaporation rate from the comet nucleus^{17} . If we assume that the dust production rate is proportional to the comet nucleus area (*i.e.*, Ref.^{23}), we have a production rate ({dot{M}}) proportional to the squared of the nucleus radius *R*. Finally, with *AD* proportional to (R^2), we find that the differential number of exocomets as a function of the nucleus size is given by (dN(R) propto R^{-gamma } dR), with (gamma = 2alpha -1).

With the fit to the observed distribution of the absorption depths, we conclude that the differential distribution of the exocomet size must follow a power law with an index (gamma =3.6 pm 0.8). This distribution is notably similar to the size distribution of comets in the Solar system (Fig. 2) and the distribution predicted in Ref.^{24} for a collisionally relaxed population ((gamma _}= 3.5)).

For the plot of the size distribution (Fig. 2), we used the cometary radii estimated following the derivation described in the Method section. The conclusion on the similarity of the size distributions in (beta) Pic and the Solar system is independent of these absolute size estimates. Nonetheless, it is remarkable that not only the distribution but also the sizes of the (beta) Pic comets nuclei are found to be similar to sizes of the Solar system comets.

The size distribution of the nucleus of Solar system comets has been estimated for various locations. The size distribution of the Jupiter family comets and Oort cloud comets are found to be similar but not exactly the same.

For the Oort cloud, the size distribution has been estimated using a cometary activity model with a survey simulation and application to 150 long-period comets (LPC) detected over 7 years by the Pan-STARRS1 near-Earth object survey^{25}. For objects with diameters above 1 km, the distribution is found to be (gamma = 3.6pm 0.4) (stat.) (pm 0.7) (sys.), which is in remarkable agreement with our value for (beta) Pic exocomets. For smaller long-period comets, a shallower distribution is found with (gamma = 0.35) for diameters between 100 m and 1 km (Ref.^{25}). Note that in practice, in a similar manner as we have used the transit absorption depth as a proxy for the size estimate of the (beta) Pic comets, the above estimates for the Oort cloud comets have been obtained using the absolute H magnitudes of the nuclei as the proxy for their size.

For the (beta) Pic exocomets, the comparison may be more appropriate with the size distribution of the comets in the Jupiter family, which originates from the Kuiper Belt. In Ref.^{26} a catalog of absolute nuclear magnitudes of Jupiter family comets (JFC) has been used to derive a size distribution and to find (gamma =3.7 pm 0.3) for nuclei with radius between 2 and 5.5 kilometers (see Fig. 8 of Ref.^{26}). More recent works have provided a shallower distribution : analyzing a large number of optical observations, Ref.^{27} found a lower value with (gamma =2.9 pm 0.2) for nuclei with radius larger than 1.25 kilometers, and (gamma = 1.2) for smaller objects. This last result is consistent with the result described in Ref.^{28}, where images of Jupiter family comets obtained with the Hubble space telescope and the Keck telescopes have been analyzed. With a model fit to the observations, it is concluded that the intrinsic size distribution of comets in the Jupiter family is consistent with a (gamma = 3.5) power-law but truncated at small nucleus radii below 2.0 kilometers. In Ref.^{29} a similar distribution is obtained with (gamma = 2.9) for radius between 2 and 5 kilometers, interpreting the distribution shallower than the canonical Dohnanyi’s size distribution (gamma _} = 3.5) (Ref.^{24}) as due to fragmentation of the JFC objects. In Ref.^{30} the measured (gamma = 3.3 pm 0.2) for Jupiter family comets of radius between 2 and 10 kilometers is to be compared to the (gamma =2.0 pm 0.1) for long-period comets between 1 and 20 kilometers in radius.

Finally, in the Solar system the size distribution of extinct or dormant comets can be determined through the population of asteroids in comet orbits (ACO). Assuming that the exhaustion processes have no major impact on the size distribution, the observed distribution for these asteroids can be considered as a good proxy for the distribution for the parent comets. For these objects, Ref.^{31} found (gamma =3.55 pm 0.04) for the full sample with radius between 2.8 and 7 kilometers, (gamma =3.2 pm 0.04) for near Earth objects (NEO) with radius down to 1.4 kilometers and (gamma =3.45 pm 0.04) for non-near Earth objects (non-NEO) with radius between 2.8 and 7 kilometers.

Taken all together these estimates for the Solar system comets are in general agreement with the value that we obtained for the (beta) Pic comets, with some slightly shallower distributions in some cases (Fig. 2). This points toward the importance of collisional fragmentation in shaping the size distribution of the exocomets in the younger (beta) Pic planetary system.

The measured absorption depth distribution is the result of the distribution of several parameters for each individual comet, *e.g.*, the orbital parameters, cometary activity, composition, size, etc. Here, following the result of numerical simulations^{17,18}, we assumed that the absorption depth distribution is mainly dominated by the distribution of the exocomet intrinsic dust production rate and hence their size. In other terms, although other parameters play a role for each individual comet, their diversities are expected to have a lower impact on the observed transit absorption depth than the size. In support of that idea, in spectroscopy it is observed that the (beta) Pic transiting comets present similar orbital characteristics, which allows the classification in two different families^{5}. With similar orbits, different transiting exocomets have different dust tails mainly because of different dust production rate, and hence because of different size nuclei.

The observed distribution of exocomets in the young planetary system of (beta) Pic is strikingly similar to the distribution observed in the Solar system. This distribution seems to be ubiquitous and is also consistent with the canonical Dohnanyi’s size distribution^{24} ((gamma _} = 3.5)), which corresponds to the size distribution of a collisionally relaxed population (see discussion in Ref.^{32}). This indicates that the collisional process with fragmentation cascades is likely one of the dominant processes that shape the population of kilometer-sized bodies in the (beta) Pic planetary systems.

Histogram of the number of exocomet transit events as a function of the absorption depth. The red line shows the fit with a power law function with (alpha =2.3 pm 0.4). The uncertainty on the fitted parameters of the power law function have been evaluated using a Poisson distribution for the number of events in each bin of width (d AD=1.5cdot 10^{-4}). The red squares represents the number of expected events in each bin as calculated with the fitted power law.

Plot of the cumulative size distribution of the exocomets in (beta) Pic. The cumulative size distribution is plotted with blue squares for each exocomet and the corresponding fit excluding the largest comet is plotted with the red thick line. For comparison, published size distributions measured in the Solar system are plotted with thin dashed lines for asteroids in comets orbits (ACO), on near Earth orbits (NEO) and non-near Earth orbits (non-NEO) (A06, Ref.^{31}), Jupiter family comets (JFC) (T06, Ref.^{26}; S11, Ref.^{27}; F13, Ref.^{29}; B17, Ref.^{30}), and long-period comets (LPC) (B17, Ref.^{30}; B19, Ref.^{25}). In this plot, the size distributions for the objects in the Solar system have been scaled to have a cumulative number of about 10 objects with radius above 2 kilometers. The radii of the (beta) Pic comets have been estimated using the derivation described in the Method section. The conclusion on the similarity of the size distributions in (beta) Pic and the Solar system is independent of these estimates.

(beta) Pictoris has been observed by TESS at 2-minute cadence several times from October 2018 to February 2021. The available data-set covers a total of 156 days of observations in the optical domain, divided into 14 epochs of about 12 days each. The raw data shows a flux dispersion of about (10^{-3}), which is mainly due to (delta) Scuti pulsations in the stellar atmosphere.

The data from 19 October 2018 to 1 February 2019 have already been analysed^{20}. Three photometric events have been identified and attributed to the transits of three different exocomets, with one spectacular transit at Julian Day (JD) equals to (2457000+1486).

In addition to these pioneering observations, (beta) Pic has been observed from 20 November 2020 to 8 February 2021. For the observations of Sectors 4 to 7 and Sectors 31 to 34 (Extended Data Table 1), the TESS data products were obtained from the Mikulski Archive for Space Telescopes (MAST). We used the Presearch Data Conditioning (PDC) flux time series^{33,34,35} processed by the TESS Science Processing Operations Center (SPOC) pipeline version 5.0.23-20210212^{36,37}, with the spacecraft cosmic ray mitigation algorithm enabled.

(beta) Pictoris is prone to (delta) Scuti type photometric variations^{38,39,40}. These variations have a dominant frequency of (47.44,hbox {d}^{-1}) (corresponding to a period of about 0.5 hours), with an amplitude of up to (4times 10^{-3}) (Extended Data Fig. 1). These variations are superimposed to the photometric signatures of exocomet transits and must be corrected before searching for these exocomets.

Each of the 14 epochs of continuous observations (Extended Data Table 1) have been reduced separately; for each of them we extracted the set of frequencies and amplitudes of the pulsations using the Period04 software as previously done in Ref.^{20}. The Period04 software performs a Fourier transform on the data, giving the frequencies, amplitude, phase and signal to noise ratio of each harmonic in the time series. We conducted the frequency search between 15 and (100,hbox {day}^{-1}), to avoid removing slow variations that might be caused by a cometary transit. For each iteration of the software, the highest-amplitude frequency within the search range is selected, and added to a multi-sine model. This model is then optimized over amplitude and frequency of each harmonic, then removed from the original signal. For each iteration, the software computes the signal to noise ratio of the main remaining frequencies. The exit condition of the loop was chosen when the main frequency’s signal to noise ratio went under 4. We considered that below this limit, all that is left is noise.

The process has been applied for each of the 14 blocks of continuous observations, leading to an average number of 43 different pulsations for each block. After identifying all these pulsations, we subtracted them from the photometric measurements. The resulting residuals were then rebinned from an initial time sampling of 120 seconds to a time sampling of 1800 seconds.

To validate the result, we checked that the three exocomet transits already identified in Ref.^{20} are clearly visible in the data set cleaned from the (delta) Scuti variations by this procedure (Extended Data Fig. 2). The data are clean enough that new potential exocomet transits are suspected from the resulting light curve.

After removal of the (delta) Scuti variations, the resulting light curve still shows slow variations. These variations can have various origins (residuals of the (delta) Scuti variations that were not properly eliminated, systematic correlated noise, instrumental effects, etc.). Whatever the origin of these slow variations, astrophysical or instrumental, they need to be cleaned before searching for exocomet transits. To do so, for each epoch of continuous observations we searched for a smooth function modelling these variations. First we rebinned the time series at a 1-day-long interval, in order to “protect” the shorter dimmings (such as comets), excluded the 1-day measurements where the deepest exocometary transits are already clearly identified (by using the detection procedure described in the main text), then we interpolated a function with a cubic spline to model the generic form of the time series. Finally, we normalized our time series with this model. The result is a flattened light curve, which can be directly used to search for shallow dip events due to exocomet transits. The procedure is illustrated in Extended Data Fig. 3 where a potential second shallower exocomet transit closely follows another deep exocomet transit event.

The conclusion on the size distribution given in the main text does not require to estimate the true size of the comet’s nuclei. It only relies on the assumption that the absorption depth is a good proxy for the dust production rate, which is supported by numerical simulations of exocomets light curves^{18}, and that the production rate is proportional to the area of the comet’s nucleus, which is consistent with Solar system comets’ models and observations.

Nonetheless, using these simulations and Solar system observations, we can make a step further and derive the typical sizes of the comets detected in the (beta) Pic TESS light curve. From a newly calculated library of exocomet transit light curves similar to the one of Ref.^{18}, we derive a typical scaling law for the absorption depth *AD*, which is

where (dot{M_{1,}}}) is the dust production rate of the comet when it is at 1 au from the star, *q* is the orbital periastron distance and (M_*) is the mass of the star. Because the transiting comet is seen from a distance that is extremely large compared to the star-comet distance, for a given cometary tail size, thickness and impact parameter, the percentage of the stellar disk covered by the tail and therefore the absorption depth do not depend on the comet distance to the star. The above estimate of the absorption depth is valid over a wide range of impact parameter and wide range of longitude of perisastron of the comet’s orbit, (omega).

The periastron distance of the detected comets can be estimated using the transit time (Delta t). We have (Delta t = L_{chord}/v_}), where (L_{chord}) is the length of the transit chord of the planetary trajectory on the stellar disk that is crossed by the planet at the velocity (v_}). The chord length depends on the transit impact parameter : it is twice the planet radius for a central transit with a zero impact parameter and it is zero for a tangential transit with an impact parameter equals to the stellar radius. Here we used the mean chord length obtained assuming a random uniform distribution of the impact parameter from zero to one stellar radius, that is (overline{L_{chord}}=pi R_{*}/2). The transit velocity is given by (v_}=sqrt{GM_* /2q}(cos omega +1)sim sqrt{GM_*/q}). With a (beta) Pic radius of (R_*=1.7 R_{odot }) (Ref.^{41}) and a mass of (M_*=1.75 M_{odot }), we find (Delta t simeq 13 (sqrt{q/1,au})) hours. The best fits values of (Delta t) correspond to distances ranging from 0.03 to 1.3 au, in good agreement for the distances expected for the comet evaporation. The mean value of the estimated periastron distances is about 0.18 au. Using this mean distance, we obtain the following relationship for the observed absorption depth and the dust production rate :

Finally, the relation between the evaporation rate and the comet’s nucleus size can be derived by scaling the observation of the Hale-Bopp comet. Using a radius of about 30 kilometers^{23,42,43} and a dust production rate of (2cdot 10^6,hbox {kg},hbox {s}^{-1}) at 1 au (Ref.^{23}) for this well-observed dusty comet, we find (dot{M_{1,}}}simeq 2cdot 10^6,hbox {kg},hbox {s}^{-1}, (R/30,hbox {km})^2(L_*/L_{odot })). With a (beta) Pic luminosity of (8.7L_{odot }), we find

All together, we conclude that the radius of the (beta) Pic comets nuclei can be estimated using the photometric transit absorption depth with

Using this relationship, we derive a size of 1.5 km for the smallest detected comets ((ADsimeq 10^{-4})), and 6.7 km for the largest comet ((ADsimeq 20cdot 10^{-4})). These sizes are remarkably similar to the sizes of comets in the Solar system.

The observational data used in this work are publicly available in the Mikulski Archive for Space Telescope (MAST). The data in the tables and the final cleaned light curves are publicly available on GitHub at https://github.com/lecaveli/BetaPic_TESS.

Vidal-Madjar, A., Hobbs, L. M., Ferlet, R., Gry, C. & Albert, C. E. The circumstellar gas cloud around Beta Pictoris. II. *Astron. Astrophys.* **167**, 325–332 (1986).

ADS CAS Google Scholar

Ferlet, R., Hobbs, L. M. & Vidal-Madjar, A. The beta Pictoris circumstellar disk. V. Time variations of the Ca II-K line. *Astron. Astrophys.* **185**, 267–270 (1987).

ADS CAS Google Scholar

Beust, H., Lagrange-Henri, A. M., Vidal-Madjar, A. & Ferlet, R. The beta Pictoris circumstellar disk. X. Numerical simulations of infalling evaporating bodies. *Astron. Astrophys.* **236**, 202 (1990).

ADS CAS Google Scholar

Vidal-Madjar, A. *et al.* HST-GHRS observations of (beta) Pictoris: Additional evidence for infalling comets. *Astron. Astrophys.* **290**, 245–258 (1994).

ADS CAS Google Scholar

Kiefer, F. *et al.* Two families of exocomets in the (beta) Pictoris system. *Nature* **514**, 462–464. https://doi.org/10.1038/nature13849 (2014).

ADS CAS Article PubMed Google Scholar

Strøm, P. A. *et al.* Exocomets from a solar system perspective. *Publ. Astron. Soc. Pac.* **132**, 101001. https://doi.org/10.1088/1538-3873/aba6a0 (2020).

ADS Article Google Scholar

Miret-Roig, N. *et al.* Dynamical traceback age of the (beta) Pictoris moving group. *Astron. Astrophys.* **642**, A179. https://doi.org/10.1051/0004-6361/202038765 (2020).

Article Google Scholar

Smith, B. A. & Terrile, R. J. A circumstellar disk around (beta) Pictoris. *Science* **226**, 1421–1424. https://doi.org/10.1126/science.226.4681.1421 (1984).

ADS CAS Article PubMed Google Scholar

Kalas, P., Larwood, J., Smith, B. A. & Schultz, A. Rings in the Planetesimal Disk of (beta) Pictoris. *Astrophys. J. Lett.* **530**, L133–L137. https://doi.org/10.1086/312494 (2000).

ADS CAS Article Google Scholar

Roberge, A., Feldman, P. D., Weinberger, A. J., Deleuil, M. & Bouret, J.-C. Stabilization of the disk around (beta)Pictoris by extremely carbon-rich gas. *Nature* **441**, 724–726. https://doi.org/10.1038/nature04832 (2006).

ADS CAS Article PubMed Google Scholar

Apai, D. *et al.* The inner disk structure, disk-planet interactions, and temporal evolution in the (beta) Pictoris system: A two-epoch HST/STIS coronagraphic study. *Astrophys. J.* **800**, 136. https://doi.org/10.1088/0004-637X/800/2/136 (2015).

ADS CAS Article Google Scholar

Brandeker, A. *et al.* Herschel detects oxygen in the (beta) Pictoris debris disk. *Astron. Astrophys.* **591**, A27. https://doi.org/10.1051/0004-6361/201628395 (2016).

CAS Article Google Scholar

Lagrange, A. M. *et al.* A giant planet imaged in the disk of the young star (beta) Pictoris. *Science* **329**, 57. https://doi.org/10.1126/science.1187187 (2010).

ADS CAS Article PubMed Google Scholar

Snellen, I. A. G. & Brown, A. G. A. The mass of the young planet Beta Pictoris b through the astrometric motion of its host star. *Nat. Astron.* **2**, 883–886. https://doi.org/10.1038/s41550-018-0561-6 (2018).

ADS Article Google Scholar

Lagrange, A. M. *et al.* Evidence for an additional planet in the (beta) Pictoris system. *Nat. Astron.* **3**, 1135–1142. https://doi.org/10.1038/s41550-019-0857-1 (2019).

ADS Article Google Scholar

Lacour, S. *et al.* The mass of Beta Pictoris c from Beta Pictoris b orbital motion. *arXiv e-prints*arXiv:2109.10671 (2021).

Lecavelier des Etangs, A., Vidal-Madjar, A. & Ferlet, R. Photometric stellar variation due to extra-solar comets. *Astron. Astrophys.* **343**, 916–922 (1999).

ADS Google Scholar

Lecavelier des Etangs, A. A library of stellar light variations due to extra-solar comets. *Astron. Astrophys. Suppl. Ser.* **140**, 15–20. https://doi.org/10.1051/aas:1999114 (1999).

ADS Article Google Scholar

Rappaport, S. *et al.* Likely transiting exocomets detected by Kepler. *Mon. Not. R. Astron. Soc.* **474**, 1453–1468. https://doi.org/10.1093/mnras/stx2735 (2018).

ADS CAS Article PubMed Google Scholar

Zieba, S., Zwintz, K., Kenworthy, M. A. & Kennedy, G. M. Transiting exocomets detected in broadband light by TESS in the (beta) Pictoris system. *Astron. Astrophys.* **625**, L13. https://doi.org/10.1051/0004-6361/201935552 (2019).

ADS CAS Article Google Scholar

Ricker, G. R. *et al.* Transiting Exoplanet Survey Satellite (TESS). *J. Astron. Telesc. Instrum. Syst.* **1**, 014003. https://doi.org/10.1117/1.JATIS.1.1.014003 (2015).

ADS Article Google Scholar

Lecavelier des Etangs, A. *et al.* Observation of the central part of the beta Pictoris disk with an anti-blooming CCD. *Astron. Astrophys.* **274**, 877–882 (1993).

ADS Google Scholar

Jewitt, D. & Matthews, H. Particulate mass loss from comet hale-bopp. *Astron. J.* **117**, 1056–1062. https://doi.org/10.1086/300743 (1999).

ADS Article Google Scholar

Dohnanyi, J. S. Collisional model of asteroids and their debris. *J. Geophys. Res.* **74**, 2531–2554. https://doi.org/10.1029/JB074i010p02531 (1969).

ADS Article Google Scholar

Boe, B. *et al.* The orbit and size-frequency distribution of long period comets observed by Pan-STARRS1. *Icarus* **333**, 252–272. https://doi.org/10.1016/j.icarus.2019.05.034 (2019).

ADS Article Google Scholar

Tancredi, G., Fernández, J. A., Rickman, H. & Licandro, J. Nuclear magnitudes and the size distribution of Jupiter family comets. *Icarus* **182**, 527–549. https://doi.org/10.1016/j.icarus.2006.01.007 (2006).

ADS Article Google Scholar

Snodgrass, C., Fitzsimmons, A., Lowry, S. C. & Weissman, P. The size distribution of Jupiter Family comet nuclei. *Mon. Not. R. Astron. Soc.* **414**, 458–469. https://doi.org/10.1111/j.1365-2966.2011.18406.x (2011).

ADS Article Google Scholar

Meech, K. J., Hainaut, O. R. & Marsden, B. G. Comet nucleus size distributions from HST and Keck telescopes. *Icarus* **170**, 463–491. https://doi.org/10.1016/j.icarus.2004.03.014 (2004).

ADS Article Google Scholar

Fernández, Y. R. *et al.* Thermal properties, sizes, and size distribution of Jupiter-family cometary nuclei. *Icarus* **226**, 1138–1170. https://doi.org/10.1016/j.icarus.2013.07.021 (2013).

ADS Article Google Scholar

Bauer, J. M. *et al.* Debiasing the NEOWISE cryogenic mission comet populations. *Astron. J.* **154**, 53. https://doi.org/10.3847/1538-3881/aa72df (2017).

ADS Article Google Scholar

Alvarez-Candal, A. & Licandro, J. The size distribution of asteroids in cometary orbits and related populations. *Astron. Astrophys.* **458**, 1007–1011. https://doi.org/10.1051/0004-6361:20064971 (2006).

ADS Article Google Scholar

O’Brien, D. P. & Greenberg, R. The collisional and dynamical evolution of the main-belt and NEA size distributions. *Icarus* **178**, 179–212. https://doi.org/10.1016/j.icarus.2005.04.001 (2005).

ADS Article Google Scholar

Smith, J. C. *et al.* Kepler presearch data conditioning II: A Bayesian approach to systematic error correction. *Publ. Astron. Soc. Pac.* **124**, 1000. https://doi.org/10.1086/667697 (2012).

ADS Article Google Scholar

Stumpe, M. C. *et al.* Kepler presearch data conditioning I-Architecture and algorithms for error correction in Kepler light curves. *Publ. Astron. Soc. Pac.* **124**, 985. https://doi.org/10.1086/667698 (2012).

ADS Article Google Scholar

Stumpe, M. C. *et al.* Multiscale systematic error correction via wavelet-based bandsplitting in Kepler data. *Publ. Astron. Soc. Pac.* **126**, 100. https://doi.org/10.1086/674989 (2014).

ADS Article Google Scholar

Jenkins, J. M. *et al.* The TESS science processing operations center. In *Software and Cyberinfrastructure for Astronomy IV*, Vol. 9913 of *Proceedings of the SPIE*, 99133E. https://doi.org/10.1117/12.2233418 (2016).

Caldwell, D. A. *et al.* TESS science processing operations center FFI target list products. *Res. Notes Am. Astron. Soc.* **4**, 201. https://doi.org/10.3847/2515-5172/abc9b3 (2020).

ADS Article Google Scholar

Koen, C. (delta) Scuti pulsations in (beta) Pictoris. *Mon. Not. R. Astron. Soc.* **341**, 1385–1387. https://doi.org/10.1046/j.1365-8711.2003.06509.x (2003).

ADS Article Google Scholar

Mékarnia, D. *et al.* The (delta) Scuti pulsations of (beta) Pictoris as observed by ASTEP from Antarctica. *Astron. Astrophys.* **608**, L6. https://doi.org/10.1051/0004-6361/201732121 (2017).

ADS CAS Article Google Scholar

Zwintz, K. *et al.* Revisiting the pulsational characteristics of the exoplanet host star (beta) Pictoris. *Astron. Astrophys.* **627**, A28. https://doi.org/10.1051/0004-6361/201834744 (2019).

CAS Article Google Scholar

Kervella, P. *et al.* VINCI/VLTI Observations of Main Sequence Stars. In *Stars as Suns : Activity, Evolution and Planets*, (Eds. Dupree, A. K. & Benz, A. O. ) Vol. 219, 80 (2004). arXiv:astro-ph/0309784.

Fernández, Y. R. *et al.* The inner coma and nucleus of comet Hale–Bopp: Results from a stellar occultation. *Icarus* **140**, 205–220. https://doi.org/10.1006/icar.1999.6127 (1999).

ADS Article Google Scholar

Bair, A. N., Schleicher, D. G. & Farnham, T. The Extremely Active Comet C/Hale-Bopp (1995 O1): Production Rates from Nearly Five Years of Narrowband Photometry. In *AAS/Division for Planetary Sciences Meeting Abstracts #50*, Vol. 50 of *AAS/Division for Planetary Sciences Meeting Abstracts* 210.06 (2018).

Download references

We thank Paul Wiegert for enlightening discussion on the size distribution of comets in the Solar system. This paper includes data collected with the TESS mission, obtained from the MAST data archive at the Space Telescope Science Institute (STScI). Funding for the TESS mission is provided by the NASA Explorer Program. STScI is operated by the Association of Universities for Research in Astronomy, Inc., under NASA contract NAS 5-26555. We thank the TESS Team members for making available the extremely accurate photometric data they obtained. ALdE, LC, GH, EM, FK, and AVM acknowledge support from the CNES (Centre national d’études spatiales, France).

Institut d’astrophysique de Paris, CNRS, UMR 7095, Sorbonne Université, 98 bis bd Arago, 75014, Paris, France

Alain Lecavelier des Etangs, Lucie Cros, Guillaume Hébrard, Eder Martioli, Flavien Kiefer & Alfred Vidal-Madjar

Ecole nationale supérieure des mines de Paris, Université PSL, 60 boulevard Saint-Michel, 75272, Paris, France

Lucie Cros

Observatoire de Haute-Provence, 04870, St Michel l’Observatoire, France

Guillaume Hébrard

Laboratório Nacional de Astrofísica, Rua Estados Unidos 154, 37504-364, Itajubá, MG, Brazil

Eder Martioli

LESIA, Observatoire de Paris, Université PSL, CNRS, Sorbonne Université, Univ. Paris Diderot, Sorbonne Paris Cité, 5 place Jules Janssen, 92195, Meudon, France

Marc Duquesnoy, Flavien Kiefer, Sylvestre Lacour & Anne-Marie Lagrange

Leiden Observatory, Leiden University, Leiden, The Netherlands

Matthew A. Kenworthy

Univ. Grenoble Alpes, CNRS, IPAG, 38000, Grenoble, France

Nadège Meunier

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

A.L.d.E. initiated and directed the project, analysed the cleaned light curves, interpreted the result and wrote the bulk of the manuscript. L.C. analysed the raw TESS data, developed and applied the algorithms, produced the cleaned light curves, extracted the first list of exocomet transits, participated to the interpretation and wrote parts of the manuscript. G.H. contributed to the project design, definition of the scientific goal, the derivation of Fig. 2. and to the writing. E.M. contributed to the TESS data extraction and interpretation, and to the writing. M.D., A.M.L. and N.M. developed an independent analysis of the data for a cross-check of the result, extracted a confirmation list of exocomet transits, and participated to the writing. M.K., F.K., S.L. and A.V.M. participated to the definition of the project, scientific interpretation of the results, clarification of the presentation and to the writing. All authors discussed the scientific result, edited the manuscript and contributed to the final version.

Correspondence to Alain Lecavelier des Etangs.

The authors declare no competing interests.

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and Permissions

Lecavelier des Etangs, A., Cros, L., Hébrard, G. *et al.* Exocomets size distribution in the (beta) Pictoris planetary system. *Sci Rep* **12, **5855 (2022). https://ift.tt/Jr6R2lo

Download citation

Received: 06 January 2022

Accepted: 16 March 2022

Published: 28 April 2022

DOI: https://ift.tt/Jr6R2lo

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.

Advertisement

Advanced search

Scientific Reports (*Sci Rep*) ISSN 2045-2322 (online)

© 2022 Springer Nature Limited

Sign up for the *Nature Briefing* newsletter — what matters in science, free to your inbox daily.

Tech

via Inferse.com https://www.inferse.com

May 12, 2022 at 09:00PM