ديناميكا جزيئية
الديناميكا الجزيئية (بالإنجليزية: Molecular dynamics (MD))، هي علم قائم على محاكاة حاسوبية للجزيئات والمركبات معتمدة على قوانين نيوتن للحركة في توقع حركة الجزيئات وسرعة تحركها واماكن وجودها بعد فترات زمنية قصيره جدا وهو عمل يتطلب حلول لمعادلات نيوتن لاعداد هائلة من الجزيئات ولفترات زمنية قصيره جدا ولضمان حساب هذه التوقعات بشكل صحيح أو اقرب للصحة يجب أخذ عدة امور عند حساب قوانين نيوتن مثل قوى التجاذب بين الجزيئات واهتزاز روابط الجزي وانثناء الرووابط حول محورها وتذبذب الزاوية بين ذرات الجزي وايضا تذبذب الزوايا بين مستويات الجزئ أو ما يعرف بالزاوية الزوجية وامور أخرى كثيرة

ولاننا في هذه العملية نقوب بحساب قوانين نيوتن وقوانين أخرى تتعلق بالمحاكاة لاعداد كبيرة جده من الجزيئات هذا يجعل العملية معقدة جدا وشبه مستحيلة لذلك يتم الاستعانة ببعض الفرضيات وطرق احصائية لتساعدنا في جعل هذه الحسابات ممكنة من الناحية العملية.

الديناميكات الجزيئية هي طريقة محاكاة حاسوبية لتحليل الحركات الفيزيائية للذرات والجزيئات. يُسمح للذرات والجزيئات بالتفاعل لفترة زمنية محددة، مما يعطي نظرة على «التطور» الديناميكي للنظام. في الإصدار الأكثر شيوعًا، يتم تحديد مسارات الذرات والجزيئات عن طريق الحل العددي لمعادلات نيوتن للحركة لنظام من الجسيمات المتفاعلة، حيث تُحسب القوى بين الجسيمات وطاقاتها الكامنة غالبًا باستخدام الكمون بين الذرات أو مجالات القوة الجزيئية. يتم تطبيق الطريقة في الغالب في الفيزياء الكيميائية وعلوم المواد والفيزياء الحيوية.

نظرًا لأن الأنظمة الجزيئية تتكون عادةً من عدد كبير من الجسيمات، فمن المستحيل تحديد خصائص هذه الأنظمة المعقدة بشكل تحليلي؛ تتغلب محاكاة الديناميكا الجزيئية على هذه المشكلة باستخدام الطرق العددية. ومع ذلك، فإن عمليات محاكاة الديناميكا الجزيئية الطويلة غير مشروطة رياضيًا، مما يؤدي إلى حدوث أخطاء تراكمية في التكامل العددي يمكن تقليلها من خلال الاختيار المناسب للخوارزميات والمعلمات، ولكن لا يتم التخلص منها تمامًا.

بالنسبة للأنظمة التي تتبع فرضية إرغوديك، يمكن استخدام تطور محاكاة ديناميكية جزيئية واحدة لتحديد الخصائص الديناميكية الحرارية العيانية للنظام: المتوسطات الزمنية لنظام ergodic تتوافق مع متوسطات المجموعة الدقيقة. كما أطلق على الديناميكا الجزيئية أيضًا اسم «الميكانيكا الإحصائية بالأرقام» و «رؤية لابلاس للميكانيكا النيوتونية» للتنبؤ بالمستقبل من خلال تحريك قوى الطبيعة  والسماح بإلقاء نظرة ثاقبة على الحركة الجزيئية على نطاق ذري.

تاريخ 
تم تطوير الديناميكا الجزيئية في الأصل في أوائل الخمسينيات من القرن الماضي، بعد النجاحات السابقة في محاكاة مونت كارلو، والتي يعود تاريخها إلى القرن الثامن عشر، في مشكلة إبرة بوفون على سبيل المثال، ولكن تم نشرها في الميكانيكا الإحصائية في مختبر لوس ألاموس الوطني بواسطة روزنبلوث ومتروبوليس فيما يعرف اليوم بخوارزمية متروبوليس-هاستينغز. يعود الاهتمام بالتطور الزمني لأنظمة الجسم N إلى القرن السابع عشر، بدءًا من نيوتن، واستمر في القرن التالي إلى حد كبير مع التركيز على الميكانيكا السماوية وقضايا مثل استقرار النظام الشمسي. تم تطوير العديد من الأساليب العددية المستخدمة اليوم خلال هذه الفترة الزمنية، والتي سبقت استخدام أجهزة الكمبيوتر؛ على سبيل المثال، خوارزمية التكامل الأكثر شيوعًا المستخدمة اليوم، خوارزمية تكامل فيرليت، تم استخدامها في وقت مبكر من عام 1791 بواسطة جان باتيست جوزيف ديلامبر. يمكن اعتبار الحسابات العددية باستخدام هذه الخوارزميات الديناميكا الجزيئية «يدويًا».

في وقت مبكر من عام 1941، تم تنفيذ تكامل معادلات الحركة للعديد من الأجسام مع أجهزة الكمبيوتر التناظرية. قام البعض بعمل كثيف العمالة لنمذجة الحركة الذرية من خلال بناء نماذج فيزيائية، على سبيل المثال، باستخدام المجالات العيانية. كان الهدف هو ترتيبها بطريقة تماثل بنية السائل واستخدامه لفحص سلوكه. برنال قال في عام 1962: «. . . أخذت عددًا من الكرات المطاطية وغرستهم معًا بقضبان بأطوال مختلفة تتراوح من 2.75 إلى 4 بوصات. حاولت القيام بذلك في المقام الأول بشكل عرضي قدر الإمكان، حيث أعمل في مكتبي الخاص، وأتعرض للمقاطعة كل خمس دقائق أو نحو ذلك ولا أتذكر ما فعلته قبل المقاطعة.» 

بعد اكتشاف الجسيمات المجهرية وتطوير أجهزة الكمبيوتر، توسع الاهتمام إلى ما وراء إثبات أنظمة الجاذبية إلى الخصائص الإحصائية للمادة. في محاولة لفهم أصل اللارجعة، اقترح فيرمي في عام 1953، ونشر في عام 1955،  استخدام مانياك، أيضًا في مختبر لوس ألاموس الوطني، لحل التطور الزمني لمعادلات الحركة للعديد من: يخضع نظام الجسم للعديد من خيارات قوانين القوة؛ اليوم، يُعرف هذا العمل الأساسي بمشكلة فيرمي - باستا - أولام - تسينغو. يظهر التطور الزمني للطاقة من العمل الأصلي في الشكل إلى اليمين.

في عام 1957، استخدم بيرني ألدر ووينرايت  حاسوب آي بي إم 704 لمحاكاة التصادمات المرنة تمامًا بين المجالات الصلبة. في عام 1960، ربما في أول محاكاة واقعية للمادة، جيبسون وآخرون. محاكاة التلف الإشعاعي للنحاس الصلب باستخدام نوع ولد ماير من التفاعل البغيض مع قوة سطحية متماسكة. في عام 1964، نشر عبد الرحمن  محاكاة للأرجون السائل تستخدم كمون لينارد جونز. مقارنة حسابات خصائص النظام، مثل معامل الانتشار الذاتي، بشكل جيد مع البيانات التجريبية.

مجالات التطبيق وحدودها 
تم استخدام طريقة الديناميكا الجزيئية لأول مرة في الفيزياء النظرية، وقد اكتسبت شعبية في علم المواد بعد ذلك بوقت قصير، ومنذ السبعينيات أصبحت شائعة أيضًا في الكيمياء الحيوية والفيزياء الحيوية. يستخدم الديناميكا الجزيئية بشكل متكرر لتحسين الهياكل ثلاثية الأبعاد للبروتينات والجزيئات الكبيرة الأخرى بناءً على القيود التجريبية من علم البلورات بالأشعة السينية أو التحليل الطيفي بالرنين المغناطيسي النووي. في الفيزياء، يتم استخدام الديناميكا الجزيئية لفحص ديناميكيات الظواهر على المستوى الذري التي لا يمكن ملاحظتها مباشرة، مثل نمو الأغشية الرقيقة وزرع الأيونات الفرعية، وكذلك لفحص الخصائص الفيزيائية للأجهزة النانوية التي لم يتم إنشاؤها أو لا يمكن إنشاؤها بعد. في الفيزياء الحيوية والبيولوجيا الهيكلية، يتم تطبيق الطريقة بشكل متكرر لدراسة حركات الجزيئات الكبيرة مثل البروتينات والأحماض النووية، والتي يمكن أن تكون مفيدة لتفسير نتائج بعض التجارب الفيزيائية الحيوية ونمذجة التفاعلات مع الجزيئات الأخرى، كما هو الحال في الإرساء الترابطي. من حيث المبدأ، يمكن استخدام الديناميكا الجزيئية للتنبؤ من البداية بهيكل البروتين عن طريق محاكاة طي سلسلة البولي ببتيد من ملف عشوائي.

يمكن اختبار نتائج محاكاة الديناميكا الجزيئية من خلال المقارنة مع التجارب التي تقيس الديناميكيات الجزيئية، والتي تعد طريقة شائعة منها هي التحليل الطيفي بالرنين المغناطيسي النووي. يمكن اختبار تنبؤات الهيكل المشتقة من الديناميكا الجزيئية من خلال تجارب على مستوى المجتمع في التقييم النقدي لتوقع بنية البروتين (CASP)، على الرغم من أن الطريقة قد حققت نجاحًا محدودًا تاريخيًا في هذا المجال. كتب مايكل ليفيت، الذي شارك في جائزة نوبل جزئياً لتطبيق الديناميكا الجزيئية على البروتينات، في عام 1999 أن المشاركين في التقدير التوقعي لهيكل البروتين لم يستخدموا الطريقة عادةً بسبب «... إحراج مركزي للميكانيكا الجزيئية، أي تقليل الطاقة أو الديناميات الجزيئية بشكل عام يؤدي إلى نموذج أقل شبهاً بالهيكل التجريبي.»  أدت التحسينات في الموارد الحسابية التي تسمح بمسارات الديناميكا الجزيئية أكثر وأطول، جنبًا إلى جنب مع التحسينات الحديثة في جودة معلمات مجال القوة، إلى بعض التحسينات في كل من تنبؤ الهيكل وصقل نموذج التماثل، دون الوصول إلى نقطة الفائدة العملية في هذه المجالات؛ يحدد العديد معلمات مجال القوة كمجال رئيسي لمزيد من التطوير.

تم الإبلاغ عن الديناميكا الجزيئية محاكاة ل الخاصة الدوائية تطوير وتصميم الأدوية. على سبيل المثال، بينتو وآخرون نفذوا عمليات محاكاة الديناميكا الجزيئية لمجمعات (Bcl-Xl) لحساب متوسط مواضع الأحماض الأمينية الحرجة المشاركة في ربط الترابط. من ناحية أخرى، كارلسون وآخرون. نفذت محاكاة الديناميكيات الجزيئية لتحديد المركبات التي تكمل المستقبل مع التسبب في الحد الأدنى من اضطراب التشكل والمرونة في الموقع النشط. تم تراكب لقطات من البروتين على فترات زمنية ثابتة أثناء المحاكاة لتحديد مناطق الربط المحفوظة (المحفوظة في ثلاثة إطارات على الأقل من أحد عشر إطارًا) لتطوير عالم الأدوية. سبيراكيس وآخرون يعتمد على سير عمل محاكاة الديناميكا الجزيئية، وبصمات الأصابع للروابط والبروتينات (FLAP) وتحليل التمييز الخطي لتحديد أفضل توافق بروتين يجند للعمل كقوالب صيدلانية استنادًا إلى تحليل (ROC) بأثر رجعي لعقاقير الصيدلانية الناتجة. في محاولة لتحسين نمذجة اكتشاف الأدوية القائمة على الهيكل، مقابل الحاجة إلى العديد من المركبات النموذجية، اقترح هاتمال مزيجًا من محاكاة الديناميكا الجزيئية وتحليل الاتصالات بين الجزيئات لمستقبلات ليجيند لتمييز الاتصالات الحرجة بين الجزيئات (تفاعلات الربط) من تلك الزائدة عن الحاجة في مركب واحد بروتين يجند. يمكن بعد ذلك تحويل جهات الاتصال الحرجة إلى نماذج صيدلانية يمكن استخدامها للفحص الافتراضي.

ترتبط حدود الطريقة بمجموعات المعلمات المستخدمة، ومجالات قوة الميكانيكا الجزيئية الأساسية. تشغيل واحد لمحاكاة الديناميكا الجزيئية يعمل على تحسين الطاقة الكامنة، بدلاً من الطاقة المجانية للبروتين، مما يعني أن جميع المساهمات الحتمية في الاستقرار الديناميكي الحراري لبنية البروتين مهملة، بما في ذلك الانتروبيا المطابقة لسلسلة البولي ببتيد (العامل الرئيسي الذي يزعزع استقرار بنية البروتين) والتأثيرات الكارهة للماء (القوى الدافعة الرئيسية لطي البروتين). عامل مهم آخر هو روابط الهيدروجين داخل الجزيئية،  والتي لم يتم تضمينها بشكل صريح في مجالات القوة الحديثة، ولكنها توصف بتفاعلات كولوم لشحنات النقطة الذرية. هذا تقدير تقريبي لأن الروابط الهيدروجينية لها طبيعة ميكانيكية وكيميائية جزئية. علاوة على ذلك، تُحسب التفاعلات الكهروستاتيكية عادةً باستخدام ثابت العزل الكهربائي للفراغ، على الرغم من أن المحلول المائي المحيط به ثابت عازل أعلى بكثير. باستخدام العيانية ثابت العزل الكهربائي على مسافات قصيرة بين الذرات غير مشكوك فيه. أخيرًا، عادةً ما يتم وصف تفاعلات قوى فان دير فالس في الديناميكا الجزيئية بواسطة إمكانات جهد لينارد-جونز استنادًا إلى نظرية فريتز لندن التي لا تنطبق إلا في الفراغ. ومع ذلك، فإن جميع أنواع قوى فان دير فال هي في نهاية المطاف من أصل إلكتروستاتيكي وبالتالي تعتمد على الخصائص العازلة للبيئة. يوضح القياس المباشر لقوى الجذب بين المواد المختلفة (مثل ثابت هامكر) أن «التفاعل بين الهيدروكربونات عبر الماء يبلغ حوالي 10٪ من ذلك عبر الفراغ». تم إهمال الاعتماد على البيئة لقوى فان دير فال في عمليات المحاكاة القياسية، ولكن يمكن تضمينها من خلال تطوير مجالات قوة قابلة للاستقطاب.

قيود التصميم 
يجب أن يأخذ تصميم محاكاة الديناميكيات الجزيئية في الاعتبار القوة الحسابية المتاحة. يجب تحديد حجم المحاكاة (n = عدد الجسيمات)، والوقت، والمدة الزمنية الإجمالية بحيث يمكن إنهاء الحساب في غضون فترة زمنية معقولة. ومع ذلك، يجب أن تكون عمليات المحاكاة طويلة بما يكفي لتكون ذات صلة بالمقاييس الزمنية للعمليات الطبيعية قيد الدراسة. للحصول على استنتاجات صحيحة إحصائيًا من عمليات المحاكاة، يجب أن يتطابق المدى الزمني الذي تمت محاكاته مع حركية العملية الطبيعية. خلافًا لذلك، يكون الأمر مشابهًا للتوصل إلى استنتاجات حول كيفية سير الإنسان عندما ينظر فقط إلى أقل من خطوة واحدة. تستخدم معظم المنشورات العلمية حول ديناميكيات البروتينات والحمض النووي  بيانات من عمليات المحاكاة التي تمتد نانوثانية (10 −9 ثوانٍ) إلى ميكروثانية (10 6 ثوانٍ). للحصول على هذه المحاكاة، هناك حاجة إلى عدة أيام من وحدة المعالجة المركزية إلى سنوات وحدة المعالجة المركزية. تسمح الخوارزميات المتوازية بتوزيع الحمل بين وحدات المعالجة المركزية؛ مثال على ذلك هو خوارزمية التحلل المكاني أو القوة.

