Trajectory-based dust evolution in disks: first results from the RAPID simulation code

Abstract The rapid depletion of dust particles in protoplanetary disks limits the time available for planetesimal formation, as solids are typically accreted onto the central star before dust particles can undergo substantial growth. Dust traps formed at sharp viscosity transitions—such as at the edges of the accretionally inactive dead zones—can halt radial drift and enhance dust coagulation. In this study, dust dynamics is investigated using , a one-dimensional Lagrangian–Eulerian simulation code that tracks representative particle trajectories over time. In order to explore the effect of physical parameters on dust evolution, a grid of 243 models was run. The simulation grid covers a range of parameters such as viscosity, width of the transition region at the edges of the dead zone, disk surface density exponent, and the collisional fragmentation velocity of the dust particles. The computational domain extends from 1 to 50 AU and covers $$5\times 10^5$$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:mrow> <mml:mn>5</mml:mn> <mml:mo>×</mml:mo> <mml:msup> <mml:mn>10</mml:mn> <mml:mn>5</mml:mn> </mml:msup> </mml:mrow> </mml:math> years of disk evolution, assuming a disk mass of $$\sim 0.005\,M_\odot $$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:mrow> <mml:mo>∼</mml:mo> <mml:mn>0.005</mml:mn> <mml:mspace/> <mml:msub> <mml:mi>M</mml:mi> <mml:mo>⊙</mml:mo> </mml:msub> </mml:mrow> </mml:math> . The results show that pressure maxima can trap up to $$3-10\,M_\oplus $$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:mrow> <mml:mn>3</mml:mn> <mml:mo>-</mml:mo> <mml:mn>10</mml:mn> <mml:mspace/> <mml:msub> <mml:mi>M</mml:mi> <mml:mo>⊕</mml:mo> </mml:msub> </mml:mrow> </mml:math> of dust, depending on the local disk conditions. However, increasing the fragmentation velocity, decreasing the viscosity, or widening the dead zone transition width tend to reduce the effectiveness of dust trapping. The simulation results with reveal that dust evolution is highly sensitive to the physical conditions of the disk, which governs the early stages of planetesimal growth.

Authors

Institutions

Publication Details

Journal
Celestial Mechanics and Dynamical Astronomy
Published
2026-02-01
DOI
https://doi.org/10.1007/s10569-026-10278-2
Primary Topic
Astrophysics and Star Formation Studies
Type
article
Field-Weighted Citation Impact
0.00

Funders

Controls
|||
ALL TIME
JAN
FEB
MAR
APR
MAY
JUN
JUL
AUG
SEP
OCT
article

Trajectory-based dust evolution in disks: first results from the RAPID simulation code

Dóra Tarczay-Nehéz
Celestial Mechanics and Dynamical Astronomy
Astrophysics and Star Formation Studies
article

Trajectory-based dust evolution in disks: first results from the RAPID simulation code

Dóra Tarczay-Nehéz
article en

Abstract

Abstract The rapid depletion of dust particles in protoplanetary disks limits the time available for planetesimal formation, as solids are typically accreted onto the central star before dust particles can undergo substantial growth. Dust traps formed at sharp viscosity transitions—such as at the edges of the accretionally inactive dead zones—can halt radial drift and enhance dust coagulation. In this study, dust dynamics is investigated using , a one-dimensional Lagrangian–Eulerian simulation code that tracks representative particle trajectories over time. In order to explore the effect of physical parameters on dust evolution, a grid of 243 models was run. The simulation grid covers a range of parameters such as viscosity, width of the transition region at the edges of the dead zone, disk surface density exponent, and the collisional fragmentation velocity of the dust particles. The computational domain extends from 1 to 50 AU and covers $$5\times 10^5$$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:mrow> <mml:mn>5</mml:mn> <mml:mo>×</mml:mo> <mml:msup> <mml:mn>10</mml:mn> <mml:mn>5</mml:mn> </mml:msup> </mml:mrow> </mml:math> years of disk evolution, assuming a disk mass of $$\sim 0.005\,M_\odot $$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:mrow> <mml:mo>∼</mml:mo> <mml:mn>0.005</mml:mn> <mml:mspace/> <mml:msub> <mml:mi>M</mml:mi> <mml:mo>⊙</mml:mo> </mml:msub> </mml:mrow> </mml:math> . The results show that pressure maxima can trap up to $$3-10\,M_\oplus $$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:mrow> <mml:mn>3</mml:mn> <mml:mo>-</mml:mo> <mml:mn>10</mml:mn> <mml:mspace/> <mml:msub> <mml:mi>M</mml:mi> <mml:mo>⊕</mml:mo> </mml:msub> </mml:mrow> </mml:math> of dust, depending on the local disk conditions. However, increasing the fragmentation velocity, decreasing the viscosity, or widening the dead zone transition width tend to reduce the effectiveness of dust trapping. The simulation results with reveal that dust evolution is highly sensitive to the physical conditions of the disk, which governs the early stages of planetesimal growth.

Celestial Mechanics and Dynamical AstronomyVol. 138(1)
Konkoly Observatory (HU)
Magyar Tudományos Akadémia, Nemzeti Kutatási Fejlesztési és Innovációs Hivatal
Openalex Percentile: Top 67%
Astrophysics and Star Formation Studies
AI Navigator

Ask Laika to Summarize, Analyze, and Connect papers live on the map.

Summarize Papers & Methodologies

Extract key findings, datasets, and comparative methods across publications.

Benchmark Rankings & Visual Analytics

Rank top research institutions, authors, funders, topics, and journals by Field-Weighted Citation Impact (FWCI) and paper volume with instant charts.

Connect Distant Disciplines

Bridge topological clusters on the map to find hidden collaborative intersections.