A Unified Geometric Model of Ring Strain: From Carbon to Nitrogen and the Role of Lone-Pair Coulomb Repulsion

तकनीकी सारांश: वलय तनाव (Ring Strain) का एक एकीकृत ज्यामितीय मॉडल: कार्बन से नाइट्रोजन तक और एकाकी-युग्म कूलम्ब प्रतिकर्षण (Lone-Pair Coulomb Repulsion) की भूमिका

समस्या विवरण
वलय तनाव ऊर्जा (Ring strain energy) चक्रीय यौगिकों की स्थिरता और अभिक्रियाशीलता का एक मौलिक निर्धारक है। वलय तनाव के लिए एक एकीकृत मात्रात्मक सिद्धांत जो विभिन्न मुख्य-समूह तत्वों में तनाव की व्याख्या कर सके, अभी भी दुर्लभ बना हुआ है। पिछले दृष्टिकोणों में सीमाएँ हैं: आणविक यांत्रिकी बल क्षेत्र (molecular mechanics force fields) क्वांटम यांत्रिक अवलोकन से सीधे जुड़े नहीं होते हैं; प्रारंभिक इलेक्ट्रॉनिक मॉडल (जैसे, "बनाना बॉन्ड्स") ने अंतर्दृष्टि प्रदान की लेकिन वे सामान्यीकरण योग्य भविष्यवक्ता संकेतकों (descriptors) से रहित थे; और आधुनिक मशीन लर्निंग मॉडल अक्सर भौतिक व्याख्यात्मकता के बिना "ब्लैक बॉक्स" के रूप में कार्य करते हैं। यह अध्ययन इन अंतरालों को संबोधित करता है, जिसका उद्देश्य स्थिर तत्वों के लिए वलय तनाव की भविष्यवाणी करने हेतु भौतिक रूप से व्याख्यात्मक संकेतकों का एक न्यूनतम सेट खोजना और स्थिर द्वितीय-आवर्त होमोटमैटिक (homoatomic) वलयों के स्थायली परिदृश्य को मानचित्रित करना है।

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

• स्थिरता स्क्रीनिंग (Stability Screening): Be, B, C, N, और O वलयों ($n=3–10$) के प्रारंभिक ज्यामिति को B3LYP-D3(BJ)/def2-SVP स्तर पर अनुकूलित (optimize) किया गया। स्थिरता को बिना विघटन के एक चक्रीय स्थानीय न्यूनतम (local minimum) तक अभिसरण (convergence) के रूप में परिभाषित किया गया।
• कॉन्फ़ॉर्मेशनल सैंपलिंग (Conformational Sampling): कार्बन (C3–C14) और नाइट्रोजन (N3–N14) वलयों के लिए, CREF/GFN2-xTB का उपयोग करके कॉन्फ़ॉर्मर एन्सेम्बल बनाए गए। एन्सेम्बल-औसत गिब्स मुक्त ऊर्जा प्राप्त करने के लिए बहु-स्तरीय DFT फनल (SVP $\to$ TZVPP $\to$ QZVPP) के माध्यम से कॉन्फ़ॉर्मर्स को परिष्कृत किया गया।
• उच्च-परिशुद्धता ऊर्जा गणना: अंतिम गिब्स मुक्त ऊर्जाएं एक संयुक्त योजना का उपयोग करके प्राप्त की गईं: B3LYP-D3(BJ)/def2-QZVPP से थर्मल/एंट्रोपिक योगदान द्वारा सुधारा गया DLPNO-CCSD(T)/def2-TZVPP सिंगल-पॉइंट ऊर्जा। वलय तनाव ऊर्जा ($E_{strain}$) को आइसोडैस्मिक हाइड्रोजन-स्थानांतरण अभिक्रियाओं के माध्यम से गणना किया गया।
• डिस्क्रिप्टर निष्कर्षण (Descriptor Extraction): अनुकूलित ज्यामिति से ज्यामितीय डिस्क्रिप्टर ($\Delta\theta^2$, $\Delta R^2$) प्राप्त किए गए। इलेक्ट्रॉनिक डिस्क्रिप्टर (ADCH चार्ज, मल्टीपोल मोमेंट्स) Multiwfn द्वारा संगणित किए गए। नाइट्रोजन वलयों के लिए एक विशिष्ट एकाकी-युग्म कूलम्ब प्रतिकर्षण ($E_{coul}$) पद को सूत्रबद्ध किया गया।
• सांख्यिकीय विश्लेषण: CCSD(T) बेंचमार्क ऊर्जाओं के विरुद्ध पूर्ण-उपसमुच्चय ऑर्डिनरी लीस्ट स्क्वायर्स (OLS) प्रतिगमन किया गया। मॉडल चयन के लिए $R^2$, लिव-वन-आउट क्रॉस-वैलिडेशन (LOOCV) RMSE, वेरिएंस इन्फ्लेशन फैक्टर (VIF), और AIC/BIC मानदंडों का उपयोग किया गया।
• इलेक्ट्रॉनिक लक्षण वर्णन: N5 के असामान्य व्यवहार की जांच CASSCF(6,5) गणनाओं और $\sigma$-होमोएरोमैटिसिटी ( $\sigma$-homoaromaticity) की जांच के लिए न्यूक्लियस-स्वतंत्र केमिकल शिफ्ट (NICS) गणनाओं द्वारा की गई।

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