أثناء محاكاة الديناميكا الجزيئية الكلاسيكية، فإن المهمة الأكثر كثافة لوحدة المعالجة المركزية هي تقييم الإمكانات كدالة للإحداثيات الداخلية للجسيمات. ضمن تقييم الطاقة هذا، يكون الجزء الأغلى هو الجزء غير المترابط أو غير التساهمي. في الكبير O التدوين، وديناميات الجزيئية المشتركة المحاكاة على نطاق و من قبل  إذا كان كل الزوج الحكيم كهرباء وفان دير فال التفاعلات يجب أن تمثل بشكل واضح. يمكن تقليل هذه التكلفة الحسابية عن طريق استخدام طرق الكهرباء الساكنة مثل شبكة الجسيمات بتلخيص ايوالد () أو الجسيمات - الجسيمات - الجسيمات - الشبكة (P3M)، أو طرق القطع الكروية الجيدة ().

هناك عامل آخر يؤثر على إجمالي وقت وحدة المعالجة المركزية الذي تحتاجه المحاكاة وهو حجم الخطوة الزمنية للتكامل. هذه هي المدة الزمنية بين تقييمات الإمكانات. يجب اختيار الخطوة الزمنية صغيرة بما يكفي لتجنب أخطاء التقدير (أي أصغر من الفترة المتعلقة بأسرع تردد اهتزازي في النظام). الخطوات الزمنية النموذجية لـ الديناميكا الجزيئية الكلاسيكي بترتيب 1 فيمتوثانية (10 - 15 ثانية). يمكن توسيع هذه القيمة باستخدام الخوارزميات مثل خوارزمية قيود شيك، والتي تثبت اهتزازات الذرات الأسرع (على سبيل المثال، الهيدروجين) في مكانها. تم أيضًا تطوير طرق متعددة للمقياس الزمني، والتي تتيح فترات ممتدة بين تحديثات القوى البعيدة المدى الأبطأ.

لمحاكاة الجزيئات في مذيب، يجب الاختيار بين المذيب الصريح والضمني. يجب حساب جسيمات المذيبات الصريحة (مثل  و SPC / E وSPC-f) بشكل مكلف عن طريق مجال القوة، بينما تستخدم المذيبات الضمنية نهج المجال المتوسط. يعد استخدام مذيب صريح مكلفًا من الناحية الحسابية، ويتطلب تضمين ما يقرب من عشرة أضعاف الجسيمات في المحاكاة. لكن حبيبات ولزوجة المذيب الصريح ضروريان لإعادة إنتاج خصائص معينة للجزيئات الذائبة. هذا مهم بشكل خاص لإعادة إنتاج الخواص الحركية الكيميائية.

في جميع أنواع محاكاة الديناميكيات الجزيئية، يجب أن يكون حجم صندوق المحاكاة كبيرًا بما يكفي لتجنب آثار حالة الحدود. غالبًا ما يتم التعامل مع شروط الحدود عن طريق اختيار القيم الثابتة عند الحواف (والتي قد تسبب تشوهات)، أو عن طريق استخدام شروط الحدود الدورية حيث يعود جانب واحد من حلقات المحاكاة إلى الجانب الآخر، ومحاكاة مرحلة مجمعة (والتي قد تسبب آثارًا أيضًا).

مجموعة متناهية الصغر (NVE) 
في المجموعة المتناهية الصغر، يتم عزل النظام عن التغيرات في الشامات (N) والحجم (V) والطاقة (E). يتوافق مع عملية ثابتة الحرارة بدون تبادل حراري. يمكن النظر إلى مسار الديناميات الجزيئية المتناهية الصغر على أنه تبادل للطاقة الكامنة والحركية، مع الحفاظ على إجمالي الطاقة. لنظام من الجسيمات N مع الإحداثيات  والسرعات ، يمكن كتابة الزوج التالي من المعادلات التفاضلية من الدرجة الأولى في تدوين نيوتن كـ

 
 

وظيفة الطاقة الكامنة  النظام هو دالة إحداثيات الجسيمات  . يشار إليه ببساطة على أنه الإمكانات في الفيزياء، أو مجال القوة في الكيمياء. تأتي المعادلة الأولى من قوانين نيوتن للحركة. القوة  يمكن حساب التأثير على كل جسيم في النظام على أنه التدرج السالب لـ  .

لكل خطوة زمنية، موضع كل جسيم  والسرعة  قد تتكامل مع طريقة تكامل سيمبلكتيك مثل تكامل فيرلت. التطور الزمني ل  و  يسمى المسار. بالنظر إلى المواقف الأولية (على سبيل المثال، من المعرفة النظرية) والسرعات (على سبيل المثال، غاوسي العشوائي)، يمكننا حساب جميع المواقف والسرعات المستقبلية (أو السابقة).

أحد مصادر الالتباس المتكررة هو معنى درجة الحرارة في الديناميكا الجزيئية. عادة لدينا خبرة في درجات الحرارة العيانية، والتي تنطوي على عدد كبير من الجسيمات. لكن درجة الحرارة هي كمية إحصائية. إذا كان هناك عدد كبير بما يكفي من الذرات، فيمكن تقدير درجة الحرارة الإحصائية من درجة الحرارة اللحظية، والتي يتم العثور عليها من خلال مساواة الطاقة الحركية للنظام بـ nk B T / 2 حيث n هو عدد درجات الحرية للنظام.

تنشأ ظاهرة مرتبطة بدرجة الحرارة بسبب قلة عدد الذرات المستخدمة في محاكاة الديناميكا الجزيئية. على سبيل المثال، ضع في اعتبارك محاكاة نمو فيلم نحاسي بدءًا من ركيزة تحتوي على 500 ذرة وطاقة ترسيب تبلغ 100 فولت. في العالم الحقيقي، سيتم نقل 100 فولت من الذرة المودعة بسرعة ومشاركتها بين عدد كبير من الذرات ( أو أكثر) دون تغيير كبير في درجة الحرارة. عندما يكون هناك 500 ذرة فقط، فإن الركيزة تتبخر على الفور تقريبًا بواسطة الترسب. يحدث شيء مشابه في عمليات المحاكاة الفيزيائية الحيوية. يتم رفع درجة حرارة النظام في (NVE) بشكل طبيعي عندما تخضع الجزيئات الكبيرة مثل البروتينات لتغيرات وترابط طارد للحرارة.

المجموعة المتعارف عليها (NVT) 
في المجموعة المتعارف عليها، يتم حفظ كمية المادة (N) والحجم (V) ودرجة الحرارة (T). كما يطلق عليه أحيانًا الديناميكيات الجزيئية ذات درجة الحرارة الثابتة (CTالديناميكا الجزيئية). في (NVT)، يتم تبادل طاقة العمليات الماصة للحرارة والطاردة للحرارة مع منظم الحرارة.

تتوفر مجموعة متنوعة من خوارزميات الترموستات لإضافة الطاقة وإزالتها من حدود محاكاة الديناميكا الجزيئية بطريقة أكثر أو أقل واقعية، بالقرب من المجموعة المتعارف عليها. تشمل الأساليب الشائعة للتحكم في درجة الحرارة إعادة قياس السرعة، ترموستات نوزيه هوفر، وسلاسل نوزيه هوفر، وترموستات بيرندسن، وديناميكيات أندرسن ولانجفين. قد يقدم منظم الحرارة بيرندسن تأثير مكعب الثلج الطائر، والذي يؤدي إلى ترجمات غير مادية وتناوب للنظام المحاكي.

ليس من التافه الحصول على توزيع متعارف عليه للمطابقة والسرعات باستخدام هذه الخوارزميات. كيف يعتمد ذلك على حجم النظام واختيار منظم الحرارة ومعلمات منظم الحرارة وخطوة الوقت والمتكامل هو موضوع العديد من المقالات في هذا المجال.

مجموعة متساوية الحرارة متساوية الضغط (NPT) 
يتم الحفاظ على كمية المادة (N) والضغط (P) ودرجة الحرارة (T) في المجموعة متساوي الحرارة - متساوي الضغط. بالإضافة إلى منظم الحرارة، هناك حاجة إلى باروستات. يتوافق بشكل وثيق مع ظروف المختبر مع وجود دورق مفتوح لدرجة الحرارة والضغط المحيطين.

في محاكاة الأغشية البيولوجية، فإن التحكم في الضغط الخواص غير مناسب. بالنسبة لطبقات الدهون الثنائية، يحدث التحكم في الضغط تحت منطقة غشاء ثابتة (NPAT) أو توتر سطحي ثابت «جاما» (NPγT).

المجموعات المعممة 
طريقة تبادل النسخ المتماثلة هي مجموعة معممة. تم إنشاؤه في الأصل للتعامل مع الديناميكيات البطيئة لأنظمة الدوران المضطربة. يطلق عليه أيضًا التقسية الموازية. تحاول صياغة الديناميكا الجزيئية للتبادل المتماثل  التغلب على مشكلة الحدود الدنيا المتعددة من خلال تبادل درجة حرارة النسخ المتماثلة غير المتفاعلة للنظام التي تعمل في درجات حرارة متعددة.

الإمكانات في محاكاة الديناميكا الجزيئية 
تتطلب محاكاة الديناميكيات الجزيئية تعريف وظيفة محتملة، أو وصفًا للمصطلحات التي تتفاعل بها الجسيمات في المحاكاة. في الكيمياء والبيولوجيا، يُشار إلى هذا عادةً باسم مجال القوة وفي فيزياء المواد على أنه جهد بين الذرات. يمكن تحديد الإمكانات على مستويات عديدة من الدقة المادية؛ تلك الأكثر شيوعًا في الكيمياء تعتمد على الميكانيكا الجزيئية وتجسد العلاج الميكانيكي الكلاسيكي لتفاعلات الجسيمات والجسيمات التي يمكن أن تعيد إنتاج التغييرات الهيكلية والتكوينية ولكن لا يمكنها عادةً إعادة إنتاج التفاعلات الكيميائية.

يستلزم الاختزال من الوصف الكمي الكامل إلى الإمكانات الكلاسيكية تقريبيين رئيسيين. الأول هو تقريب بورن أوبنهايمر، الذي ينص على أن ديناميكيات الإلكترونات سريعة جدًا بحيث يمكن اعتبارها تتفاعل على الفور مع حركة نواتها. نتيجة لذلك، يمكن معاملتهم بشكل منفصل. الثاني يتعامل مع النوى، وهي أثقل بكثير من الإلكترونات، على أنها جسيمات نقطية تتبع الديناميكيات النيوتونية الكلاسيكية. في الديناميكيات الجزيئية الكلاسيكية، يتم تقريب تأثير الإلكترونات كسطح طاقة كامن، وعادة ما يمثل الحالة الأرضية.

عندما تكون هناك حاجة إلى مستويات أدق من التفاصيل، يتم استخدام الإمكانات القائمة على ميكانيكا الكم ؛ تحاول بعض الطرق إنشاء إمكانات كلاسيكية / كمومية هجينة حيث يتم معالجة الجزء الأكبر من النظام بشكل كلاسيكي ولكن يتم التعامل مع منطقة صغيرة كنظام كمي، وعادةً ما تخضع لعملية تحول كيميائي.

الإمكانات التجريبية 
كثيرا ما تسمى الإمكانات التجريبية المستخدمة في الكيمياء مجالات القوة، في حين أن تلك المستخدمة في فيزياء المواد تسمى الجهود بين الذرية.

معظم مجالات القوة في الكيمياء تجريبية وتتكون من مجموع القوى المرتبطة المرتبطة بالروابط الكيميائية وزوايا الرابطة وثنائيات السندات والقوى غير المرتبطة المرتبطة بقوى فان دير فال والشحنة الكهروستاتيكية. تمثل الإمكانات التجريبية تأثيرات ميكانيكا الكم بطريقة محدودة من خلال التقريبات الوظيفية المخصصة. تحتوي هذه الإمكانات على معلمات مجانية مثل الشحنة الذرية، ومعلمات فان دير فالس التي تعكس تقديرات نصف القطر الذري، وطول رابطة التوازن، والزاوية، وثنائي السطوح ؛ يتم الحصول عليها عن طريق التوافق مع الحسابات الإلكترونية التفصيلية (المحاكاة الكيميائية الكمومية) أو الخصائص الفيزيائية التجريبية مثل الثوابت المرنة والمعلمات الشبكية والقياسات الطيفية.

بسبب الطبيعة غير المحلية للتفاعلات غير المرتبطة، فإنها تنطوي على تفاعلات ضعيفة على الأقل بين جميع الجسيمات في النظام. عادة ما يكون حسابه هو عنق الزجاجة في سرعة محاكاة الديناميكا الجزيئية. لخفض التكلفة الحسابية، تستخدم حقول القوة تقديرات تقريبية رقمية مثل نصف قطر القطع المغير، أو خوارزميات مجال التفاعل، أو تجميع ايوالد لشبكة الجسيمات، أو شبكة الجسيمات والجسيمات والجسيمات الأحدث (P3M).

تستخدم حقول القوة الكيميائية عادةً ترتيبات ارتباط محددة مسبقًا (استثناء هو ديناميكيات طرق كيمياء الكم)، وبالتالي فهي غير قادرة على نمذجة عملية كسر الرابطة الكيميائية والتفاعلات بشكل صريح. من ناحية أخرى، يمكن للعديد من الإمكانات المستخدمة في الفيزياء، مثل تلك القائمة على شكلية ترتيب السندات، أن تصف عدة تنسيقات مختلفة للنظام وكسر الرابطة. تتضمن أمثلة هذه الإمكانات إمكانات برينر  للهيدروكربونات وتطوراتها الإضافية لأنظمة C-Si-H  و COH. يمكن اعتبار إمكانات  هجينًا كامل التفاعل بين إمكانات ترتيب الرابطة ومجالات القوة الكيميائية.

قم بإقران الإمكانات مقابل إمكانات الأجسام المتعددة 
تتم صياغة الوظائف المحتملة التي تمثل الطاقة غير المرتبطة كمجموع على التفاعلات بين جسيمات النظام. أبسط خيار، مستخدم في العديد من مجالات القوة الشائعة، هو «جهد الزوج»، حيث يمكن حساب إجمالي الطاقة الكامنة من مجموع مساهمات الطاقة بين أزواج الذرات. لذلك، تسمى حقول القوة هذه أيضًا «حقول القوة المضافة». مثال على مثل هذا الزوج المحتمل هو احتمال جهد لينارد-جونز غير المرتبط (يُطلق عليه أيضًا احتمال 6-12)، المستخدم لحساب قوى فان دير فال.

 

