L9 · 32 دقيقة
NVT وNPT ومنظمات الحرارة والضغط
بيرندسن ليس للقياس الإنتاجي بلا وعي.
أهداف التعلّم
- تعريف N,V,T,P فيزيائياً وتمييز NVT عن NPT وعن NVE.
- مقارنة Berendsen وNosé–Hoover وLangevin وv-rescale، وBerendsen وParrinello–Rahman وMC للضغط.
ماذا تعني الحروف؟
N عدد الجسيمات (الذرات أو المواقع في النموذج): في محاكاة مغلقة لا يُخلق بروتين فجأة. V حجم الصندوق: يتغير إن سمح الباروستات. T درجة الحرارة: مقياس لطاقة الحركة عبر معدّل يرتبط بـ kB، تُضبط بمنظم حرارة لا بـ «ضبط الرقم في الإعداد ونسيانه». P الضغط: قوة على وحدة مساحة من الاصطدامات والتفاعلات مع جدران وهمية للصندوق؛ يُضبط بتغيير الحجم (أو محاولات MC للحجم). E الطاقة الكلية تُحفظ في NVE (الميكروكانوني) إن أحسن المكامل. H الإنثالبي ذات صلة بـ NPT. المختبر البيولوجي أقرب إلى T وP ثابتين (حمام مائي، 1 bar) لذا الإنتاج غالباً NPT بعد اتزان. NVE أداة تشخيص للمكامل.
| NVT | NPT | |
|---|---|---|
| الثوابت | عدد الجسيمات، الحجم، درجة الحرارة. | عدد الجسيمات، الضغط، درجة الحرارة. |
| ما يتحرك | السرعات تُعدَّل لتطابق T؛ الحجم ثابت. | الحجم يتقلب ليطابق P، والسرعات لتطابق T. |
| الاستخدام النموذجي | اتزان حراري بعد التسخين؛ أحياناً إنتاج إن كان الحجم مُعايَراً. | اتزان كثافة ثم إنتاج يشبه المختبر. |
منظمات الحرارة
- Berendsen: إعادة تدريج ضعيفة للسرعات نحو T الهدف. يُخمِد التقلبات، لا يولّد ensemble قانوني صارم. مقبول لاتزان قصير، مرفوض كإنتاج دون وعي.
- Nosé–Hoover: درجة حرية ممتدة (حمام وهمي). يقارب NVT الصحيح، قد يتذبذب، وحذر مع أنظمة ضعيفة الارتباط.
- Langevin: احتكاك + ضجيج عشوائي يحقق توازن تذبذب-تبديد. يعاين NVT، ويُخمِد الديناميكا الهيدروديناميكية (الانتشار يتأثر).
- v-rescale (Bussi–Donadio–Parrinello): إعادة تدريج مع حد عشوائي يُصلح توزيع الطاقة الحركية. شائع في GROMACS للإنتاج.
منظمات الضغط
- Berendsen: تدريج حجم ضعيف. يُقرّب الكثافة بسرعة، يُخمِد تقلبات الضغط. اتزان لا إنتاج صارم.
- Parrinello–Rahman: ديناميكا ممتدة للصندوق (قد تسمح بشكل متغير). مناسب لإنتاج NPT بعد استقرار. إن وُضع قبل اتزان الحجم فقد يتذبذب بعنف.
- MC barostat: محاولات عشوائية لتغيير الحجم تُقبل بـ Metropolis. شائع في OpenMM وبعض بروتوكولات AMBER. لا ديناميكا حجم ممتدة.
لماذا؟ — لماذا لا يبدأ الإنتاج بـ Parrinello–Rahman مباشرة بعد التقليل؟
لأن الكثافة بعيدة عن الاتزان والضغط لحظياً هائل التقلب. الباروستات الممتد يتفاعل مع ذلك بذبذبة حجم قد تمزق النظام. بروتوكول شائع: تقليل → تسخين NVT بقيود → NPT بباروستات لطيف (Berendsen) حتى تستقر الكثافة → إنتاج بـ v-rescale + Parrinello–Rahman أو ما يعادل. الترتيب حماية للفيزياء.
ماذا لو؟ — ماذا لو استخدمت Langevin باحتكاك عالٍ ثم فسّرت انتشار الليجند أو زمن ارتباط كقيمة فيزيائية؟
الاحتكاك يُبطئ النقل الهيدروديناميكي عمداً. الأزمنة لم تعد أزمنة النموذج المحافظ ولا أزمنة التجربة. Langevin جيد للمعاينة الثرموديناميكية (NVT)؛ سيئ للادعاء الحركي ما لم يُعاير الاحتكاك ويُعلن. للثرموديناميكا يمكنك أن تدفع المعاينة؛ للحركية أنت تشتري نموذجاً حركياً آخر.
تمرين
اقترح سلسلة ensemble لمحاكاة بروتين–ليجند في ماء عند 1 bar و310 K، وسمِّ أين Berendsen مقبول وأين ترفضه.
مواضيع مرتبطة