Solving the Fokker-Planck equation of discretized Dean-Kawasaki models with functional hierarchical tensor

तकनीकी सारांश: कार्यात्मक पदानुक्रमित टेंसर (Functional Hierarchical Tensor) के साथ विविक्त डीन-कावाकी मॉडल के फॉकर-प्लांक समीकरण को हल करना

समस्या विवरण
यह शोध पत्र विविक्त डीन-कावाकी (DK) मॉडल से जुड़े समय-निर्भर फॉकर-प्लांक समीकरण (FPE) को हल करने की चुनौती को संबोधित करता है, जो परस्पर क्रिया करने वाले कणों के घनत्व उतार-चढ़ाव का वर्णन करने वाला एक स्टोकेस्टिक पार्शियल डिफरेंशियल इक्वेशन (SPDE) है। यद्यपि DK मॉडल गैर-संतुलन सांख्यिकीय भौतिकी, फ्लक्चुएटिंग हाइड्रोडायनामिक्स और एजेंट-आधारित मॉडलिंग के लिए मौलिक है, लेकिन इसके संबद्ध FPE को संख्यात्मक रूप से हल करना अत्यधिक गणनात्मक रूप से कठिन है। प्राथमिक कठिनाई अवस्था स्थान (state space) की उच्च विमा (dimensionality) से उत्पन्न होती है: $d$ ग्रिड सेल में स्थानिक विवक्तीकरण (spatial discretization) के बाद, अवस्था $\pi$ कण घनत्वों का एक वेक्टर है जो एक $(d-1)$-आयामी प्रायिकता सिम्प्लेक्स (probability simplex) द्वारा बाधित है (जहाँ प्रविष्टियाँ गैर-ऋणात्मक हैं और योग एकता के बराबर है)। FPE के लिए पारंपरिक ग्रिड-आधारित विधियाँ 'कर्स ऑफ डाइमेंशनैलिटी' से ग्रस्त होती हैं, और मानक घनत्व अनुमान तकनीकें सिम्प्लेक्स बाधा और DK मॉडल में अंतर्निहित मजबूत स्थानिक सहसंबंधों (spatial correlations) के साथ संघर्ष करती हैं।

कार्यप्रणाली
लेखक एक कण-आधारित वर्कफ़्लो प्रस्तावित करते हैं जो वेवलेट आधार (wavelet basis) पर एक कार्यात्मक पदानुक्रमित टेंसर नेटवर्क (FHT-W) का उपयोग करते हुए, विविक्त DK मॉडल के समय-निर्भर नियम का घनत्व अनुमान के माध्यम से सन्निकटन करते हैं। यह कार्यप्रणाली चार विशिष्ट चरणों में आगे बढ़ती है:

1. स्थानिक और टेम्पोरल विवक्तीकरण (Spatial and Temporal Discretization):
   एक आवधिक डोमेन (periodic domain) पर एक परिमित-आयतन योजना (finite-volume scheme) का उपयोग करके निरंतर DK SPDE को विविक्त किया जाता है। यह SPDE को सेल-औसत घनत्व $\Pi(t)$ को नियंत्रित करने वाले स्टोकेस्टिक डिफरेंशियल इक्वेशंस (SDEs) की एक प्रणाली में बदल देता है। विवक्तीकरण बाहरी और युग्मवार विभव (external and pairwise potentials) से जुड़े नियतात्मक प्रवाह (deterministic fluxes) के लिए एक रो स्कीम (Roe scheme) और एक विशिष्ट स्टोकेस्टिक फ्लक्स फॉर्मूलेशन का उपयोग करता है, जो संख्यात्मक स्थिरता सुनिश्चित करता है और धनात्मकता (positivity) को बनाए रखता है। परिणामी SDEs को de facto मॉन्टे कार्लो प्रक्षेप पथों (Monte Carlo trajectories) के माध्यम से सिम्युलेट किया जाता है ताकि नमूना पथों का एक संग्रह बनाया जा सके।

