Pith. sign in

REVIEW 3 major objections 5 minor 93 references

This paper claims that gas-phase atoms of every abundant non-inert element stick to amorphous carbon dust with probability at least 0.2 across all interstellar gas and dust temperatures, making gas-phase accretion a fast enough route to exp

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

T0 review · deepseek-v4-flash

2026-08-02 02:40 UTC pith:R4FDU6PZ

load-bearing objection A credible new simulation of dust sticking coefficients, but the central high-sticking result for C and O leans on an adopted conversion factor from the authors' own prior work that needs independent verification. the 3 major comments →

arxiv 2607.14237 v1 pith:R4FDU6PZ submitted 2026-07-15 astro-ph.GA

Efficient Interstellar Grain Growth from High Sticking Coefficients on Amorphous Carbon Dust

Clarke J. Esmerian , Duncan Bossion , Francois Dulieu , Saoud Baouche , Alexey Potapov , W. M. C. Sameera , Tom J. L. C. Bakx , Susanne Aalto
show 3 more authors
Kirsten K. Knudsen Gunnar Nyman Wouter Vlemmings
This is my paper
classification astro-ph.GA
keywords interstellar dustgrain growthsticking coefficientamorphous carbongas-phase accretionmolecular dynamicsinterstellar mediumdust evolution
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper aims to remove a major uncertainty in how interstellar dust grows: whether atoms from the surrounding gas actually stick to grain surfaces often enough to add mass on astrophysical timescales. Using reactive molecular dynamics simulations of an amorphous carbon surface, the authors compute sticking coefficients for the cosmically abundant elements H, C, O, Al, Si, S, Fe, and Ni over gas temperatures of 3–10^4 K and dust temperatures of 3–300 K. They find that the total chemisorption sticking coefficient is above 0.2 for every non-inert element, and above 0.5 for most metals, across the whole parameter range. Putting these coefficients into the standard accretion formula gives growth timescales for a 50 Å carbon grain of under 100 million years in all phases of the interstellar medium, and as short as 10^3–10^4 years in dense gas. If correct, gas-phase accretion is efficient enough to be a central process in the cosmic dust life cycle, not a minor correction to stellar dust production.

Core claim

The central discovery is that the total chemisorption sticking coefficient – the Maxwell-Boltzmann-averaged probability that an incoming atom eventually forms a covalent bond with the grain – is predicted to exceed 0.2 for all non-inert elements considered, at every relevant gas and grain temperature. For hydrogen, carbon, and oxygen, direct chemisorption is supplemented by a physisorption-then-conversion pathway: atoms first held by weak van der Waals forces later convert to covalent bonds on longer timescales, using conversion probabilities adopted from the authors' earlier simulation work. For aluminum, silicon, sulfur, iron, and nickel, the direct chemisorption sticking coefficient is al

What carries the argument

The carrying object is the sticking coefficient S_i(T_g, T_D), defined as the velocity-distribution-averaged probability that species i hitting a grain at gas temperature T_g and dust temperature T_D sticks. It is computed with reactive force-field molecular dynamics on five amorphous carbon surfaces, using 50 independent collision trajectories for each element and temperature pair; collisions are classified at the end of the simulation as chemisorbed (covalent bond formed), physisorbed (van der Waals bound), or not stuck. The growth timescale τ = (√(2π)/3) a_D ρ_D / [n_g √(k_B T_g) Σ_i S_i Z_i/(μ_i √m_i)] converts the coefficients into astrophysical predictions, assuming compact spherical g

Load-bearing premise

The load-bearing premise is that a loosely held (physisorbed) carbon or oxygen atom on an amorphous carbon surface converts to a strongly bonded (chemisorbed) state at the rate assumed from the authors' earlier simulations; if that conversion rate is too high, the claimed fast growth in cold interstellar gas collapses.

What would settle it

A laboratory measurement of the sticking probability of neutral carbon atoms (or oxygen atoms) onto an amorphous carbon surface at gas temperatures around 10–100 K and dust temperatures around 10 K would settle the central claim. If the total chemisorption sticking coefficient comes out below 0.2 in that regime, the growth timescales in the paper's Figure 3 are too short and the conclusion that accretion is efficient in cold phases fails. Alternatively, a converged quantum-accurate molecular dynamics calculation of C and O on a-C, without the adopted conversion probability, would test the same

Watch this falsifier — get emailed when new claim-graph text bears on it.

If this is right

  • A 50 Å carbon grain doubles in mass in under 100 million years in every interstellar phase, and in roughly 10^3–10^4 years in dense molecular gas.
  • Even a 0.1 μm grain grows within a billion years in most phases, meaning accretion competes with stellar sources even for larger carbonaceous grains.
  • Carbon grains should naturally accumulate O, Fe, Si, S, Al, and Ni, so realistic dust models must allow mixed/composite carbon-rich compositions rather than only separate carbon and silicate populations.
  • The supplied Bernstein-polynomial fits give dust evolution models a consistent, temperature-dependent prescription for sticking coefficients that was previously unavailable for most elements.
  • For most elements, sticking coefficients rise again at gas temperatures above ~1000 K, implying that accretion is not entirely shut off in warm or shocked gas.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Because the predicted sticking coefficients barely depend on dust temperature, galaxy and cloud simulations could treat S_i as a function of gas temperature alone, which would make dust-growth subgrid models substantially cheaper to evaluate.
  • If similar high sticking holds on silicate and ice surfaces, gas-phase accretion could be a universal grain-growth channel, which would strengthen the case that dust masses in high-redshift galaxies are built by accretion rather than only by supernovae and AGB stars.
  • The predicted trace inclusions of iron and silicon in carbon grains yield a concrete observational target: high-sensitivity mid-infrared spectra of dense clouds may reveal spectral signatures of these composite grains, which current two-population dust models would not expect.
  • A direct quantum-accurate simulation of C and O atoms on amorphous carbon at low gas temperature would test the adopted physisorption-to-chemisorption conversion probability; until that is done, the low-temperature branch of the growth timescale remains the least certain part of the prediction.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

3 major / 5 minor

Summary. This paper uses ReaxFF molecular dynamics to compute sticking coefficients for H, C, O, Al, Si, S, Fe, and Ni on amorphous carbon surfaces as functions of gas temperature (3–10^4 K) and dust temperature (3–300 K). The authors report total chemisorption sticking coefficients above ~0.2 for all non-inert elements at all studied temperatures, compare their results with earlier Tg=TD calculations and with new laboratory H-recombination lower limits, and use Eq. (1) to estimate grain-growth timescales for a 50 Å carbon grain. They find growth timescales broadly below 100 Myr in all ISM phases, with dense-phase timescales as short as 10^3–10^4 yr. The paper discusses three caveats: ice mantles, Coulomb repulsion in diffuse gas, and the fact that simulations use bare a-C rather than hydrogenated a-C:H.

Significance. If the central claim holds, the paper materially strengthens the case that gas-phase accretion is a dominant dust-growth channel and provides concrete, directly usable sticking-coefficient functions for dust evolution models. Strengths include a transparent simulation protocol, provision of polynomial fitting coefficients (Table 2), explicit discussion of caveats, and a new experimental H constraint that is consistent with the simulations. However, the headline 'S>0.2 for all non-inert elements at all relevant temperatures' is not equally supported for C and O at low gas temperatures, where the estimate depends on an adopted physisorption-to-chemisorption probability from the authors' prior work rather than on direct simulation. The H laboratory data validate only the H channel, not the C/O conversion. The paper is likely publishable after the low-temperature C/O claim is either quantified with the adopted conversion values and uncertainties or supported by additional simulation/experiment.

major comments (3)
  1. [Fig. 1; Table 2; Methods A] The central claim that total chemisorption S>0.2 for C and O at low Tg is not directly supported by the raw MD results. Table 2 gives C chem. β0=1.36×10^-9 and O chem. β0=0.237, while C and O est. tot. chem. β0 are 0.399 and 0.496. The gap for C comes from multiplying the large measured physisorption probability by a physisorption-to-chemisorption probability 'adopted from our previous work [9]' (Methods A; Fig. 1 caption). The manuscript gives neither the numerical values nor the uncertainties of these conversion probabilities, nor evidence that values derived at Tg=TD in [9] apply at all (Tg,TD). Because C and O dominate the mass-weighted sum in Eq. (1) (Z_C/μ_C^(3/2) and Z_O/μ_O^(3/2) exceed the heavier metals by orders of magnitude), a factor-of-2–5 overestimate would move τ in the CNM or dense H2 phases from below 100 Myr to ≳1 Gyr, undermining the abstract's efficiency claim. Pleas
  2. [Fig. 2; Methods B] The experimental validation is limited to H recombination on porous carbonaceous dust. The measured lower limits (R=0.25±0.10 at 100–300 K; R=0.56±0.07 at 10 K) are consistent with the H simulations, but they do not constrain the C and O physisorption-to-chemisorption conversion probabilities that carry the low-temperature central claim. The statement that consistency with laboratory data gives confidence in the numerical results should be explicitly restricted to H; the C/O extrapolation remains uncalibrated by experiment. A direct laboratory measurement of C or O sticking on carbonaceous surfaces, or at least a clear acknowledgment that such data are absent, is needed before the 'all non-inert elements' claim can be accepted at face value.
  3. [Fig. 1 caption] The caption states that the chemisorption values shown are 'lower limits' because some physisorbed atoms may chemisorb on longer timescales. This is logically in tension with the star symbols, which add an adopted conversion probability and produce total chemisorption estimates above the direct chemisorption points. Direct chemisorption counts at 108.75 ps are indeed lower limits if all physisorbed atoms eventually chemisorb, but the adopted conversion probability is not a lower limit; it is a model-dependent extrapolation. Please rephrase the caption to distinguish measured direct chemisorption, measured physisorption, and estimated asymptotic total chemisorption, and explicitly state the assumptions entering the conversion.
minor comments (5)
  1. [Methods A] The term 'molecular guns' should be defined or replaced with a more standard description, e.g., 'projectile atoms are launched from randomized positions and angles'.
  2. [Third caveat paragraph] Typo: 'extimated' should be 'estimated'.
  3. [Methods D] The phrase 'calculated in Si,log10(Tg) space' is awkward; should read 'in (S_i, log10 Tg) space' or similar.
  4. [Table 2] A column listing the adopted physisorption-to-chemisorption conversion probabilities for H, C, and O, or a direct reference to the specific values in [9], would greatly aid reproducibility and make the sensitivity of the central claim transparent.
  5. [Fig. 3 caption] The 'maximal duty-cycle of de-hydrogenation' factor of 10^2 is introduced in the caption but defined only in the text; please make the caption self-contained.

Circularity Check

1 steps flagged

C and O 'total chemisorption' >0.2 rests on conversion probabilities adopted from the same authors' previous ReaxFF paper, not on the present MD outputs.

specific steps
  1. self citation load bearing [Figure 1 caption / main text, paragraph after Fig. 1]
    "Stars indicate estimated asymptotic total chemisorption values taking the physisorption-to-chemisorption probabilities for H, C, and O from [9]."

    The paper's central prediction that 'the total chemisorption sticking coefficient is predicted to be greater than 0.2 for all gas and dust temperatures' is, for C and partly O, not a raw MD output. Table 2 gives C chem. beta0 = 1.36e-9 while C est. tot. chem. beta0 = 0.399, and O chem. beta0 = 0.237 vs. O est. tot. chem. beta0 = 0.496. The low-temperature total-chemisorption values are obtained by multiplying the physisorption signal (C phys beta0 = 0.805) by conversion probabilities adopted from the same authors' previous ReaxFF paper [9]. The manuscript supplies neither the conversion values nor their uncertainty, and the only laboratory constraint in the paper is for H. Thus the 'S >= 0.2' statement for C at T_g <~ 100 K, and roughly half the O signal, reduces to a self-cited parameter

full rationale

The derivation is mostly self-contained: the MD collision counts produce raw chemisorption and physisorption coefficients (Table 2), the Bernstein-polynomial fits (Methods D) are interpolation of those outputs, and Eq. (1) is the standard definition of the growth timescale. The laboratory H recombination measurements are an external, independent constraint, albeit only for H. The circularity is localized to the 'total chemisorption' values for H, C, and O. The Fig. 1 caption and main text state that these are obtained by taking physisorption-to-chemisorption probabilities 'from [9]', a previous ReaxFF study by overlapping authors. For C the raw chemisorption beta0 is 1.36e-9, so the predicted low-temperature S_C ~ 0.4 is almost entirely the physisorption signal (0.805) multiplied by a self-cited conversion probability; for O roughly half of the 0.496 total comes from the same conversion. The manuscript does not quote the conversion probabilities or their uncertainty, and the only laboratory check is for H. Consequently the headline claim that 'the total chemisorption sticking coefficient is predicted to be greater than 0.2 for all gas and dust temperatures' rests, for C and partly O, on a same-author citation rather than on the present simulations or an external benchmark. Because C and O dominate the mass-weighted sum in Eq. (1), this also inflates the grain-growth efficiency claim. No other load-bearing step exhibits circularity.

Axiom & Free-Parameter Ledger

4 free parameters · 5 axioms · 0 invented entities

The central claim rests on: (1) the ReaxFF potential's fidelity for all eight elements on a-C, (2) an adopted physisorption-to-chemisorption conversion probability for H, C, O from the authors' previous work, and (3) ISM-phase assumptions (solar abundances, neutral atoms, compact spheres). No new physical entities are introduced.

free parameters (4)
  • physisorption-to-chemisorption conversion probability p_phys→chem (H, C, O) = Not stated numerically; adopted from Bossion et al. 2024 [9]
    Used to convert measured physisorption into the 'estimated total chemisorption' values shown as stars in Fig. 1; for C at Tg<1000 K raw chemisorption is ~0, so the S>0.2 claim for carbon rests on this parameter. No independent value or uncertainty is given in this paper.
  • Bernstein polynomial coefficients β_k (Table 2) = Listed in Table 2 for each element and sticking type
    Constrained least-squares fit to the simulated S_i(Tg,TD) data; used to evaluate S_i at arbitrary temperatures in Eq. 1. They are interpolants of the simulations, not physically measured constants, and carry no error bars into the τ estimates.
  • Selection of 5 a-C structures and 50 collision trajectories = None
    The 5 structures out of Deringer et al. 512-atom set and 10 atoms per structure were chosen without a stated selection criterion; the average sticking coefficient is sensitive to this finite sample.
  • Maximal de-hydrogenation duty-cycle factor = 100
    Used to rescale τ in Fig. 3 (dot-dashed line) to account for repeated HII-region de-hydrogenation cycles. The factor of 100 is an order-of-magnitude assumption ('Even if this duty cycle extends the grain-growth timescale by a factor of 10^2').
axioms (5)
  • domain assumption ReaxFF force fields (CHO, AlCHO, SiC, HCONSB, CHFe, NiCH) accurately describe atom-surface chemisorption barriers and energy transfer for the 8 elements on amorphous carbon.
    The sticking coefficients are entirely determined by these empirical potentials; no direct validation against DFT or experiment is provided for most elements. Invoked throughout Methods A/Table 1.
  • domain assumption Classical nuclear dynamics suffices; quantum tunneling and zero-point effects do not change sticking probabilities.
    MD propagates nuclei classically, standard for these systems; no assessment of tunneling for H at low temperatures, which could raise or lower sticking. Methods A.
  • domain assumption A Gaussian distribution of initial velocities approximates the velocity-weighted Maxwell-Boltzmann distribution well enough to define S_i(Tg,TD).
    Methods A: 'chosen at random from a Gaussian distribution matching as closely as possible a velocity-weighted Maxwell-Boltzmann distribution'; the mismatch is systematic and unquantified.
  • ad hoc to paper The pure a-C surface represents interstellar carbonaceous grains; hydrogenation only slows, not prevents, growth.
    The paper's own third caveat: a-C:H may reduce sticking; they argue H binding is weaker and de-hydrogenation cycles occur, and apply a factor-100 duty-cycle correction, which is an order-of-magnitude assumption rather than a derived result.
  • domain assumption All non-H, non-noble elements in the gas are assumed to be neutral atoms with solar abundances when computing growth timescales; no molecular depletion (CO, etc.) or ionization is included.
    Eq. 1 sums atomic species with solar mass fractions; in dense ISM, C and O are largely in CO and ices; ionization in diffuse gas changes Coulomb interactions. The paper acknowledges ice and Coulomb caveats but the τ map in Fig. 3 uses atomic neutral sticking.

