## L etter to the E ditor

## Close-in giant-planet formation via in-situ gas accretion and their natal disk properties

Yasuhiro Hasegawa 1 , Tze Yeung Mathew Yu, 2 and Bradley M. S. Hansen 2

- 1 Jet Propulsion Laboratory, California Institute of Technology, Pasadena, CA 91109, USA e-mail: yasuhiro.hasegawa@jpl.nasa.gov
- 2 Department of Physics &amp; Astronomy, University of California Los Angeles, Los Angeles, CA 90095, USA

Received; accepted

## ABSTRACT

Aims. The origin of close-in Jovian planets is still elusive. We examine the in-situ gas accretion scenario as a formation mechanism of these planets.

Methods. We reconstruct natal disk properties from the occurrence rate distribution of close-in giant planets, under the assumption that the occurrence rate may reflect the gas accretion e GLYPH&lt;14&gt; ciency onto cores of these planets.

Results. We find that the resulting gas surface density profile becomes an increasing function of the distance from the central star with some structure at r ' 0 1 au. This profile is quite di : GLYPH&lt;11&gt; erent from the standard minimum-mass solar nebula model, while our profile leads to better reproduction of the population of observed close-in super-Earths based on previous studies. We compute the resulting magnetic field profiles and find that our profiles can be fitted by stellar dipole fields ( / r GLYPH&lt;0&gt; 3 ) in the vicinity of the central star and large-scale fields ( / r GLYPH&lt;0&gt; 2 ) at the inner disk regions, either if the isothermal assumption breaks down or if nonideal MHD e GLYPH&lt;11&gt; ects become important. For both cases, the transition between these two profiles occurs at r ' 0 1 au, which corresponds to the period : valley of giant exoplanets.

Conclusions. Our work provides an opportunity to test the in-situ gas accretion scenario against disk quantities, which may constrain the gas distribution of the minimum-mass extra solar nebula.

Key words. accretion, accretion disks - magnetic fields - turbulence - planets and satellites: formation - planets and satellites: gaseous planets - protoplanetary disks

## 1. Introduction

The dawn of exoplanetary science has come with the remarkable discovery of the Jovian-mass planet orbiting in the vicinity of its host, solar-type star (Mayor &amp; Queloz 1995). The extrasolar planets referred to as hot Jupiters have stimulated a number of follow-up studies to understand their origins (Dawson &amp; Johnson 2018, for a recent review). The rapid accumulation of exoplanet populations has revealed that close-in Jovian planets including hot Jupiters are statistically rare (Winn &amp; Fabrycky 2015, also see Figure 1), and their origin is still unclear.

Three modes of forming close-in Jupiters are currently offered. The first two modes are extensions of the canonical model: (proto)Jovian planets form beyond the snow line via the core accretion scenario (Pollack et al. 1996). The new ingredients are radial movement of planets during or after the process of forming. The first mode considers that inward planetary migration caused by planet-disk interaction plays the dominant role in reproducing the current orbital distribution (Lin et al. 1996; Kley &amp; Nelson 2012). The second one attributes such short orbital periods to planet-planet interaction coupled with stellar tides (Rasio &amp; Ford 1996; Matsumura et al. 2010). The third mode is called insitu gas accretion because it does not require large-scale radial movement (Bodenheimer et al. 2000; Batygin et al. 2016; Boley et al. 2016). Identification of which mode(s) would be most responsible for shaping the occurrence rate of close-in Jupiters has not yet been achieved.

Fig. 1. Occurrence rates of close-in giant planets as a function of orbital period. The horizontal bars represent the bin sizes used in estimating the occurrence rates, while the vertical ones denote the error bars. The data are taken from Santerne et al. (2016).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00647v1-figure-28b44376-2a02-4f31-a0ae-8ef639dd184e.jpg)

The lack of migration in the in-situ scenario implies that the formation e GLYPH&lt;14&gt; ciency of giant planets must be determined by the local properties of their natal disks. Given that the occurrence rate of observed exoplanets should directly correlate with the formation e GLYPH&lt;14&gt; ciency, it may be possible to invert exoplanet observa-

tions to infer the properties of the planet-forming disks. Here, we attempt to conduct such a study, using the observational results of Santerne et al. (2016). These authors employed Kepler transit observations combined with radial velocity follow-up and obtained the occurrence rate distribution of close-in Jovian planets. Figure 1 shows that the distribution possesses some intriguing structure. A similar structure was seen by Petigura et al. (2018) where only the Kepler data were used. 1 The fundamental hypothesis of this work is that the planet formation e GLYPH&lt;14&gt; ciency and hence the occurrence rate distribution are determined by the ability of giant planet cores to accrete the local disk gas. Under this hypothesis, we show below that the occurrence rate distribution can be converted to disk properties such as the gas surface density profile and magnetic field profile threading disks. Thus, our work provides not only disk quantities for examining the in-situ gas accretion scenario, but also a test bed for the gas distribution of the minimum mass extra solar nebula model (cf., Chiang &amp; Laughlin 2013).

