Precision three-Dimensional Atmospheric Neutrino Flux Calculation Based on Honda Flux Model

होंडा फ्लक्स मॉडल पर आधारित सटीक त्रिविमीय (3D) वायुमंडलीय न्यूट्रिनो फ्लक्स गणना का तकनीकी सारांश

समस्या विवरण (Problem Statement)

वायुमंडलीय न्यूट्रिनो ($\nu_{\text{atm}}$) न्यूट्रिनो दोलन अध्ययनों, विसरित सुपरनोवा न्यूट्रिनो बैकग्राउंड (DSNB) खोजों, डार्क मैटर डिटेक्शन और न्यूक्लियॉन क्षय (nucleon decay) की खोजों के लिए महत्वपूर्ण पृष्ठभूमि (backgrounds) और प्रोब हैं। हालांकि HKKMS, Bartol और FLUKA जैसे त्रिविमीय (3D) फ्लक्स गणना ढांचे मौजूद हैं, वे वर्तमान और भविष्य के प्रयोगों के संदर्भ में तीन विशिष्ट सीमाओं का सामना करते हैं:

1. पुराने इनपुट मॉडल: पिछले गणनाओं में अक्सर पुराने प्राथमिक ब्रह्मांडीय किरण (primary cosmic ray) मॉडलों (जैसे AMS01-आधारित) और भू-चुंबकीय क्षेत्र मॉडलों (जैसे IGRF2010) पर भरोसा किया जाता है जो नवीनतम प्रयोगात्मक मापों को प्रतिबिंबित नहीं करते हैं।
2. लुप्त निम्न-ऊर्जा भौतिकी: किसी भी पिछले 3D फ्लक्स गणना में पृथ्वी के भीतर वायुमंडलीय म्यूऑन ($\mu_{\text{atm}}$) के प्रसार, उनके बाद के रुकने (stopping), और उनके क्षय (decay) या परमाणु कैप्चर (nuclear capture) को शामिल नहीं किया गया है। यह कमी 100 MeV से कम ऊर्जा सीमा के लिए महत्वपूर्ण है, जो DSNB और डार्क मैटर खोजों के लिए अत्यंत आवश्यक है।
3. अपूर्ण साइट तुलना: एक एकीकृत, अपडेटेड ढांचे का उपयोग करके विविध भू-चुंबकीय वातावरण (दक्षिण ध्रुव से लेकर निम्न-अक्षांश स्थलों तक) में फ्लक्स के व्यवस्थित अंतरों की व्यापक तुलना नहीं की गई है।

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

लेखक स्थापित HKKMS15 ढांचे पर आधारित एक व्यापक 3D वायुमंडलीय न्यूट्रिनो फ्लक्स गणना प्रस्तुत करते हैं, जिसमें कई प्रमुख पद्धतिगत अपडेट शामिल हैं:

• पृथ्वी के भीतर म्यूऑन का प्रसार: पहली बार, यह सिमुलेशन उन आवेशित द्वितीयक कणों (मुख्य रूप से म्यूऑन) को ट्रैक करता है जो पृथ्वी की पपड़ी (crust) या महासागर में प्रवेश करते हैं। यह ढांचा उनके ऊर्जा नुकसान, रुकने और उसके बाद के भाग्य को मॉडल करता है:
  • $\mu^+$ स्वतंत्र रूप से क्षय होता है ($\mu^+ \to e^+ \nu_e \bar{\nu}_\mu$)।
  • $\mu^-$ या तो कक्षा में क्षय होता है या परमाणु कैप्चर का अनुभव करता ($\mu^- + p \to \nu_\mu + n$)।
  • रुकने वाले माध्यम (stopping medium) की पहचान CRUST1.0 ग्लोबल क्रस्टल मॉडल का उपयोग करके की जाती है।
  • यह प्रक्रिया निम्न-ऊर्जा न्यूट्रिनो के पूरक फ्लक्स को उत्पन्न करती है, जहाँ परमाणु कैप्चर $\sim95$ MeV तक $\nu_\mu$ उत्पन्न करता है।
• अपडेटेड इनपुट मॉडल:
  • प्राथमिक ब्रह्मांडीय किरणें (Primary Cosmic Rays): AMS01-आधारित मॉडल को AMS02-आधारित मॉडल से बदल दिया गया है, जिसमें AMS02, BESS-polar, PAMELA, JACEE, RUNJOB और CREAM से प्राप्त डेटा शामिल है।
  • भू-चुंबकीय क्षेत्र: साइट-विशिष्ट रिजिडिटी कटऑफ (rigidity cutoffs) और म्यूऑन बेंडिंग की सटीक गणना के लिए इसे IGRF2020 मॉडल में अपडेट किया गया है।
  • हेड्रोनिक इंटरैक्शन (Hadronic Interactions): JAM और संशोधित DPMJET-III मॉडल्स के म्यूऑन-कैलिब्रेटेड संयोजन का उपयोग किया गया है, जिसे 1–100 GeV पर $\mu_{\text{atm}}$ फ्लक्स विसंगतियों को 5% के भीतर कम करने के लिए ट्यून किया गया है।
  • वायुमंडल: साइट-और-समय-निर्भर वायु घनत्व प्रोफाइल के लिए NRLMSISE-00 ग्लोबल वायुमंडलीय मॉडल का उपयोग किया गया है।
