तकनीकी सारांश: रीमान ज़ेटा फलन के सरल और सटीक सन्निकटन (Approximations)

समस्या विवरण

रीमान ज़ेटा फलन, $\zeta(s)$, के लिए काल्पनिक भाग $t = \text{Im}(s)$ के बड़े मानों की गणना करना महत्वपूर्ण चुनौतियाँ प्रस्तुत करता है। जबकि यूलर-मैकलॉरिन (Euler-Maclaurin) योग विधि कठोर त्रुटि सीमाएँ (error bounds) प्रदान करती है, इसकी कम्प्यूटेशनल जटिलता $O(|s|)$ के रूप में बढ़ती है, जो बड़े $t$ के लिए इसे अव्यावहारिक बनाती है। बड़े $t$ के लिए मानक दृष्टिकोण रीमान-सीगल (Riemann-Siegel) सूत्र है, जो एक मुख्य योग (main sum) और शेष पद (remainder term) के लिए एक स्पर्शोन्मुख विस्तार (asymptotic expansion) का उपयोग करके जटिलता को $O(t^{1/2})$ तक कम कर देता है। हालाँकि, रीमान-सीगल सूत्र में उच्च-क्रम के सुधार (higher-order corrections) में विशेष फलनों (विशेष रूप से फलन $F(z)$) के जटिल अवकलज (derivatives) शामिल होते हैं, जिन्हें सिंगुलैरिटीज़ (singularities) और जटिल सूत्रों के कारण संख्यात्मक रूप से मूल्यांकन करना कठिन होता है। इसके अलावा, गैलवे (Galway) जैसे मौजूदा उच्च-परिशुद्धता एल्गोरिदम के लिए इष्टतम सटीकता प्राप्त करने हेतु कई मापदंडों (parameters) को ट्यून करने की आवश्यकता होती है।

यह शोध पत्र ऐसे सन्निकटन (approximations) की आवश्यकता को संबोधित करता है जो कम्प्यूटेशनल रूप से कुशल हों (कम पदों की आवश्यकता हो) और उच्च-क्रम रीमान-सीगल सुधारों से जुड़ी कार्यान्वयन संबंधी कठिनाइयों के बिना, क्रिटिकल स्ट्रिप सहित किसी भी ऊर्ध्वाधर स्ट्रिप (vertical strips) में अत्यधिक सटीक हों।

कार्यप्रणाली (Methodology)

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

1. इंटीग्रल प्रतिनिधित्व (Integral Representation): सटीक रीमान-सीगल सूत्र के शेष पद को एक समाकल (integral) $I_M(s)$ के रूप में व्यक्त किया जाता है, जिसमें एक जटिल-मान वाला माप (measure) $\eta(dx) = \frac{e^{-\pi x^2}}{\cosh(\pi \theta x)} dx$ शामिल है, जहाँ $\theta = e^{-i\pi/4}$ है।
2. गौसियन क्वाड्रैचर (Gaussian Quadrature): स्पर्शोन्मुख श्रेणी विस्तार (asymptotic series expansions) के बजाय, लेखक इस समाकल को एक गौसियन क्वाड्रैचर नियम का उपयोग करके सन्निकट करते हैं। लक्ष्य नोड्स $\{x_{p,j}\}$ और भार (weights) $\{\omega_{p,j}\}$ का एक सेट खोजना है ताकि क्वाड्रैचर $4p+2$ परीक्षण फलनों (test functions) के एक विशिष्ट सेट के लिए सटीक हो।
3. गुणांक व्युत्पत्ति (Coefficient Derivation): गुणांकों $\{\omega_{p,j}\}$ और $\{\lambda_{p,j}\}$ (जहाँ $\lambda_{p,j} = \theta x_{p,j}$) को रैखिक समीकरणों के एक निकाय को हल करके निकाला जाता है। यह निकाय इस शर्त के साथ बनाया गया है कि क्वाड्रैचर $f_k(x) = e^{-\pi x^2 + 2\pi \theta y_{p,k} x}$ प्रकार के फलनों के लिए सटीक हो, जहाँ $y_{p,k}$ अंतराल $[-1, 1]$ में समान रूप से वितरित बिंदु हैं।
   • समस्या को एक रैखिक फलन $L$ के लिए नोड्स और भार खोजने में परिवर्तित किया जाता है, जो मोर्डेल इंटीग्रल (Mordell integral) $H(y)$ के मोमेंट्स (moments) द्वारा परिभाषित है।
   • ये नोड्स और भार एक तीन-पद पुनरावृत्ति संबंध (three-term recurrence relation) के माध्यम से उत्पन्न ऑर्थोगोनल बहुपदों (orthogonal polynomials) का उपयोग करके संगणित किए जाते हैं।
