Free Access
Issue
A&A
Volume 629, September 2019
Article Number A106
Number of page(s) 8
Section Planets and planetary systems
DOI https://doi.org/10.1051/0004-6361/201936202
Published online 12 September 2019

© ESO 2019

1 Introduction

The gaseous planets in the solar system, Jupiter and Saturn, have large regular satellites whose radii are larger than 1000 km, and their orbits are on the equatorial plane of the central planets. It is widely accepted that the gaseous planets were surrounded by a circumplanetary disk and that large satellites formed in the disks (e.g., Lunine & Stevenson 1982; Mosqueira & Estrada 2003; Canup & Ward 2006). Several hydrodynamic simulations revealed that circumplanetary disks can be formed as byproducts of the formation and growth of gaseous planets via mass accretion from the parent protoplanetary disk (e.g., Machida et al. 2010; Tanigawa et al. 2012).

Canup & Ward (2002) proposed the gas-starved disk model whose disk mass is regulated by the steady-state viscous accretion onto the central planet. Canup & Ward (2006) then performed N-body simulations of satellite formation in the gas-starved disk and revealed that the total mass of the satellites in circumplanetary disks is successfully reproduced. Moreover, Sasaki et al. (2010) found that the number, masses, and orbits of the Galilean satellites can be reproduced for the formation era when the disk had an inner cavity.

Jupiter has four large regular satellites called the Galilean satellites: Io, Europa, Ganymede, and Callisto. Saturn, in contrast, has only one large regular satellite, Titan, and many 100–1000 km-sized satellites. The inner three Galilean satellites orbit in a 4:2:1 mean motion resonance. One of the mechanisms for reproducing this orbital configuration is the stopping of the migration of Io near the bump in the surface density distribution and the subsequent resonant trapping of Europa and Ganymede (e.g., Sasaki et al. 2010; Fujii et al. 2017).

There have been two hypotheses to explain the formation of the bump in the surface density distribution of circumplanetary disks. Fujii et al. (2017) revealed that in the early stage of disk evolution, the disk temperature is high and the sublimation of silicate dust particles occurs within a certain orbital radius. This absence of dust particles changes the opacity (e.g., Bell & Lin 1994), causing the formation of the bump in the surface density distribution of the circumjovian disk. One problem in this scenario is whether the satellites can survive against the radial drift in the late stage of the circumjovian disk when the temperature of the disk edge is sufficiently lower than the sublimation temperature. Fujii et al. (2017) noted that the fate of satellites depends on the time evolution of the mass accretion rate onto Jupiter after the formation of planetary gas in the solar nebula (e.g., Dürmann & Kley 2015; Kanagawa et al. 2018). Another problem is that icy satellites, at a minimum, cannot survive longer than the disk dissipation timescale in such a high-temperature disk, because the temperature at the current location of icy satellites is about 1000 K at this stage (see Appendix A). We found that icy satellites at the ambient temperature of 500 K might be evaporated within 0.1 Myr, even if the satellite is as large as Ganymede.

On the other hand, Sasaki et al. (2010) proposed that the inner cavity of the circumjovian disk could be formed by the magnetic coupling between the disk and young Jupiter. This magnetic coupling can also explain why the spin rate of Jupiter is sufficiently slower than the break-up spin rate (e.g., Takata & Stevenson 1996; Batygin 2018). Nevertheless, one problem exists for this scenario: the spin period of young Jupiter should be longer than 40 h at the end of the Galilean satellite formation, as the innermost satellite should have formed at the location of the inner cavity (Ogihara & Ida 2012; Shibaike et al. 2019). This spin period is four times longer than the current spin period (approximately 10 h). Therefore, the planetary radius of young Jupiter may be at least twice as large as the current planetary radius (see Appendix B); although several evolutionary calculations have indicated that the planetary radius of young Jupiter might be approximately 1.5 times the current radius at 1 Myr after the formation of the planet (e.g., Fortney et al. 2011).

In this study, we propose another mechanism to create a bump in the surface density distribution of the circumjovian disk – the depletion of fine dust particles by photophoresis. Photophoresis is the momentum transfer mechanism whereby gas molecules collide with photo-irradiatied dust particles that exhibit a temperature gradient (e.g., Krauss & Wurm 2005; Loesche & Wurm 2012). The temperature gradient within dust particles is thought to form via direct irradiation from the central planet because the inner region of the circumjovian disk could be sufficiently depleted of dust particles. The inner region is therefore optically thin. The resulting dust-depleted inner region would have a higher ionization fraction, and would therefore admit increased accretion stress driven by magnetorotational instability (MRI; e.g., Balbus & Hawley 1991). This creates a bump in the gas surface density distribution and a pressure maximum at the boundary of the regions. We show that the photophoretic force in the circumjovian disk would overcome the gas drag, and the inward migration of dust particles would be stopped. We also reveal that this could form a pressure maximum and the bump in the surface density distribution of the circumjovian disk. In our model, the location of the bump is near the orbital radius of Io,and we do not need to consider the significant radius anomaly of young Jupiter. In addition, the temperatureof the disk at the current location of icy satellites could be lower than the sublimation temperature of H2O ice, and the resonant-trapped icy satellites may survive in the disk. Consequently, photophoresis in the circumjovian disk has the potential to explain the formation of the Galilean satellites.

2 Models

2.1 Photophoresis and gas drag

In the circumplanetary disk with gas and dust, whether or not dust particles drift inward is determined by the balance of the residual gravity caused by the gas drag, FD, and the photopho- retic force, FPh. In gaseous disks, other light-induced forces, such as the radiation pressure and the Poynting–Robertson effect, are negligible for μm- or mm-sized dust particles (Wurm & Krauss 2006, see also Appendix C).

For perfectly absorbing spherical particles whose radius is smaller than the mean free path of the gas molecules lm, the photophoretic force FPh is given by (Beresnev et al. 1993; Wurm & Krauss 2006) FPh=0.5πa2Ip3kthT/a+12σSBT4+p18kBT/(πmg),\begin{equation*} F_{\textrm{Ph}} = \frac{0.5 \pi a^{2} I p}{3 k_{\textrm{th}} T / a + 12 \sigma_{\textrm{SB}} T^{4} + p \sqrt{18 k_{\textrm{B}} T / {( \pi m_{\textrm{g}} )}}},\end{equation*}(1)

where a is the radius of dust particles, I is the intensity of the radiative flux from the planet, p is the gas pressure, kth is the thermal conductivity of dust particles, T is the temperature of the gas, σSB = 5.67 × 10−5 erg cm−2 K−4 s−1 is the Stefan–Boltzmann constant, kB = 1.38 × 10−16 erg K−1 is the Boltzmann constant, and mg = 3.9 × 10−24 g is the mean mass of gas molecules. We set a = 3 μm, which is consistent with the result of recent in-situ measurements of dust particles of the comet 67P/Churyumov-Gerasimenko (e.g., Bentley et al. 2016) and the estimation of the radius of constituent particles called monomers in the protoplanetary disk around the young protostar HL Tau (e.g., Okuzumi et al. 2016). The thermal conductivity of compact dust particles with no voids is about a few 10−1   W m−1 K−1 for organic polymers (e.g., Choy 1977) and ~1 W m−1 K−1 for glass beads (e.g., Sakatani et al. 2017). Therefore we set kth = 3 × 104 erg s−1 cm−1 K−1 (= 0.3 W m−1 K−1).

