मैं परीक्षण प्रयोजनों के लिए पानी के आणविक गतिशीलता सिमुलेशन चला रहा हूं। यदि आप एक डीएफटी आदमी से पूछते हैं, तो बॉक्स काफी छोटा है, अगर आप शास्त्रीय एमडी से चलने वाले व्यक्ति और अपेक्षाकृत बड़े से पूछते हैं, तो मेरे पास आवधिक सीमा की स्थिति में 58 पानी के अणु हैं।
CPU समय बचाने के लिए, मैं ab inioio MD चलाने से पहले अपने सेल को शास्त्रीय बल क्षेत्र के साथ अनुकूलित कर रहा हूं। मैं सिस्टम को क्लासिक रूप से 1K के लिए 300K पर संतुलित करता हूं, फिर अंतिम स्नैपशॉट लें और इसे ab initio MD के लिए इनपुट के रूप में उपयोग करें। मेरा एब इनिटियो एमडी नियमित रूप से डीएफटी-आधारित बोर्न-ओपेनहाइमर एमडी है जिसमें एक विमान तरंग आधार सेट और पीएडब्ल्यू (छद्म) क्षमता (वीएएसपी कोड है)। शास्त्रीय और ab initio दोनों प्रकार के सिमुलेशन में मैं तापमान-वेग को 300K पर वेग-रिस्क्लोज़िंग थर्मोस्टेट का उपयोग करके रख रहा हूँ।
मैं शास्त्रीय और ab initio के बीच संक्रमण करने के लिए दो अलग-अलग तरीकों का सर्वेक्षण कर रहा हूं:
- शास्त्रीय प्रक्षेपवक्र से प्रारंभिक वेगों और पदों को लें और उन्हें अब इनिटियो सिमुलेशन के लिए प्रारंभिक कॉन्फ़िगरेशन के रूप में आयात करें
- शास्त्रीय स्थितियों को ध्यान में रखते हुए तापमान को शून्य करने के लिए सिस्टम को फ्रीज करें, उसे डीएफटी कोड में आयात करें, फिर जल्दी से (मैं इसे फिलहाल 0.5 पीएस में कर रहा हूं) 300K तक गरम करें
मैं उम्मीद कर रहा था कि दोनों रणनीतियाँ संक्षिप्त (10 पीएस) संतुलन अवधि के बाद समान औसत ऊर्जा का नेतृत्व करेंगी, विशेष रूप से यह देखते हुए कि शुरुआती कॉन्फ़िगरेशन बिल्कुल समान है (समान प्रारंभिक स्थिति) उल्लेखित तापमान चाल को छोड़कर (प्रारंभिक वेग भिन्न होते हैं) । यह मामला नहीं है। नीचे दिए गए आंकड़े से पता चलता है कि सिमुलेशन जहां सिस्टम जमे हुए है और फिर जल्दी से गर्म हो जाता है, एक ऊर्जा क्षेत्र को दूसरे ईवी के मुकाबले ऊर्जा में लगभग 1 ईवी कम पाता है, जहां वेग जहां शास्त्रीय एमडी से आयात किया जाता है।
मेरे प्रश्न हैं:
- क्या यह उम्मीद की जानी है;
- शास्त्रीय से अब तक initio MD में संक्रमण को अनुकूलित करने के लिए सफल रणनीति ज्ञात हैं;
- और क्या आप मुझे इस मामले पर प्रासंगिक साहित्य की ओर इशारा कर सकते हैं?
संपादित करें:
मैं इस समय कुछ और परीक्षण चला रहा हूं और मेरे पास मौजूद सीमित डेटा के बावजूद- ऐसा लगता है कि यह एक सिस्टम-विशिष्ट समस्या हो सकती है। एक ही आकार के एक बॉक्स में पानी के बजाय मेथनॉल के साथ एक परीक्षण से पता चला है कि दो अलग-अलग प्रारंभिक वेग योजनाएं एक ही औसत ऊर्जा में जल्दी से जुट जाती हैं। हालांकि, मेथनॉल के मामले में शास्त्रीय विन्यास क्वांटम एक के बहुत करीब था, अर्थात, टी = 0 पर ऊर्जा अभिसरण के बाद औसत ऊर्जा के बहुत करीब थी। पानी एक कुख्यात कठिन प्रणाली है, इसलिए शायद यह समस्या कमोबेश पानी की है। यदि कोई उत्तर नहीं जोड़ा जाता है, तो मैं सभी परिणामों के साथ एक बार अपने परिणामों के आधार पर कोशिश करूंगा और पोस्ट करूंगा।