व्यवरोध (कम्प्यूटेशनल रसायन विज्ञान)

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

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

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

गणितीय पृष्ठभूमि
N कणों के एक सेट की गति को दूसरे क्रम के साधारण अंतर समीकरणों, न्यूटन के दूसरे नियम के एक सेट द्वारा वर्णित किया जा सकता है, जिसे आव्यूह रूप में लिखा जा सकता है



\mathbf{M} \cdot \frac{d^{2}\mathbf{q}}{dt^{2}} = \mathbf{f} = -\frac{\partial V}{\partial \mathbf{q}} $$ जहां M एक द्रव्यमान आव्यूह (मास मैट्रिक्स ) है और q सामान्यीकृत निर्देशांक का सदिश (ज्यामितीय) है जो कणों की स्थिति का वर्णन करता है। उदाहरण के लिए, सदिश q कण स्थितियों rk का 3N कार्टेशियन निर्देशांक हो सकता है, जहां k 1 से N तक चलता है; व्यवरोधओं की अनुपस्थिति में, 'M ' कण द्रव्यमान का 3Nx3N विकर्ण वर्ग आव्यूह होगा। सदिश 'f' सामान्यीकृत बलों का प्रतिनिधित्व करता है और अदिश V('q ') संभावित ऊर्जा का प्रतिनिधित्व करता है, ये दोनों सामान्यीकृत निर्देशांक 'q ' के कार्य हैं।

यदि M व्यवरोधएं उपलब्ध हैं, तो निर्देशांक को M समय-स्वतंत्र बीजगणितीय समीकरणों को भी संतुष्ट करना होगा



g_{j}(\mathbf{q}) = 0 $$ जहां सूचकांक j 1 से M तक चलता है। संक्षिप्तता के लिए, ये फलन gi हैं नीचे M-आयामी सदिश 'g ' में समूहीकृत किया गया है। कार्य न्यूटन के दूसरे नियम के सामान्य अंतर समीकरणों (ओडीई) के बजाय अंतर-बीजगणितीय (डीएई) समीकरणों के संयुक्त सेट को हल करना है।

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

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

तीसरा दृष्टिकोण व्यवरोधओं को पूरा करने के लिए आवश्यक समन्वय समायोजन निर्धारित करने के लिए लैग्रेंज मल्टीप्लायरों या व्यवरोध मैनिफोल्ड के प्रक्षेपण जैसी विधि का उपयोग करना है।

अंत में, विभिन्न संकर दृष्टिकोण हैं जिनमें व्यवरोधओं के विभिन्न सेटों को विभिन्न तरीकों से संतुष्ट किया जाता है, उदाहरण के लिए, आंतरिक निर्देशांक, स्पष्ट बल और अंतर्निहित-बल समाधान।

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

आंतरिक निर्देशांक में कुशल पुनरावर्ती ऊर्जा न्यूनीकरण के लिए मूल तरीके Gō और सहकर्मियों द्वारा विकसित किए गए थे।

कुशल पुनरावर्ती, आंतरिक-समन्वय व्यवरोध सॉल्वर को आण्विक गतिकी तक बढ़ाया गया था। एनालॉग पद्धतियां बाद में अन्य प्रणालियों में लागू की गईं।

लैग्रेंज गुणक-आधारित विधियाँ
व्यवरोध एल्गोरिदम का उपयोग करने वाले अधिकांश आण्विक गतिकी सिमुलेशन में, लैग्रेंज मल्टीप्लायरों की विधि का उपयोग करके व्यवरोधओं को लागू किया जाता है। समय t पर n रैखिक (होलोनोमिक व्यवरोधएं) व्यवरोधओं का एक सेट दिया गया है,


 * $$\sigma_k(t) := \| \mathbf x_{k\alpha}(t) - \mathbf x_{k\beta}(t) \|^2 - d_k^2 = 0, \quad k=1 \ldots n$$

