Free Access
Issue
A&A
Volume 648, April 2021
Article Number L1
Number of page(s) 6
Section Letters to the Editor
DOI https://doi.org/10.1051/0004-6361/202140464
Published online 07 April 2021

© ESO 2021

1. Introduction

HD 149026b is considered a hot Jupiter with a mass of M = 117M (Sato et al. 2005). This planet likely has a large metal (or heavy element) mass of > 40 M (Dodson-Robinson & Bodenheimer 2009; Carter et al. 2009). Recently, the TESS mission discovered a planet, TOI-849b, with M = 39 M, which is mostly composed of heavy elements (Armstrong et al. 2020). The heavy element mass in giant planets has also been estimated for other planets, whose masses and radii are both measured. Thorngren et al. (2016) derived the metal mass of giant planets (mostly warm Jupiters) and found that the typical metal mass range is approximately 10−100 M; they also found that the metal-mass fraction is typically higher than 0.1.

Although giant planets with large heavy-element masses (> 30 M) appear to be relatively common in observations, the accretion of large heavy-element masses is difficult when planetary cores are being assembled. When cores grow by accreting surrounding planetesimals, their growth ceases at the planetesimal isolation mass (= 2πaΔaΣd). Assuming Δa = 10 RH, where RH is the Hill radius, an extremely high solid surface density of Σd ≃ 1000 g cm−2 is required for the formation of cores larger than 30 M at a semimajor axis of a = 0.3 au. Even if pebble accretion is considered, the pebble isolation mass is usually smaller than 25 M (e.g., Lambrechts et al. 2014). It is possible that giant planets that formed in wide orbits (a > a few tens of au) accrete planetesimals of several tens of Earth masses during migration (Shibata et al. 2020); however, the formation of giant planets in such wide orbits is challenging (Bitsch et al. 2019; Johansen & Bitsch 2019).

Additionally, even if massive planetary cores are available in protoplanetary disks, subsequent runaway gas accretion typically does not yield giant planets with metal fractions higher than 10%. Giant planets are formed by the accretion of the gaseous components of the protoplanetary disk. Since the accretion process proceeds in a runaway fashion (e.g., Pollack et al. 1996), giant planets accrete massive envelopes, in which case the metal-mass fraction is low. For the formation of giant planets with high metal-mass fractions, the photoevaporation of accreted envelopes may play a role. The origin of HD 149026b and TOI-849b, however, is unlikely to be explained by the photoevaporation process because massive cores of ≳30 M should suppress the envelope loss by the energy-limited hydrodynamic escape (e.g., Kurokawa & Nakamoto 2014). The interior modeling suggests that HD 149026b has an envelope of ≳30% of its total mass (Fortney et al. 2006; Ikoma et al. 2006), whereas the envelope fraction of TOI-849b is no more than 3.9wt%. If TOI-849b had a larger radius (about 1.5 times the Jupiter radius), the Roche-lobe overflow driven by tidal disruption would be a possible mechanism for the origin of the large metal fraction (e.g., Armstrong et al. 2020). As an alternative way, the envelope mass loss during giant impacts can play an important role. For example, Ikoma et al. (2006) claim that a substantial number (∼90%) of envelopes can be lost by a high-speed (∼2.5 times the escape velocity) collision between giant planets, which is considered a possible mechanism for the formation of HD 149026b. For the origin of TOI-849b, Armstrong et al. (2020) also suggest that the envelope can be removed during the giant planet collision, enabling the formation of planets with large metal-mass fractions.

Meanwhile, Ogihara et al. (2020) investigated the formation of super-Earths and sub-Neptunes (SENs) and demonstrated that SENs with low-mass atmospheres were formed in a rapidly dissipating disk by photoevaporation. We expect that such a disk would also help explain the origin of giant planets with high metal fractions.

2. Model description

We investigate the planet formation using a similar approach to that in Ogihara et al. (2020), to which we refer the readers for full details. Several salient points are briefly described here.

