🔬 applied physics

Improved capabilities of the TurboGAP code for radiation induced cascade simulations: an illustration with silicon

यह शोध पत्र सिलिकॉन में विकिरण-प्रेरित कैस्केड सिमुलेशन के लिए टर्बो गैप (TurboGAP) कोड की उन्नत क्षमताओं को प्रस्तुत करता है, जिसमें दोष क्लस्टरिंग और दस लाख परमाणुओं तक की प्रणालियों में आयन-बीम मिक्सिंग का सटीक पूर्वानुमान लगाने के लिए टू-टेम्परेचर मॉडल, एडेप्टिव टाइमस्टेप्स और यथार्थवादी इलेक्ट्रॉनिक स्टॉपिंग पावर शामिल हैं।

मूल लेखक: Uttiyoarnab Saha, Ali Hamedani, Miguel A. Caro, Andrea E. Sand

प्रकाशित 2026-07-08

📖 5 मिनट में पढ़ें🧠 गहराई से पढ़ें

CC BY 4.0

मूल लेखक: Uttiyoarnab Saha, Ali Hamedani, Miguel A. Caro, Andrea E. Sand

मूल पेपर CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/) के तहत लाइसेंस किया गया है। ✨ नीचे दिए गए पेपर की यह व्याख्या AI से तैयार की गई है। इसे लेखकों ने न तो लिखा है, न इसका समर्थन किया है। तकनीकी सटीकता के लिए मूल पेपर देखें। पूरा डिस्क्लेमर पढ़ें

एक कंप्यूटर सिमुलेशन की कल्पना एक विशाल, हाई-स्पीड मूवी कैमरे के रूप में करें जो यह फिल्मा रहा है कि क्या होता है जब एक अकेला, अत्यंत तीव्र बुलेट (एक आयन) सिलिकॉन के एक ब्लॉक से टकराता है। यह "बुलेट" वास्तव में विकिरण (रेडिएशन) का एक कण है, और जब यह सिलिकॉन से टकराता है, तो यह परमाणुओं को बिलियर्ड बॉल्स की तरह इधर-उधर धकेलते हुए एक अराजक श्रृंखला अभिक्रिया (चेन रिएक्शन) पैदा करता है। इसे "रेडिएशन कैस्केड" कहा जाता है।

लंबे समय तक, वैज्ञानिक इस फिल्म को सटीक रूप से फिल्माने के लिए संघर्ष करते रहे हैं। उनके पास एक विकल्प था:

सुपर-डिटेल्ड कैमरा: अत्यंत सटीक लेकिन इतना धीमा कि यह केवल सिलिकॉन के एक बहुत ही छोटे से हिस्से को पल भर के लिए फिल्मा सकता था।
फास्ट कैमरा: यह लंबे समय तक सिलिकॉन के एक बड़े हिस्से को फिल्मा सकता था, लेकिन भौतिकी अक्सर बहुत सरल थी, जिससे परमाणुओं के बीच होने वाली बातचीत के महत्वपूर्ण विवरण छूट जाते थे।

यह पेपर TurboGAP नामक एक सॉफ्टवेयर टूल के अपग्रेड को पेश करता है। TurboGAP को एक हाई-परफॉर्मेंस इंजन के रूप में समझें जो वैज्ञानिकों को "सुपर-डिटेल्ड" भौतिकी (जिसे मशीन लर्निंग इंटरएटोमिक पोटेंशियल कहा जाता है) का उपयोग करने की अनुमति देता है, लेकिन इसे "फास्ट कैमरा" की गति पर चलाने में सक्षम बनाता है। लेखकों ने इस इंजन को और बेहतर बनाने के लिए इसमें तीन नई विशेषताएं जोड़ी हैं।

यहाँ इन तीन नई विशेषताओं को सरल उपमाओं (analogies) के माध्यम से समझाया गया है:

1. "इलेक्ट्रॉनिक ब्रेक" (इलेक्ट्रॉनिक स्टॉपिंग)

जब एक तेज़ गति वाला कण किसी पदार्थ के माध्यम से गुजरता है, तो वह केवल अन्य परमाणुओं से टकराता ही नहीं है; बल्कि वह इलेक्ट्रॉनों को भी अपने साथ खींचता है, जिससे ऊर्जा का नुकसान होता है—ठीक वैसे ही जैसे कोई कार गाढ़े कीचड़ से होकर गुजर रही हो।

