{"id":"815d2fff-b5f2-4780-838b-2dd2ae296ecf","arxiv_id":"1908.00647","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Under the assumption that occurrence rates reflect gas accretion efficiency, the paper inverts Kepler transit data to infer inner disk gas surface density and magnetic field profiles that switch character at the hot-Jupiter period valley.","lead":"This paper turns the observed occurrence rate of close-in giant planets into a reconstruction of the gas and magnetic field structure of the inner planet-forming disk. The resulting surface-density profile rises with distance, and the inferred magnetic field switches from a stellar-dipole-like r^-3 to a large-scale r^-2 profile at about 0.1 au, the location of the observed period valley.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Uniform 10-M⊕ core formation frequency is the load-bearing premise; Kepler data contradict it inside 10 days, so the 0.1 au break may be an artifact of the assumption.","rationale":"The paper's strongest claim is that fOR can be inverted to recover natal disk properties under the in-situ gas-accretion hypothesis. That inversion is only meaningful if the mapping from fOR to gas accretion efficiency is one-to-one; the paper acknowledges in Section 3 that the formation frequency of 10-M⊕ cores is assumed uniform, and that Kepler observations of smaller planets contradict this within 10 days. Because the period valley and the 0.1 au feature are extracted from the same fOR data that is being explained, the claimed 'correspondence' is not an independent test. The reader's CONDITIONAL verdict is appropriate and I agree with the reader's weakest_assumption. A test using super-Earth occurrence as f_core(r) would settle whether the reconstructed profiles are robust. Even though α and Bz magnitudes scale with the unknown kOR and adopted Mdot_d, the paper's headline claims are shape-based; the decisive issue is the radial uniformity of core formation, not parameter calibration. The paper deserves credit for explicitly flagging this limitation in Section 3, which is why the verdict should remain conditional rather than move to rejection.","tokens_in":9624,"tokens_out":12639,"duration_ms":124343,"concrete_test":"Use the observed super-Earth/sub-Neptune occurrence rate f_SE(r) from Petigura et al. (2018) as a proxy for the core formation frequency f_core(r). Replace Eq. (3) with fOR(r) = kOR f_core(r) (Mdot_p/Mdot_d) and re-run the inversion in §2.2-2.4, normalizing f_core to unity at 1 au. If the reconstructed Σg(r) and Bz(r) retain their monotonic increase and a break at r ≃ 0.1 au, the uniform-core assumption is not responsible for the central result. If the 0.1 au break shifts, disappears, or flattens, the central claim fails because the inferred disk structure is degenerate with core-formation efficiency.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central inference (Eq. 3) sets fOR(r) = kOR * Mdot_p/Mdot_d and then inverts to Σg, α, and Bz. This requires that the local formation frequency of 10-M⊕ cores is radially uniform, because any radial dependence f_core(r) enters the same equation as fOR and is degenerate with disk properties. The paper's own Section 3 flags this as 'the most critical one' and concedes that Kepler super-Earth/sub-Neptune occurrence rates increase monotonically from 1 to 10 days before flattening, so the assumption is 'justified only for longer periods.' Since the claimed r ≃ 0.1 au break sits exactly in the 1-10 day region where f_core is known to vary, the reconstructed increasing Σg profile and the r^-3 to r^-2 magnetic-field transition may simply be inherited from the assumed uniformity rather than from disk physics. kOR = 1 and Mdot_d = 1e-8 M⊙/yr also set the absolute normalization without independent calibration, so only shape claims survive even absent this degeneracy. The claim that occurrence surveys 'constrain disk quantities' therefore rests on an untested, possibly false premise.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes that the observed occurrence-rate distribution of close-in giant planets (Santerne et al. 2016) can be inverted, under the in-situ gas-accretion hypothesis, to recover the properties of the natal protoplanetary disk. The authors assume fOR(r) = kOR * Mdot_p/Mdot_d, where Mdot_p is the disk-limited gas accretion rate onto a 10 M_Earth core (Eq. 1) and Mdot_d is the steady-state disk accretion rate (Eq. 2). Solving Eqs. (1)-(4) self-consistently with viscous heating and two opacity regimes, they obtain T_d(r), Sigma_g(r), and alpha(r); they then convert alpha into B_z(r) using published ideal-MHD and nonideal-MHD alpha-beta relations. They find that Sigma_g increases with r with structure near 0.1 au, that alpha decreases with r, and that B_z can be represented by r^-3 and r^-2 power laws with a transition at roughly 0.1 au, which they associate with the giant-planet period valley. They conclude that occurrence-rate surveys can constrain inner-disk gas distributions and magnetic field geometries and provide a test of the in-situ scenario.","tokens_in":9915,"tokens_out":20847,"duration_ms":188746,"significance":"If valid, the inversion would provide a novel, observationally grounded route to inner-disk quantities (Sigma_g, alpha, B_z) that are otherwise largely unconstrained, and the explicit connection to the period valley is appealing. The paper is commendably transparent: the input data and formulas are clearly stated, the algebra is simple enough to check, and the most critical caveat (uniform core formation frequency) is acknowledged in Section 3. The strength of the paper is this clarity: a reader can immediately see that the derived Sigma_g profile \"traces the shape of fOR\" (Section 2.3) and that the B_z switch \"originates from fOR\" (Section 2.4). However, this transparency also reveals that the main results are inherited from the input occurrence-rate distribution through the assumed proportionality, and that the inversion is degenerate with the radial dependence of the formation frequency of 10 M_Earth cores. The significance is therefore conditional; the paper is best read as an illustrative inversion under a stated hypothesis rather than a demonstrated empirical reconstruction.","major_comments":[{"comment":"There is an apparent factor-of-100 normalization inconsistency between Eq. (1) and Eq. (3). Using the steady-state relation Mdot_d = 3*pi*nu*Sigma_g with nu = alpha*c_s*H_g, Eq. (1) can be rewritten as Mdot_p/Mdot_d = [0.29/(3*pi)] * (M_p/M_*)^(4/3) * (H_g/r_p)^(-4) * alpha^(-1). For the canonical parameters quoted in the text (alpha = 10^-2, H_g/r_p = 0.05, M_p = 10 M_Earth, M_* = M_Sun), this ratio equals 0.46, not 4.6*10^-3. Therefore Eq. (3), which states fOR = 0.46 kOR ... %, is too small by a factor of 100 unless kOR is correspondingly reduced. With kOR = 1 the predicted occurrence rate is about 46%, whereas the observed values in Fig. 1 are at the percent level. The absolute normalization enters alpha, Sigma_g, and B_z through Eqs. (3)-(8); in particular, the B_z values in Fig. 3 and their comparison with the dipole and large-scale field estimates in Eqs. (9) and (10) would shift by a large factor if the normalization is corrected. Please check the numerical coefficient and either correct Eq. (3) or specify the effective kOR that makes the normalization consistent with the observed rates.","section":"Section 2.1, Eqs. (1) and (3)"},{"comment":"The main inference is degenerate with the radial dependence of the formation frequency of ~10 M_Earth cores. Equation (3) should formally be fOR = kOR * f_core(r) * Mdot_p/Mdot_d, where f_core(r) is the local formation frequency of cores; any radial variation in f_core enters the equation in exactly the same way as a variation in disk properties. The manuscript itself states (Section 3) that Kepler super-Earth and sub-Neptune occurrence rates increase from 1 to 10 days and flatten only beyond 10 days, so the assumption of uniform core formation is \"justified only for longer periods.\" Since r ~ 0.1 au corresponds to orbital periods of about 10-11 days for a solar-mass star, the region that produces the claimed break in Sigma_g and B_z is precisely where the uniformity assumption is least secure. This is not merely a formal caveat: without an independent model for f_core(r), the reconstruction cannot distinguish low Sigma_g from a low core abundance inside 10 days. The abstract and Section 2.4 should be qualified accordingly, or the analysis should be repeated with an explicit f_core(r) term to show which conclusions are robust.","section":"Section 3 and Eq. (3)"},{"comment":"The paper does not propagate the observational uncertainties of the occurrence-rate measurements into the derived disk quantities. Figure 1 shows substantial error bars, and while the right panel of Fig. 2 displays error bars for alpha, no confidence intervals are given for Sigma_g or B_z, and the text does not discuss how the error bars were computed. The central claim of a \"structure\" at r ~ 0.1 au and a magnetic-field transition at the same location is therefore not yet supported against the noise in the input data. Please add a propagation of the fOR uncertainties (or a sensitivity analysis in which the bin values are varied within their error bars) and state explicitly whether the 0.1 au feature survives.","section":"Section 2.3, Fig. 2"}],"minor_comments":[{"comment":"The displayed equation contains an extra comma after W_zphi; please remove it.","section":"Section 2.4, Eq. (6)"},{"comment":"The right-panel axis label appears truncated (\"Effective \"); it should read \"Effective alpha\", and the caption should define alpha for completeness.","section":"Figure 2 caption"},{"comment":"The footnote states that the sample covers radii 4-24 R_Earth, \"equivalently 3 M_Earth - 80 M_Jupiter\"; this equivalence is dimensionally unclear and should be reworded or supported by a citation.","section":"Section 2.1, footnote 1"},{"comment":"The r-to-P conversion in the discussion of the period valley assumes a solar-mass host star, but the Santerne et al. (2016) sample includes a range of stellar masses; please state this assumption explicitly.","section":"Section 3"},{"comment":"The reference to Mayor et al. (2011) is given as an arXiv preprint; please update to the published version if one exists.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The paper is transparent and the algebra is easy to check, which is a genuine strength. My main concerns are the factor-100 normalization issue in Eq. (3) and the fact that the central 0.1 au feature lies exactly in the region where the uniform-core assumption is most questionable. These are fixable by a careful revision, but if the quantitative checks reveal that the absolute B_z profiles are no longer comparable to the stellar/large-scale field estimates, the authors should be prepared to weaken the conclusions. I do not see grounds for rejection on novelty; the inversion idea is worth publishing in revised form."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nHere's my read on Hasegawa et al. 1908.00647. The paper inverts the occurrence rate of close-in giant planets, under the in-situ gas accretion scenario, to reconstruct natal disk properties: gas surface density, effective alpha, and vertical magnetic field. The algebra is straightforward and the writing is honest.\n\nWhat's new: the explicit reconstruction of Bz from the occurrence rate is new, and the suggestion that the r^-3 to r^-2 magnetic-field transition falls at the giant-planet period valley is a striking idea. The rising Sigma_g profile was already proposed by Ogihara et al. 2015, and the authors say so. Credit where due: they use standard disk physics, two opacity regimes, and published MHD alpha relations, and they clearly state the fundamental hypothesis.\n\nThe soft spots are not hidden. The load-bearing premise is that fOR(r) is proportional to the gas accretion efficiency onto 10-M_Earth cores, with the formation frequency of such cores radially uniform. The paper itself labels this the most critical issue, and notes that Kepler data show super-Earth/sub-Neptune occurrence rising from 1 to 10 days, so the assumption is only justified at longer periods. But the period valley sits around 0.1 au, inside the region where the assumption is weakest. Since fOR is the input, and the paper says the switch in Bz profiles originates from fOR, the claimed coincidence with the period valley is inherited from the data, not an independent test. That is circular in effect: the valley is in the input occurrence rate, and it appears in the output fields.\n\nAlso, the error bars are not propagated, so the derived shapes are qualitative. The absolute normalizations depend on kOR, Mdot_d, and Mp, all set without calibration; the authors admit kOR is unknown. For a letter, these are not fatal if the claims are framed as a reconstruction rather than a prediction.\n\nWho should read this: researchers in planet formation and inner disk physics. It is a useful illustration of how occurrence-rate inversions might constrain disk structure, but it should be read with the circularity in mind.\n\nMy recommendation: engage with the paper. It deserves a serious referee who will probe the uniformity assumption and the error propagation. I would not desk-reject it, but I would insist that the abstract and conclusions stop treating the period-valley match as confirmation. Conditional acceptance, not a clean endorsement.","headline":"A transparent but ultimately circular inversion of giant-planet occurrence rates into inner-disk gas and magnetic-field profiles; the authors flag the key assumption, but it undermines the period-valley connection.","tokens_in":10422,"tokens_out":3453,"would_cite":true,"duration_ms":35055,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The distribution of close-in giant planets can be inverted, under in-situ gas accretion, to recover the gas surface density and magnetic-field structure of their natal disks.","keywords":["close-in giant planets","in-situ gas accretion","occurrence rate inversion","gas surface density","magnetic fields","period valley","minimum-mass extrasolar nebula","protoplanetary disks"],"falsifier":"A finer-binned occurrence-rate survey, or a direct measurement of inner-disk gas that showed the surface density decreasing with radius between 0.05 and 1 au, would contradict the reconstructed profile and undermine the inversion; likewise, showing that the formation frequency of ten-Earth-mass cores varies strongly within 0.1 au would remove the basis for equation (3).","tokens_in":9413,"feed_emoji":"🪐","tokens_out":8259,"duration_ms":71685,"temperature":0.7,"pith_summary":"This paper claims that exoplanet occurrence-rate surveys can be read as a map of the inner regions of planet-forming disks. Under the in-situ gas-accretion scenario, close-in giant planets form where they are observed, so the rate at which they appear at a given orbital radius should trace how efficiently roughly ten-Earth-mass cores accrete local disk gas. Inverting the observed occurrence rates this way yields a gas surface density that increases with distance from the star and has a feature near 0.1 au, in contrast to the standard minimum-mass solar nebula, which decreases outward. The same inversion reconstructs magnetic-field profiles that switch from a stellar-dipole scaling ($\\propto r^{-3}$) close to the star to a large-scale-field scaling ($\\propto r^{-2}$) farther out, with the switch at about 0.1 au, the period valley of giant exoplanets. If the claim is right, planet-occurrence statistics become a probe of the gas and magnetic environment of inner disks.","feed_headline":"Hot-Jupiter census maps the gas and magnetic fields of natal disks","feed_subtitle":"Inverting close-in giant-planet rates yields rising gas density and a magnetic switch at 0.1 au, the period valley.","key_machinery":"The load-bearing identity is equation (3), $f_{\\rm OR}(r) \\equiv k_{\\rm OR}\\,\\dot{M}_p/\\dot{M}_d$, which says the occurrence rate tracks the disk-limited gas-accretion efficiency onto planetary cores. The gas-accretion rate $\\dot{M}_p$ uses the Tanigawa-Ikoma formula; the disk accretion rate uses the steady-state relation $\\dot{M}_d = 3\\pi\\nu\\Sigma_g$ with Shakura-Sunyaev viscosity $\\nu=\\alpha c_s H_g$. A self-consistent temperature from viscous heating (equation 4) and Bell et al. opacities closes the system, so the shape of $f_{\\rm OR}$ is imprinted directly on $\\Sigma_g$, $\\alpha$, and $T_d$. The magnetic-field reconstruction then runs through two published scalings: Salvesen et al. for ideal-MHD turbulence and Bai for nonideal-MHD disk winds, both of which relate $\\alpha$ to the vertical magnetic flux $B_z$; the resulting $B_z$ profiles carry the same 0.1 au transition.","core_discovery":"The central claim is that the occurrence-rate distribution $f_{\\rm OR}(r)$ of close-in Jovian planets, taken from Kepler transit data with radial-velocity follow-up, can be inverted to recover the properties of the disks in which they formed. The inversion assumes that $f_{\\rm OR}(r)$ is proportional to the disk-limited gas-accretion efficiency onto $10\\,M_\\oplus$ cores; combining this with steady-state viscous disk accretion and self-consistent viscous heating gives the gas surface density $\\Sigma_g(r)$, effective viscosity $\\alpha(r)$, and disk temperature $T_d(r)$. The reconstructed $\\Sigma_g$ increases with radius and carries a feature at $r \\simeq 0.1$ au, unlike the $r^{-3/2}$ decline of the minimum-mass solar nebula. Converting the $\\alpha$ profile into magnetic fields through published ideal-MHD and nonideal-MHD scalings, the paper finds vertical fields that follow a stellar dipole ($B_z \\propto r^{-3}$) close to the star and large-scale fields ($B_z \\propto r^{-2}$) farther out, with the transition again at about 0.1 au, coincident with the period valley. The paper concludes that the in-situ scenario is testable against disk quantities and that occurrence data can constrain the gas distribution of the minimum-mass extrasolar nebula.","pith_inferences":["A direct test of the core-uniformity assumption would compare the occurrence rates of close-in super-Earths and sub-Neptunes, which trace where cores form, with the reconstructed gas profile; if cores are scarce inside 10 days, the recovered $\\Sigma_g$ is biased.","The same inversion applied to independent occurrence-rate data sets, such as pure radial-velocity surveys, would tell whether the rising $\\Sigma_g$ and the 0.1 au break are robust features or artifacts of the Kepler sample.","Future observations of inner-disk gas, for example line or continuum emission probing 0.05-1 au, could check directly whether the gas surface density really increases outward as predicted."],"forward_implications":["Occurrence-rate surveys of close-in giants become a direct empirical constraint on inner-disk gas surface density and temperature, effectively building a minimum-mass extrasolar nebula for $r \\lesssim 1$ au.","The predicted change from stellar-dipole to large-scale-field scaling at about 0.1 au links the hot-Jupiter period valley to a magnetic transition, implying stellar fields shape hot-Jupiter formation while large-scale fields matter for giant planets beyond 0.1 au.","The increasing surface-density profile supports earlier work showing that such a profile reproduces the observed close-in super-Earth population better than the decreasing minimum-mass solar nebula profile.","In-situ gas accretion appears viable for planets forming beyond 0.1 au; inside that radius the model needs very massive cores or very high viscosity, so other mechanisms such as high-eccentricity migration and tidal circularization may explain the innermost hot Jupiters."],"supporting_citations":[{"why":"supplies the occurrence-rate distribution of close-in giant planets that is inverted to recover disk properties.","marker":"Santerne et al. (2016)"},{"why":"provides the disk-limited gas-accretion formula used in equation (1), the core of the inversion.","marker":"Tanigawa & Ikoma (2007)"},{"why":"supplies the alpha-viscosity parametrization used to connect accretion rate and surface density.","marker":"Shakura & Sunyaev (1973)"},{"why":"gives the opacity regimes that determine the self-consistent disk temperature profiles.","marker":"Bell et al. (1997)"},{"why":"provides the ideal-MHD alpha-beta scaling used to convert the reconstructed alpha profile into magnetic field strength.","marker":"Salvesen et al. (2016)"},{"why":"provides the nonideal-MHD wind stress scaling used for the alternate magnetic-field reconstruction.","marker":"Bai (2013)"},{"why":"gives the large-scale magnetic-field estimate that the reconstructed fields are compared with.","marker":"Okuzumi et al. (2014)"},{"why":"defines the minimum-mass solar nebula profile that the reconstructed surface density is contrasted against.","marker":"Hayashi (1981)"},{"why":"establishes the in-situ gas-accretion scenario that this paper tests.","marker":"Bodenheimer et al. (2000)"}],"fun_headline_variants":["Planet census reveals rising gas density in natal disks","Close-in giants trace disk gas and magnetic switch","Inverting hot-Jupiter rates maps disk magnetic fields","Gas density rises with radius from hot-Jupiter occurrence","Natal disk fields mapped from close-in giant planet rates"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The argument stands or falls on the assumption that the occurrence rate of close-in giant planets at each radius is proportional to the gas-accretion efficiency of uniformly formed ten-Earth-mass cores, so that every bump and dip in the occurrence rate reflects gas properties rather than where cores happen to form.","fun_headline_variants_meta":{"raw":{"variants":["Planet census reveals rising gas density in natal disks","Close-in giants trace disk gas and magnetic switch","Inverting hot-Jupiter rates maps disk magnetic fields","Gas density rises with radius from hot-Jupiter occurrence","Natal disk fields mapped from close-in giant planet rates"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000235,"raw_usage":{"total_tokens":1569,"prompt_tokens":1080,"completion_tokens":489,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":696,"completion_tokens_details":{"reasoning_tokens":413}},"tokens_in":696,"tokens_out":489,"duration_ms":4690,"temperature":1.0,"reasoning_tokens":413,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T15:40:57.611818+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A finer-binned occurrence-rate survey, or a direct measurement of inner-disk gas that showed the surface density decreasing with radius between 0.05 and 1 au, would contradict the reconstructed profile and undermine the inversion; likewise, showing that the formation frequency of ten-Earth-mass cores varies strongly within 0.1 au would remove the basis for equation (3).","supporting_citations":[{"cited_title":"& Ikoma, M","cited_arxiv_id":null,"evidence_quote":"provides the disk-limited gas-accretion formula used in equation (1), the core of the inversion."},{"cited_title":"2014, ApJ, 785, 127","cited_arxiv_id":null,"evidence_quote":"gives the large-scale magnetic-field estimate that the reconstructed fields are compared with."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"establishes the in-situ gas-accretion scenario that this paper tests."}],"review_version":1}