• सिमुलेशन का दायरा: यह गणना 10 MeV से $10^4$ GeV की ऊर्जा सीमा में सात डिटेक्टर साइटों (JUNO, Super-Kamiokande, CJPL, KM3NeT/ORCA, IceCube, DUNE, और TRIDENT) को कवर करती है। यह ट्रांसवर्स मोमेंटम और पृथ्वी की वक्रता को बनाए रखने के लिए पूर्ण गोलाकार ज्यामिति (spherical geometry) और 3D भू-चुंबकीय क्षेत्रों का उपयोग करती है।

प्रमुख योगदान (Key Contributions)

1. पृथ्वी के भीतर प्रसारित होने वाले म्यूऑन का प्रथम समावेश: यह शोध पत्र पृथ्वी के भीतर रुके हुए वायुमंडलीय म्यूऑन के योगदान को न्यूट्रिनो फ्लक्स में शामिल करता है, जो एक गुणात्मक रूप से नया स्रोत है जो सभी पूर्व 3D गणनाओं (HKKMS11/15, Bartol, FLUKA) में अनुपस्थित था।
2. ग्लोबल मल्टी-साइट फ्रेमवर्क: यह सात प्रमुख प्रयोगात्मक स्थानों के लिए सटीक, साइट-विशिष्ट फ्लक्स भविष्यवाणियां प्रदान करता है, जिससे दक्षिण ध्रुव (IceCube) से लेकर निम्न-अक्षांश स्थलों (TRIDENT) तक भू-चुंबकीय प्रभावों की सीधी तुलना संभव होती है।
3. निम्न-ऊर्जा सटीकता: यह इन साइटों के लिए 100 MeV से नीचे के पहले सटीक फ्लक्स परिणाम प्रदान करता है, एक ऐसा क्षेत्र जहाँ पहले विस्तृत 3D मॉडलिंग का अभाव था।
4. पुनः अनुमानित अनिश्चितताएं: म्यूऑन-प्रतिबंधित विधि (muon-constrained method) का उपयोग करके हेड्रोनिक अनिश्चितता का अनुमान अपडेट किया गया है, जो पिछले गणनाओं की तुलना में व्यवस्थित त्रुटियों (systematic errors) को काफी कम करता है।

परिणाम (Results)

1. $E_\nu > 100$ MeV के लिए फ्लक्स की विशेषताएं:

• साइट निर्भरता: भू-चुंबकीय कटऑफ के कारण $E_\nu < 10$ GeV पर महत्वपूर्ण साइट-निर्भर अंतर दिखाई देते हैं। उदाहरण के लिए, 1 GeV से नीचे IceCube पर $\nu_\mu$ फ्लक्स JUNO की तुलना में लगभग दोगुना है। 10 GeV से ऊपर, फ्लक्स सभी साइटों पर अभिसरित (converge) हो जाते हैं।
• जेनिथ और अज़ीमुथ (Zenith and Azimuth): फ्लक्स उत्पादन ज्यामिति और भू-चुंबकीय सुरक्षा द्वारा संचालित तीव्र जेनिथ-कोण निर्भरता (up/down और horizontal/vertical अनुपात) प्रदर्शित करता है। अज़ीमुथ विषमताएं (East-West प्रभाव) सब-GeV ऊर्जाओं पर प्रमुख हैं लेकिन उच्च ऊर्जाओं पर समाप्त हो जाती हैं।
• HKKMS15 के साथ तुलना: HKKMS15 की तुलना में 2%–10% का विचलन देखा गया है, जिसका मुख्य कारण अपडेटेड प्राथमिक ब्रह्मांडीय किरण स्पेक्ट्रल शेप (40 GeV से नीचे उच्च फ्लक्स, ऊपर कम) और पुन: कैलिब्रेटेड हेड्रोनिक इंटरैक्शन मॉडल है।

2. $E_\nu < 100$ MeV के लिए फ्लक्स की विशेषताएं (म्यूऑन प्रोपेगेशन योगदान):