مثال آخر هو نموذج Born (الأيوني) للشبكة الأيونية. المصطلح الأول في المعادلة التالية هو قانون كولوم لزوج من الأيونات، والمصطلح الثاني هو التنافر قصير المدى الذي شرحه مبدأ استبعاد باولي والمصطلح الأخير هو مصطلح تفاعل التشتت. عادةً ما تتضمن المحاكاة المصطلح ثنائي القطب فقط، على الرغم من تضمين المصطلح رباعي الأقطاب في بعض الأحيان. عندما ن ل = 6، وتسمى هذه الإمكانية أيضا إمكانية كولوم-باكنغهام.

 

في إمكانات الجسم المتعددة، تشتمل الطاقة الكامنة على تأثيرات ثلاثة أو أكثر من الجسيمات التي تتفاعل مع بعضها البعض. في عمليات المحاكاة ذات الإمكانات الزوجية، توجد أيضًا تفاعلات عالمية في النظام، ولكنها تحدث فقط من خلال المصطلحات الزوجية. في إمكانات الأجسام المتعددة، لا يمكن العثور على الطاقة الكامنة بمجموع على أزواج من الذرات، حيث يتم حساب هذه التفاعلات بشكل صريح على أنها مزيج من المصطلحات ذات الترتيب الأعلى. في العرض الإحصائي، لا يمكن التعبير عن التبعية بين المتغيرات بشكل عام باستخدام المنتجات المزدوجة فقط من درجات الحرية. على سبيل المثال، إمكانات تيرسوف،  التي كانت تُستخدم في الأصل لمحاكاة الكربون، والسيليكون، والجرمانيوم، ومنذ ذلك الحين تم استخدامها لمجموعة واسعة من المواد الأخرى، تتضمن مجموع مجموعات من ثلاث ذرات، مع وجود زوايا بين الذرات عامل مهم في الإمكانات. ومن الأمثلة الأخرى طريقة الذرة المدمجة (EAM)،  (EDIP)،  وإمكانات تقريب اللحظات الثانية المحكم الربط (TBSMA)،  حيث تكون كثافة الإلكترون للحالات في منطقة الذرة هي محسوبة من مجموع المساهمات من الذرات المحيطة، وبالتالي فإن مساهمة الطاقة الكامنة هي دالة لهذا المجموع.

الإمكانات شبه التجريبية 
تستفيد الإمكانات شبه التجريبية من تمثيل المصفوفة من ميكانيكا الكم. ومع ذلك، يتم العثور على قيم عناصر المصفوفة من خلال الصيغ التجريبية التي تقدر درجة التداخل في المدارات الذرية المحددة. يتم بعد ذلك مائل المصفوفة لتحديد شغل المدارات الذرية المختلفة، ويتم استخدام الصيغ التجريبية مرة أخرى لتحديد مساهمات الطاقة في المدارات.

هناك مجموعة متنوعة من الإمكانات شبه التجريبية، تسمى إمكانات الربط المحكم، والتي تختلف وفقًا للذرات التي يتم نمذجتها.

إمكانات الاستقطاب 
تتضمن معظم مجالات القوة الكلاسيكية ضمنيًا تأثير الاستقطاب، على سبيل المثال، عن طريق زيادة الشحنات الجزئية التي تم الحصول عليها من الحسابات الكيميائية الكمومية. هذه الشحنات الجزئية ثابتة بالنسبة إلى كتلة الذرة. ولكن يمكن لمحاكاة الديناميات الجزيئية أن تكون نموذجًا صريحًا لقابلية الاستقطاب مع إدخال ثنائيات الأقطاب المستحثة من خلال طرق مختلفة، مثل الجسيمات الدقيقة أو الشحنات المتقلبة. هذا يسمح بإعادة التوزيع الديناميكي للشحنة بين الذرات والتي تستجيب للبيئة الكيميائية المحلية.

لسنوات عديدة، تم وصف عمليات محاكاة الديناميكا الجزيئية القابلة للاستقطاب على أنها الجيل التالي. بالنسبة للسوائل المتجانسة مثل الماء، تم تحقيق دقة متزايدة من خلال إدراج قابلية الاستقطاب. كما تم تحقيق بعض النتائج الواعدة للبروتينات. ومع ذلك، لا يزال من غير المؤكد أفضل طريقة لتقريب الاستقطاب في المحاكاة.

الإمكانات في طرق كيمياء الكم 
في الديناميكيات الجزيئية الكلاسيكية، يتم تمثيل سطح طاقة كامنة (عادةً الحالة الأرضية) في مجال القوة. هذا نتيجة لتقريب بورن أوبنهايمر. في الحالات المثارة أو التفاعلات الكيميائية أو عند الحاجة إلى تمثيل أكثر دقة، يمكن الحصول على السلوك الإلكتروني من المبادئ الأولى باستخدام طريقة ميكانيكا الكم، مثل نظرية الكثافة الوظيفية. يسمى هذا الديناميكيات الجزيئية لكيمياء الكم. نظرًا لتكلفة معالجة درجات الحرية الإلكترونية، فإن العبء الحسابي لهذه المحاكاة أعلى بكثير من الديناميكيات الجزيئية التقليدية. لهذا السبب، عادةً ما يقتصر نظام الديناميكا الجزيئية على أنظمة أصغر وأوقات أقصر.

يمكن استخدام الطرق الميكانيكية والكيميائية الكمومية للمبتدئين لحساب الطاقة الكامنة لنظام أثناء الطيران، حسب الحاجة للتوافق في المسار. يتم إجراء هذا الحساب عادةً في الجوار القريب لإحداثيات التفاعل. على الرغم من أنه يمكن استخدام تقديرات تقريبية مختلفة، إلا أنها تستند إلى اعتبارات نظرية، وليس على ملاءمة تجريبية. تنتج حسابات كيمياء الكم قدرًا هائلاً من المعلومات غير المتوفرة من الأساليب التجريبية، مثل كثافة الحالات الإلكترونية أو الخصائص الإلكترونية الأخرى. من المزايا المهمة لاستخدام طرق كيمياء الكم القدرة على دراسة التفاعلات التي تتضمن كسر أو تكوين روابط تساهمية، والتي تتوافق مع حالات إلكترونية متعددة. علاوة على ذلك، تسمح طرق كيمياء الكم أيضًا باستعادة التأثيرات التي تتجاوز تقريب بورن-أوبنهايمر باستخدام مناهج مثل الديناميكيات الكمومية الكلاسيكية المختلطة.

هجين QM / MM 
تعتبر طرق QM (ميكانيكا الكم) قوية جدًا. ومع ذلك، فهي مكلفة من الناحية الحسابية، في حين أن طرق MM (الميكانيكا الكلاسيكية أو الجزيئية) سريعة ولكنها تعاني من عدة حدود (تتطلب معايير واسعة النطاق ؛ تقديرات الطاقة التي تم الحصول عليها ليست دقيقة للغاية ؛ لا يمكن استخدامها لمحاكاة التفاعلات حيث يتم كسر / تشكيل الروابط التساهمية ؛ ومحدودي قدراتهم على تقديم تفاصيل دقيقة فيما يتعلق بالبيئة الكيميائية). ظهرت فئة جديدة من الأساليب التي تجمع بين النقاط الجيدة في حسابات QM (الدقة) و MM (السرعة). تسمى هذه الطرق طرق ميكانيكا الكم الميكانيكية والجزيئية المختلطة أو المختلطة (الهجين QM / MM).

إن أهم ميزة لطريقة إدارة الجودة (QM / MM) الهجينة هي السرعة. تكلفة عمل الديناميكيات الجزيئية الكلاسيكية (MM) في مقاييس الحالة الأكثر وضوحًا O (n 2)، حيث n هو عدد الذرات في النظام. هذا يرجع أساسًا إلى مصطلح التفاعلات الكهروستاتيكية (يتفاعل كل جسيم مع كل جسيم آخر). ومع ذلك، فإن استخدام نصف قطر القطع والتحديثات الدورية لقائمة الأزواج ومؤخرًا الاختلافات في طريقة شبكة الجسيمات ايوالد (PME) قد قلل هذا إلى ما بين O (n) إلى O (n 2). بمعنى آخر، إذا تمت محاكاة نظام به ضعف عدد الذرات، فسيستغرق ما بين ضعفين إلى أربعة أضعاف قوة الحوسبة. من ناحية أخرى، فإن أبسط حسابات كيمياء الكم مقياس نموذجي O (n 3) أو أسوأ (تم اقتراح حسابات هارتري فوك المقيدة لتوسيع نطاق ~ O (n 2.7)). للتغلب على هذا الحد، يتم التعامل مع جزء صغير من النظام ميكانيكيًا (عادةً موقع نشط من الإنزيم) ويتم معالجة النظام المتبقي بشكل كلاسيكي.

في التطبيقات الأكثر تعقيدًا، توجد طرق QM / MM لمعالجة كل من النوى الخفيفة المعرضة للتأثيرات الكمومية (مثل الهيدروجين) والحالات الإلكترونية. هذا يسمح بتوليد وظائف موجات الهيدروجين (على غرار وظائف الموجة الإلكترونية). كانت هذه المنهجية مفيدة في التحقيق في ظواهر مثل نفق الهيدروجين. أحد الأمثلة حيث قدمت طرق QM / MM اكتشافات جديدة هو حساب نقل الهيدريد في إنزيم نازع هيدروجين الكحول في الكبد. في هذه الحالة، نفق الكم مهم للهيدروجين، لأنه يحدد معدل التفاعل.

الخشنة الحبيبات والتمثيلات المخفضة 
في الطرف الآخر من مقياس التفاصيل توجد نماذج ذات حبيبات خشنة ونماذج شعرية. بدلاً من التمثيل الصريح لكل ذرة في النظام، يستخدم المرء «ذرات زائفة» لتمثيل مجموعات من الذرات. قد تتطلب عمليات محاكاة الديناميكا الجزيئية على أنظمة كبيرة جدًا موارد كمبيوتر كبيرة بحيث لا يمكن دراستها بسهولة من خلال طرق كل الذرة التقليدية. وبالمثل، فإن محاكاة العمليات على نطاقات زمنية طويلة (تتجاوز حوالي 1 ميكرو ثانية) باهظة التكلفة، لأنها تتطلب الكثير من الخطوات الزمنية. في هذه الحالات، يمكن للمرء أحيانًا معالجة المشكلة باستخدام تمثيلات مختصرة، والتي تسمى أيضًا نماذج الحبيبات الخشنة.

من الأمثلة على طرق الحبيبات الخشنة (CG) الديناميكيات الجزيئية المتقطعة (CG-Dالديناميكا الجزيئية)  ونماذج Go. يتم إجراء الحبيبات الخشنة أحيانًا بأخذ ذرات زائفة أكبر. تم استخدام مثل هذه التقريبات الذرية الموحدة في محاكاة الديناميكا الجزيئية للأغشية البيولوجية. قد يكون تنفيذ مثل هذا النهج على الأنظمة التي تكون فيها الخصائص الكهربائية ذات أهمية أمرًا صعبًا نظرًا لصعوبة استخدام توزيع الشحنة المناسب على الذرات الزائفة. يتم تمثيل ذيول الليبيدات الأليفاتية بواسطة عدد قليل من الذرات الزائفة عن طريق تجميع 2 إلى 4 مجموعات ميثيلين في كل ذرة زائفة.

يجب أن يتم تحديد معلمات هذه النماذج الخشنة للغاية بشكل تجريبي، من خلال مطابقة سلوك النموذج مع البيانات التجريبية المناسبة أو محاكاة جميع الذرة. من الناحية المثالية، يجب أن تأخذ هذه المعلمات في الاعتبار المساهمات الحرارية والنتروبية في الطاقة الحرة بطريقة ضمنية. عندما يتم إجراء الحبيبات الخشنة على مستويات أعلى، قد تكون دقة الوصف الديناميكي أقل موثوقية. ولكن تم استخدام النماذج ذات الحبيبات الخشنة بنجاح لفحص مجموعة واسعة من الأسئلة في البيولوجيا الهيكلية وتنظيم البلورات السائلة وزجاج البوليمر.

أمثلة على تطبيقات الحبيبات الخشنة:

 غالبًا ما تُجرى دراسات طي البروتين ودراسات التنبؤ ببنية البروتين باستخدام ذرة زائفة واحدة أو بضع ذرات لكل حمض أميني.
 تم فحص انتقالات الطور البلوري السائل في الأشكال الهندسية المحصورة و / أو أثناء التدفق باستخدام إمكانات جاي بيرن، التي تصف الأنواع متباينة الخواص.
 تمت دراسة زجاج البوليمر أثناء التشوه باستخدام نوابض توافقية بسيطة أو نوابض فين لتوصيل الكرات الموصوفة بواسطة إمكانات لينارد جونز.
 تم فحص الالتفاف الفائق للحمض النووي باستخدام 1-3 ذرات زائفة لكل زوج أساس، وبدقة أقل.
 تم فحص تغليف الحمض النووي المزدوج الحلزوني في العاثية بنماذج حيث تمثل الذرة الزائفة منعطفًا واحدًا (حوالي 10 أزواج أساسية) من الحلزون المزدوج.
 تم تصميم بنية الحمض النووي الريبي في الريبوسوم والأنظمة الكبيرة الأخرى باستخدام ذرة زائفة واحدة لكل نوكليوتيد.
 محاكاة افتراضية للخلية لدراسة تفاعل الخلايا والركائز المختلفة.

أبسط شكل من أشكال الحبيبات الخشنة هو الذرة المتحدة (تسمى أحيانًا الذرة الممتدة)، وقد تم استخدامها في معظم عمليات محاكاة الديناميكا الجزيئية المبكرة للبروتينات والدهون والأحماض النووية. على سبيل المثال، بدلاً من معالجة جميع الذرات الأربع لمجموعة ميثيل CH 3 بشكل صريح (أو جميع الذرات الثلاث لمجموعة الميثيلين CH 2)، يمثل المرء المجموعة بأكملها بذرة زائفة واحدة. يجب بالطبع تحديد معلماتها بشكل صحيح بحيث تتمتع تفاعلاتها مع المجموعات الأخرى بالاعتماد المناسب على المسافة. تنطبق اعتبارات مماثلة على الروابط والزوايا والالتواءات التي تشارك فيها الذرة الزائفة. في هذا النوع من التمثيل الذري الموحد، يزيل المرء عادةً جميع ذرات الهيدروجين الصريحة باستثناء تلك التي لديها القدرة على المشاركة في الروابط الهيدروجينية (الهيدروجين القطبي). مثال على ذلك هو مجال القوة تشارم 19.

عادةً ما يتم الاحتفاظ بالهيدروجين القطبي في النموذج، لأن المعالجة المناسبة للروابط الهيدروجينية تتطلب وصفًا دقيقًا بشكل معقول للاتجاه والتفاعلات الكهروستاتيكية بين مجموعات المتبرع والمستقبل. يمكن لمجموعة الهيدروكسيل، على سبيل المثال، أن تكون مانحًا لرابطة الهيدروجين، ومتقبلًا لرابطة الهيدروجين، وسيكون من المستحيل معالجة ذلك باستخدام ذرة زائفة (OH) واحدة. حوالي نصف الذرات الموجودة في البروتين أو الحمض النووي هي هيدروجين غير قطبي، لذا فإن استخدام الذرات الموحدة يمكن أن يوفر توفيرًا كبيرًا في وقت الكمبيوتر.

