🔬 mesoscale physics

Silicate cosmic dust grain collisions in the interstellar medium: A molecular dynamics study

Using molecular dynamics simulations of colliding silicate grains, this study reveals that shattering velocity thresholds are approximately 6 km/s—significantly higher than previously assumed due to a correction in earlier theoretical models—and demonstrates that existing models fail to accurately predict the resulting mass fractions and size distributions of shattered products.

Original authors: C. J. Esmerian, S. R. Hashemi, W. M. C. Sameera, W. Vlemmings, S. Andersson, T. J. L. C. Bakx, K. K. Knudsen, S. Aalto, G. Nyman

Published 2026-05-22

📖 4 min read☕ Coffee break read

CC BY 4.0

Original authors: C. J. Esmerian, S. R. Hashemi, W. M. C. Sameera, W. Vlemmings, S. Andersson, T. J. L. C. Bakx, K. K. Knudsen, S. Aalto, G. Nyman

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine the space between stars (the Interstellar Medium) not as an empty void, but as a bustling, invisible highway filled with tiny specks of dust. These aren't just random dirt; they are cosmic dust grains, mostly made of silicates (think of them as microscopic, rocky sand). They are crucial because they act like tiny factories where new molecules form and like shields that protect those molecules from being blasted apart by starlight.

For a long time, astronomers had a "rulebook" for what happens when two of these dust grains crash into each other. They believed that if two grains hit each other at a speed of about 2.7 kilometers per second (roughly 6,000 miles per hour), they would shatter into tiny pieces, like dropping a ceramic plate on the floor. If they hit even faster, they would vaporize, turning instantly into gas.

The New Experiment: A High-Speed Crash Test
In this paper, a team of scientists decided to test that old rulebook using a super-powerful computer simulation. Instead of dropping real dust grains (which are too small and fast to catch in a lab), they built digital models of these grains atom-by-atom.

Think of it like a video game crash test, but instead of cars, they are smashing together two perfect spheres of "digital sand." They simulated collisions at speeds ranging from a gentle 0.1 km/s up to a blistering 20 km/s. They tested two types of "sand": pure silica (like glass) and a more complex mix called "astrodust" (which contains iron and magnesium, like the rocks in our solar system).

The Big Surprise: The Dust is Tougher Than We Thought
The results were a shock to the system. The old rulebook said the dust would break at 2.7 km/s. The new computer experiments showed that the dust grains are actually much tougher. They didn't start shattering until they hit speeds of about 6 km/s.

Why the Old Rulebook Was Wrong
The authors found that the old rulebook wasn't just slightly off; it had a math error in its foundation. It was like a recipe that said "add 2 cups of flour" when it actually meant "add 4 cups." When they fixed the math in the old theory, the predicted breaking speed jumped up to about 7.9 km/s. This new, corrected number is much closer to what their computer simulations actually showed (around 6 km/s).

So, the main takeaway is: Cosmic dust is more durable than we previously thought. It can survive much faster crashes than we assumed.

What Happens When They Do Break?
When the grains finally did break in the simulations, the results didn't look like the old theories predicted either.

The Old Theory: Predicted that broken dust would follow a neat, predictable pattern (like a smooth slide where you get a specific number of big pieces and a specific number of tiny pieces).
The Reality: The broken pieces were messy and chaotic. The size of the fragments depended heavily on exactly how fast they were going and how big the original grains were. There was no single "perfect pattern."

Also, the old theory guessed that a certain amount of dust would turn into gas (vaporize) at high speeds. The simulations showed that the old theory was way too optimistic about shattering and way too pessimistic about vaporization. In reality, the grains held together longer, and when they did break, they didn't turn into gas as easily as the old models suggested.

Why Does This Matter?
This changes how we understand the "lifecycle" of dust in the universe.

Resilience: Because the dust is tougher, it survives longer in the harsh environment of space. It doesn't get destroyed as quickly by collisions.
Growth: Since the grains don't shatter as easily, they might be more likely to stick together (coagulate) to form bigger grains, rather than being smashed into dust.
The Math: Astronomers who build models of how galaxies evolve will need to update their calculations. They can't use the old "2.7 km/s" breaking point anymore; they need to use the new, higher speed limits to get accurate pictures of how dust behaves in the universe.

In short, this paper is a "crash test" for the universe's smallest building blocks. It tells us that cosmic dust is a lot more resilient than we gave it credit for, and it fixes a decades-old math mistake that has been used in astronomy textbooks.

Technical Summary: Silicate Cosmic Dust Grain Collisions in the Interstellar Medium

