ऊर्जा अपवाह

यांत्रिक प्रणालियों के कंप्यूटर सिमुलेशन में समय के साथ बंद प्रणाली की कुल ऊर्जा में क्रमिक परिवर्तन ऊर्जा बहाव के रूप में है। यांत्रिकी के नियमों के अनुसार ऊर्जा गतिमान स्थिर रूप में बनी रहनी चाहिए और उसे परिवर्तित नहीं होनी चाहिए। चूंकि, सिमुलेशन में ऊर्जा कम समय के पैमाने पर उतार-चढ़ाव कर सकती है और इस प्रकार संख्यात्मक साधारण अंतर समीकरण की कलाकृतियों के कारण बहुत लंबे समय के पैमाने पर बढ़ या घट सकती है, जो परिमित समय Δt चरण के उपयोग के साथ उत्पन्न होती है। यह कुछ सीमा तक फ्लाइंग आइस क्यूब समस्या के समान होती है, जिसके द्वारा ऊर्जा के समविभाजन पर नियंत्रण में संख्यात्मक त्रुटियां कंपन ऊर्जा में बदल सकती हैं।

और इस प्रकार विशेष रूप से ऊर्जा में तेजी से वृद्धि होने की प्रवृत्ति को अंतःबोध के द्वारा सहज रूप से समझा जा सकता है क्योंकि प्रत्येक चरण में एक छोटे से पर्टर्बेशन δv को वास्तविक वेग vtrue के रूप में प्रस्तुत किया जाता है, जो v के साथ असंबंधित है और जो सरल समाकलन विधियों के लिए सही रूप में होता है और इस प्रकार ऊर्जा में द्वितीय क्रम में वृद्धि होती है।


 * $$E = \sum m \mathbf{v}^{2} = \sum m \mathbf{v}_\mathrm{true}^{2} + \sum m \ \delta \mathbf{v}^{2}$$

क्रॉस टर्म में v · δv शून्य रूप में होते है, क्योंकि इनमे कोई संबंध नहीं है।

ऊर्जा बहाव सामान्यतः अवमंदन - संख्यात्मक समाकलन योजनाओं के लिए पर्याप्त है जो सहानुभूतिपूर्ण इंटीग्रेटर नहीं हैं, जैसे कि रनगे-कुट्टा परिवार। आमतौर पर आणविक गतिकी में उपयोग किए जाने वाले सहानुभूतिपूर्ण इंटीग्रेटर्स, जैसे कि वेरलेट  समाकलन परिवार, बहुत लंबे समय के पैमाने पर ऊर्जा में वृद्धि प्रदर्शित करते हैं, हालांकि त्रुटि लगभग स्थिर रहती है। ये इंटीग्रेटर्स वास्तव में सिस्टम के हैमिल्टनियन यांत्रिकी को पुन: उत्पन्न नहीं करते हैं; इसके बजाय, वे बारीकी से संबंधित छाया हैमिल्टनियन को पुन: उत्पन्न करते हैं जिसका मूल्य वे परिमाण के कई आदेशों को अधिक बारीकी से संरक्षित करते हैं। सच्चे हैमिल्टनियन के लिए ऊर्जा संरक्षण की सटीकता समय कदम पर निर्भर है। रेफ नाम = गन्स > एक सहानुभूतिपूर्ण इंटीग्रेटर के संशोधित हैमिल्टनियन से गणना की गई ऊर्जा है $$\mathcal{O}\left(\Delta t^{p}\right)$$ सच्चे हैमिल्टनियन से।

एनर्जी ड्रिफ्ट पैरामीट्रिक ऑसिलेटर # पैरामीट्रिक रेजोनेंस के समान है, जिसमें एक परिमित, असतत टाइमस्टेपिंग योजना के परिणामस्वरूप वेग अपडेट की आवृत्ति के करीब आवृत्ति के साथ गति के गैर-भौतिक, सीमित नमूने होंगे। इस प्रकार अधिकतम चरण आकार पर प्रतिबंध जो किसी दिए गए सिस्टम के लिए स्थिर होगा, सिस्टम की गति के सबसे तेज़ मौलिक मोड की अवधि के समानुपाती होता है। एक प्राकृतिक आवृत्ति ω के साथ एक गति के लिए, कृत्रिम अनुनाद पेश किए जाते हैं जब वेग की आवृत्ति अद्यतन होती है, $$\frac{2\pi}{\Delta t}$$ ω से संबंधित है
 * $$\frac{n}{m}\omega = \frac{2\pi}{\Delta t}$$

