Higher-order interactions for controlling time-delayed Kuramoto model

तकनीकी सारांश: समय-विलंबित कुरोमोटो मॉडल (Kuramoto model) के नियंत्रण के लिए उच्च-क्रम अंतःक्रियाएं

समस्या विवरण
युग्मित ऑसिलेटर प्रणालियों का सामूहिक सिंक्रोनाइज़ेशन (synchronization) न्यूरोनल डायनेमिक्स, पावर ग्रिड और स्वार्म रोबोटिक्स जैसे विभिन्न इंजीनियरिंग और जैविक अनुप्रयोगों के लिए मौलिक है। हालांकि, वास्तविक दुनिया के नेटवर्क में सिग्नल प्रसार की सीमित गति के कारण समय-विलंब (time delays) शामिल होते हैं। ये विलंब प्रणाली को अनंत-आयामी (delay-differential equations द्वारा वर्णित) बना देते हैं, जिससे विश्लेषण और नियंत्रण काफी चुनौतीपूर्ण हो जाता है। जबकि समय-विलंब द्वि-स्थिरता (bistability) और हिस्टेरेसिस (hysteresis) जैसे जटिल व्यवहार उत्पन्न कर सकते हैं, मानक गणितीय मॉडलिंग अक्सर इन्हें नजरअंदाज करती है या प्रथम-क्रम चरण-लैग सन्निकटन (जैसे कि साकागुची-कुरोमोटो मॉडल) पर निर्भर करती है। ये पारंपरिक युग्म-आधारित (pairwise) सन्निकटन विलंबित प्रणालियों की विशिष्ट गतिक विशेषताओं, विशेष रूप से सिंक्रोनाइज्ड और इनकोहेरेंट अवस्थाओं के बीच द्वि-स्थिरता को पुनरुत्पादित करने में विफल रहते हैं।

कार्यप्रणाली
लेखक समय-विलंबित कुरोमोटो मॉडल के लिए एक विलंब-मुक्त नियंत्रण ढांचे का प्रस्ताव करते हैं, जो समय-विलंबित युग्म अंतःक्रियाओं (pairwise interactions) को प्रभावी उच्च-क्रम (गैर-युग्म) अंतःक्रियाओं के रूप में व्याख्यायित करता है। कार्यप्रणाली निम्नलिखित चरणों के माध्यम से आगे बढ़ती है:

1. उच्च-क्रम सन्निकटन (Higher-Order Approximation): समय-विलंबित युग्म अंतःक्रियाओं वाले एक वैश्विक रूप से युग्मित $N$ फेज ऑसिलेटर सिस्टम और एक नियंत्रण इनपुट से शुरू करते हुए, लेखक छोटे कपलिंग स्ट्रेंथ $\epsilon$ के लिए विलंबित पदों का सन्निकटन करते हैं। इससे एक विलंब-मुक्त कुरोमोटो मॉडल प्राप्त होता जिसमें युग्म (pairwise) और प्रभावी तीन-शरीर (three-body/higher-order) दोनों प्रकार की अंतःक्रियाएं शामिल होती हैं, जहाँ समय विलंब $\tau$ चरण अंतराल (phase lags) के रूप में दिखाई देता है।
2. आयामी न्यूनीकरण (Dimensional Reduction): लेखक ऑसिलेटर के फेज के प्रायिकता घनत्व फलन (probability density function) पर ओट्ट-एंटोनसें (Ott–Antonsen - OA) एंसेट लागू करते हैं। प्राकृतिक आवृत्तियों के लॉरेंट्ज़ वितरण (जो समान ऑसिलेटर्स के मामले में सीमित हो जाता है) को मानते हुए, वे जटिल ऑर्डर पैरामीटर $z(t) = R(t)e^{i\Psi(t)}$ के लिए एक द्वि-आयामी गतिक प्रणाली व्युत्पन्न करते हैं।
3. औसत निकालना (Averaging): सिंक्रोनाइज़ेशन डिग्री $R$ की धीमी गतिकी को अलग करने के लिए, लेखक द्वितीय-क्रम औसत पद्धति (second-order averaging method) का उपयोग करते हैं। यह तीव्र चरण गतिकी को समाप्त करता है, जिसके परिणामस्वरूप ऑर्डर पैरामीटर $R$ के लिए एक एक-आयामी न्यूनीकृत साधारण अवकल समीकरण (ODE) प्राप्त होता है। यह न्यूनीकृत समीकरण स्पष्ट रूप से $O(\tau\epsilon^2)$ के पदों को बनाए रखता है, जो उच्च-क्रम अंतःक्रियाओं के प्रभावों को पकड़ते हैं।
4. नियंत्रण डिजाइन (Control Design): न्यूनीकृत समीकरण के आधार पर दो मीन-फील्ड फीडबैक नियंत्रण स्कीमें डिजाइन की गई हैं:
   • रैखिक फीडबैक (Linear Feedback): $u(t) = C R(t) \cos \Psi(t)$।
   • गैर-रैखिक फीडबैक (Nonlinear Feedback): $u(t) = C R(t) (R(t) - E) \cos \Psi(t)$, जो विशिष्ट सिंक्रोनाइज़ेशन स्तरों को लक्षित करने के लिए एक अतिरिक्त पैरामीटर $E$ पेश करता है।
