🔬 condensed matter

The Uhlenbeck-Ford model in two dimensions: Reference system for fluid-phase free-energy calculations

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

मूल लेखक: Samuel Cajahuaringa, Rodolfo Paula Leite, Maurice de Koning

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

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

CC BY 4.0

मूल लेखक: Samuel Cajahuaringa, Rodolfo Paula Leite, Maurice de Koning

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

कल्पना कीजिए कि आप एक नन्हे खोजकर्ता हैं जो एक द्वि-आयामी (two-dimensional) दुनिया में रहते हैं, जैसे कि ग्राफ पेपर की एक शीट पर रहने वाला कोई पात्र। हमारी वास्तविक, त्रि-आयामी (three-dimensional) दुनिया में, जब आप बर्फ के एक ठोस टुकड़े को गर्म करते हैं, तो वह एक ही नाटकीय छलांग में पानी में पिघल जाता है। लेकिन इस सपाट दुनिया में, चीजें अजीब होती हैं। कभी-कभी, बर्फ केवल पिघलती नहीं है; यह एक अजीब, डगमगाती हुई मध्यवर्ती अवस्था में बदल सकती है जिसे "हेक्साटिक" (hexatic) चरण कहा जाता है, जहाँ कण एक कठोर आकार बनाए रखने में सक्षम नहीं होते हैं, लेकिन फिर भी वे अपने पड़ोसियों के साथ एक विशिष्ट क्रम बनाए रखते हैं। वैज्ञानिक इन सपाट दुनियाओं को समझने के लिए जुनूनी हैं क्योंकि वे आधुनिक तकनीक में हर जगह मौजूद हैं, जैसे कि आपके फोन में ग्रेफीन की अति-पतली परतें और जिस तरह से परमाणु छोटी बूंदों में व्यवस्थित होते हैं। इन दुनियाओं के बारे में सटीक रूप से समझने के लिए, वैज्ञानिकों को "मुक्त ऊर्जा" (free energy) नामक कुछ गणना करने की आवश्यकता होती है। मुक्त ऊर्जा को एक स्कोरकार्ड की तरह समझें जो यह बताता है कि कणों का समूह किस अवस्था को पसंद करता है: क्या वे एक कठोर क्रिस्टल चाहते हैं, एक अस्त-व्यस्त तरल, या वह अजीब हेक्साटिक मध्यवर्ती अवस्था? जितना कम स्कोर होगा, कण उतने ही अधिक खुश होंगे। मुक्त ऊर्जा की गणना करना अविश्वसनीय रूप से कठिन है, जैसे कि ज्वार आने के दौरान समुद्र तट पर रेत के प्रत्येक कण को गिनने की कोशिश करना। गणित को आसान बनाने के लिए, वैज्ञानिक एक "संदर्भ प्रणाली" (reference system) का उपयोग करते हैं—एक सरल, काल्पनिक दुनिया जहाँ वे पहले से ही स्कोर जानते हैं, और फिर वे यह मापते हैं कि उनकी वास्तविक, जटिल दुनिया को हल करना उनके सरल वाले की तुलना में कितना कठिन है।

यह शोध पत्र इस बारे में है कि उन द्वि-आयामी तरल पदार्थों के लिए एक संदर्भ के रूप में उपयोग करने हेतु एकदम सही "सरल दुनिया" कैसे खोजी जाए। लेखकों ने "उहलेनबेक-फोर्ड" (Uhlenbeck-Ford - UF) मॉडल नामक एक विशिष्ट गणितीय मॉडल की जांच की। कल्पना कीजिए कि यह मॉडल कणों की एक ऐसी भीड़ है जो पास आने पर एक-दूसरे को बहुत धीरे से धकेलते हैं, लेकिन यदि वे बिल्कुल एक ही स्थान घेरने की कोशिश करते हैं, तो वह धक्का अनंत रूप से शक्तिशाली हो जाता है (हालांकि तुरंत नहीं, बल्कि लॉगरिदमिक रूप से)। इस मॉडल की सुंदरता यह है कि इसका व्यवहार मुख्य रूप से इस बात पर निर्भर करता है कि कण कितने घने हैं, जबकि तापमान धक्के की शक्ति को ऊपर या नीचे करने के लिए एक साधारण डायल की तरह कार्य करता है। शोधकर्ता जानना चाहते थे: क्या यह UF मॉडल 2D तरल पदार्थों के लिए एक अच्छा "संदर्भ सिस्टम" है? दूसरे शब्दों में, क्या हम वास्तविक, जटिल 2-आयामी तरल पदार्थों की मुक्त ऊर्जा की गणना करने के लिए इसे एक विश्वसनीय आधार के रूप में उपयोग कर सकते हैं?