The key ingredient of our model is the detailed disk evolution. Suzuki et al. (2016) developed a 1D disk-evolution model that considers mass loss due to magnetohydrodynamic (MHD) disk winds (Suzuki & Inutsuka 2009) and wind-driven accretion (Blandford & Payne 1982), in addition to viscous evolution. Subsequently, Kunitomo et al. (2020) improved this model by considering the mass loss due to photoevaporative winds (PEWs) (e.g., Alexander et al. 2006). The temporal evolution of the gas-surface density is shown in Fig. 1 of Ogihara et al. (2020). We observe that the surface density profile is flattened in the close-in region inside r ∼ 1 au, where r is the radial distance from the star due to the effects of the disk winds. The bottom panels of Fig. 1 in this Letter show the disk accretion rate of r = 0.3 au for models with and without PEWs. The disk-accretion rate is as high as 100 M Myr−1 even at t = 10 Myr for models with no PEWs, while for models with PEWs, the disk-accretion rate rapidly decreases during the disk-dissipation phase (t > 1 Myr).

thumbnail Fig. 1.

Evolution of the semimajor axis a, metal mass Mz, and envelope mass fraction Menv/M for our fiducial model in disks with and without PEW effects. The dashed line in panel a indicates that the growth of the core ceases at the pebble isolation mass if mutual collisions are ignored. Bottom panels: gas accretion rate at r = 0.3 au. (a) Without PEW, Fpb, fid, (b) with PEW, Fpb, fid.

Another major physical ingredient is the envelope mass loss during giant impacts. Stewart et al. (2014) derived a scaling law for the mass loss during giant impacts (see Sect. 2.1.6 of Ogihara & Hori 2020); we utilized an updated model of Kegerreis et al. (2020). Although the mass loss was estimated by the ground speed in Stewart et al. (2014), the envelope was directly included in the simulations of Kegerreis et al. (2020). The envelope loss fraction, X1, from the planet with Mcore, 1 in a collision with a planet with Mcore, 2 is given by

X 1 = 0.64 [ ( v imp v esc ) 2 ( M core , 1 M tot ) 1 / 2 f M ] 0.65 , $$ \begin{aligned} X_1 = 0.64 \left[\left(\frac{{ v}_{\rm imp}}{{ v}_{\rm esc}}\right)^{2} \left(\frac{M_{\rm core,1}}{M_{\rm tot}}\right)^{1/2} f_{\rm M}\right]^{0.65}, \end{aligned} $$(1)

where vimp, vesc, and Mtot are the impact velocity, two-body escape velocity, and total core mass (=Mcore, 1 + Mcore, 2), respectively. The fractional interacting mass, fM, is given by Eq. (B.1) of Kegerreis et al. (2020). Here, we assume that the bulk core densities of colliding planets are equal. The envelope loss from the planet with Mcore, 2 was calculated using the same method.

Here, we summarize other conditions. The simulations were initiated with 40 bare embryos with M = 0.01 M, which were placed between a = 0.2 − 1 au with orbital separations of 15 mutual Hill radii. Cores can grow via pebble accretion (e.g., Ormel & Klahr 2010; Lambrechts & Johansen 2012) until they reach the pebble isolation mass (e.g., Bitsch et al. 2018). As mentioned in Sect. 1, the pebble isolation mass is smaller than 25 M; thus, it is difficult to explain the massive metal mass of the observed giant planets. The pebble flux is assumed to be 3 × 10−4exp(−t/Myr) M yr−1, which is larger than that used in the simulations of super-Earth formation (Ogihara et al. 2020), to form more massive planets. The pebble flux was reduced in accordance with the number of pebbles accreted on outer planets.

Cores can start accreting massive envelopes when they reach the critical core mass (Ogihara & Hori 2020):

M crit = 13 M ( M ˙ acc , pb 10 6 M yr 1 ) 0.23 , $$ \begin{aligned} M_{\rm crit} = 13\,M_\oplus \left(\frac{\dot{M}_{\rm acc,pb}}{10^{-6}\,M_\oplus \,\mathrm{yr^{-1}}}\right)^{0.23}, \end{aligned} $$(2)