पुराना तरीका (फ्रिक्शन मॉडल): कल्पना करें कि आप एक फिसलते हुए पक (puck) पर ब्रेक लगा रहे हैं, लेकिन ब्रेक तभी लगाते हैं जब पक एक विशिष्ट गति से तेज चल रहा हो। यदि पक उस गति से नीचे धीमा हो जाता है, तो ब्रेक जादुई रूप से बंद हो जाता है। समस्या यह थी कि वैज्ञानिकों को उस "जादुई गति" का अनुमान लगाना पड़ता था। यदि उनका अनुमान गलत होता, तो सिमुलेशन गलत परिणाम देता।
नया तरीका (EPH मॉडल): लेखकों ने एक स्मार्ट सिस्टम जोड़ा है। एक साधारण ऑन/ऑफ ब्रेक के बजाय, उन्होंने एक "स्मार्ट थर्मोस्टेट" बनाया है जो चलते हुए परमाणुओं को इलेक्ट्रॉनों के एक विशाल, अदृश्य हीट बाथ (ऊष्मा कुंड) से जोड़ता है। यह प्रणाली स्वचालित रूप से उस स्थान पर इलेक्ट्रॉनों की भीड़ के आधार पर ब्रेकिंग फोर्स को एडजस्ट करती है। इसे किसी "जादुई गति" वाले कटऑफ की आवश्यकता नहीं है; यह स्वाभाविक रूप से काम करता है, जिससे परमाणु अधिक यथार्थवादी तरीके से धीमे हो जाते हैं।

2. "स्मार्ट टाइमर" (एडैप्टिव टाइमस्टेप)

एक सिमुलेशन में, कंप्यूटर यह गणना करने के लिए कि परमाणु आगे कहाँ जाएंगे, समय के छोटे "कदम" (steps) लेता है।

समस्या: जब परमाणु शांति से उड़ रहे होते हैं, तो कंप्यूटर बड़े कदम ले सकता है। लेकिन जब दो परमाणु उच्च गति से आपस में टकराते हैं, तो बल (forces) तुरंत बदल जाते हैं। यदि कंप्यूटर बड़े कदम लेना जारी रखता है, तो वह टक्कर को मिस कर देगा या उसे गलत तरीके से कैलकुलेट करेगा।
समाधान: नए TurboGAP में एक "स्मार्ट टाइमर" है। यह परमाणुओं पर नज़र रखता है। यदि वे धीरे चल रहे हैं, तो यह समय बचाने के लिए बड़े कदम लेता है। जैसे ही परमाणु टकराने के करीब आते हैं, टाइमर तुरंत कदमों को सूक्ष्म (microscopic) बना देता है, जिससे यह सुनिश्चित होता है कि टक्कर को पूरी तरह से कैप्चर किया गया है। एक बार जब अराजकता शांत हो जाती है, तो यह कदमों की गति को फिर से बढ़ा देता है।

3. "कूलिंग बॉर्डर" (परमाणुओं का समूहीकरण)

जब एक रेडिएशन कैस्केड होता है, तो यह सिमुलेशन के केंद्र में भारी मात्रा में गर्मी पैदा करता है। यदि यह गर्मी फंसी रह जाती है, तो कंप्यूटर में सिलिकॉन का पूरा ब्लॉक पिघल जाएगा, जो अधिकांश प्रयोगों के लिए यथार्थवादी नहीं है।

समाधान: लेखकों ने सिमुलेशन बॉक्स के बिल्कुल किनारे पर परमाणुओं को समूहबद्ध करने का एक तरीका जोड़ा है। वे एक कूलिंग जैकेट की तरह कार्य करते हैं, जो लगातार केंद्र से अतिरिक्त गर्मी को सोखते रहते हैं और किनारों को एक स्थिर तापमान पर रखते हैं। यह टक्कर से उत्पन्न होने वाली "हीट वेव" को स्वाभाविक रूप से बाहर बहने देता है, ठीक वैसे ही जैसे गर्मी एक गर्म पैन से हवा में बाहर निकलती है।

उन्होंने क्या पाया?

इन नए उपकरणों का उपयोग करते हुए, लेखकों ने दस लाख परमाणुओं (इस प्रकार की भौतिकी के लिए एक बड़ी संख्या) वाले सिलिकॉन ब्लॉक का सिमुलेशन किया और इसमें 10,000 इलेक्ट्रॉन-वोल्ट की ऊर्जा वाला एक कण छोड़ा।