इसका उत्तर देने के लिए, टीम ने दो मुख्य कार्य किए। सबसे पहले, उन्होंने "विरियल गुणांकों" (virial coefficients) की गणना करने के लिए भारी-भरकम गणित का उपयोग किया (जो कि इस तरह के सूत्र के शुरुआती कुछ पद हैं जिससे गैस के दबाव के घनत्व के साथ परिवर्तन का पता चलता है) दसवें क्रम तक। उन्होंने पाया कि हालांकि गणित बहुत कम घनत्व के लिए पूरी तरह से काम करता है, लेकिन जैसे-जैसे कण घने होते जाते हैं, यह बहुत जल्दी अव्यवधर और गलत हो जाता है। इसलिए, उन्होंने अपनी रणनीति बदली और घनत्व की एक विस्तृत श्रृंखला में दबाव और ऊर्जा के व्यवहार को मैप करने के लिए शक्तिशाली कंप्यूटर सिमुलेशन (मॉलिक्यूलर डायनेमिक्स) का उपयोग किया। उन्होंने एक अत्यधिक सटीक "इक्वेशन ऑफ स्टेट" (equation of state) बनाया जो UF मॉडल के व्यवहार को पूरी तरह से वर्णित करता है।

इसके बाद, उन्होंने इस मॉडल के "फेज डायग्राम" (phase diagram) की खोज की ताकि यह देखा जा सके कि घनत्व और "तापमान डायल" (जिसे वे स्केलिंग पैरामीटर $p$ कहते हैं) को बदलने पर क्या होता है। वे उन सीमाओं की तलाश कर रहे थे जहाँ तरल पदार्थ एक ठोस या उस रहस्यमय हेक्साटिक चरण में बदल जाता है। उन्होंने एक महत्वपूर्ण विवरण की खोज की: जबकि कम घनत्व पर स्केलिंग पैरामीटर $p \approx 70$ तक तरल अवस्था ही एकमात्र ऊष्मागतिक रूप से स्थिर अवस्था है, यह कम घनत्व के लिए सभी स्थितियों में सत्य नहीं है। उच्च घनत्व पर, यह मॉडल एक ठोस चरण और एक मध्यवर्ती हेक्साटिक चरण प्रदर्शित करता है। हालांकि, मुख्य निष्कर्ष यह है कि जब तक स्केलिंग पैरामीटर $p$ लगभग 70 से नीचे रहता है, तरल पदार्थ उन ठोस/हेक्साटिक सीमाओं तक पहुँचने से पहले पूरे घनत्व रेंज में स्थिर अवस्था बना रहता है। इसका अर्थ है कि UF मॉडल उस विशिष्ट पैरामीटर रेंज के भीतर तरल गणनाओं के लिए एक अत्यंत सुदृढ़ संदर्भ प्रणाली है। उन्होंने यह भी पुष्टि की कि उच्च घनत्व के लिए, यह मॉडल 2D पिघलने की समृद्ध, जटिल भौतिकी को पकड़ता है, जिसमें KTHNY सिद्धांत द्वारा अनुमानित दो-चरणीय पिघलने की प्रक्रिया शामिल है, जहाँ कण अपना कठोर ग्रिड खो देते हैं लेकिन अपनी दिशात्मक व्यवस्था बनाए रखते हैं (हेक्साटिक चरण) और अंततः एक अराजक तरल में पिघल जाते हैं।

