सहानुभूतिपूर्ण इंटीग्रेटर: Difference between revisions

From Vigyanwiki
No edit summary
No edit summary
Line 1: Line 1:
गणित में, सिम्पलेक्टिक इंटीग्रेटर (एसआई) हैमिल्टनियन प्रणालियों के लिए एक संख्यात्मक साधारण अवकल समीकरण है। सिम्पलेक्टिक इंटीग्रेटर्स [[ज्यामितीय इंटीग्रेटर]]्स के उपवर्ग का निर्माण करते हैं, जो परिभाषा के अनुसार, [[ विहित परिवर्तन ]] हैं। वे व्यापक रूप से गैर-रैखिक गतिशीलता, आणविक गतिशीलता, [[असतत तत्व विधि]]यों, [[कण त्वरक]], [[प्लाज्मा भौतिकी]], [[क्वांटम भौतिकी]] और [[आकाशीय यांत्रिकी]] में उपयोग किए जाते हैं। '''सिम्पलेक्टिक इंटीग्रेटर्स [[ज्यामितीय इंटीग्रेटर]]्स के उपवर्ग का निर्माण करते हैं, जो परिभाषा के अनुसार, [[ विहित परिवर्तन | विहित परिवर्तन]] हैं। वे व्यापक रूप से गैर-रैखिक गतिशीलता, आणविक गतिशीलता, [[असतत तत्व विधि]]यों, [[कण त्वरक]], [[प्लाज्मा भौतिकी]], [[क्वांटम भौतिकी]] और [[आकाशीय यांत्रिकी]] में उपयोग किए जाते हैं।'''
गणित में, सिम्पलेक्टिक इंटीग्रेटर (एसआई) हैमिल्टनियन प्रणालियों के लिए संख्यात्मक साधारण अवकल समीकरण है। सिम्पलेक्टिक इंटीग्रेटर्स [[ज्यामितीय इंटीग्रेटर]] के उपवर्ग का निर्माण करते हैं, जो परिभाषा के अनुसार, [[ विहित परिवर्तन ]] हैं। वे व्यापक रूप से गैर-रैखिक गतिशीलता, आणविक गतिशीलता, [[असतत तत्व विधि]]यों, [[कण त्वरक]], [[प्लाज्मा भौतिकी]], [[क्वांटम भौतिकी]] और [[आकाशीय यांत्रिकी]] में उपयोग किए जाते हैं।  


== परिचय ==
== परिचय ==
Line 9: Line 9:
स्थिति और संवेग का सेट निर्देशांक <math>(q,p)</math> [[विहित निर्देशांक]] कहलाते हैं। (अधिक पृष्ठभूमि के लिए [[हैमिल्टनियन यांत्रिकी]] देखें।)
स्थिति और संवेग का सेट निर्देशांक <math>(q,p)</math> [[विहित निर्देशांक]] कहलाते हैं। (अधिक पृष्ठभूमि के लिए [[हैमिल्टनियन यांत्रिकी]] देखें।)


हैमिल्टन के समीकरणों का समय विकास एक [[sympletomorphism|सिम्प्टोमोर्फिज़्म]] है, जिसका अर्थ है कि यह सिम्पलेक्टिक [[2-प्रपत्र]] का संरक्षण करता है <math>dp \wedge dq</math>. एक संख्यात्मक योजना एक सिम्पलेक्टिक इंटीग्रेटर है यदि यह इस 2-फॉर्म को भी संरक्षित करती है।
हैमिल्टन के समीकरणों का समय विकास [[sympletomorphism|सिम्प्टोमोर्फिज़्म]] है, जिसका अर्थ है कि यह सिम्पलेक्टिक [[2-प्रपत्र]] का संरक्षण करता है <math>dp \wedge dq</math>. संख्यात्मक योजना सिम्पलेक्टिक इंटीग्रेटर है यदि यह इस 2-फॉर्म को भी संरक्षित करती है।


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


आदिम [[यूलर एकीकरण]] और क्लासिकल रनगे-कुट्टा योजना जैसी अधिकांश सामान्य संख्यात्मक विधियां, सिम्पलेक्टिक समाकलक नहीं हैं।
आदिम [[यूलर एकीकरण]] और क्लासिकल रनगे-कुट्टा योजना जैसी अधिकांश सामान्य संख्यात्मक विधियां, सिम्पलेक्टिक समाकलक नहीं हैं।
Line 19: Line 19:
=== वियोज्य हैमिल्टन के लिए विभाजन विधियाँ ===
=== वियोज्य हैमिल्टन के लिए विभाजन विधियाँ ===


बंटवारे की विधियों से सिम्पलेक्टिक इंटीग्रेटर्स का एक व्यापक रूप से उपयोग किया जाने वाला वर्ग बनता है।
बंटवारे की विधियों से सिम्पलेक्टिक इंटीग्रेटर्स का व्यापक रूप से उपयोग किया जाने वाला वर्ग बनता है।


मान लें कि हैमिल्टन वियोज्य है, जिसका अर्थ है कि इसे रूप में लिखा जा सकता है
मान लें कि हैमिल्टन वियोज्य है, जिसका अर्थ है कि इसे रूप में लिखा जा सकता है
Line 33: Line 33:
{{NumBlk|:|<math>\dot{z}=\{z,H(z)\},</math>|{{EquationRef|2}}}}
{{NumBlk|:|<math>\dot{z}=\{z,H(z)\},</math>|{{EquationRef|2}}}}


जहाँ <math>\{\cdot, \cdot\}</math> पॉसों कोष्ठक है। इसके अतिरिक्त, एक ऑपरेटर का प्रारंभ करके <math>D_H \cdot = \{\cdot, H\}</math>, जो हेमिल्टनियन यांत्रिकी गणितीय औपचारिकता के साथ ऑपरेंड का पॉइसन ब्रैकेट लौटाता है, हैमिल्टन के समीकरण की अभिव्यक्ति को और सरल बनाया जा सकता है
जहाँ <math>\{\cdot, \cdot\}</math> पॉसों कोष्ठक है। इसके अतिरिक्त, ऑपरेटर का प्रारंभ करके <math>D_H \cdot = \{\cdot, H\}</math>, जो हेमिल्टनियन यांत्रिकी गणितीय औपचारिकता के साथ ऑपरेंड का पॉइसन ब्रैकेट लौटाता है, हैमिल्टन के समीकरण की अभिव्यक्ति को और सरल बनाया जा सकता है


:<math>\dot{z}=D_H z.</math>
:<math>\dot{z}=D_H z.</math>
Line 45: Line 45:
{{NumBlk|:|<math>z(\tau) = \exp[\tau (D_T + D_V)]z(0).</math>|{{EquationRef|4}}}}
{{NumBlk|:|<math>z(\tau) = \exp[\tau (D_T + D_V)]z(0).</math>|{{EquationRef|4}}}}


एसआई योजना समय-विकास ऑपरेटर का अनुमान लगाती है <math>\exp[\tau (D_T + D_V)]</math> औपचारिक समाधान में ({{EquationNote|4}}) ऑपरेटरों के एक गुणांक के रूप में
एसआई योजना समय-विकास ऑपरेटर का अनुमान लगाती है <math>\exp[\tau (D_T + D_V)]</math> औपचारिक समाधान में ({{EquationNote|4}}) ऑपरेटरों के गुणांक के रूप में


{{NumBlk|:|<math>
{{NumBlk|:|<math>
Line 54: Line 54:
</math>|{{EquationRef|5}}}}
</math>|{{EquationRef|5}}}}