"मैजिक स्पीड" का मुद्दा: उन्होंने पाया कि पुराना "फ्रिक्शन मॉडल" (जिसमें अनुमान लगाने योग्य स्पीड लिमिट थी) अलग-अलग स्पीड लिमिट चुनने पर अलग-अलग परिणाम देता था। नया "स्मार्ट थर्मोस्टेट" (EPH मॉडल) बिना किसी अनुमान के सुसंगत परिणाम देता है।
डिफेक्ट क्लस्टर्स: जब रेडिएशन के बाद परमाणु स्थिर होते हैं, तो वे "डिफेक्ट्स" (छेद या अतिरिक्त परमाणु) बनाते हैं। नए मॉडल ने दिखाया कि ये डिफेक्ट्स पुराने मॉडल की तुलना में बड़े, अधिक यथार्थवादी समूहों में इकट्ठा होते हैं।
मिक्सिंग: उन्होंने मापा कि टक्कर द्वारा परमाणुओं को कितना "मिश्रित" किया गया। उनके नए, यथार्थवादी मॉडल के परिणाम पिछले सिमुलेशन की तुलना में वास्तविक प्रयोगशाला प्रयोगों से कहीं बेहतर मेल खाते हैं।

मुख्य निष्कर्ष (The Bottom Line)

लेखकों ने केवल सॉफ्टवेयर को तेज़ नहीं बनाया; उन्होंने इसे स्मार्ट बनाया है। इलेक्ट्रॉनों को ऊर्जा हानि संभालने के एक यथार्थवादी तरीके, टकराव के लिए एक स्मार्ट टाइमर, और एक बेहतर कूलिंग सिस्टम को जोड़कर, वे अब उच्च सटीकता के साथ सिलिकॉन में बड़े रेडिएशन इवेंट्स का सिमुलेशन कर सकते हैं। यह वैज्ञानिकों को यह समझने में मदद करता है कि विकिरण के तहत सामग्रियां कैसे खराब होती हैं, बिना खेल के नियमों का अनुमान लगाए।

तकनीकी सारांश: विकिरण-प्रेरित कैस्केड सिमुलेशन के लिए टर्बोगैप (TurboGAP) कोड की उन्नत क्षमताएं

समस्या विवरण
विकिरण क्षति (radiation damage) के सटीक परमाणु स्तर के सिमुलेशन के लिए पोटेंशियल एनर्जी सरफेस (PES) की भौतिक सटीकता और गणनात्मक दक्षता के बीच संतुलन बनाना आवश्यक है। जबकि मशीन लर्निंग इंटरएटोमिक पोटेंशियल (MLIPs), जैसे कि गॉसियन एप्रोक्सिमेशन पोटेंशियल (GAP), एम्पिरिकल पोटेंशियल की तुलना में बेहतर सटीकता प्रदान करते हैं, उच्च प्राथमिक प्रहार परमाणु (PKA) ऊर्जाओं के लिए उनका अनुप्रयोग अक्सर गणनात्मक लागत के कारण सीमित होता है। इसके अलावा, मानक मॉलिक्यूलर डायनेमिक्स (MD) कोड में अक्सर वास्तविक विकिरण क्षति मॉडलिंग के लिए आवश्यक आवश्यक विशेषताएं नहीं होती हैं: उच्च-ऊर्जा टकरावों के दौरान अनुकूली टाइमस्टेप (adaptive timesteps) को संभालने की क्षमता, सिमुलेशन सेल सीमाओं का प्रभावी कूलिंग, और इलेक्ट्रॉनिक ऊर्जा अपव्यय (इलेक्ट्रॉनिक स्टॉपिंग) का सटीक मॉडलिंग। पारंपरिक घर्षण-आधारित इलेक्ट्रॉनिक स्टॉपिंग (FES) मॉडल अक्सर मनमाने गतिज ऊर्जा कटऑफ (kinetic energy cutoffs) पर निर्भर करते हैं, जो दोष सूक्ष्म संरचनाओं (defect microstructures) को महत्वपूर्ण रूप से प्रभावित कर सकते हैं। एक ऐसे सॉफ्टवेयर फ्रेमवर्क की आवश्यकता है जो $10^6$ परमाणुओं तक के सिस्टम में 10 keV तक की PKA ऊर्जाओं के लिए कैस्केड को सिम्युलेट करने हेतु टर्बोगैप की दक्षता को उन्नत अपव्यय मॉडलों, जैसे कि इलेक्ट्रॉन-फोनोन (EPH) कपलिंग मॉडल के साथ जोड़ सके।

