Модель плотноупакованных случайных сфер¶

Модель плотноупакованных случайных сфер (также известная как модель случайной плотной упаковки, Random Close Packing, RCP) — это концептуальная и вычислительная модель, описывающая конфигурацию твёрдых непересекающихся сфер одинакового радиуса, которые уложены в пространстве с максимально возможной плотностью при условии отсутствия дальнего порядка (кристаллической решётки). В отличие от правильных кристаллических упаковок (например, гранецентрированной кубической или гексагональной), в RCP сферы расположены хаотично, но при этом не могут быть сдвинуты без увеличения общего объёма. Данная модель является ключевой для понимания структуры аморфных твёрдых тел, жидкостей с высокой вязкостью, гранулированных материалов и коллоидных систем.

¶История

Исследования плотной упаковки сфер восходят к XVII веку, когда Иоганн Кеплер в трактате «О шестиугольных снежинках» (1611) выдвинул гипотезу о том, что максимальная плотность упаковки равных сфер достигается в гранецентрированной кубической решётке (плотность π/√18 ≈ 0,74048). Однако интерес к случайным упаковкам возник лишь в середине XX века в связи с развитием физики конденсированного состояния и материаловедения.

В 1950-х годах британский физик Джон Бернал (John Bernal) провёл серию экспериментов с пластилиновыми шариками, сжимая их в резиновых мешках. Он обнаружил, что при случайном заполнении объёма плотность упаковки составляет около 0,64. Бернал ввёл понятие «случайной плотной упаковки» и описал её как модель для структуры простых жидкостей. В 1970-х годах Джордж Скотт (George Scott) и другие исследователи уточнили это значение, получив плотность 0,6366 ± 0,0005 для бесконечно большого объёма.

С развитием компьютерного моделирования в 1980–1990-х годах стало возможным численно воспроизводить случайные упаковки. Методы, такие как алгоритм Лубачевского–Стиллинджера (Lubachevsky–Stillinger), позволяют генерировать конфигурации с заданной плотностью, имитируя рост сфер из точечных зародышей. Современные исследования показывают, что RCP не является единственным состоянием, а представляет собой часть более широкого континуума аморфных упаковок.

¶Определение и ключевые характеристики

¶Плотность упаковки

Основной параметр модели — плотность упаковки φ, определяемая как отношение суммарного объёма всех сфер к объёму вмещающего контейнера. Для RCP φ ≈ 0,64. Это значение значительно ниже плотности кристаллических упаковок (0,74), но выше плотности случайной рыхлой упаковки (Random Loose Packing, RLP), где φ ≈ 0,56. Разница между RCP и RLP обусловлена степенью механической стабильности: в RCP частицы не могут быть сдвинуты без увеличения объёма, тогда как в RLP они могут быть уплотнены при встряхивании.

¶Координационное число

Среднее число контактов каждой сферы с соседями (координационное число) для RCP составляет около 6. Это меньше, чем в кристаллических решётках (12), но больше, чем в рыхлых упаковках. Распределение контактов в RCP широкое: от 4 до 10, что отражает структурную неоднородность.

¶Радиальная функция распределения

Характерной особенностью RCP является радиальная функция распределения g(r), которая показывает вероятность найти соседнюю сферу на расстоянии r от центра данной. Для RCP g(r) имеет ярко выраженный первый пик на расстоянии 2R (где R — радиус сферы), затем несколько меньших пиков, соответствующих вторым и третьим соседям, и при r → ∞ стремится к 1. Отсутствие дальних пиков свидетельствует об отсутствии кристаллического порядка.

¶Методы моделирования

¶Экспериментальные методы

Классические эксперименты Бернала включали ручное сжатие шариков в эластичных контейнерах. Современные подходы используют:

• Сжатие коллоидных суспензий: полимерные или кремниевые микросферы в жидкости осаждаются под действием центрифугирования, формируя случайную упаковку.
• Гранулированные материалы: песок, стальные шарики или пластиковые гранулы уплотняются вибрацией или механическим давлением.

¶Компьютерные алгоритмы