4. सन्निकटन सूत्र (Approximation Formula): अंतिम सन्निकटन को इस प्रकार परिभाषित किया गया है:
   $$\zeta_p(s) = \sum_{n=1}^{N_t} n^{-s} + \chi(s) \sum_{n=1}^{N_t} n^{s-1} - \frac{(-1)^{N_t}}{2} \left[ I_{N_t+1/2, p}(s) + \chi(s) I_{N_t+1/2, p}(1-s) \right]$$
   जहाँ $I_{M,p}(s)$ समाकल पद का क्वाड्रैचर सन्निकटन है, जिसमें केवल प्रारंभिक फलन (elementary functions - घातांक, लघुगणक) और पूर्व-संगणित जटिल गुणांक शामिल हैं।

मुख्य योगदान

• सरलीकृत कार्यान्वयन (Simplified Implementation): प्रस्तावित सन्निकटन एक एकल पूर्णांक पैरामीटर $p$ और पूर्व-संगणित गुणांकों $\{\omega_{p,j}\}$ और $\{\lambda_{p,j}\}$ के एक सेट पर निर्भर करता है। रीमान-सीगल सूत्र के विपरीत, जिसमें उच्च-क्रम पदों के लिए विशेष फलनों के जटिल अवकलजों के मूल्यांकन की आवश्यकता होती है, $\zeta_p(s)$ के लिए केवल प्रारंभिक फलनों और गामा फलन (Gamma function) की आवश्यकता होती है।
• पूर्व-संगणित गुणांक (Precomputed Coefficients): लेखक $1 \le p \le 30$ के लिए पूर्व-संगणित गुणांक और $p=150$ तक चयनित मान प्रदान करते हैं। ये गुणांक मध्यम परिमाण के दिखाए गए हैं और इनमें तीव्र क्षय (rapid decay) प्रदर्शित होता है, जो स्थिर संख्यात्मक मूल्यांकन की सुविधा प्रदान करता है।
• अवकलज सन्निकटन (Derivative Approximation): शोध पत्र $\zeta'(s)$ को सन्निकट करने के लिए पद्धति का विस्तार करता है, जो $s$ के सापेक्ष सन्निकटन सूत्र का अवकलन करता है (इसे निश्चित $N_t$ पर मूल्यांकित किया जाता है), जिससे $N_t$ के फ्लोर फंक्शन (floor function) के अवकलन से जुड़ी विच्छिन्नता (discontinuities) से बचा जा सके।
• हाइब्रिड एल्गोरिदम (Hybrid Algorithm): $s$ के क्षेत्र के आधार पर $\zeta(s)$ की गणना के लिए एक रणनीति प्रस्तावित की गई है, जो प्रत्यक्ष योग (direct summation), यूलर-मैकलॉरिन योग और $\zeta_p(s)$ सन्निकटन के बीच स्विच करती है।

परिणाम

शोध पत्र $\zeta_p(s)$ की सटीकता के संबंध में व्यापक संख्यात्मक साक्ष्य प्रस्तुत करता है:

• क्रिटिकल स्ट्रिप में सटीकता: क्रिटिकल स्ट्रिप ($0 \le \sigma \le 1$) में संख्यात्मक प्रयोगों से पता चलता है कि जैसे-जैसे $p$ बढ़ता है, त्रुटि $|\zeta_p(s) - \zeta(s)|$ काफी कम हो जाती है।
  • $p=10$ के लिए, $t > 250$ होने पर त्रुटियाँ $< 10^{-15}$ हैं।
  • $p=50$ के लिए, $t > 4000$ होने पर त्रुटियाँ $< 10^{-100}$ हैं।
  • $p=150$ के लिए, $t > 6900$ होने पर त्रुटियाँ $< 10^{-300}$ हैं।
