วิธีการมอนติคาร์โลและพลังงานอิสระ

วิธีการมอนติคาร์โลสุ่มตัวอย่างการจัดเรียงโครงสร้างแบบสุ่ม และเทคนิคพลังงานอิสระที่สร้างขึ้นจากวิธีการเหล่านี้จะคำนวณปริมาณทางเทอร์โมไดนามิกที่ควบคุมการจับกัน การละลาย และสมดุล

Definition

วิธีการสุ่มตัวอย่างการจัดเรียงโครงสร้างแบบสุ่มและตัวประมาณค่าพลังงานอิสระที่ได้จากวิธีการเหล่านี้ ซึ่งใช้ในการคำนวณความแตกต่างของพลังงานอิสระทางเทอร์โมไดนามิกในระบบโมเลกุล

Scope

ครอบคลุมการสุ่มตัวอย่างแบบ Metropolis Monte Carlo, การสุ่มตัวอย่างแบบความสำคัญ (importance sampling) และสมดุลละเอียด (detailed balance) รวมถึงเทคนิคพลังงานอิสระหลัก ได้แก่ การรบกวนพลังงานอิสระ (free-energy perturbation), การรวมทางเทอร์โมไดนามิก (thermodynamic integration), การสุ่มตัวอย่างแบบร่ม (umbrella sampling) และวิธีการสุ่มตัวอย่างแบบเสริมประสิทธิภาพที่สามารถเอาชนะอุปสรรคสูงได้ เน้นการคำนวณพลังงานอิสระสัมพัทธ์และสัมบูรณ์ที่มีประโยชน์ในทางปฏิบัติ

Core questions

• การสุ่มตัวอย่างแบบ Metropolis สร้างการจัดเรียงโครงสร้างที่มีน้ำหนัก Boltzmann ที่ถูกต้องได้อย่างไร?
• เหตุใดพลังงานอิสระจึงคำนวณยากกว่าพลังงาน และวิธีการรบกวนและการรวมจัดการกับปัญหานี้อย่างไร?
• การสุ่มตัวอย่างแบบร่มกู้คืนโปรไฟล์พลังงานอิสระข้ามอุปสรรคได้อย่างไร?
• กลยุทธ์ใดบ้างที่ช่วยเพิ่มประสิทธิภาพการสุ่มตัวอย่างเหตุการณ์ที่หายาก?

Key theories

การสุ่มตัวอย่างแบบ Metropolis

ยอมรับหรือปฏิเสธการเคลื่อนที่ทดลองด้วยความน่าจะเป็นที่สอดคล้องกับสมดุลละเอียด (detailed balance) เพื่อให้การจัดเรียงโครงสร้างที่สร้างขึ้นเป็นไปตามการกระจายแบบ Boltzmann ซึ่งช่วยให้สามารถหาค่าเฉลี่ยทางเทอร์โมไดนามิกที่ไม่ลำเอียงได้

การรบกวนพลังงานอิสระ

แสดงความแตกต่างของพลังงานอิสระระหว่างสองสถานะเป็นค่าเฉลี่ยของกลุ่ม (ensemble average) ของความแตกต่างของพลังงานแบบเอกซ์โพเนนเชียล ซึ่งเป็นรากฐานของการคำนวณพลังงานอิสระแบบแปรธาตุ (alchemical free-energy calculations)

Clinical relevance

วิธีการพลังงานอิสระเป็นพื้นฐานของการทำนายเชิงปริมาณของความสัมพันธ์ในการจับกัน (binding affinities), พลังงานอิสระของการละลาย (solvation free energies) และสัมประสิทธิ์การแบ่งส่วน (partition coefficients) ทำให้เป็นหัวใจสำคัญของการค้นพบยาด้วยคอมพิวเตอร์และการประมาณค่าคุณสมบัติทางกายภาพ

History

อัลกอริทึม Metropolis ในปี 1953 ได้นำเสนอการสุ่มตัวอย่างแบบความสำคัญ (importance-sampling Monte Carlo); สูตรการรบกวนของ Zwanzig ในปี 1954 และแผนการรวมทางเทอร์โมไดนามิกและการสุ่มตัวอย่างแบบร่มในภายหลังได้สร้างชุดเครื่องมือพลังงานอิสระสมัยใหม่ที่ปัจจุบันเป็นมาตรฐานในการจำลองโมเลกุล

Key figures

• Nicholas Metropolis
• Marshall Rosenbluth
• Robert Zwanzig
• Daan Frenkel

Seminal works

• metropolis1953
• zwanzig1954

Frequently asked questions

Monte Carlo และพลวัตโมเลกุลแตกต่างกันอย่างไร?

Monte Carlo สุ่มตัวอย่างการจัดเรียงโครงสร้างแบบสุ่มโดยไม่มีเวลาจริงหรือพลวัต ในขณะที่พลวัตโมเลกุลสร้างวิถีการเคลื่อนที่ตามลำดับเวลา ทั้งสองวิธีสามารถให้ค่าเฉลี่ยสมดุลได้ แต่มีเพียงพลวัตเท่านั้นที่ให้คุณสมบัติที่ขึ้นกับเวลา

เหตุใดการคำนวณพลังงานอิสระจึงเป็นเรื่องยาก?

พลังงานอิสระขึ้นอยู่กับปริมาตรทั้งหมดของปริภูมิเฟสที่เข้าถึงได้ แทนที่จะเป็นเพียงการจัดเรียงโครงสร้างที่มีพลังงานต่ำสุด ดังนั้นจึงต้องมีการสุ่มตัวอย่างที่ระมัดระวังและมีการลู่เข้าที่ดี รวมถึงตัวประมาณค่าเฉพาะทาง