## 2. Occurrence rate as a probe of inner disk structures

## 2.1. Basic hypothesis

We begin with the basic hypothesis that can be drawn from the assumption that the occurrence rate distribution ( f OR( )) of r close-in giant exoplanets may be viewed as a proxy for the probability of giant planet formation at r . Under the in-situ gasaccretion scenario, this probability should scale with the e GLYPH&lt;14&gt; -ciency for planetary cores to accrete the local disk gas. Consequently, one can obtain the relationship between f OR and inner disk properties (see below).

Gas accretion onto planets ( ˙ Mp ) can be divided into three stages in the core accretion scenario (e.g., Mordasini et al. 2012; Hasegawa et al. 2019). The first stage is envelope contraction around planetary cores characterized by the Kelvin-Helmholtz timescale. Given that the disk gas contained within the Hill radius of the cores with &gt; GLYPH&lt;24&gt; 10 M GLYPH&lt;8&gt; is a few Neptune masses at most at small orbital periods, this stage would end quickly and is not important for f OR. In the second stage, proto-Jovian planets keep accreting the disk gas, but the resulting accretion rate becomes too high and should eventually be limited by the local mass supply via disk evolution. This stage is the so-called disklimited gas accretion. Such e GLYPH&lt;14&gt; cient gas accretion is needed for the formation of giant planets, some of which would constitute the samples analyzed in Santerne et al. (2016). When planets are massive enough to open up gaps in their gas disks, ˙ Mp is determined by the gas coming from the polar direction (Morbidelli et al. 2014). While the final mass of planets should be determined during this stage, the value does not represent the local disk properties. Instead, the mass supply from outer disks over the gas disk lifetime of GLYPH&lt;24&gt; 10 6 yr controls this quantity (Bodenheimer et al. 2000). 2 Thus, the disk-limited gas accretion is the most likely determinant of the gas giant formation e GLYPH&lt;14&gt; ciency and hence f OR under the in-situ scenario. Mathematically, we formulate the disk-limited gas accretion rate as (Tanigawa &amp; Ikoma

1 Such a distribution is derived from planets that have radii of GLYPH&lt;24&gt; 4 GLYPH&lt;0&gt; 24 R GLYPH&lt;8&gt; (equivalently GLYPH&lt;24&gt; 3 M GLYPH&lt;8&gt; GLYPH&lt;0&gt; 80 M Jupiter ), while planets residing in the lower and upper mass ends are rare.

2 Under the steady state disk accretion model, ˙ Mp in this stage becomes an increasing function of r (e.g., see equations (8) and (9) in Hasegawa et al. 2019), which may be related to a radial distribution of the final mass of planets.

2007)