जहाँ $$\scriptstyle \mathbf x_{k\alpha}(t)$$ और $$\scriptstyle\mathbf x_{k\beta}(t)$$ समय t और पर kth व्यवरोध में सम्मिलित दो कणों की स्थिति हैं $$d_k$$ निर्धारित अंतर-कण दूरी है।

इन व्यवरोधओं के कारण बलों को गति के समीकरणों में जोड़ा जाता है, जिसके परिणामस्वरूप, सिस्टम में प्रत्येक N कण के लिए


 * $$\frac{\partial^2 \mathbf x_i(t)}{\partial t^2} m_i = -\frac{\partial}{\partial \mathbf x_i} \left[ V(\mathbf x_i(t)) - \sum_{k=1}^n \lambda_k \sigma_k(t) \right], \quad i=1 \ldots N.$$

व्यवरोध बलों को जोड़ने से कुल ऊर्जा में परिवर्तन नहीं होता है, क्योंकि व्यवरोध बलों (कणों के समूह पर लिया गया जिन पर व्यवरोधएं कार्य करती हैं) द्वारा किया गया शुद्ध कार्य शून्य है। ध्यान दें कि साइन ऑन है $$\lambda_k$$ स्वच्छंद है और कुछ संदर्भ एक विपरीत चिन्ह है.

समय के संबंध में समीकरण के दोनों पक्षों को एकीकृत करने से, उस समय कणों के बाधित निर्देशांक, $$t + \Delta t$$, दिया जाता है,


 * $$\mathbf x_i(t + \Delta t) = \hat{\mathbf x}_i(t + \Delta t) + \sum_{k=1}^n \lambda_k \frac{\partial\sigma_k(t)}{\partial \mathbf x_i}\left(\Delta t\right)^2m_i^{-1}, \quad i=1 \ldots N$$

जहाँ $$\hat{\mathbf x}_i(t + \Delta t)$$ गति के अप्रतिबंधित समीकरणों को एकीकृत करने के बाद iवें कण की अप्रतिबंधित (या असंशोधित) स्थिति है।

व्यवरोधओं को पूरा करने के लिए $$\sigma_k(t + \Delta t)$$ अगले समय चरण में, लैग्रेंज गुणक को निम्नलिखित समीकरण के रूप में निर्धारित किया जाना चाहिए,


 * $$\sigma_k(t + \Delta t) := \left\| \mathbf x_{k\alpha}(t+\Delta t) - \mathbf x_{k\beta}(t+\Delta t)\right\|^2 - d_k^2 = 0.$$

इसका तात्पर्य एक प्रणाली को हल करना है $$n$$ गैर-रैखिक समीकरण


 * $$\sigma_j(t + \Delta t) := \left\| \hat{\mathbf x}_{j\alpha}(t+\Delta t) - \hat{\mathbf x}_{j\beta}(t+\Delta t) +   \sum_{k=1}^n \lambda_k \left(\Delta t\right)^2 \left[ \frac{\partial\sigma_k(t)}{\partial \mathbf x_{j\alpha}}m_{j\alpha}^{-1} -  \frac{\partial\sigma_k(t)}{\partial \mathbf x_{j\beta}}m_{j\beta}^{-1}\right] \right\|^2 - d_j^2 = 0, \quad j = 1 \ldots n$$

के लिए एक साथ $$n$$ अज्ञात लैग्रेंज गुणक $$\lambda_k$$.

की यह व्यवस्था $$n$$ गैर-रैखिक समीकरण $$n$$ अज्ञात को सामान्यतः न्यूटन की विधि| न्यूटन-रेफसन विधि का उपयोग करके हल किया जाता है जहां समाधान सदिश होता है $$\underline{\lambda}$$ का उपयोग कर अद्यतन किया जाता है


 * $$\underline{\lambda}^{(l+1)} \leftarrow \underline{\lambda}^{(l)} - \mathbf J_\sigma^{-1} \underline{\sigma}(t+\Delta t)$$

