2014 dxdy logo

Научный форум dxdy

Математика, Физика, Computer Science, Machine Learning, LaTeX, Механика и Техника, Химия,
Биология и Медицина, Экономика и Финансовая Математика, Гуманитарные науки




 
 Расчет равновесных характеристик методом Монте-Карло
Сообщение28.04.2016, 17:14 
Аватара пользователя
Добрый день,
Допустим, что необходимо посчитать значение какой-либо величины в случае NVT ансамбля. Пусть мы просемплировали пространство состояний при помощи алгоритма Метрополиса-Гастингса с вероятностью перехода в новую конфигурацию $p=\min\{1, \exp(\frac{\Delta E}{kT})\}$, где $\Delta E $ - разница в энергиях между новой и старой конфигурацией. Теперь вопрос: можно ли потом еще раз перевзвесить по Больцману конфигурации в полученной траектории, если это самое Больцмановское распределение и так участвовало в получении этой траектории, и если да, то на что можно в этом случае сослаться?

 
 
 
 Re: Расчет равновесных характеристик методом Монте-Карло
Сообщение10.05.2016, 20:03 
Аватара пользователя
К предыдущему сообщению: так делать, видимо, не стоит...

Но теперь новый вопрос. Методом Метрополиса (с вышеуказанной вероятностью принятия шага) считается определенного вида система (молекула). Сэмплируется наблюдаемые в эксперименте свойства при NVT ансамбле. Начало траектории не выкидывается (система стартует с самой выгодной конфигурации и предполагается существенная локализация системы в определенном многомерном "колодце" при данной температуре). Далее на основе этой траектории строится теоретическая модель и сравнивается с экспериментом. Предполагается, что до определенного момента (пока система все же не перескочит в соседний "колодец") рост длины траектории должен давать улучшение согласие с экспериментом. Но, это улучшение наблюдается далекооооо не всегда.
В чем может быть причина такого поведения, и что нужно делать, чтобы от этого избавиться?

 
 
 [ Сообщений: 2 ] 


Powered by phpBB © 2000, 2002, 2005, 2007 phpBB Group