अंत में, लेखकों ने अपने नए संदर्भ सिस्टम का परीक्षण किया। उन्होंने एक अलग, अधिक यथार्थवादी तरल जैसे कि "लेनार्ड-जोन्स" (Lennard-Jones) तरल (जो आकर्षित और प्रतिकर्षित करने वाले परमाणुओं के लिए एक मानक मॉडल है) की मुक्त ऊर्जा की गणना करने के लिए UF मॉडल का उपयोग किया। उन्होंने "नॉन-इक्विलिब्रियम थर्मोडायनामिक इंटीग्रेशन" नामक एक तकनीक का उपयोग किया, जो धीरे-धीरे सरल UF तरल को जटिल लेनार्ड-जोन्स तरल में बदलने जैसा है, जबकि इसके लिए आवश्यक कार्य को मापा जाता है। परिणाम आश्चर्यजनक थे: उनके द्वारा गणना की गई मुक्त ऊर्जा अन्य अत्यधिक सटीक विधियों के साथ लगभग पूरी तरह मेल खाती थी, जिसमें कुछ मामलों में अंतर 0.04% जितना छोटा था।

संक्षेप में, यह शोध पत्र सिद्ध करता है कि 2D उहलेनबेक-फोर्ड मॉडल द्वि-आयामी तरल पदार्थों की मुक्त ऊर्जा की गणना करने के लिए एक शानदार, विश्वसनीय उपकरण है, बशर्ते कि स्केलिंग पैरामीटर $p$ 70 के बराबर या उससे कम हो। यह वैज्ञानिकों को गणितीय उलझनों में फंसे बिना, सपाट सामग्रियों के जटिल पिघलने के व्यवहार और चरण संक्रमणों का अध्ययन करने के लिए एक ठोस, सटीक आधार प्रदान करता है। लेखकों ने केवल यह सुझाव नहीं दिया कि यह काम कर सकता है; उन्होंने इसका सिमुलेशन किया, इसके फेज डायग्राम को मैप किया (यह पहचानते हुए कि तरल कब एकमात्र स्थिर अवस्था नहीं रह जाता), और अन्य विधियों के विरुद्ध इसका सत्यापन किया, जिससे यह सिद्ध हुआ कि यह 2D सामग्री विज्ञान के भविष्य के लिए एक व्यावहारिक और सटीक संदर्भ प्रणाली है।

तकनीकी सारांश: मुक्त-ऊर्जा गणनाओं के लिए एक संदर्भ प्रणाली के रूप में 2D में उहलेनबेक-फोर्ड (Uhlenbeck–Ford) मॉडल

समस्या विवरण
त्रि-आयामी (3D) प्रणालियों की तुलना में द्वि-आयामी (2D) प्रणालियों का सैद्धांतिक अध्ययन अद्वितीय चुनौतियाँ प्रस्तुत करता है, विशेष रूप से चरण व्यवहार (phase behavior) के संबंध में जिसमें दोष-मध्यस्थ पिघलना (defect-mediated melting) और मध्यवर्ती हेक्साटिक (hexatic) चरण शामिल हैं। इन प्रणालियों के लक्षण वर्णन में एक केंद्रीय कठिनाई सटीक मुक्त ऊर्जा (free energy) की गणना करना है, जो प्रतिस्पर्धी चरणों की सापेक्ष स्थिरता निर्धारित करने के लिए आवश्यक है। यद्यपि थर्मोडायनामिक इंटीग्रेशन (TI) विधि ज्ञात मुक्त ऊर्जा वाले एक उपयुक्त संदर्भ तंत्र की उपलब्धता पर अत्यधिक निर्भर करती है, यह मुक्त ऊर्जा अंतरों के लिए एक सटीक मार्ग प्रदान करती है। हालांकि उहलेनबेक-फोर्ड (UF) मॉडल को 3D तरल पदार्थों के लिए एक कुशल संदर्भ प्रणाली के रूप में स्थापित किया गया है, लेकिन इसकी ऊष्मागतिक (thermodynamic) गुणधर्मों, चरण व्यवहार और 2D तरल पदार्थों के लिए एक संदर्भ प्रणाली के रूप में इसकी उपयुक्तता का व्यवस्थित रूप से लक्षण वर्णन नहीं किया गया है।

कार्यप्रणाली
लेखक 2D UF मॉडल की जांच करने के लिए एक बहु-आयामी कम्प्यूटेशनल दृष्टिकोण का उपयोग करते हैं:

सटीक विरियल विस्तार (Exact Virial Expansion): अध्ययन 2D UF मॉडल के लिए दसवें क्रम तक सटीक विरियल गुणांकों की गणना करता है। इसमें दिए गए क्रम $n$ के लिए सभी अनलेबल बाइकोनेक्टेड ग्राफ उत्पन्न करना और मेयर क्लस्टर एक्सपेंशन फॉर्मलिज्म का उपयोग करके उनके भार (weights) की गणना करना शामिल है।
आणविक गतिकी (Molecular Dynamics - MD) सिमुलेशन: लैंजवेन थर्मोस्टेट (Langevin thermostat) का उपयोग करते हुए LAMMPS कोड का उपयोग करके व्यापक MD सिमुलेशन किए गए हैं। ये सिमुलेशन समीकरण अवस्था (equation of state - EOS) के लिए डेटा उत्पन्न करने और चरण संक्रमणों का पता लगाने के लिए घनत्व और स्केलिंग मापदंडों ( $p$ ) की एक विस्तृत श्रृंखला को कवर करते हैं।
समीकरण अवस्था (Equation of State) निर्माण: उच्च घनत्व पर विरियल विस्तार की सीमाओं को दूर करने के लिए, लेखक MD दबाव डेटा के क्यूबिक स्प्लाइन इंटरपोलेशन (cubic spline interpolation) का उपयोग करके EOS के संख्यात्मक निरूपण तैयार करते हैं। इसके बाद सीधे ऊष्मागतिक संबंध के एकीकरण के माध्यम से अतिरिक्त हेल्महोल्ट्ज़ मुक्त ऊर्जा (excess Helmholtz free energy) प्राप्त की जाती है।
गैर-संतुलन ऊष्मागतिक एकीकरण (Non-Equilibrium Thermodynamic Integration - NETI): मुक्त-ऊर्जा गणनाओं को मान्य करने और मॉडल की उपयोगिता प्रदर्शित करने के लिए, लेखक एक NETI दृष्टिकोण का उपयोग करते हैं। इसमें सिस्टम ऑफ इंटरेस्ट और UF संदर्भ प्रणाली के बीच एक गैर-संतुलन स्विचिंग प्रक्रिया शामिल है। परिमित स्विचिंग समय से उत्पन्न होने वाली व्यवस्थित त्रुटियों को समाप्त करने के लिए फॉरवर्ड और बैकवर्ड स्विचिंग प्रक्रियाओं के औसत द्वारा रिवर्सिबल वर्क (reversible work) का अनुमान लगाया जाता है।
चरण आरेख मानचित्रण (Phase Diagram Mapping): चरण व्यवहार का विश्लेषण संरचनात्मक ऑर्डर मापदंडों की निगरानी करके किया जाता है, जिसमें स्टैटिक स्ट्रक्चर फैक्टर $S(k)$ , ट्रांसलेशनल ऑर्डर मापदंड $\Psi_T$ , और बॉन्ड-ओरिएंटेशनल ऑर्डर मापदंड $\Psi_6$ , तथा बॉन्ड-ओरिएंटेशनल सहसंबंध फलन $g_6(r)$ शामिल हैं। इन मेट्रिक्स का उपयोग KTHNY सिद्धांत के अनुसार ठोस, हेक्साटिक और तरल चरणों के बीच अंतर करने के लिए किया जाता है।

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