1. द्वितीय-आवर्त वलयों का स्थिरता द्वैत (Stability Dichotomy):
   कम्प्यूटेशनल स्क्रीनिंग एक तीव्र स्थिरता सीमा प्रकट करती है। केवल कार्बन (C) और नाइट्रोजन (N) ही एक विस्तृत आकार सीमा (C3–C30+ और N3–N14+) में स्थिर एकल-बंधित होमोटमैटिक वलयों का समर्थन करते हैं।

   • बेरिलियम (Be): इलेक्ट्रॉन की कमी के कारण वलय इलेक्ट्रॉन-न्यूनता वाले पिंजरों (cages) में ढह जाते हैं।
   • बोरोन (B): B4 से बड़े वलय बहु-केंद्र बंधन अस्थिरता के कारण अस्थिर होते हैं।
   • ऑक्सीजन (O): अत्यधिक एकाकी-युग्म प्रतिकर्षण द्वारा संचालित पेरोक्सी श्रृंखला निर्माण के कारण O4 से परे के वलय विघटित हो जाते हैं।
     इस स्थिरता को एक समरूपता-मिलान सिद्धांत (symmetry-matching principle) द्वारा तर्कसंगत बनाया गया है, जहाँ C और N ऐसी इलेक्ट्रॉनिक विन्यास (टेट्राहेड्रल/पिरामिडल $sp^3$) रखते हैं जो चक्रीय टोपोलॉजी के अनुकूल हैं, जबकि Be, B, और O वैकल्पिक संरचनाओं को प्राथमिकता देते हैं।

2. कार्बन वलयों के लिए एकीकृत ज्यामितीय मॉडल:
   कार्बन वलयों (C3–C14) के लिए, बॉन्ड कोण विचलन ($\Delta\theta^2$) और बॉन्ड लंबाई विरूपण ($\Delta R^2$) पर आधारित एक शुद्ध ज्यामितीय दो-पैरामीटर मॉडल वलय तनाव के विचरण का 98.2% ($R^2 = 0.982$, RMSE = 0.36 kcal mol$^{-1}$ atom$^{-1}$) कैप्चर करता है। इलेक्ट्रॉनिक डिस्क्रिप्टर को शामिल करने से नगण्य सुधार होता है, जो पुष्टि करता है कि कार्बन वलय तनाव मुख्य रूप से ज्यामितीय कारकों द्वारा नियंत्रित होता है।

