A priori error analysis of a mass-lumped midpoint finite element method with a structure-preserving solver for the Landau--Lifshitz--Gilbert equation with Dzyaloshinskii--Moriya interaction

तकनीकी सारांश: बल्क DMI के साथ LLG समीकरण के लिए एक मास-लम्पड मिडपॉइंट परिमित तत्व विधि का ए प्रियोरी त्रुटि विश्लेषण

1. समस्या विवरण
यह शोध पत्र फेरोमैग्नेटिक सामग्रियों में चुंबकत्व गतिशीलता (magnetisation dynamics) को नियंत्रित करने वाले लैंडौ-लिफ़शिट्ज़-गिल्बर्ट (LLG) समीकरण के संख्यात्मक सन्निकटन (numerical approximation) को संबोधित करता है। विशेष रूप से, यह अध्ययन बल्क डज़ालोशिंस्की-मोरिया इंटरैक्शन (DMI) के समावेश पर केंद्रित है, जो प्रभावी क्षेत्र (effective field) में एक प्रथम-क्रम पद पेश करता है और एक स्वाभाविक काइरल सीमा स्थिति (chiral boundary condition) उत्पन्न करता है। मुख्य गणितीय चुनौतियाँ हैं: बिंदुवार इकाई-लंबाई प्रतिबंध ($|m|=1$), समाधान स्थान की गैर-उत्तल ज्यामिति (non-convex geometry), और समीकरण की अत्यधिक गैर-रैखिक संरचना। जबकि केवल एक्सचेंज मॉडल का व्यापक अध्ययन किया जा चुका है, DMI का समावेश वेरिएशनल संरचना और सीमा व्यवहार को संशोधित करता है, जिससे स्थिर, सटीक और संरचना-संरक्षण (structure-preserving) परिमित तत्व योजनाओं का निर्माण अत्यंत सूक्ष्म हो जाता है।

2. कार्यप्रणाली
लेखक एक पूर्णतः विविक्त (fully discrete), मास-लम्पड मिडपॉइंट परिमित तत्व योजना प्रस्तावित करते हैं और उसका विश्लेषण करते हैं। इस कार्यप्रणाली में निम्नलिखित घटक शामिल हैं:

• विविक्तकरण (Discretization): स्थानिक डोमेन को निम्नतम-क्रम के अनुरूप लैग्रेंज परिमित तत्वों ($P_1$) का उपयोग करके विविक्त किया गया है। समय एकीकरण के लिए निहित मिडपॉइंट नियम (implicit midpoint rule) का उपयोग किया जाता है।
• मास-लम्पिंग (Mass-Lumping): एक न्यूनीकरण-एकीकृत (reduced-integration/mass-lumped) आंतरिक उत्पाद (inner product) का उपयोग किया जाता है। यह तकनीक महत्वपूर्ण है क्योंकि यह प्रत्येक नोड पर नोडल इकाई-लंबाई प्रतिबंध को सटीक रूप से लागू करने की अनुमति देती है, एक ऐसा गुण जिसे मानक $L^2$-आंतरिक उत्पाद अक्सर गैर-रैखिक योजनाओं में बनाए रखने में विफल रहते हैं।
• प्रभावी क्षेत्र ऑपरेटर (Effective Field Operator): विविक्त प्रभावी क्षेत्र में एक्सचेंज स्टिफनेस और बल्क DMI दोनों शामिल हैं। DMI पद को एक विविक्त द्विपद रूप (bilinear form) के माध्यम से संभाला जाता है जो स्वाभाविक रूप से काइरल सीमा स्थिति का सम्मान करता है। विश्लेषण को सुगम बनाने के लिए, एक शिफ्टेड प्रभावी क्षेत्र ऑपरेटर पेश किया गया है ताकि विविक्त द्विपद रूप की समान कोएर्सिविटी (uniform coercivity) सुनिश्चित की जा सके।
• सॉल्वर रणनीति (Solver Strategy): चूंकि मिडपॉइंट योजना प्रत्येक समय चरण में एक गैर-रैखिक बीजगणितीय प्रणाली का परिणाम देती है, इसलिए यह शोध पत्र एक संरचना-संरक्षण फिक्स्ड-पॉइंट इटरेशन प्रस्तावित करता है। सामान्य सॉल्वरों के विपरीत, यह इटरेशन प्रत्येक इटरेशन पर नोडल इकाई-लंबाई प्रतिबंध को बनाए रखने के लिए डिज़ाइन किया गया है।
• त्रुटि विश्लेषण ढांचा (Error Analysis Framework): विश्लेषण सटीक समाधान के एलिप्टिक प्रोजेक्शन के विरुद्ध संख्यात्मक समाधान की तुलना करता है। प्रमाण मिडपॉइंट विधि में अंतर्निहित ज्यामितीय निरस्तीकरण (geometric cancellations), प्रथम-क्रम DMI योगदान के सावधानीपूर्वक अनुमान, और एक विविक्त ऊर्जा नियम से प्राप्त स्थिरता अनुमानों पर आधारित है।

3. मुख्य योगदान
यह शोध पत्र कम्प्यूटेशनल माइक्रोजेMagnetics के साहित्य में तीन प्राथमिक योगदान देता है:

1. फॉर्मुलेशन और स्थिरता: लेखक एक्सचेंज और बल्क DMI के साथ LLG समीकरण के लिए एक पूर्णतः विविक्त मास-लम्पड मिडपॉइंट योजना तैयार करते हैं। वे सिद्ध करते हैं कि यह योजना नोडल इकाई-लंबाई प्रतिबंध को सटीक रूप से संरक्षित करती है और एक विविक्त ऊर्जा पहचान (energy identity) को संतुष्ट करती है, जो अन्य संरचना-संरक्षण विधियों में पाए जाने वाले सख्त CFL-प्रकार के प्रतिबंधों या ज्यामितिक मेश प्रतिबंधों (जैसे न्यून कोण) की आवश्यकता के बिना स्थिरता की नींव प्रदान करती है।
2. इष्टतम ए प्रियोरी त्रुटि अनुमान (आदर्श योजना): उपयुक्त नियमितता धारणाओं के तहत, यह पत्र आदर्श योजना (मान लीजिए कि गैर-रैखिक प्रणाली का सटीक समाधान प्राप्त है) के लिए एक इष्टतम-क्रम अभिसरण अनुमान सिद्ध करता है। त्रुटि ऊर्जा नॉर्म ($H^1$) में $O(h + k^2)$ के क्रम के साथ अभset होती है, जो स्थान में प्रथम-क्रम और समय में द्वितीय-क्रम की सटीकता को दर्शाती है। लेखकों के ज्ञान तक, मास-लम्पड मिडपॉइंट विधि के लिए, यहाँ तक कि DMI की अनुपस्थिति में भी, यह पहला ऐसा ए प्रियोरी त्रुटि अनुमान है।
3. अपूर्ण सॉल्वर (Inexact Solvers) का विश्लेषण: यह पहचानते हुए कि व्यवहार में गैर-रैखिक प्रणाली को पुनरावृत्ति (iteratively) द्वारा हल किया जाना चाहिए, यह पत्र सॉल्वर त्रुटियों के प्रभाव का विश्लेषण करता है। यह सिद्ध करता है कि यदि फिक्स्ड-पॉइंट इटरेशन को एक विशिष्ट अवशेष (residual) स्टॉपिंग मानदंड के आधार पर समाप्त किया जाता है, तो परिणामी अपूर्ण योजना आदर्श योजना के समान ही अभिसरण दर बनाए रखती है, बशर्ते कि सॉल्वर टॉलरेंस विविक्त मापदंडों के सापेक्ष पर्याप्त रूप से छोटा हो।

4. परिणाम

• सैद्धांतिक सीमाएँ (Theoretical Bounds): मुख्य सैद्धांतिक परिणाम (थ्योरम 3.7 और थ्योरम 4.8) यह स्थापित करता है कि ऊर्जा नॉर्म में त्रुटि $C(h + k^2)$ और सॉल्वर टॉलरेंस पर निर्भर एक पद द्वारा सीमित है। विशेष रूप रूप से, एक समान सॉल्वर टॉलरेंस $\varepsilon_{fp}$ के लिए, स्थिति $\varepsilon_{fp} = O(h)$ इष्टतम अभिसरण दर बनाए रखने के लिए पर्याप्त है।
• संख्यात्मक प्रयोग:
  • अभिसरण दरें: एक सुचारू समाधान पर संख्यात्मक परीक्षण सैद्धांतिक अभिसरण दरों की पुष्टि करते हैं: स्थानिक परिशोधन (refinement) के लिए स्थानिक रूप से $H^1$ में लगभग प्रथम-क्रम और $L^2$ में द्वितीय-क्रम, और सामयिक परिशोधन के लिए द्वितीय-क्रम।
  • सॉल्वर टॉलरेंस: फिक्स्ड-पॉइंट टॉलरेंस को बदलते हुए प्रयोग प्रदर्शित करते हैं कि त्रुटि $\varepsilon_{fp}/h$ के साथ रैखिक रूप से स्केल करती है, जो सॉल्वर-प्रेरित त्रुटियों के संबंध में सैद्धांतिक भविष्यवाणी को मान्य करती है।
  • संरचना संरक्षण: एक ब्लॉक-प्रकार के स्काईर्मियन (Bloch-type skyrmion) और एक उच्च-ग्रेडिएंट "बबलिंग" प्रयोग वाले सिमुलेशन पुष्टि करते हैं कि विधि मशीन परिशुद्धता तक नोडल इकाई-लंबाई प्रतिबंध को संरक्षित करती है और सही ढंग से विसरित ऊर्जा व्यवहार (dissipative energy behavior) को कैप्चर करती है।
  • DMI प्रभाव: प्रयोग DMI चिह्न ($\kappa_d$) के असममित प्रभाव को दर्शाते हैं। नकारात्मक $\kappa_d$ ग्रेडिएंट एकाग्रता (बबल फॉर्मेशन) को तेज करता है, जबकि सकारात्मक $\kappa_d$ एकाग्रता को दबाने और विश्रांति (relaxation) को बढ़ावा देने की प्रवृत्ति रखता है, जो काइरल इंटरैक्शन के अनुरूप है।

5. महत्व और दावे
यह शोध पत्र दावा करता है कि यह बल्क DMI के साथ LLG समीकरण पर लागू मास-लम्पड मिडपॉइंट परिमित तत्व विधि के लिए पहला कठोर ए प्रियोरी त्रुटि विश्लेषण प्रदान करता है। इसका महत्व इस तथ्य में निहित है कि यह इस विधि के ज्ञात स्थिरता गुणों (सटीक प्रतिबंध संरक्षण और ऊर्जा अपव्यय) और मात्रात्मक अभिसरण दरों के बीच के अंतर को पाटता है।

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