विरियल गुणांक और EOS: लेखक $p=1$ के लिए $n=10$ तक सटीक संख्यात्मक विरियल गुणांक प्रदान करते हैं। वे प्रदर्शित करते हैं कि जबकि कम घनत्व ( $x \lesssim 1$ ) पर विरियल विस्तार सटीक है, उच्च घनत्व पर गुणांकों के वैकल्पिक संकेतों और धीमी गिरावट के कारण यह तेजी से विचलन (diverge) करता है। फलस्वरूप, वे MD डेटा के स्प्लाइन इंटरपोलेशन का उपयोग करके पूर्णांक मानों (1, 25, 50, और 70) के लिए EOS और अतिरिक्त हेल्महोल्ट्ज़ मुक्त ऊर्जा के अत्यधिक सटीक संख्यात्मक निरूपण स्थापित करते हैं।
मुक्त ऊर्जा का सत्यापन: स्प्लाइन-आधारित EOS से प्राप्त अतिरिक्त हेल्महोल्ट्ज़ मुक्त ऊर्जा को NETI गणनाओं के विरुद्ध स्वतंत्र रूप से मान्य किया गया है। परिणाम पूरे घनत्व रेंज में $2 \times 10^{-4}$ से कम के सापेक्ष अवशेषों (relative residuals) के साथ उत्कृष्ट सहमति दिखाते हैं, जो निर्मित EOS की विश्वसनीयता की पुष्टि करते हैं।
चरण आरेख और स्थिरता सीमाएँ: अध्ययन 2D UF मॉडल के चरण आरेख का मानचित्रण करता है।
- कम घनत्व ( $x < 0.1$ ) के लिए, सिस्टम स्केलिंग मापदंड की परवाह किए बिना तरल चरण में रहता है।
- $x \geq 0.225$ के लिए, मॉडल विशिष्ट पिघलने के व्यवहार को प्रदर्शित करता है। निचले घनत्व ( $x=0.5$ ) पर, पिघलना ठोस से तरल के सीधे प्रथम-क्रम संक्रमण (first-order transition) के माध्यम से होता है। उच्च घनत्व ( $x=1.0$ ) पर, सिस्टम KTHNY सिद्धांत के अनुरूप दो-चरणीय पिघलने की स्थिति प्रदर्शित करता है, जिसमें एक मध्यवर्ती हेक्साटिक चरण होता है जहाँ ट्रांसलेशनल ऑर्डर लुप्त हो जाता है जबकि बॉन्ड-ओरिएंटेशनल ऑर्डर बना रहता है।
- महत्वपूर्ण रूप से, पूरे जांचे गए घनत्व रेंज में स्केलिंग मापदंड $p \lesssim 70$ तक तरल चरण को एकमात्र ऊष्मागतिक रूप से स्थिर चरण पाया गया है। इस सीमा के ऊपर, ठोस और हेक्साटिक चरण उभरते हैं।
लैनार्ड-जोन्स (Lennard–Jones - LJ) तरल पर अनुप्रयोग: 2D लैनार्ड-जोन्स तरल के हेल्महोल्ट्ज़ मुक्त ऊर्जा की गणना करके 2D UF मॉडल की व्यावहारिक उपयोगिता प्रदर्शित की गई है। UF मॉडल को संदर्भ के रूप में उपयोग करते हुए NETI गणनाएं, संशोधित बेनेडिक्ट-वेब-रूबिन (MBWR) समीकरण अवस्था और टू-फेज थर्मोडायनामिक (2PT) विधि से प्राप्त मूल्यों के साथ उत्कृष्ट सहमति (विचलन $< 1\%$ ) दर्शाती हैं।

महत्व
यह शोध पत्र 2D उहलेनबेक-फोल्ड मॉडल को 2D तरल पदार्थों के लिए मुक्त-ऊर्जा गणनाओं हेतु एक व्यावहारिक, सटीक और कम्प्यूटेशनल रूप से कुशल संदर्भ प्रणाली के रूप में स्थापित करता है। तरल चरण की ऊष्मागतिक स्थिरता सीमाओं (विशेष रूप से $p \lesssim 70$ के लिए) को परिभाषित करके, यह कार्य उस अनुप्रयोग की सीमा की पहचान करता जहाँ मॉडल का विश्वसनीय रूप से थर्मोडायमिक इंटीग्रेशन में उपयोग किया जा सकता है। 2D लैनार्ड-जोन्स तरल पर सफल अनुप्रयोग इस मॉडल की जटिल 2D प्रणालियों के लिए एक मजबूत संदर्भ के रूप में कार्य करने की क्षमता को प्रमाणित करता है, जो निम्न-आयामी सामग्रियों के लिए सैद्धांतिक विधियों के अंतराल को भरता है। लेखक भविष्य के अध्ययनों में इस मॉडल के उपयोग को सुगम बनाने के लिए आवश्यक संख्यात्मक निरूपण (EOS और मुक्त ऊर्जा अभिव्यक्ति) और स्क्रिप्ट प्रदान करते हैं।

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