जहाँ n और m अनुनाद क्रम का वर्णन करने वाले पूर्णांक हैं। वेरलेट समाकलन के लिए, चौथे क्रम तक अनुनाद $$\left(\frac{n}{m} = 4\right)$$ अक्सर संख्यात्मक अस्थिरता का कारण बनता है, जिससे टाइमस्टेप के आकार पर प्रतिबंध लग जाता है
 * $$\Delta t < \frac{\sqrt{2}}{\omega} \approx 0.225p$$

जहां ω प्रणाली में सबसे तेज गति की आवृत्ति है और पी इसकी अवधि है। अधिकांश जैव-आण्विक प्रणालियों में सबसे तेज़ गति में हाइड्रोजन परमाणुओं की गति शामिल होती है; इस प्रकार हाइड्रोजन गति को प्रतिबंधित करने के लिए बाधा एल्गोरिदम का उपयोग करना आम है और इस प्रकार सिमुलेशन में उपयोग किए जा सकने वाले अधिकतम स्थिर समय कदम को बढ़ाता है। हालांकि, क्योंकि भारी-परमाणु गतियों के समय के पैमाने हाइड्रोजन गतियों से व्यापक रूप से भिन्न नहीं होते हैं, व्यवहार में यह समय चरण में केवल दो गुना वृद्धि की अनुमति देता है। ऑल-एटम बायोमोलेक्युलर सिमुलेशन में सामान्य अभ्यास अनियंत्रित सिमुलेशन के लिए 1 गुजरने  (एफएस) के समय चरण का उपयोग करना है और प्रतिबंधित सिमुलेशन के लिए 2 एफएस है, हालांकि कुछ प्रणालियों या पैरामीटर के विकल्पों के लिए बड़े समय के चरण संभव हो सकते हैं।

कार्य (गणित) के मूल्यांकन में खामियों के कारण ऊर्जा बहाव भी हो सकता है, आमतौर पर सिमुलेशन मापदंडों के कारण जो कम्प्यूटेशनल गति के लिए सटीकता का त्याग करते हैं। उदाहरण के लिए, इलेक्ट्रोस्टैटिक बलों के मूल्यांकन के लिए कटऑफ योजनाएं हर बार कदम के साथ ऊर्जा में व्यवस्थित त्रुटियों का परिचय देती हैं, क्योंकि यदि पर्याप्त चिकनाई का उपयोग नहीं किया जाता है तो कण कटऑफ त्रिज्या में आगे और पीछे चलते हैं। कण जाल इवाल्ड  समेशन इस प्रभाव का एक समाधान है, लेकिन अपनी खुद की कलाकृतियों का परिचय देता है। सिम्युलेटेड की जा रही प्रणाली में त्रुटियां विस्फोटक के रूप में विशेषता वाले ऊर्जा बहाव को भी प्रेरित कर सकती हैं जो कलाकृतियां नहीं हैं, लेकिन प्रारंभिक स्थितियों की अस्थिरता को दर्शाती हैं; यह तब हो सकता है जब उत्पादन गतिशीलता शुरू करने से पहले सिस्टम को पर्याप्त संरचनात्मक न्यूनीकरण के अधीन नहीं किया गया हो। व्यवहार में, ऊर्जा बहाव को समय के साथ प्रतिशत वृद्धि के रूप में मापा जा सकता है, या सिस्टम में दी गई मात्रा में ऊर्जा जोड़ने के लिए आवश्यक समय के रूप में।

ऊर्जा बहाव के व्यावहारिक प्रभाव सिमुलेशन स्थितियों पर निर्भर करते हैं, थर्मोडायनामिक पहनावा सिम्युलेटेड किया जा रहा है, और अध्ययन के तहत सिमुलेशन का इरादा उपयोग; उदाहरण के लिए, जहां तापमान स्थिर रखा जाता है, वहां कैनोनिकल पहनावा की तुलना में माइक्रोकैनोनिकल पहनावा के सिमुलेशन के लिए ऊर्जा बहाव के बहुत अधिक गंभीर परिणाम होते हैं। हालांकि, यह दिखाया गया है कि लंबे समय तक माइक्रोकैनोनिकल पहनावा सिमुलेशन महत्वहीन ऊर्जा बहाव के साथ किया जा सकता है, जिसमें लचीले अणु शामिल हैं जो बाधाओं और इवाल्ड योगों को शामिल करते हैं। ऊर्जा बहाव को अक्सर सिमुलेशन की गुणवत्ता के एक उपाय के रूप में उपयोग किया जाता है, और प्रोटीन डाटा बैंक के अनुरूप आणविक गतिशीलता प्रक्षेपवक्र डेटा के एक बड़े भंडार में नियमित रूप से रिपोर्ट किए जाने के लिए एक गुणवत्ता मीट्रिक के रूप में प्रस्तावित किया गया है।

अग्रिम पठन

 * Sanz-Serna JM, Calvo MP. (1994). Numerical Hamiltonian Problems. Chapman & Hall, London, England.