$$\text{Here, we } \int _ { \real } 2 0 0 / \real ) \\ \text{I results} \\ \text{and ob-} \\ \text{and ob-} \\ \text{triguing} \\ \int _ { \real } \int _ { \real } \int _ { \real } \int _ { \real } \int _ { \real } \int _ { \real } \int _ { \real } \int _ { \real } \int _ { \real } \int _ { \real } \int _ { \real } \int _ { \real } \int _ { \real } \int _ { \real } \int _ { \real } \int _ { \real } \int _ { \real } \int _ { \real } \int _ { \real } \int _ { \real } \int _ { \real } \int _ { \real } \int _ { \real } \int _ { \real } \int _ { \real } \int _ { \real } \int _ { \real } \int _ { \real } \int _ { \real } \int _ \real } \int _ { \real } \int _ { \real } \int _ { \real } \int _ { \real } \int _ { \real } \int _ { \real } \int _ { \real } \int _ { \real } \int _ { \real } \int _ { \real } \int _ { \real } \int _ { \real } \int _ { \real } \int _ { \real } \int _ { \real } \int _ { \real } \int _ { \real } \int _ { \real } \int _ { \real } \int _ { \real } \int _ { \real } \int _ { \real } \int _ { \real } \int _ { \real } \int _ { \real } \int _ { \real } \int _ { \real } \int _ { \real } \\ \text{mat by-} \\ \text{ncy and} \\ \text{the abil-} \\ \text{der this} \\ \text{This formula is derived from hvdrodynamical simulations (Tani-}$$

This formula is derived from hydrodynamical simulations (Tanigawa &amp; Watanabe 2002) and reliably reproduces the simulation results that were run independently by di GLYPH&lt;11&gt; erent groups (Hasegawa et al. 2018, references therein). We note that the steady state disk accretion model is adopted in equation (1), where the disk accretion rate ( ˙ M d ) is given as (Frank et al. 2002)

$$\dot { M } _ { d } & = 3 \pi \nu \Sigma _ { g },$$

where GLYPH&lt;6&gt; g is the gas surface density, GLYPH&lt;23&gt; = GLYPH&lt;11&gt; c H s g is the e GLYPH&lt;11&gt; ective viscosity, c s ( / T 1 2 = d ) is the sound speed, T d is the disk temperature, H g = c s = GLYPH&lt;10&gt; is the pressure scale height, and GLYPH&lt;10&gt; is the Keplerian angular velocity around a solar mass star. The GLYPH&lt;11&gt; -parameterization is used for quantifying the e GLYPH&lt;14&gt; ciency of angular momentum transport in disks (Shakura &amp; Sunyaev 1973).

Consequently, our hypothesis can be written as

$$\overset { \nu } { \text{itu gas} } ^ { m } \text{$\nu$}, \\ \overset { \text{itu eff} } { \text{Cons} } ^ { m } \text{$\nu$}, \\ \overset { \text{Cons} } { \text{Cons} } ^ { m } \text{$\nu$}, \\ \text{nd inner} ^ { m }$$

where k OR is an unknown factor related to the giant planet formation e GLYPH&lt;14&gt; ciency. This equation shows that a higher ˙ Mp leads to a higher f OR.

Thus, one can use f OR for specifying disk properties (either GLYPH&lt;11&gt; or T d ) under our hypothesis.

## 2.2. Disk model

Determination of T d becomes fundamental for exploring natal disk structures from f OR (see equation (3)). Here we describe how T d is computed self-consistently.

Protoplanetary disks can be heated by two sources (Armitage 2011): viscous heating and stellar irradiation. Viscous heating originates from turbulent accretion stress and generally becomes important in the inner region of the disk (D'Alessio et al. 1998). On the other hand, stellar irradiation regulates the disk temperature in the outer region of disk (Chiang &amp; Goldreich 1997). This work focuses on the inner region of disk and hence on viscous heating. Accordingly, T d is computed as under the optically thick assumption (Ruden &amp; Lin 1986; Nakamoto &amp; Nakagawa 1994):

$$\underset { \substack { e f o r m u - \\ \& \text{ Ikoma} } } { \text{ency and} } \, T _ { d } ^ { 4 } = \frac { 2 7 \tau } { 1 2 8 \sigma _ { S B } } \Sigma _ { g } \nu \Omega ^ { 2 } = \frac { 2 7 \kappa } { 1 2 8 \sigma _ { S B } } \left ( \frac { \dot { M } _ { d } } { 3 \pi } \right ) ^ { 2 } \frac { \Omega ^ { 3 } } { \alpha c _ { s } ^ { 2 } },$$

where GLYPH&lt;28&gt; = GLYPH&lt;20&gt; GLYPH&lt;6&gt; g is the optical depth, GLYPH&lt;20&gt; is the opacity, and GLYPH&lt;27&gt; SB is the Stefan-Boltzmann constant. We also assume in equation (4) that the vertical heat flux is constant and the temperature at the disk surface is negligible compared to that at the disk midplane (Hasegawa &amp; Pudritz 2011). Thus, T d can be computed self-consistently for a given value of GLYPH&lt;20&gt; , using equations (3) and (4).

We follow Bell et al. (1997) to specify GLYPH&lt;20&gt; that are functions of GLYPH&lt;6&gt; g and T d . In the preliminary study, we considered all the opacity regimes summarized in Table 2 of Bell et al. (1997, from n = 1 to n = 9, where n is the opacity regime index) and found that the results of all the cases are bracketed by those of n = 7 and n = 8 within the water snow line. The former case corresponds to the opacity originating from metal grains, while the latter corresponds to the opacity regime where metal grains sublimate due to high temperatures. Here, we focus only on these two opacity regimes, which are hereafter referred to as the metal grain and the sublimation cases, respectively.

In summary, there are three parameters in our model: k OR, ˙ M d , and Mp . We simply assume that k OR = 1. For ˙ M d , we adopt 10 GLYPH&lt;0&gt; 8 M GLYPH&lt;12&gt; yr GLYPH&lt;0&gt; 1 based on disk observations (Williams &amp; Cieza 2011). We set Mp to 10 M GLYPH&lt;8&gt; because this value may represent the typical core mass (see Section 3 for more discussion). We note that variation of these parameters does not a GLYPH&lt;11&gt; ect the profiles of the disk properties; variation only causes shifts in the absolute magnitude.

## 2.3. Resulting disk properties

We present all the disk properties ( T d , GLYPH&lt;6&gt; g , and GLYPH&lt;11&gt; ) that are computed self-consistently. For comparison purposes, we also consider two given temperature profiles: the isothermal case ( T d = 10 3 K), and the flat disk case: T d ' (2 = (3 GLYPH&lt;25&gt; )) 1 4 = ( R GLYPH&lt;3&gt; = r ) 3 4 = T GLYPH&lt;3&gt; , where R GLYPH&lt;3&gt; = 1 5 : R GLYPH&lt;12&gt; and T GLYPH&lt;3&gt; = 5780 K are the radius and temperature of the central star, respectively (Chiang &amp; Goldreich 1997). We note that T d is determined by stellar irradiation for the flat disk case.

Figure 2 shows our results for T d , GLYPH&lt;6&gt; g , and GLYPH&lt;11&gt; . We find that the radial variation of T d is small for the sublimation case. Interestingly, this behavior is very similar to that seen in the isothermal case. For the metal grain case, T d becomes a roughly decreasing function of r . This trend is broadly consistent with the results of the flat disk case. We confirm similar features for GLYPH&lt;6&gt; g . We can therefore conclude that for the opacity regime currently considered, the isothermal assumption works well for the sublimation case and the flat disk model provides reasonable results for the metal grain case. In the following, we consider only the isothermal and the flat disk cases. Our results also show that GLYPH&lt;6&gt; g traces the shape of f OR for all cases considered here (see Figure 1). The in-situ gas accretion hypothesis can therefore reproduce the occurrence rate of observed close-in planets if the gas surface density of disks has a similar distribution. We now discuss the GLYPH&lt;11&gt; profile (see right panel). We observe that the resulting GLYPH&lt;11&gt; profile is characterized by 1 = GLYPH&lt;6&gt; g , which is expected under the steady state disk model. We note that these disk properties should be viewed as the mean values due to plausible di GLYPH&lt;11&gt; erences in core mass and formation timing of close-in Jovian planets.

In the following, we use the above disk properties to compute the radial distribution of magnetic fields threading disks.

## 2.4. Magnetic field profile

Magnetic fields threading disks play a key role in disk evolution (Turner et al. 2014, for a review). Recent studies make it evident that the angular momentum transport of disks and the resulting GLYPH&lt;11&gt; are regulated by both nonideal MHD e GLYPH&lt;11&gt; ects and vertical magnetic flux. As shown in Figure 2, however, the inner region of disks becomes hot ( &gt; GLYPH&lt;24&gt; 800 K) enough for thermal and collisional ionization to be e GLYPH&lt;11&gt; ective (Gammie 1996). This indicates that the ideal MHD limit would still work well for most parts of inner disks. We therefore consider both ideal MHD and nonideal MHD cases here.

We use the results of detailed MHD simulations. For the ideal MHD case, the angular momentum of the disk is transported radially by magnetorotational instability (MRI) and the resulting MHD turbulence (Balbus &amp; Hawley 1998). We adopt the GLYPH&lt;11&gt; -expression given by Salvesen et al. (2016):

$$\intertext { \mathbb { m m m } } \, \text{app} _ { \cdot } \, \text{ \quad } \alpha _ { S 1 6 } = 1 1 \beta _ { z } ^ { - 0. 5 3 },$$

where GLYPH&lt;12&gt; z = ( GLYPH&lt;26&gt; g c 2 s ) = ( Bz 2 = 8 GLYPH&lt;25&gt; ) is the plasma beta, GLYPH&lt;26&gt; g = GLYPH&lt;6&gt; g = ( p 2 GLYPH&lt;25&gt; H g ) is the gas density at the disk midplane, and Bz is the net vertical magnetic flux threading disks.

When nonideal MHD e GLYPH&lt;11&gt; ects (Ohmic and ambipolar di GLYPH&lt;11&gt; usion and Hall term) come into play, the picture of disk evolution changes considerably. In fact, recent nonideal MHD simulations reveal that inclusion of Ohmic and ambipolar di GLYPH&lt;11&gt; usion with large-scale magnetic fields can lead to the quench of MRIinduced turbulence and the launch of disk winds that remove the angular momentum of the disk in the vertical direction (Bai &amp; Stone 2013; Bai 2013). When all three nonideal MHD terms are included, the results become sensitive to the direction of magnetic fields (aligned or anti-aligned with the rotation axis of the disk) while disks tend to be laminar and disk winds still play a considerable role in disk evolution (Bai 2014; Lesur et al. 2014). The magnetic field configuration relative to the rotation axis of the disk is unknown, and the results with Ohmic and ambipolar di GLYPH&lt;11&gt; usion provide the intermediate behavior between the results of full nonideal MHD simulations with the aligned and antialigned cases (Bai 2013, 2014). Therefore, here we use the results of nonideal MHD simulations with Ohmic and ambipolar di GLYPH&lt;11&gt; usion. Assuming that angular momentum is removed largely by disk winds vertically in almost laminar disks, equation (2) can be modified as (Hasegawa et al. 2017)

$$\underset { \mathfrak { W } \, \equiv \, \infty } { \text{creating} } \\ \underset { \mathfrak { W } \, \equiv \, \infty } { \text{creating} } \\ \underset { \mathfrak { W } \, \equiv \, \infty } { \text{$\mathfrak { W } } } \, \overset { \cdot } { M _ { d } } \simeq 4 \sqrt { 2 \pi } r \Sigma _ { g } c _ { s } W _ { z \phi }, \,,$$

where W z GLYPH&lt;30&gt; is the normalized accretion stress originating from disk winds. We adopt the functional form of W z GLYPH&lt;30&gt; obtained by Bai (2013):

$$\colon \text{ison} ^ { - } _ { g } \text{ traces} \quad W _ { z \phi } ^ { B 1 3 } = 0. 2 3 \left ( \frac { r } { 1 \, \text{au} } \right ) ^ { 0. 4 6 } \beta _ { z } ^ { - 0. 6 6 }.$$

Consequently, the e GLYPH&lt;11&gt; ective GLYPH&lt;11&gt; can be written as

$$\underset { \text{div state} } { \text{uss the $\alpha$} } } \quad \alpha _ { B 1 3 } \equiv \frac { 4 r } { 3 \sqrt { \pi } H _ { g } } W _ { z \phi } ^ { B 1 3 }.$$

We now compute the magnetic field profiles ( Bz ) that are reconstructed from the GLYPH&lt;11&gt; profiles (Figure 2), using equations (5) and (8). Figure 3 shows our results. We find that there are plateaus in the resulting magnetic field profiles around r ' 0 1 au : for all cases. Importantly, these locations correspond to the transition between hot Jupiters and giant planets in the period valley. Our results also show that except for the isothermal case in the ideal MHD limit, Bz switches from r GLYPH&lt;0&gt; 3 to r GLYPH&lt;0&gt; 2 as r increases. It is clear that the switch in profiles originates from f OR.

Therefore, the occurrence rate distribution can be converted to disk properties ( GLYPH&lt;6&gt; g , GLYPH&lt;11&gt; , and Bz ) under our hypothesis.

## 3. Discussion

We first point out that the resulting profiles of GLYPH&lt;6&gt; g and GLYPH&lt;11&gt; are quite di GLYPH&lt;11&gt; erent from those generated by the standard model: our GLYPH&lt;6&gt; g

Fig. 2. Resulting structures of T d , GLYPH&lt;6&gt; g , and GLYPH&lt;11&gt; as a function of the distance from the central star in the left, central, and right panels, respectively. The results for the isothermal case are denoted by the red solid line, those for the sublimation case by the green dotted line, those for the flat disk case by the black dashed line, and those for the metal grain case by the orange dash-dotted line. As in Figure 1, the horizontal bars are displayed while error bars are omitted for the purpose of clear presentation in the left and central panels. Both the horizontal bars and error bars are displayed in the right panel.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00647v1-figure-9d92fdf1-36e6-4577-8f38-cd1d3c976877.jpg)

