L9 · 26 دقيقة
المكاملات والخطوة الزمنية
Verlet وقيود الهيدروجين و2 fs.
أهداف التعلّم
- شرح فكرة Verlet: مواضع من قوى دون اشتقاق سرعات كخطوة أولى ضرورية.
- ربط الخطوة الزمنية بأسرع حركة متبقية بعد القيود.
المكامل هو حساب التفاضل للمواضع
نعرف a = F/m في الزمن t، ونريد r عند t+Δt. مكامل Verlet (وصيغه: leap-frog، velocity-Verlet) يستخدم أن التسارع يدخل التوسعة التايلورية كحد Δt²، بينما حد السرعة الفردي يُلغى بين خطوة أمامية وخلفية. الناتج مستقر طويلًا إن كانت Δt صغيرة أمام أسرع اهتزاز، وحافظ للطاقة في NVE أفضل من Euler الساذج الذي ينفجر. اختيار المكامل في الممارسة غالباً ما يفرضه المحرّك (GROMACS leap-frog شائع). فهمك يُظهر في اختيار Δt والقيود، لا في إعادة اشتقاق الورقة الأصلية.
فكرة Verlet للمواضع
الموضع اللاحق من الموضع الحالي والسابق والتسارع الحالي. السرعة ليست صريحة في هذه الصيغة.
- \mathbf{r}(t)متجه موضع الذرة عند الزمن t · nm
- \Delta tالخطوة الزمنية · ps أو fs
- \mathbf{a}(t)التسارع من F/m · nm ps⁻²
خطأ محلي من رتبة Δt⁴ في هذه العائلة (حسب الصيغة)، لكن الخطأ يتراكم. Δt كبيرة تجعل الحدود المهملة في تايلور غير مهملة: النظام يسخن أو ينهار. القيود (SHAKE، LINCS) تُحل بعد أو ضمن الخطوة لتعيد روابط H إلى طولها.
تطبيق في التصميم الدوائي: اختبار صحة الخطوة: محاكاة NVE قصيرة ومراقبة انحراف الطاقة الكلية. انحراف منتظم متصاعد يعني Δt أخشن من الفيزياء.
صيغ تلتقي في الممارسة
- Verlet الأصلي: مواضع؛ السرعات تُستنتج عند الحاجة.
- leap-frog: سرعات نصف خطوة بعيداً عن المواضع. شائع في GROMACS.
- velocity-Verlet: مواضع وسرعات متزامنة، مريح لمنظمات حرارة معيّنة.
- Euler: لا تستخدمه لإنتاج MD. غير مستقر على هذا السطح.
لماذا؟ — لماذا لا نأخذ Δt = 5 fs «لنوفر الوقت» بعد تقييد الهيدروجين؟
لأن الماء وزوايا الثني ودوران مجموعات الميثيل ما زالت سريعة. المكامل يخطئ في تلك الحدود فيضخ طاقة (تسخين اصطناعي) أو يكسر قيوداً. التوفير الوهمي يُشترى بفيزياء تالفة. إن رغبت خطوة أطول، فالحجة تجريبية: استقرار NVE، لا شعار الإنتاجية.
ماذا لو؟ — ماذا لو كانت الطاقة في NVE تتذبذب حول متوسط ثابت، فانتقلت إلى إنتاج NPT بخطوة أكبر مرتين؟
اختبار NVE كان لـ Δt الأصلية. مضاعفة الخطوة تُعيد إدخال خطأ التكامل. منظم الحرارة قد يُخفي التسخين بامتصاص الطاقة الزائدة، فتبدو المحاكاة «مستقرة» وهي تُدفأ رقمياً. لا تُغيّر Δt بين الاختبار والإنتاج دون إعادة الاختبار.
تمرين
صف تجربة ذهنية تُظهر لماذا Euler ينفجر بينما Verlet يبقى على بندول توافقي تقريباً.
مواضيع مرتبطة