พีชคณิตเชิงเส้นเชิงตัวเลขและปัญหาค่าลักษณะเฉพาะในฟิสิกส์

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

Definition

พีชคณิตเชิงเส้นเชิงตัวเลขในฟิสิกส์คือชุดของอัลกอริทึมสำหรับการแก้สมการเมทริกซ์และปัญหาค่าลักษณะเฉพาะที่เกิดขึ้นเมื่อตัวดำเนินการทางฟิสิกส์แบบต่อเนื่องถูกแทนด้วยฐานจำกัดหรือบนกริด

Scope

หัวข้อนี้ครอบคลุมการคำนวณเมทริกซ์ที่เป็นหัวใจสำคัญของฟิสิกส์: การแก้ระบบเชิงเส้นด้วยวิธีตรงและวิธีวนซ้ำ และการคำนวณค่าลักษณะเฉพาะและเวกเตอร์ลักษณะเฉพาะของเมทริกซ์เฮอร์มิเชียนขนาดใหญ่ ซึ่งมักจะเป็นเมทริกซ์แบบสปาร์ส ผ่านอัลกอริทึม QR, Jacobi, Lanczos และ conjugate-gradient โดยเน้นโครงสร้างของเมทริกซ์ทางฟิสิกส์ เช่น ความเป็นสปาร์สและความเป็นเฮอร์มิเชียน

Core questions

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

Key theories

ตัวแก้ระบบเชิงเส้นแบบตรงและแบบวนซ้ำ

ระบบเชิงเส้นถูกแก้โดยการแยกตัวประกอบโดยตรง เช่น LU และ Cholesky ซึ่งแม่นยำจนถึงค่าปัดเศษ หรือโดยวิธีการวนซ้ำแบบ Krylov เช่น conjugate gradients ที่ใช้ประโยชน์จากความเป็นสปาร์สและลู่เข้าสู่ค่าความคลาดเคลื่อนที่ยอมรับได้

อัลกอริทึมค่าลักษณะเฉพาะ

ค่าลักษณะเฉพาะและเวกเตอร์ลักษณะเฉพาะถูกคำนวณโดยอัลกอริทึม QR และการหมุน Jacobi สำหรับเมทริกซ์หนาแน่น ซึ่งให้สเปกตรัมแบบไม่ต่อเนื่องของตัวดำเนินการทางฟิสิกส์ที่แสดงในฐานจำกัด

วิธีการ Lanczos และปริภูมิ Krylov

อัลกอริทึม Lanczos สร้างการฉายภาพแบบ tridiagonal ขนาดเล็กของเมทริกซ์เฮอร์มิเชียนแบบสปาร์สขนาดใหญ่ในปริภูมิ Krylov ซึ่งช่วยให้สามารถค้นหาค่าลักษณะเฉพาะและเวกเตอร์ลักษณะเฉพาะที่ปลายสุดบางค่าได้โดยไม่ต้องจัดเก็บเมทริกซ์ทั้งหมด

Clinical relevance

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

History

การคำนวณค่าลักษณะเฉพาะของเมทริกซ์ในทางปฏิบัติพัฒนาขึ้นในช่วงกลางศตวรรษที่ 20 ด้วยการวนซ้ำของ Lanczos ในปี 1950 และอัลกอริทึม QR ในช่วงต้นทศวรรษ 1960 การเพิ่มขึ้นของปัญหาแบบสปาร์สขนาดใหญ่ในฟิสิกส์ทำให้วิธีการปริภูมิ Krylov กลายเป็นเครื่องมือหลักสำหรับการวิเคราะห์สเปกตรัมของแฮมิลโทเนียนที่มีมิติสูง

Key figures

• Cornelius Lanczos
• Gene H. Golub
• James H. Wilkinson

Seminal works

• golub2013
• lanczos1950

Frequently asked questions

เหตุใดจึงใช้วิธีการวนซ้ำแทนที่จะทำการทำให้เมทริกซ์ทั้งหมดเป็นแนวทแยง?

แฮมิลโทเนียนทางฟิสิกส์อาจมีมิติเป็นพันล้านแต่เป็นแบบสปาร์ส ดังนั้นการจัดเก็บหรือแยกตัวประกอบทั้งหมดจึงเป็นไปไม่ได้ วิธีการวนซ้ำแบบ Krylov เช่น Lanczos ต้องการเพียงการกระทำของเมทริกซ์บนเวกเตอร์และสามารถสกัดสถานะพื้นฐานที่ต่ำที่สุดไม่กี่สถานะที่ฟิสิกส์มักจะสนใจได้

เหตุใดความเป็นเฮอร์มิเชียนของเมทริกซ์ทางฟิสิกส์จึงมีความสำคัญเชิงตัวเลข?

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