5. सत्यापन (Validation): उच्च-क्रम न्यूनीकृत समीकरण से प्राप्त स्थिरता आरेखों (stability diagrams) की तुलना एक पारंपरिक युग्म सन्निकटन (जो $O(\tau\epsilon^2)$ पदों को छोड़ देता है) से की जाती है। इन भविष्यवाणियों को मूल समय-विलंबित डिले-डिफरेंशियल समीकरणों के संख्यात्मक सिमुलेशन के माध्यम से सत्यापित किया जाता है।

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

• श्रेष्ठ भविष्य कहनेवाला सटीकता (Superior Predictive Accuracy): संख्यात्मक सिमुलेशन प्रदर्शित करते हैं कि उच्च-क्रम सन्निकटन मूल समय-विलंबित प्रणाली की गतिकी को पारंपरिक युग्म सन्निकटन की तुलना में अधिक सटीक रूप से भविष्यवाणी करता है। इसे ग्लोबल रूट मीन स्क्वेर्ड एरर (RMSE) का उपयोग करके मापा गया है, जो विभिन्न विलंब मानों ($\tau$) के लिए उच्च-क्रम मॉडल के लिए लगातार कम रहता है।
• द्वि-स्थिरता का क्रियान्वयन (Realization of Bistability): यह अध्ययन रैखिक फीडबैक नियंत्रण का उपयोग करके इनकोहेरेंट अवस्था ($R=0$) और पूर्णतः सिंक्रोनाइज्ड अवस्था ($R=1$) के बीच द्वि-स्थिरता को सफलतापूर्वक पुनरुत्पादित करता है। महत्वपूर्ण रूप से, यह द्वि-स्थिर क्षेत्र केवल उच्च-क्रम न्यूनीकृत समीकरण द्वारा ही पकड़ा जाता है; युग्म सन्निकटन इस घटना की भविष्यवाणी करने में विफल रहता है, और केवल मोनोस्टेबल क्षेत्र दिखाता है।
• मध्यवर्ती अवस्थाओं का स्थिरीकरण (Stabilization of Intermediate States): गैर-रैखिक फीडबैक नियंत्रण योजना किसी भी वांछित मध्यवर्ती सिंक्रोनाइज़ेशन स्तर ($0 < R < 1$) को स्थिर करने में सक्षम बनाती है। $C$ और $E$ मापदंडों को ट्यून करके, सिस्टम को एक वांछित संदर्भ मान $R_{ref}$ पर स्थिर किया जा सकता है।
• जटिल द्वि-स्थिर शासन (Complex Bistable Regimes): यह ढांचा एक ऐसे द्वि-स्थिर शासन को प्रकट करता है जहाँ पूर्णतः सिंक्रोनाइज्ड अवस्था ($R=1$) एक मध्यवर्ती सिंक्रोनाइज्ड अवस्था ($R=R_{ref}$) के साथ सह-अस्तित्व में रहती है। यह जटिल गतिक संरचना उच्च-क्रम मॉडल द्वारा सटीक रूप से पकड़ी जाती है, लेकिन युग्म सन्निकटन में अनुपस्थित है।

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