A Novel Compartmental Model for the Mathematical Epidemiology of Addiction: Backward Bifurcation, Bistability, Persistence, and Treatment-Capacity Policy Implications

तकनीकी सारांश: व्यसन के गणितीय महामारी विज्ञान के लिए एक नवीन कम्पार्टमेंटल मॉडल

समस्या विवरण
व्यसन (Addiction) के महामारी विज्ञान के मौजूदा गणितीय मॉडल अक्सर रैखिक स्थानांतरण दरों (linear transfer rates) पर निर्भर करते हैं और वास्तविक दुनिया की उपचार प्रणालियों की दो अनुभवजन्य रूप से महत्वपूर्ण विशेषताओं को पकड़ने में विफल रहते हैं: (1) सीमित उपचार क्षमता, जहाँ बढ़ती हुई व्यसनी आबादी के कारण उपचार में प्रवेश की दर संतृप्त हो जाती है (वेटिंग लिस्ट और निश्चित बजट के कारण); और (2) प्रत्यक्ष रिलैप्स (Direct relapse), जहाँ स्वस्थ हुए व्यक्ति पूरी तरह से सुग्राह्य (susceptible) अवस्था में लौटने के बजाय सीधे सक्रिय व्यसन की ओर लौट आते हैं। मानक मॉडलों में इन तंत्रों की अनुपस्थिति 'बैकवर्ड बिफर्केशन' (backward bifurcation) की संभावना को छिपा देती है, जो एक ऐसी घटना है जहाँ मूल प्रजनन संख्या ($R_0$) को इकाई से कम करने के बावजूद व्यसन को समाप्त करना अपर्याप्त होता है, क्योंकि सिस्टम एक स्थिर, उच्च-प्रचलन वाले स्थानिक (endemic) अवस्था में फंसा रह सकता है।

कार्यप्रणाली
लेखक एक चार-कम्पार्टमेंट वाला नियतात्मक कोहोर्ट मॉडल ($S, A, T, R$) विकसित करते हैं जो सुग्राह्य (Susceptible), सक्रिय रूप से व्यसनी (Actively Addicted), उपचार में (In Treatment), और स्वस्थ (Recovered) आबादी का प्रतिनिधित्व करता है। मॉडल तीन नवीन गैर-रैखिक/संरचनात्मक घटकों को शामिल करता है:

1. संतृप्त उद्दीपन (Saturating incidence): $\frac{\beta SA}{1+\alpha A}$, जो सीमित एक्सपोजर क्षमता को दर्शाता है।
2. क्षमता-सीमित उपचार प्रवेश: $\frac{\sigma A}{1+\gamma A}$, जहाँ $\gamma$ क्षमता बाधाओं की गंभीरता को मापता है।
3. प्रत्यक्ष रिलैप्स: दर $\psi$ पर $R$ से सीधे $A$ की ओर एक मार्ग।

अध्ययन एक कठोर गणितीय टूलकिट का उपयोग करता है:

• साम्यावस्था विश्लेषण (Equilibrium Analysis): एडिक्शन-फ्री इक्विलिब्रियम (AFE) और स्थानिक (endemic) साम्यावस्थाओं का व्युत्पन्न, जो मूल (root) के अस्तित्व का विश्लेषण करने के लिए सिस्टम को एक द्विघात समीकरण में कम करता है।
• स्थिरता सिद्धांत (Stability Theory): नेक्स्ट-जेनरेशन मैट्रिक्स विधि ($R_0$) के माध्यम से स्थानीय स्थिरता विश्लेषण और मेट्ज़लर-मैट्रिक्स तुलना तर्क तथा लियापुनोव कार्यों (Lyapunov functions) का उपयोग करके वैश्विक स्थिरता प्रमाण।
• बिफर्केशन सिद्धांत (Bifurcation Theory): बैकवर्ड बिफर्केशन के लिए एक सटीक, क्लोज्ड-फॉर्म मानदंड प्राप्त करने के लिए कास्टिलो-चावेज़-सोंग सेंटर-मैनिफ़ोल्ड प्रमेय का अनुप्रयोग।
• संवेदनशीलता विश्लेषण (Sensitivity Analysis): $R_0$ और बिफर्केशन गुणांक पर पैरामीटर के प्रभाव को निर्धारित करने के लिए स्थानीय (सामान्यीकृत संवेदनशीलता सूचकांक) और वैश्विक (सोबोल सूचकांक और PRCC) विश्लेषण।
• विस्तार (Extensions): जनसांख्यिकीय शोर (demographic noise) के प्रभावों का आकलन करने के लिए स्टोकेस्टिक डिफ्यूजन एप्रोक्सिमेशन (Itô SDE) का उपयोग, और उपचार क्षमता एवं रोकथाम निवेश के लिए अनुकूलतम नियंत्रण (optimal control) सूत्रीकरण (पोंट्रियागिन के अधिकतम सिद्धांत द्वारा)।
• अनुभवजन्य अंशांकन (Empirical Calibration): सैद्धांतिक निष्कर्षों की प्रासंगिकता का परीक्षण करने के लिए अमेरिकी NSDUH हेरोइन-उपयोग डेटा (2013-2022) के साथ मॉडल को फिट करना।

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

