Расчет критического размера и свойств микрокапель методом Монте-Карло в обобщенном ансамбле

Предлагаются два варианта метода Монте-Карло для расчета свойств микрокластеров с использованием обобщенного $\mu pT$-ансамбля. Первый способ предполагает генерирование отрезков неравновесных марковских цепей при заданных давлении, температуре и химическом потенциале и начальном числе частиц. Проведены расчеты для двух начальных условий в системе леннард-джонсовских частиц. Второй (стационарный) вариант метода отличается от первого тем, что число частиц в системе в каждой цепи может меняться только в пределах $n$, $n+1$ и, следовательно, генерируется стационарная марковская цепь. Рассчитывая отношение вероятностей нахождения системы в состоянии с тем и другим числом частиц, можно вычислить химический потенциал системы, а затем термодинамический потенциал Гиббса. Аналогичным образом, используя $\mu VT$- ансамбль, можно рассчитать свободную энергию Гельмгольца. Результаты, полученные стационарным методом, сравниваются с численными экспериментами более ранних работ, где используются другие способы расчета.