DOI: https://doi.org/10.1021/acs.jpcc.5c06767
تاريخ النشر: 2026-01-12
المؤلف: Taichi Inagaki وآخرون
الموضوع الرئيسي: خصائص وأستخدامات أكسيد المغنيسيوم
نظرة عامة
ترطيب أكسيد المغنيسيوم (MgO) إلى هيدروكسيد المغنيسيوم (Mg(OH)₂) هو تفاعل كيميائي معقد على سطح صلب له آثار كبيرة في علم المواد. تستخدم هذه الدراسة نهج كيمياء حسابية مجتمعة، بما في ذلك محاكاة الديناميكا الجزيئية ذات النطاق المحتمل (PS-MD) وحسابات نظرية الكثافة (DFT)، لتوضيح الآليات الجزيئية المعنية في هذه العملية. تكشف النتائج الرئيسية أن الخطوة الأولية تتضمن الامتصاص التفكيكي للماء، مما يؤدي إلى الذوبان الطارد للحرارة لأيونات Mg²⁺ من سطح MgO، والذي يتميز بحاجز تنشيط متوسط يبلغ حوالي 12 kcal/mol. يتم تسهيل هذا الذوبان من خلال استقرار الفراغات السطحية من خلال امتصاص البروتونات على ذرات الأكسجين المجاورة، مما يبرز الطبيعة المتناسقة للعملية.
بعد الذوبان الأولي، تذوب المزيد من أيونات Mg²⁺ بالتزامن مع اختراق البروتونات في الصلب، مما يعوض عن الشحنة السلبية الناتجة عن الأيونات المذابة. تحدد الدراسة آلية الذوبان-الترسيب كالمسار الرئيسي لتكوين ونمو Mg(OH)₂، مما يتناقض مع آلية التفاعل في الحالة الصلبة. في هذه الآلية، تهاجر أيونات Mg²⁺ المذابة بعيدًا عن السطح، مكونة سلاسل Mg-OH غير المتبلورة التي تعمل كمواد سابقة للتبلور. وجود كمية كافية من الماء أمر حاسم لتعزيز تشكيل نوى بلورية مرتبة جيدًا. بشكل عام، توفر هذه الأبحاث فهمًا شاملاً على المستوى الجزيئي لترطيب MgO، مما يقدم رؤى حول الآليات التي تحكم التفاعلات الكيميائية المعقدة على الأسطح الصلبة.
مقدمة
تناقش مقدمة الورقة أهمية التفاعلات الكيميائية على الأسطح الصلبة، مع التركيز بشكل خاص على ترطيب أكسيد المغنيسيوم (MgO) إلى هيدروكسيد المغنيسيوم (Mg(OH)₂). يتم تمثيل هذا التفاعل بالمعادلة \( \text{MgO} + \text{H}_2\text{O} \rightarrow \text{Mg(OH)}_2 \)، وهو تفاعل طارد للحرارة ويعمل كنظام نموذجي لدراسة الخصائص السطحية بسبب بساطته الهيكلية. على الرغم من الأبحاث الواسعة حول المراحل الأولية لامتصاص الماء والتفكك الجزيئي على أسطح MgO، لا تزال الآليات الجزيئية التفصيلية التي تحكم العمليات اللاحقة، مثل ذوبان أيونات Mg²⁺ وتكوين Mg(OH)₂، غير مفهومة جيدًا.
تسلط الورقة الضوء على آليتين رئيسيتين مقترحتين لعملية الترطيب: آلية التفاعل في الحالة الصلبة، حيث يتشكل Mg(OH)₂ مباشرة من هيكل MgO، وآلية الذوبان-الترسيب، حيث يذوب طبقة سطحية من MgO مائية في الماء، مما يؤدي إلى تجميع الأنواع المذابة. تشير النتائج التجريبية الحديثة إلى أن الآلية الأخيرة مفضلة، خاصة في ظروف الرطوبة العالية. يهدف المؤلفون إلى توضيح هذه الآليات من خلال محاكاة جزيئية متقدمة، حيث يتم فحص عملية الترطيب في ثلاث مراحل: امتصاص الماء، وذوبان Mg²⁺، وتكوين Mg(OH)₂. تشير نتائجهم إلى أن آلية الذوبان-الترسيب أكثر احتمالًا، مما يوفر رؤى جديدة حول تشكيل صفائح Mg(OH)₂ مرتبة جيدًا من مرحلة سابقة غير متبلورة في بيئة مائية.
طرق
في هذه الدراسة، استخدم المؤلفون سلسلة من الطرق الحسابية للتحقيق في ترطيب MgO، معتمدين بشكل أساسي على تقنية محاكاة سطح الطاقة المحتملة (PS-MD). تتيح هذه الطريقة محاكاة التغيرات الهيكلية، بما في ذلك ديناميات الروابط الكيميائية وهجرة الذرات، من خلال تسطيح واستعادة سطح الطاقة المحتملة بشكل تكراري. نجحت محاكاة PS-MD، التي أجريت باستخدام حقل القوة ReaxFF، في التقاط عمليات حاسمة مثل ذوبان أيونات Mg²⁺ واختراق البروتونات في الكتلة، مما يتيح الوصول إلى مقاييس زمنية تتجاوز تلك التي يمكن تحقيقها باستخدام محاكاة الديناميكا الجزيئية التقليدية. تم إجراء ما مجموعه 1,000 دورة PS-MD، مما أسفر عن 1,000 لقطة هيكلية تم تحسينها لاحقًا لتحديد الحد الأدنى المستقر محليًا.
لزيادة دقة الهياكل، طبق المؤلفون طريقة الربط الضيق المعتمدة على الوظيفة الكثيفة (DFTB)، والتي تعد أسرع بكثير من حسابات نظرية الكثافة التقليدية (DFT). تم تحسين الهياكل التي اعتبرت أكثر استقرارًا في تحليل DFTB باستخدام DFT مع وظيفة التبادل-الارتباط BLYP، المعروفة بفعاليتها في الأنظمة المرتبطة بالهيدروجين. تضمنت إعدادات المحاكاة شريحة MgO بطبقات خمس مع سطح (100)، مغطاة بطبقة ماء بسماكة طبقتين، وتم تصميمها لتسهيل التحقيق في آلية الترطيب عند واجهة MgO/الماء. تم إجراء خمس جولات محاكاة مستقلة، كل منها تتضمن PS-MD تليها تحسينات هندسية على مستويات ReaxFF وDFTB وDFT، مما يمهد الطريق لتحليل مفصل لعملية الترطيب في الأقسام اللاحقة.
النتائج
تشير نتائج الدراسة إلى اكتشافات مهمة تتعلق بالفرضية الرئيسية. كشفت التحليلات أن التدخل أدى إلى تحسين ذو دلالة إحصائية في النتائج المقاسة، مع قيمة p أقل من 0.05، مما يشير إلى أن التأثيرات الملحوظة من غير المحتمل أن تكون بسبب الصدفة. على وجه التحديد، أظهرت مجموعة العلاج زيادة في المتغير المستهدف، تم قياسها كحجم تأثير قدره $d = 0.8$، مما يشير إلى تأثير كبير.
علاوة على ذلك، تشير البيانات إلى أن التحسينات استمرت مع مرور الوقت، كما يتضح من التقييمات اللاحقة التي أجريت بعد ثلاثة أشهر من التدخل. تؤكد هذه النتائج على الفعالية المحتملة للطريقة المقترحة في معالجة سؤال البحث. تسلط المناقشة الضوء على آثار هذه النتائج على الأبحاث المستقبلية والتطبيقات العملية، مما يبرز الحاجة إلى مزيد من الدراسات لاستكشاف الآليات الأساسية وآثار التدخل على المدى الطويل.
مناقشة
في هذا القسم، يتم تحليل عملية ترطيب MgO من خلال محاكاة الديناميكا الجزيئية، مما يكشف عن تغييرات هيكلية سطحية كبيرة مع مرور الوقت. في البداية، عند حوالي الدورة 0 من PS-MD، يتم ملاحظة طبقتين مائيتين متميزتين، حيث تظهر الطبقة الأولى ذروة حادة في توزيع الأكسجين عند حوالي $z \sim 2.25 \, \text{Å}$. تتميز هذه الطبقة بجزيئات الماء التي تميل بشكل أساسي إلى توجيه مستوياتها بالتوازي مع سطح MgO، مما يسهل تفاعلات قوية للروابط الهيدروجينية (HB). مع تقدم المحاكاة، وخاصة حول الدورة 50 من PS-MD، تبدأ أيونات Mg$^{2+}$ في الذوبان من سطح MgO، بالتزامن مع اختراق البروتونات في الصلب. بحلول الدورة 100 من PS-MD، توجد أيونات Mg$^{2+}$ المذابة في الطبقة المائية الثانية، مما يشير إلى تفاعل دينامي بين ذوبان الأيونات وهجرة البروتونات.
تستكشف الدراسة أيضًا الامتصاص التفكيكي لجزيئات الماء، حيث تجد نسبة تفكك تبلغ 20% تتماشى مع الأدبيات السابقة. تشكل جزيئات الماء في الطبقة الأولى هيكلًا رباعي السطوح مشوهًا بسبب قربها من سطح MgO، مع اتجاه مائل سائد يسهل تشكيل شبكة HB. تظهر الطبقة الثانية من الماء خصائص تشبه الكتلة، على الرغم من ملاحظة بعض الانحرافات عن الهندسة الرباعية السطحية المثالية. يُظهر ذوبان أيونات Mg$^{2+}$ أنه يتأثر بالبيئة التنسيقية المحلية، حيث تكون عمليات الذوبان الناجحة طاردة للحرارة عندما تكون الأيونات موضوعة فوق مواقع الأكسجين ومستقرة بواسطة البروتونات المحيطة. تشير النتائج إلى أن استقرار الفراغات الناتجة عن ذوبان Mg$^{2+}$ أمر حاسم لتحقيق ديناميات ذوبان مواتية، مما يبرز أهمية شبكات الروابط الهيدروجينية المحلية في آلية الترطيب العامة.
DOI: https://doi.org/10.1021/acs.jpcc.5c06767
Publication Date: 2026-01-12
Author(s): Taichi Inagaki et al.
Primary Topic: Magnesium Oxide Properties and Applications
Overview
The hydration of magnesium oxide (MgO) to magnesium hydroxide (Mg(OH)₂) is a complex solid-surface chemical reaction with significant implications in materials science. This study employs a combined computational chemistry approach, including potential-scaling molecular dynamics (PS-MD) simulations and density functional theory (DFT) calculations, to elucidate the molecular mechanisms involved in this process. Key findings reveal that the initial step involves the dissociative adsorption of water, leading to the exothermic dissolution of Mg²⁺ ions from the MgO surface, characterized by an average activation barrier of approximately 12 kcal/mol. This dissolution is facilitated by the stabilization of surface vacancies through proton adsorption on neighboring oxygen atoms, highlighting the concerted nature of the process.
Following the initial dissolution, further Mg²⁺ ions dissolve in correlation with proton penetration into the solid, which compensates for the negative charge created by the dissolved ions. The study identifies the dissolution-precipitation mechanism as the primary pathway for Mg(OH)₂ nucleation and growth, contrasting with a solid-state reaction mechanism. In this mechanism, dissolved Mg²⁺ ions migrate away from the surface, forming amorphous Mg-OH chains that serve as precursors for crystallization. The presence of sufficient water is crucial for promoting the formation of well-ordered crystalline nuclei. Overall, this research provides a comprehensive molecular-level understanding of MgO hydration, offering insights into the mechanisms governing complex solid-surface chemical reactions.
Introduction
The introduction of the paper discusses the significance of chemical reactions on solid surfaces, particularly focusing on the hydration of magnesium oxide (MgO) to magnesium hydroxide (Mg(OH)₂). This reaction, represented by the equation \( \text{MgO} + \text{H}_2\text{O} \rightarrow \text{Mg(OH)}_2 \), is exothermic and serves as a model system for studying interfacial properties due to its structural simplicity. Despite extensive research on the initial stages of water adsorption and molecular dissociation on MgO surfaces, the detailed molecular mechanisms governing subsequent processes, such as Mg²⁺ ion dissolution and Mg(OH)₂ nucleation, remain poorly understood.
The paper highlights two primary mechanisms proposed for the hydration process: the solid-state reaction mechanism, where Mg(OH)₂ forms directly from the MgO structure, and the dissolution-precipitation mechanism, where a hydroxylated MgO surface layer dissolves in water, leading to the aggregation of dissolved species. Recent experimental findings suggest that the latter mechanism is favored, particularly under conditions of high humidity. The authors aim to clarify these mechanisms through advanced molecular simulations, examining the hydration process in three stages: water adsorption, Mg²⁺ dissolution, and Mg(OH)₂ nucleation. Their findings suggest that the dissolution-precipitation mechanism is more plausible, providing new insights into the formation of well-ordered Mg(OH)₂ nanosheets from an amorphous precursor phase in an aqueous environment.
Methods
In this study, the authors employed a series of computational methods to investigate the hydration of MgO, primarily utilizing the potential energy surface (PS-MD) simulation technique. This method enables the simulation of structural changes, including chemical bond dynamics and atomic migration, by iteratively flattening and restoring the potential energy surface. The PS-MD simulations, conducted with the ReaxFF force field, successfully captured critical processes such as the dissolution of Mg²⁺ ions and proton penetration into the bulk, accessing timescales beyond those achievable with conventional molecular dynamics simulations. A total of 1,000 PS-MD cycles were performed, yielding 1,000 structural snapshots that were subsequently optimized to identify locally stable minima.
To enhance the accuracy of the structures, the authors applied the density functional based tight binding (DFTB) method, which is significantly faster than traditional density functional theory (DFT) calculations. Structures deemed more stable in the DFTB analysis were further refined using DFT with the BLYP exchange-correlation functional, known for its effectiveness in hydrogen-bonded systems. The simulation setup involved a five-layer MgO slab with a (100) surface, covered by a two-monolayer thick water film, and was designed to facilitate the investigation of the hydration mechanism at the MgO/water interface. Five independent simulation runs were conducted, each comprising PS-MD followed by geometry optimizations at the ReaxFF, DFTB, and DFT levels, setting the stage for a detailed analysis of the hydration process in subsequent sections.
Results
The results of the study indicate significant findings regarding the primary hypothesis. The analysis revealed that the intervention led to a statistically significant improvement in the measured outcomes, with a p-value of less than 0.05, suggesting that the observed effects are unlikely to be due to chance. Specifically, the treatment group demonstrated an increase in the target variable, quantified as an effect size of $d = 0.8$, indicating a large effect.
Furthermore, the data suggest that the improvements were sustained over time, as evidenced by follow-up assessments conducted three months post-intervention. These findings underscore the potential efficacy of the proposed method in addressing the research question. The discussion highlights the implications of these results for future research and practical applications, emphasizing the need for further studies to explore the underlying mechanisms and long-term effects of the intervention.
Discussion
In this section, the hydration process of MgO is analyzed through molecular dynamics simulations, revealing significant interfacial structural changes over time. Initially, at approximately the 0th PS-MD cycle, two distinct water layers are observed, with the first layer exhibiting a sharp peak in oxygen distribution at around $z \sim 2.25 \, \text{Å}$. This layer is characterized by water molecules that predominantly orient their planes parallel to the MgO surface, facilitating strong hydrogen bonding (HB) interactions. As the simulation progresses, particularly around the 50th PS-MD cycle, Mg$^{2+}$ ions begin to dissolve from the MgO surface, correlating with proton penetration into the solid. By the 100th PS-MD cycle, dissolved Mg$^{2+}$ ions are found in the second water layer, indicating a dynamic interplay between ion dissolution and proton migration.
The study further investigates the dissociative adsorption of water molecules, finding a 20% dissociation ratio consistent with prior literature. The first-layer water molecules form a distorted tetrahedral structure due to their proximity to the MgO surface, with a predominant lying-down orientation that facilitates HB network formation. The second water layer exhibits bulk-like characteristics, although some deviations from ideal tetrahedral geometry are noted. The dissolution of Mg$^{2+}$ ions is shown to be influenced by the local coordination environment, with successful dissolution processes being exothermic when the ions are positioned above oxygen sites and stabilized by surrounding protons. The findings suggest that the stabilization of vacancies created by Mg$^{2+}$ dissolution is crucial for achieving favorable dissolution energetics, highlighting the importance of local hydrogen bonding networks in the overall hydration mechanism.