• ग्लोबल एब्सोल्यूट इंक्रीमेंट: पृथ्वी के भीतर म्यूऑन प्रसार का योगदान एक वैश्विक प्रभाव है। फ्लक्स में पूर्ण वृद्धि ($\Delta \phi_\nu$) साइट-स्वतंत्र है, जो 10–50 MeV रेंज में चरम पर होती है।
• सापेक्ष वृद्धि (Relative Enhancement): चूंकि बेसलाइन फ्लक्स साइट के अनुसार भिन्न होता है (भू-चुंबकीय कटऑफ के कारण), इसलिए सापेक्ष वृद्धि विपरीत रूप से भिन्न होती है। कम बेसलाइन फ्लक्स वाली साइटें (जैसे JUNO) उच्च-फ्लक्स साइटों (जैसे IceCube) की तुलना में अधिक भिन्नात्मक (fractional) वृद्धि दिखाती हैं।
• **फ्लेवर पदानुक्रम (Flavor Hierarchy):ity: $\nu_e \approx \bar{\nu}_\mu > \nu_\mu > \bar{\nu}_e$ का पालन करता है। यह रुकने वाले म्यूऑन ($\mu^+ > \mu^-$) की चार्ज विषमता और परमाणु कैप्चर प्रक्रिया द्वारा संचालित है, जो $\sim95$ MeV तक ऊर्जा वाले $\nu_\mu$ का उत्पादन करता है, जो मिशेल एंडपॉइंट ($\sim52.8$ MeV) से काफी ऊपर है।

3. अनिश्चितता अनुमान (Uncertainty Estimation):

• म्यूऑन-प्रतिबंधित विधि का उपयोग करते हुए, हेड्रोनिक अनिश्चितता ($\delta_{\text{had}}$) काफी कम हो गई है। 1–10 GeV पर क्षैतिज दिशा के लिए, $\delta_{\text{had}}$ $\nu_\mu/\bar{\nu}_\mu$ के लिए $\sim3\%$ और $\nu_e/\bar{\nu}_e$ के लिए $\sim4\%$ तक कम हो जाता है।
• कुल सापेक्ष अनिश्चितता ($\delta_{\text{tot}}$) 1–10 GeV रेंज (क्षैतिज) में $\sim5\%$–$7\%$ तक कम हो जाती है, जो पिछले अनुमानों की तुलना में एक उल्लेखनीय सुधार है।

महत्व (Significance)

यह शोध दावा करता है कि ये परिणाम व्यापक भौतिकी कार्यक्रमों के लिए आवश्यक, सटीक इनपुट प्रदान करते हैं:

• न्यूट्रिनो दोलन (Neutrino Oscillations): 3D फ्लक्स भविष्यवाणियां और कम हुई अनिश्चितताएं वायुमंडलीय न्यूट्रिनो प्रयोगों (JUNO, Super-Kamiokande, Hyper-Kamiokande, DUNE, KM3NeT/ORCA, IceCube) के लिए बेहतर इनपुट के रूप में कार्य करती हैं।
• DSNB खोज: निम्न-ऊर्जा फ्लक्स गणनाएं (म्यूऑन-प्रसार घटक सहित) 10–30 MeV विंडो में DSNB खोजों के लिए प्रमुख अपरिहार्य बैकग्राउंड (irreducible background) को परिभाषित करती हैं।
• डार्क मैटर खोज:
  • अप्रत्यक्ष (Indirect): अपडेटेड फ्लक्स WIMP विलोपन (IceCube) और हल्के डार्क मैटर विलोपन (JUNO, Borexino, Super-K) की खोजों के लिए बैकग्राउंड लक्षण वर्णन में सुधार करते हैं।
  • प्रत्यक्ष (Direct): सटीक निम्न-ऊर्जा भविष्यवाणियां उस "न्यूट्रिनो फॉग" (neutrino fog) को सीमित करती हैं जो प्रत्यक्ष डार्क मैटर डिटेक्टरों (LZ, XENONnT, PandaX-4T) की संवेदनशीलता को सीमित करता है।
• न्यूक्लियॉन क्षय (Nucleon Decay): बेहतर-नियंत्रित वायुमंडलीय न्यूट्रिनो बैकग्राउंड वर्तमान और अगली पीढ़ी के डिटेक्टरों में न्यूक्लियॉन क्षय खोजों के लिए संवेदनशीलता सीमाओं को सीधे प्रभावित करता है।

लेखक निष्कर्ष निकालते हैं कि जबकि वर्तमान गणना एक मजबूत आधार प्रदान करती है, भविष्य के सुधार—जैसे स्थानीय पर्वत प्रोफाइल, समय-निर्भर सौर मॉड्यूलेशन, और NA61/SHINE तथा LHCf डेटा का उपयोग करके अपडेटेड हेड्रोनिक मॉडल को शामिल करना—विशेष रूप से निम्न-से-मध्यम ऊर्जा क्षेत्रों में व्यवस्थित अनिश्चितताओं को और कम करेंगे।