2. निर्देशांक रूपांतरण (सिम्प्लेक्स से यूक्लिडियन स्पेस):
   मानक टेंसर नेटवर्क घनत्व अनुमान लागू करने के लिए, लेखक अवस्था को प्रायिकता सिम्प्लेक्स $\Delta$ से एक अनियंत्रित यूक्लिडियन स्पेस में मैप करते हैं। यह दो-चरणीय रूपांतरण के माध्यम से प्राप्त किया जाता है:

• सेंटर्ड लॉग-रेशियो (CLR) रूपांतरण: घनत्व वेक्टर $\pi$ को $s$ में मैप किया जाता है जहाँ $s_i = \log(\pi_i) - \frac{1}{d}\sum_j \log(\pi_j)$। यह सिम्प्लेक्स के आंतरिक भाग को $\mathbb{R}^d$ में एक हाइपरप्लेन में मैप करता है जहाँ $\sum s_i = 0$, जो प्रभावी रूप से सिम्प्लेक्स बाधा को हटा देता है।
• वेवलेट रूपांतरण: $s$-निर्देशांकों को फिर एक डिस्क्रीट वेवलेट ट्रांसफॉर्म (DWT) का उपयोग करके वेवलेट गुणांक $c$ में परिवर्तित किया जाता है। यह चरण जालक मॉडल (lattice model) की बहु-पैमाने (multi-scale) संरचना का लाभ उठाता। चूंकि प्रसार (diffusion) के कारण DK मॉडल में आसन्न ग्रिड सेल अत्यधिक सहसंबंधित होते हैं, इसलिए वेवलेट रूपांतरण विभिन्न पैमानों पर चरों को डिकोरिलेट (decorrelate) करता है, जिससे परिणामी प्रायिकता वितरण $P_c(c, t)$ की संख्यात्मक रैंक काफी कम हो जाती है।

3. कार्यात्मक पदानुक्रमित टेंसर (FHT-W) एंसेट (Ansatz):
   वेवलेट समन्वय स्थान में प्रायिकता घनत्व का सन्निकटन एक कार्यात्मक पदानुक्रमित टेंसर (FHT) नेटवर्क का उपयोग करके किया जाता है। विशेष रूप से, लेखक FHT-W संस्करण का उपयोग करते हैं, जो एक Hierarchical Tucker (HT) ट्री टेंसर नेटवर्क है जहाँ पत्तियाँ (leaves) वेवलेट गुणांकों के अनुरूप होती हैं।

• ट्री संरचना वेवलेट रूपांतरण के डायडिक अपघटन (dyadic decomposition) को दर्शाती है।
• घनत्व को एक निश्चित एकचर आधार (लेजेंड्रे बहुपदों) के विस्तार के रूप में दर्शाया जाता है जिसके गुणांक एक निम्न-रैंक टेंसर प्रारूप में व्यवस्थित होते हैं।
• यह प्रतिनिधित्व यह सुनिश्चित करता है कि भंडारण जटिलता (storage complexity) आयाम $d$ के साथ रैखिक रूप से स्केल करे, न कि घातांकीय रूप से, बशर्ते कि पदानुक्रमित रैंक कम बनी रहें।

4. इंटरपोलेशन के माध्यम से अवलोकन अनुमान (Observable Estimation via Interpolation):
   सीखे गए घनत्व से भौतिक अवलोकनों (जैसे एन्ट्रॉपी) की गणना करने के लिए, लेखक एक इंटरपोलेशन योजना प्रस्तावित करते हैं। चूंकि मूल $\pi$-स्पेस में परिभाषित अवलोकन वेवलेट $c$-स्पेस में जटिल फलन बन जाते हैं, उन्हें क्रॉस-इंटरपोलेशन विधि का उपयोग करके FHT-W प्रारूप में संकुचित किया जाता है। इसमें अवलोकन फलन के बिंदु मूल्यांकन (point evaluations) से प्राप्त संरचित रैखिक समीकरणों को हल करना शामिल है, जिससे घनत्व टेंसर के साथ कुशल संकुचन (contraction) संभव होता है।

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

