Trefftz DG Approximation of the T-Matrix for Scattering by Periodic Layered Structures

तकनीकी सारांश: आवधिक स्तरित संरचनाओं द्वारा प्रकीर्णन के लिए ट्रेफ़्त्ज़ डीजी (Trefftz DG) टी-मैट्रिक्स सन्निकटन

समस्या विवरण
यह शोध पत्र आवधिक स्तरित ग्रेटिंग्स (periodic layered gratings) द्वारा समय-हार्मोनिक विद्युत चुम्बकीय तरंग प्रकीर्णन के संख्यात्मक सिमुलेशन को संबोधित करता है, जिसे द्वि-आयामी हेल्महोल्ट्ज़ समीकरण (Helmholtz equation) द्वारा मॉडल किया गया है। भौतिक डोमेन पारगम्य डाइइलेक्ट्रिक मीडिया और संभावित रूप से अभेद्य परम विद्युत चालक (PEC) बाधाओं वाले सीमित संख्या में स्टैक्ड परतों ($N$) से बना है। सापेक्ष पारगम्यता $\varepsilon$ खंडवार-स्थिर (piecewise-constant) है और क्षैतिज दिशा ($x_1$) में आवधिकता $L$ के साथ है। प्रकीर्णन समस्या को एक एकल आवधिक सेल पर अर्ध-आवधिक (quasi-periodic) सीमा स्थितियों और ऊर्ध्वाधर दिशाओं में विकिरण स्थितियों को लागू करने के लिए डिरिचलेट-टू-नेमन (DtN) ऑपरेटरों का उपयोग करते हुए तैयार किया गया है।

मुख्य चुनौती मानक मोनोलिथिक (monolithic) विधियों का उपयोग करके बड़े, जटिल बहु-परत डोमेन को विविक्त (discretize) करने से जुड़ी कम्प्यूटेशनल लागत है। जैसे-जैसे परतों की संख्या बढ़ती है, वैश्विक रैखिक प्रणाली का आकार आनुपातिक रूप से बढ़ता है, जिससे महत्वपूर्ण कम्प्यूटेशनल बोझ उत्पन्न होता है।

कार्यप्रणाली
लेखक एक हाइब्रिड संख्यात्मक ढांचे का प्रस्ताव करते हैं जो मल्टीपल-स्कैटरिंग थ्योरी से टी-मैट्रिक्स (T-matrix) विधि को ट्रेफ़्त्ज़ डिस्कंटीन्यूअस गैलेरकिन (TDG) विधि के साथ जोड़ता है। यह दृष्टिकोण वैश्विक प्रकीर्णन समस्या को छोटी, स्वतंत्र स्थानीय समस्याओं के अनुक्रम में विभाजित करता है।

1. डोमेन अपघटन और टी-मैट्रिक्स सूत्रीकरण:
   वैश्विक डोमेन $\Omega$ को $N$ गैर-अतिव्यापी परतों $\Omega_j$ में विभाजित किया गया है। प्रत्येक परत के लिए, एक स्थानीय सीमा मान समस्या (BVP) को परिभाषित किया जाता है जहाँ आपतित क्षेत्र (incident field) अर्ध-आवधिक मोड विस्तार (ऊपर और नीचे की ओर प्रसारित तरंगों) से बना होता है। इस विधि का मूल तत्व प्रत्येक परत के लिए एक टी-मैट्रिक्स ($T_j$) की गणना करना है। यह मैट्रिक्स आने वाले तरंग मोड के फूरियर गुणांकों को बाहर जाने वाले (प्रकीर्णित और संचारित) तरंग मोड के गुणांकों से मैप करता है।

• टी-मैट्रिक्स का व्युत्पन्न अर्ध-आवधिकता द्वारा निर्देशित प्राकृतिक फूरियर विस्तार का उपयोग करके किया जाता है।
• वैश्विक समाधान को परतों के बीच ऊर्ध्वाधर दूरी को ध्यान में रखने वाले प्रोपेगेशन मैट्रिसेस के माध्यम से स्थानीय टी-मैट्रिक्स को जोड़कर पुनर्गठित किया जाता है। इसके परिणामस्वरूप एक वैश्विक रैखिक प्रणाली प्राप्त होती है जो बैंडेड और ब्लॉक-ट्राइडिगोनल (block-tridiagonal) होती है, जिसे बैंड गॉसियन एलिमिनेशन का उपयोग करके $O(N)$ रैखिक जटिलता के साथ हल किया जा सकता है।

