[i] <50) or (y [i]> 350) vy [i]: =-vy [i]; {від стінок} j: = 1 to n do (i <> j) and (abs (x [i]-x [j]) <= 2) and (abs (y [i]-y [j]) <= 2) [i]: = x [i] + vx [i] * dt;
vy [j]: =-vy [j];;;; (5000), (3); (48, 46, 352, 46); (48, 48, 48, 352) ; (352, 46, 352, 352); (352, 352, 48, 352); (1,11); (48,48,352,352); (0); i: = 1 to N do circle (round (x [ i]), round (y [i]), 2); KeyPressed;;.
5. Потенціал міжмолекулярної взаємодії
Насамперед нам необхідно визначити модель системи, яку ми бажаємо моделювати. Оскільки ми хочемо зрозуміти якісні властивості систем багатьох частинок, підемо на спрощення завдання, припускаючи, що динаміку можна вважати класичною, а молекули - хімічно інертними кульками. Ми припускаємо також, що сила взаємодії будь-яких двох молекул залежить тільки від відстані між ними. У цьому випадку повна потенційна енергія U визначається сумою двох часткових взаємодій:
U = V (r12) + V (r13) + ... + V (r23) + ... == (5.1)
де V (r ij ) залежить тільки від абсолютної величини відстані r ij між частинками i і j. Парне взаємодія виду (5.1) відповідає В«простимВ» рідин, наприклад рідкий Аргон.
В принципі форму V (r) для електрично нейтральних атомів можна побудувати шляхом детального розрахунку, що базується на основних законах квантової механіки. Такий розрахунок дуже важкий і, крім того, зазвичай буває цілком достатньо в якості V (r) вибрати просту феноменологічну формулу. Найбільш важливими особливостями V (r) для простих рідин є сильне відштовхування для малих r і слабке тяжіння на великих відстанях. Відштовхування при малих r обумовлено правилом заборони. Інакше кажучи, якщо електронні хмари двох атомів перекриваються, деякі електрони повинні збільшити свою кінетичну енергію, щоб перебувати в різних квантових станах. Сумарний ефект виражається в відштовхуванні між електронами, званому відштовхуванням кора . Слабке тяжіння при великих r обумовлено головним чином взаємної поляризацією кожного атома; результуюча сила тяжіння називається силою Ван-дер-Ваальса .
Однією з найбільш уживаних феноменологічних формул для V (r) є потенціал Леннарда - Джонса
(5.2)
Графік потенціалу Ленарда-Джонса зображений на рис. 6.1. Хоча залежність r -6 у формулі (6.2) отримана теоретично, залежність r -12
В
Рис. 5.1 Графік потенціалу Ленарда-Джонса. Зазначимо, що r вимірюється в одиницях і V (r)-в одиницях
...