Problem Statement
Cosmic dust grains are fundamental components of the interstellar medium (ISM), influencing gas thermodynamics, chemistry, and star formation. A critical process governing the evolution of the dust grain size distribution is grain-grain collisions, which can result in coagulation, shattering, or vaporization. Existing astrophysical models rely heavily on analytical frameworks established by Tielens et al. (1994) and Jones et al. (1996). These models predict shattering thresholds for silicate grains at approximately 2.7 km/s and assume continuum fluid-dynamical behavior with plane-parallel collision geometries. However, these assumptions are increasingly questionable for the smallest interstellar grains (radii $\sim$ 0.5–5 nm), where the continuum approximation breaks down. Furthermore, laboratory experiments mimicking these high-velocity, small-grain dynamics are extremely challenging. Consequently, there is a lack of atomistic data to validate or refine the canonical thresholds and outcome distributions used in ISM dust evolution models.

Methodology
This study employs classical Molecular Dynamics (MD) simulations to model head-on collisions between silicate dust grains for the first time at astrophysically relevant velocities (0.1–20 km/s). The authors simulate grains with radii $a \in [5, 50]$ Å, representing the lower end of the expected ISM grain size distribution.

Two distinct material compositions are investigated:

Amorphous Silica ( $SiO_2$ ): Modeled using the reactive force field (ReaxFF) potential.
"Astrodust" Silicates: Modeled with the composition $Mg_{1.3}Fe_{0.3}Si_1O_{3.6}$ , inferred from the Draine & Hensley (2021) model. These are simulated using the semi-empirical quantum mechanical method GFN1-xTB due to the complexity of the multi-component system.

The simulations utilize the AMS driver program to initiate collisions between spherical, non-porous, amorphous grains. The study focuses on head-on collisions (impact parameter $b=0$ ) to isolate the dependence of collision outcomes on relative velocity. Post-simulation analysis categorizes products into "largest grain," "shattered" (grain-sized fragments), and "vaporized" (gas-phase atoms/molecules) based on a geometric radius threshold of 4.5 Å.

Key Results

Shattering Thresholds: The simulations reveal that the velocity threshold for grain shattering is approximately 6 km/s for both amorphous $SiO_2$ and astrodust materials. This is roughly a factor of two higher than the canonical value of 2.7 km/s cited in Jones et al. (1996).
Resolution of Discrepancy: The authors demonstrate that the discrepancy arises primarily from an algebraic error in the derivation of the velocity threshold expression in Tielens et al. (1994). When corrected (Eq. 5 in the paper), the theoretical prediction for silicates rises to $\sim$ 7.9 km/s, which is in moderate agreement with the simulation results (over-predicting by $\sim$ 30%).
Mass Fractions: The fraction of mass that is shattered versus vaporized does not follow the predictions of Tielens et al. (1994) or Hirashita & Kobayashi (2013). The continuum-based expressions over-predict the shattered fraction and under-predict the vaporized fraction for the small grains simulated.
Size Distributions: The size distributions of shattered products do not conform to the single power-law distributions ( $dn/da \propto a^{-3.3}$ ) predicted by Jones et al. (1996). Instead, the distributions are complex, velocity-dependent, and often multimodal, likely due to the breakdown of continuum assumptions and changes in grain density/shape during collision.
Mass Ratio Effects: Collisions between grains with significantly different masses (e.g., a 100:1 ratio) do not disrupt the larger grain at ISM-relevant velocities, suggesting shattering and vaporization are primarily relevant for collisions between grains of comparable mass (within a factor of 10).

Significance and Claims
The authors claim to provide the first atomistic validation of grain collision outcomes for silicate nanoparticles in the ISM. The primary contribution is the provision of updated shattering velocity thresholds ( $\sim$ 6 km/s) for standard candidate grain materials.

The paper argues that these updated thresholds suggest astrophysical dust grains, particularly silicates, are more robust to shattering in the ISM than previously assumed. While the authors acknowledge that the precise impact on global dust population models is complex due to the interplay of multiple physical processes (turbulence, gas drag, etc.), they note that the higher thresholds could increase the predicted maximum grain sizes in certain ISM phases (e.g., the Warm Ionized Medium) by up to an order of magnitude, depending on local gas density assumptions.

The study concludes that while the continuum dynamical framework of Tielens et al. (1994) remains conceptually useful for defining thresholds, its specific expressions for mass fractions and size distributions fail to describe the atomistic reality of small grain collisions. The authors recommend adopting the corrected threshold velocity expression (Eq. 5) for future astrophysical models but caution that prescriptions for shattered mass fractions and size distributions for small grains remain an open area requiring further simulation and theoretical development.

More like this