1. क्लोज्ड-फॉर्म बिफर्केशन मानदंड: शोध पत्र बिफर्केशन गुणांक $a$ के लिए एक सटीक अभिव्यक्ति व्युत्पन्न करता है। विश्लेषण प्रकट करता है कि उपचार-क्षमता संतृप्ति ($\gamma$) ही द्विरूपता (bistability) का एकमात्र स्रोत है (बैकवर्ड बिफर्केशन), जबकि उद्दीपन संतृप्ति ($\alpha$) पूरी तरह से स्थिर करने वाली है। गुणांक इस प्रकार है:
   $$a = 2\beta^*v_2w_1w_2 - \frac{2\Lambda\alpha\beta^*}{\mu}v_2w_2^2 + 2\gamma\sigma w_2^2(v_2 - v_3)$$
   यह सिद्ध करता है कि बैकवर्ड बिफर्केशन तभी होता है जब $\gamma > 0$ हो और रिलैप्स का दबाव पर्याप्त ($v_2 > v_3$) हो।

2. क्षमता की $R_0$ के प्रति "अदृश्यता": एक केंद्रीय सैद्धांतिक निष्कर्ष यह है कि मूल प्रजनन संख्या $R_0$ उपचार क्षमता पैरामीटर $\gamma$ से स्वतंत्र है ($\partial R_0 / \partial \gamma \equiv 0$)। फलस्वरूप, मानक निगरानी प्रणाली $R_0$ का उपयोग करके क्षमता-जनित द्विरूपता के जोखिम का पता नहीं लगा सकती है। $R_0$ संचरण और जनसांख्यिकीय दरों द्वारा नियंत्रित होता है, जबकि बिफर्केशन का प्रकार (फॉरवर्ड बनाम बैकवर्ड) $\gamma$ और $\psi$ द्वारा नियंत्रित होता है।

3. वैश्विक स्थिरता और निरंतरता (Global Stability and Persistence):

   • AFE की वैश्विक स्थिरता के लिए एक पर्याप्त स्थिति एक मेट्ज़लर-मैट्रिक्स तुलना तर्क के माध्यम से स्थापित की गई है, जो दिखाती है कि रैखिक सीमाएँ ठीक उसी क्षेत्र में विफल हो जाती हैं जहाँ द्विरूपता होती है।
   • $R_0 > 1$ के लिए एक यूनिफॉर्म पर्सिस्टेंस प्रमेय सिद्ध किया गया है, जो सुनिश्चित करता है कि आंतरिक प्रारंभिक स्थितियों से व्यसन समाप्त नहीं हो सकता।
   • फॉरवर्ड बिफर्केशन शासन ($a < 0$) के लिए, अद्वितीय स्थानिक साम्यावस्था के लिए एक वैश्विक स्थिरता परिणाम प्रदान किया गया है, जो उद्दीपन घटक के लिए एक हैंड-प्रूवन लियापुनोव फंक्शन और उपचार/रिलैप्स चक्र के लिए कम्प्यूटेशनल सर्टिफिकेशन द्वारा समर्थित है।

4. सैडल-नोड "घोस्ट" तंत्र (Saddle-Node "Ghost" Mechanism): अमेरिकी हेरोइन डेटा के अंशांकन से एक ऐसी स्थिति का पता चलता है जहाँ $R_0 < 1$ है लेकिन $R_0 > R_0^c$ (सैडल-नोड थ्रेशोल्ड) है। इस शासन में, स्थानिक साम्यावस्थाएँ लुप्त हो चुकी हैं, फिर भी सिस्टम सैडल-नोड बिफर्केशन के एक "घोस्ट" (ghost) को प्रदर्शित करता है। इसके परिणामस्वरूप, सदियों लंबे क्वासी-स्थानिक ट्रांजिएंट्स (quasi-endemic transients) उत्पन्न होते हैं (जैसे, 150 वर्षों में व्यापकता में 20% की कमी), जिससे जनसंख्या किसी भी यथार्थवादी नीति क्षितिज पर स्थिर और स्थानिक प्रतीत होती है, भले ही सैद्धांतिक अट्रैक्टर 'रोग-मुक्त अवस्था' हो।