Fig. 3. Resulting magnetic field profiles as a function of r ; as in the right panel of Figure 2. The results for the ideal MHD case are shown in the left panel, and those for the nonideal MHD case are shown in the right panel. For reference, the profiles of r GLYPH&lt;0&gt; 2 and r GLYPH&lt;0&gt; 3 are denoted by the blue and the purple solid lines, respectively.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00647v1-figure-8091cc4c-4310-4519-b28f-cae8d6901e22.jpg)

is a roughly increasing function of r with some structures (Figure 2), while the minimum-mass solar nebula (MMSN) model goes to r GLYPH&lt;0&gt; 3 2 = (Hayashi 1981). The gas distribution in the inner region of disks is poorly constrained by the MMSN model due to the absence of gas giants there. It is worth mentioning that the importance of increasing surface density profiles with r has already been recognized in the literature. Ogihara et al. (2015) suggest that with such a profile, the population of observed close-in super-Earths can be better reproduced. We note that they assume that the origin of the profile is mass loss originating from disk winds without considering magnetic field profiles. Our results would be reasonable within &lt; GLYPH&lt;24&gt; 1 au because the opacity regimes currently considered are e GLYPH&lt;11&gt; ective there; other opacity regimes become important beyond this region, which changes the gas surface density and disk temperature profiles, and the importance of nonideal MHD e GLYPH&lt;11&gt; ects. The recent analysis of radial velocity observations suggests that the occurrence rate of giant planets may have a peak at r = 2 GLYPH&lt;0&gt; 3 au (Fernandes et al. 2019). This implies that other formation mechanisms would operate there as well.

