طريقة متروبوليس مونت كارلو في الفيزياء
تُعد خوارزمية متروبوليس الأداة الأساسية في محاكاة الفيزياء الإحصائية: فمن خلال قبول أو رفض التحركات المقترحة بناءً على تكلفتها الطاقوية، تبني سلسلة ماركوف تُنمذج التكوينات باحتمالية بولتزمان الصحيحة.
Definition
خوارزمية متروبوليس هي طريقة مونت كارلو لسلسلة ماركوف تُنشئ تسلسلاً من التكوينات التي يكون توزيعها الحدي هو التجمع القانوني (canonical ensemble)، وذلك باقتراح تغييرات موضعية وقبولها باحتمالية يحددها عامل بولتزمان لتغير الطاقة.
Scope
يغطي هذا الموضوع خوارزميات متروبوليس ومتروبوليس-هاستينغز كما تُطبق على الأنظمة الفيزيائية: قاعدة القبول، التوازن التفصيلي والشمولية (ergodicity)، الاستقرار (equilibration) والارتباط الذاتي، وتقدير المتوسطات الحرارية وأخطائها الإحصائية. إنها طريقة النمذجة الأساسية التي تقوم عليها منطقة مونت كارلو الأوسع.
Core questions
- كيف تعتمد احتمالية القبول على تغير الطاقة للحركة المقترحة؟
- لماذا يضمن التوازن التفصيلي التوزيع الثابت الصحيح؟
- كيف يتم تشخيص أوقات الاستقرار والارتباط الذاتي وأخذها في الاعتبار؟
- كيف يُقدر الخطأ الإحصائي لمتوسط مونت كارلو من العينات المترابطة؟
Key theories
- التوازن التفصيلي والثبات
- يضمن اختيار احتمالات القبول التي تحقق التوازن التفصيلي بالنسبة لتوزيع بولتزمان أن يكون هذا التوزيع ثابتًا ضمن سلسلة ماركوف، بحيث تتقارب المتوسطات طويلة المدى مع قيم التوقع الحرارية.
- تعميم متروبوليس-هاستينغز
- عمم هاستينغز قاعدة القبول لتشمل توزيعات الاقتراح غير المتماثلة، مما وسع نطاق الخوارزمية لتتجاوز التحركات الموضعية المتماثلة مع الحفاظ على التوزيع الثابت المستهدف.
- الارتباط الذاتي وتقدير الخطأ
- تكون عينات متروبوليس المتتالية مترابطة، لذا يُقلل العدد الفعال للعينات المستقلة بواسطة زمن الارتباط الذاتي، والذي يجب قياسه لتحديد هوامش خطأ دقيقة للمتوسطات الحرارية.
Clinical relevance
تُستخدم طريقة متروبوليس لحساب الكميات الديناميكية الحرارية لنماذج اللف المغزلي الشبكي (lattice spin models)، والسوائل والبوليمرات، وتحديد التحولات الطورية، وتعمل كمحرك أساسي ضمن محاكاة مونت كارلو الجزيئية والعديد من مخططات مونت كارلو الكمومية.
History
قُدمت في عام 1953 لحساب معادلة الحالة لسائل قرص صلب ثنائي الأبعاد على حاسوب MANIAC في لوس ألاموس، ثم عُممت الخوارزمية بواسطة هاستينغز في عام 1970 وأصبحت طريقة المحاكاة الأكثر استخدامًا في الفيزياء الإحصائية، ولاحقًا في الإحصاء البايزي.
Key figures
- Nicholas Metropolis
- Arianna Rosenbluth
- W. Keith Hastings
Related topics
Seminal works
- metropolis1953
- hastings1970
Frequently asked questions
- لماذا تُقبل دائمًا التحركات التي تُخفض الطاقة؟
- تزيد الحركة التي تُخفض الطاقة من وزن بولتزمان، لذا فإن قبولها دائمًا يُحرك السلسلة نحو حالات أكثر احتمالية؛ تُقبل التحركات الصاعدة أحيانًا فقط، باحتمالية يحددها زيادة الطاقة، وهذا ما يسمح للسلسلة باستكشاف التوزيع الحراري الكامل بدلاً من الانحدار فقط.
- لماذا يجب التخلص من العينات في بداية التشغيل؟
- تبدأ السلسلة من تكوين عشوائي لا يمثل بعد توزيع التوازن. تُهمل فترة الاستقرار الأولية أو فترة الإحماء (burn-in) بحيث تعكس المتوسطات المقاسة التجمع الحراري الحقيقي بدلاً من التحيز الأولي.