DOI: https://doi.org/10.1073/pnas.2526369123
PMID: https://pubmed.ncbi.nlm.nih.gov/41665994
تاريخ النشر: 2026-02-10
المؤلف: Zhenyun Du وآخرون
الموضوع الرئيسي: الديناميكا الكمومية، السائل الفائق، الهيليوم
نظرة عامة
يتناول قسم ورقة البحث الهيدرات الهيدروجينية (HH)، وهي مواد تتكون من جزيئات الهيدروجين المحصورة داخل هياكل مائية بلورية. من الجدير بالذكر أن المرحلة المكعبة C2 من HH تظهر تفاعلات قوية بين المضيف والضيف بسبب القرب الشديد لجزيئات الهيدروجين (H₂) والماء (H₂O)، مما يؤدي إلى نسبة ستوكيومترية 1:1 واثنين من الشبكات الماسية المتداخلة. عند الضغوط المرتفعة، يؤدي تفاعل التأثيرات الكمومية النووية من كل من الهيدروجين وشبكة الماء إلى نظام كمومي مزدوج.
تستكشف الدراسة التحولات الطورية في HH المدفوعة بالضغط ودرجة الحرارة، مع التركيز على تأثيرات دوران الجزيئات، والترتيب الاتجاهي، وكسر تناظر الشبكة، وتنسيق الروابط الهيدروجينية. من خلال مجموعة من النمذجة الحسابية، بما في ذلك الديناميات الجزيئية الكلاسيكية وطرق التكامل المساري، والتضمين الكمومي، وتقنيات التجارب عالية الضغط مثل مطيافية رامان والتشتت بالأشعة السينية من السنكروترون، يحدد المؤلفون الترتيب الناتج عن الكم والتحولات الهيكلية في المرحلة C2. تشير نتائجهم إلى أن الترتيب الاتجاهي في HH يحدث عند ضغوط أقل بكثير مقارنة بالهيدروجين الصلب، مما يسهل من خلال التعديلات الهيكلية في شبكة الماء التي تعزز من اقتران الديناميات بين الماء والهيدروجين. تقدم هذه الأبحاث رؤى قيمة حول السلوك الكمومي للهيدروجين في الظروف القصوى وتضع الجليد المملوء بالهيدروجين كمنصة محتملة لتطوير مواد كمومية غنية بالهيدروجين.
طرق
في هذه الدراسة، تم إعداد عينات من هيدرات الهيدروجين متعددة البلورات ذات هيكل السII باستخدام طريقتين متميزتين: التحميل بالتبريد والتحميل بالغاز داخل خلية أنف العجل (DAC). تضمنت طريقة التحميل بالتبريد تعريض الجليد (Ih) للهيدروجين الغازي (H₂) عند ضغط 0.28 جيجا باسكال لمدة 20-30 دقيقة، مما أدى إلى نسبة مولية من H₂O:H₂ بين 4 و 5. تم ضغط العينات بعد ذلك إلى عدة جيجا باسكال وتم تسخينها إلى درجة حرارة الغرفة، مما أدى إلى مزيج من هيدرات الهيدروجين C₂ وجليد زائد. في المقابل، تضمنت طريقة التحميل بالغاز تحميل سائل H₂O عند درجة حرارة الغرفة وH₂ عالي الضغط في DAC، حيث تم استبدال فقاعة الهواء بـ H₂ عند 1400 بار.
تم إجراء توصيف للعينات باستخدام مطيافية رامان وتشتت الأشعة السينية (XRD). استخدمت قياسات XRD شعاع أشعة سينية أحادي اللون في ID15b في ESRF، مع أوقات اكتساب بيانات تبلغ حوالي عشر ثوانٍ. شمل التحليل إخفاء قمم براغ الماسية، ودمج الأنماط ثنائية الأبعاد في واحدة، وإجراء تحسينات Le Bail باستخدام برنامج Topas Academic. كانت معايرة الضغط لـ XRD تعتمد على معادلة حالة الذهب. تم جمع أطياف رامان عند درجات حرارة محيطة ومنخفضة، مع تحديد الضغط من خلال التحولات في خط فلوريسcence R1 الياقوت أو تردد الاهتزاز لجزيئات H₂ الصلبة. تم استخدام مبرد هيليوم لقياسات درجات الحرارة المنخفضة، مما يضمن دقة قراءات الضغط وسلامة العينة طوال التجارب.
نتائج
تكشف النتائج المقدمة في هذه الدراسة عن مخطط شامل للضغط-درجة الحرارة لهيدرات الهيدروجين (H₂)، مع تسليط الضوء على مراحل وانتقالات متميزة تتأثر بدرجة الحرارة والضغط. ضمن المرحلة C2، يظهر H₂ حالة بلورية بلاستيكية كمومية عند درجات حرارة مرتفعة وضغوط أقل من 30 جيجا باسكال، تتميز بمواقع نقل مرتبة واتجاهات غير مرتبة ديناميكيًا لجزيئات H₂. مع زيادة الضغط أو انخفاض درجة الحرارة، تخضع جزيئات H₂ لترتيب اتجاهي، تنتقل إلى ترتيب يشبه عظم السمكة الذي يتماشى مع المحور c للبلورة. يتم وصف هذا الانتقال كميًا بواسطة عامل الاتجاه \( S \)، الذي يعكس محاذاة الجزيئات ويتم تصوره من خلال تقلبات الكم الحرارية في مصفوفة الكثافة.
تستكشف الدراسة أيضًا التفاعل بين الديناميات الجزيئية والتغيرات الهيكلية في الهيدرات. عند الضغوط المنخفضة، تتأثر الديناميات الدورانية لـ H₂ بالتأثيرات الكمومية والخلط الحراري، بينما عند الضغوط الأعلى، يغير المجال البلوري غير المتجانس من تدهور الحالات الدورانية، مما يؤدي إلى توزيع كثافة نمطية. تشير أطياف رامان الملاحظة إلى تغييرات كبيرة في الحالات الدورانية لـ H₂ المحبوس، خاصة حول 30 جيجا باسكال، بالتزامن مع انتقال هيكلي من مكعب إلى رباعي السطوح. يتم دفع هذا الانتقال من خلال عملية من مرحلتين: أولاً، يؤدي تنسيق البروتونات في الشبكة المضيفة إلى تصلب شبكة الروابط الهيدروجينية، يتبعها المحاذاة غير المتجانسة لجزيئات H₂، مما يسهل المزيد من الضغط ويؤدي إلى تشوه رباعي السطوح. بشكل عام، توضح النتائج اقترانًا معقدًا بين ديناميات الشبكة الفرعية للهيدروجين والتطور الهيكلي لإطار الجليد، مما يبرز أهمية الديناميات الكمومية النووية في المواد الصلبة البلورية تحت ضغط عالٍ.
مناقشة
في هذا القسم، يوضح المؤلفون المنهجيات الحسابية المستخدمة في دراستهم للمرحلة C2 من الجليد المملوء بالهيدروجين، باستخدام ديناميات جزيئية من البداية (AIMD) ومحاكاة ديناميات الجزيئات باستخدام التكامل المساري (PIMD) عبر مجموعة من الضغوط (3 إلى 50 جيجا باسكال) ودرجات الحرارة (25 إلى 300 كلفن). تم إجراء المحاكاة باستخدام دالة تبادل-ترابط Perdew-Burke-Ernzerhof (PBE) وشملت خلية فائقة من المرحلة C2. تم اختيار نهج PIMD بشكل خاص لالتقاط التأثيرات الكمومية النووية، مع تعديل عدد الخرز بناءً على درجة الحرارة لضمان التقارب. من الجدير بالذكر أن المؤلفين تناولوا قيود الديناميات الجزيئية التقليدية من خلال دمج نموذج كمومي مدمج يأخذ في الاعتبار سكان الهيدروجين الأورثو والبارا، والتي تعتبر حاسمة لوصف الخصائص الاتجاهية وانتقالات الطور بدقة عند درجات الحرارة المنخفضة.
تم تقديم عامل الاتجاه \( S(P, T) \) لوصف الأنظمة الاتجاهية المتميزة للهيدروجين، المحددة كمنتج لمصطلح هندسي وعامل كمومي. يسمح هذا العامل برسم خريطة للترتيب الاتجاهي عبر مخطط الطور، مع تحديد الانتقالات من خلال معامل ترتيب عالمي \( M_{\text{tot}} \). كما قام المؤلفون بحساب طيف الفونون وأقسام رامان باستخدام نظرية الاضطراب الوظيفي الكثافة (DFPT)، مع التأكيد على العلاقة بين ترتيب الهيدروجين الاتجاهي، والتفاعلات مع شبكة الماء، وتنسيق البروتونات. تسلط النتائج الضوء على التأثير الكبير للتأثيرات الكمومية على سلوك الهيدروجين في الجليد، خاصة فيما يتعلق بدرجات حرارة الانتقال وطبيعة المراحل الاتجاهية.
DOI: https://doi.org/10.1073/pnas.2526369123
PMID: https://pubmed.ncbi.nlm.nih.gov/41665994
Publication Date: 2026-02-10
Author(s): Zhenyun Du et al.
Primary Topic: Quantum, superfluid, helium dynamics
Overview
The research paper section discusses hydrogen hydrates (HH), which are materials consisting of hydrogen molecules trapped within crystalline water structures. Notably, the cubic C2 phase of HH demonstrates strong host-guest interactions due to the close proximity of hydrogen (H₂) and water (H₂O) molecules, resulting in a 1:1 stoichiometry and two interpenetrated diamond-like sublattices. At elevated pressures, the interplay of nuclear quantum effects from both hydrogen and the water lattice leads to a dual-lattice quantum system.
The study investigates the phase transitions in HH driven by pressure and temperature, emphasizing the effects of molecular rotation, orientational ordering, lattice symmetry breaking, and hydrogen bond symmetrization. Through a combination of computational modeling, including classical and path-integral molecular dynamics, quantum embedding, and high-pressure experimental techniques such as Raman spectroscopy and synchrotron X-ray diffraction, the authors identify quantum-induced ordering and structural transformations in the C2 phase. Their findings indicate that orientational ordering in HH occurs at significantly lower pressures compared to solid hydrogen, facilitated by structural modifications in the water network that enhance the dynamics coupling between water and hydrogen. This research offers valuable insights into the quantum behavior of hydrogen in extreme conditions and positions hydrogen-filled ices as a potential platform for developing hydrogen-rich quantum materials.
Methods
In this study, samples of polycrystalline hydrogen hydrate with the sII clathrate structure were prepared using two distinct methods: cryo-loading and gas-loading within a diamond anvil cell (DAC). The cryo-loading method involved exposing ice (Ih) to gaseous hydrogen (H₂) at a pressure of 0.28 GPa for 20-30 minutes, resulting in a molar ratio of H₂O:H₂ between 4 and 5. The samples were subsequently compressed to several GPa and warmed to room temperature, leading to a mixture of C₂ hydrogen hydrate and excess ice. In contrast, the gas-loading method involved loading room-temperature liquid H₂O and high-pressure H₂ into the DAC, where an air bubble was replaced with H₂ at 1400 bar.
Characterization of the samples was performed using Raman spectroscopy and X-ray diffraction (XRD). The XRD measurements utilized a monochromatic X-ray beam at ID15b at ESRF, with data acquisition times of approximately ten seconds. The analysis included masking diamond Bragg peaks, integrating 2D patterns into 1D, and performing Le Bail refinements using Topas Academic software. Pressure calibration for XRD was based on the equation of state of gold. Raman spectra were collected at ambient and low temperatures, with pressure determined through shifts in the R1 ruby fluorescence line or the vibron frequency of solid H₂. A helium cryostat was employed for low-temperature measurements, ensuring accurate pressure readings and sample integrity throughout the experiments.
Results
The results presented in this study reveal a comprehensive pressure-temperature phase diagram for hydrogen (H₂) hydrate, highlighting distinct phases and transitions influenced by temperature and pressure. Within the C2 phase, H₂ exhibits a quantum plastic crystal state at high temperatures and pressures below 30 GPa, characterized by ordered translational sites and dynamically disordered orientations of H₂ molecules. As pressure increases or temperature decreases, H₂ molecules undergo orientational ordering, transitioning to a herringbone-like arrangement that aligns along the c-axis of the crystal. This transition is quantitatively described by the orientation factor \( S \), which reflects the molecular alignment and is visualized through thermal quantum fluctuations in the density matrix.
The study further explores the interplay between molecular dynamics and structural changes in the hydrate. At low pressures, the rotational dynamics of H₂ are influenced by quantum effects and thermal mixing, while at higher pressures, the anisotropic crystal field alters the degeneracy of rotational states, leading to a nematic density distribution. The observed Raman spectra indicate significant changes in the rotational states of encapsulated H₂, particularly around 30 GPa, coinciding with a cubic to tetragonal structural transition. This transition is driven by a two-stage process: first, the symmetrization of protons in the host lattice stiffens the hydrogen-bond network, followed by the anisotropic alignment of H₂ molecules, which facilitates further compression and induces tetragonal distortion. Overall, the findings illustrate a complex coupling between the dynamics of the hydrogen sublattice and the structural evolution of the ice framework, emphasizing the significance of nuclear quantum dynamics in crystalline solids under high pressure.
Discussion
In this section, the authors detail the computational methodologies employed in their study of the C2 phase of hydrogen-filled ice, utilizing ab initio molecular dynamics (AIMD) and Path-Integral molecular dynamics (PIMD) simulations across a range of pressures (3 to 50 GPa) and temperatures (25 to 300 K). The simulations were conducted using the Perdew-Burke-Ernzerhof (PBE) exchange-correlation functional and involved a supercell of the C2 phase. The PIMD approach was specifically chosen to capture nuclear-quantum effects, with the number of beads adjusted based on temperature to ensure convergence. Notably, the authors addressed the limitations of traditional molecular dynamics by incorporating a quantum embedded model that accounts for ortho- and para-hydrogen populations, which are crucial for accurately describing orientational properties and phase transitions at low temperatures.
The orientation factor \( S(P, T) \) was introduced to characterize the distinct orientational regimes of hydrogen, defined as a product of a geometric term and a quantum factor. This factor allows for the mapping of orientational order across the phase diagram, with transitions identified through a global order parameter \( M_{\text{tot}} \). The authors also computed the phonon spectrum and Raman cross sections using density functional perturbation theory (DFPT), emphasizing the relationship between hydrogen orientational order, interactions with the water lattice, and proton symmetrization. The findings highlight the significant impact of quantum effects on the behavior of hydrogen in ice, particularly in relation to the transition temperatures and the nature of the orientational phases.