For GLYPH&lt;11&gt; , our profile broadly predicts a decreasing function of r with some features around r ' 0 1 au (Figure 2). In particular, : one would notice that GLYPH&lt;11&gt; &gt; GLYPH&lt;24&gt; 1 in the vicinity of the central star. We remind the readers that the absolute values of disk quantities are scaled by model parameters ( k OR, Mp , and ˙ M d ) and that k OR is totally unknown. For reference, we have computed the plasma GLYPH&lt;12&gt; z value with equations (5) and (7) and confirmed that under the current setup, the values of GLYPH&lt;11&gt; at r &gt; GLYPH&lt;24&gt; 0 04 au and at : r &gt; GLYPH&lt;24&gt; 0 4 : au are achieved by recent ideal and nonideal MHD simulations, respectively (Salvesen et al. 2016; Bai 2013).

The most intriguing quantity would be magnetic field profiles. While it is very di GLYPH&lt;14&gt; cult to infer the magnetic field profile observationally, a number of theoretical studies are available. For instance, it is well recognized that disk inner edges may be threaded by stellar magnetic fields and the resulting magnetospheres would be important for disk evolution in the vicinity of classical T Tauri stars (Ghosh &amp; Lamb 1979; Koenigl 1991). If the dipole approximation is reasonable, the stellar vertical magnetic flux ( B s ) threading disks can be written as (e.g., Armitage 2016)