pith-pipeline@v1.3.0-alltime-deepseek · 16788 in / 16587 out tokens · 169903 ms · 2026-08-02T02:40:40.012884+00:00 · methodology

0 comments
read the original abstract

Cosmic dust is the solid phase of the interstellar medium (ISM), classically assumed to be composed of carbonaceous and silicate grains with size distributions spanning $\sim 5~\AA$ to $\sim 1~\mu$m (Weingartner & Draine 2001, Draine & Li 2007, Hensley & Draine 2023). While it constitutes at most order-of-magnitude $\mathbf{1\%}$ of the ISM mass, dust is second only to stars in importance for the observable properties of galaxies (Zavala et al. 2021). Large uncertainties in the efficiency of grain growth obfuscate the relative contribution of the two dominant sources of dust in the Universe: direct production from evolved stars versus gas-phase accretion in the ambient ISM (Feldmann 2015, Esmerian & Gnedin 2022, Esmerian & Gnedin 2024). Advances in supercomputers have only recently allowed us to move beyond simple, idealized predictions of dust grain growth efficiencies (Leitch-Devlin & Williams 1985) with atomistic dynamical calculations (Bossion et al. 2024). We show that small carbon dust grains can grow significantly on timescales much shorter than the age of the universe and, in some ISM phases, comparable to the lifetimes of giant molecular clouds. Specifically, we perform molecular dynamics simulations of an amorphous carbon (a-C) grain surface impacted by gas-phase atoms of cosmologically abundant elements with realistic interstellar conditions, finding high ($\gtrsim 0.2$) sticking coefficients for all non-inert elements at all relevant gas and grain temperatures. We present the results of experiments conducted on similar dust candidate materials that support our theoretical calculations. Our results therefore confirm that the process of gas-phase accretion onto grains is likely an efficient mechanism for the growth of interstellar dust mass on astrophysical timescales, and plausibly central to the evolutionary life-cycle of interstellar grains at all cosmic epochs. (abridged)

Figures

Figures reproduced from arXiv: 2607.14237 by Alexey Potapov, Clarke J. Esmerian, Duncan Bossion, Francois Dulieu, Gunnar Nyman, Kirsten K. Knudsen, Saoud Baouche, Susanne Aalto, Tom J. L. C. Bakx, W. M. C. Sameera, Wouter Vlemmings.

