Представлено результати динамічного моделювання процесу формування пилового хвоста комети C/1995 O1 (Хейла—Боппа). Для побудови моделі пилового хвоста були простежені траєкторії 2·10⁶ пилових частинок. Їхні розміри, напрями, моменти вильоту та початкові швидкості були задані за допомогою алгоритму Монте-Карло. Отриманий тривимірний модельний хвіст був спроектований на небесну сферу для порівняння з даними спостережень. Згідно з нашими модельними дослідженнями спостережувані хвости могли бути сформовані частинками розмірами від 0.3 до 8.0 мкм, швидкості вильоту яких становили від 0.155 до 0.670 км/с з експоненціальним розподілом по розмірах з показником степеня –3.6…–3.7. Показано, що врахування процесу фрагментації частинок призводить до помітного покращення результатів моделювання.