जहाँ $$\mathbf J_\sigma$$ जैकोबियन आव्यूह और समीकरणों का निर्धारक है σk:


 * $$\mathbf J = \left( \begin{array}{cccc}

\frac{\partial\sigma_1}{\partial\lambda_1} & \frac{\partial\sigma_1}{\partial\lambda_2} & \cdots & \frac{\partial\sigma_1}{\partial\lambda_n} \\[5pt] \frac{\partial\sigma_2}{\partial\lambda_1} & \frac{\partial\sigma_2}{\partial\lambda_2} & \cdots & \frac{\partial\sigma_2}{\partial\lambda_n} \\[5pt] \vdots & \vdots & \ddots & \vdots \\[5pt] \frac{\partial\sigma_n}{\partial\lambda_1} & \frac{\partial\sigma_n}{\partial\lambda_2} & \cdots & \frac{\partial\sigma_n}{\partial\lambda_n} \end{array}\right).$$ चूँकि सभी कण सभी व्यवरोधओं में योगदान नहीं करते हैं, $$\mathbf J_\sigma$$ एक ब्लॉक आव्यूह है और इसे आव्यूह की ब्लॉक-यूनिट में व्यक्तिगत रूप से हल किया जा सकता है। दूसरे शब्दों में, $$\mathbf J_\sigma$$ प्रत्येक अणु के लिए व्यक्तिगत रूप से हल किया जा सकता है।

सदिश को लगातार अपडेट करने के बजाय $$\underline{\lambda}$$, से पुनरावृत्ति प्रारंभ की जा सकती है $$\underline{\lambda}^{(0)} = \mathbf 0$$, जिसके परिणामस्वरूप सरल अभिव्यक्तियाँ प्राप्त होती हैं $$\sigma_k(t)$$ और $$\frac{\partial \sigma_k(t)}{\partial \lambda_j}$$. इस मामले में


 * $$ J_{ij} = \left.\frac{\partial\sigma_j}{\partial\lambda_i}\right|_{\mathbf \lambda=0} = 2\left[\hat{x}_{j\alpha} - \hat{x}_{j\beta}\right]\left[\frac{\partial \sigma_i}{\partial x_{j\alpha}} m_{j\alpha}^{-1} - \frac{\partial \sigma_i}{\partial x_{j\beta}} m_{j\beta}^{-1} \right]. $$

तब $$\lambda$$ को अद्यतन किया गया है
 * $$ \mathbf \lambda_j = - \mathbf J^{-1}\left[ \left\| \hat{\mathbf x}_{j\alpha}(t+\Delta t) - \hat{\mathbf x}_{j\beta}(t+\Delta t)\right\|^2 - d_j^2\right].$$

प्रत्येक पुनरावृत्ति के बाद, अप्रतिबंधित कण स्थितियों का उपयोग करके अद्यतन किया जाता है


 * $$\hat{\mathbf x}_i(t+\Delta t) \leftarrow \hat{\mathbf x}_i(t+\Delta t) + \sum_{k=1}^n \lambda_k\frac{\partial \sigma_k}{\partial \mathbf x_i} \left(\Delta t\right)^2m_i^{-1}.$$

फिर सदिश को रीसेट कर दिया जाता है


 * $$\underline{\lambda} = \mathbf 0.$$

उपरोक्त प्रक्रिया व्यवरोध समीकरणों के समाधान होने तक दोहराई जाती है, $$\sigma_k(t+\Delta t)$$, एक संख्यात्मक त्रुटि की निर्धारित सहनशीलता में परिवर्तित हो जाता है।

हालाँकि लैग्रेंज मल्टीप्लायरों की गणना करने के लिए कई एल्गोरिदम हैं, लेकिन ये अंतर केवल समीकरणों की प्रणाली को हल करने के तरीकों पर निर्भर करते हैं। इस विधि के लिए सामान्यतः अर्ध-न्यूटन विधियों का उपयोग किया जाता है।

