розрахунку, яка відповідає температурі і щільності, близьким до необхідних. Наступна підпрограма запам'ятовує останню конфігурацію пропускається варіанту і може бути використана в якості останньої підпрограми, що викликається програмою md. p> Наступні інструкції можна додати в підпрограму initiаl, з тим щоб можна було скористатися якоюсь колишньою рівноважної конфігурацією в якості початкової конфігурації нового розрахунку.
Необхідно тільки один раз провести моделювання молекулярної динаміки, вирушаючи від системи з довільною конфігурацією. Наступні розрахунки для різних станів можна починати, модифікуючи умови, сформовані на момент завершення відповідного попереднього розрахунку. Питання про те, як змінювати температуру системи при незмінній щільності, розглянуто в задачі. p> ЗАВДАННЯ 1. Якісні властивості рідини і газу
а. У всіх пунктах даного завдання для отримання своєї початкової конфігурації використовуйте наступні інструкції dаtа і rеаd. Інструкція RЕАD читає елементи даних з інструкцій DАТА і привласнює їх значення відповідним змінним. Кожна інструкція READ читає черговий елемент з інструкції DАТА. Напишіть коротку програму, в якій використовуються інструкції RЕАD і DАТА, і переконайтеся, що ви правильно розумієте роботу цих інструкцій. p align="justify"> Виберіть N = 16, Lх = Lу = 6, rsса1е = 1 і? t = 0.02. Чому дорівнює наведена щільність? Виберіть nsnар = 20 і проведіть моделювання протягом принаймні 200 кроків за часом. Чому дорівнює початкова температура системи? Який характер ваших знімків? Як довго система досягає рівноваги? Який ваш критерій рівноваги? Обчисліть повну енергію системи, температуру і тиск. Чи не враховуйте при усередненні за часом нерівноважні конфігурації. p align="justify"> б. Як було знайдено в п. В«аВ», повна енергія негативна. Який сенс має знак повної енергії? Повторіть моделювання з тими ж початковими умовами (rscа1е = 1), але з Lх = Ly = 30. Чому дорівнює повна енергія в цьому випадку? Тому ж, що й у п. В«аВ»? Якщо ні, то чому? Опишіть характер отриманих знімків. Заповнюють Чи частинки ящик рівномірно за час порядку 200 кроків або у них спостерігається тенденція до утворення В«крапельокВ»? p align="justify"> в. Збільште початкова відстань між частинками. Використовуйте Lx = Lу = 30 і rsсаlе = 2. Чому дорівнює наведена щільність системи? Оцініть початкову щільність крапельки. Яка вийшла повна енергія-негативна або позитивна? Якщо величину rsсаlе зробити досить великий, крапелька повинна зрештою В«випаруватисяВ», навіть якщо повна енергія негативна. У реальній системі крапелька знаходився б в рівновазі зі своїм паром і обидві фази існували б спільно. У міру зменшення щільності початковій крапельки все більше частинок усередині крапельки буде випаровуватися і розмір крапельки буде скорочуватися. Зрозуміло, потрібно з обережністю ставитися до будь-яких висновків про структуру системи, заснованим на моделюванні тільки з 16 частинками. p align="justify"> р. ...
