Розрахунок середньоквадратичних зміщень презентация

N

M

2. Для кожного початкового моменту розраховується 3N різниць DX=X(I)-X(0), DY=Y(I)-Y(0), DZ=Z(I)-Z(0) для I=1,N і для кожної перевіряються граничні умови.

3. Результат нормується на кількість усереднень (число частинок сорту A,B та число різних функцій для кожного I)

Модель води SPC/E:
Жорсткі молекули
rOH = 1 A
HOH=109.47o
qO = -0.8476 |e|
qH = +0.4238 |e|
Потенціали: LJ+Coulomb

Ο

Η

Η

Введемо Лагранжіан системи з α різними звязками

- невідомі множники Лагранжа

Рівняння руху:

Додаткова сила через наявність звязку:

Для знаходження λ використовуємо :

Рівняння руху:

Множник Лагранжа:

Тобто, додаткова сила що діє внаслідок накладеного звязку:

Застосовується не одночасно до всіх звязків, а послідовно – звязок за звязком

Тому: