Finite-Slab Reflectance and Transmittance for Henyey-Greenstein Scattering via a First-Passage Transfer Operator

तकनीकी सारांश: प्रथम-प्रवेश स्थानांतरण ऑपरेटर (First-Passage Transfer Operator) के माध्यम से हेनी-ग्रीनस्टीन स्कैटरिंग हेतु परिमित-स्लैब परावर्तन और संचरण

समस्या विवरण
यह शोध पत्र हेनी-ग्रीनस्टीन (HG) स्कैटरिंग वाले एक प्लेन-पैरेलल स्लैब के लिए परावर्तन ($R$) और संचरण ($T$) की गणना को संबोधित करता है। यह प्रणाली विषमता पैरामीटर (asymmetry parameter) $g$, एकल-स्कैटरिंग एल्बेडो (single-scattering albedo) $a$, और आपतन कोण (incidence angle) $\theta_0$ द्वारा परिभाषित है। जबकि मानक दृष्टिकोण मोंटे कार्लो (MC) सैंपलिंग या डिस्क्रीट-ऑर्डिनेट/H-फंक्शन मशीनरी पर निर्भर करते हैं, यह कार्य एक ऐसे विधि की खोज करता है जो अपने स्वरूप में सटीक है (स्वयं ट्रांसपोर्ट मॉडल से परे), यादृच्छिक सैंपलिंग से मुक्त है, और एकल ऑपरेटर मूल्यांकन से संपूर्ण एल्बेडो अक्ष के लिए परिणाम देने में सक्षम है।

कार्यप्रणाली: प्रथम-प्रवेश स्थानांतरण ऑपरेटर (The First-Passage Transfer Operator)
विधि का मूल एक प्रथम-प्रवेश स्थानांतरण ऑपरेटर है जो गहराई-दिशा अवस्था (depth–direction state) $(z, \mu)$ पर कार्य करता है, जहाँ $z$ ऑप्टिकल थिकनेस (optical depth) है और $\mu$ दिशा कोसाइन (direction cosine) है। यह दृष्टिकोण पूर्ण 3D रैंडम वॉक को $(z, \mu)$ में एक बंद मार्कोव चेन (Markov chain) में कम कर देता है, बिना अज़ीमुथल सहसंबंधों (azimuthal correlations) को त्याग दिए या डिफ्यूजन सन्निकटन (diffusion approximations) का आह्वान किए।