सेटल एल्गोरिदम
सेटल (SETTLE) एल्गोरिथम गैर-रैखिक समीकरणों की प्रणाली को विश्लेषणात्मक रूप से हल करता है $$n=3$$ निरंतर समय में व्यवरोधएँ. यद्यपि यह बड़ी संख्या में व्यवरोधओं को मापता नहीं है, इसका उपयोग प्रायः दृढ़  जल के अणुओं को बाधित करने के लिए किया जाता है, जो लगभग सभी जैविक सिमुलेशन में उपलब्ध होते हैं और सामान्यतः तीन व्यवरोधओं (जैसे एसपीसी/ई और टीआईपी3पी जल मॉडल) का उपयोग करके तैयार किए जाते हैं।

शेक एल्गोरिदम
शेक (SHAKE) एल्गोरिथ्म को पहली बार आण्विक गतिकी सिमुलेशन के दौरान एक बंधन ज्यामिति व्यवरोध को संतुष्ट करने के लिए विकसित किया गया था। किसी भी होलोनोमिक व्यवरोध को संभालने के लिए विधि को सामान्यीकृत किया गया था, जैसे कि निरंतर बंधन कोण, या आणविक दृढता को बनाए रखने के लिए आवश्यक।

शेक एल्गोरिथ्म में, गैर-रैखिक व्यवरोध समीकरणों की प्रणाली को गॉस-सीडेल विधि का उपयोग करके हल किया जाता है जो न्यूटन पुनरावृत्तिl न्यूटन-रेफसन विधि का उपयोग करके समीकरणों की रैखिक प्रणाली के समाधान का अनुमान लगाता है;


 * $$\underline{\lambda} = -\mathbf J_\sigma^{-1} \underline{\sigma}.$$

यह ऐसा मानने के बराबर है $$\mathbf J_\sigma$$ विकर्ण रूप से प्रभावशाली है और हल कर रहा है $$k$$th के लिए समीकरण $$k$$ अज्ञात है। व्यवहार में, हम गणना करते हैं



\begin{align} \lambda_k & \leftarrow \frac{\sigma_k(t)}{\partial \sigma_k(t)/\partial \lambda_k}, \\[5pt] \mathbf x_{k\alpha} & \leftarrow \mathbf x_{k\alpha} + \lambda_k \frac{\partial \sigma_k(t)}{\partial \mathbf x_{k\alpha}}, \\[5pt] \mathbf x_{k\beta} & \leftarrow \mathbf x_{k\beta} + \lambda_k \frac{\partial \sigma_k(t)}{\partial \mathbf x_{k\beta}}, \end{align} $$ सभी के लिए $$k=1\ldots n$$ व्यवरोध समीकरणों तक पुनरावर्ती रूप से $$\sigma_k(t+\Delta t)$$ एक निश्चित सहिष्णुता के अनुसार हल किया जाता है।

प्रत्येक पुनरावृत्ति की गणना लागत है $$\mathcal O(n)$$, और पुनरावृत्तियाँ स्वयं रैखिक रूप से अभिसरण होती हैं।

बाद में शेक का एक अपुनरावृत्तीय रूप विकसित किया गया।

शेक एल्गोरिथम के कई प्रकार उपलब्ध हैं। यद्यपि वे स्वयं व्यवरोधओं की गणना या लागू करने के तरीके में भिन्न हैं, फिर भी व्यवरोधओं को लैग्रेंज मल्टीप्लायरों का उपयोग करके तैयार किया जाता है जिनकी गणना गॉस-सीडेल विधि का उपयोग करके की जाती है।

