Предположим, у меня есть симуляция группы молекул азота (N_2). Просто взглянув на положение группы молекул в пространстве, как я могу определить, что изображено на картинке: твердое тело, жидкость или газ?
На каком расстоянии они должны быть друг от друга, чтобы я сделал вывод, что это газ, а не жидкость? Насколько близко они должны быть друг к другу, прежде чем я сделаю вывод, что это твердое тело?
Мотивация: я хотел бы проводить симуляции различной сложности для получения таких величин, как точка замерзания и точка кипения при различных давлениях. Для субрегионов, например вокселей, я хочу иметь возможность задать вопрос: Это твердое тело, жидкость или газ?
Качественные ли это различия или количественные?
В газе расстояние между молекулами явно больше. Но между аморфным твердым телом и жидкостью порядок и расстояния эквивалентны (хуже того: помните, что лед менее плотен, чем жидкая вода), так что вы не сможете сказать. Для кристаллического твердого тела порядок будет легче обнаружить, но полярные жидкости также имеют некоторый ближний порядок.
Конечно, очень близко к переходу (или, что еще хуже, к тройной точке), вы не сможете определить фазу только по позициям. Заметим также, что фаза не ограничивается твердой/жидкой/фазой: помимо экстремальных состояний существуют также различные перестройки кристаллов в зависимости от температуры и давления.
Идентифицировать твердую фазу проще всего. вы увидите кристалл, организованную структуру с шестью соседями для каждой частицы. (теорема Мермина – Вагнера здесь не применима, потому что взаимодействие имеет достаточно большой диапазон. Моделирование и снимки твердой фазы см. В этой статье ). вот как выглядит твердая фаза:
будет гексатическая фаза (жидкокристаллическая, снэпшоты не могу найти).
Трудная часть - различить жидкость и газ. вы должны заметить, что при достаточно высокой температуре нет различия между газом и жидкостью (выше критической точки), поэтому, если вы просто посмотрите на расстояние между частицами, вы не сможете сказать, что такое жидкость, а что газ.
Если вы посмотрите на температуру ниже критической температуры, то в зависимости от прилагаемого давления вы можете увидеть газ, жидкость или капли жидкости внутри газа. когда я говорю о каплях жидкости в газе, я имею в виду, что у вас будут скопления частиц и несколько свободных частиц, движущихся между ними — см. рисунок в области жидкого газа.
если вы находитесь выше критической температуры, нет различия между газом и жидкостью, мы называем это жидкостью (см. однофазную область на рисунке).
Если вы хотите различать газ\жидкость\флюид - есть способ:
что вы можете сделать, это посмотреть на распределение плотности. разделите снимок вашей симуляции на прямоугольники (вам нужно будет поиграть с размерами, в зависимости от количества частиц и плотности), как на этой картинке:
спросите, какова плотность в каждом ящике, и поместите результаты в гистограмму. сделайте это для всех шагов в вашей симуляции и усредните результаты. см. эту статью для более подробной информации.
если у вас есть жидкость, смешанная с газом, вы увидите две плотности:
по мере увеличения температуры две плотности будут сближаться, пока не станут одной. это означает, что вы преодолели критическую температуру и теперь находитесь выше критической температуры.
еще одно замечание: для определения точки замерзания и точки испарения могут потребоваться разные методы моделирования. по мере приближения к критической точке вам потребуется гораздо больше шагов, чтобы уравновесить симуляцию, поэтому вы, вероятно, захотите использовать ансамбль Гиббса, чтобы избежать этого.
K7PEH
Джабаву Адамс
K7PEH