DOI: https://doi.org/10.1038/s41467-024-47879-0
PMID: https://pubmed.ncbi.nlm.nih.gov/38693110
تاريخ النشر: 2024-04-30
المؤلف: Joseph P. Heindel وآخرون
الموضوع الرئيسي: التطبيقات الطيفية والدراسات الكيميائية الكمومية
طرق
في هذه الدراسة، يستخدم المؤلفون نموذج الكتلة-الاستمرارية للتحقيق في تفاعلات الأيون-الماء من خلال محاكاة الديناميكا الجزيئية لقطرات الماء بقطر 4 نانومتر، باستخدام حقل القوة التفاعلي ReaxFF/CGeM. يقومون بإنشاء 100 تكوين لـ \(X(H_2O)_{35}\) من محاكاة قطرات محايدة ومشحونة تحتوي على أعداد متباينة من أيونات الهيدروكسيد (\(OH^-\)) وأيونات الهيدروجين (\(H^+\)). يتم تحديد الشروط الأولية للقطرات باستخدام PACKMOL، تليها توازن عند 293 كلفن لمدة 1 نانوثانية، مع 250 بيكوثانية إضافية لجني التكوين. للمقارنة، يتم أيضًا إنشاء 200 تكوين لقطرات ماء نقية. يتم التحقق من توتر السطح للنموذج، مما يعطي قيمًا تبلغ 60 مللي نيوتن/متر و64 مللي نيوتن/متر، والتي تتماشى بشكل معقول مع القيمة التجريبية البالغة 72 مللي نيوتن/متر.
يتم حساب مساهمة الطاقة المرتبطة بالارتباط في الإنثالبي عند مستوى نظرية ωB97M-V/aug-cc-pVDZ باستخدام Q-Chem، مع تعريف الطاقة المرتبطة بالارتباط بواسطة المعادلة:
\[
\Delta H(n)_{gas}[X] = \langle E[X(H_2O)_n] \rangle – \langle E[(H_2O)_n] \rangle – E[X]
\]
تلتقط هذه الطريقة المساهمات الاهتزازية، بما في ذلك عدم التناسق، وتأخذ في الاعتبار طاقة التجويف، مما يميزها عن نماذج الكتلة-الاستمرارية التقليدية التي تركز فقط على الهياكل الدنيا العالمية. يتم حساب طاقة الذوبان باستخدام مزيج من حسابات الاستمرارية القابلة للاستقطاب عبر نموذج COSMO ودمج الشحنات لمعالجة تأثير الشحنات الزائدة في البيئة، مع تمثيل تحول طاقة الذوبان كفرق في طاقات الذوبان لحالات الأيون الفردي والمتعدد.
نتائج
يقدم قسم “النتائج” في ورقة البحث النتائج الرئيسية المستمدة من التجارب أو التحليلات التي تم إجراؤها. يسلط الضوء على النتائج المهمة التي تدعم الفرضيات أو أسئلة البحث المطروحة في الدراسة. عادةً ما يتم توضيح النتائج من خلال أشكال مختلفة من تمثيل البيانات، مثل الجداول، الرسوم البيانية، أو التحليلات الإحصائية، التي توفر فهمًا بصريًا واضحًا للنتائج.
قد يتضمن القسم أيضًا مقارنات بين المجموعات التجريبية، تدابير التحكم، وأي ارتباطات أو اتجاهات تم ملاحظتها. غالبًا ما يتم الإبلاغ عن الدلالة الإحصائية، مما يشير إلى موثوقية النتائج. بشكل عام، يخدم هذا القسم لنقل الأدلة التجريبية التي تدعم استنتاجات الدراسة، مع التأكيد على آثار النتائج ضمن السياق الأوسع لمجال البحث.
مناقشة
في هذا القسم، يناقش المؤلفون الديناميكا الحرارية لنقل الإلكترون في القطرات المشحونة، خاصة في سياق تفاعلات الأكسدة والاختزال في البيئات المائية. يبرزون أن الأكسدة تتطلب عادةً الجذور الهيدروكسيلية (OH•)، بينما يعتمد الاختزال على ظروف نقل الإلكترون المواتية. تستخدم الدراسة نموذجًا ميكانيكيًا كمومي-إحصائي لتقدير مكونات الطاقة الحرة للترطيب لأيونات الهيدونيوم (H⁺) والهيدروكسيد (OH⁻) في القطرات المشحونة. تكشف محاكاة الديناميكا الجزيئية أن إنتروبيا الترطيب تتأثر بشكل أساسي بمساحة السطح القابلة للوصول للمذيب (SASA) للأيونات، مع كون إنثالبي الترطيب هو العامل السائد في الطاقة الحرة الديناميكية الحرارية.
تشير النتائج إلى أنه مع زيادة الشحنة في القطرات، تنخفض إنثالبيات الترطيب لكل من أيونات H⁺ وOH⁻ بشكل كبير، مما يشير إلى أن نقل الإلكترون يمكن أن يصبح مواتيًا ديناميكيًا حراريًا حتى عند تركيزات الشحنة أقل من حد رايلي. على وجه التحديد، يجد المؤلفون أن حوالي 12 أيون OH⁻ أو 16 أيون H⁺ يمكن أن تتجاوز حاجز 107 كيلو كالوري/مول لنقل الإلكترون التلقائي. كما يستكشفون استقرار الجذور OH• وH• في البيئات المشحونة، مشيرين إلى أن طاقة التأين العمودية (VIE) لـ OH⁻ تنخفض مع زيادة الشحنة، بينما تصبح الألفة الإلكترونية العمودية (VEA) لـ H⁺ أكثر ملاءمة. تمتد آثار هذه النتائج إلى القطرات الدقيقة الأكبر، حيث تنطبق مبادئ الديناميكا الحرارية المماثلة، مما يوفر فهمًا شاملاً لتفاعلات الأكسدة والاختزال المعجلة في أنظمة القطرات الدقيقة.
DOI: https://doi.org/10.1038/s41467-024-47879-0
PMID: https://pubmed.ncbi.nlm.nih.gov/38693110
Publication Date: 2024-04-30
Author(s): Joseph P. Heindel et al.
Primary Topic: Spectroscopy and Quantum Chemical Studies
Methods
In this study, the authors employ a cluster-continuum model to investigate ion-water interactions through molecular dynamics simulations of water droplets with a radius of 4 nm, utilizing the ReaxFF/CGeM reactive force field. They generate 100 configurations of \(X(H_2O)_{35}\) from simulations of neutral and charged droplets containing varying numbers of hydroxide (\(OH^-\)) and hydrogen ions (\(H^+\)). Initial conditions for the droplets are established using PACKMOL, followed by equilibration at 293 K for 1 ns, with an additional 250 ps for configuration harvesting. For comparison, 200 configurations of pure water droplets are also generated. The surface tension of the model is validated, yielding values of 60 mN/m and 64 mN/m, which align reasonably with the experimental value of 72 mN/m.
The binding energy contribution to the enthalpy is calculated at the ωB97M-V/aug-cc-pVDZ level of theory using Q-Chem, with the binding energy defined by the equation:
\[
\Delta H(n)_{gas}[X] = \langle E[X(H_2O)_n] \rangle – \langle E[(H_2O)_n] \rangle – E[X]
\]
This approach captures vibrational contributions, including anharmonicity, and accounts for cavitation energy, distinguishing it from traditional cluster-continuum models that focus solely on global minimum structures. The solvation energy is computed using a combination of polarizable continuum calculations via the COSMO model and charge embedding to address the influence of excess charges in the environment, with the solvation energy shift represented as the difference in solvation energies for single and multiple ion cases.
Results
The “Results” section of the research paper presents the key findings derived from the conducted experiments or analyses. It highlights the significant outcomes that support the hypotheses or research questions posed in the study. The results are typically illustrated through various forms of data representation, such as tables, graphs, or statistical analyses, which provide a clear visual understanding of the findings.
The section may also include comparisons between experimental groups, control measures, and any observed correlations or trends. Statistical significance is often reported, indicating the reliability of the results. Overall, this section serves to convey the empirical evidence that underpins the study’s conclusions, emphasizing the implications of the findings within the broader context of the research field.
Discussion
In this section, the authors discuss the thermodynamics of electron transfer in charged droplets, particularly in the context of redox reactions in aqueous environments. They highlight that oxidation typically requires hydroxyl radicals (OH•), while reduction is contingent on favorable electron transfer conditions. The research employs a quantum mechanical-statistical mechanical model to estimate the hydration-free energy components of hydronium (H⁺) and hydroxide (OH⁻) ions in charged droplets. Molecular dynamics simulations reveal that the hydration entropy is primarily influenced by the solvent accessible surface area (SASA) of the ions, with the hydration enthalpy being the dominant factor in the thermodynamic free energy.
The findings indicate that as the charge in the droplets increases, the hydration enthalpies for both H⁺ and OH⁻ ions decrease significantly, suggesting that electron transfer can become thermodynamically favorable even at charge concentrations below the Rayleigh limit. Specifically, the authors find that approximately 12 OH⁻ or 16 H⁺ ions can overcome the 107 kcal/mol barrier for spontaneous electron transfer. They also explore the stabilization of OH• and H• radicals in charged environments, noting that the vertical ionization energy (VIE) of OH⁻ decreases with increasing charge, while the vertical electron affinity (VEA) of H⁺ becomes more favorable. The implications of these results extend to larger microdroplets, where similar thermodynamic principles apply, thereby providing a comprehensive understanding of accelerated redox reactions in microdroplet systems.