3. नाइट्रोजन वलयों में एकाकी-युग्म प्रतिकर्षण की भूमिका:
   नाइट्रोजन वलयों के लिए, केवल ज्यामितीय ढांचा ही 92.5% विचरण की व्याख्या करता है। हालांकि, एक भौतिक रूप से व्युत्पन्न एकाकी-युग्म कूलम्ब प्रतिकर्षण पद ($E_{coul}$) को शामिल करने से मॉडल में $R^2 = 0.968$ (RMSE = 0.64 kcal mol$^{-1}$ atom$^{-1}$) तक महत्वपूर्ण सुधार होता है। यह पुष्टि करता है कि नाइट्रोजन वलयों में निरंतर अवशिष्ट तनाव (जो कार्बन की तुलना में अधिक रहता है) नाइट्रोजन के एकाकी-युग्मों के बीच प्रतिकर्षण अंतःक्रियाओं द्वारा संचालित है।

4. N5 विसंगति और $\sigma$-होमोएरोमैटिसिटी:
   N5 को एक सांख्यिकीय आउटलियर के रूप में पहचाना गया है जिसका अवशिष्ट (residual) $-4.0$ kcal mol$^{-1}$ है, जो अप्रत्याशित स्थिरीकरण का संकेत देता है। इसे $\sigma$-होमोएरोमैटिसिटी के कारण माना जाता है: पाँच इन-प्ले नाइट्रोजन एकाकी-युग्म $\sigma$-फ्रेमवर्क के माध्यम से विस्थानीकृत (delocalize) होते हैं, जिससे एक डायट्रोपिक रिंग करंट (N-I-C-S(0) = $-15.18$ ppm) बनता है। यह प्रभाव एकाकी-युग्म प्रतिकर्षण का मुकाबला करता है, जो बेंजीन की $\pi$-एरोमैटिसिटी से भिन्न है, जैसा कि आणविक तल के ऊपर NICS(1) के माध्यम से रिंग करंट के तीव्र क्षीणन से सिद्ध होता है।

महत्व और प्रस्तावित ढांचा
यह शोध एक वलय तनाव के स्तरित सिद्धांत (layered theory of ring strain) का प्रस्ताव करता है जो ज्यामितीय, इलेक्ट्रोस्टैटिक और क्वांटम एरोमैटिक प्रभावों को एक एकल मात्रात्मक ढांचे में एकीकृत करता है:
$$E_{strain} = E_{geom}(\Delta\theta^2, \Delta R^2) + E_{LP}(E_{coul}) + E_{\sigma-aromatic}(NICS)$$

• $E_{geom}$: सभी मुख्य-समूह होमोटमैटिक वलयों के लिए एक सार्वभौमिक ज्यामितीय ढांचा, जो बायर के कोण तनाव (Baeyer's angle strain) को पुनः प्राप्त करता है और बॉन्ड-लंबाई विरूपण को जोड़ता है।
• $E_{LP}$: एकाकी-युग्मों वाले तत्वों (पनिकोजेन्स) के लिए आवश्यक एक एकाकी-युग्म कूलम्ब सुधार, जो नाइट्रोजन के लिए प्रमुख गैर-ज्यामितीय कारक बन जाता है।
• $E_{\sigma-aromatic}$: विशिष्ट मामलों (जैसे N5) के लिए एक गैर-योगात्मक क्वांटम सुधार जहाँ ज्यामितीय और इलेक्ट्रॉनिक स्थितियाँ थ्रू-बॉन्ड एकाकी-युग्म विस्थानीकरण को सक्षम करती हैं।

यह कार्य यह स्थापित करता है कि जबकि वलय तनाव कार्बन के लिए मौलिक रूप से ज्यामितीय है, इसके लिए नाइट्रोजन के लिए एकाकी-युग्म प्रतिकर्षण का स्पष्ट उपचार आवश्यक है, जिसमें $\sigma$-एरोमैटिसिटी छोटे वलयों के लिए एक महत्वपूर्ण सुधार के रूप में कार्य करती है। यह अध्ययन अनुभवजन्य बल क्षेत्रों और क्वांटम यांत्रिक वास्तविकता के बीच के अंतर को पाटता है, जो CCSD(T) ऊर्जाओं से सीधे जुड़े और भौतिक रूप से ठोस डिस्क्रिप्टर प्रदान करता है।