В Метод распределения с перекрытием был представлен Чарльз Х. Беннетт[1] для оценки химический потенциал.
Теория
Для двух N системы частиц 0 и 1 с функция распределения
и
,
из ![{displaystyle F (N, V, T) = - k_ {B} Tln Q}](https://wikimedia.org/api/rest_v1/media/math/render/svg/ad2c9ecb61b922bc5d6c6dbb76dd8dd9d230da5d)
получить термодинамическая свободная энергия разница ![{displaystyle Delta F = -k_ {B} Tln (Q_ {1} / Q_ {0}) = - k_ {B} Tln ({frac {int ds ^ {N} exp [- eta U_ {1} (s ^ {N})]} {int ds ^ {N} exp [- eta U_ {0} (s ^ {N})]}})}](https://wikimedia.org/api/rest_v1/media/math/render/svg/0609baade70017808cdf8055ed4e24fb6fe6c80a)
Для каждой конфигурации, посещенной во время этой выборки системы 1, мы можем вычислить потенциальную энергию U как функцию пространства конфигурации, а разность потенциальной энергии равна
![{displaystyle Delta U = U_ {1} (s ^ {N}) - U_ {0} (s ^ {N})}](https://wikimedia.org/api/rest_v1/media/math/render/svg/8e272b734ad27430a396b1c6b844f179ae92dddc)
Теперь постройте плотность вероятности потенциальной энергии из приведенного выше уравнения:
![{displaystyle p_ {1} (Delta U) = {frac {int ds ^ {N} exp (- eta U_ {1}) delta (U_ {1} -U_ {0} -Delta U)} {Q_ {1} }}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/59ee9385ddef0b7358dffa1b1b1a8380545ae323)
в которой
является конфигурационной частью статистической суммы
![{displaystyle = {frac {Q_ {0}} {Q_ {1}}} exp (- eta Delta U) {frac {int ds ^ {N} exp (- eta U_ {0}) delta (U_ {1} - U_ {0} -Delta U)} {Q_ {0}}} = {frac {Q_ {0}} {Q_ {1}}} exp (- eta Delta U) p_ {0} (Delta U)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/bd842ccfb5b5afc1defdaf3a99aed0edf44f2363)
поскольку
![{displaystyle Delta F = -k_ {B} Tln (Q_ {1} / Q_ {0})}](https://wikimedia.org/api/rest_v1/media/math/render/svg/7c822ec7352aa4003d2faa70857eca524de52b65)
![{displaystyle ln p_ {1} (Delta U) = eta (Delta F-Delta U) + ln p_ {0} (Delta U)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/2e3d02bccee397896b97fb5a2e3d2461a115622d)
теперь определите две функции:
![{displaystyle f_ {0} (Delta U) = ln p_ {0} (Delta U) - {frac {eta Delta U} {2}} f_ {1} (Delta U) = ln p_ {1} (Delta U) + {frac {eta Delta U} {2}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/46d383e74ead281e5992b321dc566c2a3b3fbc5d)
таким образом, что
![{displaystyle f_ {1} (Delta U) = f_ {0} (Delta U) + eta Delta F}](https://wikimedia.org/api/rest_v1/media/math/render/svg/aba17ffaf6475e7674334740bda262be8f4d25c9)
и
можно получить, установив
и ![f _ {{0}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/6423b30a4c5770c59b5ab92dcb4ce378755440ab)
Рекомендации