ЗАДАЧА 16.16. Моделирование твердых стержней методом Моите-Карло

а.         Модифицируйте программу hard_disk так, чтобы она годилась длясистемы твердых стержней. Задайте L = 12 и Л/ = 10. Как согласу-ется плотность с максимальной возможной плотностью? Выберите на-чальные положения так, чтобы оии образовывали одномерную сетку.Чему равен коэффициент принятия для варианта dxmax = 0.1? Сколь-ко приблизительно требуется шагов Монте-Карло на частицу длядостижения равновесия? Вычислите парную корреляционную функциюg(r) с г = \х\. Не забудьте, что г измеряется в единицах о*.

б.         Постройте график g(r) как функцию от г. Почему при значенияхг < 1 функция g(r) = 0? Почему значения g(r) для г > L/2 не име-ют смысла? Какой физический смысл имеют пики на кривой g(r)lЭкстраполируя полученные для g(r) результаты на точку г = 1 ииспользуя формулы (16.31), определите давление.

в.         Получите функцию g(r) для нескольких меньших значений плот-ности, взяв за основу какую-нибудь равновесную конфигурацию изпредыдущего варианта и увеличивая L. Как меняются размеры и рас-положение пиков на кривой g(r)?