मूल शेक एल्गोरिदम दृढ़ और लचीले दोनों अणुओं (जैसे पानी, बेंजीन और बाइफिनाइल) को नियंत्रित करने में सक्षम है और आण्विक गतिकी सिमुलेशन में नगण्य त्रुटि या ऊर्जा बहाव पेश करता है। शेक के साथ एक मुद्दा यह है कि अभिसरण के एक निश्चित स्तर तक पहुंचने के लिए आवश्यक पुनरावृत्तियों की संख्या बढ़ जाती है क्योंकि आणविक ज्यामिति अधिक जटिल हो जाती है। 64 बिट कंप्यूटर सटीकता (सापेक्ष सहनशीलता) तक पहुंचने के लिए $$\approx 10^{-16}$$) 310K के तापमान पर एक विशिष्ट आण्विक गतिकी सिमुलेशन में, आणविक ज्यामिति को बनाए रखने के लिए 3 व्यवरोधओं वाले 3-साइट जल मॉडल को औसतन 9 पुनरावृत्तियों की आवश्यकता होती है (जो प्रति साइट प्रति समय-चरण 3 है)। 5 व्यवरोधओं वाले 4-साइट ब्यूटेन मॉडल को 17 पुनरावृत्तियों (22 प्रति साइट) की आवश्यकता होती है, 12 व्यवरोधओं वाले 6-साइट बेंजीन मॉडल को 36 पुनरावृत्तियों (72 प्रति साइट) की आवश्यकता होती है, जबकि 29 व्यवरोधओं वाले 12-साइट बाइफिनाइल मॉडल को 92 पुनरावृत्तियों की आवश्यकता होती है ( 229 प्रति साइट प्रति समय-चरण)। इसलिए शेक एल्गोरिदम की सीपीयू आवश्यकताएं महत्वपूर्ण हो सकती हैं, खासकर अगर आणविक मॉडल में उच्च स्तर की  दृढता हो।

विधि का एक बाद का विस्तार, QSHAKE (चार का समुदाय शेक) को कठोर इकाइयों से बने अणुओं के लिए एक तेज़ विकल्प के रूप में विकसित किया गया था, लेकिन यह सामान्य उद्देश्य के रूप में नहीं है। यह सुगंधित अंगूठी सिस्टम जैसे कठोर लूप के लिए संतोषजनक ढंग से काम करता है लेकिन QSHAKE लचीले लूप के लिए विफल रहता है, जैसे कि जब प्रोटीन में डाइसल्फ़ाइड बॉन्ड होता है।

जबकि रैटल शेक की तरह ही काम करता है, फिर भी वेलोसिटी वेरलेट समय एकीकरण योजना का उपयोग करते हुए, WIGGLE लैग्रेंज मल्टीप्लायरों के लिए प्रारंभिक अनुमान का उपयोग करके SHAKE और RATTLE का विस्तार करता है $$\lambda_k$$ कण वेग के आधार पर. उल्लेखनीय है कि MSHAKE बेहतर अभिसरण प्राप्त करने के लिए व्यवरोध बलों पर सुधार की गणना करता है।

SHAKE एल्गोरिथम का अंतिम संशोधन P-SHAKE एल्गोरिथम है जिसे बहुत दृढ या अर्ध-दृढ अणुओं पर लागू किया जाता है। P-SHAKE एक प्री-कंडीशनर की गणना और अद्यतन करता है जो SHAKE पुनरावृत्ति से पहले व्यवरोध ग्रेडिएंट्स पर लागू होता है, जिससे जैकोबियन होता है $$\mathbf J_\sigma$$ विकर्ण या दृढ़ता से विकर्ण रूप से प्रभावशाली बनना। इस प्रकार वियुग्मित व्यवरोधएं बहुत तेजी से (रैखिक रूप से विपरीत द्विघात रूप से) एकाग्र होती हैं $$\mathcal O(n^2)$$.

M-शेक एल्गोरिदम
M-शेक एल्गोरिदम सीधे न्यूटन की विधि का उपयोग करके समीकरणों की गैर-रेखीय प्रणाली को हल करता है। प्रत्येक पुनरावृत्ति में, समीकरणों की रैखिक प्रणाली


 * $$\underline{\lambda} = -\mathbf J_\sigma^{-1} \underline{\sigma}$$