$$\underset { 1 9 } { \text{$t$ panels} } \text{$B_{s}$} \simeq 1 0 ^ { 3 } \left ( \frac { r } { 1. 5 R _ { \odot } } \right ) ^ { - 3 } \, G \simeq 4 1 \left ( \frac { r } { 2 \times 1 0 ^ { - 2 } \, \text{au} } \right ) ^ { - 3 } \, G.$$

We note that the averaged stellar surface magnetic field strength is adopted in equation (9), which is on the order of 10 3 G for classical T Tauri stars (e.g., Johns-Krull 2007). Such strength is probably needed for truncating the disk inner edge and making magnetospheric accretion e GLYPH&lt;11&gt; ective (Koenigl 1991). How-

ever, spectropolarimetry observations suggest that quadrupolar and higher-order moments provide the dominant contribution to the strength (e.g., Donati et al. 2008), implying that the dipole field may be much weaker (Hartmann et al. 2016, for a review). In fact, Russo &amp; Thompson (2015b) adopt a weak stellar field and compute the disk properties in detail (also see Russo &amp; Thompson 2015a). Therefore, equation (9) should be viewed as an upper limit for stellar dipole fields. In addition, protoplanetary disks would be threaded by large-scale magnetic fields ( B d ) that originate from interstellar fields that are the consequence of disk formation and evolution inherited from star formation (Lubow et al. 1994; Okuzumi et al. 2014). For this case, B d is determined by two competing processes: advection and di GLYPH&lt;11&gt; usion of magnetic flux in the radial direction. Adopting the steady state solution, B d can be given as (Okuzumi et al. 2014)

$$B _ { d } & \simeq 0. 1 \left ( \frac { r } { 1 \, \text{au} } \right ) ^ { - 2 } \, G. & ( 1 0 ) & \text{Armitage}, \\ & \text{Armitage},$$

This should be viewed as an upper limit for the large-scale magnetic fields of disks because it is derived from the assumption of highly conducting accretion disks. Interestingly, our results of Bz are comparable with these two upper limits (Figure 3). This work therefore implies that stellar fields might a GLYPH&lt;11&gt; ect the formation of hot Jupiters while large-scale fields might be important for giant planet formation beyond r ' 0 1 au. Furthermore, switching : these two fields might lead to the period valley in f OR.