We checked the validity of the assumption that dust particles can be treated as perfectly absorbing particles for the radiation from the planet. The effective temperature of young Jupiter is approximately 1000 K (e.g., Fortney et al. 2011), and the characteristic wavelength of the radiation from the young Jupiter is λ ≃ 3 μm. Therefore, we can apply Eq. (1) for μm-sized dust particles whose size parameter x ≡ 2πaλ is larger than unity, that is, a ≳ 0.5 μm. We also confirmed the validity of the assumption of a < lm (see Appendix D).

The rotation of dust particles is not considered in this study, although fast rotation results in a decrease in the photophoretic force. For μm-sized dust particles, the main source of random rotation is Brownian motion. We confirmed that the thermal rotation timescale is several orders of magnitude longer than the heat-conduction timescale and that the effect of rotation is negligible (see Appendix E).

The gas in a circumplanetary disk usually rotates slower than the Keplerian velocity owing to the negative pressure gradient in the accretion disk. In contrast, dust particles in the gas disk are not supported by the pressure gradient, and their stable orbital velocity is Keplerian; however, fine dust particles do rotate with the gas. Therefore, there is a gas drag on dust particles, and the residual gravity FD is given by (Weidenschilling 1977; Wurm & Krauss 2006) FD=4πa33ρpρgdpdr,\begin{equation*} F_{\textrm{D}} = - \frac{4 \pi a^{3}}{3} \frac{\rho_{\textrm{p}}}{\rho_{\textrm{g}}} \frac{\textrm{d}p}{\textrm{d}r}, \end{equation*}(2)

where ρp = 2 g cm−3 is the density of dust particles, ρg is the gas density, and dp∕dr is the pressure gradient.

In the outer region of the circumplanetary disk, the residual gravity FD is larger than the photophoretic force FPh, and the inward drift dominates, even in the presence of direct irradiation. However, there is a turnover distance where FD = FPh, and the inward drift is prevented. Therefore, the inner region of the circumplanetary disk would be depleted of dust particles. Wurm & Krauss (2006) named this transition location the light barrier (see Fig. 1).

thumbnail Fig. 1

Schematic of the circumplanetary disk with a dust-poor inner region formed by photophoresis. Both the pressure maximum and the bump in the surface density distribution are at the light barrier.

2.2 Circumplanetary disk model