कार्यप्रणाली
लेखकों ने तीन महत्वपूर्ण मॉड्यूल को लागू करके टर्बोगैप सॉफ्टवेयर पैकेज को उन्नत किया है:

इलेक्ट्रॉनिक ऊर्जा अपव्यय मॉडल:
- घर्षण-आधारित इलेक्ट्रॉनिक स्टॉपिंग (FES): इसे पूर्व-गणना की गई इलेक्ट्रॉनिक स्टॉपिंग पावर (ESP) तालिकाओं (जैसे, SRIM से) से प्राप्त वेग-निर्भर घर्षण बल के रूप में लागू किया गया है। एक उपयोगकर्ता-निर्धारित गतिज ऊर्जा कटऑफ ( $E_{el\_cut}$ ) यह निर्धारित करती है कि किन परमाणुओं पर यह बल लगेगा।
- इलेक्ट्रॉन-फोनोन (EPH) मॉडल: इसे एक सामान्यीकृत टू-टेम्परेचर मॉलिक्यूलर डायनेमिक्स (TTMD) फ्रेमवर्क में एकीकृत किया गया है। यह मॉडल इलेक्ट्रॉन घनत्व-निर्भर मापदंडों के माध्यम से परमाणु और इलेक्ट्रॉनिक उपप्रणालियों को जोड़ता है। यह परमाणु गतिशीलता के साथ इलेक्ट्रॉनिक तापमान ( $T_e$ ) के लिए ऊष्मा प्रसार समीकरण (heat diffusion equation) को हल करता है। परमाणु गति एडियाबेटिक बलों, सापेक्ष वेगों पर निर्भर घर्षण बलों और फ्लक्चुएशन-डिसिपेशन प्रमेय से प्राप्त स्टोकेस्टिक रैंडम बलों द्वारा नियंत्रित होती है। यह दृष्टिकोण मनमाने गतिज ऊर्जा कटऑफ से बचता है।
अनुकूली टाइमस्टेप गणना (Adaptive Timestep Calculation): निकट टकरावों के दौरान तीव्र पोटेंशियल एनर्जी ग्रेडिएंट्स को पकड़ने के लिए, टाइमस्टेप को दो मानदंडों के आधार पर गतिशील रूप से पुनर्गणना किया जाता है: किसी भी परमाणु के अधिकतम अनुमत विस्थापन ( $x_{max}$ ) और अधिकतम गतिज ऊर्जा में परिवर्तन ( $e_{max}$ )। इन मानदंडों से प्राप्त टाइमस्टेप्स में से न्यूनतम को चुना जाता है।
परमाणु समूहन और सीमा कूलिंग (Atom Grouping and Border Cooling): एक समूहन तंत्र (grouping mechanism) थर्मोस्टेट लागू करने के लिए स्थानिक मानदंडों (जैसे, एक बॉर्डर लेयर) के आधार पर परमाणुओं के चयन की अनुमति देता है। यह कैस्केड से अतिरिक्त ऊष्मा को बल्क सामग्री में विसर्जित करने में सक्षम बनाता है, जिससे वांछित सिमुलेशन तापमान बना रहता है।

पोटेंशियल विकास
सिलिकॉन (Si) के लिए एक विशिष्ट GAP MLIP को 2-body, 3-body, और turboSOAP डिस्क्रिप्टर (बेहतर दक्षता के लिए एक अनुकूलित वेरिएंट) का उपयोग करके पुन: प्रशिक्षित किया गया था। पोटेंशियल को विभिन्न Si संरचनाओं के डेटासेट पर प्रशिक्षित किया गया था, जिसमें तरल और अमोर्फस चरणों के लिए विशिष्ट रेगुलराइजेशन पैरामीटर शामिल थे। टकरावों के दौरान उच्च-ऊर्जा रिपल्सिव इंटरैक्शन को सटीक रूप से मॉडल करने के लिए, यूनिवर्सल ज़िगलर-बियर्सैक-लिटमर (ZBL) पोटेंशियल को GAP मॉडल में एक बाहरी पोटेंशियल के रूप में जोड़ा गया था, जो कम अंतर-परमाणु दूरियों पर सक्रिय होता है।