1. अवस्था न्यूनीकरण (State Reduction): लेखक प्रदर्शित करते हैं कि पलायन की संभावना (escape probability) केवल गहराई चरणों $\{z_k\}$ और दिशा कोसाइन $\{\mu_k\}$ के अनुक्रम पर निर्भर करती है। अज़ीमुथल कोण विकसित होता है लेकिन उन ट्रांसवर्स निर्देशांकों को प्रभावित नहीं करता जो प्लेन-पैरेलल पलायन के लिए प्रासंगिक हैं, जिससे इस प्रक्रिया को केवल $(z, \mu)$ में एक मार्कोव चेन के रूप में मॉडल किया जा सकता है।
2. ऑपरेटर निर्माण:
   • मुक्त उड़ान (Free Flight): फोटॉन निश्चित $\mu$ पर गहराई में मुक्त उड़ान चलते हैं, जो पथ लंबाई के घातांकीय वितरण (exponential distribution) द्वारा नियंत्रित होते हैं।
   • कोणीय पुनर्वितरण (Angular Redistribution): टकराव होने पर, दिशा को HG फेज फंक्शन के माध्यम से अपडेट किया जाता है। यूनिफॉर्म स्कैटरिंग अज़ीमुथ पर औसत निकालकर, लेखक एक पुनर्वितरण कर्नेल $\Psi(\mu \to \mu')$ प्राप्त करते हैं जो केवल आने वाले और जाने वाले दिशा कोसाइन पर निर्भर करता है।
   • अवशोषक सीमाएँ (Absorbing Boundaries): ऑपरेटर को स्लैब $(0, \tau)$ के भीतर सीमित किया गया है जिसमें अवशोषक सीमाएँ हैं। सीमा तक पहुँचने वाले द्रव्यमान (mass) को हटा दिया जाता है, जो ठीक $n$ टकरावों के बाद बाहर निकलने वाले फोटॉनों के क्रम-अनुरूप परावर्तन ($P_R(n)$) और संचरण ($P_T(n)$) नियमों को परिभाषित करता है।
3. एल्बेडो वेटिंग (Albedo Weighting): किसी भी एल्बेडो $a$ के लिए कुल परावर्तन और संचरण को क्रम-अनुरूप नियमों को वेटिंग देकर प्राप्त किया जाता है:
   $$R(\tau, g, a, \theta_0) = \sum_{n \ge 1} P_R(n) a^n$$
   $$T(\tau, g, a, \theta_0) = e^{-\tau/\mu_0} + \sum_{n \ge 1} P_T(n) a^n$$
   यह एकल डिस्क्रीटाइज्ड ऑपरेटर मूल्यांकन से सभी $a \in [0, 1]$ के लिए $R$ और $T$ की गणना करने की अनुमति देता है।
4. संख्यात्मक कार्यान्वयन: विधि एक समान गहराई ग्रिड और कोणों के लिए गॉस-लेजेंड्रे क्वाड्रचर (Gauss–Legendre quadrature) का उपयोग करती है। ऑपरेटर को एक कटऑफ ऑर्डर $n_{max}$ तक इटरेट किया जाता है जो कंजर्वेटिव स्लैब्स ($a=1$) के लिए द्रव्यमान को कैप्चर करने के लिए पर्याप्त है।

प्रमुख परिणाम और सत्यापन

• सटीकता: ऑपरेटर एक विस्तृत पैरामीटर स्पेस ($g \in [0, 0.95]$, $\tau \in [0.5, 16]$, $a \in [0.5, 1]$) में पूर्ण 3D मोंटे कार्लो परिणामों को $\le 1.6 \times 10^{-3}$ के पूर्ण त्रुटि के भीतर पुनरुत्पादित करता है। ऊर्जा संरक्षण ($R+T=1$) $a=1$ के लिए $< 10^{-4}$ तक संतुष्ट है।
• उभरते वितरण (Emergent Distributions): विधि बिना किसी अतिरिक्त गणना लागत के उभरते कोणीय वितरण $R(\mu)$ और $T(\mu)$ प्रदान करती है, जो परावर्तन में MC बिनिंग के कारण $<1\%$ त्रुटि और संचरण में $\sim 2\%$ त्रुटि (तीव्र फॉरवर्ड पीक्स के MC बिनिंग द्वारा सीमित) के साथ MC डेटा से मेल खाती है।
• मोटाई और अनिसोट्रॉपी (Thickness and Anisotropy): परिणाम पुष्टि करते हैं कि फॉरवर्ड स्कैटरिंग ($g > 0$) परावर्तन को दबाती है और संचरण को बढ़ाती है। जैसे-जैसे मोटाई बढ़ती है, कंजर्वेटिव मीडिया के लिए $R$ इकाई की ओर बढ़ता है, जिससे डिफ्यूसिव $T \sim 1/\tau$ टेल रिकवर होती है।
• अवशोषक मीडिया: $\tau \to \infty$ की सीमा में, ऑपरेटर अवशोषक मीडिया के लिए सटीक अर्ध-अनंत एल्बेडो को रिकवर करता है, जो $\le 6 \times 10^{-4}$ तक सटीक H-फंक्शन समाधान के अनुरूप है।

संरचनात्मक योगदान और महत्व
यह शोध पत्र इस बात पर जोर देता है कि इसका प्राथमिक योगदान संरचनात्मक है न कि केवल कम्प्यूटेशनल। केंद्रीय परिणाम परिमित-स्लैब परावर्तन का सटीक गुणनखंड (exact factorization) है:
$$R_\tau(a) = \sum_{n \ge 1} P_\infty(n) S(n, \tau) a^n$$
जहाँ:

• $P_\infty(n)$ एक हाफ-स्पेस फर्स्ट-रिटर्न लॉ (अर्ध-अनंत माध्यम में $n$-वें टकराव पर पलायन की संभावना) है।
• $S(n, \tau)$ एक डेप्थ-सर्वाइवल फैक्टर है, जो यह दर्शाता है कि $n$ क्रम के रिटर्न पाथ के $\tau$ से कम गहरा रहने की संभावना है।

यह गुणनखंड स्थापित करता है कि परिमित स्लैब समस्या और हाफ-स्पेस रिटर्न सांख्यिकी "एक ही ऑपरेटर के दो चेहरे" हैं। जैसे-जैसे $\tau \to \infty$ होता है, स्लैब परावर्तन नियम $P_R(n, \tau)$ टर्म-दर-टर्म हाफ-स्पेस लॉ $P_\infty(n)$ की ओर अभिसरित होता है, जिससे स्पारे-एंडरसन/कॉची-कर्नेल (Sparre–Andersen/Cauchy-kernel) टेल व्यवहार प्राप्त होता है।

भाग II: सर्वाइवल फैक्टर पर वर्किंग नोट्स
शोध पत्र में $S(n, \tau)$ के विश्लेषणात्मक रूप के संबंध में प्रारंभिक, अनुभवजन्य नोट्स (भाग II) शामिल हैं।

• स्केलिंग: अनुभवजन्य डेटा बताता है कि $S(n, \tau)$ ट्रांसपोर्ट मीन फ्री पाथ $\ell^* = 1/(1-g)$ और एक एडिटिव डेप्थ ऑफसेट $c_g$ के आधार पर एक स्केलिंग कोलैप्स (scaling collapse) का पालन करता है।
• ब्राउनियन एनालॉजी: सर्वाइवल फैक्टर एक ब्राउनियन एक्सकर्शन मैक्सिमम (Brownian excursion maximum) के संचयी वितरण फलन (CDF) का अनुसरण करता प्रतीत होता है, जिसे $g$-निर्भर ऑफसेट और एक हल्के आकार सुधार द्वारा संशोधित किया गया है।
• स्थिति: ये निष्कर्ष स्पष्ट रूप से "वर्किंग नोट्स" और "अनुभवजन्य" के रूप में लेबल किए गए हैं, जो प्रथम सिद्धांतों से व्युत्पत्ति की प्रतीक्षा कर रहे हैं। इन्हें भाग I के फॉरवर्ड मॉडल के सत्यापित परिणामों के रूप में दावा नहीं किया गया है।

सीमाएं और दायरा
लेखक नोट करते हैं कि इस विधि का उद्देश्य एकल रेडिएंस की गणना के लिए अनुकूलित डिस्क्रीट-ऑर्डिनेट कोड (जैसे DISORT) से बेहतर प्रदर्शन करना नहीं है। इसके बजाय, यह पूरक है, जो क्रम-अनुरूप कानूनों और हाफ-स्पेस सांख्यिकी के साथ स्पष्ट संबंध प्रदान करता है जिसे अन्य सॉल्वर सीधे उजागर नहीं करते हैं। वर्तमान कार्यान्वयन अज़ीमुथल औसत कोणीय वितरण प्रदान करता है; तिरछी आपतन (oblique incidence) के लिए पूर्ण BRDF/BTDF निर्भरता को एक सरल विस्तार के रूप में नोट किया गया है। शोध पत्र इस बात के साथ समाप्त होता है कि स्लैब ऑप्टिक्स और हाफ-स्पेस रिटर्न सांख्यिकी का प्रथम-प्रवेश ऑपरेटर के माध्यम से एकीकरण ही मुख्य वैचारिक प्रगति है।