where acc,pb denotes the pebble accretion rate. While pebble accretion results in core growth, it can also heat the envelope, leading to a significant increase in the critical core mass. When the pebble accretion ceases after the cores reach the pebble isolation mass, the critical core mass decreases to a value lower than 10 M (Hori & Ikoma 2010). The actual envelope accretion rate was calculated by the Kelvin-Helmholtz contraction and regulated by the disk-accretion rate, disk.

Orbits of planets with masses larger than that of Mars change due to the gravitational interaction with the disk gas. We consider the damping of eccentricity and inclination (e.g., Cresswell & Nelson 2008), as well as orbital migration (e.g., Paardekooper et al. 2011; Kanagawa et al. 2018).

3. Results

We present the typical simulation results. Figure 1 shows the temporal evolution of the semimajor axis, the metal mass, and the envelope mass fraction for simulations with our fiducial pebble flux. Panel a indicates the results of simulations for a disk without PEWs, while panel b corresponds to those of simulations for a disk with PEWs. As mentioned in Sect. 2, the bottom panels show the disk-accretion rate at r = 0.3 au. First, we examined the results of simulations in a disk with no PEWs (panel a). As shown in the evolution of the semimajor axis, planets did not undergo significant migration1. This is because the gas-surface density profile was altered from simple power-law distributions due to the effects of mass loss and mass accretion caused by MHD disk winds. Further, we find that planets exhibit orbital instability (hence the giant impact events) during the disk dissipation (t ≳ 1 Myr). After giant impact events and the ejection of some planets, including the two innermost planets, a few giant planets eventually formed.

The second row of Fig. 1 highlights the planets that can have a large mass of heavy elements (> 10 M). The gray dashed curve indicates the metal mass evolution of planets at a = 0.3 au for a reference simulation in which the mutual interactions and collisions between planets were ignored. This shows that the growth of cores is terminated when they reach the pebble isolation mass (≃10 − 20 M). Conversely, in the results of N-body simulations that consider mutual interactions, planets undergo giant impacts after they reach the pebble isolation mass. Therefore, here, we demonstrate that the high metal mass of giant planets can be achieved by giant impacts. Similar results promoting the formation of massive cores have also been demonstrated in N-body simulations of cold-Jupiter formation (Wimarsson et al. 2020). A related analytical argument has also been reported by Ginzburg & Chiang (2020) recently.

Regarding the envelope mass fraction, planets acquire massive envelopes (> 80wt%) in a disk with no PEW effects. While planets grow via pebble accretion, the critical core mass can increase to a value larger than 10 M due to the heating of the envelope by the accreting pebbles (Eq. (2)). After planets reach the pebble isolation mass, the critical core mass decreases to ≃3 M or smaller (Hori & Ikoma 2010). At t ≃ 1 Myr, the cores exceed the critical core mass and start accreting massive envelopes. The gas-accretion rate in the disk is sufficiently high for envelope accretion (> 100 M Myr−1) until t ≃ 10 Myr (see bottom panel), leading to the accretion of massive envelopes. Moreover, we find that even if a portion of the envelope is lost during giant impacts, the envelope mass fraction is high at the end. This is because planets can reaccrete envelopes after envelope erosion during impacts, which is discussed in Sect. 4.

Here, we examine the results of simulations in a disk with PEW effects (panel b). The evolutions of the semimajor axis and metal mass are similar to those in panel a. In other words, planets do not undergo significant migration and have a high mass of heavy elements after giant impacts. In the evolution of the envelope mass fraction, we observe a notable difference. Planets can avoid accreting massive envelopes, and the final envelope mass fraction is less than 0.5. This is because the gas disk is rapidly cleared out during the disk-dissipation phase (t ≳ 1 Myr), and the envelope accretion onto planets is quenched. Consequently, we find that both the high metal mass and the high metal fraction in the observed giant planets can be explained in a rapidly dissipating disk with PEWs2.