• त्रुटि पैटर्न (Error Patterns): त्रुटि एक विशिष्ट पैटर्न प्रदर्शित करती है जहाँ यह अंतरालों $[t_n, t_{n+1}]$ (जहाँ $t_n = 2\pi n^2$) के सिरों पर और इन अंतरालों के भीतर $4p$ समान रूप से वितरित बिंदुओं पर शून्य हो जाती है। यह व्यवहार विश्लेषणात्मक रूप से इस तथ्य द्वारा समझाया गया है कि प्रमुख त्रुटि पद $H_p(B(t)) - H(B(t))$ के समानुपाती है, जो क्वाड्रैचर नोड्स पर शून्य होता है।
• रीमान-सीगल के साथ तुलना: मध्यम $t$ (जैसे, $10^2 \le t \le 10^4$) के लिए, $\zeta_7(s)$ का सन्निकटन $K=5$ वाले रीमान-सीगल सन्निकटन ($\zeta^{RS}_5$) की तुलना में अधिक सटीक दिखाया गया है। हालांकि $\zeta^{RS}_5$ की स्पर्शोन्मुख क्षय दर ($O(t^{-13/4})$) $\zeta_7$ ($O(t^{-1/4})$) की तुलना में तेज़ है, फिर भी $\zeta_7$ के लिए बहुत कम पूर्व-संगठित गुणांकों (उच्च परिशुद्धता के लिए 15 बनाम ~100–170) की आवश्यकता होती है और इसे लागू करना आसान है।
• परिशुद्धता सीमाएँ (Precision Limits): डबल-प्रिसिजन कार्यान्वयन में, राउंडिंग एरर (rounding errors) $t > 200$ के लिए सन्निकटन त्रुटि पर हावी हो जाते हैं। हालांकि, क्वाड्रपल प्रिसिजन (quadruple precision) में, सन्निकटन बहुत बड़े $t$ (जैसे, $t > 10^4$ के साथ त्रुटियाँ $< 10^{-29}$) तक उच्च सटीकता बनाए रखते हैं।
• प्रदर्शन (Performance): हाइब्रिड एल्गोरिदम के MATLAB कार्यान्वयन को बिल्ट-इन MATLAB zeta फंक्शन की तुलना में सैकड़ों से हजारों गुना तेज़ बताया गया है।

महत्व और दावे

शोध पत्र का दावा है कि ये सन्निकटन बड़े $t$ के लिए क्रिटिकल स्ट्रिप और अन्य वर्टिकल स्ट्रिप्स में $\zeta(s)$ और $\zeta'(s)$ की उच्च-परिशुद्धता गणना को सक्षम बनाते हैं। इसका महत्व कार्यान्वयन की जटिलता और सटीकता के बीच के संतुलन में निहित है: यह विधि रीमान-सीगल सूत्र के एक "सरल" विकल्प के रूप में कार्य करती है जो उच्च-क्रम सुधार पदों के कम्प्यूटेशनल बोझ से बचते हुए प्रतिस्पर्धी या बेहतर सटीकता बनाए रखती है।

लेखक नोट करते हैं कि हालांकि इस कार्य में $\zeta_p(s)$ के लिए कठोर त्रुटि सीमाएँ प्रदान नहीं की गई हैं, लेकिन संख्यात्मक साक्ष्य व्यापक हैं। इस पद्धति को उच्च-परिशुद्धता गणनाओं के लिए एक व्यावहारिक उपकरण के रूप में प्रस्तुत किया गया है, जो मर्टेंस अनुमान (Mertens conjecture) के खंडन जैसी गणितीय खोजों के लिए महत्वपूर्ण हैं। शोध पत्र MATLAB, Python और Fortran में ओपन-सोर्स कार्यान्वयन प्रदान करके समाप्त होता है, जिससे अनुसंधान समुदाय द्वारा तत्काल उपयोग की सुविधा मिलती है।