How do our disk properties a GLYPH&lt;11&gt; ect planetary migration? As shown by previous work, type I migration is halted at the disk inner edge (Masset et al. 2006; Ogihara et al. 2015). For type II migration, its speed can be reduced significantly if the planet mass at r = rp exceeds the local disk mass ( GLYPH&lt;1&gt; M d ' 2 GLYPH&lt;25&gt; GLYPH&lt;6&gt; g r 2 p , e.g., Hasegawa &amp; Ida 2013). At rp ' 1 au, GLYPH&lt;6&gt; g GLYPH&lt;24&gt; 10 2 g cm GLYPH&lt;0&gt; 2 and GLYPH&lt;1&gt; M d GLYPH&lt;24&gt; 24 M GLYPH&lt;8&gt; for our current setup. Thus, both type I and II migration would be negligible in our model, which justifies the assumption of the in-situ scenario.

We now list the potential issues contained in this work. The most critical one is that we have adopted Mp = 10 M GLYPH&lt;8&gt; in equation (3). The presence of such massive cores in the inner region of disks is fundamental for the in-situ gas accretion scenario. Nonetheless, equation (3) implicitly assumes that the formation frequency of such cores may be uniform there. Since Kepler observations show that the populations of both super-Earths and sub-Neptunes increase monotonically from 1 day to 10 days and become relatively flat beyond 10 days (e.g., Petigura et al. 2018), our assumption would be justified only for longer periods. In other words, low f OR within 10 days can be caused by either low GLYPH&lt;6&gt; g or low abundance of planetary cores. We also note that the Kepler data lead to a similar occurrence rate distribution even with di GLYPH&lt;11&gt; erent analyses (Santerne et al. 2016; Petigura et al. 2018). On the contrary, the analysis of only radial velocity data tends to gain a higher value of the occurrence rate for close-in Jovian planets (Mayor et al. 2011). More observations and or system-/ atic analyses are needed to examine the di GLYPH&lt;11&gt; erence.

We finally discuss how close-in giant planets form. Our results clearly show that the in-situ gas accretion can work for planets located beyond 0.1 au. For planets within 0.1 au, it may not be impossible to form some fraction of hot Jupiters via the in-situ gas accretion scenario as suggested by Bodenheimer et al. (2000) and this work. However, the scenario would entail certain problems, such as the need to populate massive cores and or un-/ reasonably high values of GLYPH&lt;11&gt; there. Instead, other scenarios may provide reasonable explanations for such planets. For instance, the peak value of f OR at GLYPH&lt;24&gt; 4 days can be understood by the combination of high-eccentricity migration and tidal circularization (Wu &amp; Lithwick 2011), and other orbital features such as obliquity can be explained by the Lidov-Kozai e GLYPH&lt;11&gt; ect (e.g., Naoz 2016).

In the near future, better modeling and observations of inner disk regions will become available to directly test our prediction and eventually to reveal the appearance of natal protoplanetary disks and how planet formation takes place in the vicinity of central stars.

Acknowledgements. The authors thank an anonymous referee for useful comments on our manuscript. Y.H. thanks Konstantin Batygin and Trevor David for stimulating discussions. This research was carried out at JPL Caltech under a / contract with NASA. Y.H. is supported by JPL Caltech. /

## References

Armitage, P. J. 2011, ARA&amp;A, 49, 195 Armitage, P. J. 2016, ApJ, 833, L15 Bai, X.-N. 2013, ApJ, 772, 96 Bai, X.-N. 2014, ApJ, 791, 137 Bai, X.-N. &amp; Stone, J. M. 2013, ApJ, 769, 76 Balbus, S. A. &amp; Hawley, J. F. 1998, Reviews of Modern Physics, 70, 1 Batygin, K., Bodenheimer, P. H., &amp; Laughlin, G. P. 2016, ApJ, 829, 114 Bell, K. R., Cassen, P. M., Klahr, H. H., &amp; Henning, T. 1997, ApJ, 486, 372 Bodenheimer, P., Hubickyj, O., &amp; Lissauer, J. J. 2000, Icarus, 143, 2 Boley, A. C., Granados Contreras, A. P., &amp; Gladman, B. 2016, ApJ, 817, L17 Chiang, E. &amp; Laughlin, G. 2013, MNRAS, 431, 3444 Chiang, E. I. &amp; Goldreich, P. 1997, ApJ, 490, 368 D'Alessio, P., Cantó, J., Calvet, N., &amp; Lizano, S. 1998, ApJ, 500, 411 Dawson, R. I. &amp; Johnson, J. A. 2018, ARA&amp;A, 56, 175 Donati, J. F., Morin, J., Petit, P., et al. 2008, MNRAS, 390, 545 Fernandes, R. B., Mulders, G. D., Pascucci, I., Mordasini, C., &amp; Emsenhuber, A. 2019, ApJ, 874, 81

- Frank, J., King, A., &amp; Raine, D. J. 2002, Accretion Power in Astrophysics: Third Edition
- Gammie, C. F. 1996, ApJ, 457, 355
- Ghosh, P. &amp; Lamb, F. K. 1979, ApJ, 232, 259
- Hartmann, L., Herczeg, G., &amp; Calvet, N. 2016, ARA&amp;A, 54, 135
- Hasegawa, Y., Bryden, G., Ikoma, M., Vasisht, G., &amp; Swain, M. 2018, ApJ, 865, 32
- Hasegawa, Y., Hansen, B. M. S., &amp; Vasisht, G. 2019, ApJ, 876, L32
- Hasegawa, Y. &amp; Ida, S. 2013, ApJ, 774, 146
- Hasegawa, Y., Okuzumi, S., Flock, M., &amp; Turner, N. J. 2017, ApJ, 845, 31
- Hasegawa, Y. &amp; Pudritz, R. E. 2011, MNRAS, 417, 1236
- Hayashi, C. 1981, Progress of Theoretical Physics Supplement, 70, 35
- Johns-Krull, C. M. 2007, ApJ, 664, 975
- Kley, W. &amp; Nelson, R. P. 2012, ARA&amp;A, 50, 211
- Koenigl, A. 1991, ApJ, 370, L39
- Lesur, G., Kunz, M. W., &amp; Fromang, S. 2014, A&amp;A, 566, A56
- Lin, D. N. C., Bodenheimer, P., &amp; Richardson, D. C. 1996, Nature, 380, 606
- Lubow, S. H., Papaloizou, J. C. B., &amp; Pringle, J. E. 1994, MNRAS, 267, 235
- Masset, F. S., Morbidelli, A., Crida, A., &amp; Ferreira, J. 2006, ApJ, 642, 478
- Matsumura, S., Peale, S. J., &amp; Rasio, F. A. 2010, ApJ, 725, 1995
- Mayor, M. &amp; Queloz, D. 1995, Nature, 378, 355
- Mayor, M. et al. 2011, preprint (astro-ph arXiv:1109.2497v1) /
- Morbidelli, A., Szulágyi, J., Crida, A., et al. 2014, Icarus, 232, 266
- Mordasini, C., Alibert, Y., Klahr, H., &amp; Henning, T. 2012, A&amp;A, 547, A111
- Nakamoto, T. &amp; Nakagawa, Y. 1994, ApJ, 421, 640
- Naoz, S. 2016, ARA&amp;A, 54, 441
- Ogihara, M., Morbidelli, A., &amp; Guillot, T. 2015, A&amp;A, 584, L1
- Okuzumi, S., Takeuchi, T., &amp; Muto, T. 2014, ApJ, 785, 127
- Petigura, E. A., Marcy, G. W., Winn, J. N., et al. 2018, AJ, 155, 89
- Pollack, J. B., Hubickyj, O., Bodenheimer, P., et al. 1996, Icarus, 124, 62
- Rasio, F. A. &amp; Ford, E. B. 1996, Science, 274, 954
- Ruden, S. P. &amp; Lin, D. N. C. 1986, ApJ, 308, 883

Russo, M. &amp; Thompson, C. 2015a, ApJ, 815, 38

- Russo, M. &amp; Thompson, C. 2015b, ApJ, 813, 81
- Salvesen, G., Simon, J. B., Armitage, P. J., &amp; Begelman, M. C. 2016, MNRAS, 457, 857
- Santerne, A., Moutou, C., Tsantaki, M., et al. 2016, A&amp;A, 587, A64
- Shakura, N. I. &amp; Sunyaev, R. A. 1973, A&amp;A, 24, 337
- Tanigawa, T. &amp; Ikoma, M. 2007, ApJ, 667, 557
- Tanigawa, T. &amp; Watanabe, S.-i. 2002, ApJ, 580, 506
- Turner, N. J., Fromang, S., Gammie, C., et al. 2014, Protostars and Planets VI, 411
- Williams, J. P. &amp; Cieza, L. A. 2011, ARA&amp;A, 49, 67

Winn, J. N. &amp; Fabrycky, D. C. 2015, ARA&amp;A, 53, 409

Wu, Y. &amp; Lithwick, Y. 2011, ApJ, 735, 109
Yasuhiro Hasegawa Tze Yeung Mathew Yu Bradley M. S. Hansen
Close-in giant-planet formation via in-situ gas accretion and their natal disk properties
The origin of close-in Jovian planets is still elusive. We examine the
in-situ gas accretion scenario as a formation mechanism of these planets. We
reconstruct natal disk properties from the occurrence rate distribution of
close-in giant planets, under the assumption that the occurrence rate may
reflect the gas accretion efficiency onto cores of these planets. We find that
the resulting gas surface density profile becomes an increasing function of the
distance from the central star with some structure at $r \simeq 0.1$ au. This
profile is quite different from the standard minimum-mass solar nebula model,
while our profile leads to better reproduction of the population of observed
close-in super-Earths based on previous studies. We compute the resulting
magnetic field profiles and find that our profiles can be fitted by stellar
dipole fields ($\propto r^{-3}$) in the vicinity of the central star and
large-scale fields ($\propto r^{-2}$) at the inner disk regions, either if the
isothermal assumption breaks down or if nonideal MHD effects become important.
For both cases, the transition between these two profiles occurs at $r \simeq
0.1$ au, which corresponds to the period valley of giant exoplanets. Our work
provides an opportunity to test the in-situ gas accretion scenario against disk
quantities, which may constrain the gas distribution of the minimum-mass {\it
extra}solar nebula.