We model the structure of the circumplanetary disk as Fujii et al. (2014) and Shibaike et al. (2017). In this study, we assume the circumplanetary disk is a steady-state viscous accretion disk with a supply from the protoplanetary disk. The gas surface density of the circumplanetary disk Σ is given by Σ=M˙2πrbr3/2ν(29r1/2+23rbr3/2),\begin{equation*} \Sigma = \frac{\dot{M}}{2 \pi r_{\textrm{b}}} \frac{r^{3/2}}{\nu} {\left( - \frac{2}{9} r^{-1/2} &#x002B; \frac{2}{3} r_{\textrm{b}} r^{-3/2} \right)}, \end{equation*}(3)

where is the mass accretion rate, rb is the radius of the region where the gas falls in, r is the distance from the planet, and ν is the kinematic viscosity. The numerical simulation by Tanigawa et al. (2012) found that rb ≃ 20 RJ when the planet has a mass of 0.4MJ, where RJ = 7.15 × 104 km and MJ = 1.9 × 1030 g are the planetary radius and mass of Jupiter (see also Fujii et al. 2014), respectively. Assuming that rb scales with the Hill radius of the planet, Shibaike et al. (2017) suggested that rb = 27 RJ when the planet has a mass of 1MJ.

The gas temperature at the disk midplane T is given by (Nakamoto & Nakagawa 1994) T=[9νΣΩK28σSB(3τR8+12τP)]1/4,\begin{equation*} T = {\left[ \frac{9 \nu \Sigma {\Omega_{\textrm{K}}}^{2}}{8 \sigma_{\textrm{SB}}} {\left( \frac{3 \tau_{\textrm{R}}}{8} &#x002B; \frac{1}{2 \tau_{\textrm{P}}} \right)} \right]}^{1/4}, \end{equation*}(4)

where ΩK is the Keplerian angular velocity, and τR and τP are the Rosseland and Planck mean optical depths, respectively. The Keplerian angular velocity is given by ΩK=GM/r3$\Omega_{\textrm{K}} = \sqrt{{G M} / {r^{3}}}$, where G = 6.67 × 10−8 dyn cm2 g−2 is the gravitational constant, and M is the mass of the central planet. We set M = MJ in this study. The sound velocity cs is cs=kBT/mg$c_{\textrm{s}} = \sqrt{k_{\textrm{B}} T / m_{\textrm{g}}}$, the gas scale height of the disk hg is hg = cs∕ΩK, the gas density at the disk midplane is ρg=Σ/(2πhg)$\rho_{\textrm{g}} = \Sigma / {( \sqrt{2 \pi} h_{\textrm{g}} )}$, and the pressure at the disk midplane is p = ρgkBTmg. We employ the standard α-prescription for the kinematic viscosity ν as ν=αcs2/ΩK$\nu = \alpha {c_{\textrm{s}}}^{2} / \Omega_{\textrm{K}}$ (Shakura & Sunyaev 1973).

The Rosseland mean optical depth τR is given by τR=Σ2(fdκR,d+κR,g),\begin{equation*} \tau_{\textrm{R}} = \frac{\Sigma}{2} {\left( f_{\textrm{d}} \kappa_{\textrm{R, d}} &#x002B; \kappa_{\textrm{R, g}} \right)}, \end{equation*}(5)

where fd is the scaling factor of the dust abundance, and κR, d and κR, g are the contributions of dust and gas to the Rosseland mean opacity of a disk with solar metallicity. Similarly, the Planck mean optical depth τP is given by τP=Σ2(fdκP,d+κP,g),\begin{equation*} \tau_{\textrm{P}} = \frac{\Sigma}{2} {\left( f_{\textrm{d}} \kappa_{\textrm{P, d}} &#x002B; \kappa_{\textrm{P, g}} \right)}, \end{equation*}(6)

where κP, d and κP, g are the contributions of dust and gas to the Planck mean opacity of a disk with solar metallicity. The dust and gas opacities are taken from Henning & Stognienko (1996) and Freedman et al. (2014), respectively (see Appendix F).

Because the dust particles in the solar nebula tend to settle down toward the midplane, the dust-to-gas mass ratio of the inflow gas is thought to be lower than the dust-to-gas mass ratio of the solar nebula (e.g., Tanigawa et al. 2012). In addition, the dust-to-gas mass ratio in the inner region and the outer region must differ by orders of magnitude by the light barrier. Thus, we assumed that fd = 10−1 for the outer region, and fd = 10−6 for the inner region.

In this study, we assume that the abundance of nm- or sub-μm-sized dust particles is negligibly small, and κR, d and κP, d are determined by the abundance of μm-sized dust particles which affect the photophoretic force. Thus, we expect a drastic change of fd at the light barrier. We also neglect the turbulent diffusion of dust particles and consider the boundary at the light barrier to be sharp for simplicity. We note that the diffusion of dust particles may be important if the diffused region becomes optically thick and the intensity of the radiative flux at the light barrier decreases. Hence, future studies on this point are essential, although a detailed discussion is beyond the scope of this study.

The depletion of dust also affects the turbulent viscosity of the disk. Fujii et al. (2011, 2014) found that the midplane of circumplanetary disks is in general MRI-inactive because the diffusion timescale of the magnetic field is short. In the MRI-inactive region, it is likely that hydrodynamic instabilities, for example, the convective overstability (e.g., Klahr & Hubbard 2014) and subcritical baroclinic instability (e.g., Lyra & Klahr 2011), could dominate the angular momentum transport with values of α ~ 10−3 (e.g., Lyra 2014)1. In contrast, the disk might be (partly) MRI-active if it is significantly depleted of dust grains (e.g., Turner et al. 2014; Fujii et al. 2014), because the lower dust abundance leads to a higher ionization fraction of the disk gas. This is true when the dominant fraction of the total surface area of dust particles (where free charges can stick) is in μm-sized dust particles which will be stopped at the light barrier. In this case, the effective turbulent viscosity could be larger than the viscosity of the completely MRI-inactive disks. Thus, we assumed α = 10−3 for the outer region and α = 10−2 for the inner region. We note that the choice of absolute value of α does not change the results drastically; the increase in α at the light barrier is the key to generating the bump in the surface density distribution (see Appendix G).

3 Results

3.1 Location of the light barrier

First, we calculate the location of the light barrier in the steady-state accretion disk. At the location of the light barrier, the intensity of the radiative flux from the planet I should be I = L∕(4πr2), where L is the luminosity of the planet. As the rapid gas accretion was turned off and the planetary mass reached approximately 1 MJ, the luminosity of young Jupiter rapidly collapsed to 10−5L (e.g., Marley et al. 2007), where L = 3.83 × 1033 erg s−1 is the solar luminosity. Therefore, the luminosity of young Jupiter was on the order of 10−5 L to 10−6 L within the first 10 Myr (e.g., Fortney et al. 2011), and we set L = 4 × 10−6L in our calculation.

We calculated the force ratio of the photophoretic force to the gas drag, FPhFD, by using both the normal model (fd = 10−1 and α = 10−3) and the dust-poor model (fd = 10−6 and α = 10−2). The force ratio calculated from the normal model is (FPh/FD)n${\left( F_{\textrm{Ph}} / F_{\textrm{D}} \right)}_{\textrm{n}}$, and the force ratio calculated from the dust-poor model is (FPh/FD)p${\left( F_{\textrm{Ph}} / F_{\textrm{D}} \right)}_{\textrm{p}}$. We then defined the location of the light barrier rlb as the location where min{(FPh/FD)n,(FPh/FD)p}=1$\min {\left\{ {\left( F_{\textrm{Ph}} / F_{\textrm{D}} \right)}_{\textrm{n}}, {\left( F_{\textrm{Ph}} / F_{\textrm{D}} \right)}_{\textrm{p}} \right\}} = 1$, and we applied the normal model for the outer region whose distance is larger than rlb and applied the dust-poor model for the inner region whose distance is smaller than rlb.

Figure 2 shows the force ratio FPhFD under the assumption that the intensity of the radiative flux is I = L∕(4πr2). This estimate of I is true at the light barrier. We found that rlb = 6.9 RJ, 6.1 RJ, and 4.2 RJ for the cases of = 10−7MJ yr−1, 10−8MJ yr−1, and 10−9MJ yr−1, respectively.Therefore, the photophoretic force stops the inward migration of dust particles, and the inner region of the circumjovian disk is depleted of fine dust particles. Moreover, the location of this light barrier is close to the current orbit of the innermost satellite, Io (r = 5.9 RJ).

We note that for the case of = 10−9MJ yr−1, the force ratio is FPhFD = 1 not only at the light barrier (r = 4.2 RJ) but also at r = 10.9 RJ. We speculate that the true location of the light barrier is between 4.2 RJ and 10.9 RJ, although a detailed calculation of the dust and gas disk structure around the light barrier is needed to determine its location.

Behind the light barrier, the intensity of the radiative flux would decrease exponentially due to the increase in the optical depth from the light barrier. Therefore, the photophoretic force might only work around the light barrier, and the structure of the outer region might not be affected. However, the structure of the inner region must be affected by photophoresis, and our results do not change significantly. We also note that rlb does not significantly depend on the radius of dust particles when the two conditions, x ≫ 1 and alm, are satisfied (see Appendix H).

thumbnail Fig. 2

Force ratio FPhFD under the assumption that I = L∕(4πr2). The solid lines correspond to the force ratio calculated from the normal model (fd = 10−1 and α = 10−3), and the dashed lines are the dust-poor model (fd = 10−6 and α = 10−2). The vertical lines represent the location of the light barrier. The line colors represent the mass accretion rate ( = 10−7MJ yr−1 for blue, = 10−8MJ yr−1 for pink, and = 10−9MJ yr−1 for green). As a reference, the current orbit of Io is r = 5.9 RJ.

thumbnail Fig. 3

Pressure distribution of the disk. The solid lines correspond to the pressure p calculated from the normal model, and the dashed lines are the dust-poor model. The vertical lines represent the location of the light barrier, and a pressure maximum is formed at the light barrier.

3.2 Pressure and surface density distributions

Figure 3 shows the pressure distribution of the circumjovian disk. We found that a pressure maximum was formed at the light barrier. This structure originates from the change in the value of α at the light barrier in the steady-state accretion disk. We speculate that the radial drift of pebbles and satellitesimals might stop at the pressuremaximum (e.g., Haghighipour & Boss 2003) and help the inner region to maintain its dust-poor environment.

Figure 4 shows the surface density distribution. We found that the bump in the surface density distribution is formed at the light barrier when the mass accretion rate is lower than 10−7MJ yr−1. The migration direction and speed depend on the mass of the satellite and the disk structure (e.g., Tanaka et al. 2002; Paardekooper et al. 2011; Romanova et al. 2019). Fujii et al. (2017) found that the innermost satellite is stopped by the convergent migration to the bump. From our calculation, a steep bump in the surface density distribution is formed at the light barrier. Thus, the innermost satellite could be trapped around the light barrier.

Moreover, we confirmed the possibility of the formation of the resonant-trapped Galilean satellites by the following migration of Europa and Ganymede (see Appendix I). We found that the migration timescale is sufficiently longer than the critical timescale for resonance capture, and the adjacent satellites (i.e., Io–Europa and Europa–Ganymede) can be captured in the 2:1 mean motion resonance in the circumjovian accretion disk.

thumbnail Fig. 4

Surface density distribution of the disk. The solid lines correspond to the surface density Σ calculated from the normal model, and the dashed lines are the dust-poor model. The vertical lines represent the location of the light barrier, and a bump in the surface density distribution is formed at the light barrier.

thumbnail Fig. 5

Temperature distribution of the disk. The solid lines correspond to the temperature T calculated from the normal model, and the dashed lines are the dust-poor model.

3.3 Disk temperature

Figure 5 shows the temperature distribution. We found that the temperature at the light barrier is lower than the sublimation temperature of silicate dust (≳1500 K; Lodders 2003) when ≲ 10−7MJ yr−1. Therefore, photophoresis works before the sublimation of fine dust particles. Moreover, the temperature at the current location of Ganymede (r = 15 RJ) becomes lower than the sublimation temperature of H2O ice (≃200 K) when ≲ 10−8MJ yr−1. This result indicates that the Galilean satellites formed in the circumjovian disk with a mass accretion rate of ≲ 10−8MJ yr−1.

We note, however, that the disk temperature depends on the assumed opacity (e.g., Canup & Ward 2002), and the opacity should be calculated from the size distribution and the true abundance of the dust particles. The disk temperature is also dependent on the vertical structure of the accretion disk (e.g., Mori et al. 2019), and the vertical distribution of dust particles might be altered by photophoresis itself (e.g., McNally & Hubbard 2015; McNally & McClure 2017). Thus, future studies in these areas are essential.

4 Conclusions

The inner three Galilean satellites are in a 4:2:1 mean motion resonance, and their orbital configuration may be originated from the stopping of the migration of Io around the bump in surface density distribution, although the formation mechanism of the bump is under debate. We found that photophoresis in the circumjovian disk stopped the inward migration offine dust particles. The resulting dust-depleted inner region would have a higher ionization degree, and thus admit increased turbulent viscosity. This also forms a bump in the surface density distribution, and the migration of the innermost large satellite terminated around the bump. The location of the bump is near the current orbit of Io. Therefore, photophoresis is a strong candidate for the origin of the resonant-trapped Galilean satellites.

Acknowledgements

We thank Taishi Nakamoto for meaningful comments. S.A. is supported by the Grant-in-Aid for JSPS Research Fellow (JP17J06861). Y.S. is supported by a JSPS KAKENHI Grant (JP15H02065). Parts of this work have been carried out within the frame of the National Centre for Competence in Research PlanetS supported by the Swiss National Science Foundation (SNSF). Y.S. acknowledges the financial support of the SNSF.

Appendix A Evaporation of icy satellites

If the ambient temperature is much higher than the sublimation temperature of H2 O ice, then the icysatellites might be evaporated within the dissipation timescale of the circumplanetary disk. Here, we estimate the evaporation timescale. The mass and radius of the largest satellite, Ganymede, are MGa = 1.48 × 1026 g and RGa = 2.63 × 108 cm, and the evaporation timescale τevap is evaluated as τevapLevapMsat4πRsat2σSBT44.1×104(MsatMGa)(RsatRGa)2(T500 K)4 yr, \begin{equation*} \begin{split} \hspace*{-2pt}\tau_{\textrm{evap}} &\simeq \frac{L_{\textrm{evap}} M_{\textrm{sat}}}{4 \pi {R_{\textrm{sat}}}^{2} \sigma_{\textrm{SB}} T^{4}} \\ \hspace*{-2pt}&\simeq 4.1 \,{\times}\, 10^{4} {\left( \frac{M_{\textrm{sat}}}{M_{\textrm{Ga}}} \right)} {\left( \frac{R_{\textrm{sat}}}{R_{\textrm{Ga}}} \right)}^{-2} {\left( \frac{T}{500\ \textrm{K}} \right)}^{-4}\ \textrm{yr}, \end{split} \end{equation*}(A.1)

where Levap = 2.7 × 1010 erg g−1 is the latent heat of sublimation of H2O ice (Tanaka et al. 2013), Msat is the mass of the satellite, and Rsat is the radius of the satellite. Therefore, the evaporation timescale is shorter than the dissipation timescale of the circumplanetary disk (a few Myr). Moreover, the evaporation timescale is proportional to the radius Rsat and inversely proportional to the fourth power of the ambient temperature T, such that not only large satellites but also icy pebbles and satellitesimals cannot survive in the high-temperature environment.

Appendix B Inner cavity and corotation radius

If the inner cavity is located around the orbit of Io (r = 5.9 RJ), the orbital configuration of the Galilean satellites can be reproduced (Ogihara & Ida 2012; Shibaike et al. 2019). If the planetary magnetic field can be coupled effectively to the inner region of a circumplanetary accretion disk, then the disk could be disrupted at the radius rin (e.g., Königl 1991): rin=(μ42GMM˙2)1/75.2RJ(B150 G)4/7(M˙108MJyr1)2/7, \begin{equation*} \begin{split} r_{\textrm{in}} &= {\left( \frac{\mu^{4}}{2 G M {\dot{M}}^{2}} \right)}^{1/7} \\ \hspace*{-2pt}&\simeq 5.2 R_{\textrm{J}} {\left( \frac{B}{150\ \textrm{G}} \right)}^{4/7} {\left( \frac{\dot{M}}{10^{-8} M_{\textrm{J}}\ \textrm{yr}^{-1}} \right)}^{- 2/7}, \end{split} \end{equation*}(B.1)

where μ ~ BR3 is the dipole moment of the planet, B is the surface magnetic field of the planet, and R = RJ is the planetary radius. Although the current magnetic field of Jupiter is 7.8 G, according to the observations from the Juno spacecraft (Bolton et al. 2017; Connerney et al. 2017), young Jupiter might have a stronger magnetic field, on the order of 100 G (e.g., Sánchez-Lavega 2004; Reiners & Christensen 2010). Thus, we can expect that rin ≃ 5.9 RJ could be achieved if the ionization degree of the circumjovian disk is high enough, and the planetary magnetic field is coupled to the disk.

However, there is another requirement to make an inner cavity in a steady-state disk. To maintain its steady-state accretion, the corotation radius rco should not be smaller than the radius of the inner cavity rin (e.g., Königl 1991). The corotation radius is given by rco=(GMωspin2)1/35.3RJ(2π/ωspin36 h)2/3,\begin{equation*} r_{\textrm{co}} = {\left( \frac{G M}{{\omega_{\textrm{spin}}}^{2}} \right)}^{1/3} \,{\simeq}\, 5.3\, R_{\textrm{J}} {\left( \frac{2 \pi / \omega_{\textrm{spin}}}{36\ \textrm{h}} \right)}^{2/3}, \end{equation*}(B.2)

where ωspin is the angular velocity of the planetary spin. The rotation period of current Jupiter is 10 h; thus, the angular velocity of young Jupiter should have been about four times lower than that of the current Jupiter if the inner cavity were created by the magnetic coupling between young Jupiter and the circumjovian disk. In contrast, if rco > rin, a steady-state accretion disk cannot form, but a disk with no accretion might form. Instead, the spin rate of the central planet would decrease by propeller-driven outflows (e.g., Romanova et al. 2005). Then, the spin-down of the planet increases the corotation radius, reaching a stable state with rcorin, as Sasaki et al. (2010) predicted.

If we model the interior structure of the young Jupiter as a polytropic body with a polytropic index of 1.5, the spin angular momentum Lspin is given by (Batygin 2018) Lspin0.21MR2ωspin.\begin{equation*} L_{\textrm{spin}} \simeq 0.21 M R^{2} \omega_{\textrm{spin}}. \end{equation*}(B.3)

Consequentially, the planetary radius of young Jupiter would be as large as twice that of the current Jupiter if we assumed conservationof spin angular momentum. Evolutionary calculations of the planetary interior demonstrated, however, that the planetary radius of young Jupiter should have been about 1.5 times the current radius at 1 Myr after the formation of Jupiter (e.g., Fortney et al. 2011). Therefore, we need to consider the significant radius anomaly of young Jupiter in this scenario.

We note that Sasaki et al. (2010) argued that Io could form at its current position, even if the inner cavity was at the current corotation radius (r = 2.25 RJ); however, subsequent N-body simulations by Ogihara & Ida (2012) revealed that the positions of the cavity and the innermost satellite must have been the same.

Appendix C Radiation pressure

The radiation pressure is caused by the direct momentum transfer of photons. The radiation force Frad is given by Frad=Icπa2,\begin{equation*} F_{\textrm{rad}} = \frac{I}{c} \pi a^{2}, \end{equation*}(C.1)

where c = 3.00 × 1010 cm s−1 is the speed of light. Here, we compare the photophoretic force FPh and the radiation force Frad. The force ratio of the radiation force to the photophoretic force is FradFPh6kthTcap6×102(T300 K)(a3μm)1(p102 dyn cm2)1. \begin{equation*} \begin{split} \hspace*{-2pt}\frac{F_{\textrm{rad}}}{F_{\textrm{Ph}}} &\simeq \frac{6 k_{\textrm{th}} T}\textrm{c a p} \\ \hspace*{-2pt}&\simeq 6 \times 10^{-2} {\left( \frac{T}{300\ \textrm{K}} \right)} {\left( \frac{a}{3\ {\mu}\textrm{m}} \right)}^{-1} {\left( \frac{p}{10^{2}\ \textrm{dyn}\ \textrm{cm}^{-2}} \right)}^{- 1}. \end{split} \end{equation*}(C.2)

We then found that the radiation force is negligibly weaker than the photophoretic force in the gaseous circumplanetary disks (see also Loesche & Wurm 2012).

Appendix D Mean free path of gas molecules

The mean free path of the gas molecules lm is given by lm=mgσmolρg,\begin{equation*} l_{\textrm{m}} = \frac{m_{\textrm{g}}}{\sigma_{\textrm{mol}} \rho_{\textrm{g}}}, \end{equation*}(D.1)

where σmol = 2 × 10−15 cm2 is the collisional cross section of gas molecules (Adachi et al. 1976). In the circumjovian disk, lm is larger than 10−2 cm at the location of the light barrier (see Fig. D.1), and the photophoretic force for small dust particles of a = 3 μm is given by Eq. (1), as we assumed.

For spherical particles whose radius is smaller than the mean free path of the gas molecules, a < lm, the photophoretic force FPh is given by (Beresnev et al. 1993; Wurm & Krauss 2006) FPhπa3Ip6kthT,\begin{equation*} F_{\textrm{Ph}} \simeq \frac{\pi a^{3} I p}{6 k_{\textrm{th}} T}, \end{equation*}(D.2)

and therefore FPh does not depend on lm. On the other hand, when the dust radius is larger than the mean free path of the gas molecules, a > lm, the photophoretic force is approximately given by (Takeuchi & Krauss 2008) FPhπa3Ip6kthT×4.5(lma)2,\begin{equation*} F_{\textrm{Ph}} \,{\simeq}\, \frac{\pi a^{3} I p}{6 k_{\textrm{th}} T} \times 4.5 {\left( \frac{l_{\textrm{m}}}{a} \right)}^{2}, \end{equation*}(D.3)

and FPh decreases with increasing of a. In contrast to FPh, the residual gravity caused by gas drag FD is independent of alm (Weidenschilling 1977): FD=4πa33ρpρgdpdr.\begin{equation*} F_{\textrm{D}} = - \frac{4 \pi a^{3}}{3} \frac{\rho_{\textrm{p}}}{\rho_{\textrm{g}}} \frac{\textrm{d}p}{\textrm{d}r}. \end{equation*}(D.4)

thumbnail Fig. D.1

Mean free path of gas molecules in the disk. The solid lines correspond to the mean free path lm calculated from the normal model, and the dashed lines are the dust-poor model.

Therefore, FPhFD is independent of a for the case of a < lm, while FPhFD is inversely proportional to the square of the particle radius for the case of a > lm.

This means that the radial drift of cm-sized compact pebbles could not be prevented by photophoresis itself, although the radial drift would be stopped because the pressure maximum formed at the light barrier. We note, however, that the thermal conductivityof dust aggregates is exceedingly lower than that of compact dust particles (e.g., Arakawa et al. 2017, 2019), and the radial drift of aggregated pebbles might be stopped by photophoresis, even if their radius is larger than the mean free path of the gas molecules.

Appendix E Rotationof dust particles

If particles rotate faster than the heat conduction timescale, the temperature gradient would be reduced, and the photophoretic force might be suppressed. For μm-sized dust particles, the main source of particle rotation is the Brownian motion. The thermal rotation timescale of a spherical particle, τrot, is (Krauss & Wurm 2005; Loesche & Wurm 2012) τrot=8π2a5ρp45kBT1.4×102(a3μm)5/2(T300 K)1/2 s.\begin{equation*} \tau_{\textrm{rot}} = \sqrt{\frac{8 {\pi}^{2} a^{5} \rho_{\textrm{p}}}{45 k_{\textrm{B}} T}} \simeq 1.4 \times 10^{-2} {\left( \frac{a}{3\ {\mu}\textrm{m}} \right)}^{5/2} {\left( \frac{T}{300\ \textrm{K}} \right)}^{- 1/2}\ \textrm{s}. \end{equation*}(E.1)

The heat conduction timescale within a dust particle, τcond, is (Krauss & Wurm 2005; Loesche & Wurm 2012) τcond=a2kthρpcp6.0×105(a3μm)2 s,\begin{equation*} \tau_{\textrm{cond}} = \frac{a^{2}}{k_{\textrm{th}}} \rho_{\textrm{p}} c_{\textrm{p}} \simeq 6.0 \times 10^{-5} {\left( \frac{a}{3\ {\mu}\textrm{m}} \right)}^{2}\ \textrm{s}, \end{equation*}(E.2)

where cp = 107 erg g−1 K−1 is the specific heat (e.g., Ciesla et al. 2004). We found that τrotτcond, and the effect of particle rotation is negligible for μm-sized dust particles.

In addition to Brownian motion, dust-dust collisions spin up dust particles. However, this effect is negligibly small for a dust-to-gas mass ratio lower than unity (see McNally & McClure 2017).

Appendix F Dust and gas opacities

We calculated the opacity of the disk using standard temperature-dependent opacity models (Henning & Stognienko 1996; Freedman et al. 2014).

thumbnail Fig. F.1

Planck mean opacity of gas κP, g calculated by Freedman et al. (2014) (gray solid lines) and the analytical fit used in this study (blue dashed line).

The contribution of dust to the Rosseland mean opacity of a disk with solar metallicity, κR,d, is given by Henning & Stognienko (1996), and we use the iron-poor model of the study. The contribution of dust to the Rosseland mean opacity of a disk with solar metallicity, κP, d, is given by κP, d = 2.4κR, d (Nakamoto & Nakagawa 1994).

The contribution of gas to the Rosseland mean opacity of a disk with solar metallicity, κR,g, is given by Freedman et al. (2014), and we use their analytical fit of κR,g for the solar metallicity gas. We also evaluated κP, g, using Table 3 of Freedman et al. (2014) as κP,g=[15(T/100 K)8+0.15(T/100 K)0.8]1cm2g1.\begin{equation*} \kappa_{\textrm{P, g}} = {\left[ 15 {\left( T / {100\ \textrm{K}} \right)}^{-8} &#x002B; 0.15 {\left( T / {100\ \textrm{K}} \right)}^{0.8} \right]}^{-1}\ \textrm{cm}^{2}\ \textrm{g}^{-1}. \end{equation*}(F.1)

Within the temperature range of 100–2000 K, the pressure dependence of κP, g is small. Therefore, we assume κP, g is only the function of the temperature (see Fig. F.1).

We note that the dominant contributions to the Rosseland mean opacity of gas are in opacity windows, while the dominant contributions to the Planck mean opacity of gas are strong bands (see Fig. 2 of Freedman et al. 2014). Therefore, κP, g is several orders of magnitude larger than κR, g in general.

Appendix G Dependence on the turbulent viscosity

In the main text, we assumed α = 10−3 for the outerregion and α = 10−2 for the inner region. However, the choice of absolute values of α does not change the results significantly,but the increase in α at the light barrier is the key to generate the bump in the surface density distribution.

Figure G.1 shows the force ratio FPhFD for various viscosity parameters. Here, we assume that the value of α for the inner region is ten times higher than that for the outer region, and the mass accretion rate is = 10−8MJ yr−1. The line colors represent the value of α for the outerregion (α = 10−3 for pink, α = 5 × 10−4 for orange, and α =2.5 × 10−4 for blue). We found that the location of the light barrier is rlb = 6.1 RJ, 7.5 RJ, and 9.0 RJ for the cases of α =10−3, 5 × 10−4, and 2.5 × 10−4, respectively. The location of the light barrier is not strongly dependent on α and we can adjust rlb = 5.9 RJ by choosing appropriate values of the luminosity L and mass accretion rate .

thumbnail Fig. G.1

Force ratio FPhFD. The solid lines correspond to the force ratio calculated from the normal model, and the dashed lines are the dust-poor model. The vertical lines represent the location of the light barrier. The line colors represent the value of α for the outer region (α = 10−3 for pink, α = 5 × 10−4 for orange, and α = 2.5 × 10−4 for blue). Weassume that the value of α for the inner region is ten times higher than that for the outer region, and the mass accretion rate is = 10−8MJ yr−1.

Figure G.2 shows the surface density distribution of the disk for various viscosity parameters. We confirmed that a bump in the surface density distribution forms at the light barrier in all cases. Therefore, the general view of the disk structure hardly depends on the choice of absolute values of α, provided that the turbulence of the inner region is stronger than that of the outer region.

thumbnail Fig. G.2

Surface density distribution of the disk for various viscosity parameters. The solid lines correspond to the surface density Σ calculated from the normal model, and the dashed lines are the dust-poor model. The vertical lines represent the location of the light barrier, and a bump in the surface density distribution is formed at the light barrier.

Appendix H Dependence on the particle radius

For μm-sized dust particles in the circumjovian disk, we can simplify Eq. (1) as follows: FPhπa3Ip6kthT.\begin{equation*} F_{\textrm{Ph}} \simeq \frac{\pi a^{3} I p}{6 k_{\textrm{th}} T}. \end{equation*}(H.1)

Thus the force ratio FPhFD is given by FPhFDρgIp8ρpkthT(dpdr)1,\begin{equation*} \frac{F_{\textrm{Ph}}}{F_{\textrm{D}}} \,{\simeq}\, \frac{\rho_{\textrm{g}} I p}{8 \rho_{\textrm{p}} k_{\textrm{th}} T} {\left( \frac{\textrm{d}p}{\textrm{d}r} \right)}^{-1}, \end{equation*}(H.2)

and the location of the light barrier rlb is independent of the particle radius a, because both FPh and FD are proportional to the cubic of a.

Appendix I Condition for resonance capture

The Galilean satellites are believed to have migrated in the circumjovian disk by exchanging their angular momentum with the disk. If the migration timescale, τmig, is longer than the critical migration timescale for resonance capture, τcrit, the satellites are captured in the 2:1 mean-motion resonance. Ogihara & Kobayashi (2013) performed N-body simulations and revealed that the critical migration timescale for resonance capture τcrit is given by τcrit2×105(MinnerMIo)4/3ΩK1,\begin{equation*} \tau_{\textrm{crit}} \,{\simeq}\, 2 \times 10^{5} {\left( \frac{M_{\textrm{inner}}}{M_{\textrm{Io}}} \right)}^{- 4/3} {\Omega_{\textrm{K}}}^{-1}, \end{equation*}(I.1)

where MIo = 8.9 × 1025 g is the mass of Io and Minner is the mass of the inner satellite. The migration timescale τmig is given by (e.g., Paardekooper et al. 2011; Fujii et al. 2017) τmig15(MsatM)1(r2ΣM)1(csrΩK)2ΩK17×107(T300 K)(Σ300 g cm2)1(r9RJ)1(MsatMEu)1ΩK1, \begin{equation*} \begin{split} \tau_{\textrm{mig}} &\simeq \frac{1}{5} {\left( \frac{M_{\textrm{sat}}}{M} \right)}^{-1} {\left( \frac{r^{2} \Sigma}{M} \right)}^{-1} {\left( \frac{c_{\textrm{s}}}{r \Omega_{\textrm{K}}} \right)}^{2} {\Omega_{\textrm{K}}}^{-1} \\[3pt] \hspace*{-2pt}&\simeq 7 \times 10^{7} {\left( \frac{T}{300\ \textrm{K}} \right)} {\left( \frac{\Sigma}{300\ \textrm{g}\ \textrm{cm}^{-2}} \right)}^{-1} {\left( \frac{r}{9\, R_{\textrm{J}}} \right)}^{-1} {\left( \frac{M_{\textrm{sat}}}{M_{\textrm{Eu}}} \right)}^{-1} {\Omega_{\textrm{K}}}^{-1}, \end{split} \end{equation*}(I.2)

where MEu = 4.8 × 1025 g is the mass of Europa. We found that τmigτcrit and, therefore, the adjacent satellites (i.e., Io–Europa and Europa–Ganymede) can be captured in the 2:1 mean motion resonance.

The migration timescale τmig is on the order of a few 105 years for the Galilean satellites, and they can migrate in the circumjovian disk. We note, however, that τmig might be longer than the dissipation timescale of the disk when the mass of the satellite is far smaller than 10−2MEu. This fact indicates that Europa and Ganymede should be larger than 10−2 times the current mass at the time of resonance capture.

References

  1. Adachi, I., Hayashi, C., & Nakazawa, K. 1976, Prog. Theor. Phys., 56, 1756 [NASA ADS] [CrossRef] [Google Scholar]
  2. Arakawa, S., Tanaka, H., Kataoka, A., & Nakamoto, T. 2017, A&A, 608, L7 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  3. Arakawa, S., Tatsuuma, M., Sakatani, N., & Nakamoto, T. 2019, Icarus, 324, 8 [NASA ADS] [CrossRef] [Google Scholar]
  4. Balbus, S. A., & Hawley, J. F. 1991, ApJ, 376, 214 [NASA ADS] [CrossRef] [Google Scholar]
  5. Batygin, K. 2018, AJ, 155, 178 [NASA ADS] [CrossRef] [Google Scholar]
  6. Bell, K. R., & Lin, D. N. C. 1994, ApJ, 427, 987 [NASA ADS] [CrossRef] [Google Scholar]
  7. Bentley, M. S., Schmied, R., Mannel, T., et al. 2016, Nature, 537, 73 [NASA ADS] [CrossRef] [Google Scholar]
  8. Beresnev, S., Chernyak, V., & Fomyagin, G. 1993, Phys. Fluids A, 5, 2043 [NASA ADS] [CrossRef] [Google Scholar]
  9. Bolton, S. J., Adriani, A., Adumitroaie, V., et al. 2017, Science, 356, 821 [NASA ADS] [CrossRef] [Google Scholar]
  10. Canup, R. M., & Ward, W. R. 2002, AJ, 124, 3404 [NASA ADS] [CrossRef] [Google Scholar]
  11. Canup, R. M., & Ward, W. R. 2006, Nature, 441, 834 [NASA ADS] [CrossRef] [Google Scholar]
  12. Choy, C. L. 1977, Polymer, 18, 984 [CrossRef] [Google Scholar]
  13. Ciesla, F. J., Hood, L. L., & Weidenschilling, S. J. 2004, Meteorit. Planet. Sci., 39, 1809 [NASA ADS] [CrossRef] [Google Scholar]
  14. Connerney, J. E. P., Adriani, A., Allegrini, F., et al. 2017, Science, 356, 826 [NASA ADS] [CrossRef] [Google Scholar]
  15. Dürmann, C., & Kley, W. 2015, A&A, 574, A52 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  16. Fortney, J. J., Ikoma, M., Nettelmann, N., Guillot, T., & Marley, M. S. 2011, ApJ, 729, 32 [NASA ADS] [CrossRef] [Google Scholar]
  17. Freedman, R. S., Lustig-Yaeger, J., Fortney, J. J., et al. 2014, ApJS, 214, 25 [NASA ADS] [CrossRef] [Google Scholar]
  18. Fujii, Y. I., Okuzumi, S., & Inutsuka, S.-i. 2011, ApJ, 743, 53 [NASA ADS] [CrossRef] [Google Scholar]
  19. Fujii, Y. I., Okuzumi, S., Tanigawa, T., & Inutsuka, S.-i. 2014, ApJ, 785, 101 [NASA ADS] [CrossRef] [Google Scholar]
  20. Fujii, Y. I., Kobayashi, H., Takahashi, S. Z., & Gressel, O. 2017, AJ, 153, 194 [NASA ADS] [CrossRef] [Google Scholar]
  21. Haghighipour, N., & Boss, A. P. 2003, ApJ, 583, 996 [NASA ADS] [CrossRef] [Google Scholar]
  22. Henning, T., & Stognienko, R. 1996, A&A, 311, 291 [NASA ADS] [Google Scholar]
  23. Kanagawa, K. D., Tanaka, H., & Szuszkiewicz, E. 2018, ApJ, 861, 140 [NASA ADS] [CrossRef] [Google Scholar]
  24. Klahr, H., & Hubbard, A. 2014, ApJ, 788, 21 [Google Scholar]
  25. Königl, A. 1991, ApJ, 370, L39 [NASA ADS] [CrossRef] [Google Scholar]
  26. Krauss, O., & Wurm, G. 2005, ApJ, 630, 1088 [NASA ADS] [CrossRef] [Google Scholar]
  27. Lodders, K. 2003, ApJ, 591, 1220 [NASA ADS] [CrossRef] [Google Scholar]
  28. Loesche, C., & Wurm, G. 2012, A&A, 545, A36 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  29. Lunine, J. I., & Stevenson, D. J. 1982, Icarus, 52, 14 [NASA ADS] [CrossRef] [Google Scholar]
  30. Lyra, W. 2014, ApJ, 789, 77 [Google Scholar]
  31. Lyra, W., & Klahr, H. 2011, A&A, 527, A138 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  32. Machida, M. N., Kokubo, E., Inutsuka, S.-i., & Matsumoto, T. 2010, MNRAS, 405, 1227 [NASA ADS] [Google Scholar]
  33. Marley, M. S., Fortney, J. J., Hubickyj, O., Bodenheimer, P., & Lissauer, J. J. 2007, ApJ, 655, 541 [NASA ADS] [CrossRef] [Google Scholar]
  34. McNally, C. P., & Hubbard, A. 2015, ApJ, 814, 37 [NASA ADS] [CrossRef] [Google Scholar]
  35. McNally, C. P., & McClure, M. K. 2017, ApJ, 834, 48 [NASA ADS] [CrossRef] [Google Scholar]
  36. Mori, S., Bai, X.-N., & Okuzumi, S. 2019, ApJ, 872, 98 [NASA ADS] [CrossRef] [Google Scholar]
  37. Mosqueira, I., & Estrada, P. R. 2003, Icarus, 163, 198 [NASA ADS] [CrossRef] [Google Scholar]
  38. Nakamoto, T., & Nakagawa, Y. 1994, ApJ, 421, 640 [NASA ADS] [CrossRef] [Google Scholar]
  39. Ogihara, M., & Ida, S. 2012, ApJ, 753, 60 [NASA ADS] [CrossRef] [Google Scholar]
  40. Ogihara, M., & Kobayashi, H. 2013, ApJ, 775, 34 [NASA ADS] [CrossRef] [Google Scholar]
  41. Okuzumi, S., Momose, M., Sirono, S.-i., Kobayashi, H., & Tanaka, H. 2016, ApJ, 821, 82 [NASA ADS] [CrossRef] [Google Scholar]
  42. Paardekooper, S. J., Baruteau, C., & Kley, W. 2011, MNRAS, 410, 293 [NASA ADS] [CrossRef] [Google Scholar]
  43. Reiners, A., & Christensen, U. R. 2010, A&A, 522, A13 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  44. Romanova, M. M., Ustyugova, G. V., Koldoba, A. V., & Lovelace, R. V. E. 2005, ApJ, 635, L165 [NASA ADS] [CrossRef] [Google Scholar]
  45. Romanova, M. M., Lii, P. S., Koldoba, A. V., et al. 2019, MNRAS, 485, 2666 [NASA ADS] [CrossRef] [Google Scholar]
  46. Sakatani, N., Ogawa, K., Iijima, Y., et al. 2017, AIP Adv., 7, 015310 [NASA ADS] [CrossRef] [Google Scholar]
  47. Sánchez-Lavega, A. 2004, ApJ, 609, L87 [NASA ADS] [CrossRef] [Google Scholar]
  48. Sasaki, T., Stewart, G. R., & Ida, S. 2010, ApJ, 714, 1052 [NASA ADS] [CrossRef] [Google Scholar]
  49. Shakura, N. I., & Sunyaev, R. A. 1973, A&A, 24, 337 [NASA ADS] [Google Scholar]
  50. Shibaike, Y., Okuzumi, S., Sasaki, T., & Ida, S. 2017, ApJ, 846, 81 [NASA ADS] [CrossRef] [Google Scholar]
  51. Shibaike, Y., Ormel, C. W., Ida, S., Okuzumi, S., & Sasaki, T. 2019, ApJ, submitted [arXiv:1909.00285] [Google Scholar]
  52. Takata, T., & Stevenson, D. J. 1996, Icarus, 123, 404 [NASA ADS] [CrossRef] [Google Scholar]
  53. Takeuchi, T., & Krauss, O. 2008, ApJ, 677, 1309 [NASA ADS] [CrossRef] [Google Scholar]
  54. Tanaka, H., Takeuchi, T., & Ward, W. R. 2002, ApJ, 565, 1257 [NASA ADS] [CrossRef] [Google Scholar]
  55. Tanaka, K. K., Yamamoto, T., Tanaka, H., et al. 2013, ApJ, 764, 120 [NASA ADS] [CrossRef] [Google Scholar]
  56. Tanigawa, T., Ohtsuki, K., & Machida, M. N. 2012, ApJ, 747, 47 [NASA ADS] [CrossRef] [Google Scholar]
  57. Turner, N. J., Lee, M. H., & Sano, T. 2014, ApJ, 783, 14 [NASA ADS] [CrossRef] [Google Scholar]
  58. Weidenschilling, S. J. 1977, MNRAS, 180, 57 [NASA ADS] [CrossRef] [Google Scholar]
  59. Wurm, G., & Krauss, O. 2006, Icarus, 180, 487 [NASA ADS] [CrossRef] [Google Scholar]

1

We acknowledge, however, that these studies (Klahr & Hubbard 2014; Lyra & Klahr 2011) are for protoplanetary disks and we need to confirm whether or not hydrodynamic instabilities could dominate the angular momentum transport in the MRI-inactive region in the future.

All Figures

thumbnail Fig. 1

Schematic of the circumplanetary disk with a dust-poor inner region formed by photophoresis. Both the pressure maximum and the bump in the surface density distribution are at the light barrier.

In the text
thumbnail Fig. 2

Force ratio FPhFD under the assumption that I = L∕(4πr2). The solid lines correspond to the force ratio calculated from the normal model (fd = 10−1 and α = 10−3), and the dashed lines are the dust-poor model (fd = 10−6 and α = 10−2). The vertical lines represent the location of the light barrier. The line colors represent the mass accretion rate ( = 10−7MJ yr−1 for blue, = 10−8MJ yr−1 for pink, and = 10−9MJ yr−1 for green). As a reference, the current orbit of Io is r = 5.9 RJ.

In the text
thumbnail Fig. 3

Pressure distribution of the disk. The solid lines correspond to the pressure p calculated from the normal model, and the dashed lines are the dust-poor model. The vertical lines represent the location of the light barrier, and a pressure maximum is formed at the light barrier.

In the text
thumbnail Fig. 4

Surface density distribution of the disk. The solid lines correspond to the surface density Σ calculated from the normal model, and the dashed lines are the dust-poor model. The vertical lines represent the location of the light barrier, and a bump in the surface density distribution is formed at the light barrier.

In the text
thumbnail Fig. 5

Temperature distribution of the disk. The solid lines correspond to the temperature T calculated from the normal model, and the dashed lines are the dust-poor model.

In the text
thumbnail Fig. D.1

Mean free path of gas molecules in the disk. The solid lines correspond to the mean free path lm calculated from the normal model, and the dashed lines are the dust-poor model.

In the text
thumbnail Fig. F.1

Planck mean opacity of gas κP, g calculated by Freedman et al. (2014) (gray solid lines) and the analytical fit used in this study (blue dashed line).

In the text
thumbnail Fig. G.1

Force ratio FPhFD. The solid lines correspond to the force ratio calculated from the normal model, and the dashed lines are the dust-poor model. The vertical lines represent the location of the light barrier. The line colors represent the value of α for the outer region (α = 10−3 for pink, α = 5 × 10−4 for orange, and α = 2.5 × 10−4 for blue). Weassume that the value of α for the inner region is ten times higher than that for the outer region, and the mass accretion rate is = 10−8MJ yr−1.

In the text
thumbnail Fig. G.2

Surface density distribution of the disk for various viscosity parameters. The solid lines correspond to the surface density Σ calculated from the normal model, and the dashed lines are the dust-poor model. The vertical lines represent the location of the light barrier, and a bump in the surface density distribution is formed at the light barrier.

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.