Figure 2 compares the estimates of the metal mass of the observed giant planets with our simulation results. The gray circles represent the estimates derived by Thorngren et al. (2016). Both the metal mass and the metal-mass fraction can be fitted by simple power laws. HD 149026b and TOI-849b are also plotted with purple and brown symbols, respectively. Panels a and b indicate the results of simulations in disks without and with PEW effects, respectively. The pebble flux was varied from the fiducial flux of Fpb, fid = 3 × 10−4exp(−t/Myr) M yr−1, in which the results of simulations with the same pebble flux are indicated with the same color. The values of the pebble flux were chosen such that the final masses of the planets are between ≃20 and 1000 M. For the cases without PEWs, the pebble flux was varied between 2/3 and 4/3 Fpb, fid (the time-integrated mass flux is between 200 and 400 M). For the cases with PEWs, the flux was varied between 1 and 3 Fpb, fid (the time-integrated mass flux is between about 230 and 690 M). Although the total mass flux is high, some observations suggest a high solid mass reservoir in protoplanetary disks (e.g., Carrasco-González et al. 2016). Five runs with different initial locations of embryos were performed for each pebble flux. Regardless of whether the PEW effects are considered or not, it is shown that planets with a large metal mass of more than 30 M can form. This is due to the promoted growth by giant impacts between the planets, which stop their growth once the pebble isolation mass is achieved.

thumbnail Fig. 2.

Comparison of the metal mass Mz and metal-mass fraction Mz/M between the observed transiting giant planets and the results of our simulations. Gray circles indicate the estimates by Thorngren et al. (2016), which were fitted by power-law distributions (dashed lines). HD 149026b and TOI-849b are plotted with purple and brown symbols, respectively. The simulation results are plotted with cresses and circles. The pebble flux was varied, and five runs were performed for each pebble flux. (a) Without PEW, (b) with PEW.

Regarding the metal-mass fraction, the fraction tends to be lower than 0.1 after the accretion of massive envelops in a disk with no PEW effects, which is generally lower than the estimates of Thorngren et al. (2016). For the simulations in a disk with PEW effects, on the other hand, the metal-mass fraction tends to be higher than 0.1 because the envelope accretion is terminated by the rapid disk clearing. The metal-mass fraction is consistent with the estimates for observed giant planets. The trend that massive planets tend to have lower metal-mass fractions is interpreted as follows. For relatively large pebble fluxes, cores reach the pebble isolation mass quite early and start the runaway gas accretion. In addition, large cores can form after giant impacts of large pebble fluxes.

Furthermore, we find that the origins of the distinctive planets (i.e., HD 149026b and TOI-849b) can be explained by our simulations. In the simulations with no PEW effects, no planets that are similar to the distinctive planets form. However, we find that, in a disk with rapid disk clearing by PEWs, planets with a high metal mass and metal fraction similar to HD 149026b and TOI-849b can be formed. Although the origins of the large metal fraction of these two planets are not well understood (see Sect. 1), our model can explain their origins. We note that the formation of hot Jupiters (a < 0.1 au) was not directly treated in our simulations because of the extreme computational cost required for N-body simulations of close-in planets. According to previous simulations that investigate the formation of closer-in planets (Ogihara & Hori 2020), it is most likely that orbital locations of hot Jupiters can also be explained by our model.

More interestingly, the relationship between the eccentricity and the metal mass suggested by Thorngren et al. (2016) also seems to be better explained with PEWs. Figure 13 of Thorngren et al. (2016) suggests that planets with a high eccentricity tend to have a high heavy-element mass, although this is not statistically significant because the number of planets is not high. We do not go into detail here; however, we find that this trend can be better explained in the case with PEWs.

4. Effects of envelope erosion during giant impacts

For the origin of HD 149026b and TOI-849b, it has been suggested that the impact erosion of accreted envelopes during a high-speed impact ( > 2.5 vesc) between planets plays an important role in increasing the metal-mass fraction. In the present study, however, we find that the impact erosion of envelopes alone does not explain the formation of giant planets with low envelope mass fractions. We instead show that the cutoff in disk accretion by PEWs plays a crucial role. We discuss this point thoroughly.

