ऊर्जा अपवाह

यांत्रिक प्रणालियों के कंप्यूटर सिमुलेशन में समय के साथ बंद प्रणाली की कुल ऊर्जा में क्रमिक परिवर्तन ऊर्जा बहाव के रूप में है। यांत्रिकी के नियमों के अनुसार ऊर्जा गतिमान स्थिर रूप में बनी रहनी चाहिए और उसे परिवर्तित नहीं होनी चाहिए। चूंकि, सिमुलेशन में ऊर्जा कम समय के पैमाने पर उतार-चढ़ाव कर सकती है और इस प्रकार संख्यात्मक साधारण अंतर समीकरण की कलाकृतियों के कारण बहुत लंबे समय के पैमाने पर बढ़ या घट सकती है, जो परिमित समय Δ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.