प्रमुख परिणाम

प्रदर्शन और स्केलिंग: नए मॉड्यूल का परीक्षण हाई-परफॉर्मेंस कंप्यूटिंग (HPC) सेटअप पर किया गया। $10^6$ परमाणुओं के सिस्टम के लिए 600 प्रोसेसरों के साथ स्केलिंग परीक्षणों ने प्रदर्शित किया कि लागू किए गए मॉड्यूल मूल टर्बोगैप कोड की गणनात्मक दक्षता बनाए रखते हैं।
कोड सत्यापन: LAMMPS (USER-EPH प्लगइन और मानक GAP कार्यान्वयन के साथ) के साथ तुलनाओं ने 100 eV PKA कैस्केड के लिए ऊष्मप्रवैगिकी गुणों (तापमान, दबाव, ऊर्जा) और संचयी इलेक्ट्रॉनिक ऊर्जा हानि के विकास में उत्कृष्ट सहमति दिखाई।
दोष उत्पादन और क्लस्टरिंग (2 keV PKA):
- EPH मॉडल ने FES मॉडलों (10 eV और 1 eV कटऑफ के साथ) की तुलना में अधिक जीवित दोषों (~33) की भविष्यवाणी की।
- EPH मॉडल के परिणामस्वरूप बड़े दोष क्लस्टर (कुछ 100 इंटरस्टिशियल्स से अधिक) बने, जबकि FES मॉडल ने आम तौर पर 100 इंटरस्टिशियल्स से कम के क्लस्टर उत्पन्न किए।
- FES मॉडल के परिणाम गतिज ऊर्जा कटऑफ के मनमाने चुनाव के प्रति संवेदनशील थे, जबकि EPH मॉडल के परिणाम विभिन्न इलेक्ट्रॉनिक उपप्रणाली सेटिंग्स के बावजूद सुदृढ़ (robust) थे।
- EPH मॉडल ने कैस्केड के ठंडा होने के बाद भी इलेक्ट्रॉनिक सिस्टम में निरंतर ऊर्जा हस्तांतरण की अनुमति दी, जबकि Fस (FES) मॉडल में ऊर्जा हस्तांतरण तब रुक गया जब परमाणु कटऑफ ऊर्जा से नीचे गिर गए।
आयन-बीम मिक्सिंग:
- EPH मॉडल और GAP MLIP का उपयोग करके गणना किए गए मिक्सिंग मान ( $Q$ ) Ne, Ar, Kr, और Xe आयनों के लिए Si के प्रायोगिक डेटा के साथ अच्छी सहमति दिखाते हैं।
- GAP MLIP की भविष्यवाणियां क्लासिकल स्टिलिंगर-वेबर (SW) पोटेंशियल का उपयोग करने वालों की तुलना में काफी भिन्न (2-4 गुना) थीं, जो सटीक PES प्रतिनिधित्व के महत्व को उजागर करती हैं।
- EPH मॉडल के मिक्सिंग मान FES परिणामों के तुलनीय थे, लेकिन यह बिना किसी मनमाने गतिज ऊर्जा कटऑफ के प्राप्त किया गया।

महत्व और दावे
यह शोध दावा करता है कि टर्बोगैप में EPH मॉडल, अनुकूली टाइमस्टेप और बॉर्डर कूलिंग का एकीकरण विकिरण क्षति सिमुलेशन की क्षमता को महत्वपूर्ण रूप से बढ़ाता है। टर्बोगैप की गणनात्मक दक्षता को GAP MLIPs की सटीकता के साथ जोड़कर, लेखक $10^6$ परमाणुओं तक के सिस्टम में 10 keV तक की PKA ऊर्जाओं के साथ कैस्केड को सफलतापूर्वक सिम्युलेट करने में सफल रहे हैं।

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

1. "इलेक्ट्रॉनिक ब्रेक" (इलेक्ट्रॉनिक स्टॉपिंग)

2. "स्मार्ट टाइमर" (एडैप्टिव टाइमस्टेप)

3. "कूलिंग बॉर्डर" (परमाणुओं का समूहीकरण)

उन्होंने क्या पाया?

मुख्य निष्कर्ष (The Bottom Line)

इसी तरह के और पेपर