ऊर्जा अपवाह

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

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

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

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

अग्रिम पठन

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