आणविक यांत्रिकी

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

सभी-परमाणु आणविक यांत्रिकी विधियों में निम्नलिखित गुण हैं:
 * प्रत्येक परमाणु को एक कण के रूप में अनुकरण किया जाता है
 * प्रत्येक कण को ​​एक अर्धव्यास (सामान्य तौर पर वैन डेर वाल्स त्रिज्या), ध्रुवीकरण, और एक स्थिर शुद्ध आवेश (सामान्य तौर पर क्वांटम गणना) सौंपा गया है।
 * जुडात्मक प्रभाव को प्रयोगात्मक या गणना की गई जुडात्मक लंबाई के बराबर संतुलन दूरी के साथ स्प्रिंग्स के रूप में माना जाता है।

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

कार्यात्मक रूप
निम्नलिखित कार्यात्मक सारांश, जिसे रसायन विज्ञान में संभावित अंतर-परमाणु कार्य या बल क्षेत्र (रसायन विज्ञान) कहा जाता है, आणविक प्रणाली की संभावित ऊर्जा (E) की व्यक्तिगत ऊर्जा के योग अंतर्गत दी गई रचना में गणना करता है।

$$\ E = E_\text{covalent} + E_\text{noncovalent} \, $$

जहां सहसंयोजक और गैर-सहसंयोजक योगदान के घटक निम्नलिखित योगों द्वारा दिए गए हैं:

$$\ E_\text{covalent} = E_\text{bond} + E_\text{angle} + E_\text{dihedral}$$

$$\ E_\text{noncovalent} = E_\text{electrostatic} + E_\text{van der Waals} $$

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

पूर्ण रूप से गणना करने के लिए गैर-बंधित शर्तें बहुत अधिक कम्प्यूटेशनल रूप से महंगी हैं, क्योंकि एक विशिष्ट परमाणु अपने कुछ पड़ोसियों से ही जुड़ा होता है, लेकिन अणु में हर दूसरे परमाणु के साथ संपर्क करता है। सौभाग्य से वैन डेर वाल्स बल की अवधि तेजी से गिरती है। इसे आम तौर पर 6-12 लेनार्ड-जोन्स क्षमता का उपयोग करके तैयार किया जाता है, जिसका अर्थ है कि आकर्षक बल आर के रूप में दूरी के साथ गिर जाते हैं।−6 और प्रतिकारक बल r के रूप में−12, जहाँ r दो परमाणुओं के बीच की दूरी को दर्शाता है। प्रतिकारक भाग आर−12 हालांकि अभौतिक है, क्योंकि प्रतिकर्षण तेजी से बढ़ता है। लेनार्ड-जोन्स 6–12 क्षमता द्वारा वैन डेर वाल्स बलों का विवरण अशुद्धि का परिचय देता है, जो कम दूरी पर महत्वपूर्ण हो जाता है। आम तौर पर एक कटऑफ त्रिज्या का उपयोग गणना को गति देने के लिए किया जाता है ताकि परमाणु जोड़े जो कि कटऑफ से अधिक दूरी पर हैं, शून्य की वैन डेर वाल्स इंटरैक्शन ऊर्जा हो।

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

प्रत्येक ऊर्जा अवधि के कार्यात्मक रूप के अलावा, एक उपयोगी ऊर्जा फ़ंक्शन को बल स्थिरांक, वैन डेर वाल्स मल्टीप्लायरों और अन्य स्थिर शर्तों के लिए पैरामीटर निर्दिष्ट किया जाना चाहिए। ये शब्द, संतुलन बंधन, कोण, और डायहेड्रल मान, आंशिक आवेश मान, परमाणु द्रव्यमान और त्रिज्या, और ऊर्जा कार्य परिभाषाओं के साथ सामूहिक रूप से एक बल क्षेत्र (रसायन विज्ञान) कहलाते हैं। प्राचलीकरण आमतौर पर प्रयोगात्मक मूल्यों और सैद्धांतिक गणना परिणामों के साथ समझौते के माध्यम से किया जाता है। पिछले MM4 संस्करण में नॉर्मन एल. एलींगर के बल क्षेत्र की गणना 0.35 kcal/mol की rms त्रुटि के साथ हाइड्रोकार्बन के गठन की ऊष्मा, 24 सेमी की rms त्रुटि के साथ कंपन स्पेक्ट्रा के लिए की जाती है-1, 2.2 की rms त्रुटि के साथ घूर्णी अवरोध°, सी-सी बांड लंबाई 0.004 ए के भीतर और सी-सी-सी कोण 1 के भीतर°. बाद में MM4 संस्करणों में एलिफैटिक एमाइन जैसे हेटेरोएटम्स के साथ यौगिक भी शामिल हैं। प्रत्येक बल क्षेत्र को आंतरिक रूप से सुसंगत होने के लिए पैरामीटरकृत किया जाता है, लेकिन पैरामीटर आमतौर पर एक बल क्षेत्र से दूसरे में स्थानांतरित नहीं होते हैं।

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

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

बाध्यकारी स्थिरांक की गणना के लिए आणविक यांत्रिकी संभावित ऊर्जा कार्यों का उपयोग किया गया है,    प्रोटीन तह कैनेटीक्स, प्रोटोनेशन संतुलन, डॉकिंग (आणविक), और प्रोटीन डिजाइन के लिए।

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

सॉफ्टवेयर पैकेज
मुख्य लेख: आणविक यांत्रिकी मॉडलिंग के लिए सॉफ्टवेयर की तुलना

यह एक सीमित सूची है; कई और पैकेज उपलब्ध हैं।


 * ऐबालोन
 * एसीईएमडी - जीपीयू एमडी [14]
 * एम्बर
 * एस्कलाफ डिजाइनर [15]
 * मालिक
 * आकर्षण
 * कॉसमॉस [16]
 * CP2K
 * घेमिकल
 * GROMACS
 * ग्रोमोस
 * आंतरिक समन्वय यांत्रिकी (आईसीएम)
 * लैम्प्स
 * मैक्रोमॉडल
 * MDynaMix
 * आणविक परिचालन पर्यावरण (एमओई)
 * NAMD
 * क्यू
 * क्यू केम
 * परहेज़गार
 * स्ट्रुएमएम3डी (STR3DI32) [17]
 * टिन से मढ़नेवाला
 * एक्स-प्लोर
 * यासरा
 * राशि चक्र [

यह भी देखें

 * आणविक ग्राफिक्स
 * आणविक गतिकी
 * अणु संपादक
 * बल क्षेत्र (रसायन विज्ञान)
 * बल क्षेत्र कार्यान्वयन की तुलना
 * आणविक डिजाइन सॉफ्टवेयर
 * GPU पर आणविक मॉडलिंग
 * आणविक यांत्रिकी मॉडलिंग के लिए सॉफ्टवेयर की तुलना
 * मोंटे कार्लो आण्विक मॉडलिंग के लिए सॉफ्टवेयर की सूची

संदर्भ

 * 1) ज़गरबोवा एम, एट अल। (2010)। "जोड़ीवार-योगात्मक अनुभवजन्य बल क्षेत्रों में त्रुटियों का बड़े पैमाने पर मुआवजा: कठोर डीएफटी-एसएपीटी गणनाओं के साथ एम्बर इंटरमॉलिक्यूलर शर्तों की तुलना"। भौतिक। रसायन। रसायन। भौतिकी । 12 (35): 10476-10493। बिबकोड : 2010PCCP...1210476Z । डीओआई : 10.1039/सी002656ई . पीएमआईडी  20603660 ।
 * 2) एलींगर, एनएल; चेन, के.; ली, जे.-एच। जे कम्प्यूट। रसायन। 1996, 17 , 642 3बी2-यू
 * 3) कुओ-सियांग चेन, जेन-हुई ली, यी फैन, नॉर्मन एल. एलिंगर जे. कंप्यूट। रसायन। 2007, 28, 2391 https://onlinelibrary.wiley.com/doi/full/10.1002/jcc.20737
 * 4) कुह्न बी, कोलमैन पीए (अक्टूबर 2000)। "एविडिन और स्ट्रेप्टाविडिन के लिए लिगैंड्स के एक विविध सेट का बंधन: आणविक यांत्रिकी और निरंतर विलायक मॉडल के संयोजन द्वारा उनके सापेक्ष समानता की सटीक मात्रात्मक भविष्यवाणी"। जर्नल ऑफ़ मेडिसिनल केमिस्ट्री । 43 (20): 3786–91। डीओआई : 10.1021/jm000241h . पीएमआईडी  11020294 ।
 * 5) हुओ एस, मासोवा I, कोलमैन पीए (जनवरी 2002)। "1:1 मानव विकास हार्मोन-रिसेप्टर कॉम्प्लेक्स की कम्प्यूटेशनल एलेनिन स्कैनिंग"। जे कंप्यूट केम । 23 (1): 15–27। डीओआई : 10.1002/जेसीसी.1153 । पीएमआईडी  11913381 । एस 2 सीआईडी ​​10381457 .
 * 6) मोब्ले डीएल, ग्रेव्स एपी, चोदेरा जेडी, मैकरेनॉल्ड्स एसी, शोइचेत बीके, दिल केए (अगस्त 2007)। "एक साधारण मॉडल साइट के लिए पूर्ण लिगैंड बाइंडिंग फ्री एनर्जी की भविष्यवाणी" । जे मोल बायोल । 371 (4): 1118–34। डीओआई : 10.1016/जे.जेएमबी.2007.06.002 । पीएमसी  2104542 । पीएमआईडी  17599350 ।
 * 7) वांग जे, कांग एक्स, कुंट्ज़ आईडी, कोलमैन पीए (अप्रैल 2005)। "फार्माकोफोर मॉडल, कठोर डॉकिंग, सॉल्वेशन डॉकिंग और एमएम-पीबी/एसए का उपयोग करके एचआईवी-1 रिवर्स ट्रांसक्रिपटेस के लिए पदानुक्रमित डेटाबेस स्क्रीनिंग"। जर्नल ऑफ़ मेडिसिनल केमिस्ट्री । 48 (7): 2432-44। डीओआई : 10.1021/jm049606e . पीएमआईडी  15801834 ।
 * 8) कोलमैन पीए, मासोवा I, रेयेस सी, एट अल। (दिसंबर 2000)। "संरचनाओं की गणना और जटिल अणुओं की मुक्त ऊर्जा: आणविक यांत्रिकी और सातत्य मॉडल का संयोजन"। एसीसी केम रेस । 33 (12): 889–97। साइटसेरएक्स  10.1.1.469.844 । डीओआई : 10.1021/ar000033j । पीएमआईडी  11123888 ।
 * 9) स्नो सीडी, गुयेन एच, पांडे वीएस, ग्रुबेले एम (नवंबर 2002)। "सिम्युलेटेड और प्रायोगिक प्रोटीन-फोल्डिंग डायनेमिक्स की पूर्ण तुलना"। प्रकृति । 420 (6911): 102–6। बिबकोड : 2002Natur.420..102S . डीओई : 10.1038/प्रकृति 01160 . पीएमआईडी 12422224 . एस 2 सीआईडी ​​1061159.
 * 10) बर्थ पी, अल्बर्ट टी, हारबरी पीबी (मार्च 2007)। "प्रोटीन आयनीकरण स्थिरांक पर विलायक प्रभावों की सटीक, संरचना-निर्भर भविष्यवाणियां" । प्रोक नेटल एकेड साइंस यूएसए । 104 (12): 4898–903। बिबकोड : 2007पीएनएएस..104.4898बी । डीओआई : 10.1073/पीएनएएस.0700188104 । पीएमसी  1829236 । पीएमआईडी  17360348.
 * 11) चक्रवर्ती आर, क्लिबानोव एएम, फ्रेज़नर आरए (जुलाई 2005)। "देशी प्रोटीन लिगैंड-बाइंडिंग और एंजाइम सक्रिय साइट अनुक्रमों की कम्प्यूटेशनल भविष्यवाणी" । प्रोक नेटल एकेड साइंस यूएसए । 102 (29): 10153–8। बिबकोड : 2005PNAS..10210153C । डीओआई : 10.1073/पीएनएएस.0504023102 । पीएमसी  1177389 । पीएमआईडी  15998733.
 * 12) बोस एफई, हारबरी पीबी (जुलाई 2008)। "आणविक-यांत्रिकी ऊर्जा मॉडल के आधार पर प्रोटीन-लिगैंड बाइंडिंग का डिज़ाइन" । जे मोल बायोल । 380 (2): 415-24। डीओआई : 10.1016/जे.जेएमबी.2008.04.001 । पीएमसी  2569001 । पीएमआईडी  18514737 ।
 * 13) क्रैमर, क्रिस्टोफर जे। (2004)। कम्प्यूटेशनल रसायन विज्ञान की अनिवार्यता: सिद्धांत और मॉडल (दूसरा संस्करण)। चिचेस्टर, वेस्ट ससेक्स, इंग्लैंड: विली। आईएसबीएन 0-470-09182-7. ओसीएलसी  55887497 ।
 * 14) एसीईएमडी - जीपीयू एमडी
 * 15) एस्कलाफ
 * 16) कॉसमॉस
 * 17) स्ट्रूएमएम3डी (STR3DI32)
 * 18) वेबैक मशीन पर राशि संग्रह 2009-12-16
 * एलींगर एनएल, बुर्कर्ट यू (1982)। आणविक यांत्रिकी । एक अमेरिकन केमिकल सोसायटी प्रकाशन। आईएसबीएन 978-0-8412-0885-8.
 * बॉक्स वीजी (मार्च 1997)। "मात्राबद्ध वैलेंस बांड के आणविक यांत्रिकी"। जे मोल मॉडल । 3 (3): 124–41। डीओई : 10.1007/एस 008940050026 . एस 2 सीआईडी ​​93821090.
 * बॉक्स वीजी (12 नवंबर 1998)। "मोनोसेकेराइड और उनके डेरिवेटिव का विसंगतिपूर्ण प्रभाव। नए क्यूवीबीएमएम आण्विक यांत्रिकी बल क्षेत्र से अंतर्दृष्टि" । हेट्रोसायकल । 48 (11): 2389–417। डीओआई : 10.3987/आरईवी-98-504 . 31 मई 2012 को मूल से संग्रहीत । 12 नवंबर 2008 को पुनःप्राप्त ।
 * बॉक्स वीजी (2004)। "पोलिन्यूक्लियोटाइड्स और उनके दोहरे हेलिकॉप्टरों में स्टीरियो-इलेक्ट्रॉनिक प्रभाव"। जे मोल। संरचना । 689 (1–2): 33–41। बिबकोड : 2004JMoSt.689...33B . डीओआई : 10.1016/जे.मोलस्ट्रक.2003.10.019 ।
 * बेकर ओम (2001)। कम्प्यूटेशनल बायोकैमिस्ट्री और बायोफिजिक्स । न्यूयॉर्क, एनवाई: मार्सेल डेकर। आईएसबीएन 978-0-8247-0455-1.
 * मैकेरल एडी (अक्टूबर 2004)। "जैविक मैक्रोमोलेक्यूल्स के लिए अनुभवजन्य बल क्षेत्र: सिंहावलोकन और मुद्दे"। जे कंप्यूट केम । 25 (13): 1584–604। डीओआई : 10.1002/जेसीसी.20082 । पीएमआईडी  15264253 । एस 2 सीआईडी ​​9162620 .
 * श्लिक टी (2002)। आणविक मॉडलिंग और सिमुलेशन: एक अंतःविषय गाइड । बर्लिन: स्प्रिंगर. आईएसबीएन 978-0-387-95404-2.
 * कृष्णन नम्बूरी; रामचंद्रन, के.एस.; दीपा गोपाकुमार (2008)। कम्प्यूटेशनल रसायन विज्ञान और आणविक मॉडलिंग: सिद्धांत और अनुप्रयोग । बर्लिन: स्प्रिंगर. आईएसबीएन 978-3-540-77302-3.

इस पेज में लापता आंतरिक लिंक की सूची

 * मोलेकुलर
 * अंतर-परमाणु क्षमता
 * प्रोटीन की संभावित ऊर्जा
 * गाऊसी (सॉफ्टवेयर)
 * ढतला हुआ वंश
 * निहित समाधान
 * पानी का मॉडल

बाहरी संबंध

 * आण्विक गतिशीलता सिमुलेशन विधियों को संशोधित किया गया
 * आणविक यांत्रिकी - यह सरल है