ديناميكا براونية

يمكن أن تُستعمل الديناميكا البراونية أو الديناميك البراوني قالب:إنج (و يسمى اختصاراً BD) لوصف حركة الجزيئات في المحاكاة الجزيئة. هي نسخة مبسطة لديناميكا لانجيفين Langevin dynamics وتنسجم مع الحد التي لا يأخذ فيها متوسط التسارع مكاناً ثابتاً عند تشغيل المحاكاة. يمكن أن تُوصف هذا التقريب على أنه ديناميكا لانجيفين 'مفرطة التخامد overdamped', أو, كديناميكا لانجيفين من دون عطالة.

في ديناميكا لانجيفين, تكون معادلة الحركة هي

$M\ddot{X}=-\mathrm{\nabla }U(X)-\gamma M\dot{X}+\sqrt{2\gamma k_{B}TM}R(t)$

حيث أن $U(X)$ هي كمون تفاعل الجسيم; و$\mathrm{\nabla }$ هي مؤثر التدرج gradient operator عندما تكون $\mathrm{\nabla }U(X)$ هي القوة المحسوبة من كمونات تفاعل الجسيم; والنقطة هي عملية الضرب وهي مشتقة عندما تكون $\dot{X}$ هي السرعة و$\ddot{X}$ هي التسارع; T هي الحرارة, k_B هي ثابت بولتزمان; و$R(t)$ هي التفاعل الغاوسي المستقر ذات الارتباط-دلتا بالمعنى-الصفري، وتكون

$⟨R(t)⟩=0$

$⟨R(t)R(t^{\prime })⟩=\delta (t-t^{\prime })$

في الديناميكا البراونية, لا يوجد تسارع مفترض يمكن أخذها في الاعتبار. وهكذا يُهمل الحد $M\ddot{X}(t)$, ومجموع هذه الحدود هو الصفر.

$0=-\mathrm{\nabla }U(X)-\gamma M\dot{X}+\sqrt{2\gamma k_{B}TM}R(t)$

بتعريف $\zeta =\gamma M$, وباستعمال علاقة أينشتاين Einstein relation, $D=k_{B}T/\zeta$, تكون أسهل كتابة المعادلة بهذه الطريقة

$\dot{X}(t)=-\mathrm{\nabla }U(X)/\zeta +\sqrt{2D}R(t)$

انظر أيضاً

• حركة براونية

وصلات خارجية

• Course on Brownian and Langevin Dynamics