2. DtN-TDG के माध्यम से स्थानीय सन्निकटन:
   प्रत्येक परत $\Omega_j$ के लिए टी-मैट्रिक्स के प्रविष्टियों (entries) की गणना करने के लिए, लेखक DtN-TDG विधि का उपयोग करते हैं।

• विविक्तीकरण (Discretization): परत को एक उत्तल-बहुभुज (convex-polygonal) परिमित तत्व विभाजन के साथ मेष (mesh) किया गया है। समाधान स्थान एक ट्रेफ़्त्ज़ स्पेस है, जहाँ आधार फलन (basis functions) स्थानीय रूप से हेल्महोल्ट्ज़ समीकरण को संतुष्ट करते हैं।
• आधार फलन (Basis Functions): प्लेन वेव्स (plane waves) का उपयोग आधार कार्यों के रूप में किया जाता है। यह विकल्प स्टिफनेस मैट्रिक्स प्रविष्टियों के स्पष्ट विश्लेषणात्मक अभिव्यक्तियों की अनुमति देता है, जिससे संख्यात्मक क्वाड्रैचर त्रुटियों से बचा जा सके।
• सीमा स्थितियाँ: विकिरण स्थिति और अर्ध-आवधिकता को लागू करने के लिए क्षैतिज इंटरफेस ($\Gamma^\pm_j$) पर DtN ऑपरेटर लागू किए जाते हैं। टी-मैट्रिक्स के कॉलम को भरने के लिए इस सूत्र को मनमाने अर्ध-आवधिक आपतित क्षेत्रों को संभालने के लिए विस्तारित किया गया है।
• ट्रंकेशन (Truncation): अनंत टी-मैट्रिक्स और DtN ऑपरेटरों को $M$ पैरामीटर (रिटेन फूरियर मोड की संख्या) द्वारा निर्धारित सीमित आकार तक ट्रंकेट किया जाता है।

3. असेंबली और समाधान:
   एक परत के लिए टी-मैट्रिक्स को $4M+2$ स्थानीय BVPs (आधार आपतित मोड के अनुरूप) को हल करके असेंबल किया जाता है। इसके बाद, स्कैटरड फील्ड गुणांकों को हल करने के लिए इन स्थानीय टी-मैट्रिक्स को जोड़ा जाता है।

प्रमुख योगदान

• स्केल्स का पृथक्करण (Decoupling of Scales): सूत्रण वैश्विक युग्मन (global coupling) की जटिलता को स्थानीय विवक्तीकरण से अलग करता है। वैश्विक प्रणाली का आकार केवल परतों की संख्या ($N$) और ट्रंकेशन ऑर्डर ($M$) पर निर्भर करता है, न कि परतों के ऊर्ध्वाधर आकार, उनकी भौतिक विशेषताओं या आंतरिक मेश रिज़ॉल्यूशन पर। यह वैश्विक युग्मन चरण की लागत बढ़ाए बिना परतों के भीतर उच्च-क्रम या सूक्ष्म विवक्तीकरण की अनुमति देता है।
• पुन: प्रयोज्यता और समानांतरकरण (Reusability and Parallelization): यदि परतें समान हैं (या एक-दूसरे के प्रतिबिंब हैं), तो उनके टी-मैट्रिक्स समान होते हैं। यह संरचनाओं के लिए टी-मैट्रिक्स के पूर्व-गणना और पुन: उपयोग की अनुमति देता है जो दोहराए गए पैटर्न रखते हैं। विभिन्न परतों के लिए टी-मैट्रिक्स के संयोजन को समानांतर में किया जा सकता है।
• जटिल ज्यामिति को संभालना: अक्ष-संरेखित (axis-aligned) विशेषताओं तक सीमित विधियों के विपरीत, प्लेन वेव बेस के साथ TDG विधि परतों के भीतर घूर्णित समावेशन (rotated inclusions) और मनमाने बहुभुज आकारों, जिसमें अभेद्य (PEC) बाधाएं भी शामिल हैं, को संभाल सकती है।
• अर्ध-आवधिक क्षेत्रों का विस्तार: कार्य मौजूदा DtN-TDG ढांचे को मनमाने अर्ध-आवधिक आपतित क्षेत्रों को समायोजित करने के लिए विस्तारित करता है, जो स्तरित सेटिंग में टी-मैट्रिक्स के निर्माण के लिए एक पूर्व शर्त है।

संख्यात्मक परिणाम
लेखक कई संख्यात्मक प्रयोगों के माध्यम से अपनी विधि को मान्य करते हैं:

• अभिसरण (Convergence): समरूप सामग्रियों के बीच फ्लैट इंटरफेस वाले एक परीक्षण मामले में, विधि प्लेन वेव आधार कार्यों ($p$) की संख्या के संबंध में घातांकीय अभिसरण (exponential convergence) प्रदर्शित करती है, जो सैद्धांतिक अपेक्षाओं के अनुरूप है।
• विलक्षणता (Singularities): विधि जटिल सामग्री इंटरफेस वाले स्तरित संरचनाओं में कॉर्नर सिंगुलैरिटी को सफलतापूर्वक संभालती है। यह दिखाया गया है कि त्रुटि TDG विवक्तीकरण द्वारा हावी है, एक बार जब फूरियर ट्रंकेशन ऑर्डर $M$ एक निश्चित सीमा को पार कर जाता है।
• स्केलेबिलिटी (Scalability): पूर्ण डोमेन पर एक मोनोलिथिक DtN-TDG सॉल्वर के साथ प्रस्तावित टी-मैट्रिक्स दृष्टिकोण के बीच तुलना प्रस्तावित विधि की स्केलेबिलिटी को प्रदर्शित करती है। 45 परतों वाली संरचना के लिए, टी-मैट्रिक्स दृष्टिकोण लगभग स्थिर कम्प्यूटेशनल लागत बनाए रखता है (जो स्थानीय टी-मैट्रिक्स गणना द्वारा नियंत्रित है), जबकि मोनोलिथिक दृष्टिकोण बढ़ते वैश्विक रैखिक प्रणाली के आकार के कारण समय में तीव्र वृद्धि दिखाता है।
• बेंचमार्किंग (RCWA): विधि की तुलना एक डिस्ट्रीब्यूटेड ब्रैग रिफ्लेक्टर (DBR) संरचना पर रिगोरस कपल्ड-वेव एनालिसिस (RCWA) के विरुद्ध की गई है, जिसमें क्वैसी-बाउंड-स्टेट्स-इन-द-कंटीनम (qBIC) गाइडेड-मोड रेजोनेंस शामिल हैं। परिणाम बैंडगैप और ट्रांसमिटेंस ऑसिलेशन की भविष्यवाणी करने में उत्कृष्ट सहमति दिखाते हैं। टी-मैट्रिक्स TDG विधि RCWA के समान संकीर्ण रेजोनेंस को कैप्चर करती है, जिसमें रेजोनेंस पीक्स के बाहर सामान्य सापेक्ष त्रुटियां $10^{-2}$ से कम होती हैं।

महत्व और दावे
लेखक दावा करते हैं कि प्रस्तावित ढांचा आवधिक स्तरित संरचनाओं के लिए एक मजबूत और कुशल विकल्प प्रदान करता है। इसका प्राथमिक महत्व परतों की संख्या के संबंध में इसके रैखिक कम्प्यूटेशनल जटिलता में निहित है, जो इसे विशेष रूप से बहु-परत विन्यास के लिए फायदेमंद बनाता है, खासकर दोहराए गए ज्यामिति वाले मामलों में।

लेखक इस बात पर जोर देते हैं कि जबकि RCWA एक स्थापित और अक्सर मानक आवधिक संरचनाओं के लिए तेज़ विधि है, टी-मैट्रिक्स TDG दृष्टिकोण अधिक व्यापकता (generality) प्रदान करता है। विशेष रूप से, यह निम्नलिखित को संभाल सकता है:

1. घूर्णित समावेशन (Rotated inclusions) जिन्हें समन्वय-संरेखित (coordinate-aligned) मेश की आवश्यकता नहीं होती।
2. अभेद्य (PEC) बाधाएं, जिन्हें मानक RCWA सूत्रीकरणों के साथ संभालना कठिन या असंभव है।

यह कार्य एक व्यावहारिक कार्यान्वयन (ऑनलाइन उपलब्ध) प्रदान करता है जो उच्च-गुणवत्ता-कारक रेजोनेंस (high-quality-factor resonances) सहित जटिल फोटोनिक संरचनाओं के लिए सटीकता और दक्षता के बीच सफलतापूर्वक संतुलन बनाता है। वे रेजोनेंस पीक्स (qBIC-GMRs) पर अपने मेथड के प्रदर्शन के संबंध में विनम्र रहते हैं, यह नोट करते हुए कि इन बिंदुओं पर BVP समाधान की गैर-विशिष्टता (non-uniqueness) अभिसरण संबंधी समस्याओं का कारण बन सकती है, जो कि अन्य संख्यात्मक विधियों के साथ साझा की गई एक ज्ञात सीमा है।