• एक व्यापक कम्प्यूटेशनल फ्रेमवर्क: यह शोध पत्र DK मॉडल के अध्ययन के लिए एक पूर्ण पाइपलाइन प्रदान करता है, जो घनत्व अनुमान के माध्यम से स्टोकेस्टिक कण सिमुलेशन और उच्च-आयामी PDE समाधान के बीच के अंतर को पाटता है।
• सिम्प्लेक्स बाधाओं का प्रबंधन: सेंटर्ड लॉग-रेशियो रूपांतरण के साथ वेवलेट निर्देशांकों का परिचय प्रायिकता सिम्प्लेक्स समस्याओं पर टेंसर नेटवर्क विधियों को लागू करने में सक्षम बनाता है, जो एक ऐसा क्षेत्र है जहाँ मानक यूक्लिडियन टेंसर नेटवर्क विफल हो जाते हैं।
• उच्च आयामों तक स्केलेबिलिटी: यह पद्धति $d=64$ स्वतंत्रता की डिग्री वाले सिस्टम (1D और 2D दोनों मामलों में) को सफलतापूर्वक संभालती है, जो एक ऐसी विमा है जिसे प्रत्यक्ष FPE सॉल्वर या साधारण घनत्व अनुमान के लिए दुर्गम माना जाता है।
• पोटेंशियल और इंटरैक्शन का एकीकरण: यह ढांचा दोनों बाहरी विभव ($V_1$) और युग्मवार इंटरैक्शन विभव ($V_2$) को समाहित करता है, जो जटिल साइट-वाइज सहसंबंधों और इंटरैक्शन प्रभावों को कैप्चर करता है।

परिणाम
1D और 2D में संख्यात्मक प्रयोगों से पता चलता है कि FHT-W दृष्टिकोण प्रायिकता नियम के समय विकास को सटीक रूप से कैप्चर करता है।

• सटीकता: विधि वास्तविक मॉडल के साइट-वाइज सहसंबंधों और अन्य सांख्यिकीय अवलोकनों को उच्च निष्ठा (fidelity) के साथ पुनरुत्पादित करती है।
• स्थिरता: यह वर्कफ़्लो "स्ट्रॉन्ग-डिफ्यूजन" शासन (regime) के धारणा के तहत स्थिर रहता है, जिसमें दुर्लभ घटनाओं (rare events) में विस्फोट को रोकने के लिए शोर (noise) पर क्लैम्पिंग तकनीक लागू की गई है।
• दक्षता: वितरण की प्रभावी रैंक को कम करने के लिए वेवलेट रूपांतरण का उपयोग करके, विधि प्रबंधनीय टेंसर रैंक (उदाहरण के लिए, बॉन्ड डायमेंशन $r=20$) के साथ $d=64$ के लिए सटीक परिणाम प्राप्त करती है, जिससे पूर्ण टेंसर प्रतिनिधित्व से जुड़ी घातांकीय भंडारण लागतों से बचा जा सकता है।

महत्व और दावे
लेखक इस कार्य को स्वयं एक नवीन टेंसर नेटवर्क आर्किटेक्चर के रूप में नहीं, बल्कि मौजूदा ट्री-आधारित टेंसर नेटवर्क (विशेष रूप से HT प्रारूप) के एक चुनौतीपूर्ण वर्ग के मजबूत अनुप्रयोग के रूप में प्रस्तुत करते हैं। प्राथमिक योगदान DK मॉडल के लिए एक व्यावहारिक, स्केलेबल फ्रेमवर्क का निर्माण है। शोध पत्र दावा करता है कि यह दृष्टिकोण फ्लक्चुएटिंग हाइड्रोडायनामिक्स के अध्ययन के लिए पारंपरिक कण विधियों के लिए एक व्यवहार्य विकल्प प्रदान करता है, विशेष रूप से उन शासन (regimes) में जहाँ केवल नमूना पथों के बजाय पूर्ण प्रायिकता वितरण को समझना आवश्यक होता है। यह कार्य रेखांकित करता है कि ज्यामितीय रूपांतरण (लॉग-रेशियो और वेवलेट) के साथ टेंसर नेटवर्क संपीड़न का संयोजन, परस्पर क्रिया करने वाले कण प्रणालियों में निहित आयामी बाधाओं को पार करने के लिए आवश्यक है।