एलयू अपघटन का उपयोग करके बिल्कुल हल किया जाता है। प्रत्येक पुनरावृत्ति की लागत होती है $$\mathcal O(n^3)$$ संचालन, फिर भी समाधान द्विघात अभिसरण को अभिसरण करता है, जिसके लिए SHAKE की तुलना में कम पुनरावृत्तियों की आवश्यकता होती है।

यह समाधान पहली बार 1986 में जियोवन्नी सिस्कोटी और रेकैर्ट द्वारा प्रस्तावित किया गया था शीर्षक के तहत आव्यूह विधि, फिर भी समीकरणों की रैखिक प्रणाली के समाधान में भिन्नता है। सिस्कोटी और रेकैर्ट आव्यूह को उलटने का सुझाव देते हैं $$\mathbf J_\sigma$$ प्रत्यक्ष रूप से, फिर भी ऐसा केवल एक बार, पहली पुनरावृत्ति में। पहले पुनरावृत्ति की लागत होती है $$\mathcal O(n^3)$$ संचालन, जबकि निम्नलिखित पुनरावृत्तियों की लागत केवल है $$\mathcal O(n^2)$$ संचालन (आव्यूह-सदिश गुणन के लिए)। हालाँकि यह सुधार एक लागत पर आता है, क्योंकि जैकोबियन अब अद्यतन नहीं है, अभिसरण केवल रैखिक अभिसरण है, भले ही SHAKE एल्गोरिथ्म की तुलना में बहुत तेज़ दर पर हो।

विरल आव्यूह तकनीकों पर आधारित इस दृष्टिकोण के कई प्रकारों का अध्ययन बार्थ एट अल द्वारा किया गया था।

आकार (SHAPE) एल्गोरिथ्म
आकार एल्गोरिथ्म तीन या अधिक केंद्रों के दृढ़ पिंडों को बाधित करने के लिए SHAKE का एक बहुकेंद्रीय एनालॉग है। SHAKE की तरह, एक अनियंत्रित कदम उठाया जाता है और फिर सीधे   रिजिड बॉडी रोटेशन मैट्रिक्स (दृढ़ पिंड परिक्रमण आव्यूह) की गणना और लागू करके सही किया जाता है जो संतुष्ट करता है:


 * $$ L^\text{rigid} \left( t + \frac{\Delta t} 2 \right) = L^\text{nonrigid} \left( t + \frac{\Delta t} 2 \right)$$

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

लिंक्स एल्गोरिदम
एक वैकल्पिक व्यवरोध विधि, LINCS (रैखिक व्यवरोध सॉल्वर) 1997 में हेस, बेकर, बेरेन्डसेन और फ्रैजे द्वारा विकसित की गई थी। और यह एडबर्ग, इवांस और मॉरिस (ईईएम) की 1986 पद्धति पर आधारित था।

LINCS लैग्रेंज मल्टीप्लायरों को व्यवरोध बलों पर लागू करता है और जैकोबियन के व्युत्क्रम का अनुमान लगाने के लिए श्रृंखला विस्तार का उपयोग करके मल्टीप्लायरों का समाधान करता है। $$\mathbf J_\sigma$$:


 * $$(\mathbf I - \mathbf J_\sigma)^{-1} = \mathbf I + \mathbf J_\sigma + \mathbf J_\sigma^2 + \mathbf J_\sigma^3 + \cdots$$

न्यूटन पुनरावृत्ति के प्रत्येक चरण में। यह सन्निकटन केवल 1 से छोटे eigenvalues ​​​​वाले आव्यूह के लिए काम करता है, जिससे LINCS एल्गोरिदम केवल कम कनेक्टिविटी वाले अणुओं के लिए उपयुक्त हो जाता है।

बताया गया है कि LINCS, SHAKE से 3-4 गुना तेज़ है।

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

यह भी देखें

 * आण्विक गतिकी
 * आणविक यांत्रिकी मॉडलिंग के लिए सॉफ्टवेयर की सूची