دمج تأثيرات المذيبات 
في العديد من عمليات محاكاة نظام المذيب المذاب، ينصب التركيز الرئيسي على سلوك المذاب مع القليل من الاهتمام بسلوك المذيب خاصة في جزيئات المذيب الموجودة في مناطق بعيدة عن الجزيء المذاب. قد تؤثر المذيبات على السلوك الديناميكي للمذابات من خلال الاصطدامات العشوائية وفرض مقاومة احتكاكية على حركة المذاب عبر المذيب. يمكن أن يساعد استخدام شروط الحدود الدورية غير المستطيلة والحدود العشوائية وأغلفة المذيبات في تقليل عدد جزيئات المذيبات المطلوبة وتمكين نسبة أكبر من وقت الحوسبة ليتم إنفاقها بدلاً من ذلك في محاكاة المذاب. من الممكن أيضًا دمج تأثيرات المذيب دون الحاجة إلى وجود أي جزيئات مذيب واضحة. أحد الأمثلة على هذا النهج هو استخدام القوة المتوسطة المحتملة (PMF) التي تصف كيفية تغير الطاقة الحرة كتنوع إحداثيات معينة. يحتوي تغيير الطاقة الحرة الموصوف بواسطة (PMF) على متوسط تأثيرات المذيب.

قوى بعيدة المدى 
التفاعل طويل المدى هو تفاعل لا يسقط فيه التفاعل المكاني أسرع من  أين  هي أبعاد النظام. تشمل الأمثلة تفاعلات الشحنة والشحنة بين الأيونات وتفاعلات ثنائي القطب بين الجزيئات. تمثل نمذجة هذه القوى تحديًا كبيرًا لأنها كبيرة على مسافة قد تكون أكبر من نصف طول الصندوق مع محاكاة لعدة آلاف من الجسيمات. على الرغم من أن أحد الحلول سيكون زيادة حجم طول الصندوق بشكل كبير، إلا أن أسلوب القوة الغاشمة هذا أقل من مثالي لأن المحاكاة ستصبح باهظة الثمن من الناحية الحسابية. إن اقتطاع الإمكانات كرويًا هو أيضًا غير وارد حيث يمكن ملاحظة السلوك غير الواقعي عندما تكون المسافة قريبة من مسافة القطع.

الديناميكات الجزيئية الموجهة (SMD) 
تقوم محاكاة الديناميكيات الجزيئية (Sالديناميكا الجزيئية)، أو محاكاة مسبار القوة، بتطبيق قوى على بروتين من أجل التلاعب بهيكله عن طريق سحبه على طول درجات الحرية المرغوبة. يمكن استخدام هذه التجارب للكشف عن التغيرات الهيكلية في البروتين على المستوى الذري. غالبًا ما يستخدم Sالديناميكا الجزيئية لمحاكاة أحداث مثل التكشف أو التمدد الميكانيكي.

هناك نوعان من البروتوكولات النموذجية لـ الديناميكا الجزيئية: أحدهما تكون فيه سرعة السحب ثابتة، والآخر تكون فيه القوة المطبقة ثابتة. بشكل نموذجي، يتم تقييد جزء من النظام المدروس (على سبيل المثال، ذرة في بروتين) بواسطة جهد توافقي. ثم يتم تطبيق القوى على ذرات معينة إما بسرعة ثابتة أو بقوة ثابتة. يتم استخدام أخذ العينات الشامل لتحريك النظام على طول تنسيق التفاعل المطلوب عن طريق تغيير، على سبيل المثال، القوى والمسافات والزوايا التي تم التلاعب بها في المحاكاة. من خلال أخذ العينات المظلة، يتم أخذ عينات من جميع تكوينات النظام - عالية الطاقة ومنخفضة الطاقة - بشكل كافٍ. بعد ذلك، يمكن حساب تغيير كل تكوين في الطاقة الحرة على أنه احتمال متوسط القوة. إحدى الطرق الشائعة لحساب (PMF) هي من خلال طريقة تحليل الرسم البياني الموزون (WHAM)، والتي تحلل سلسلة من محاكاة أخذ العينات الشاملة.

يوجد الكثير من التطبيقات المهمة لـ Sالديناميكا الجزيئية في مجال اكتشاف الأدوية وعلوم الجزيئات الحيوية. على سبيل المثال تم استخدام Sالديناميكا الجزيئية للتحقيق في استقرار ليفية الزهايمر الأولية،  لدراسة تفاعل بروتين يجند في كيناز المعتمد على السيكلين 5  وحتى لإظهار تأثير المجال الكهربائي على مركب الثرومبين (البروتين) والأبتامر (النيوكليوتيد)  بين العديد من الدراسات الأخرى المثيرة للاهتمام.

أمثلة على التطبيقات 

تستخدم الديناميات الجزيئية في العديد من مجالات العلوم.

 تم نشر أول محاكاة الديناميكا الجزيئية لعملية الطي البيولوجي المبسطة في عام 1975. مهدت المحاكاة التي نُشرت في دورية نتشر الطريق لمساحة واسعة من طي البروتين الحسابي الحديث.
 تم نشر أول محاكاة الديناميكا الجزيئية لعملية بيولوجية في عام 1976. مهدت المحاكاة التي نُشرت في دورية نتشر الطريق لفهم حركة البروتين باعتبارها ضرورية في الوظيفة وليس مجرد ملحق.
 الديناميكا الجزيئية هي الطريقة القياسية لمعالجة شلالات الاصطدام في نظام ارتفاع الحرارة، أي التأثيرات التي تحدثها إشعاعات النيوترونات والأيونات النشطة على المواد الصلبة والأسطح الصلبة.

توضح الأمثلة الفيزيائية الحيوية التالية الجهود الملحوظة لإنتاج محاكاة لأنظمة ذات حجم كبير جدًا (فيروس كامل) أو أوقات محاكاة طويلة جدًا (تصل إلى 1.112 مللي ثانية):

 محاكاة الديناميكا الجزيئية لفيروس فسيفساء التبغ الساتلي الكامل (STMV) (2006، الحجم: 1 مليون ذرة، وقت المحاكاة: 50 نانوثانية، البرنامج: NAالديناميكا الجزيئية) هذا الفيروس عبارة عن فيروس نباتي صغير عشروني الوجوه يؤدي إلى تفاقم أعراض الإصابة بفيروس التبغ الفسيفسائي (TMV). تم استخدام محاكاة الديناميكيات الجزيئية لاستكشاف آليات التجمع الفيروسي. يتكون جسيم STMV بأكمله من 60 نسخة متطابقة من بروتين واحد يتكون من القفيصة الفيروسية (الغلاف)، وجينوم RNA أحادي النوكليوتيد 1063 تقطعت به السبل. تتمثل إحدى النتائج الرئيسية في أن الكابسيد يكون غير مستقر للغاية عند عدم وجود RNA بداخله. ستستغرق المحاكاة جهاز كمبيوتر مكتبي واحدًا في عام 2006 حوالي 35 عامًا حتى تكتمل. وهكذا تم القيام به في العديد من المعالجات بالتوازي مع التواصل المستمر بينهم.
 محاكاة قابلة للطي لـ فيلين هيدبيس بتفاصيل شاملة (2006، الحجم: 20000 ذرة ؛ وقت المحاكاة: 500 ميكرو ثانية = 500000 نانوثانية، البرنامج: فولدنغ@هوم) تم تشغيل هذه المحاكاة في 200000 وحدة معالجة مركزية من أجهزة الكمبيوتر الشخصية المشاركة حول العالم. تم تثبيت برنامج فولدنغ@هوم على أجهزة الكمبيوتر هذه، وهو جهد حوسبي موزع على نطاق واسع بتنسيق من فيجاي باندي في جامعة ستانفورد. تم فحص الخصائص الحركية لبروتين فيلين هيدبيس باستخدام العديد من المسارات المستقلة والقصيرة التي تديرها وحدة المعالجة المركزية دون اتصال مستمر في الوقت الحقيقي. كانت إحدى الطرق المستخدمة هي تحليل قيمة بفولد، والتي تقيس احتمالية الطي قبل الكشف عن تشكيل بداية محدد. يعطي بفولد معلومات حول هياكل الحالة الانتقالية وترتيب المطابقة على طول مسار الطي. يمكن أن يكون كل مسار في حساب بفولد قصيرًا نسبيًا، ولكن هناك حاجة إلى العديد من المسارات المستقلة.
 تم إجراء عمليات محاكاة طويلة للمسار المستمر على انطون، وهو كمبيوتر عملاق متوازي بشكل كبير تم تصميمه وبنائه حول دوائر متكاملة خاصة بالتطبيقات المخصصة (ASIC) ووصلات متصلة بواسطة دي شو للأبحاث. أطول نتيجة تم نشرها لمحاكاة تم إجراؤها باستخدام أنطون هي محاكاة NTL9 تبلغ 1.112 مللي ثانية عند 355 كلفن ؛ تم أيضًا إجراء محاكاة ثانية مستقلة تبلغ 1.073 مللي ثانية لهذا التكوين (والعديد من عمليات المحاكاة الأخرى لأكثر من 250 ميكرو ثانية من الوقت الكيميائي المستمر). في كيف تطوى البروتينات سريعة الطي، يناقش الباحثون كريستين ليندورف لارسن وستيفانو بيانا ورون أو درور وديفيد إي. تكشف عن مجموعة من المبادئ المشتركة الكامنة وراء طي 12 بروتينًا متنوعًا هيكليًا ". يسمح فحص هذه المسارات الطويلة المتنوعة، التي تم تمكينها بواسطة أجهزة مخصصة ومخصصة، باستنتاج أنه "في معظم الحالات، يتبع الطي مسارًا واحدًا مهيمنًا تظهر فيه عناصر الهيكل الأصلي بترتيب مرتبط بشدة بميلها للتشكل في تكشفت الدولة ". في دراسة منفصلة، تم استخدام أنطون لإجراء محاكاة 1.013 مللي ثانية لديناميات الحالة الأصلية لمثبط التربسين البنكرياس البقري (BPTI) عند 300 كلفن.

يستفيد تطبيق مهم آخر لطريقة الديناميكا الجزيئية من قدرتها على التوصيف ثلاثي الأبعاد وتحليل التطور المجهرية على النطاق الذري.

 تُستخدم محاكاة الديناميكا الجزيئية في توصيف تطور حجم الحبوب، على سبيل المثال، عند وصف التآكل والاحتكاك لمواد Al و Al (Zr) النانوية. يتم تحليل تطور الخلع وتطور حجم الحبيبات أثناء عملية الاحتكاك في هذه المحاكاة. نظرًا لأن طريقة الديناميكا الجزيئية قدمت المعلومات الكاملة للبنية المجهرية، فقد تم حساب تطور حجم الحبيبات في صورة ثلاثية الأبعاد باستخدام مطابقة القالب متعدد السطوح، وتجزئة الحبوب، وتجميع الرسم البياني . في مثل هذه المحاكاة، قدمت طريقة الديناميكا الجزيئية قياسًا دقيقًا لحجم الحبوب. من خلال الاستفادة من هذه المعلومات، تم استخراج الهياكل الفعلية للحبوب وقياسها وتقديمها. مقارنة بالطريقة التقليدية لاستخدام الديناميكات الجزيئية الموجهة (SMD) مع شريحة واحدة ثنائية الأبعاد من المادة، يوفر الديناميكا الجزيئية طريقة ثلاثية الأبعاد ودقيقة لتوصيف تطور البنية المجهرية على النطاق الذري.

خوارزميات الديناميات الجزيئية 

 فحص إمكانات كولوم نموذج المذيب الضمني.

التكامل 

 تكامل عاطفي.
 التكامل بين فيرلت و ستويرمر.
 تكامل طريقة رونج-كوتا.
 خوارزمية بيمان.
 خوارزميات القيد (للأنظمة المقيدة).

خوارزميات التفاعل قصير المدى 

 قوائم الخلايا.
 قائمة فيرلت.
 التفاعلات المستعبدة.

خوارزميات التفاعل بعيد المدى 

 تلخيص ايوالد.
 تجميع ايوالد لشبكة الجسيمات (PME).
 جسيمات - جسيمات - جسيمات - شبكة (P3M).
 طريقة القوة المتغيرة.

استراتيجيات الموازاة 

 طريقة تحليل المجال (توزيع بيانات النظام للحوسبة المتوازية).

الديناميات الجزيئية لكيمياء الكم 

 الديناميكيات الجزيئية لسيارة بارينيلو.

أجهزة متخصصة لمحاكاة الديناميكا الجزيئية 

 أنطون - كمبيوتر عملاق متخصص ومتوازي بشكل كبير مصمم لتنفيذ عمليات محاكاة الديناميكا الجزيئية.
 MDGRAPE - نظام ذو غرض خاص مصمم لمحاكاة الديناميكيات الجزيئية، وخاصة التنبؤ ببنية البروتين.

بطاقة الرسومات كأداة لمحاكاة الديناميكا الجزيئية 

 النمذجة الجزيئية على (GPU).

انظر أيضًا 

 النمذجة الجزيئية.
 الكيمياء الحسابية.
 مجال القوة (كيمياء).
 مقارنة تطبيقات مجال القوة.
 طريقة مونت كارلو.
 برمجيات التصميم الجزيئي.
 ميكانيكا جزيئية.
 متعدد النطاقات وظيفة جرين.
 طريقة كار-بارينيلو.
 مقارنة برامج نمذجة الميكانيكا الجزيئية.
 كيمياء الكم.
 طريقة العنصر المنفصل.
 مقارنة بين برامج محاكاة الحمض النووي.
 محرر جزيء.
 ديناميكيات الكم الكلاسيكية المختلطة.

مراجع

مراجع عامة 

 M. P. Allen, D. J. Tildesley (1989) Computer simulation of liquids. Oxford University Press. .
 J. A. McCammon, S. C. Harvey (1987) Dynamics of Proteins and Nucleic Acids. Cambridge University Press.  (hardback).
 D. C. Rapaport (1996) The Art of Molecular Dynamics Simulation. .
 
 
 J. M. Haile (2001) Molecular Dynamics Simulation: Elementary Methods. 
 R. J. Sadus, Molecular Simulation of Fluids: Theory, Algorithms and Object-Orientation, 2002, 
 Oren M. Becker, Alexander D. Mackerell, Jr., Benoît Roux, Masakatsu Watanabe (2001) Computational Biochemistry and Biophysics. Marcel Dekker. .
 Andrew Leach (2001) Molecular Modelling: Principles and Applications. (2nd Edition) Prentice Hall. .
 Tamar Schlick (2002) Molecular Modeling and Simulation. Springer. .
  (1991) Computational Statistical Mechanics, Elsevier, .
 D. J. Evans and G. P. Morriss (2008) Statistical Mechanics of Nonequilibrium Liquids, Second Edition, Cambridge University Press, .