जहाँ <math>c_i</math> और <math>d_i</math> वास्तविक संख्याएँ हैं, <math>k</math> एक पूर्णांक है, जिसे इंटीग्रेटर का क्रम कहा जाता है, और जहाँ <math display="inline">\sum_{i=1}^k c_i = \sum_{i=1}^k d_i = 1</math>. ध्यान दें कि प्रत्येक ऑपरेटर <math>\exp(c_i \tau D_T)</math> और <math>\exp(d_i \tau D_V)</math> एक समानता प्रदान करता है, इसलिए उनका गुणांक दाईं ओर प्रदर्शित होता है ({{EquationNote|5}}) भी एक सिम्पलेक्टिक मानचित्र बनाता है।
जहाँ <math>c_i</math> और <math>d_i</math> वास्तविक संख्याएँ हैं, <math>k</math> पूर्णांक है, जिसे इंटीग्रेटर का क्रम कहा जाता है, और जहाँ <math display="inline">\sum_{i=1}^k c_i = \sum_{i=1}^k d_i = 1</math>. ध्यान दें कि प्रत्येक ऑपरेटर <math>\exp(c_i \tau D_T)</math> और <math>\exp(d_i \tau D_V)</math> समानता प्रदान करता है, इसलिए उनका गुणांक दाईं ओर प्रदर्शित होता है ({{EquationNote|5}}) भी सिम्पलेक्टिक मानचित्र बनाता है।


तब से <math>D_T^2 z = \{\{z,T\},T\} = \{(\dot{q}, 0),T\} = (0,0)</math> सभी के लिए <math>z</math>, हम यह निष्कर्ष निकाल सकते हैं कि
तब से <math>D_T^2 z = \{\{z,T\},T\} = \{(\dot{q}, 0),T\} = (0,0)</math> सभी के लिए <math>z</math>, हम यह निष्कर्ष निकाल सकते हैं कि
Line 64: Line 64:
{{NumBlk|:|<math>\exp(a D_T) = \sum_{n=0}^\infty \frac{(a D_T)^n}{n!},</math>|{{EquationRef|7}}}}
{{NumBlk|:|<math>\exp(a D_T) = \sum_{n=0}^\infty \frac{(a D_T)^n}{n!},</math>|{{EquationRef|7}}}}


जहाँ <math>a</math> एक इच्छानुसार वास्तविक संख्या है। संयोजन ({{EquationNote|6}}) और ({{EquationNote|7}}), और के लिए समान तर्क का उपयोग करके <math>D_V</math> जैसा कि हमने <math>D_T</math> उपयोग किया है, हम पाते हैं
जहाँ <math>a</math> इच्छानुसार वास्तविक संख्या है। संयोजन ({{EquationNote|6}}) और ({{EquationNote|7}}), और के लिए समान तर्क का उपयोग करके <math>D_V</math> जैसा कि हमने <math>D_T</math> उपयोग किया है, हम पाते हैं