5. स्टोकेस्टिक मजबूती (Stochastic Robustness): दिखाया गया है कि जनसांख्यिकीय शोर के तहत रोग-मुक्त और स्थानिक अवस्थाओं के बीच का नियततात्मक बेसिन बाउंड्री एक सुचारू प्रायिकता प्रवणता (probability gradient) बन जाता है। सीमित आबादी में, स्थानिक शाखा तक पहुँचने की संभावना प्रारंभिक स्थितियों का एक निरंतर फलन (continuous function) है, न कि एक कठोर थ्रेशोल्ड।

6. अनुकूलतम नियंत्रण (Optimal Control): एक अनुकूलतम नियंत्रण समस्या तैयार की गई है जो सक्रिय व्यसन और हस्तक्षेप लागतों को कम करती है। परिणाम संकेत देते हैं कि जबकि अनुकूलतम नियंत्रण प्रकोप के बोझ को कम कर सकता है, यह स्वयं बिफर्केशन शासन को नहीं बदलता है, जब तक कि नियंत्रण स्पष्ट रूप से संतृप्ति पैरामीटर $\gamma$ को लक्षित न करे (जैसे, केवल स्लॉट्स जोड़ने के बजाय कतार संबंधी घर्षणों को हटाना)।

महत्व और नीतिगत निहितार्थ
शोध पत्र तर्क देता है कि व्यसन नियंत्रण के लिए $R_0 < 1$ को कम करने का मानक नीतिगत मीट्रिक अपर्याप्त है जब उपचार क्षमता संतृप्त होती है। लेखक दावा करते हैं:

• $R_0 < 1$ पर्याप्त नहीं है: द्विरूपता विंडो ($R_0^c < R_0 < 1$) में, सिस्टम एक उच्च-व्यसन अवस्था में फंसा रह सकता है। उन्मूलन के लिए या तो व्यापकता में भारी एकमुश्त कमी (सिस्टम को बेसिन बाउंड्री के नीचे धकेलना) या उन संरचनात्मक परिवर्तनों की आवश्यकता है जो $R_0$ को महत्वपूर्ण थ्रेशोल्ड $R_0^c$ से नीचे ले जाएं।
• क्षमता नीति $R_0$ के लिए अदृश्य है: उपचार क्षमता का विस्तार करना ($\gamma$ को कम करना) $R_0$ को कम नहीं करता है; इसका प्रभाव पूरी तरह से बिफर्केशन आरेख के आकार पर पड़ता है। इसलिए, केवल $R_0$ को ट्रैक करने वाली निगरानी प्रणालियाँ क्षमता-जनित पतन के करीब आने का पता लगाने में विफल रहेंगी।
• जागरूकता के दोहरे लाभ: जागरूकता बढ़ाना ($\alpha$ बढ़ाना) दोहरा लाभ प्रदान करता है: यह $R_0$ को कम करता है और साथ ही सिस्टम को बैकवर्ड बिफर्केशन के विरुद्ध स्थिर करता है।
• संरचनात्मक बनाम मात्रात्मक हस्तक्षेप: केवल उपचार चाहने की दर ($\sigma$) को बढ़ाना बिना वास्तविक क्षमता ($\gamma$) के विस्तार के, सिस्टम को बैकवर्ड ब्रांच के साथ आगे बढ़ाता है लेकिन द्विरूपता ट्रैप को समाप्त नहीं करता है। केवल $\gamma$ को कम करना ही फॉरवर्ड-बिफर्केशन शासन में प्रवेश कर सकता है।

अध्ययन निष्कर्ष निकालता है कि प्रभावी व्यसन नीति को केवल संचरण में कमी या उपचार चाहने को बढ़ावा देने पर निर्भर रहने के बजाय, सुरक्षा मार्जिन $R_0/R_0^c$ की निगरानी करनी चाहिए और संरचनात्मक क्षमता बाधाओं को सीधे संबोधित करना चाहिए।