Figure 3 shows the envelope mass loss fraction of each collision for all simulations of the fiducial flux in disks with and without PEW effects. It is shown that although the typical impact velocity is less than 2 vesc, there exist some high-speed impacts with > 2 vesc. In some high-speed impacts, more than 90% of the envelope can be lost. Notably, both the impact speed and the envelope loss fraction in the simulations in a disk without PEW effects are not lower than those in the simulations in a disk with PEW effects. Nevertheless, there is a large difference in the final envelope mass fraction between the simulations, as shown in Fig. 2. This is because some accretions exist in the disk after the giant impacts (t = 1 − 10 Myr) (bottom panel of Fig. 1) in the disk with no PEW effects, and planets can reaccrete envelopes after giant impacts. Consequently, the final envelope mass fraction increases. Therefore, we can conclude that the rapid disk clearing after giant impacts plays a crucial role in the formation of giant planets with low envelope mass fractions. We note that planets with low envelope mass fractions may form in the simulations with no PEW effects if further giant impacts occur after the disk gas is depleted (t > 10 Myr). However, this is highly unlikely because the orbital separation between the planets is already large after the “first” giant impact phase (t = 1 − 10 Myr). The system shown in Fig. 1a is Hill stable (Gladman 1993), and no further impact events would occur.

thumbnail Fig. 3.

Envelope mass-loss fraction for all the giant impacts in the simulations with the fiducial pebble flux. Circles with blue edges represent collisions in simulations without PEW effects, while those with orange edges are in simulations with PEW effects. Each size is proportional to mass, and the color corresponds to its envelope fraction before impacts. A summary of all mass loss events is presented in Fig. A.1.

Next, we performed additional simulations to confirm the effect of impact erosion. Figure 4 shows a summary of the simulations in which the mass loss during impacts was ignored. The other conditions were maintained as in the simulations shown in Fig. 2b. By comparing Fig. 4 with Fig. 2b, we find that the final metal-mass fraction is almost the same. In other words, even if the impact erosion is not considered, giant planets with high metal-mass fractions can form in a disk with PEWs. In contrast, by comparing the metal mass fraction in Fig. 4 with that in Fig. 2a, it is clear that the metal-mass fraction in Fig. 4 is significantly larger. Therefore, this result confirmed that the rapid disk clearing by PEWs was the main reason for the high metal-mass fraction.

thumbnail Fig. 4.

Results of additional simulations in which the envelope loss during giant impacts is ignored.

Here, we comment on the model for the mass loss during impacts. First, we do not consider the loss of metal materials during collisions. However, solid cores cannot merge completely when high-speed collisions occur (e.g., Denman et al. 2020). Even if such effects are considered, the planetary envelopes can be lost significantly, and the metal-mass fractions are expected to be large. Thus, our conclusions are likely to remain unchanged. Second, it is not clear if the model used in this study is a good approximation because the model has not been confirmed for collisions between planets with high core masses (> 10 M) and for hit-and-run collisions yet. Although these would not change our conclusions, detailed studies on the impact-induced envelope escape are still important.

5. Conclusions

We demonstrate that the runaway gas accretion triggered by the formation of massive cores through giant impacts can be quenched by a rapid disk clearing due to PEWs, resulting in the formation of giant planets with high metal fractions. Further, we find that the metal mass in giant planets can increase after giant impacts between planets with the pebble isolation mass. The results reveal that giant planet collisions alone cannot explain the observed properties of transiting gas planets. The suppression of envelope accretion by the rapid-disk dissipation plays an important role in the increase in the metal-mass fraction of the giant planets. As a result, the high metal mass and metal-mass fraction of the observed giant planets, including HD 149026b and TOI-849b, can be naturally explained without relying on the impact erosion after the substantial depletion of the disk gas (t > 10 Myr) and other envelope-loss mechanisms.


1

According to the suppression of type-I migration (e.g., Ogihara et al. 2018, 2020), it is possible to prevent super-Earth cores from having an icy composition, which is consistent with the inference that the majority of super-Earths and sub-Neptunes have rocky cores (e.g., Agol et al. 2021; Rogers & Owen 2021).

2

There is a slight difference in the number of remaining planets between Figs. 1a and b. This difference may be due to the presence of runaway gas accretion. In all of our simulations, however, there is no clear trend that the number of remaining planets depends on whether or not PEWs are included.