• Алгоритм Лубачевского–Стиллинджера: сферы начинают как точки, которые равномерно растут в объёме с заданной скоростью, при этом их центры отталкиваются при соприкосновении. Процесс останавливается при достижении целевой плотности.
• Метод Монте-Карло: случайные сдвиги сфер принимаются или отвергаются по критерию Метрополиса, что позволяет достичь равновесной конфигурации.
• Метод молекулярной динамики: сферы взаимодействуют через потенциал отталкивания (например, потенциал Леннарда-Джонса), и система охлаждается до застывания в случайной конфигурации.

¶Применение

¶Физика конденсированного состояния

RCP служит моделью для структуры аморфных твёрдых тел, таких как стекло, аморфный кремний и полимеры. В отличие от кристаллов, аморфные материалы не имеют дальнего порядка, и их свойства (вязкость, теплопроводность, прочность) определяются локальной структурой, близкой к RCP.

¶Гранулированные материалы

В промышленности (строительство, фармацевтика, пищевая промышленность) RCP описывает укладку сыпучих веществ: песка, цемента, порошков. Знание плотности упаковки позволяет оптимизировать дозирование, транспортировку и хранение.

¶Коллоидные системы

Коллоидные суспензии, содержащие микросферы, при высокой концентрации образуют случайные плотные упаковки. Это явление используется в производстве фотонных кристаллов, сенсоров и материалов с управляемой пористостью.

¶Геофизика

RCP применяется для моделирования структуры осадочных пород, таких как песчаники, где зёрна минералов упакованы случайным образом. Плотность упаковки влияет на пористость, проницаемость и механическую прочность горных пород.

¶Критика и ограничения

¶Неоднозначность определения

Термин «плотноупакованные случайные сферы» критикуется за отсутствие строгого математического определения. В отличие от кристаллических решёток, где плотность фиксирована, RCP может варьироваться в зависимости от метода получения и граничных условий. Некоторые исследователи (например, Стивен Торквато) утверждают, что RCP — это не уникальное состояние, а часть континуума «джемминга» (jamming), где система может быть как плотной, так и рыхлой в зависимости от давления.

¶Влияние размера контейнера

В экспериментах с конечным числом сфер (например, 10⁴–10⁶ частиц) плотность RCP может отклоняться от теоретического значения 0,64 из-за краевых эффектов. Для бесконечно большого объёма плотность стремится к 0,6366, но в реальных системах это значение достигается лишь приближённо.

¶Альтернативные модели

В 2000-х годах была предложена модель «случайной упаковки с максимальной плотностью» (Maximally Random Jammed State, MRJ), которая, по мнению её авторов, является более точной альтернативой RCP. MRJ предполагает, что система является не только плотной, но и максимально неупорядоченной, что достигается при определённых условиях (например, при медленном охлаждении). Плотность MRJ также близка к 0,64, но её определение строже.

¶Интересные факты

• Парадокс Кеплера: долгое время считалось, что случайная упаковка не может быть плотнее 0,64, но в 2010-х годах были найдены конфигурации с плотностью до 0,68, однако они не являются полностью случайными и содержат элементы локального порядка.
• Связь с жидкостями: RCP часто рассматривается как модель для «замороженной» жидкости, так как радиальная функция распределения RCP напоминает таковую для простых жидкостей при высокой плотности.
• Применение в косметологии: в производстве пудр и тональных кремов используется принцип RCP для достижения оптимальной плотности пигментов без образования комков.

¶Источники

• Bernal, J. D. (1960). «Geometry of the Structure of Monatomic Liquids». Nature.
• Scott, G. D. (1962). «Packing of Spheres». Nature.
• Torquato, S. (2002). «Random Heterogeneous Materials: Microstructure and Macroscopic Properties». Springer.
• Lubachevsky, B. D., Stillinger, F. H. (1990). «Geometric Properties of Random Disk Packings». Journal of Statistical Physics.
• Donev, A. et al. (2004). «Improving the Density of Jammed Disordered Packings». Physical Review Letters.