Figure 2
Figure 2. Figure 2: The only existing theoretical estimates for similar [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: The coefficients βk for each fitting function are provided in [PITH_FULL_IMAGE:figures/full_fig_p013_3.png] view at source ↗

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

93 extracted references · 10 canonical work pages

  1. [1]

    Investisse- ments d’Avenir

    To ensure that the structure is fully thermalized, the MD is run for 5 ps. In the second step, molecular guns are used to shoot ten colliding atoms, one every 1.25 ps, from different positions and at different colliding angles around the a-C structure in order to obtain sticking coefficients that are averaged over many possible surface sites, and not bias...

    arXiv 2020

  2. [2]

    Draine, B. T. & Li, A. Infrared Emission from Interstel- lar Dust. IV . The Silicate-Graphite-PAH Model in the Post-Spitzer Era.Astrophys. J.657, 810–837 (2007). astro-ph/0608003

    Pith/arXiv arXiv 2007

  3. [3]

    Hensley, B. S. & Draine, B. T. The Astrodust+PAH Model: A Unified Description of the Extinction, Emis- sion, and Polarization from Dust in the Diffuse Interstel- lar Medium.Astrophys. J.948, 55 (2023).2208.123 65

    2023

  4. [4]

    A.et al.The Evolution of the IR Luminos- ity Function and Dust-obscured Star Formation over the Past 13 Billion Years.Astrophys

    Zavala, J. A.et al.The Evolution of the IR Luminos- ity Function and Dust-obscured Star Formation over the Past 13 Billion Years.Astrophys. J.909, 165 (2021). 2101.04734

    Pith/arXiv arXiv 2021

  5. [5]

    The equilibrium view on dust and metals in galaxies: Galactic outflows drive low dust-to-metal ratios in dwarf galaxies.Mon

    Feldmann, R. The equilibrium view on dust and metals in galaxies: Galactic outflows drive low dust-to-metal ratios in dwarf galaxies.Mon. Not. R. Astron. Soc.449, 3274–3292 (2015).1412.2755

    Pith/arXiv arXiv 2015

  6. [6]

    Esmerian, C. J. & Gnedin, N. Y . Modeling Dust Pro- duction, Growth, and Destruction in Reionization-era Galaxies with the CROC Simulations: Methods and Parameter Exploration.Astrophys. J.940, 74 (2022). 2208.02277

    Pith/arXiv arXiv 2022

  7. [7]

    Esmerian, C. J. & Gnedin, N. Y . Modeling Dust Pro- duction, Growth, and Destruction in Reionization-era Galaxies with the CROC Simulations. II. Predicting the Dust Content of High-redshift Galaxies.Astrophys. J. 968, 113 (2024).2308.11723

    Pith/arXiv arXiv 2024

  8. [8]

    Leitch-Devlin, M. A. & Williams, D. A. Sticking co- efficients for atoms and molecules at the surfaces of in- terstellar dust grains.Mon. Not. R. Astron. Soc.213, 295–306 (1985)

    1985

  9. [9]

    Astrophys.692, A249 (2024).2411.06125

    Bossion, D.et al.Accurate sticking coefficient calcu- lation for carbonaceous dust growth through accretion and desorption in astrophysical environments.Astron. Astrophys.692, A249 (2024).2411.06125

    Pith/arXiv arXiv 2024

  10. [10]

    & Rotundi, A

    Mennella, V ., Colangeli, L., Bussoletti, E., Palumbo, P. & Rotundi, A. A new approach to the puzzle of the ul- traviolet interstellar extinction bump.Astrophys. J.507, L177 (1998). URLhttps://doi.org/10.1086/ 311693

    1998

  11. [11]

    & Mutschke, H

    Henning, T., J ¨ager, C. & Mutschke, H. Laboratory Stud- ies of Carbonaceous Dust Analogs. In Witt, A. N., Clay- ton, G. C. & Draine, B. T. (eds.)Astrophysics of Dust, vol. 309 ofAstronomical Society of the Pacific Confer- ence Series, 603 (2004)

    2004

  12. [12]

    Hensley, B. S. & Draine, B. T. Detection of PAH Absorption and Determination of the Mid-infrared Dif- fuse Interstellar Extinction Curve from the Sight Line toward Cyg OB2-12.Astrophys. J.895, 38 (2020). 2002.02457

    Pith/arXiv arXiv 2020

  13. [13]

    J., Jim ´enez-Redondo, M., Pel ´aez, R

    Herrero, V . J., Jim ´enez-Redondo, M., Pel ´aez, R. J., Mat´e, B. & Tanarro, I. Structure and evolution of inter- stellar carbonaceous dust. insights from the laboratory. Frontiers in Astronomy and Space SciencesVolume 9 - 2022(2022). URLhttps://www.frontiersin. org/journals/astronomy-and-space-sci ences/articles/10.3389/fspas.2022.10 83288

    work page doi:10.3389/fspas.2022.10 2022

  14. [14]

    Astrophys.684, A34 (2024).2401.07739

    Ysard, N.et al.THEMIS 2.0: A self-consistent model for dust extinction, emission, and polarisation.Astron. Astrophys.684, A34 (2024).2401.07739. 6

    Pith/arXiv arXiv 2024

  15. [15]

    Hecht, J. H. The Nature of the Dust around R Coronae Borealis Stars: Isolated Amorphous Carbon or Graphite Fractals?Astrophys. J.367, 635 (1991)

    1991

  16. [16]

    L.et al.A comprehensive study of the dust declines in r coronae borealis stars.Mon

    Crawford, C. L.et al.A comprehensive study of the dust declines in r coronae borealis stars.Mon. Not. R. Astron. Soc.537, 2635–2646 (2025). URLhttps:// doi.org/10.1093/mnras/staf215.https: //academic.oup.com/mnras/article-pdf /537/3/2635/61813744/staf215.pdf

    work page doi:10.1093/mnras/staf215.https: 2025

  17. [17]

    Elyajouri, M.et al.PDRs4All. V . Modelling the dust evolution across the illuminated edge of the Orion Bar. Astron. Astrophys.685, A76 (2024).2401.01221

    Pith/arXiv arXiv 2024

  18. [18]

    & Micelotta, E

    Sarangi, A., Matsuura, M. & Micelotta, E. R. Dust in Supernovae and Supernova Remnants I: Formation Sce- narios.Space Sci. Rev.214, 63 (2018)

    2018

  19. [19]

    R., Matsuura, M

    Micelotta, E. R., Matsuura, M. & Sarangi, A. Dust in Supernovae and Supernova Remnants II: Processing and Survival. In Bykov, A.et al.(eds.)Supernovae, vol. 68, 361–418 (2019)

    2019

  20. [20]

    & Olofsson, H

    H ¨ofner, S. & Olofsson, H. Mass loss of stars on the asymptotic giant branch. Mechanisms, models and mea- surements.A&A Rev.26, 1 (2018)

    2018

  21. [21]

    D.et al.Carbon-rich Dust Injected into the Interstellar Medium by Galactic WC Binaries Sur- vives for Hundreds of Years.Astrophys

    Richardson, N. D.et al.Carbon-rich Dust Injected into the Interstellar Medium by Galactic WC Binaries Sur- vives for Hundreds of Years.Astrophys. J.987, 160 (2025).2505.11616

    Pith/arXiv arXiv 2025

  22. [22]

    Draine, B. T. Interstellar Dust Models and Evolutionary Implications. In Henning, T., Gr ¨un, E. & Steinacker, J. (eds.)Cosmic Dust - Near and Far, vol. 414 ofAs- tronomical Society of the Pacific Conference Series, 453 (2009).0903.1658

    Pith/arXiv arXiv 2009

  23. [23]

    Dust Growth and Evolution in Protoplan- etary Disks.Annual Review of Astronomy and Astro- physics62, 157–202 (2024).2312.13287

    Birnstiel, T. Dust Growth and Evolution in Protoplan- etary Disks.Annual Review of Astronomy and Astro- physics62, 157–202 (2024).2312.13287

    Pith/arXiv arXiv 2024

  24. [24]

    T.Physics of the Interstellar and Intergalac- tic Medium(2011)

    Draine, B. T.Physics of the Interstellar and Intergalac- tic Medium(2011)

    2011

  25. [25]

    Draine, B. T. Evolution of interstellar dust. In Blitz, L. (ed.)The Evolution of the Interstellar Medium, vol. 12 ofAstronomical Society of the Pacific Conference Se- ries, 193–205 (1990)

    1990

  26. [26]

    The Evolution of the Elemental Abundances in the Gas and Dust Phases of the Galaxy.Astrophys

    Dwek, E. The Evolution of the Elemental Abundances in the Gas and Dust Phases of the Galaxy.Astrophys. J. 501, 643 (1998).astro-ph/9707024

    Pith/arXiv arXiv 1998

  27. [27]

    Dust Growth Timescale and Mass Func- tion of Molecular Clouds in the Galaxy.Publ

    Hirashita, H. Dust Growth Timescale and Mass Func- tion of Molecular Clouds in the Galaxy.Publ. Astron. Soc. Jpn.52, 585–588 (2000)

    2000

  28. [28]

    Savage, B. D. & Sembach, K. R. Interstellar Abun- dances from Absorption-Line Observations with the Hubble Space Telescope.Annual Review of Astronomy and Astrophysics34, 279–330 (1996)

    1996

  29. [29]

    Jenkins, E. B. A Unified Representation of Gas-Phase Element Depletions in the Interstellar Medium.Astro- phys. J.700, 1299–1348 (2009).0905.3173

    Pith/arXiv arXiv 2009

  30. [30]

    Roman-Duval, J.et al.METAL: The Metal Evolu- tion, Transport, and Abundance in the Large Magellanic Cloud Hubble Program. II. Variations of Interstellar De- pletions and Dust-to-gas Ratio within the LMC.Astro- phys. J.910, 95 (2021).2101.09399

    Pith/arXiv arXiv 2021

  31. [31]

    Watson, D.et al.A dusty, normal galaxy in the epoch of reionization.Nature519, 327–330 (2015).1503.0 0002

    2015

  32. [32]

    Michałowski, M. J. Dust production 680-850 million years after the Big Bang.Astron. Astrophys.577, A80 (2015).1503.08210

    Pith/arXiv arXiv 2015

  33. [33]

    Chenoweth, K., van Duin, A. C. T. & Goddard, W. A. Reaxff reactive force field for molecular dynamics sim- ulations of hydrocarbon oxidation.J. Phys. Chem. A 112, 1040–1053 (2008). URLhttps://doi.org/ 10.1021/jp709896w. PMID: 18197648,https: //doi.org/10.1021/jp709896w

    doi:10.1021/jp709896w 2008

  34. [34]

    van Duin, A. C. T., Dasgupta, S., Lorant, F. & Goddard, W. A. Reaxff: A reactive force field for hydrocarbons. The Journal of Physical Chemistry A105, 9396–9409 (2001). URLhttps://doi.org/10.1021/jp00 4368u

    work page doi:10.1021/jp00 2001

  35. [35]

    Planck Collaborationet al.Planck 2018 results. VI. Cosmological parameters.Astron. Astrophys.641, A6 (2020).1807.06209

    Pith/arXiv arXiv 2018

  36. [36]

    S., Rood, R

    Balser, D. S., Rood, R. T., Bania, T. M. & Anderson, L. D. H II Region Metallicity Distribution in the Milky Way Disk.Astrophys. J.738, 27 (2011).1106.1660

    Pith/arXiv arXiv 2011

  37. [37]

    Astrophys.707, A54 (2026).2508.19339

    Hansson, K.et al.The binding energies of atoms on amorphous silicate dust: A computational study.Astron. Astrophys.707, A54 (2026).2508.19339

    arXiv 2026

  38. [38]

    J.760, 149 (2012).1210.3631

    Paladini, R.et al.Spitzer and Herschel Multiwavelength Characterization of the Dust Content of Evolved H II Regions.Astrophys. J.760, 149 (2012).1210.3631

    Pith/arXiv arXiv 2012

  39. [39]

    Binder, B. A. & Povich, M. S. A Multiwavelength Look at Galactic Massive Star-forming Regions.Astrophys. J. 864, 136 (2018).1808.00454

    Pith/arXiv arXiv 2018

  40. [40]

    & Dulieu, F

    Grieco, F., Theul ´e, P., De Looze, I. & Dulieu, F. En- hanced star formation through the high-temperature for- mation of H 2 on carbonaceous dust grains.Nature As- tronomy7, 541–545 (2023). 7

    2023

  41. [41]

    Pirronello, V ., Liu, C., Roser, J. E. & Vidali, G. Mea- surements of molecular hydrogen formation on carbona- ceous grains.Astron. Astrophys.344, 681–686 (1999)

    1999

  42. [42]

    V ., Field, D

    Hornekær, L., Baurichter, A., Petrunin, V . V ., Field, D. & Luntz, A. C. Importance of Surface Morphology in Interstellar H 2 Formation.Science302, 1943–1946 (2003)

    1943

  43. [43]

    H., Dulieu, F., Momeni, A

    Amiaud, L., Fillion, J. H., Dulieu, F., Momeni, A. & Lemaire, J. L. Probing amorphous solid water by temperature programmed desorption experiments with molecular hydrogen: from compact to porous ices. In Lemaire, J. L. & Combes, F. (eds.)Molecules in Space and Laboratory, 71 (2007)

    2007

  44. [44]

    Kwok, S.Physics and Chemistry of the Interstellar Medium(2007)

    2007

  45. [45]

    & Troland, T

    Heiles, C. & Troland, T. H. The Millennium Arecibo 21 Centimeter Absorption-Line Survey. II. Properties of the Warm and Cold Neutral Media.Astrophys. J.586, 1067–1093 (2003).astro-ph/0207105

    Pith/arXiv arXiv 2003

  46. [46]

    Chevance, M.et al.The lifecycle of molecular clouds in nearby star-forming disc galaxies.Mon. Not. R. Astron. Soc.493, 2872–2909 (2020).1911.03479

    Pith/arXiv arXiv 2020

  47. [47]

    Peters, T.et al.The turbulent life of dust grains in the supernova-driven, multiphase interstellar medium.Mon. Not. R. Astron. Soc.467, 4322–4342 (2017).1610.0 6579

    2017

  48. [48]

    & Maiolino, R

    Schneider, R. & Maiolino, R. The formation and cosmic evolution of dust in the early Universe: I. Dust sources. A&A Rev.32, 2 (2024).2310.00053

    Pith/arXiv arXiv 2024

  49. [49]

    L., Whittet, D

    Gibb, E. L., Whittet, D. C. B., Boogert, A. C. A. & Tie- lens, A. G. G. M. Interstellar Ice: The Infrared Space Observatory Legacy.Astrophys. J. Suppl.151, 35–73 (2004)

    2004

  50. [50]

    Boogert, A. C. A.et al.Ice and Dust in the Quiescent Medium of Isolated Dense Cores.Astrophys. J.729, 92 (2011).1101.2486

    Pith/arXiv arXiv 2011

  51. [51]

    Boogert, A. C. A., Gerakines, P. A. & Whittet, D. C. B. Observations of the icy universe.Annual Review of Astronomy and Astrophysics53, 541–581 (2015). 1501.05317

    Pith/arXiv arXiv 2015

  52. [52]

    Decleir, M.et al.A First Taste of MEAD (Measuring Extinction and Abundances of Dust). I. Diffuse Milky Way Interstellar Dust Extinction Features in JWST In- frared Spectra.AJ169, 99 (2025).2412.14378

    Pith/arXiv arXiv 2025

  53. [53]

    Astrophys.512, A30 (2010).0903.3 120

    Dulieu, F.et al.Experimental evidence for water forma- tion on interstellar dust grains by hydrogen and oxygen atoms.Astron. Astrophys.512, A30 (2010).0903.3 120

    2010

  54. [54]

    & Ceccarelli, C

    Ferrara, A., Viti, S. & Ceccarelli, C. The problematic growth of dust in high-redshift galaxies.Mon. Not. R. Astron. Soc.463, L112–L116 (2016).1606.07214

    Pith/arXiv arXiv 2016

  55. [55]

    & Henning, T

    Potapov, A., J ¨ager, C. & Henning, T. Ice Coverage of Dust Grains in Cold Astrophysical Environments. Phys. Rev. Lett.124, 221103 (2020).2005.00757

    Pith/arXiv arXiv 2020

  56. [56]

    Tielens, A. G. G. M., McKee, C. F., Seab, C. G. & Hol- lenbach, D. J. The Physics of Grain-Grain Collisions and Gas-Grain Sputtering in Interstellar Shocks.Astro- phys. J.431, 321 (1994)

    1994

  57. [57]

    P., Tielens, A

    Jones, A. P., Tielens, A. G. G. M. & Hollenbach, D. J. Grain Shattering in Shocks: The Interstellar Grain Size Distribution.Astrophys. J.469, 740 (1996)

    1996

  58. [58]

    J.et al.Silicate cosmic dust grain col- lisions in the interstellar medium: A molecular dy- namics study.arXiv e-printsarXiv:2605.22225 (2026)

    Esmerian, C. J.et al.Silicate cosmic dust grain col- lisions in the interstellar medium: A molecular dy- namics study.arXiv e-printsarXiv:2605.22225 (2026). 2605.22225

    Pith/arXiv arXiv 2026

  59. [59]

    Whittet, D. C. B.et al.Infrared Spectroscopy of Dust in the Diffuse Interstellar Medium toward Cygnus OB2 No. 12.Astrophys. J.490, 729–734 (1997)

    1997

  60. [60]

    Weingartner, J. C. & Draine, B. T. Photoelectric Emis- sion from Interstellar Dust: Grain Charging and Gas Heating.Astrophys. J. Suppl.134, 263–281 (2001). astro-ph/9907251

    Pith/arXiv arXiv 2001

  61. [61]

    C., Draine, B

    Weingartner, J. C., Draine, B. T. & Barr, D. K. Pho- toelectric Emission from Dust Grains Exposed to Ex- treme Ultraviolet and X-Ray Radiation.Astrophys. J. 645, 1188–1197 (2006).astro-ph/0601296

    Pith/arXiv arXiv 2006

  62. [62]

    Weingartner, J. C. & Draine, B. T. Interstellar Depletion onto Very Small Dust Grains.Astrophys. J.517, 292– 298 (1999).astro-ph/9812469

    Pith/arXiv arXiv 1999

  63. [63]

    C.et al.Dust charge distribution in the interstellar medium.Mon

    Ib ´a˜nez-Mej´ıa, J. C.et al.Dust charge distribution in the interstellar medium.Mon. Not. R. Astron. Soc.485, 1220–1247 (2019).1812.08281

    Pith/arXiv arXiv 2019

  64. [64]

    A., Kravtsov, A

    Semenov, V . A., Kravtsov, A. V . & Gnedin, N. Y . The Physical Origin of Long Gas Depletion Times in Galax- ies.Astrophys. J.845, 133 (2017).1704.04239

    Pith/arXiv arXiv 2017

  65. [65]

    & Svoboda, B

    Battersby, C., Bally, J. & Svoboda, B. The Lifetimes of Phases in High-mass Star-forming Regions.Astrophys. J.835, 263 (2017).1702.02199

    Pith/arXiv arXiv 2017

  66. [66]

    Bonne, L.et al.The SOFIA FEEDBACK [CII] Legacy Survey: Rapid molecular cloud dispersal in RCW 79. Astron. Astrophys.679, L5 (2023).2310.01657

    Pith/arXiv arXiv 2023

  67. [67]

    Faerber, T.et al.Expansion Signatures in 35 H II Re- gions Traced by SOFIA [C II] Emission.Astrophys. J. 990, 30 (2025).2506.16700. 8

    Pith/arXiv arXiv 2025

  68. [68]

    Cautun, M.et al.The milky way total mass profile as inferred from Gaia DR2.Mon. Not. R. Astron. Soc.494, 4291–4313 (2020).1911.04557

    Pith/arXiv arXiv 2020

  69. [69]

    Astrodust

    Draine, B. T. & Hensley, B. S. The Dielectric Function of “Astrodust” and Predictions for Polarization in the 3.4 and 10µm Features.Astrophys. J.909, 94 (2021).20 09.11314

    2021

  70. [70]

    D.et al.Milky Way Mid-Infrared Spitzer Spectroscopic Extinction Curves: Continuum and Sili- cate Features.Astrophys

    Gordon, K. D.et al.Milky Way Mid-Infrared Spitzer Spectroscopic Extinction Curves: Continuum and Sili- cate Features.Astrophys. J.916, 33 (2021).2105.0 5087

    2021

  71. [71]

    Stephan, T.et al.The Presolar Grain Database. I. Silicon Carbide.Astrophys. J. Suppl.270, 27 (2024)

    2024

  72. [72]

    Astrophys.462, 667–676 (2007)

    Min, M.et al.The shape and composition of interstellar silicate grains.Astron. Astrophys.462, 667–676 (2007). astro-ph/0611329

    Pith/arXiv arXiv 2007

  73. [73]

    Roman-Duval, J.et al.METAL: The Metal Evolu- tion, Transport, and Abundance in the Large Magel- lanic Cloud Hubble Program. III. Interstellar Deple- tions, Dust-to-Metal, and Dust-to-Gas Ratios versus Metallicity.Astrophys. J.928, 90 (2022).2202.047 65

    2022

  74. [74]

    Whittet, D. C. B. Oxygen Depletion in the Interstellar Medium: Implications for Grain Models and the Distri- bution of Elemental Oxygen.Astrophys. J.710, 1009– 1016 (2010).0912.3298

    Pith/arXiv arXiv 2010

  75. [75]

    A., Whittet, D

    Poteet, C. A., Whittet, D. C. B. & Draine, B. T. The Composition of Interstellar Grains towardζOphiuchi: Constraining the Elemental Budget near the Diffuse- dense Cloud Transition.Astrophys. J.801, 110 (2015). 1501.02810

    Pith/arXiv arXiv 2015

  76. [76]

    J., Hua, L

    Yang, X. J., Hua, L. & Li, A. Where Have All the Sul- fur Atoms Gone? Polycyclic Aromatic Hydrocarbon as a Possible Sink for the Missing Sulfur in the Interstel- lar Medium. I. The C—S Band Strengths.Astrophys. J. 974, 30 (2024).2408.04857

    Pith/arXiv arXiv 2024

  77. [77]

    L.et al.Computational surface chemistry of tetrahedral amorphous carbon by combining machine learning and density functional theory.Chem

    Deringer, V . L.et al.Computational surface chemistry of tetrahedral amorphous carbon by combining machine learning and density functional theory.Chem. Mater.30, 7438–7445 (2018). URLhttps://doi.org/10.1 021/acs.chemmater.8b02410.https://do i.org/10.1021/acs.chemmater.8b02410

    work page doi:10.1021/acs.chemmater.8b02410 2018

  78. [78]

    J.et al.The amsterdam modeling suite.J

    Baerends, E. J.et al.The amsterdam modeling suite.J. Chem. Phys.162, 162501 (2025). URLhttps://do i.org/10.1063/5.0258496.https://pubs .aip.org/aip/jcp/article-pdf/doi/10. 1063/5.0258496/20497318/162501_1_5.0 258496.pdf

    work page doi:10.1063/5.0258496.https://pubs 2025

  79. [79]

    van Duin, A.et al.Reaxff 2025.1, scm, theoreti- cal chemistry, vrije universiteit, amsterdam, the nether- lands, http://www.scm.com

    2025

  80. [80]

    & van Duin, A

    Hong, S. & van Duin, A. C. T. Atomistic-scale anal- ysis of carbon coating and its effect on the oxidation of aluminum nanoparticles by reaxff-molecular dynam- ics simulations.The Journal of Physical Chemistry C 120, 9464–9474 (2016). URLhttps://doi.or g/10.1021/acs.jpcc.6b00786.https: //doi.org/10.1021/acs.jpcc.6b00786

    work page doi:10.1021/acs.jpcc.6b00786.https: 2016

Showing first 80 references.