Acknowledgments

We thank the anonymous referee whose comments improved this manuscript. Numerical computations were in part carried out on PC cluster at Center for Computational Astrophysics of the National Astronomical Observatory of Japan. M.O. is supported by JSPS KAKENHI Grant Numbers 18K13608 and 19H05087; Y.H. by 18H05439; K.K. by 20J01258. This study was inspired by discussions especially with Shigeru Ida and Daniel Thorngren at Ringberg Castle Workshop held in September 2019.

References

  1. Agol, E., Dorn, C., Grimm, S. L., et al. 2021, Planet. Sci. J., 2, 1 [Google Scholar]
  2. Alexander, R. D., Clarke, C. J., & Pringle, J. E. 2006, MNRAS, 369, 216 [Google Scholar]
  3. Armstrong, D. J., Lopez, T. A., Adibekyan, V., et al. 2020, Nature, 583, 39 [CrossRef] [Google Scholar]
  4. Bitsch, B., Morbidelli, A., Johansen, A., et al. 2018, A&A, 612, A30 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  5. Bitsch, B., Izidoro, A., Johansen, A., et al. 2019, A&A, 623, A88 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  6. Blandford, R. D., & Payne, D. G. 1982, MNRAS, 199, 883 [NASA ADS] [CrossRef] [Google Scholar]
  7. Carrasco-González, C., Henning, T., Chandler, C. J., et al. 2016, ApJ, 821, L16 [NASA ADS] [CrossRef] [Google Scholar]
  8. Carter, J. A., Winn, J. N., Gilliland, R., & Holman, M. J. 2009, ApJ, 696, 241 [Google Scholar]
  9. Cresswell, P., & Nelson, R. P. 2008, A&A, 482, 677 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  10. Denman, T. R., Leinhardt, Z. M., Carter, P. J., & Mordasini, C. 2020, MNRAS, 496, 1166 [Google Scholar]
  11. Dodson-Robinson, S. E., & Bodenheimer, P. 2009, ApJ, 695, L159 [Google Scholar]
  12. Fortney, J. J., Saumon, D., Marley, M. S., Lodders, K., & Freedman, R. S. 2006, ApJ, 642, 495 [Google Scholar]
  13. Ginzburg, S., & Chiang, E. 2020, MNRAS, 498, 680 [Google Scholar]
  14. Gladman, B. 1993, Icarus, 106, 247 [Google Scholar]
  15. Hori, Y., & Ikoma, M. 2010, ApJ, 714, 1343 [NASA ADS] [CrossRef] [Google Scholar]
  16. Ikoma, M., Guillot, T., Genda, H., Tanigawa, T., & Ida, S. 2006, ApJ, 650, 1150 [Google Scholar]
  17. Johansen, A., & Bitsch, B. 2019, A&A, 631, A70 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  18. Kanagawa, K. D., Tanaka, H., & Szuszkiewicz, E. 2018, ApJ, 861, 140 [NASA ADS] [CrossRef] [Google Scholar]
  19. Kegerreis, J. A., Eke, V. R., Catling, D. C., et al. 2020, ApJ, 901, L31 [Google Scholar]
  20. Kunitomo, M., Suzuki, T. K., & Inutsuka, S.-I. 2020, MNRAS, 492, 3849 [Google Scholar]
  21. Kurokawa, H., & Nakamoto, T. 2014, ApJ, 783, 54 [Google Scholar]
  22. Lambrechts, M., & Johansen, A. 2012, A&A, 544, A32 [CrossRef] [EDP Sciences] [Google Scholar]
  23. Lambrechts, M., Johansen, A., & Morbidelli, A. 2014, A&A, 572, A35 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  24. Ogihara, M., & Hori, Y. 2020, ApJ, 892, 124 [Google Scholar]
  25. Ogihara, M., Kokubo, E., Suzuki, T. K., & Morbidelli, A. 2018, A&A, 615, A63 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  26. Ogihara, M., Kunitomo, M., & Hori, Y. 2020, ApJ, 899, 91 [NASA ADS] [CrossRef] [Google Scholar]
  27. Ormel, C. W., & Klahr, H. H. 2010, A&A, 520, A43 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  28. Paardekooper, S. J., Baruteau, C., & Kley, W. 2011, MNRAS, 410, 293 [Google Scholar]
  29. Pollack, J. B., Hubickyj, O., Bodenheimer, P., et al. 1996, Icarus, 124, 62 [NASA ADS] [CrossRef] [Google Scholar]
  30. Rogers, J. G., & Owen, J. E. 2021, MNRAS, 503, 1526 [Google Scholar]
  31. Sato, B., Fischer, D. A., Henry, G. W., et al. 2005, ApJ, 633, 465 [Google Scholar]
  32. Shibata, S., Helled, R., & Ikoma, M. 2020, A&A, 633, A33 [Google Scholar]
  33. Stewart, S. T., Lock, S. J., & Mukhopadhyay, S. 2014, Lunar and Planetary Science Conference, Lunar and Planetary Science Conference, 2869 [Google Scholar]
  34. Suzuki, T. K., & Inutsuka, S.-I. 2009, ApJ, 691, L49 [NASA ADS] [CrossRef] [Google Scholar]
  35. Suzuki, T. K., Ogihara, M., Morbidelli, A. R., Crida, A., & Guillot, T. 2016, A&A, 596, A74 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  36. Thorngren, D. P., Fortney, J. J., Murray-Clay, R. A., & Lopez, E. D. 2016, ApJ, 831, 64 [NASA ADS] [CrossRef] [Google Scholar]
  37. Wimarsson, J., Liu, B., & Ogihara, M. 2020, MNRAS, 496, 3314 [Google Scholar]

