REVIEW 4 major objections 5 minor 44 references
Gas Pressure Driven Screening Forces and Pebble Aggregation: A Pathway for Growth in Planet Formation
T0 review · 4 major / 5 minor · reviewed 2026-08-06 · deepseek-v4-flash
Pith's one-line read 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.
desk verdict A plausible screening-force mechanism for pebble growth, but the headline 10^5 yr growth timescale is not supported by the calculation as written. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
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.
What would settle it
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.
Extended reading notes
Core claim
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.
Load-bearing premise
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.
Editorial extensions
If this is right
- 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.
Reading between the lines
- 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.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
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.
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 (4)
- [§3.2.2, Eq. (14), Fig. 4] 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.
- [§2.7, Eqs. (18)-(20)] 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.
- [§2.1, Eq. (2), Caution] 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.
- [§3.2.1-3.2.2] 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.
minor comments (5)
- [Figure 3 and §3.1] 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.
- [Eq. (8)] 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.
- [Throughout] 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'.
- [§3.1 and Conclusions] 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.
- [Conclusions] 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.
Circularity Check
No significant circularity: the screening-force derivation is self-contained and its disk inputs are stated standard values rather than fitted targets.
full rationale
The paper's derivation chain is self-contained. The screening force (Eq. 2) follows from integrating the pressure anisotropy over a shadowed solid angle; the critical velocity (Eq. 4) follows from energy balance; the relative velocity (Eq. 8) is an adopted turbulent-drag prescription; and the binding probability (Eq. 18) is a product of geometric, Maxwellian, and collision-rate factors. The disk inputs (T0 = 200 K, Sigma0 = 10^4 g cm^-2, alpha = 10^-3, Z = 0.01) are stated as standard representative values and are not fitted to the predicted ring locations or the 10^5 yr growth timescale. The ALMA comparisons in Sections 1 and 4 are qualitative and do not feed back into the model. The 10^5 yr growth claim is an output of the stated collision and drag timescale calculation, not a quantity used to define any input. The main weaknesses are non-circular: Eq. (14) normalizes the projectile number density as if all dust mass were in the accreting size bin, which is inconsistent with the paper's own q_N = 3 size distribution, and Qs is unspecified for Figure 4, making the growth timescale non-reproducible as written. These are correctness, consistency, and reproducibility concerns, not equivalences between outputs and inputs. No self-citation chain is load-bearing.
Assumptions & free parameters
free parameters (10)
- Sigma0 =
10^4 g cm^-2
- T0 =
200 K (100 K in a variant)
- qT =
0.5
- alpha =
10^-3 (10^-2 in a variant)
- Z =
0.01
- rho_p =
1 g cm^-3
- Qs =
not specified (presumably 1 in Figure 4)
- qN =
3
- pebble radius range =
0.01 to 10 cm
- ya =
1.6
assumptions (8)
- domain assumption Isotropic, noninteracting thermal gas with straight-line molecular trajectories between pebbles when s < lambda.
- ad hoc to paper Perfect shadow with zero transmitted flux and a step-function cutoff at s = lambda.
- domain assumption Pebbles are smooth spheres.
- domain assumption Two-body interaction only; multi-body shielding is neglected.
- domain assumption Relative pebble velocities follow a Maxwellian distribution with dispersion vrel/sqrt(3), where vrel comes from Eq. 8.
- domain assumption Every bound pair merges via gas drag.
- ad hoc to paper Small particles with R2 < lambda are continuously available at density n_peb = Qs Z rho / M2.
- domain assumption Gas temperature and density follow the adopted Shakura-Sunyaev power-law disk with vertical stratification.
invented entities (1)
-
Screening (shadow) force in neutral gas
Cite this review
Pith. "Pith review of Gas Pressure Driven Screening Forces and Pebble Aggregation: A Pathway for Growth in Planet Formation." pith.science (2026). https://pith.science/paper/YOD54BD7
@misc{pith2026250702570,
author = {Pith},
title = {Pith review of: Gas Pressure Driven Screening Forces and Pebble Aggregation: A Pathway for Growth in Planet Formation},
year = {2026},
howpublished = {\url{https://pith.science/paper/YOD54BD7}},
note = {Machine review of arXiv:2507.02570}
}
abstract
The formation of planetesimals from cm-sized pebbles in protoplanetary disks faces significant barriers, including fragmentation and radial drift. We identify a previously unaccounted screening force, arising from mutual shielding of thermal gas particles between pebbles when their separation falls below the gas mean free path. This force facilitates pebble binding, overcoming key growth barriers under turbulent disk conditions. Unlike conventional mechanisms, screening forces operate independently of surface adhesion and complement streaming instability and pressure traps by enhancing aggregation in high-density regions. Our analysis predicts that screening interactions are most effective in the {middle disk regions ($ \sim 0.3$ to few AU),} consistent with ALMA observations (e.g., TW Hya) of enhanced dust concentrations. {Furthermore, we find that screening-induced pebble growth from centimeter to kilometer scales can occur on timescales significantly shorter than the disk lifetime ($\sim 10^5$ years). Importantly, this growth naturally terminates when particles smaller than the local gas mean free path are depleted, thereby avoiding runaway accretion.} Beyond planetary science, the screening forces have {potential} implications for high-energy astrophysics, dusty plasmas, confined particle suspensions and other relevant areas, suggesting a broader fundamental significance.
Figures
Reference graph
Works this paper leans on
-
[1]
Andrews, S. M., Wilner, D. J., Hughes, A. M., et al. 2012, ApJ, 744, 162, doi: 10.1088/0004-637X/744/2/162
-
[2]
Andrews, S. M., Wilner, D. J., Zhu, Z., et al. 2016, ApJL, 820, L40, doi: 10.3847/2041-8205/820/2/L40
-
[3]
Armitage, P. J. 2020, Astrophysics of planet formation, Second Edition
2020
-
[4]
2018, A&A, 609, L2, doi: 10.1051/0004-6361/201732034
Banzatti, A., Garufi, A., Kama, M., et al. 2018, A&A, 609, L2, doi: 10.1051/0004-6361/201732034
-
[5]
1995, A&A, 295, L1, doi: 10.48550/arXiv.astro-ph/9501050
Barge, P., & Sommeria, J. 1995, A&A, 295, L1, doi: 10.48550/arXiv.astro-ph/9501050
-
[6]
2016, SSRv, 205, 41, doi: 10.1007/s11214-016-0256-1
Birnstiel, T., Fang, M., & Johansen, A. 2016, SSRv, 205, 41, doi: 10.1007/s11214-016-0256-1
-
[7]
2012, A&A, 539, A148, doi: 10.1051/0004-6361/201118136
Birnstiel, T., Klahr, H., & Ercolano, B. 2012, A&A, 539, A148, doi: 10.1051/0004-6361/201118136
-
[8]
Birnstiel, T., Ormel, C. W., & Dullemond, C. P. 2011, A&A, 525, A11, doi: 10.1051/0004-6361/201015228
Show all 44 references
-
[9]
2015, A&A, 582, A112, doi: 10.1051/0004-6361/201526463
Bitsch, B., Lambrechts, M., & Johansen, A. 2015, A&A, 582, A112, doi: 10.1051/0004-6361/201526463
2015 doi
-
[10]
2008, ARA&A, 46, 21, doi: 10.1146/annurev.astro.46.060407.145152 Carrasco-Gonz´alez, C., Henning, T., Chandler, C
Blum, J., & Wurm, G. 2008, ARA&A, 46, 21, doi: 10.1146/annurev.astro.46.060407.145152 Carrasco-Gonz´alez, C., Henning, T., Chandler, C. J., et al. 2016, ApJL, 821, L16, doi: 10.3847/2041-8205/821/1/L16 Carrasco-Gonz´alez, C., Sierra, A., Flock, M., et al. 2019, ApJ, 883, 71, d...
2008
-
[11]
Chamberlin, T. C. 1916, JRASC, 10, 473
1916
-
[12]
Chapman, S., & Cowling, T. G. 1970, The mathematical theory of non-uniform gases. an account of the kinetic theory of viscosity, thermal conduction and diffusion in gases
1970
-
[13]
Chatterjee, S., & Tan, J. C. 2014, ApJ, 780, 53, doi: 10.1088/0004-637X/780/1/53 —. 2015, ApJL, 798, L32, doi: 10.1088/2041-8205/798/2/L32
2014 doi
-
[14]
Dohnanyi, J. S. 1969, J. Geophys. Res., 74, 2531, doi: 10.1029/JB074i010p02531
1969 doi
-
[15]
1949, Reviews of Modern Physics, 21, 367, doi: 10.1103/RevModPhys.21.367
Gamow, G. 1949, Reviews of Modern Physics, 21, 367, doi: 10.1103/RevModPhys.21.367
1949 doi
-
[16]
C., Zhu, Z., et al
Hu, X., Tan, J. C., Zhu, Z., et al. 2018, ApJ, 857, 20, doi: 10.3847/1538-4357/aaad08
2018 doi
-
[17]
C., & Chatterjee, S
Hu, X., Zhu, Zhaohuan an d Tan, J. C., & Chatterjee, S. 2016, ApJ, 816, 19, doi: 10.3847/0004-637X/816/1/19
2016 doi
-
[18]
Hogerheijde, M. R. 2011, ApJ, 727, 85, doi: 10.1088/0004-637X/727/2/85
2011 doi
-
[19]
Ignatov, A. M. 1996, Plasma Physics Reports, 22, 585
1996
-
[20]
2010, MNRAS, 404, 475, doi: 10.1111/j.1365-2966.2010.16309.x
Johansen, A., & Lacerda, P. 2010, MNRAS, 404, 475, doi: 10.1111/j.1365-2966.2010.16309.x
2010
-
[21]
2017, Annual Review of Earth and Planetary Sciences, 45, 359, doi: 10.1146/annurev-earth-063016-020226
Johansen, A., & Lambrechts, M. 2017, Annual Review of Earth and Planetary Sciences, 45, 359, doi: 10.1146/annurev-earth-063016-020226
2017 doi
-
[22]
S., Mac Low, M.-M., et al
Johansen, A., Oishi, J. S., Mac Low, M.-M., et al. 2007, Nature, 448, 1022, doi: 10.1038/nature06086
2007 doi
-
[23]
2007, ApJ, 662, 627, doi: 10.1086/516730
Johansen, A., & Youdin, A. 2007, ApJ, 662, 627, doi: 10.1086/516730
2007 doi
-
[24]
2009, ApJL, 704, L75, doi: 10.1088/0004-637X/704/2/L75
Johansen, A., Youdin, A., & Mac Low, M.-M. 2009, ApJL, 704, L75, doi: 10.1088/0004-637X/704/2/L75
2009 doi
-
[25]
2009, PhRvL, 103, 215502, doi: 10.1103/PhysRevLett.103.215502
Kelling, T., & Wurm, G. 2009, PhRvL, 103, 215502, doi: 10.1103/PhysRevLett.103.215502
2009 doi
-
[26]
2016, ApJ, 827, 110, doi: 10.3847/0004-637X/827/2/110 Le Sage, G
Kruss, M., Demirci, T., Koester, M., Kelling, T., & Wurm, G. 2016, ApJ, 827, 110, doi: 10.3847/0004-637X/827/2/110 Le Sage, G. 1782, Nouveaux M´emoires de L’Acad´emie Royale des Sciences et Belles Lettres, Berlin Le Sage, G.-L. 2022, The Le Sage theory of gravitation (DigiCat)
2016 doi
-
[27]
Lissauer, J. J. 1987, Icarus, 69, 249, doi: 10.1016/0019-1035(87)90104-7 —. 1993, ARA&A, 31, 129, doi: 10.1146/annurev.aa.31.090193.001021
1987
-
[28]
J., et al
Long, F., Pinilla, P., Herczeg, G. J., et al. 2018, ApJ, 869, 17, doi: 10.3847/1538-4357/aae8e1
2018 doi
-
[29]
2013, ApJ, 775, 17, doi: 10.1088/0004-637X/775/1/17
Lyra, W., & Lin, M.-K. 2013, ApJ, 775, 17, doi: 10.1088/0004-637X/775/1/17
2013 doi
-
[30]
1980, Progress of Theoretical Physics, 64, 544, doi: 10.1143/PTP.64.544
Mizuno, H. 1980, Progress of Theoretical Physics, 64, 544, doi: 10.1143/PTP.64.544
1980 doi
-
[31]
W., & Cuzzi, J
Ormel, C. W., & Cuzzi, J. N. 2007, A&A, 466, 413, doi: 10.1051/0004-6361:20066899
2007 doi
-
[32]
B., & Murray-Clay, R
Perets, H. B., & Murray-Clay, R. A. 2011, ApJ, 733, 56, doi: 10.1088/0004-637X/733/1/56
2011 doi
-
[33]
B., Hubickyj, O., Bodenheimer, P., et al
Pollack, J. B., Hubickyj, O., Bodenheimer, P., et al. 1996, Icarus, 124, 62, doi: 10.1006/icar.1996.0190
1996
-
[34]
I., & Sunyaev, R
Shakura, N. I., & Sunyaev, R. A. 1973, A&A, 24, 337
1973
-
[35]
1941, ApJ, 94, 232, doi: 10.1086/144328 12
Spitzer, Jr., L. 1941, ApJ, 94, 232, doi: 10.1086/144328 12
1941 doi
-
[36]
S., & Wilner, D
Testi, L., Natta, A., Shepherd, D. S., & Wilner, D. J. 2003, A&A, 403, 323, doi: 10.1051/0004-6361:20030362
2003 doi
-
[37]
Weidenschilling, S. J. 1977, MNRAS, 180, 57, doi: 10.1093/mnras/180.2.57 —. 2000, SSRv, 92, 295, doi: 10.1023/A:1005259615299
1977 doi
- [38]
-
[39]
2005, ApJL, 626, L109, doi: 10.1086/431757
Hartmann, L. 2005, ApJL, 626, L109, doi: 10.1086/431757
2005 doi
-
[40]
2021, Nature Reviews Physics, 3, 405, doi: 10.1038/s42254-021-00312-7
Wurm, G., & Teiser, J. 2021, Nature Reviews Physics, 3, 405, doi: 10.1038/s42254-021-00312-7
2021 doi
- [41]
-
[42]
N., & Lithwick, Y
Youdin, A. N., & Lithwick, Y . 2007, Icarus, 192, 588, doi: 10.1016/j.icarus.2007.07.012
2007 doi
-
[43]
J., Booth, R
Zagaria, F., Clarke, C. J., Booth, R. A., Facchini, S., & Rosotti, G. P. 2023, ApJL, 959, L15, doi: 10.3847/2041-8213/ad0c54
2023 doi
-
[44]
W., G¨uttler, C., Blum, J., & Dullemond, C
Zsom, A., Ormel, C. W., G¨uttler, C., Blum, J., & Dullemond, C. P. 2010, A&A, 513, A57, doi: 10.1051/0004-6361/200912976
2010 doi
Reviewed August 6, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.