{{NumBlk|:|<math>\begin{cases}
{{NumBlk|:|<math>\begin{cases}
Line 104: Line 104:
जहाँ <math>F(x)</math> पर बल सदिश है <math>x</math>, <math>a(x)</math> त्वरण वेक्टर है <math>x</math>, और <math>m</math> द्रव्यमान की अदिश राशि है।
जहाँ <math>F(x)</math> पर बल सदिश है <math>x</math>, <math>a(x)</math> त्वरण वेक्टर है <math>x</math>, और <math>m</math> द्रव्यमान की अदिश राशि है।


कई सिम्पलेक्टिक इंटीग्रेटर्स नीचे दिए गए हैं। उनका उपयोग करने का एक उदाहरण तरीका स्थिति के साथ एक कण पर विचार करना है <math>q</math> और गति <math>p</math> है
कई सिम्पलेक्टिक इंटीग्रेटर्स नीचे दिए गए हैं। उनका उपयोग करने का उदाहरण विधियाँ स्थिति के साथ कण पर विचार करना है <math>q</math> और गति <math>p</math> है


मूल्यों के साथ एक टाइमस्टेप प्रयुक्त करने के लिए <math>c_{1,2,3}, d_{1,2,3}</math> कण के लिए, निम्न चरणों का पालन करें: क्रमिक रूप से:
मूल्यों के साथ टाइमस्टेप प्रयुक्त करने के लिए <math>c_{1,2,3}, d_{1,2,3}</math> कण के लिए, निम्न चरणों का पालन करें: क्रमिक रूप से:
* स्थिति अद्यतन करें <math>i</math> इसके (पहले अद्यतन) वेग को जोड़कर कण का <math>i</math> से गुणा <math>c_i</math>
* स्थिति अद्यतन करें <math>i</math> इसके (पहले अद्यतन) वेग को जोड़कर कण का <math>i</math> से गुणा <math>c_i</math>
* वेग अद्यतन करें <math>i</math> कण का इसमें त्वरण (अद्यतन स्थिति में) जोड़कर गुणा किया जाता है <math>d_i</math>
* वेग अद्यतन करें <math>i</math> कण का इसमें त्वरण (अद्यतन स्थिति में) जोड़कर गुणा किया जाता है <math>d_i</math>
Line 114: Line 114:
[[सहानुभूतिपूर्ण यूलर विधि|सिम्पलेक्टिक यूलर विधि]] प्रथम-क्रम समाकलक है <math>k=1</math> और गुणांक
[[सहानुभूतिपूर्ण यूलर विधि|सिम्पलेक्टिक यूलर विधि]] प्रथम-क्रम समाकलक है <math>k=1</math> और गुणांक
:<math>c_1 = d_1 = 1.</math>
:<math>c_1 = d_1 = 1.</math>
ध्यान दें कि समय-प्रतिवर्तीता की आवश्यकता होने पर उपरोक्त एल्गोरिदम काम नहीं करता है। एल्गोरिथ्म को दो भागों में प्रयुक्त किया जाना है, एक सकारात्मक समय चरणों के लिए, एक नकारात्मक समय चरणों के लिए।
ध्यान दें कि समय-प्रतिवर्तीता की आवश्यकता होने पर उपरोक्त एल्गोरिदम काम नहीं करता है। एल्गोरिथ्म को दो भागों में प्रयुक्त किया जाना है, सकारात्मक समय चरणों के लिए नकारात्मक समय चरणों के लिए।


==== दूसरे क्रम का उदाहरण ====
==== दूसरे क्रम का उदाहरण ====
Line 123: Line 123:
==== तीसरे क्रम का उदाहरण ====
==== तीसरे क्रम का उदाहरण ====


एक तीसरा क्रम सिम्पलेक्टिक इंटीग्रेटर (के साथ <math>k=3</math>) की खोज 1983 में रोनाल्ड रुथ ने की थी।<ref>{{cite journal
तीसरा क्रम सिम्पलेक्टिक इंटीग्रेटर (के साथ <math>k=3</math>) की खोज 1983 में रोनाल्ड रुथ ने की थी।<ref>{{cite journal
|last = Ruth
|last = Ruth
|first = Ronald D.
|first = Ronald D.
Line 146: Line 146:


==== चौथे क्रम का उदाहरण ====
==== चौथे क्रम का उदाहरण ====
एक चौथे क्रम के इंटीग्रेटर (के साथ <math>k=4</math>) भी 1983 में रूथ द्वारा खोजा गया था और उस समय कण-त्वरक समुदाय को सामान्यतः वितरित किया गया था। वन द्वारा जीवंत समीक्षा लेख में इसका वर्णन किया गया था।<ref>{{cite journal
चौथे क्रम के इंटीग्रेटर (के साथ <math>k=4</math>) भी 1983 में रूथ द्वारा खोजा गया था और उस समय कण-त्वरक समुदाय को सामान्यतः वितरित किया गया था। वन द्वारा जीवंत समीक्षा लेख में इसका वर्णन किया गया था।<ref>{{cite journal
|last = Forest
|last = Forest
|first = Etienne
|first = Etienne
Line 199: Line 199:
\end{align}
\end{align}
</math>
</math>
इन गुणांकों को निर्धारित करने के लिए, बेकर-कैंपबेल-हॉसडॉर्फ सूत्र का उपयोग किया जा सकता है। योशिदा, विशेष रूप से, उच्च-क्रम के समाकलकों के लिए गुणांकों की एक सुंदर व्युत्पत्ति देता है। बाद में, ब्लेन्स और मून<ref>{{cite journal| last=Blanes|first=S.| author2=Moan, P. C. |title=Practical symplectic partitioned Runge–Kutta and Runge–Kutta–Nyström methods| journal=Journal of Computational and Applied Mathematics|date=May 2002| volume=142| issue=2|pages=313–330| doi=10.1016/S0377-0427(01)00492-7| bibcode = 2002JCoAM.142..313B |doi-access=free}}</ref> बहुत कम त्रुटि स्थिरांक वाले वियोज्य हैमिल्टनियन के साथ प्रणाली के एकीकरण के लिए आगे विकसित विभाजित रनगे-कुट्टा विधियाँ।
इन गुणांकों को निर्धारित करने के लिए, बेकर-कैंपबेल-हॉसडॉर्फ सूत्र का उपयोग किया जा सकता है। योशिदा, विशेष रूप से, उच्च-क्रम के समाकलकों के लिए गुणांकों की सुंदर व्युत्पत्ति देता है। बाद में, ब्लेन्स और मून<ref>{{cite journal| last=Blanes|first=S.| author2=Moan, P. C. |title=Practical symplectic partitioned Runge–Kutta and Runge–Kutta–Nyström methods| journal=Journal of Computational and Applied Mathematics|date=May 2002| volume=142| issue=2|pages=313–330| doi=10.1016/S0377-0427(01)00492-7| bibcode = 2002JCoAM.142..313B |doi-access=free}}</ref> बहुत कम त्रुटि स्थिरांक वाले वियोज्य हैमिल्टनियन के साथ प्रणाली के एकीकरण के लिए आगे विकसित विभाजित रनगे-कुट्टा विधियाँ।


=== === सामान्य अविभाज्य हैमिल्टनियन === के लिए विभाजन विधियाँ ===
=== सामान्य अविभाज्य हैमिल्टनियन के लिए विभाजन विधियाँ ===
'''सामान्य अविभाज्य हैमिल्टनियन भी स्पष्ट और सिम्पलेक्टिक रूप से एकीकृत हो सकते हैं।'''
'''सामान्य अविभाज्य हैमिल्टनियन भी स्पष्ट और सिम्पलेक्टिक रूप से एकीकृत हो सकते हैं।'''


ऐसा करने के लिए, ताओ ने एक अवरोध प्रस्तुत किया जो इस तरह के प्रणाली के स्पष्ट विभाजन को सक्षम करने के लिए चरण स्थान की दो प्रतियों को एक साथ बांधता है।<ref>{{cite journal
ऐसा करने के लिए, ताओ ने अवरोध प्रस्तुत किया जो इस तरह के प्रणाली के स्पष्ट विभाजन को सक्षम करने के लिए चरण स्थान की दो प्रतियों को साथ बांधता है।<ref>{{cite journal
|last = Tao
|last = Tao
|first = Molei
|first = Molei
Line 221: Line 221:




इसके अतिरिक्त विचार है <math>H(Q,P)</math>, एक अनुकरण करता है <math>\bar{H}(q,p,x,y) = H(q,y) + H(x,p) + \omega \left(\left\|q-x\right\|_2^2/2 + \left\|p-y\right\|_2^2/2\right)</math>, जिसका समाधान इससे सहमत है <math>H(Q,P)</math> इस अर्थ में कि <math>q(t)=x(t)=Q(t),p(t)=y(t)=P(t)</math>.
इसके अतिरिक्त विचार है <math>H(Q,P)</math>, अनुकरण करता है <math>\bar{H}(q,p,x,y) = H(q,y) + H(x,p) + \omega \left(\left\|q-x\right\|_2^2/2 + \left\|p-y\right\|_2^2/2\right)</math>, जिसका समाधान इससे सहमत है <math>H(Q,P)</math> इस अर्थ में कि <math>q(t)=x(t)=Q(t),p(t)=y(t)=P(t)</math>.


नया हैमिल्टनियन स्पष्ट सिम्पलेक्टिक एकीकरण के लिए लाभदायक है, क्योंकि इसे तीन उप-हैमिल्टनियों के योग में विभाजित किया जा सकता है, <math>H_A=H(q,y)</math>, <math>H_B=H(x,p)</math>, और <math>H_C = \omega \left(\left\|q-x\right\|_2^2/2 + \left\|p-y\right\|_2^2/2\right)</math>. सभी तीन उप-हैमिल्टनियों का सटीक समाधान स्पष्ट रूप से प्राप्त किया जा सकता है: दोनों <math>H_A, H_B</math> समाधान बेमेल स्थिति और गति के बदलाव के अनुरूप हैं, और <math>H_C</math> एक रेखीय परिवर्तन के अनुरूप है। प्रणाली को सिम्पलेक्टिक रूप से अनुकरण करने के लिए, बस इन समाधान मानचित्रों की रचना करता है।
नया हैमिल्टनियन स्पष्ट सिम्पलेक्टिक एकीकरण के लिए लाभदायक है, क्योंकि इसे तीन उप-हैमिल्टनियों के योग में विभाजित किया जा सकता है, <math>H_A=H(q,y)</math>, <math>H_B=H(x,p)</math>, और <math>H_C = \omega \left(\left\|q-x\right\|_2^2/2 + \left\|p-y\right\|_2^2/2\right)</math>. सभी तीन उप-हैमिल्टनियों का सटीक समाधान स्पष्ट रूप से प्राप्त किया जा सकता है: दोनों <math>H_A, H_B</math> समाधान बेमेल स्थिति और गति के बदलाव के अनुरूप हैं, और <math>H_C</math> रेखीय परिवर्तन के अनुरूप है। प्रणाली को सिम्पलेक्टिक रूप से अनुकरण करने के लिए, बस इन समाधान मानचित्रों की रचना करता है।


== अनुप्रयोग ==
== अनुप्रयोग ==


=== प्लाज्मा भौतिकी में ===
=== प्लाज्मा भौतिकी में ===
हाल के दशकों में प्लाज्मा भौतिकी में सिम्पलेक्टिक इंटीग्रेटर एक सक्रिय शोध विषय बन गया है,<ref>{{cite journal| last=Qin| first=H.|author2=Guan,X.|journal=Physical Review Letters|date=2008|volume=100|page=035006|title=सामान्य चुंबकीय क्षेत्रों में लंबे समय के सिमुलेशन के लिए आवेशित कणों के मार्गदर्शक केंद्र गति के लिए एक परिवर्तनशील सहानुभूतिपूर्ण इंटीग्रेटर| issue=3| doi=10.1103/PhysRevLett.100.035006| pmid=18232993| url=https://digital.library.unt.edu/ark:/67531/metadc932923/m2/1/high_res_d/960290.pdf}}</ref> क्योंकि मानक सिम्पलेक्टिक विधियों के सीधे-सीधे अनुप्रयोग पेटा- से एक्स-स्केल कंप्यूटिंग हार्डवेयर द्वारा सक्षम बड़े पैमाने के प्लाज्मा सिमुलेशन की आवश्यकता के अनुरूप नहीं हैं। जांच के तहत भौतिकी समस्या की विशेष संरचनाओं में टैप करते हुए, विशेष सिम्पलेक्टिक एल्गोरिदम को कस्टम रूप से डिज़ाइन करने की आवश्यकता है। ऐसा ही एक उदाहरण विद्युत चुम्बकीय क्षेत्र में आवेशित कण गतिकी है। विहित सिम्पलेक्टिक संरचना के साथ, गतिकी का हैमिल्टनियन है <math display="block">H(\boldsymbol{p},\boldsymbol{x})=\frac{1}{2}\left(\boldsymbol{p}-\boldsymbol{A}\right)^{2}+\phi,</math> जिसकी <math display="inline">\boldsymbol{p}</math>-निर्भरता और <math display="inline">\boldsymbol{x}</math>निर्भरता वियोज्य नहीं हैं, और मानक स्पष्ट सिम्पलेक्टिक विधियों प्रयुक्त नहीं होते हैं। बड़े पैमाने पर समांतर समूहों पर बड़े पैमाने पर सिमुलेशन के लिए, चुकीं, स्पष्ट तरीकों को प्राथमिकता दी जाती है।
हाल के दशकों में प्लाज्मा भौतिकी में सिम्पलेक्टिक इंटीग्रेटर सक्रिय शोध विषय बन गया है,<ref>{{cite journal| last=Qin| first=H.|author2=Guan,X.|journal=Physical Review Letters|date=2008|volume=100|page=035006|title=सामान्य चुंबकीय क्षेत्रों में लंबे समय के सिमुलेशन के लिए आवेशित कणों के मार्गदर्शक केंद्र गति के लिए एक परिवर्तनशील सहानुभूतिपूर्ण इंटीग्रेटर| issue=3| doi=10.1103/PhysRevLett.100.035006| pmid=18232993| url=https://digital.library.unt.edu/ark:/67531/metadc932923/m2/1/high_res_d/960290.pdf}}</ref> क्योंकि मानक सिम्पलेक्टिक विधियों के सीधे-सीधे अनुप्रयोग पेटा- से एक्स-स्केल कंप्यूटिंग हार्डवेयर द्वारा सक्षम बड़े पैमाने के प्लाज्मा सिमुलेशन की आवश्यकता के अनुरूप नहीं हैं। जांच के तहत भौतिकी समस्या की विशेष संरचनाओं में टैप करते हुए, विशेष सिम्पलेक्टिक एल्गोरिदम को कस्टम रूप से डिज़ाइन करने की आवश्यकता है। ऐसा ही उदाहरण विद्युत चुम्बकीय क्षेत्र में आवेशित कण गतिकी है। विहित सिम्पलेक्टिक संरचना के साथ, गतिकी का हैमिल्टनियन है <math display="block">H(\boldsymbol{p},\boldsymbol{x})=\frac{1}{2}\left(\boldsymbol{p}-\boldsymbol{A}\right)^{2}+\phi,</math> जिसकी <math display="inline">\boldsymbol{p}</math>-निर्भरता और <math display="inline">\boldsymbol{x}</math>निर्भरता वियोज्य नहीं हैं, और मानक स्पष्ट सिम्पलेक्टिक विधियों प्रयुक्त नहीं होते हैं। बड़े पैमाने पर समांतर समूहों पर बड़े पैमाने पर सिमुलेशन के लिए, चुकीं, स्पष्ट तरीकों को प्राथमिकता दी जाती है।
इस कठिनाई को दूर करने के लिए, हम विशिष्ट विधियों का पता लगा सकते हैं <math display="inline">\boldsymbol{p}</math>-निर्भरता और <math display="inline">\boldsymbol{x}</math>-निर्भरता इस हैमिल्टनियन में उलझी हुई है, और केवल इस या इस प्रकार की समस्या के लिए एक सिम्पलेक्टिक एल्गोरिथम डिजाइन करने का प्रयास करें। सबसे पहले, हम ध्यान दें कि <math display="inline">\boldsymbol{p}</math>-निर्भरता द्विघात है, इसलिए पहले क्रम की सहानुभूति यूलर विधि में निहित है <math display="inline">\boldsymbol{p}</math> वास्तव में मुखर है। यह वही है जो कैनोनिकल सिम्प्लेक्टिक [[ कण-इन-सेल ]] (पीआईसी) एल्गोरिथम में उपयोग किया जाता है।<ref>{{cite journal| last=Qin|first=H.|author2=Liu, J.|author3=Xiao,J.|journal=Nuclear Fusion|date=2016|volume=56|page=014001|title=Canonical symplectic particle-in-cell method for long-term large-scale simulations of the Vlasov–Maxwell equations| issue=1| doi=10.1088/0029-5515/56/1/014001|arxiv=1503.08334|bibcode=2016NucFu..56a4001Q|s2cid=29190330 }}</ref> उच्च क्रम स्पष्ट विधियों का निर्माण करने के लिए, हम आगे ध्यान दें कि <math display="inline">\boldsymbol{p}</math>-निर्भरता और <math display="inline">\boldsymbol{x}</math>-इसमें निर्भरता <math display="inline">H(\boldsymbol{p},\boldsymbol{x})</math> गुणांक-वियोज्य हैं, दूसरे और तीसरे क्रम के स्पष्ट सिम्पलेक्टिक एल्गोरिदम का निर्माण कार्यों का उपयोग करके किया जा सकता है,<ref>{{cite journal|last=Zhang|first=R.|author2=Qin, H.|author3=Tang, Y.|journal=Physical Review E|date=2016|volume=94|title=चार्ज कण गतिकी के लिए कार्यों को उत्पन्न करने के आधार पर स्पष्ट सहानुभूतिपूर्ण एल्गोरिदम| issue=1| page=013205| doi=10.1103/PhysRevE.94.013205| pmid=27575228|arxiv=1604.02787|bibcode=2016PhRvE..94a3205Z |s2cid=2166879 }}</ref> और समय-निर्भर विद्युत चुम्बकीय क्षेत्रों के लिए मनमाने ढंग से उच्च-क्रम वाले स्पष्ट सिम्पलेक्टिक इंटीग्रेटर्स का निर्माण रनगे-कुट्टा तकनीकों का उपयोग करके भी किया जा सकता है।<ref>{{cite journal|last=Tao|first=M.|journal=Journal of Computational Physics|date=2016|volume=327|title=सामान्य विद्युतचुंबकीय क्षेत्रों में आवेशित कणों के लिए स्पष्ट उच्च-आदेश सहानुभूतिपूर्ण इंटीग्रेटर्स|page=245|doi=10.1016/j.jcp.2016.09.047|arxiv=1605.01458|bibcode=2016JCoPh.327..245T |s2cid=31262651 }}</ref>
इस कठिनाई को दूर करने के लिए, हम विशिष्ट विधियों का पता लगा सकते हैं <math display="inline">\boldsymbol{p}</math>-निर्भरता और <math display="inline">\boldsymbol{x}</math>-निर्भरता इस हैमिल्टनियन में उलझी हुई है, और केवल इस या इस प्रकार की समस्या के लिए सिम्पलेक्टिक एल्गोरिथम डिजाइन करने का प्रयास करें। सबसे पहले, हम ध्यान दें कि <math display="inline">\boldsymbol{p}</math>-निर्भरता द्विघात है, इसलिए पहले क्रम की सहानुभूति यूलर विधि में निहित है <math display="inline">\boldsymbol{p}</math> वास्तव में मुखर है। यह वही है जो कैनोनिकल सिम्प्लेक्टिक [[ कण-इन-सेल ]] (पीआईसी) एल्गोरिथम में उपयोग किया जाता है।<ref>{{cite journal| last=Qin|first=H.|author2=Liu, J.|author3=Xiao,J.|journal=Nuclear Fusion|date=2016|volume=56|page=014001|title=Canonical symplectic particle-in-cell method for long-term large-scale simulations of the Vlasov–Maxwell equations| issue=1| doi=10.1088/0029-5515/56/1/014001|arxiv=1503.08334|bibcode=2016NucFu..56a4001Q|s2cid=29190330 }}</ref> उच्च क्रम स्पष्ट विधियों का निर्माण करने के लिए, हम आगे ध्यान दें कि <math display="inline">\boldsymbol{p}</math>-निर्भरता और <math display="inline">\boldsymbol{x}</math>-इसमें निर्भरता <math display="inline">H(\boldsymbol{p},\boldsymbol{x})</math> गुणांक-वियोज्य हैं, दूसरे और तीसरे क्रम के स्पष्ट सिम्पलेक्टिक एल्गोरिदम का निर्माण कार्यों का उपयोग करके किया जा सकता है,<ref>{{cite journal|last=Zhang|first=R.|author2=Qin, H.|author3=Tang, Y.|journal=Physical Review E|date=2016|volume=94|title=चार्ज कण गतिकी के लिए कार्यों को उत्पन्न करने के आधार पर स्पष्ट सहानुभूतिपूर्ण एल्गोरिदम| issue=1| page=013205| doi=10.1103/PhysRevE.94.013205| pmid=27575228|arxiv=1604.02787|bibcode=2016PhRvE..94a3205Z |s2cid=2166879 }}</ref> और समय-निर्भर विद्युत चुम्बकीय क्षेत्रों के लिए मनमाने ढंग से उच्च-क्रम वाले स्पष्ट सिम्पलेक्टिक इंटीग्रेटर्स का निर्माण रनगे-कुट्टा तकनीकों का उपयोग करके भी किया जा सकता है।<ref>{{cite journal|last=Tao|first=M.|journal=Journal of Computational Physics|date=2016|volume=327|title=सामान्य विद्युतचुंबकीय क्षेत्रों में आवेशित कणों के लिए स्पष्ट उच्च-आदेश सहानुभूतिपूर्ण इंटीग्रेटर्स|page=245|doi=10.1016/j.jcp.2016.09.047|arxiv=1605.01458|bibcode=2016JCoPh.327..245T |s2cid=31262651 }}</ref>
समस्या के निम्नलिखित गैर-विहित सहानुभूति संरचना को देखने के लिए एक अधिक सुरुचिपूर्ण और बहुमुखी विकल्प है, <math display="block">i_{(\dot{\boldsymbol{x}},\dot{\boldsymbol{v}})}\Omega=-dH,\ \ \ \Omega=d(\boldsymbol{v}+\boldsymbol{A})\wedge d\boldsymbol{x},\ \ \ H=\frac{1}{2}\boldsymbol{v}^{2}+\phi.</math> यहाँ <math display="inline">\Omega</math> एक गैर-निरंतर गैर-विहित सिम्पलेक्टिक रूप है। गैर-निरंतर गैर-कैनोनिकल सिम्प्लेक्टिक संरचना के लिए सामान्य सिम्प्लेक्टिक इंटीग्रेटर, स्पष्ट या अंतर्निहित, उपस्थित नहीं है। चुकीं, इस विशिष्ट समस्या के लिए, हे विभाजन विधि का उपयोग करके उच्च-क्रम स्पष्ट गैर-कैनोनिकल सिम्प्लेक्टिक इंटीग्रेटर्स का एक परिवार बनाया जा सकता है।<ref>{{cite journal|last=He|first=Y.|author2=Qin, H.|author3=Sun, Y.|journal=Physics of Plasmas|date=2015|volume=22|title=वेलासोव-मैक्सवेल समीकरणों के लिए हैमिल्टनियन एकीकरण के तरीके|page=124503| doi=10.1063/1.4938034 |arxiv=1505.06076 |s2cid=118560512 }}</ref> विभाजन <math display="inline">H</math> 4 भागों में, <math display="block">\begin{aligned}
समस्या के निम्नलिखित गैर-विहित सहानुभूति संरचना को देखने के लिए अधिक सुरुचिपूर्ण और बहुमुखी विकल्प है, <math display="block">i_{(\dot{\boldsymbol{x}},\dot{\boldsymbol{v}})}\Omega=-dH,\ \ \ \Omega=d(\boldsymbol{v}+\boldsymbol{A})\wedge d\boldsymbol{x},\ \ \ H=\frac{1}{2}\boldsymbol{v}^{2}+\phi.</math> यहाँ <math display="inline">\Omega</math> गैर-निरंतर गैर-विहित सिम्पलेक्टिक रूप है। गैर-निरंतर गैर-कैनोनिकल सिम्प्लेक्टिक संरचना के लिए सामान्य सिम्प्लेक्टिक इंटीग्रेटर, स्पष्ट या अंतर्निहित, उपस्थित नहीं है। चुकीं, इस विशिष्ट समस्या के लिए, हे विभाजन विधि का उपयोग करके उच्च-क्रम स्पष्ट गैर-कैनोनिकल सिम्प्लेक्टिक इंटीग्रेटर्स का परिवार बनाया जा सकता है।<ref>{{cite journal|last=He|first=Y.|author2=Qin, H.|author3=Sun, Y.|journal=Physics of Plasmas|date=2015|volume=22|title=वेलासोव-मैक्सवेल समीकरणों के लिए हैमिल्टनियन एकीकरण के तरीके|page=124503| doi=10.1063/1.4938034 |arxiv=1505.06076 |s2cid=118560512 }}</ref> <math display="inline">H</math> का विभाजन  4 भागों में, <math display="block">\begin{aligned}
H & = H_{x} + H_{y} + H_{z} + H_{\phi},\\
H & = H_{x} + H_{y} + H_{z} + H_{\phi},\\
H_{x} &= \frac{1}{2}v_{x}^{2},\ \ H_{y} = \frac{1}{2}v_{y}^{2},\ \ H_{z} = \frac{1}{2}v_{z}^{2},\ \ H_{\phi} = \phi,
H_{x} &= \frac{1}{2}v_{x}^{2},\ \ H_{y} = \frac{1}{2}v_{y}^{2},\ \ H_{z} = \frac{1}{2}v_{z}^{2},\ \ H_{\phi} = \phi,
\end{aligned}</math> हम गंभीर रूप से पाते हैं कि प्रत्येक उपप्रणाली के लिए, उदाहरण के लिए, <math display="block">i_{(\dot{\boldsymbol{x}},\dot{\boldsymbol{v}})}\Omega=-dH_{x}</math> और <math display="block">i_{(\dot{\boldsymbol{x}},\dot{\boldsymbol{v}})}\Omega=-dH_{\phi},</math> समाधान मानचित्र को स्पष्ट रूप से लिखा जा सकता है और सटीक गणना की जा सकती है। फिर विभिन्न रचनाओं का उपयोग करके स्पष्ट उच्च-क्रम गैर-विहित सिम्पलेक्टिक एल्गोरिदम का निर्माण किया जा सकता है। माना <math display="inline">\Theta_{x},\Theta_{y},\Theta_{z}</math> और <math display="inline">\Theta_{\phi}</math> 4 उपप्रणालियों के लिए सटीक समाधान मानचित्रों को निरूपित करें। एक प्रथम-क्रम सिम्पलेक्टिक योजना है <math display="block">\begin{aligned}
\end{aligned}</math> हम गंभीर रूप से पाते हैं कि प्रत्येक उपप्रणाली के लिए, उदाहरण के लिए, <math display="block">i_{(\dot{\boldsymbol{x}},\dot{\boldsymbol{v}})}\Omega=-dH_{x}</math> और <math display="block">i_{(\dot{\boldsymbol{x}},\dot{\boldsymbol{v}})}\Omega=-dH_{\phi},</math> समाधान मानचित्र को स्पष्ट रूप से लिखा जा सकता है और सटीक गणना की जा सकती है। फिर विभिन्न रचनाओं का उपयोग करके स्पष्ट उच्च-क्रम गैर-विहित सिम्पलेक्टिक एल्गोरिदम का निर्माण किया जा सकता है। माना <math display="inline">\Theta_{x},\Theta_{y},\Theta_{z}</math> और <math display="inline">\Theta_{\phi}</math> 4 उपप्रणालियों के लिए सटीक समाधान मानचित्रों को निरूपित करें। प्रथम-क्रम सिम्पलेक्टिक योजना है <math display="block">\begin{aligned}
\Theta_{1}\left(\Delta\tau\right)=\Theta_{x}\left(\Delta\tau\right)\Theta_{y}\left(\Delta\tau\right)\Theta_{z}\left(\Delta\tau\right)\Theta_{\phi}\left(\Delta\tau\right)~.\end{aligned}</math> एक सममित द्वितीय-क्रम सिम्पलेक्टिक योजना है, <math display="block">\begin{aligned}
\Theta_{1}\left(\Delta\tau\right)=\Theta_{x}\left(\Delta\tau\right)\Theta_{y}\left(\Delta\tau\right)\Theta_{z}\left(\Delta\tau\right)\Theta_{\phi}\left(\Delta\tau\right)~.\end{aligned}</math> सममित द्वितीय-क्रम सिम्पलेक्टिक योजना है, <math display="block">\begin{aligned}
\Theta_{2}\left(\Delta\tau\right) & =\Theta_{x}\left(\Delta\tau/2\right)\Theta_{y}\left(\Delta\tau/2\right)\Theta_{z}\left(\Delta\tau/2\right)\Theta_{\phi}\left(\Delta\tau\right)\\
\Theta_{2}\left(\Delta\tau\right) & =\Theta_{x}\left(\Delta\tau/2\right)\Theta_{y}\left(\Delta\tau/2\right)\Theta_{z}\left(\Delta\tau/2\right)\Theta_{\phi}\left(\Delta\tau\right)\\
  & \Theta_{z}\left(\Delta t/2\right)\Theta_{y}\left(\Delta t/2\right)\Theta_{x}\left(\Delta t/2\right)\!,
  & \Theta_{z}\left(\Delta t/2\right)\Theta_{y}\left(\Delta t/2\right)\Theta_{x}\left(\Delta t/2\right)\!,
\end{aligned}</math> जो कि एक अनुकूलित मॉडिफाइड स्ट्रैंग स्प्लिटिंग है। A <math display="inline">2(l+1)</math>-वें आदेश योजना का निर्माण a से किया जा सकता है <math display="inline">2l</math>-वें आदेश योजना ट्रिपल जंप की विधि का उपयोग कर, <math display="block">\begin{aligned}
\end{aligned}</math> जो कि अनुकूलित मॉडिफाइड स्ट्रैंग स्प्लिटिंग है। A <math display="inline">2(l+1)</math>-वें आदेश योजना का निर्माण a से किया जा सकता है <math display="inline">2l</math>-वें आदेश योजना ट्रिपल जंप की विधि का उपयोग कर, <math display="block">\begin{aligned}
\Theta_{2(l+1)}(\Delta\tau) & =\Theta_{2l}(\alpha_{l}\Delta\tau)\Theta_{2l}(\beta_{l}\Delta\tau)\Theta_{2l}(\alpha_{l}\Delta\tau)~,\\
\Theta_{2(l+1)}(\Delta\tau) & =\Theta_{2l}(\alpha_{l}\Delta\tau)\Theta_{2l}(\beta_{l}\Delta\tau)\Theta_{2l}(\alpha_{l}\Delta\tau)~,\\
\alpha_{l} & =1/(2-2^{1/(2l+1)})~,\\
\alpha_{l} & =1/(2-2^{1/(2l+1)})~,\\
\beta_{l} & =1-2\alpha_{l}~.
\beta_{l} & =1-2\alpha_{l}~.
\end{aligned}</math> हे विभाजन विधि संरचना-संरक्षण ज्यामितीय कण-इन-सेल (पीआईसी) एल्गोरिदम में उपयोग की जाने वाली प्रमुख तकनीकों में से एक है।<ref>{{cite journal|last=Xiao|first=J.|author2=Qin, H.|author3=Liu, J.|journal=Physics of Plasmas| date=2015|volume=22|title=Vlasov-Maxwell सिस्टम के लिए स्पष्ट उच्च-क्रम गैर-कैनोनिकल सिम्प्लेक्टिक पार्टिकल-इन-सेल एल्गोरिदम| issue=11|page=112504|doi=10.1063/1.4935904|arxiv=1510.06972|bibcode=2015PhPl...22k2504X |s2cid=12893515 }}</ref><ref>{{cite journal| last=Kraus| first=M|author2=Kormann, K|author3=Morrison, P.|author4=Sonnendrucker, E|journal=Journal of Plasma Physics |date=2017| volume=83| title=GEMPIC: geometric electromagnetic particle-in-cell methods| issue=4| page=905830401| doi=10.1017/S002237781700040X| arxiv=1609.03053| bibcode=2017JPlPh..83d9001K| s2cid=8207132}}</ref><ref>{{cite journal|last=Xiao|first=J.|author2=Qin, H.|author3=Liu, J.|journal=Plasma Science and Technology | date=2018|volume=20|title=वेलासोव-मैक्सवेल सिस्टम के लिए संरचना-संरक्षण ज्यामितीय कण-इन-सेल तरीके| issue=11| page=110501|doi=10.1088/2058-6272/aac3d1|arxiv=1804.08823|bibcode=2018PlST...20k0501X |s2cid=250801157 }}</ref><ref>{{cite journal| last=Glasser| first=A.|author2=Qin, H.|journal=Journal of Plasma Physics |date=2022| volume=88| title=परिमित तत्व बाहरी कलन का उपयोग करके कण-इन-सेल सिमुलेशन के लिए एक गेज-संगत हैमिल्टन विभाजन एल्गोरिथ्म| issue=2| page=835880202| doi=10.1017/S0022377822000290| arxiv=2110.10346| bibcode=2022JPlPh..88b8302G| s2cid=239049433}}</ref>
\end{aligned}</math> ही विभाजन विधि संरचना-संरक्षण ज्यामितीय कण-इन-सेल (पीआईसी) एल्गोरिदम में उपयोग की जाने वाली प्रमुख तकनीकों में से एक है।<ref>{{cite journal|last=Xiao|first=J.|author2=Qin, H.|author3=Liu, J.|journal=Physics of Plasmas| date=2015|volume=22|title=Vlasov-Maxwell सिस्टम के लिए स्पष्ट उच्च-क्रम गैर-कैनोनिकल सिम्प्लेक्टिक पार्टिकल-इन-सेल एल्गोरिदम| issue=11|page=112504|doi=10.1063/1.4935904|arxiv=1510.06972|bibcode=2015PhPl...22k2504X |s2cid=12893515 }}</ref><ref>{{cite journal| last=Kraus| first=M|author2=Kormann, K|author3=Morrison, P.|author4=Sonnendrucker, E|journal=Journal of Plasma Physics |date=2017| volume=83| title=GEMPIC: geometric electromagnetic particle-in-cell methods| issue=4| page=905830401| doi=10.1017/S002237781700040X| arxiv=1609.03053| bibcode=2017JPlPh..83d9001K| s2cid=8207132}}</ref><ref>{{cite journal|last=Xiao|first=J.|author2=Qin, H.|author3=Liu, J.|journal=Plasma Science and Technology | date=2018|volume=20|title=वेलासोव-मैक्सवेल सिस्टम के लिए संरचना-संरक्षण ज्यामितीय कण-इन-सेल तरीके| issue=11| page=110501|doi=10.1088/2058-6272/aac3d1|arxiv=1804.08823|bibcode=2018PlST...20k0501X |s2cid=250801157 }}</ref><ref>{{cite journal| last=Glasser| first=A.|author2=Qin, H.|journal=Journal of Plasma Physics |date=2022| volume=88| title=परिमित तत्व बाहरी कलन का उपयोग करके कण-इन-सेल सिमुलेशन के लिए एक गेज-संगत हैमिल्टन विभाजन एल्गोरिथ्म| issue=2| page=835880202| doi=10.1017/S0022377822000290| arxiv=2110.10346| bibcode=2022JPlPh..88b8302G| s2cid=239049433}}</ref>





Revision as of 12:56, 26 May 2023

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

परिचय

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

जहाँ स्थिति निर्देशांक को दर्शाता है, गति निर्देशांक, और हैमिल्टनियन है।

स्थिति और संवेग का सेट निर्देशांक विहित निर्देशांक कहलाते हैं। (अधिक पृष्ठभूमि के लिए हैमिल्टनियन यांत्रिकी देखें।)

हैमिल्टन के समीकरणों का समय विकास सिम्प्टोमोर्फिज़्म है, जिसका अर्थ है कि यह सिम्पलेक्टिक 2-प्रपत्र का संरक्षण करता है . संख्यात्मक योजना सिम्पलेक्टिक इंटीग्रेटर है यदि यह इस 2-फॉर्म को भी संरक्षित करती है।

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

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

सिम्पलेक्टिक एल्गोरिदम बनाने की विधियाँ

वियोज्य हैमिल्टन के लिए विभाजन विधियाँ

बंटवारे की विधियों से सिम्पलेक्टिक इंटीग्रेटर्स का व्यापक रूप से उपयोग किया जाने वाला वर्ग बनता है।

मान लें कि हैमिल्टन वियोज्य है, जिसका अर्थ है कि इसे रूप में लिखा जा सकता है

 

 

 

 

(1)

यह हैमिल्टनियन यांत्रिकी में अधिकांशतः होता है, जिसमें T गतिज ऊर्जा और V स्थितिज ऊर्जा है।

सांकेतिक सरलता के लिए, आइए हम प्रतीक का परिचय दें; विहित निर्देशांकों को निरूपित करने के लिए जिसमें स्थिति और संवेग दोनों निर्देशांक सम्मिलित हैं। फिर, परिचय में दिए गए हैमिल्टन के समीकरणों के सेट को एकल अभिव्यक्ति में व्यक्त किया जा सकता है

 

 

 

 

(2)

जहाँ पॉसों कोष्ठक है। इसके अतिरिक्त, ऑपरेटर का प्रारंभ करके , जो हेमिल्टनियन यांत्रिकी गणितीय औपचारिकता के साथ ऑपरेंड का पॉइसन ब्रैकेट लौटाता है, हैमिल्टन के समीकरण की अभिव्यक्ति को और सरल बनाया जा सकता है

समीकरणों के इस सेट का औपचारिक समाधान आव्यूह घातांक के रूप में दिया गया है:

 

 

 

 

(3)

की सकारात्मकता पर ध्यान दें आव्यूह एक्सपोनेंशियल में।

जब हैमिल्टनियन के पास समीकरण का रूप है (1), समाधान (3) के बराबर है

 

 

 

 

(4)

एसआई योजना समय-विकास ऑपरेटर का अनुमान लगाती है औपचारिक समाधान में (4) ऑपरेटरों के गुणांक के रूप में

 

 

 

 

(5)

जहाँ और वास्तविक संख्याएँ हैं, पूर्णांक है, जिसे इंटीग्रेटर का क्रम कहा जाता है, और जहाँ . ध्यान दें कि प्रत्येक ऑपरेटर और समानता प्रदान करता है, इसलिए उनका गुणांक दाईं ओर प्रदर्शित होता है (5) भी सिम्पलेक्टिक मानचित्र बनाता है।

तब से सभी के लिए , हम यह निष्कर्ष निकाल सकते हैं कि

 

 

 

 

(6)

टेलर श्रृंखला का उपयोग करके, के रूप में व्यक्त किया जा सकता है

 

 

 

 

(7)

जहाँ इच्छानुसार वास्तविक संख्या है। संयोजन (6) और (7), और के लिए समान तर्क का उपयोग करके जैसा कि हमने उपयोग किया है, हम पाते हैं

 

 

 

 

(8)

ठोस शब्दों में, मैपिंग देता है

और देता है

ध्यान दें कि ये दोनों मानचित्र व्यावहारिक रूप से संगणनीय हैं।

उदाहरण

समीकरणों का सरलीकृत रूप (निष्पादित क्रम में) हैं:

लाग्रंगियन निर्देशांक में परिवर्तित करने के बाद:

जहाँ पर बल सदिश है , त्वरण वेक्टर है , और द्रव्यमान की अदिश राशि है।

कई सिम्पलेक्टिक इंटीग्रेटर्स नीचे दिए गए हैं। उनका उपयोग करने का उदाहरण विधियाँ स्थिति के साथ कण पर विचार करना है और गति है

मूल्यों के साथ टाइमस्टेप प्रयुक्त करने के लिए कण के लिए, निम्न चरणों का पालन करें: क्रमिक रूप से:

  • स्थिति अद्यतन करें इसके (पहले अद्यतन) वेग को जोड़कर कण का से गुणा
  • वेग अद्यतन करें कण का इसमें त्वरण (अद्यतन स्थिति में) जोड़कर गुणा किया जाता है


पहले क्रम का उदाहरण

सिम्पलेक्टिक यूलर विधि प्रथम-क्रम समाकलक है और गुणांक

ध्यान दें कि समय-प्रतिवर्तीता की आवश्यकता होने पर उपरोक्त एल्गोरिदम काम नहीं करता है। एल्गोरिथ्म को दो भागों में प्रयुक्त किया जाना है, सकारात्मक समय चरणों के लिए नकारात्मक समय चरणों के लिए।

दूसरे क्रम का उदाहरण

वेरलेट इंटीग्रेशन दूसरे क्रम का इंटीग्रेटर है और गुणांक

तब से उपरोक्त एल्गोरिदम समय में सममित है। एल्गोरिथ्म के 3 चरण हैं, और चरण 1 और 3 बिल्कुल समान हैं, इसलिए सकारात्मक समय संस्करण का उपयोग नकारात्मक समय के लिए किया जा सकता है।

तीसरे क्रम का उदाहरण

तीसरा क्रम सिम्पलेक्टिक इंटीग्रेटर (के साथ ) की खोज 1983 में रोनाल्ड रुथ ने की थी।[1]

अनेक समाधानों में से एक द्वारा दिया गया है


चौथे क्रम का उदाहरण

चौथे क्रम के इंटीग्रेटर (के साथ ) भी 1983 में रूथ द्वारा खोजा गया था और उस समय कण-त्वरक समुदाय को सामान्यतः वितरित किया गया था। वन द्वारा जीवंत समीक्षा लेख में इसका वर्णन किया गया था।[2]

यह चौथा क्रम संपूर्नकर्ता 1990 में वन और रूथ द्वारा प्रकाशित किया गया था और उसी समय के आसपास दो अन्य समूहों द्वारा स्वतंत्र रूप से खोजा गया था।[3][4][5]

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

सामान्य अविभाज्य हैमिल्टनियन के लिए विभाजन विधियाँ

सामान्य अविभाज्य हैमिल्टनियन भी स्पष्ट और सिम्पलेक्टिक रूप से एकीकृत हो सकते हैं।

ऐसा करने के लिए, ताओ ने अवरोध प्रस्तुत किया जो इस तरह के प्रणाली के स्पष्ट विभाजन को सक्षम करने के लिए चरण स्थान की दो प्रतियों को साथ बांधता है।[7]


इसके अतिरिक्त विचार है , अनुकरण करता है , जिसका समाधान इससे सहमत है इस अर्थ में कि .

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

अनुप्रयोग

प्लाज्मा भौतिकी में

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

जिसकी -निर्भरता और निर्भरता वियोज्य नहीं हैं, और मानक स्पष्ट सिम्पलेक्टिक विधियों प्रयुक्त नहीं होते हैं। बड़े पैमाने पर समांतर समूहों पर बड़े पैमाने पर सिमुलेशन के लिए, चुकीं, स्पष्ट तरीकों को प्राथमिकता दी जाती है। इस कठिनाई को दूर करने के लिए, हम विशिष्ट विधियों का पता लगा सकते हैं -निर्भरता और -निर्भरता इस हैमिल्टनियन में उलझी हुई है, और केवल इस या इस प्रकार की समस्या के लिए सिम्पलेक्टिक एल्गोरिथम डिजाइन करने का प्रयास करें। सबसे पहले, हम ध्यान दें कि -निर्भरता द्विघात है, इसलिए पहले क्रम की सहानुभूति यूलर विधि में निहित है वास्तव में मुखर है। यह वही है जो कैनोनिकल सिम्प्लेक्टिक कण-इन-सेल (पीआईसी) एल्गोरिथम में उपयोग किया जाता है।[9] उच्च क्रम स्पष्ट विधियों का निर्माण करने के लिए, हम आगे ध्यान दें कि -निर्भरता और -इसमें निर्भरता गुणांक-वियोज्य हैं, दूसरे और तीसरे क्रम के स्पष्ट सिम्पलेक्टिक एल्गोरिदम का निर्माण कार्यों का उपयोग करके किया जा सकता है,[10] और समय-निर्भर विद्युत चुम्बकीय क्षेत्रों के लिए मनमाने ढंग से उच्च-क्रम वाले स्पष्ट सिम्पलेक्टिक इंटीग्रेटर्स का निर्माण रनगे-कुट्टा तकनीकों का उपयोग करके भी किया जा सकता है।[11] समस्या के निम्नलिखित गैर-विहित सहानुभूति संरचना को देखने के लिए अधिक सुरुचिपूर्ण और बहुमुखी विकल्प है,
यहाँ गैर-निरंतर गैर-विहित सिम्पलेक्टिक रूप है। गैर-निरंतर गैर-कैनोनिकल सिम्प्लेक्टिक संरचना के लिए सामान्य सिम्प्लेक्टिक इंटीग्रेटर, स्पष्ट या अंतर्निहित, उपस्थित नहीं है। चुकीं, इस विशिष्ट समस्या के लिए, हे विभाजन विधि का उपयोग करके उच्च-क्रम स्पष्ट गैर-कैनोनिकल सिम्प्लेक्टिक इंटीग्रेटर्स का परिवार बनाया जा सकता है।[12] का विभाजन 4 भागों में,
हम गंभीर रूप से पाते हैं कि प्रत्येक उपप्रणाली के लिए, उदाहरण के लिए,
और
समाधान मानचित्र को स्पष्ट रूप से लिखा जा सकता है और सटीक गणना की जा सकती है। फिर विभिन्न रचनाओं का उपयोग करके स्पष्ट उच्च-क्रम गैर-विहित सिम्पलेक्टिक एल्गोरिदम का निर्माण किया जा सकता है। माना और 4 उपप्रणालियों के लिए सटीक समाधान मानचित्रों को निरूपित करें। प्रथम-क्रम सिम्पलेक्टिक योजना है
सममित द्वितीय-क्रम सिम्पलेक्टिक योजना है,
जो कि अनुकूलित मॉडिफाइड स्ट्रैंग स्प्लिटिंग है। A -वें आदेश योजना का निर्माण a से किया जा सकता है -वें आदेश योजना ट्रिपल जंप की विधि का उपयोग कर,
ही विभाजन विधि संरचना-संरक्षण ज्यामितीय कण-इन-सेल (पीआईसी) एल्गोरिदम में उपयोग की जाने वाली प्रमुख तकनीकों में से एक है।[13][14][15][16]


यह भी देखें

संदर्भ

  1. Ruth, Ronald D. (August 1983). "A Canonical Integration Technique". IEEE Transactions on Nuclear Science. NS-30 (4): 2669–2671. Bibcode:1983ITNS...30.2669R. doi:10.1109/TNS.1983.4332919. S2CID 5911358.
  2. Forest, Etienne (2006). "Geometric Integration for Particle Accelerators". J. Phys. A: Math. Gen. 39 (19): 5321–5377. Bibcode:2006JPhA...39.5321F. doi:10.1088/0305-4470/39/19/S03.
  3. Forest, E.; Ruth, Ronald D. (1990). "Fourth-order symplectic integration" (PDF). Physica D. 43: 105–117. Bibcode:1990PhyD...43..105F. doi:10.1016/0167-2789(90)90019-L.
  4. Yoshida, H. (1990). "Construction of higher order symplectic integrators". Phys. Lett. A. 150 (5–7): 262–268. Bibcode:1990PhLA..150..262Y. doi:10.1016/0375-9601(90)90092-3.
  5. Candy, J.; Rozmus, W (1991). "A Symplectic Integration Algorithm for Separable Hamiltonian Functions". J. Comput. Phys. 92 (1): 230–256. Bibcode:1991JCoPh..92..230C. doi:10.1016/0021-9991(91)90299-Z.
  6. Blanes, S.; Moan, P. C. (May 2002). "Practical symplectic partitioned Runge–Kutta and Runge–Kutta–Nyström methods". Journal of Computational and Applied Mathematics. 142 (2): 313–330. Bibcode:2002JCoAM.142..313B. doi:10.1016/S0377-0427(01)00492-7.
  7. Tao, Molei (2016). "Explicit symplectic approximation of nonseparable Hamiltonians: Algorithm and long time performance". Phys. Rev. E. 94 (4): 043303. arXiv:1609.02212. Bibcode:2016PhRvE..94d3303T. doi:10.1103/PhysRevE.94.043303. PMID 27841574. S2CID 41468935.
  8. Qin, H.; Guan,X. (2008). "सामान्य चुंबकीय क्षेत्रों में लंबे समय के सिमुलेशन के लिए आवेशित कणों के मार्गदर्शक केंद्र गति के लिए एक परिवर्तनशील सहानुभूतिपूर्ण इंटीग्रेटर" (PDF). Physical Review Letters. 100 (3): 035006. doi:10.1103/PhysRevLett.100.035006. PMID 18232993.
  9. Qin, H.; Liu, J.; Xiao,J. (2016). "Canonical symplectic particle-in-cell method for long-term large-scale simulations of the Vlasov–Maxwell equations". Nuclear Fusion. 56 (1): 014001. arXiv:1503.08334. Bibcode:2016NucFu..56a4001Q. doi:10.1088/0029-5515/56/1/014001. S2CID 29190330.
  10. Zhang, R.; Qin, H.; Tang, Y. (2016). "चार्ज कण गतिकी के लिए कार्यों को उत्पन्न करने के आधार पर स्पष्ट सहानुभूतिपूर्ण एल्गोरिदम". Physical Review E. 94 (1): 013205. arXiv:1604.02787. Bibcode:2016PhRvE..94a3205Z. doi:10.1103/PhysRevE.94.013205. PMID 27575228. S2CID 2166879.
  11. Tao, M. (2016). "सामान्य विद्युतचुंबकीय क्षेत्रों में आवेशित कणों के लिए स्पष्ट उच्च-आदेश सहानुभूतिपूर्ण इंटीग्रेटर्स". Journal of Computational Physics. 327: 245. arXiv:1605.01458. Bibcode:2016JCoPh.327..245T. doi:10.1016/j.jcp.2016.09.047. S2CID 31262651.
  12. He, Y.; Qin, H.; Sun, Y. (2015). "वेलासोव-मैक्सवेल समीकरणों के लिए हैमिल्टनियन एकीकरण के तरीके". Physics of Plasmas. 22: 124503. arXiv:1505.06076. doi:10.1063/1.4938034. S2CID 118560512.
  13. Xiao, J.; Qin, H.; Liu, J. (2015). "Vlasov-Maxwell सिस्टम के लिए स्पष्ट उच्च-क्रम गैर-कैनोनिकल सिम्प्लेक्टिक पार्टिकल-इन-सेल एल्गोरिदम". Physics of Plasmas. 22 (11): 112504. arXiv:1510.06972. Bibcode:2015PhPl...22k2504X. doi:10.1063/1.4935904. S2CID 12893515.
  14. Kraus, M; Kormann, K; Morrison, P.; Sonnendrucker, E (2017). "GEMPIC: geometric electromagnetic particle-in-cell methods". Journal of Plasma Physics. 83 (4): 905830401. arXiv:1609.03053. Bibcode:2017JPlPh..83d9001K. doi:10.1017/S002237781700040X. S2CID 8207132.
  15. Xiao, J.; Qin, H.; Liu, J. (2018). "वेलासोव-मैक्सवेल सिस्टम के लिए संरचना-संरक्षण ज्यामितीय कण-इन-सेल तरीके". Plasma Science and Technology. 20 (11): 110501. arXiv:1804.08823. Bibcode:2018PlST...20k0501X. doi:10.1088/2058-6272/aac3d1. S2CID 250801157.
  16. Glasser, A.; Qin, H. (2022). "परिमित तत्व बाहरी कलन का उपयोग करके कण-इन-सेल सिमुलेशन के लिए एक गेज-संगत हैमिल्टन विभाजन एल्गोरिथ्म". Journal of Plasma Physics. 88 (2): 835880202. arXiv:2110.10346. Bibcode:2022JPlPh..88b8302G. doi:10.1017/S0022377822000290. S2CID 239049433.
  • Leimkuhler, Ben; Reich, Sebastian (2005). Simulating Hamiltonian Dynamics. Cambridge University Press. ISBN 0-521-77290-7.
  • Hairer, Ernst; Lubich, Christian; Wanner, Gerhard (2006). Geometric Numerical Integration: Structure-Preserving Algorithms for Ordinary Differential Equations (2 ed.). Springer. ISBN 978-3-540-30663-4.
  • Kang, Feng; Qin, Mengzhao (2010). Symplectic geometric algorithms for Hamiltonian systems. Springer.