Determination of active forces in actomyosin systems as inverse source problems for the Stokes equation

तकनीकी सारांश: स्टोक समीकरण के लिए व्युत्क्रम स्रोत समस्याओं के रूप में एक्टोमायोसिन प्रणालियों में सक्रिय बलों का निर्धारण

समस्या विवरण
इस कार्य का केंद्रीय उद्देश्य उन सक्रिय बलों की पहचान करना है जो फिलामेंट-मोटर नेटवर्क (विशेष रूप से F-एक्टिन और मायोसिन, जिन्हें एक्टोमायोसिन के रूप में जाना जाता है) द्वारा उत्पन्न होते हैं जो जैविक प्रणालियों में तरल प्रवाह को संचालित करते हैं। जबकि परिणामी तरल वेग क्षेत्र (velocity fields) को ऑप्टिकल माइक्रोस्कोपी और पार्टिकल इमेज वेलोसिटी (PIV) के माध्यम से प्रयोगात्मक रूप से मापा जा सकता है, अंतर्निहित बल सीधे सुलभ नहीं होते हैं। लेखक इसे निम्न-रेनॉल्ड्स संख्या हाइड्रोडायनामिक्स के संदर्भ में एक व्युत्क्रम स्रोत समस्या (inverse source problem) के रूप में तैयार करते हैं: प्रेक्षित वेग क्षेत्र $v$ से बल क्षेत्र $f$ का पुनर्निर्माण करना।

अध्ययन दो अलग-अलग भौतिक परिवेशों को संबोधित करता है:

1. संयमित प्रणालियाँ (Confined Systems): पानी-में-तेल की बूंदों में समाहित एक्टोमायोसिन नेटवर्क, जिन्हें रॉबिन (स्लिप) सीमा स्थितियों (boundary conditions) के साथ मॉडल किया गया है।
2. बल्क प्रणालियाँ (Bulk Systems): बिना किसी दृश्य कैप्सूलीकरण के स्वतंत्र रूप से निलंबित नेटवर्क, जिन्हें असजातीय डिरिचलेट (inhomogeneous Dirichlet) सीमा स्थितियों के साथ मॉडल किया गया है।

एक महत्वपूर्ण चुनौती पहचानी गई है कि प्रयोगात्मक डेटा अपूर्ण है; विशेष रूप से, दबाव क्षेत्र $p$ अपरिभाषित/अमापित है, और श्यानता (viscosity) $\eta$ को आसपास के तरल के सापेक्ष फिलामेंट नेटवर्क के नगण्य आयतन अंश के आधार पर स्थिर माना गया है।

कार्यप्रणाली
लेखक स्टोक समीकरण को अग्रवर्ती मॉडल (forward model) के रूप में उपयोग करते हैं, जिसमें असंपीड्यता की स्थिति ($\nabla \cdot v = 0$) और स्थानिक रूप से परिवर्तनशील श्यानता $\eta$ शामिल है। गणितीय ढांचा निम्नलिखित चरणों पर आधारित है:

1. अग्रवर्ती समस्या विश्लेषण (Forward Problem Analysis): लेखक दोनों प्रकार की सीमा स्थितियों (रॉबिन और डिरिचलेट) के लिए स्टोक समस्या की सु-निर्धारितता (well-posedness) का कठोरता से विश्लेषण करते हैं। वे उपयुक्त फलन स्थानों (जैसे, रॉबिन स्थितियों के लिए $H^1_{\parallel}(\Omega)^d$) को परिभाषित करते हैं और इन्फ-सप (inf-sup) स्थिति और कॉर्न की असमानता (Korn's inequality) का उपयोग करके वेग-दबाव युग्म $(v, p)$ के अद्वितीय अस्तित्व और स्थिरता को सिद्ध करते हैं।
2. व्युत्क्रम समस्या सूत्रीकरण (Inverse Problem Formulation): व्युत्क्रम समस्या को ऐसे $f$ को खोजने के रूप में परिभाषित किया गया है जिससे अग्रवर्ती ऑपरेटर $A$ (बल से वेग तक का मानचित्रण) मापित डेटा $v_{meas}$ से मेल खाता हो। दबाव डेटा के अभाव के कारण, इस समस्या को केवल वेग के आधार पर बल की पहचान के रूप में माना गया है।
3. सैद्धांतिक अपघटन (Theoretical Decomposition): शोध पत्र बल के अपसरण-मुक्त (divergence-free/solenoidal) और चक्र-मुक्त (curl-free/irrotational) घटकों में हेल्महोल्ट्ज़ अपघटन का विश्लेषण करता है। यह दिखाया गया है कि स्थिर श्यानता के लिए, बल का अपसरण-मुक्त घटक स्टोक समीकरण में वेग-निर्भर पद के अनुरूप होता है, जबकि चक्र-मुक्त घटक दबाव प्रवणता (pressure gradient) के अनुरूप होता है। फलस्वरूप, दबाव डेटा के बिना, सक्रिय बल के केवल अपसरण-मुक्त घटक को ही अद्वितीय रूप से पुनर्निर्मित किया जा सकता है।
4. एडजॉइंट व्युत्पत्ति (Adjoint Derivation): संख्यात्मक समाधान को सक्षम करने के लिए, लेखक अग्रवर्ती ऑपरेटरों और उनके संबंधित एडजॉइंट्स (बनाच और हिल्बर्ट स्पेस एडजॉइंट्स) के फ्रेचे डेरिवेटिव (Fréchet derivatives) को व्युत्पन्न करते हैं। यह व्युत्पत्ति ग्रेडिएंट-आधारित नियमितीकरण (regularization) के लिए आवश्यक है।
5. संख्यात्मक नियमितीकरण (Numerical Regularization): व्युत्क्रम समस्या को लैंडवेबर पुनरावृत्ति (Landweber iteration) विधि का उपयोग करके हल किया जाता है। एल्गोरिदम मापित वेग और सिम्युलेटेड वेग के बीच विसंगति को कम करता है। शोर (noise) के विरुद्ध स्थिरता सुनिश्चित करने के लिए, पुनरावृत्ति को विसंगति सिद्धांत (discrepancy principle) का उपयोग करके 'ए पोस्टेरिओरी' (a posteriori) समाप्त किया जाता है। प्रारंभिक अनुमान एक "नाइव" (naive) पुनर्निर्माण है जो सीधे स्टोक समीकरण के वेग पद ($f_0 = -\nabla \cdot \eta D(v_{meas})$) से प्राप्त किया गया है।

मुख्य परिणाम
कार्यप्रणाली को कृत्रिम डेटा (synthetic data) और प्रयोगात्मक माप दोनों का उपयोग करके मान्य किया गया:

• कृत्रिम डेटा: पुनर्निर्माण बूंदों (रॉबिन) और बल्क (डिरिचलेट) परिदृश्यों का प्रतिनिधित्व करने वाले 2D डोमेन पर किए गए।
  • शोर-मुक्त डेटा के लिए, लैंडवेबर पुनरावृत्ति ने जटिल भंवर पैटर्न (swirl patterns) सहित वास्तविक बल (ground truth force) को सफलतापूर्वक पुनर्निर्मित किया, जिसमें निम्न सापेक्ष त्रुटि ($\tilde{\epsilon} \approx 0.3$ रॉबिन के लिए, $\approx 0.73$ डिरिचलेट के लिए) देखी गई।
  • शोर युक्त डेटा के लिए, विधि ने मजबूती प्रदर्शित की। $\tilde{\delta} = 1.37$ तक के सापेक्ष शोर स्तरों के साथ भी, सापेक्ष पुनर्निर्माण त्रुटि सीमित रही ($\tilde{\epsilon} \approx 0.29$), जिससे बल की सामान्य संरचना और परिमाण सुरक्षित रहा, हालांकि उच्च शोर स्तरों पर सूक्ष्म विवरण खो गए।
• प्रयोगात्मक डेटा: एल्गोरिदम को बूंदों और बल्क नमूनों में एक्टोमायोसिन नेटवर्क के PIV डेटा पर लागू किया गया।
  • बूंदें (Droplets): पुनर्निर्मित बल क्षेत्र ने बूंद के भीतर देखे गए प्रवाह की दिशा के अनुरूप एक अपसरण-मुक्त घटक दिखाया।
  • बल्क (Bulk): पुनर्निर्माण बिना कैप्सूलीकरण वाले नेटवर्क पर किए गए। लेखक नोट करते हैं कि बल्क परिवेश में पुनर्निर्माण अधिक चुनौतीपूर्ण है क्योंकि लुप्त चक्र-मुक्त (दबाव-संबंधित) बल घटक की संभावित महत्ता और इस बात की अस्पष्टता बनी रहती है कि देखे गए प्रवाह केवल एक्टोमायोसिन गतिविधि से उत्पन्न हुए हैं या तरल-परिवेश अंतःक्रियाओं से।

महत्व और दावे
यह शोध पत्र उन जैविक प्रणालियों में सक्रिय बलों की पहचान के लिए एक कठोर गणितीय ढांचा प्रदान करने का दावा करता है जहाँ प्रत्यक्ष मापन असंभव है। इसके प्राथमिक योगदान हैं:

1. गणितीय कठोरता (Mathematical Rigor): सक्रिय पदार्थ (active matter) के लिए प्रासंगिक विशिष्ट सीमा स्थितियों (रॉबिन और डिरिचलेट) के साथ अग्रवर्ती स्टोक समस्याओं की सु-निर्धारितता स्थापित करना और नियमितीकरण के लिए आवश्यक एडजॉइंट ऑपरेटरों को व्युत्पन्न करना।
2. अपूर्ण डेटा का प्रबंधन: स्थिर श्यानता धारणाओं के तहत यह प्रदर्शित करके दबाव डेटा की कमी को स्पष्ट रूप से संबोधित करना कि अपसरण-मुक्त घटक ही पुनर्निर्मित करने योग्य मात्रा है।
3. प्रूफ ऑफ कॉन्सेप्ट: कृत्रिम और वास्तविक प्रयोगात्मक डेटा दोनों पर व्युत्क्रम दृष्टिकोण को मान्य करके, यह दिखाना कि वेग क्षेत्रों से ही सक्रिय बलों का अनुमान लगाया जा सकता है।

लेखक इस कार्य को मौजूदा बल पहचान तकनीकों (जैसे ट्रैक्शन फोर्स माइक्रोस्कोपी) के पूरक के रूप में देखते हैं, लेकिन इसे तरल वातावरण और सक्रिय पदार्थ के लिए अनुकूलित किया गया है। वे विनम्रतापूर्वक कहते हैं कि हालांकि विधि सफलतापूर्वक प्रमुख अपसरण-मुक्त बलों को पुनर्निर्मित करती है, लेकिन सक्रिय बल (चक्र-मुक्त भाग सहित) का पूर्ण पुनर्निर्माण के लिए दबाव क्षेत्र के ज्ञान या आगे के अनुमानों की आवश्यकता होगी, जो भविष्य के अध्ययन का विषय बना हुआ है। वर्तमान कार्य इन प्रणालियों के स्थिर-अवस्था (steady-state), 2D सन्निकटन पर केंद्रित है।