فيزياء
كيمياء الجزيئات الضخمة
كيمياء حسابية
محاكاة
نمذجة جزيئية
نمذجة علمية
الديناميكا الجزيئية ( بالإنجليزية : Molecular dynamics ( MD ) ) ، هي علم قائم على محاكاة حاسوبية للجزيئات و المركبات معتمدة على قوانين نيوتن ل الحركة في توقع حركة الجزيئات و سرعة تحرك ها و أماكن وجود ها بعد فترات زمنية قصير ه جدا و هو عمل يتطلب حلول لمعادلات نيوتن ل إعداد هائلة من الجزيئات ولفترات زمنية قصير ه جدا ولضمان حساب هذه التوقعات ب شكل صحيح أو اقرب ل الصحة يجب أخذ عدة امور عند حساب قوانين نيوتن مثل قوى التجاذب بين الجزيئات واهتزاز روابط الجزي وانثناء الرووابط حول محور ها وتذبذب الزاوية بين ذرات الجزي و أيضا تذبذب الزوايا بين مستويات الجزئ أو ما يعرف بالزاوية الزوجية وامور أخرى كثيرة ولان نا في هذه العملية نقوب ب حساب قوانين نيوتن و قوانين أخرى تتعلق بالمحاكاة ل إعداد كبيرة جد ه من الجزيئات هذا يجعل العملية معقدة جدا و شبه مستحيلة ل ذٰلك يتم الاستعانة ب بعض الفرضيات و طرق احصائية ل تساعد نا في جعل هذه الحسابات ممكنة من الناحية العملية . الديناميكات الجزيئية هي طريقة محاكاة حاسوبية ل تحليل الحركات الفيزيائية للذرات والجزيئات . يُسمح للذرات والجزيئات بالتفاعل ل فترة زمنية محددة ، مما يعطي نظرة على « التطور » الديناميكي ل النظام . في ال إصدار الأكثر شيوعًا ، يتم تحديد مسارات الذرات والجزيئات عن طريق الحل العددي لمعادلات نيوتن ل الحركة ل نظام من الجسيمات المت فاعلة ، حيث تُحسب القوى بين الجسيمات وطاقات ها الكامنة غالبًا ب استخدام الكمون بين الذرات أو مجالات القوة الجزيئية . يتم تطبيق الطريقة في الغالب في الفيزياء الكيميائية وعلوم المواد والفيزياء الحيوية . نظرًا لأن الأنظمة الجزيئية تتكون عادةً من عدد كبير من الجسيمات ، ف من المستحيل تحديد خصائص هذه الأنظمة المعقدة ب شكل تحليلي ؛ تتغلب محاكاة الديناميكا الجزيئية على هذه المشكلة ب استخدام الطرق العددية . و مع ذلك ، ف إن عمليات محاكاة الديناميكا الجزيئية الطويلة غير مشروطة رياضيًا ، مما يؤدي إلى حدوث أخطاء تراكمية في التكامل العددي يمكن تقليل ها من خلال الاختيار المناسب للخوارزميات والمعلمات ، و لٰكن لا يتم ا التخلص من ها تمامًا . ب النسبة للأنظمة التي تتبع فرضية إرغوديك ، يمكن استخدام تطور محاكاة ديناميكية جزيئية واحدة ل تحديد الخصائص الديناميكية الحرارية العيانية ل النظام : المتوسطات الزمنية ل نظام ergodic تتوافق مع متوسطات المجموعة الدقيقة . كما أطلق على الديناميكا الجزيئية أيضًا اسم « الميكانيكا ال إحصائية بالأرقام » و « رؤية لابلاس للميكانيكا النيوتونية » للتنبؤ ب المستقبل من خلال تحريك قوى الطبيعة و السماح بإلقاء نظرة ثاقبة على الحركة الجزيئية على نطاق ذري . تاريخ تم تطوير الديناميكا الجزيئية في الأصل في أوائل الخمسينيات من القرن الماضي ، بعد النجاحات السابقة في محاكاة مونت كارلو ، و التي يعود تاريخ ها إلى القرن الثامن عشر ، في مشكلة إبرة بوفون على سبيل المثال ، و لٰكن تم نشر ها في الميكانيكا ال إحصائية في مختبر لوس ألاموس الوطني ب واسطة روزنبلوث ومتروبوليس فيما يعرف اليوم بخوارزمية متروبوليس - هاستينغز . يعود الاهتمام ب التطور الزمني ل أنظمة الجسم N إلى القرن السابع عشر ، بدءًا من نيوتن ، واستمر في القرن التالي إلى حد كبير مع التركيز على الميكانيكا السماوية و قضايا مثل استقرار النظام الشمسي . تم تطوير العديد من الأساليب العددية المستخدمة اليوم خلال هذه الفترة الزمنية ، و التي سبقت استخدام أجهزة الكمبيوتر ؛ على سبيل المثال ، خوارزمية التكامل الأكثر شيوعًا المستخدمة اليوم ، خوارزمية تكامل فيرليت ، تم استخدام ها في وقت مبكر من عامب واسطة جان باتيست جوزيف ديلامبر . يمكن اعتبار الحسابات العددية ب استخدام هذه الخوارزميات الديناميكا الجزيئية « يدويًا » . في وقت مبكر من عام تم تنفيذ تكامل معادلات الحركة ل العديد من الأجسام مع أجهزة الكمبيوتر التناظرية . قام البعض بعمل كثيف العمالة لنمذجة الحركة الذرية من خلال بناء نماذج فيزيائية ، على سبيل المثال ، ب استخدام المجالات العيانية . كان الهدف هو ترتيب ها ب طريقة تماثل بنية السائل و استخدام ه ل فحص سلوكه . برنال قال في عام « .\.\ . أخذت عددًا من الكرات المطاطية وغرستهم معًا بقضبان بأطوال مختلفة تتراوح من 75 إلىبوصات . حاولت القيام ب ذٰلك في المقام الأول ب شكل عرضي قدر ال إمكان ، حيث أعمل في مكتبي الخاص ، وأتعرض للمقاطعة كل خمس دقائق أو نحو ذلك و لا أتذكر ما فعلت ه قبل المقاطعة . » بعد اكتشاف الجسيمات المجهرية و تطوير أجهزة الكمبيوتر ، توسع الاهتمام إلى ما وراء إثبات أنظمة الجاذبية إلى الخصائص ال إحصائية للمادة . في محاولة لفهم أصل اللارجعة ، اقترح فيرمي في عام و نشر في عام استخدام مانياك ، أيضًا في مختبر لوس ألاموس الوطني ، ل حل التطور الزمني لمعادلات الحركة ل العديد من : يخضع نظام الجسم ل العديد من خيارات قوانين القوة ؛ اليوم ، يُعرف هذا العمل الأساسي ب مشكلة فيرمي - باستا - أولام - تسينغو . يظهر التطور الزمني ل الطاقة من العمل الأصلي في الشكل إلى اليمين . في عام استخدم بيرني ألدر ووينرايت حاسوب آي بي إملمحاكاة التصادمات المرنة تمامًا بين المجالات الصلبة . في عام ربما في أول محاكاة واقعية للمادة ، جيبسون و آخرون . محاكاة التلف الإشعاعي للنحاس الصلب ب استخدام نوع ولد ماير من التفاعل البغيض مع قوة سطحية متماسكة . في عام نشر عبد الرحمن محاكاة للأرجون السائل تستخدم كمون لينارد جونز . مقارنة حسابات خصائص النظام ، مثل معامل الانتشار الذاتي ، ب شكل جيد مع البيانات التجريبية . مجالات التطبيق وحدود ها تم استخدام طريقة الديناميكا الجزيئية ل أول مرة في الفيزياء ال نظرية ، و قد اكتسبت شعبية في علم المواد بعد ذلك ب وقت قصير ، و منذ السبعينيات أصبحت شائعة أيضًا في الكيمياء الحيوية والفيزياء الحيوية . يستخدم الديناميكا الجزيئية ب شكل متكرر ل تحسين ال هياكل ثلاثية ال أبعاد للبروتينات والجزيئات الكبيرة الأخرى بناءً على القيود التجريبية من علم البلورات بالأشعة السينية أو ال تحليل الطيفي بالرنين المغناطيسي النووي . في الفيزياء ، يتم استخدام الديناميكا الجزيئية ل فحص ديناميكيات الظواهر على المستوى الذري التي لا يمكن ملاحظة ها مباشرة ، مثل نمو الأغشية الرقيقة وزرع الأيونات الفرعية ، و كذٰلك ل فحص الخصائص الفيزيائية للأجهزة النانوية التي لم يتم إنشاء ها أو لا يمكن إنشاء ها بعد . في الفيزياء الحيوية والبيولوجيا الهي كلية ، يتم تطبيق الطريقة ب شكل متكرر ل دراسة حركات الجزيئات الكبيرة مثل البروتينات والأحماض النووية ، و التي يمكن أن تكون مفيدة ل تفسير نتائج بعض التجارب الفيزيائية الحيوية ونمذجة التفاعلات مع الجزيئات الأخرى ، كما هو الحال في الإرساء الترابطي . من حيث المبدأ ، يمكن استخدام الديناميكا الجزيئية للتنبؤ من البداية بهيكل البروتين عن طريق محاكاة طي سلسلة البولي ببتيد من ملف عشوائي . يمكن اختبار نتائج محاكاة الديناميكا الجزيئية من خلال المقارنة مع التجارب التي تقيس الديناميكيات الجزيئية ، و التي تعد طريقة شائعة من ها هي ال تحليل الطيفي بالرنين المغناطيسي النووي . يمكن اختبار تنبؤات الهيكل المشتقة من الديناميكا الجزيئية من خلال تجارب على مستوى المجتمع في التقييم النقدي لتوقع بنية البروتين ( CA SP ) ، على الرغم من أن الطريقة قد حققت نجاحًا محدودًا تاريخيًا في هذا المجال . كتب مايكل ليفيت ، الذي شارك في جائزة نوبل جزئياً ل تطبيق الديناميكا الجزيئية على البروتينات ، في عامأن المشاركين في التقدير التوقعي ل هيكل البروتين لم يستخدموا الطريقة عادةً ب سبب « . . . إحراج مركزي للميكانيكا الجزيئية ، أي تقليل الطاقة أو الديناميات الجزيئية ب شكل عام يؤدي إلى نموذج أقل شبهاً بالهيكل التجريبي . » أدت التحسينات في الموارد الحسابية التي تسمح بمسارات الديناميكا الجزيئية أكثر وأطول ، جنبًا إلى جنب مع التحسينات الحديثة في جودة معلمات مجال القوة ، إلى بعض التحسينات في كل من تنبؤ الهيكل وصقل نموذج التماثل ، دون الوصول إلى نقطة الفائدة العملية في هذه المجالات ؛ يحدد العديد معلمات مجال القوة كمجال رئيسي ل مزيد من التطوير . تم الإبلاغ عن الديناميكا الجزيئية محاكاة ل الخاصة الدوائية تطوير و تصميم الأدوية . على سبيل المثال ، بينتو و آخرون نفذوا عمليات محاكاة الديناميكا الجزيئية لمجمعات ( B cl - Xl ) ل حساب متوسط مواضع الأحماض الأمينية الحرجة المشاركة في ربط الترابط . من ناحية أخرى ، كارلسون و آخرون . نفذت محاكاة الديناميكيات الجزيئية ل تحديد المركبات التي تكمل المستقبل مع التسبب في الحد الأدنى من اضطراب ال تشكل والمرونة في الموقع النشط . تم تراكب لقطات من البروتين على فترات زمنية ثابتة أثناء المحاكاة ل تحديد مناطق الربط المحفوظة ( المحفوظة في ثلاثة إطارات على الأقل من أحد عشر إطارًا ) ل تطوير عالم الأدوية . سبيراكيس و آخرون يعتمد على سير عمل محاكاة الديناميكا الجزيئية ، وبصمات الأصابع للروابط والبروتينات ( FLAP ) و تحليل التمييز الخطي ل تحديد أفضل توافق بروتين يجند ل العمل كقوالب صيدلانية استنادًا إلى تحليل ( ROC ) بأثر رجعي لعقاقير الصيدلانية الناتجة . في محاولة ل تحسين نمذجة اكتشاف الأدوية القائمة على الهيكل ، مقابل الحاجة إلى العديد من المركبات النموذجية ، اقترح هاتمال مزيجًا من محاكاة الديناميكا الجزيئية و تحليل الاتصالات بين الجزيئات لمستقبلات ليجيند لتمييز الاتصالات الحرجة بين الجزيئات ( تفاعلات الربط ) من تلك الزائدة عن الحاجة في مركب واحد بروتين يجند . يمكن بعد ذلك تحويل جهات الاتصال الحرجة إلى نماذج صيدلانية يمكن استخدام ها للفحص الافتراضي . ترتبط حدود الطريقة بمجموعات المعلمات المستخدمة ، و مجالات قوة الميكانيكا الجزيئية الأساسية . تشغيل واحد لمحاكاة الديناميكا الجزيئية يعمل على تحسين الطاقة الكامنة ، بدلاً من الطاقة المجانية للبروتين ، مما يعني أن جميع المساهمات الحتمية في الاستقرار الديناميكي الحرار ي لبنية البروتين مهملة ، ب ما في ذلك الانتروبيا المطابقة ل سلسلة البولي ببتيد ( العامل الرئيسي الذي يزعزع استقرار بنية البروتين ) و التأثيرات الكارهة ل الماء ( القوى الدافعة الرئيسية لطي البروتين ) . عامل مهم آخر هو روابط الهيدروجين داخل الجزيئية ، و التي لم يتم تضمين ها ب شكل صريح في مجالات القوة الحديثة ، و لٰكن ها توصف بتفاعلات كولوم لشحنات النقطة الذرية . هذا تقدير تقريبي لأن الروابط الهيدروجينية ل ها طبيعة ميكانيكية وكيميائية جزئية . علاوة على ذلك ، تُحسب التفاعلات الكهروستاتيكية عادةً ب استخدام ثابت العزل الكهربائي للفراغ ، على الرغم من أن المحلول المائي المحيط ب ه ثابت عازل أعلى ب كثير . ب استخدام العيانية ثابت العزل الكهربائي على مسافات قصيرة بين الذرات غير مشكوك في ه . أخيرًا ، عادةً ما يتم وصف تفاعلات قوى ف إن دير فالس في الديناميكا الجزيئية ب واسطة إمكانات جهد لينارد - جونز استنادًا إلى نظرية فريتز لندن التي لا تنطبق إلا في الفراغ . و مع ذلك ، ف إن جميع أنواع قوى ف إن دير فال هي في نهاية المطاف من أصل إلكتروستاتيكي و ب التالي تعتمد على الخصائص العازلة ل البيئة . يوضح القياس المباشر ل قوى الجذب بين المواد المختلفة ( مثل ثابت هامكر ) أن « التفاعل بين الهيدروكربونات عبر الماء يبلغ حوالي من ذلك عبر الفراغ » . تم إهمال الاعتماد على البيئة ل قوى ف إن دير فال في عمليات المحاكاة القياسية ، و لٰكن يمكن تضمين ها من خلال تطوير مجالات قوة قابلة لل استقطاب . قيود التصميم يجب أن يأخذ تصميم محاكاة الديناميكيات الجزيئية في الاعتبار القوة الحسابية المتاحة . يجب تحديد حجم المحاكاة ( n = عدد الجسيمات ) ، و الوقت ، و المدة الزمنية الإجمالية ب حيث يمكن إنهاء الحساب في غضون فترة زمنية معقولة . و مع ذلك ، يجب أن تكون عمليات المحاكاة طويلة ب ما يكفي ل تكون ذات صلة ب المقاييس الزمنية للعمليات الطبيعية قيد الدراسة . ل الحصول على استنتاجات صحيحة إحصائيًا من عمليات المحاكاة ، يجب أن يتطابق المدى الزمني الذي تمت محاكات ه مع حركية العملية الطبيعية . خلافًا ل ذٰلك ، يكون الأمر مشابهًا ل التوصل إلى استنتاجات حول كيفية سير الإنسان عندما ينظر فقط إلى أقل من خطوة واحدة . تستخدم معظم المنشورات العلمية حول ديناميكيات البروتينات والحمض النووي بيانات من عمليات المحاكاة التي تمتد نانوثانية − 9 ثوانٍ ) إلى ميكروثانية 6 ثوانٍ ) . ل الحصول على هذه المحاكاة ، هناك حاجة إلى عدة أيام من وحدة ا المعالجة المركزية إلى سنوات وحدة ا المعالجة المركزية . تسمح الخوارزميات المتوازية ب توزيع الحمل بين وحدات ا المعالجة المركزية ؛ مثال على ذلك هو خوارزمية التحلل المكاني أو القوة . أثناء محاكاة الديناميكا الجزيئية الكلاسيكية ، ف إن المهمة الأكثر كثافة لوحدة ا المعالجة المركزية هي تقييم الإمكانات كدالة للإحداثيات الداخلية للجسيمات . ضمن تقييم الطاقة هذا ، يكون الجزء الأغلى هو الجزء غير المترابط أو غير التساهمي . في الكبير O التدوين ، وديناميات الجزيئية المشتركة المحاكاة على نطاق و من قبل إذا كان كل الزوج الحكيم كهرباء وفان دير فال التفاعلات يجب أن تمثل ب شكل واضح . يمكن تقليل هذه التكلفة الحسابية عن طريق استخدام طرق الكهرباء الساكنة مثل شبكة الجسيمات بتلخيص ايوالد ( ) أو الجسيمات - الجسيمات - الجسيمات - الشبكة ( P3M ) ، أو طرق القطع الكروية الجيدة ( ) . هناك عامل آخر يؤثر على إجمالي وقت وحدة ا المعالجة المركزية الذي تحتاج ه المحاكاة و هو حجم الخطوة الزمنية ل التكامل . هذه هي المدة الزمنية بين تقييمات الإمكانات . يجب اختيار الخطوة الزمنية صغيرة ب ما يكفي ل تجنب أخطاء التقدير ( أي أصغر من الفترة المتعلقة ب أسرع تردد اهتزازي في النظام ) . الخطوات الزمنية النموذجية لـ الديناميكا الجزيئية الكلاسيكي بترتيبفيمتوثانية - ثانية ) . يمكن توسيع هذه القيمة ب استخدام الخوارزميات مثل خوارزمية قيود شيك ، و التي تثبت اهتزازات الذرات ال أسرع ( على سبيل المثال ، الهيدروجين ) في مكان ها . تم أيضًا تطوير طرق متعددة للمقياس الزمني ، و التي تتيح فترات ممتدة بين تحديثات القوى البعيدة المدى الأبطأ . لمحاكاة الجزيئات في مذيب ، يجب الاختيار بين المذيب الصريح والضمني . يجب حساب جسيمات المذيبات ال صريحة ( مثل و SPC / E وS PC - f ) ب شكل مكلف عن طريق مجال القوة ، بينما تستخدم المذيبات الضمنية نهج المجال المتوسط . يعد استخدام مذيب صريح مكلفًا من الناحية الحسابية ، و يتطلب تضمين ما يقرب من عشرة أضعاف الجسيمات في المحاكاة . لكن حبيبات ولزوجة المذيب الصريح ضروريان ل إعادة إنتاج خصائص معينة للجزيئات الذائبة . هذا مهم ب شكل خاص ل إعادة إنتاج الخواص الحركية الكيميائية . في جميع أنواع محاكاة الديناميكيات الجزيئية ، يجب أن يكون حجم صندوق المحاكاة كبيرًا ب ما يكفي ل تجنب آثار حالة الحدود . غالبًا ما يتم التعامل مع شروط الحدود عن طريق اختيار القيم الثابتة عند الحواف ( و التي قد تسبب تشوهات ) ، أو عن طريق استخدام شروط الحدود الدورية حيث يعود جانب واحد من حلقات المحاكاة إلى الجانب الآخر ، ومحاكاة مرحلة مجمعة ( و التي قد تسبب آثارًا أيضًا ) . مجموعة متناهية الصغر ( NVE ) في المجموعة المتناهية الصغر ، يتم عزل النظام عن التغيرات في الشامات ( N ) والحجم ( V ) و الطاقة ( E ) . يتوافق مع عملية ثابتة الحرارة ب دون تبادل حرار ي . يمكن النظر إلى مسار الديناميات الجزيئية المتناهية الصغر على أن ه تبادل ل الطاقة الكامنة والحركية ، مع الحفاظ على إجمالي الطاقة . ل نظام من الجسيمات N مع الإحداثيات والسرعات ، يمكن كتابة الزوج التالي من المعادلات التفاضلية من الدرجة الأولى في تدوين نيوتن كـ وظيفة الطاقة الكامنة النظام هو دالة إحداثيات الجسيمات \ . يشار إلي ه ببساطة على أن ه الإمكانات في الفيزياء ، أو مجال القوة في الكيمياء . تأتي المعادلة الأولى من قوانين نيوتن ل الحركة . القوة يمكن حساب التأثير على كل جسيم في النظام على أن ه التدرج السالب لـ \ . ل كل خطوة زمنية ، موضع كل جسيم والسرعة قد تتكامل مع طريقة تكامل سيمبلكتيك مثل تكامل فيرلت . التطور الزمني ل و يسمى المسار . ب النظر إلى المواقف الأولية ( على سبيل المثال ، من المعرفة ال نظرية ) والسرعات ( على سبيل المثال ، غاوسي العشوائي ) ، يمكن نا حساب جميع المواقف والسرعات المستقبلية ( أو السابقة ) . أحد مصادر الالتباس المتكررة هو معنى درجة الحرارة في الديناميكا الجزيئية . عادة لدي نا خبرة في درجات الحرارة العيانية ، و التي تنطوي على عدد كبير من الجسيمات . لكن درجة الحرارة هي كمية إحصائية . إذا كان هناك عدد كبير ب ما يكفي من الذرات ، ف يمكن تقدير درجة الحرارة ال إحصائية من درجة الحرارة اللحظية ، و التي يتم العثور علي ها من خلال مساواة الطاقة الحركية ل النظام بـ nk B T / حيث n هو عدد درجات الحرية ل النظام . تنشأ ظاهرة مرتبطة ب درجة الحرارة ب سبب قلة عدد الذرات المستخدمة في محاكاة الديناميكا الجزيئية . على سبيل المثال ، ضع في اعتبارك محاكاة نمو فيلم نحاسي بدءًا من ركيزة تحتوي علىذرة و طاقة ترسيب تبلغفولت . في العالم الحقيقي ، س يتم نقلفولت من الذرة المودعة ب سرعة ومشاركة ها بين عدد كبير من الذرات ( أو أكثر ) دون تغيير كبير في درجة الحرارة . عندما يكون هناكذرة فقط ، ف إن الركيزة تتبخر على الفور تقريبًا ب واسطة الترسب . يحدث شيء مشاب ه في عمليات المحاكاة الفيزيائية الحيوية . يتم رفع درجة حرارة النظام في ( NVE ) ب شكل طبيعي عندما تخضع الجزيئات الكبيرة مثل البروتينات لتغيرات وترابط طارد للحرارة . المجموعة المتعارف علي ها ( NVT ) في المجموعة المتعارف علي ها ، يتم حفظ كمية المادة ( N ) والحجم ( V ) و درجة الحرارة ( T ) . كما يطلق علي ه أحيانًا الديناميكيات الجزيئية ذات درجة الحرارة الثابتة ( CTالديناميكا الجزيئية ) . في ( NVT ) ، يتم تبادل طاقة العمليات الماصة للحرارة والطاردة للحرارة مع منظم الحرارة . تتوفر مجموعة متنوعة من خوارزميات الترموستات ل إضافة الطاقة وإزالة ها من حدود محاكاة الديناميكا الجزيئية ب طريقة أكثر أو أقل واقعية ، ب القرب من المجموعة المتعارف علي ها . تشمل الأساليب الشائعة ل التحكم في درجة الحرارة إعادة قياس السرعة ، ترموستات نوزيه هوفر ، وسلاسل نوزيه هوفر ، وترموستات بيرندسن ، وديناميكيات أندرسن ولانجفين . قد يقدم منظم الحرارة بيرندسن تأثير مكعب الثلج الطائر ، و الذي يؤدي إلى ترجمات غير مادية و تناوب ل النظام المحاكي . ليس من التاف ه الحصول على توزيع متعارف علي ه للمطابقة والسرعات ب استخدام هذه الخوارزميات . كيف يعتمد ذلك على حجم النظام و اختيار منظم الحرارة ومعلمات منظم الحرارة وخطوة الوقت والمتكامل هو موضوع العديد من المقالات في هذا المجال . مجموعة متساوية الحرارة متساوية الضغط ( NPT ) يتم الحفاظ على كمية المادة ( N ) والضغط ( P ) و درجة الحرارة ( T ) في المجموعة متساوي الحرارة - متساوي الضغط . ب الإضافة إلى منظم الحرارة ، هناك حاجة إلى باروستات . يتوافق ب شكل وثيق مع ظروف المختبر مع وجود دورق مفتوح ل درجة الحرارة والضغط المحيطين . في محاكاة الأغشية البيولوجية ، ف إن التحكم في الضغط الخواص غير مناسب . ب النسبة لطبقات الدهون الثنائية ، يحدث التحكم في الضغط تحت منطقة غشاء ثابتة ( NPAT ) أو توتر سطحي ثابت « جاما » ( NPγT ) . المجموعات المعممة طريقة تبادل النسخ المتماثلة هي مجموعة معممة . تم إنشاء ه في الأصل ل التعامل مع الديناميكيات البطيئة ل أنظمة ال دوران المضطربة . يطلق علي ه أيضًا التقسية الموازية . تحاول صياغة الديناميكا الجزيئية ل التبادل المتماثل التغلب على مشكلة الحدود الدنيا المتعددة من خلال تبادل درجة حرارة النسخ المتماثلة غير المت فاعلة ل النظام التي تعمل في درجات حرارة متعددة . الإمكانات في محاكاة الديناميكا الجزيئية تتطلب محاكاة الديناميكيات الجزيئية تعريف وظيفة محتملة ، أو وصفًا للمصطلحات التي تتفاعل ب ها الجسيمات في المحاكاة . في الكيمياء والبيولوجيا ، يُشار إلى هذا عادةً ب اسم مجال القوة و في فيزياء المواد على أن ه جهد بين الذرات . يمكن تحديد الإمكانات على مستويات عديدة من ا الدقة المادية ؛ تلك الأكثر شيوعًا في الكيمياء تعتمد على الميكانيكا الجزيئية وتجسد العلاج الميكانيكي الكلاسيكي لتفاعلات الجسيمات والجسيمات التي يمكن أن تعيد إنتاج التغييرات الهي كلية والتكوينية و لٰكن لا يمكن ها عادةً إعادة إنتاج التفاعلات الكيميائية . يستلزم الاختزال من الوصف الكمي الكامل إلى الإمكانات الكلاسيكية تقريبيين رئيسيين . الأول هو تقريب بورن أوبنهايمر ، الذي ينص على أن ديناميكيات الإلكترونات سريعة جدًا ب حيث يمكن اعتبار ها تتفاعل على الفور مع حركة نواة ها . نتيجة ل ذٰلك ، يمكن معاملت هم ب شكل منفصل . الثاني يتعامل مع النوى ، و هي أثقل ب كثير من الإلكترونات ، على أن ها جسيمات نقطية تتبع الديناميكيات النيوتونية الكلاسيكية . في الديناميكيات الجزيئية الكلاسيكية ، يتم تقريب تأثير الإلكترونات كسطح طاقة كامن ، وعادة ما يمثل الحالة الأرضية . عندما تكون هناك حاجة إلى مستويات أدق من التفاصيل ، يتم استخدام الإمكانات القائمة على ميكانيكا الكم ؛ تحاول بعض الطرق إنشاء إمكانات كلاسيكية / كمومية هجينة حيث يتم معالجة الجزء الأكبر من النظام ب شكل كلاسيكي و لٰكن يتم التعامل مع منطقة صغيرة كنظام كمي ، وعادةً ما تخضع ل عملية تحول كيميائي . الإمكانات التجريبية كثيرا ما تسمى الإمكانات التجريبية المستخدمة في الكيمياء مجالات القوة ، في حين أن تلك المستخدمة في فيزياء المواد تسمى الجهود بين الذرية . معظم مجالات القوة في الكيمياء تجريبية و تتكون من مجموع القوى المرتبطة المرتبطة بالروابط الكيميائية وزوايا ال رابطة وثنائيات السندات و القوى غير المرتبطة المرتبطة ب قوى ف إن دير فال والشحنة الكهروستاتيكية . تمثل الإمكانات التجريبية تأثيرات ميكانيكا الكم ب طريقة محدودة من خلال التقريبات الوظيفية المخصصة . تحتوي هذه الإمكانات على معلمات مجانية مثل الشحنة الذرية ، ومعلمات ف إن دير فالس التي تعكس تقديرات نصف ال قطر الذري ، و طول رابطة التوازن ، والزاوية ، وثنائي السطوح ؛ يتم الحصول علي ها عن طريق التوافق مع الحسابات الإلكترونية التفصيلية ( المحاكاة الكيميائية الكمومية ) أو الخصائص الفيزيائية التجريبية مثل الثوابت المرنة والمعلمات الشبكية والقياسات الطيفية . ب سبب الطبيعة غير المحلية للتفاعلات غير المرتبطة ، ف إن ها تنطوي على تفاعلات ضعيفة على الأقل بين جميع الجسيمات في النظام . عادة ما يكون حساب ه هو عنق الزجاجة في سرعة محاكاة الديناميكا الجزيئية . ل خفض التكلفة الحسابية ، تستخدم حقول القوة تقديرات تقريبية رقمية مثل نصف قطر القطع المغير ، أو خوارزميات مجال التفاعل ، أو تجميع ايوالد لشبكة الجسيمات ، أو شبكة الجسيمات والجسيمات والجسيمات ال أحدث ( P3M ) . تستخدم حقول القوة الكيميائية عادةً ترتيبات ارتباط محددة مسبقًا ( استثناء هو ديناميكيات طرق كيمياء الكم ) ، و ب التالي ف هي غير قادرة على نمذجة عملية كسر ال رابطة الكيميائية والتفاعلات ب شكل صريح . من ناحية أخرى ، يمكن ل العديد من الإمكانات المستخدمة في الفيزياء ، مثل تلك القائمة على ش كلية ترتيب السندات ، أن تصف عدة تنسيقات مختلفة ل النظام و كسر ال رابطة . تتضمن أمثلة هذه الإمكانات إمكانات برينر للهيدروكربونات وتطورات ها الإضافية ل أنظمة C - Si - H و COH . يمكن اعتبار إمكانات هجينًا كامل التفاعل بين إمكانات ترتيب ال رابطة و مجالات القوة الكيميائية . قم بإقران الإمكانات مقابل إمكانات الأجسام المتعددة تتم صياغة الوظائف المحتملة التي تمثل الطاقة غير المرتبطة ك مجموع على التفاعلات بين جسيمات النظام . أبسط خيار ، مستخدم في العديد من مجالات القوة الشائعة ، هو « جهد الزوج » ، حيث يمكن حساب إجمالي الطاقة الكامنة من مجموع مساهمات الطاقة بين أزواج الذرات . ل ذٰلك ، تسمى حقول القوة هذه أيضًا « حقول القوة المضافة » . مثال على مثل هذا الزوج المحتمل هو احتمال جهد لينارد - جونز غير المرتبط ( يُطلق علي ه أيضًا احتمال المستخدم ل حساب قوى ف إن دير فال . مثال آخر هو نموذج Born ( الأيوني ) ل الشبكة الأيونية . المصطلح الأول في المعادلة التالية هو قانون كولوم لزوج من الأيونات ، والمصطلح الثاني هو التنافر قصير المدى الذي شرح ه مبدأ استبعاد باولي والمصطلح الأخير هو مصطلح تفاعل التشتت . عادةً ما تتضمن المحاكاة المصطلح ثنائي القطب فقط ، على الرغم من تضمين المصطلح رباعي الأقطاب في بعض الأحيان . عندما ن ل = وتسمى هذه ال إمكانية أيضا إمكانية كولوم - باكنغهام . في إمكانات الجسم المتعددة ، تشتمل الطاقة الكامنة على تأثيرات ثلاثة أو أكثر من الجسيمات التي تتفاعل مع بعض ها البعض . في عمليات المحاكاة ذات الإمكانات الزوجية ، توجد أيضًا تفاعلات عالمية في النظام ، و لٰكن ها تحدث فقط من خلال المصطلحات الزوجية . في إمكانات الأجسام المتعددة ، لا يمكن العثور على الطاقة الكامنة ب مجموع على أزواج من الذرات ، حيث يتم حساب هذه التفاعلات ب شكل صريح على أن ها مزيج من المصطلحات ذات الترتيب الأعلى . في العرض الإحصائي ، لا يمكن التعبير عن التبعية بين المتغيرات ب شكل عام ب استخدام المنتجات المزدوجة فقط من درجات الحرية . على سبيل المثال ، إمكانات تير سوف ، التي كانت تُستخدم في الأصل لمحاكاة الكربون ، والسيليكون ، والجرمانيوم ، و منذ ذلك الحين تم استخدام ها ل مجموعة واسعة من المواد الأخرى ، تتضمن مجموع مجموعات من ثلاث ذرات ، مع وجود زوايا بين الذرات عامل مهم في الإمكانات . و من ال أمثلة الأخرى طريقة الذرة المدمجة ( EAM ) ، ( E DIP ) ، و إمكانات تقريب اللحظات الثانية المحكم الربط ( TB S M A ) ، حيث تكون كثافة الإلكترون للحالات في منطقة الذرة هي محسوبة من مجموع المساهمات من الذرات المحيطة ، و ب التالي ف إن مساهمة الطاقة الكامنة هي دالة ل هٰذا ال مجموع . الإمكانات شبه التجريبية تستفيد الإمكانات شبه التجريبية من تمثيل المصفوفة من ميكانيكا الكم . و مع ذلك ، يتم العثور على قيم عناصر المصفوفة من خلال الصيغ التجريبية التي تقدر درجة التداخل في المدارات الذرية المحددة . يتم بعد ذلك مائل المصفوفة ل تحديد شغل المدارات الذرية المختلفة ، و يتم استخدام الصيغ التجريبية مرة أخرى ل تحديد مساهمات الطاقة في المدارات . هناك مجموعة متنوعة من الإمكانات شبه التجريبية ، تسمى إمكانات الربط المحكم ، و التي تختلف وفقًا للذرات التي يتم نمذجة ها . إمكانات ال استقطاب تتضمن معظم مجالات القوة الكلاسيكية ضمنيًا تأثير ال استقطاب ، على سبيل المثال ، عن طريق زيادة الشحنات الجزئية التي تم الحصول علي ها من الحسابات الكيميائية الكمومية . هذه الشحنات الجزئية ثابتة ب النسبة إلى كتلة الذرة . و لٰكن يمكن لمحاكاة الديناميات الجزيئية أن تكون نموذجًا صريحًا لقابلية ال استقطاب مع إدخال ثنائيات الأقطاب المستحثة من خلال طرق مختلفة ، مثل الجسيمات الدقيقة أو الشحنات المتقلبة . هذا يسمح ب إعادة التوزيع الديناميكي ل الشحنة بين الذرات و التي تستجيب ل البيئة الكيميائية المحلية . ل سنوات عديدة ، تم وصف عمليات محاكاة الديناميكا الجزيئية القابلة لل استقطاب على أن ها الجيل التالي . ب النسبة للسوائل المتجانسة مثل الماء ، تم تحقيق دقة متزايدة من خلال إدراج قابلية ال استقطاب . كما تم تحقيق بعض النتائج الواعدة للبروتينات . و مع ذلك ، لا يزال من غير المؤكد أفضل طريقة لتقريب ال استقطاب في المحاكاة . الإمكانات في طرق كيمياء الكم في الديناميكيات الجزيئية الكلاسيكية ، يتم تمثيل سطح طاقة كامنة ( عادةً الحالة الأرضية ) في مجال القوة . هذا نتيجة لتقريب بورن أوبنهايمر . في الحالات المثارة أو التفاعلات الكيميائية أو عند الحاجة إلى تمثيل أكثر دقة ، يمكن الحصول على السلوك الإلكتروني من المبادئ الأولى ب استخدام طريقة ميكانيكا الكم ، مثل نظرية الكثافة الوظيفية . يسمى هذا الديناميكيات الجزيئية لكيمياء الكم . نظرًا لتكلفة معالجة درجات الحرية الإلكترونية ، ف إن العبء الحسابي ل هٰذه المحاكاة أعلى ب كثير من الديناميكيات الجزيئية التقليدية . ل هٰذا السبب ، عادةً ما يقتصر نظام الديناميكا الجزيئية على أنظمة أصغر وأوقات أقصر . يمكن استخدام الطرق الميكانيكية والكيميائية الكمومية للمبتدئين ل حساب الطاقة الكامنة ل نظام أثناء الطيران ، حسب الحاجة للتوافق في المسار . يتم إجراء هذا الحساب عادةً في الجوار القريب لإحداثيات التفاعل . على الرغم من أن ه يمكن استخدام تقديرات تقريبية مختلفة ، إلا أن ها تستند إلى اعتبارات نظرية ، و ليس على ملاءمة تجريبية . تنتج حسابات كيمياء الكم قدرًا هائلاً من المعلومات غير المتوفرة من الأساليب التجريبية ، مثل كثافة الحالات الإلكترونية أو الخصائص الإلكترونية الأخرى . من المزايا المهمة ل استخدام طرق كيمياء الكم القدرة على دراسة التفاعلات التي تتضمن كسر أو تكوين روابط تساهمية ، و التي تتوافق مع حالات إلكترونية متعددة . علاوة على ذلك ، تسمح طرق كيمياء الكم أيضًا ب استعادة التأثيرات التي تتجاوز تقريب بورن - أوبنهايمر ب استخدام مناهج مثل الديناميكيات الكمومية الكلاسيكية المختلطة . هجين QM / MM تعتبر طرق QM ( ميكانيكا الكم ) قوية جدًا . و مع ذلك ، ف هي مكلفة من الناحية الحسابية ، في حين أن طرق MM ( الميكانيكا الكلاسيكية أو الجزيئية ) سريعة و لٰكن ها تعاني من عدة حدود ( تتطلب معايير واسعة النطاق ؛ تقديرات الطاقة التي تم الحصول علي ها ليست دقيقة ل الغاية ؛ لا يمكن استخدام ها لمحاكاة التفاعلات حيث يتم كسر / تشكيل الروابط التساهمية ؛ و محدودي قدرات هم على تقديم تفاصيل دقيقة فيما يتعلق بالبيئة الكيميائية ) . ظهرت فئة جديدة من الأساليب التي تجمع بين النقاط الجيدة في حسابات QM ( ا الدقة ) و MM ( السرعة ) . تسمى هذه الطرق طرق ميكانيكا الكم الميكانيكية والجزيئية المختلطة أو المختلطة ( الهجين QM / MM ) . إن أهم ميزة ل طريقة إدارة الجودة ( QM / MM ) الهجينة هي السرعة . تكلفة عمل الديناميكيات الجزيئية الكلاسيكية ( MM ) في مقاييس الحالة الأكثر وضوحًا O ( n حيث n هو عدد الذرات في النظام . هذا يرجع أساسًا إلى مصطلح التفاعلات الكهروستاتيكية ( يتفاعل كل جسيم مع كل جسيم آخر ) . و مع ذلك ، ف إن استخدام نصف قطر القطع والتحديثات الدورية لقائمة الأزواج ومؤخرًا الاختلافات في طريقة شبكة الجسيمات ايوالد ( P ME ) قد قلل هذا إلى ما بين O ( n ) إلى O ( n ب معنى آخر ، إذا تمت محاكاة نظام ب ه ضعف عدد الذرات ، فس يستغرق ما بين ضعفين إلى أربعة أضعاف قوة الحوسبة . من ناحية أخرى ، ف إن أبسط حسابات كيمياء الكم مقياس نموذجي O ( nأو أسوأ ( تم اقتراح حسابات هارتري فوك المقيدة ل توسيع نطاق ~ O ( n7 ) ) . ل التغلب على هذا الحد ، يتم التعامل مع جزء صغير من النظام ميكانيكيًا ( عادةً موقع نشط من الإنزيم ) و يتم معالجة النظام الم تبقي ب شكل كلاسيكي . في التطبيقات الأكثر تعقيدًا ، توجد طرق QM / MM ل معالجة كل من النوى الخفيفة المعرضة للتأثيرات الكمومية ( مثل الهيدروجين ) والحالات الإلكترونية . هذا يسمح بتوليد وظائف موجات الهيدروجين ( على غرار وظائف الموجة الإلكترونية ) . كانت هذه المنهجية مفيدة في التحقيق في ظواهر مثل نفق الهيدروجين . أحد ال أمثلة حيث قدمت طرق QM / MM اكتشافات جديدة هو حساب نقل الهيدريد في إنزيم نازع هيدروجين الكحول في الكبد . في هذه الحالة ، نفق الكم مهم للهيدروجين ، لأن ه يحدد معدل التفاعل . الخشنة الحبيبات والتمثيلات المخفضة في الطرف الآخر من مقياس التفاصيل توجد نماذج ذات حبيبات خشنة ونماذج شعرية . بدلاً من التمثيل الصريح ل كل ذرة في النظام ، يستخدم المرء « ذرات زائفة » ل تمثيل مجموعات من الذرات . قد تتطلب عمليات محاكاة الديناميكا الجزيئية على أنظمة كبيرة جدًا موارد كمبيوتر كبيرة ب حيث لا يمكن دراسة ها ب سهولة من خلال طرق كل الذرة التقليدية . وب المثل ، ف إن محاكاة العمليات على نطاقات زمنية طويلة ( تتجاوز حواليميكرو ثانية ) باهظة التكلفة ، لأن ها تتطلب الكثير من الخطوات الزمنية . في هذه الحالات ، يمكن للمرء أحيانًا معالجة المشكلة ب استخدام تمثيلات مختصرة ، و التي تسمى أيضًا نماذج الحبيبات الخشنة . من ال أمثلة على طرق الحبيبات الخشنة ( CG ) الديناميكيات الجزيئية المتقطعة ( CG - Dالديناميكا الجزيئية ) ونماذج Go . يتم إجراء الحبيبات الخشنة أحيانًا ب أخذ ذرات زائفة أكبر . تم استخدام مثل هذه التقريبات الذرية الموحدة في محاكاة الديناميكا الجزيئية للأغشية البيولوجية . قد يكون تنفيذ مثل هذا ال نهج على الأنظمة التي تكون في ها الخصائص الكهربائية ذات أهمية أمرًا صعبًا نظرًا ل صعوبة استخدام توزيع الشحنة المناسب على الذرات الزائفة . يتم تمثيل ذيول الليبيدات الأليفاتية ب واسطة عدد قليل من الذرات الزائفة عن طريق تجميعإلىمجموعات ميثيلين في كل ذرة زائفة . يجب أن يتم تحديد معلمات هذه النماذج الخشنة ل الغاية ب شكل تجريبي ، من خلال مطابقة سلوك النموذج مع البيانات التجريبية المناسبة أو محاكاة جميع الذرة . من الناحية المثالية ، يجب أن تأخذ هذه المعلمات في الاعتبار المساهمات الحرارية والنتروبية في الطاقة الحرة ب طريقة ضمنية . عندما يتم إجراء الحبيبات الخشنة على مستويات أعلى ، قد تكون دقة الوصف الديناميكي أقل موثوقية . و لٰكن تم استخدام النماذج ذات الحبيبات الخشنة ب نجاح ل فحص مجموعة واسعة من الأسئلة في البيولوجيا الهي كلية و تنظيم البلورات السائلة وزجاج البوليمر . أمثلة على تطبيقات الحبيبات الخشنة : غالبًا ما تُجرى دراسات طي البروتين ودراسات التنبؤ ببنية البروتين ب استخدام ذرة زائفة واحدة أو بضع ذرات ل كل حمض أميني . تم فحص انتقالات الطور البلوري السائل في الأشكال الهندسية المحصورة و / أو أثناء التدفق ب استخدام إمكانات جاي بيرن ، التي تصف الأنواع متباينة الخواص . تمت دراسة زجاج البوليمر أثناء التشوه ب استخدام نوابض توافقية بسيطة أو نوابض فين ل توصيل الكرات الموصوفة ب واسطة إمكانات لينارد جونز . تم فحص الالتفاف الفائق للحمض النووي باستخدامذرات زائفة ل كل زوج أساس ، وبدقة أقل . تم فحص تغليف الحمض النووي المزدوج الحلزوني في العاثية بنماذج حيث تمثل الذرة الزائفة منعطفًا واحدًا ( حواليأزواج أساسية ) من الحلزون المزدوج . تم تصميم بنية الحمض النووي الريبي في الريبوسوم والأنظمة الكبيرة الأخرى ب استخدام ذرة زائفة واحدة ل كل نوكليوتيد . محاكاة افتراضية للخلية ل دراسة تفاعل الخلايا والركائز المختلفة . أبسط شكل من أشكال الحبيبات الخشنة هو الذرة المتحدة ( تسمى أحيانًا الذرة الممتدة ) ، و قد تم استخدام ها في معظم عمليات محاكاة الديناميكا الجزيئية المبكرة للبروتينات والدهون والأحماض النووية . على سبيل المثال ، بدلاً من معالجة جميع الذرات الأربع ل مجموعة ميثيل CH ب شكل صريح ( أو جميع الذرات الثلاث ل مجموعة الميثيلين CH يمثل المرء المجموعة ب أكمل ها بذرة زائفة واحدة . يجب ب الطبع تحديد معلمات ها ب شكل صحيح ب حيث تتمتع تفاعلات ها مع المجموعات الأخرى ب الاعتماد المناسب على المسافة . تنطبق اعتبارات مماثلة على الروابط و الزوايا والالتواءات التي تشارك في ها الذرة الزائفة . في هذا النوع من التمثيل الذري الموحد ، يزيل المرء عادةً جميع ذرات الهيدروجين ال صريحة ب استثناء تلك التي لدي ها القدرة على المشاركة في الروابط الهيدروجينية ( الهيدروجين القطبي ) . مثال على ذلك هو مجال القوة تشارم عادةً ما يتم ا الاحتفاظ بالهيدروجين القطبي في النموذج ، لأن ا المعالجة المناسبة للروابط الهيدروجينية تتطلب وصفًا دقيقًا ب شكل معقول للاتجاه والتفاعلات الكهروستاتيكية بين مجموعات المتبرع و المستقبل . يمكن ل مجموعة الهيدروكسيل ، على سبيل المثال ، أن تكون مانحًا ل رابطة الهيدروجين ، ومتقبلًا ل رابطة الهيدروجين ، و س يكون من المستحيل معالجة ذلك ب استخدام ذرة زائفة ( OH ) واحدة . حوالي نصف الذرات الموجودة في البروتين أو الحمض النووي هي هيدروجين غير قطبي ، لذا ف إن استخدام الذرات الموحدة يمكن أن يوفر توفيرًا كبيرًا في وقت الكمبيوتر . دمج تأثيرات المذيبات في العديد من عمليات محاكاة نظام المذيب المذاب ، ينصب التركيز الرئيسي على سلوك المذاب مع القليل من الاهتمام بسلوك المذيب خاصة في جزيئات المذيب الموجودة في مناطق بعيدة عن الجزيء المذاب . قد تؤثر المذيبات على السلوك الديناميكي للمذابات من خلال الاصطدامات العشوائية و فرض مقاومة احتكاكية على حركة المذاب عبر المذيب . يمكن أن يساعد استخدام شروط الحدود الدورية غير المستطيلة و الحدود العشوائية وأغلفة المذيبات في تقليل عدد جزيئات المذيبات المطلوبة و تمكين نسبة أكبر من وقت الحوسبة ل يتم إنفاق ها بدلاً من ذلك في محاكاة المذاب . من الممكن أيضًا دمج تأثيرات المذيب دون الحاجة إلى وجود أي جزيئات مذيب واضحة . أحد ال أمثلة على هذا ال نهج هو استخدام القوة المتوسطة المحتملة ( PMF ) التي تصف كيفية تغير الطاقة الحرة كتنوع إحداثيات معينة . يحتوي تغيير الطاقة الحرة الموصوف ب واسطة ( PMF ) على متوسط تأثيرات المذيب . قوى بعيدة المدى التفاعل طويل المدى هو تفاعل لا يسقط في ه التفاعل المكاني أسرع من أين هي أبعاد النظام . تشمل ال أمثلة تفاعلات الشحنة والشحنة بين الأيونات وتفاعلات ثنائي القطب بين الجزيئات . تمثل نمذجة هذه القوى تحديًا كبيرًا لأن ها كبيرة على مسافة قد تكون أكبر من نصف طول الصندوق مع محاكاة ل عدة آلاف من الجسيمات . على الرغم من أن أحد الحلول س يكون زيادة حجم طول الصندوق ب شكل كبير ، إلا أن أسلوب القوة الغاشمة هذا أقل من مثالي لأن المحاكاة س تصبح باهظة الثمن من الناحية الحسابية . إن اقتطاع الإمكانات كرويًا هو أيضًا غير وارد حيث يمكن ملاحظة السلوك غير الواقعي عندما تكون المسافة قريبة من مسافة القطع . الديناميكات الجزيئية الموجهة ( SMD ) تقوم محاكاة الديناميكيات الجزيئية ( Sالديناميكا الجزيئية ) ، أو محاكاة مسبار القوة ، ب تطبيق قوى على بروتين من أجل التلاعب بهيكله عن طريق سحب ه على طول درجات الحرية المرغوبة . يمكن استخدام هذه التجارب ل الكشف عن التغيرات الهي كلية في البروتين على المستوى الذري . غالبًا ما يستخدم Sالديناميكا الجزيئية لمحاكاة أحداث مثل ال تكشف أو التمدد الميكانيكي . هناك نوعان من البروتوكولات النموذجية لـ الديناميكا الجزيئية : أحد هما تكون في ه سرعة السحب ثابتة ، و الآخر تكون في ه القوة المطبقة ثابتة . ب شكل نموذجي ، يتم تقييد جزء من النظام المدروس ( على سبيل المثال ، ذرة في بروتين ) ب واسطة جهد توافقي . ثم يتم تطبيق القوى على ذرات معينة إما ب سرعة ثابتة أو ب قوة ثابتة . يتم استخدام أخذ العينات الشامل ل تحريك النظام على طول تنسيق التفاعل المطلوب عن طريق تغيير ، على سبيل المثال ، القوى والمسافات و الزوايا التي تم التلاعب ب ها في المحاكاة . من خلال أخذ العينات المظلة ، يتم أخذ عينات من جميع تكوينات النظام - عالية الطاقة ومنخفضة الطاقة - ب شكل كافٍ . بعد ذلك ، يمكن حساب تغيير كل تكوين في الطاقة الحرة على أن ه احتمال متوسط القوة . إحدى الطرق الشائعة ل حساب ( PMF ) هي من خلال طريقة تحليل الرسم البياني الموزون ( WHAM ) ، و التي تحلل سلسلة من محاكاة أخذ العينات الشاملة . يوجد الكثير من التطبيقات المهمة لـ Sالديناميكا الجزيئية في مجال اكتشاف الأدوية وعلوم الجزيئات الحيوية . على سبيل المثال تم استخدام Sالديناميكا الجزيئية ل التحقيق في استقرار ليفية الزهايمر الأولية ، ل دراسة تفاعل بروتين يجند في كيناز المعتمد على السيكلين و حتى لإظهار تأثير المجال الكهربائي على مركب الثرومبين ( البروتين ) والأبتامر ( النيوكليوتيد ) بين العديد من الدراسات الأخرى المثيرة للاهتمام . أمثلة على التطبيقات تستخدم الديناميات الجزيئية في العديد من مجالات العلوم . تم نشر أول محاكاة الديناميكا الجزيئية ل عملية الطي البيولوجي المبسطة في عام مهدت المحاكاة التي نُشرت في دورية نتشر الطريق لمساحة واسعة من طي البروتين الحسابي الحديث . تم نشر أول محاكاة الديناميكا الجزيئية ل عملية بيولوجية في عام مهدت المحاكاة التي نُشرت في دورية نتشر الطريق لفهم حركة البروتين ب اعتبار ها ضرورية في الوظيفة و ليس مجرد ملحق . الديناميكا الجزيئية هي الطريقة القياسية ل معالجة شلالات الاصطدام في نظام ارتفاع الحرارة ، أي التأثيرات التي تحدث ها إشعاعات النيوترونات والأيونات النشطة على المواد الصلبة والأسطح الصلبة . توضح ال أمثلة الفيزيائية الحيوية التالية الجهود الملحوظة ل إنتاج محاكاة ل أنظمة ذات حجم كبير جدًا ( فيروس كامل ) أو أوقات محاكاة طويلة جدًا ( تصل إلى 112 مللي ثانية ) : محاكاة الديناميكا الجزيئية لفيروس فسيفساء التبغ الساتلي الكامل ( ST M V ) الحجم : مليون ذرة ، وقت المحاكاة : نانوثانية ، البرنامج : NAالديناميكا الجزيئية ) هذا الفيروس عبارة عن فيروس نباتي صغير عشروني الوجوه يؤدي إلى تفاقم أعراض الإصابة بفيروس التبغ الفسيفسائي ( T M V ) . تم استخدام محاكاة الديناميكيات الجزيئية ل استكشاف آليات التجمع الفيروسي . يتكون جسيم ST M V بأكمل ه مننسخة متطابقة من بروتين واحد يتكون من القفيصة الفيروسية ( الغلاف ) ، وجينوم RNA أحادي النوكليوتيدتقطعت ب ه السبل . تتمثل إحدى النتائج الرئيسية في أن الكابسيد يكون غير مستقر ل الغاية عند عدم وجود RNA بداخل ه . س تستغرق المحاكاة جهاز كمبيوتر مكتبي واحدًا في عامحواليعامًا حتى تكتمل . وهكذا تم القيام ب ه في العديد من المعالجات بالتوازي مع ا التواصل المستمر بين هم . محاكاة قابلة للطي لـ فيلين هيدبيس بتفاصيل شاملة الحجم : ذرة ؛ وقت المحاكاة : ميكرو ثانية = نانوثانية ، البرنامج : فولدنغ @ هوم ) تم تشغيل هذه المحاكاة فيوحدة معالجة مركزية من أجهزة الكمبيوتر الشخصية المشاركة حول العالم . تم تثبيت برنامج فولدنغ @ هوم على أجهزة الكمبيوتر هذه ، و هو جهد حوسبي موزع على نطاق واسع ب تنسيق من فيجاي باندي في جامعة ستانفورد . تم فحص الخصائص الحركية لبروتين فيلين هيدبيس ب استخدام العديد من المسارات المستقلة والقصيرة التي تدير ها وحدة ا المعالجة المركزية دون اتصال مستمر في الوقت الحقيقي . كانت إحدى الطرق المستخدمة هي تحليل قيمة بفولد ، و التي تقيس احتمالية الطي قبل الكشف عن تشكيل بداية محدد . يعطي بفولد معلومات حول هياكل الحالة الانتقالية و ترتيب المطابقة على طول مسار الطي . يمكن أن يكون كل مسار في حساب بفولد قصيرًا نسبيًا ، و لٰكن هناك حاجة إلى العديد من المسارات المستقلة . تم إجراء عمليات محاكاة طويلة للمسار المستمر على انطون ، و هو كمبيوتر عملاق متوازي ب شكل كبير تم تصميم ه وبناء ه حول دوائر متكاملة خاصة بالتطبيقات المخصصة ( A SIC ) ووصلات متصلة ب واسطة دي شو ل الأبحاث . أطول نتيجة تم نشر ها لمحاكاة تم إجراء ها ب استخدام أنطون هي محاكاة N TL9 تبلغ 112 مللي ثانية عندكلفن ؛ تم أيضًا إجراء محاكاة ثانية مستقلة تبلغ 073 مللي ثانية ل هٰذا التكوين ( والعديد من عمليات المحاكاة الأخرى ل أكثر منميكرو ثانية من الوقت الكيميائي المستمر ) . في كيف تطوى البروتينات سريعة الطي ، يناقش الباحثون كريستين ليندورف لارسن وستيفانو بيانا ورون أو درور وديفيد إي . تكشف عن مجموعة من المبادئ المشتركة الكامنة وراء طيبروتينًا متنوعًا هيكليًا " . يسمح فحص هذه المسارات الطويلة المتنوعة ، التي تم تمكين ها ب واسطة أجهزة مخصصة ومخصصة ، باستنتاج أن ه " في معظم الحالات ، يتبع الطي مسارًا واحدًا مهيمنًا تظهر في ه عناصر الهيكل الأصلي ب ترتيب مرتبط ب شدة بميل ها لل تشكل في ت كشفت الدولة " . في دراسة منفصلة ، تم استخدام أنطون ل إجراء محاكاة 013 مللي ثانية لديناميات الحالة الأصلية لمثبط التربسين البنكرياس البقري ( BPTI ) عندكلفن . يستفيد تطبيق مهم آخر ل طريقة الديناميكا الجزيئية من قدرة ها على التوصيف ثلاثي ال أبعاد و تحليل التطور المجهرية على النطاق الذري . تُستخدم محاكاة الديناميكا الجزيئية في توصيف تطور حجم الحبوب ، على سبيل المثال ، عند وصف التآكل والاحتكاك لمواد Al و Al ( Zr ) النانوية . يتم تحليل تطور الخلع و تطور حجم الحبيبات أثناء عملية الاحتكاك في هذه المحاكاة . نظرًا لأن طريقة الديناميكا الجزيئية قدمت المعلومات الكاملة للبنية المجهرية ، ف قد تم حساب تطور حجم الحبيبات في صورة ثلاثية ال أبعاد ب استخدام مطابقة القالب متعدد السطوح ، وتجزئة الحبوب ، و تجميع الرسم البياني \ . في مثل هذه المحاكاة ، قدمت طريقة الديناميكا الجزيئية قياسًا دقيقًا ل حجم الحبوب . من خلال الاستفادة من هذه المعلومات ، تم استخراج ال هياكل الفعلية للحبوب وقياس ها وتقديم ها . مقارنة ب الطريقة التقليدية ل استخدام الديناميكات الجزيئية الموجهة ( SMD ) مع شريحة واحدة ثنائية ال أبعاد من المادة ، يوفر الديناميكا الجزيئية طريقة ثلاثية ال أبعاد و دقيقة لتوصيف تطور البنية المجهرية على النطاق الذري . خوارزميات الديناميات الجزيئية فحص إمكانات كولوم نموذج المذيب الضمني . التكامل تكامل عاطفي . التكامل بين فيرلت و ستويرمر . تكامل طريقة رونج - كوتا . خوارزمية بيمان . خوارزميات القيد ( للأنظمة المقيدة ) . خوارزميات التفاعل قصير المدى قوائم الخلايا . قائمة فيرلت . التفاعلات المستعبدة . خوارزميات التفاعل ب عيد المدى تلخيص ايوالد . تجميع ايوالد لشبكة الجسيمات ( P ME ) . جسيمات - جسيمات - جسيمات - شبكة ( P3M ) . طريقة القوة المتغيرة . استراتيجيات الموازاة طريقة تحليل المجال ( توزيع بيانات النظام للحوسبة المتوازية ) . الديناميات الجزيئية لكيمياء الكم الديناميكيات الجزيئية لسيارة بارينيلو . أجهزة متخصصة لمحاكاة الديناميكا الجزيئية أنطون - كمبيوتر عملاق متخصص ومتوازي ب شكل كبير مصمم ل تنفيذ عمليات محاكاة الديناميكا الجزيئية . MDG RAPE - نظام ذو غرض خاص مصمم لمحاكاة الديناميكيات الجزيئية ، و خاصة التنبؤ ببنية البروتين . بطاقة ال رسومات ك أداة لمحاكاة الديناميكا الجزيئية النمذجة الجزيئية على ( GPU ) . انظر أيضًا النمذجة الجزيئية . الكيمياء الحسابية . مجال القوة ( كيمياء ) . مقارنة تطبيقات مجال القوة . طريقة مونت كارلو . برمجيات التصميم الجزيئي . ميكانيكا جزيئية . متعدد النطاقات وظيفة جرين . طريقة كار - بارينيلو . مقارنة برامج نمذجة الميكانيكا الجزيئية . كيمياء الكم . طريقة العنصر المنفصل . مقارنة بين برامج محاكاة الحمض النووي . محرر جزيء . ديناميكيات الكم الكلاسيكية المختلطة . مراجع مراجع عامة M . P . Allen , D . J . Tildesley Computer simulation of liquids . Oxford University Press.\ . J . A . Mc Camm on , S . C . HarveyDynamics of Proteins and Nucleic A cids . Cambridge University Press . ( hardback ) . D . C . RapaportThe Art of Molecular Dynamics Simulati o n .\ . J . M . HaileM olecular Dynamics Simulation : Elementary Meth ods . R . J . Sadus , Molecular Simulation of Fluids : Theory , Algorithms and Object - Orientation , Oren M . Becker , Alexander D . Macker e ll , Jr . , Benoît Roux , Masakatsu Watanabe Computatio nal Biochemistry and Bioph ysics . Marcel Dekke r.\ . A ndrew Leach M olecular Modelling : Principles and A pplicati o n s.d Edition ) Prentice H all.\ . Tamar SchlickM olecular Modeling and Simulation . Springe r .\ . Computatio nal Statistic al Mech anics , Elsevie r ,\ . D . J . Evans and G . P . MorrissStatistic al Mech anics of Nonequilibrium Liquids , Second Edition , Cambridge University Press,\ . فيزياء كيمياء الجزيئات الضخمة كيمياء حسابية محاكاة نمذجة جزيئية نمذجة علمية