{"id":"cbc7aa70-4837-4693-a5e9-83bb069705df","arxiv_id":"2507.02570","paper_version":1,"verdict":"REJECT","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"high","formal_verification":"none","parameter_count":10,"one_line_summary":"The paper derives an attractive screening force between pebbles in protoplanetary gas and claims it can grow centimeter pebbles to kilometer planetesimals in about 10^5 years, mainly between 0.3 and a few AU.","lead":"A new study proposes that pebbles in planet-forming disks can attract each other through a gas-pressure screening force when they are closer than the gas mean free path, enabling growth from centimeters to kilometers in about one hundred thousand years. The claim matters because it offers a surface-independent route past the fragmentation and drift barriers that block planetesimal formation.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. (14) mis-assigns the full dust mass to the accreting sub-λ size bin; with the paper's own q_N=3 size distribution, the usable projectile density drops by ~Rmax/λ, so the 10^5 yr growth time in Fig. 4 is unsupported.","rationale":"I approached the paper as a proposal of a physical idea: two pebbles separated by less than the gas mean free path can attract via pressure shadowing. The force law (Eq. 2) is at least dimensionally reasonable and the manuscript's own Ncol/Nrec consistency check at §3 is a useful sanity check, though it is not a derivation of the timescale. The author's list of caveats is honest, but none of them addresses the projectile-supply normalization issue. The strongest claim, growth from 1 cm to 10 km in ~10^5 yr, requires an iterative application of Eq. (20). The published description of that iteration is incomplete: no algorithm, no Qs, no size binning. More importantly, Eq. (14) treats all dust mass as if it were in the single accreting size M2, while §2.4's power-law distribution means only a small fraction of dust mass lies below λ. This overestimates collision rates by roughly Rmax/λ (~30) in the fiducial model and, combined with the unstated Qs, makes the 10^5 yr number an artifact of an inconsistent normalization rather than a derived prediction. A rerun with the differential size distribution is a decisive, low-cost check. Because this concern independently undermines the central quantitative claim, I agree with the REJECT verdict and recommend no change.","tokens_in":14781,"tokens_out":16752,"duration_ms":196204,"concrete_test":"Recompute the §3.2.2 growth track from first principles: for each survivor radius R1, take n_peb from the differential distribution dN/dR = A R^{-3} (A fixed by ∫ M(R)dN/dR = Qs Z ρ over Rmin=0.01 cm to Rmax=10 cm), set the collisional cross-section π(R1+R2)^2, integrate Eq. (19) over R2 from Rmin to λ(r), update R1 -> (R1^3+R2^3)^{1/3} per successful merger, and sum Tg. Vary Qs from 1 to 100 and report the time to reach 10 km at r=1 AU. If the time exceeds 1 Myr at Qs=1 (or requires Qs ≳ 30 to reach 10^5 yr), the headline growth-time claim is unsupported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative claim is the §3.2.2 growth track: 1 cm to 10 km in ~10^5 yr (Fig. 4). The track uses Tg = tcoll + tdrag (Eq. 20) with tcoll = 1/(n_peb σ_coll vrel) (Eq. 19). Eq. (14) sets n_peb = Qs Z ρ / M2, which places the entire dust mass density Zρ into monodisperse particles of mass M2. But §2.4 already adopts a Dohnanyi number distribution N(R) ∝ R^{-3} over 0.01–10 cm. For a body that accretes only projectiles with R2 < λ (the termination condition Eq. 23), the dust mass in the usable bins is only a fraction λ/Rmax ≈ 0.03 of Zρ at r=1 AU for the top-left Fig. 3 parameters (λ≈0.28 cm, Rmax=10 cm). Using the differential size distribution instead of Eq. (14) lowers the growth rate dR1/dt by this factor (≈30×), pushing the 10^5 yr timescale to several Myr even before depletion of the sub-λ reservoir and with unspecified Qs in Fig. 4. The paper provides no code, no Qs value, and no step-by-step iteration for Fig. 4, so the published timescale cannot be reproduced from the manuscript as written. This is the load-bearing weak point for the abstract's headline claim.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper proposes a new attractive force between pebbles in a protoplanetary disk, arising from mutual shadowing of thermal gas molecules when the surface-to-surface separation is smaller than the gas mean free path. The authors derive an inverse-square force (Eq. 2), a critical velocity for bound encounters (Eq. 4), and a binding probability that combines collision rate, Maxwellian velocity distribution, and geometric shadowing (Eq. 18). Using a power-law disk model, they predict the binding probability peaks near 0.7-0.8 AU, and they present a growth track in Figure 4 in which a 1 cm pebble grows to 10 km in about 10^5 years. They argue the process terminates naturally when particles smaller than the mean free path are exhausted, and they discuss comparisons with ALMA ring observations and possible broader applications.","tokens_in":15115,"tokens_out":11915,"duration_ms":142423,"significance":"If the screening force exists with the strength claimed and the 10^5 yr growth timescale is correct, the paper would open a new, surface-adhesion-independent pathway for planetesimal formation and would provide a concrete mechanism for enhanced dust aggregation in the inner-to-mid disk. The analytical derivation is explicit and the paper is candid about its limitations in the Caution section. However, the central quantitative claim is not reproducible from the manuscript: the key number density in Eq. (14) is inconsistent with the adopted size distribution, Qs is unspecified for Figure 4, the iterative growth calculation is not described, and the depletion that is invoked as a termination mechanism is absent from the growth model. The physical force derivation also omits a treatment of molecular re-emission from the pebble surfaces. These are not cosmetic issues; they affect the headline timescale and the existence of the mechanism.","major_comments":[{"comment":"The central growth timescale is not supported by the equations as written. Equation (14) sets n_peb = Qs Z ρ / M2, placing the entire dust mass density Zρ into monodisperse particles of the projectile mass M2. This contradicts the Dohnanyi size distribution N(R) ∝ R^{-3} adopted in §2.4 and the termination condition R2 < λ (Eq. 23). For the Figure 3 top-left parameters at 1 AU, λ ≈ 0.28 cm while the size distribution extends to Rmax = 10 cm, so only a fraction λ/Rmax ≈ 0.03 of the dust mass is in accretable projectiles. Using the differential size distribution raises t_coll (Eq. 19) by roughly a factor of 30, changing the 10^5 yr track to several Myr for Qs = 1, i.e., comparable to or longer than the disk lifetime. In addition, Figure 4 does not state the value of Qs, the projectile size/time step, or the iteration procedure, so the plotted curve cannot be reproduced from the manuscript. Because the abstract's headline claim is the 10^5 yr growth, this is a load-bearing deficiency.","section":"§3.2.2, Eq. (14), Fig. 4"},{"comment":"The growth loop ignores the binding probability PB in the effective collision rate. Equation (20) uses t_coll = 1/(n_peb σ_coll v_rel), but if only a fraction PB of encounters results in a bound state, the average time per successful merger is t_coll/PB, not t_coll. The text asserts PB ≈ 1 in the middle disk regions, but the supporting Figure 3 is ambiguous: the caption says the solid line is for the smallest pair R1 = R2 = 10^-2 cm, while §3.1 says it is for R1 = 10 cm with R2 = 10^-2 cm. These two cases have very different v_crit (Eq. 7), so the figure does not establish PB ≈ 1 for the size ratios relevant to the Figure 4 growth track. If PB is 10^-2 or lower for part of the growth sequence, the integrated timescale increases correspondingly.","section":"§2.7, Eqs. (18)-(20)"},{"comment":"Equation (2) is derived by integrating the ambient isotropic pressure over the unshadowed portion of the smaller sphere, treating the larger sphere as a perfect geometric blocker of the momentum flux. Real pebbles re-emit or reflect gas molecules, and molecules scattered from the facing surface can strike the other pebble, partially filling the shadow and reducing the net force. The magnitude of this effect depends on the accommodation coefficient and the scattering kernel, neither of which is specified. The Caution paragraph acknowledges oblique flux and a weak repulsive component but gives no quantitative estimate. Because the existence and sign of the screening force is the foundation of the paper, a free-molecular-flow calculation or particle simulation is needed before the derived timescales can be accepted.","section":"§2.1, Eq. (2), Caution"},{"comment":"The growth model and the termination mechanism are mutually inconsistent. Section 3.2.1 states that growth halts when the local sub-λ particles are depleted, and that this avoids runaway accretion. Yet §3.2.2 and Figure 4 assume n_peb = Qs Z ρ / M2 (Eq. 14) is maintained throughout the entire merger sequence from 1 cm to 10 km; no depletion term or time-dependent size distribution is included. The 10^5 yr curve therefore uses an upper bound on the projectile supply and does not account for the finite reservoir of R2 < λ particles. A self-consistent calculation must integrate over the decreasing number of available projectiles, which will lengthen the growth time and is particularly important for the kilometer-size tail.","section":"§3.2.1-3.2.2"}],"minor_comments":[{"comment":"The caption and the main text disagree on the solid-line curve: the caption states R1 = R2 = 10^-2 cm, while the text says R1 = 10 cm and R2 = 10^-2 cm; please reconcile because the two configurations have very different binding probabilities.","section":"Figure 3 and §3.1"},{"comment":"The parameter ya = 1.6 is introduced without a definition or a reference; the meaning of this parameter should be stated explicitly, and its relation to the Ormel & Cuzzi (2007) formulation should be clarified.","section":"Eq. (8)"},{"comment":"There are several typographical errors, including 'Lagragian' in §2.2 and 'Stocks regime' in §3.2.2; these should be corrected to 'Lagrangian' and 'Stokes regime'.","section":"Throughout"},{"comment":"The favorable region for screening is described as both '0.3 to a few AU' and '0.3 to 10 AU' at different points; please use one consistent radial range and explain whether the outer boundary is determined by the 90% of peak criterion or by the disk model.","section":"§3.1 and Conclusions"},{"comment":"The conclusion that the screening force 'naturally emerges from the established physics of gas-rich disks' is stronger than the derivation supports, since the derivation assumes a step-function cutoff and perfect shadowing; a more cautious statement of assumptions would be appropriate.","section":"Conclusions"}],"recommendation":"reject","confidential_remarks":"The manuscript's central quantitative claim (1 cm to 10 km in ~10^5 yr) depends on an inconsistent use of the size distribution in Eq. (14) and on an unspecified Qs in Figure 4, so the result is not reproducible as written. The physical mechanism also lacks a treatment of molecular re-emission from pebble surfaces, which could substantially alter or cancel the force. These are load-bearing issues rather than presentation issues. If the authors can provide a size-resolved coagulation calculation, specify all parameters, and support the force law with a kinetic-theory calculation, the idea might be worth revisiting, but the current manuscript does not meet the standard for publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe paper's physical idea is worth taking seriously: mutual shadowing of thermal gas can create an attractive force between pebbles separated by less than a mean free path. That is a known kinetic effect, but the application to neutral gas in protoplanetary disks is genuinely new, and the force law (Eq. 2) is internally consistent for a static pair under perfect shadowing. The paper is also honest about its ancestry, citing Le Sage, Spitzer, and Ignatov, and it lists a set of caveats about oblique fluxes and shape dependence.\n\nThe soft spot is exactly where the reader and stress-test put it. The headline claim—growth from 1 cm to 10 km in ~10^5 yr (Fig. 4)—is not backed by a reproducible calculation. Eq. (14) sets the projectile number density to Qs Z ρ / M2, which assigns the entire dust mass to monodisperse particles of mass M2. But §2.4 adopts a Dohnanyi-like distribution N(R) ∝ R^{-3} over 0.01–10 cm. The termination condition R2 < λ means a growing pebble can only use the sub-λ tail; at r=1 AU with the paper's own parameters, that tail holds only ~λ/Rmax ≈ 3% of the dust mass. Using the differential size distribution cuts the accretion rate by roughly 30x, which pushes the growth time from 10^5 yr to a few Myr—even before worrying about depletion and with Qs left unspecified in Fig. 4. The paper provides no code, no step-by-step iteration, and no error bars, so the timescale cannot be reproduced from the manuscript.\n\nThe rest of the paper is qualitative. The claim that binding is most effective in the middle disk is a consequence of the chosen disk model, not an independent prediction, and the ALMA comparison is suggestive rather than quantitative. None of this kills the underlying idea: the screening force itself may be real and testable with particle-in-cell simulations, and the paper says so explicitly.\n\nWho this is for: researchers working on planetesimal formation who want a catalogue of possible binding mechanisms. It deserves a serious referee because the central mechanism is plausible and the application is new, but the quantitative growth track needs to be redone with a proper size-dependent collision kernel. I would send it out, with a strong request for revision.\n\nRecommendation: accept for peer review; reject in current form.","headline":"A plausible screening-force mechanism for pebble growth, but the headline 10^5 yr growth timescale is not supported by the calculation as written.","tokens_in":15703,"tokens_out":3191,"would_cite":false,"duration_ms":33008,"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 paper argues that gas-pressure screening between close pebbles can bind them into aggregates and grow centimeter pebbles to kilometer planetesimals in roughly 10^5 years.","keywords":["screening force","pebble aggregation","planetesimal formation","protoplanetary disks","gas mean free path","binding probability","dust growth","streaming instability"],"falsifier":"Recalculate the Section 3.2.2 growth track while explicitly counting the number of sub-mean-free-path particles consumed at each step; if the local reservoir is exhausted before the pebble reaches 1 km, the claimed $10^5$-year growth to 10 km is refuted.","tokens_in":14476,"feed_emoji":"🪐","tokens_out":9632,"duration_ms":105683,"temperature":0.7,"pith_summary":"The paper argues that two pebbles in a protoplanetary disk attract each other when their separation is smaller than the gas mean free path: each pebble shields the other from thermal gas molecules, creating a pressure imbalance that pulls them together. This screening force is proposed as a binding mechanism that does not depend on sticky surfaces, so it could bypass the fragmentation and bouncing barriers that normally block growth beyond centimeter sizes. If the argument is right, pebbles can grow from 1 cm to 10 km in about $10^5$ years in the middle disk region, roughly 0.3 to a few AU, faster than the disk lifetime, and the process stops naturally when particles smaller than the mean free path are used up. The claim matters because it would give planetesimal formation a route that operates under turbulent conditions and complements streaming instability and pressure traps, with a peak location that matches ALMA dust concentrations.","feed_headline":"Gas shadowing can grow pebbles to 10 km in 100,000 years","feed_subtitle":"A pressure imbalance between close pebbles offers a non-sticky route past the usual growth barriers.","key_machinery":"The central object is the screening force, an attractive force between two bodies in a gas that arises from mutual shadowing of thermal molecular flux when their separation falls below the mean free path $\\lambda$. The load-bearing identity is the force law $F = -\\pi n k_B T R_1^2 R_2^2 / x^2$, which gives an inverse-square potential; equating the reduced-mass kinetic energy of an encounter with the work integral of this force defines the critical binding velocity $v_{\\mathrm{crit}}$, and a Maxwellian distribution of turbulent relative velocities converts that threshold into a binding probability. Gas drag then sets the merging timescale, and the geometric condition $R_2 < \\lambda$, derived both from the shadow length and from the collision-geometry probability $P_2$, terminates the growth sequence.","core_discovery":"On the paper's own terms, the central discovery is a previously unaccounted attractive force between spherical pebbles in a gas. When the surface-to-surface separation $s$ satisfies $s \\le \\lambda$, where $\\lambda$ is the local gas mean free path, the larger sphere blocks the thermal momentum flux reaching the smaller one over a solid angle, and the resulting anisotropic pressure gives $F = -\\pi n k_B T R_1^2 R_2^2 / x^2$, an inverse-square attraction that strengthens with gas density and temperature and vanishes for point particles. Comparing the relative kinetic energy of an encounter with the work available from this force yields a critical velocity $v_{\\mathrm{crit}}$; encounters slower than $v_{\\mathrm{crit}}$ form bound pairs, gas drag then dissipates their relative motion and merges them, and iterative merging lets a 1 cm pebble at 1 AU grow to 10 km in about $10^5$ years. Growth is self-limiting because the shadow geometry requires the accreted projectile to satisfy $R_2 < \\lambda$, so when the local reservoir of sub-mean-free-path particles is exhausted the process halts instead of running away.","pith_inferences":["The paper leaves implicit that the screening force should leave a measurable signature in grain-size distributions: a depletion of particles just above the local mean free path and an accumulation below it, which multi-wavelength ALMA and polarization maps could in principle test.","A natural testable extension is a laboratory experiment at the same ratio of particle separation to gas mean free path, where two suspended spheres should show an attractive force scaling with gas density and inverse separation squared even though disk-like conditions cannot be reproduced directly.","A direct extension not developed in the paper is a coagulation model with a screening collision kernel; such a model would show whether the $10^5$-year track survives when particle consumption and size-dependent resupply are tracked explicitly."],"forward_implications":["Planetesimal formation would no longer need surface adhesion: a purely pressure-based attraction can bind pebbles, with van der Waals forces only assisting at the final contact.","The most favorable window is the middle disk, roughly 0.3 to a few AU, where binding probabilities approach unity (peak near 0.7–0.8 AU for the fiducial disk), offering a physical explanation for enhanced dust concentrations seen by ALMA in such regions.","A 1 cm pebble at 1 AU can reach 10 km in about $10^5$ years, well below the typical disk lifetime, making the pathway dynamically competitive.","Growth stops when particles smaller than the local mean free path are depleted, so the mechanism cannot produce runaway accretion and instead allows several aggregates to grow in parallel.","The mechanism survives turbulence up to $\\alpha = 10^{-2}$ and becomes more effective as streaming instability or pressure traps raise the local pebble density, so it reinforces existing planetesimal-formation routes rather than replacing them."],"supporting_citations":[{"why":"Supplies the turbulent relative-velocity model (Equation 8) that sets the velocity dispersion entering the binding probability.","marker":"Ormel & Cuzzi 2007"},{"why":"Establishes the radial-drift barrier for centimeter pebbles that the screening mechanism is meant to overcome.","marker":"Weidenschilling 1977"},{"why":"Provides the review of fragmentation and bouncing barriers that motivates a non-adhesive binding route.","marker":"Blum & Wurm 2008"},{"why":"Introduces streaming instability, the density-enhancement context that the paper's Qs factor parameterizes and that screening complements.","marker":"Youdin & Goodman 2005"},{"why":"Gives the pebble number density formula n_peb = Qs Z rho / M2 used for collision rates.","marker":"Bitsch et al. 2015"},{"why":"Supplies the gas-drag force expression used to compute the merging timescale of bound pairs.","marker":"Perets & Murray-Clay 2011"},{"why":"Provides the disk model conventions, including sound speed, mean free path, stopping time, and surface density, used throughout.","marker":"Armitage 2020"},{"why":"Supplies the collisional size distribution power-law exponent q_N = 3 used to sample pebble pairs.","marker":"Dohnanyi 1969"}],"fun_headline_variants":["Gas shadowing attracts pebbles, enabling 10-km growth","Gas screening binds pebbles, bypassing stickiness to make planetesimals","Pebble growth via gas shadowing: 10 km in 100,000 years","Gas pressure screening forces drive pebble aggregation for planet formation","Self-limiting gas shadowing grows pebbles to kilometer scale"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The $10^5$-year growth estimate assumes that small particles below the gas mean free path remain available at the local dust density throughout the sequence and that every bound pair merges via gas drag, so the growth track never waits for resupply.","fun_headline_variants_meta":{"raw":{"variants":["Gas shadowing attracts pebbles, enabling 10-km growth","Gas screening binds pebbles, bypassing stickiness to make planetesimals","Pebble growth via gas shadowing: 10 km in 100,000 years","Gas pressure screening forces drive pebble aggregation for planet formation","Self-limiting gas shadowing grows pebbles to kilometer scale"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000665,"raw_usage":{"total_tokens":3069,"prompt_tokens":1015,"completion_tokens":2054,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":631,"completion_tokens_details":{"reasoning_tokens":1957}},"tokens_in":631,"tokens_out":2054,"duration_ms":20353,"temperature":1.0,"reasoning_tokens":1957,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T20:26:56.731236+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recalculate the Section 3.2.2 growth track while explicitly counting the number of sub-mean-free-path particles consumed at each step; if the local reservoir is exhausted before the pebble reaches 1 km, the claimed $10^5$-year growth to 10 km is refuted.","supporting_citations":[],"review_version":1}