Appendix A: A summary of all envelope mass loss events for Fig. 2

Figure 3 shows the envelope mass loss during giant impacts only for simulations with the fiducial pebble flux. Here, Fig. A.1 summarizes the envelope mass loss for impact events in all simulations shown in Fig. 2. As stated in Sect. 4, there are collisions in which a large fraction of the envelope is lost in both cases with and without PEW effects. Regarding the dependence on the planetary mass, we find that the impact velocity and/or the envelope loss fraction tend to increase when the mass is larger. Since the planetary mass and the envelope mass fraction is positively correlated (see Fig. 2), the impact velocity tends to be higher for planets with high envelope mass fractions.

thumbnail Fig. A.1.

Same as Fig. 3, but for all simulations for Fig. 2. (a) Without PEW, (b) with PEW.

All Figures

thumbnail Fig. 1.

Evolution of the semimajor axis a, metal mass Mz, and envelope mass fraction Menv/M for our fiducial model in disks with and without PEW effects. The dashed line in panel a indicates that the growth of the core ceases at the pebble isolation mass if mutual collisions are ignored. Bottom panels: gas accretion rate at r = 0.3 au. (a) Without PEW, Fpb, fid, (b) with PEW, Fpb, fid.

In the text
thumbnail Fig. 2.

Comparison of the metal mass Mz and metal-mass fraction Mz/M between the observed transiting giant planets and the results of our simulations. Gray circles indicate the estimates by Thorngren et al. (2016), which were fitted by power-law distributions (dashed lines). HD 149026b and TOI-849b are plotted with purple and brown symbols, respectively. The simulation results are plotted with cresses and circles. The pebble flux was varied, and five runs were performed for each pebble flux. (a) Without PEW, (b) with PEW.

In the text
thumbnail Fig. 3.

Envelope mass-loss fraction for all the giant impacts in the simulations with the fiducial pebble flux. Circles with blue edges represent collisions in simulations without PEW effects, while those with orange edges are in simulations with PEW effects. Each size is proportional to mass, and the color corresponds to its envelope fraction before impacts. A summary of all mass loss events is presented in Fig. A.1.

In the text
thumbnail Fig. 4.

Results of additional simulations in which the envelope loss during giant impacts is ignored.

In the text
thumbnail Fig. A.1.

Same as Fig. 3, but for all simulations for Fig. 2. (a) Without PEW, (b) with PEW.

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.