DOI: https://doi.org/10.5194/ar-4-1-2026
تاريخ النشر: 2026-01-07
المؤلف: Ivo Neefjes وآخرون
الموضوع الرئيسي: كيمياء الغلاف الجوي والهباء الجوي
نظرة عامة
في هذا البحث، قام المؤلفون بتقييم منهجي لمجموعة متنوعة من الطرق الكيميائية الكمومية لتعزيز النمذجة الحاسوبية لمجموعات الجزيئات المائية في الغلاف الجوي، وخاصة تلك التي تتضمن الأحماض والقواعد ذات الصلة بالغلاف الجوي مع ما يصل إلى خمس جزيئات ماء. وجدوا أن الطريقة المستخدمة على نطاق واسع ωB97X-D/6-31++G(d,p) المدمجة مع تصحيح الطاقة الإلكترونية DLPNO NormalPNO-CCSD(T0)/aug-cc-pVTZ فعالة في التنبؤ بدقة بالخصائص الهيكلية والطاقة لهذه المجموعات. بالإضافة إلى ذلك، تعتبر طرق الكثافة الوظيفية المركبة مثل B97-3c و r2SCAN-3c و ωB97X-3c مناسبة لفحص المجموعات الأكبر، بينما يُوصى باستخدام نهج المدار الطبيعي المحلي LNO-CCSD(T)/aug-cc-pVTZ للتعديلات الدقيقة نظرًا لمتطلبات الذاكرة المنخفضة والدقة العالية.
استكشفت الدراسة أيضًا تصحيحات حرارية كيميائية تتجاوز تقريب المذبذب التوافقي التقليدي، حيث تم تقييم توزيعات الترطيب ومعدلات تشكيل الجسيمات لمجموعات مثل حمض الكبريتيك-الأمونيا (SA-AM) وحمض الكبريتيك-ثنائي ميثيل الأمين (SA-DMA) وحمض الميثان سلفونيك-ميثيل أمين (MSA-MA). أشارت النتائج إلى أن تأثيرات الترطيب تعتمد على النظام، حيث أظهرت مجموعات MSA-MA حساسية أكبر للترطيب مقارنةً بـ SA-AM أو SA-DMA، على الرغم من انخفاض معدلات تشكيل الجسيمات لديها. تشير النتائج إلى أنه بينما يكون للترطيب تأثير متواضع على تشكيل الجسيمات الجديدة تحت الظروف الجوية النموذجية، فإنه يلعب دورًا أكثر أهمية في مراحل النمو اللاحقة. يؤكد المؤلفون على الحاجة إلى طرق حسابية متسقة عند مقارنة توزيعات الترطيب ويوصون بمزيد من التحقيق في المجموعات الأكبر وتأثير الرطوبة على تطور الهباء الجوي.
مقدمة
تسلط مقدمة ورقة البحث الضوء على التأثير الكبير لجسيمات الهباء الجوي على المناخ العالمي وصحة الإنسان، مع التأكيد على دورها المزدوج في التأثيرات المناخية المباشرة وغير المباشرة. يمكن أن تنشأ الهباء الجوي من الانبعاثات المباشرة أو من خلال تشكيل الجسيمات الجديدة (NPF)، وهي عملية تتجمع فيها جزيئات الغاز ذات الضغط المنخفض لتشكيل مجموعات تنمو لتصبح جسيمات هباء جوي. تحدد الهيئة الحكومية الدولية المعنية بتغير المناخ (IPCC) الهباء الجوي كمصدر رئيسي للغموض في نماذج المناخ، خاصة بسبب الفهم المحدود للمراحل المبكرة من NPF، حيث تتشكل مجموعات يبلغ قطرها حوالي 2 نانومتر. تشمل المساهمات الرئيسية في NPF حمض الكبريتيك ومجموعة متنوعة من الأحماض العضوية، حيث تلعب المياه أيضًا دورًا حاسمًا، على الرغم من أن تأثيرها يبقى معقدًا ويعتمد على الظروف.
تناقش الورقة التحديات في دراسة NPF، خاصة فيما يتعلق بتوصيف الجسيمات التي يقل قطرها عن 2-3 نانومتر، حيث تكافح التقنيات التجريبية التقليدية لتحليلها بسبب محدودية رؤى التركيب الكيميائي. ظهرت الكيمياء الحاسوبية كأداة قيمة لاستكشاف هذه المجموعات، ومع ذلك، غالبًا ما يتم استبعاد المياه من النماذج لتقليل التكاليف الحاسوبية. تشير الدراسات الحديثة، مع ذلك، إلى أن المياه تؤثر بشكل كبير على نقل البروتون في مجموعات الأحماض والقواعد ويمكن أن تعزز أو تثبط NPF اعتمادًا على الظروف البيئية. يهدف هذا البحث إلى تقييم طرق الكيمياء الكمومية لنمذجة المجموعات المائية بدقة والتحقيق في كيفية تأثير إدراج المياه على فهم دور الرطوبة في NPF، مع التركيز على المراحل الأولية لتشكيل الجسيمات.
طرق
في هذه الدراسة، تم إجراء تقييم شامل لمجموعة متنوعة من طرق الكيمياء الكمومية (QC)، تتراوح من الطرق شبه التجريبية إلى الطرق المعتمدة على دالة الموجة عالية الدقة. شملت الطرق التي تم تقييمها r2SCAN-3c و ωB97X-3c، المصممة لتعزيز الدقة في التفاعلات بين الجزيئات مع الحفاظ على الكفاءة الحاسوبية. بالإضافة إلى ذلك، تم تقييم الوظائف الهجينة والوظائف المعممة التقريبية (GGA) مثل ωB97X-D و M06-2X، حيث تم الإشارة إلى ωB97X-D لأدائها الموثوق في الكيمياء الحرارية وحسابات طاقة الربط لمجموعات الجزيئات الجافة. تم تضمين الطريقة المعتمدة على دالة الموجة RI-MP2 وطرق المدار المحلية (DLPNO و LNO) أيضًا لتوفير طيف واسع من الدقة والكفاءة.
استخدمت الدراسة مجموعات أساسية متنوعة مصممة خصيصًا لطرق QC المحددة، مع استخدام مجموعات أساس Pople لـ M06-2X و PW91، ومجموعات أساس متسقة مع الارتباط المعزز لـ RI-MP2 و DLPNO-CCSD(T0) و LNO-CCSD(T). تم تقييم أداء ωB97X-D عبر جميع مجموعات الأساس، بما في ذلك متغير متخصص لـ LNO-CCSD(T) الذي يستبعد الوظائف المنتشرة على ذرات الهيدروجين لتحسين الكفاءة الحاسوبية. تم تطبيق التقنية المرتبطة صراحةً (F12) على طرق CCSD(T*)-F12 و DLPNO-CCSD(T0)-F12. استكشفت الدراسة أيضًا معايير تقليم مختلفة للمدارات الطبيعية في DLPNO-CCSD(T0) و LNO-CCSD(T)، مما يؤثر على كل من الدقة والتكلفة الحاسوبية. تم إجراء حسابات QC باستخدام عدة حزم برمجية، مع ملاحظة تقليص ملحوظ في وقت الحساب في أحدث إصدار من ORCA، على الرغم من الفروق الضئيلة في حسابات طاقة الربط.
النتائج
يقدم قسم “النتائج” في ورقة البحث النتائج الرئيسية المستمدة من التجارب أو التحليلات التي تم إجراؤها. يسلط الضوء على الاتجاهات البيانية المهمة، والنتائج الإحصائية، وأي علاقات ملحوظة بين المتغيرات. عادةً ما يتم توضيح النتائج من خلال الجداول أو الرسوم البيانية أو الأشكال، مما يوفر تمثيلًا بصريًا ويعزز فهم البيانات.
في هذا القسم، قد يناقش المؤلفون أيضًا تداعيات نتائجهم بالنسبة للأدبيات الموجودة، مؤكدين كيف تساهم نتائجهم في المجال الأوسع للدراسة. بالإضافة إلى ذلك، يتم الاعتراف بأي قيود واجهت خلال عملية البحث، مما يضمن تفسيرًا شاملاً للنتائج. بشكل عام، يخدم هذا القسم لتأكيد فرضية البحث ويضع الأساس للمناقشات والاستنتاجات اللاحقة.
المناقشة
في هذا القسم، يناقش المؤلفون تقييم طرق الكيمياء الكمومية (QC) للتنبؤ بالطاقة الربط الإلكترونية، والهندسات التوازنية، والخصائص الحرارية الكيميائية للمجموعات المائية الدقيقة المكونة من الأحماض والقواعد ذات الصلة بالغلاف الجوي. تم إنشاء ما مجموعه 395 مجموعة فريدة، وتم تقييم طرق QC المختلفة مقارنةً بطريقة مرجعية، CCSD(T*)-F12/cc-pVTZ-F12، التي تعتبر المعيار الذهبي للطاقة. تسلط الدراسة الضوء على أداء الطرق المختلفة، وخاصة طريقة DLPNO NormalPNO-CCSD(T0)/aug-cc-pVTZ، التي أظهرت خطأ مطلقًا متوسطًا منخفضًا (MAE) يبلغ حوالي 0.20 كيلو كالوري مول⁻¹ وتم اختيارها للمجموعات الأكبر بسبب كفاءتها ودقتها.
كما فحص المؤلفون الهندسات التوازنية للمجموعات، مستخدمين انحرافات الجذر التربيعي (RMSDs) لتقييم دقة طرق QC المختلفة مقارنةً بمصدر موثوق. بالإضافة إلى ذلك، استكشفوا الخصائص الحرارية الكيميائية، مؤكدين أهمية الترددات الاهتزازية في حساب طاقات جيبس الحرة للتكوين. وجدت الدراسة أن التصحيحات للانحرافات غير التوافقية والترددات الاهتزازية المنخفضة ضرورية للتنبؤات الحرارية الكيميائية الدقيقة، حيث تم استخدام التقريب شبه التوافقي (QHA) ونظرية الاضطراب الاهتزازي (VPT2) لتعزيز دقة الترددات الاهتزازية. بشكل عام، تؤكد النتائج على أهمية اختيار طرق QC المناسبة لدراسة الديناميكا الحرارية لمجموعات الجزيئات المائية، خاصة في سياقات الكيمياء الجوية.
DOI: https://doi.org/10.5194/ar-4-1-2026
Publication Date: 2026-01-07
Author(s): Ivo Neefjes et al.
Primary Topic: Atmospheric chemistry and aerosols
Overview
In this research, the authors systematically evaluated various quantum-chemical methods to enhance the computational modeling of hydrated atmospheric molecular clusters, particularly those involving atmospherically relevant acids and bases with up to five water molecules. They found that the widely used method ωB97X-D/6-31++G(d,p) combined with DLPNO NormalPNO-CCSD(T0)/aug-cc-pVTZ electronic energy correction is effective for accurately predicting structural and energetic properties of these clusters. Additionally, composite density functional methods like B97-3c, r2SCAN-3c, and ωB97X-3c are suitable for pre-screening larger clusters, while the local natural orbital approach LNO-CCSD(T)/aug-cc-pVTZ is recommended for precise refinements due to its low memory requirements and high accuracy.
The study also explored thermochemical corrections beyond the conventional harmonic oscillator approximation, assessing hydration distributions and particle formation rates for clusters such as sulfuric acid-ammonia (SA-AM), sulfuric acid-dimethylamine (SA-DMA), and methanesulfonic acid-methylamine (MSA-MA). Results indicated that hydration effects are system-dependent, with MSA-MA clusters showing more sensitivity to hydration than SA-AM or SA-DMA, despite their lower particle formation rates. The findings suggest that while hydration has a modest impact on new particle formation under typical atmospheric conditions, it plays a more significant role in later growth stages. The authors emphasize the need for consistent computational methods when comparing hydration distributions and recommend further investigation into larger clusters and the influence of humidity on aerosol evolution.
Introduction
The introduction of the research paper highlights the significant impact of aerosol particles on global climate and human health, emphasizing their dual role in direct and indirect climate effects. Aerosols can originate from direct emissions or through new particle formation (NPF), a process where low-volatility gas-phase molecules aggregate to form clusters that grow into aerosol particles. The Intergovernmental Panel on Climate Change (IPCC) identifies aerosols as a major source of uncertainty in climate models, particularly due to limited understanding of the early stages of NPF, where clusters approximately 2 nm in diameter form. Key contributors to NPF include sulfuric acid and various organic acids, with water also playing a critical role, although its influence remains complex and condition-dependent.
The paper discusses the challenges in studying NPF, particularly regarding the characterization of sub-2-3 nm particles, which traditional experimental techniques struggle to analyze due to limited chemical composition insights. Computational chemistry has emerged as a valuable tool to explore these clusters, yet water is often excluded from models to reduce computational costs. Recent studies, however, indicate that water significantly influences proton transfer in acid-base clusters and can either enhance or inhibit NPF depending on environmental conditions. This research aims to benchmark quantum chemistry methods for accurately modeling hydrated clusters and to investigate how the inclusion of water affects the understanding of humidity’s role in NPF, focusing on the initial stages of particle formation.
Methods
In this study, a comprehensive benchmarking of various quantum chemistry (QC) methods was conducted, ranging from semiempirical to high-accuracy wavefunction-based approaches. The methods evaluated included r2SCAN-3c and ωB97X-3c, which are designed to enhance accuracy in intermolecular interactions while maintaining computational efficiency. Additionally, hybrid and meta-generalized gradient approximation (GGA) functionals such as ωB97X-D and M06-2X were assessed, with ωB97X-D noted for its reliable performance in thermochemistry and binding energy calculations for dry molecular clusters. The wavefunction-based method RI-MP2 and localized orbital methods (DLPNO and LNO) were also included to provide a broad spectrum of accuracy and efficiency.
The study utilized various basis sets tailored to the specific QC methods, employing Pople basis sets for M06-2X and PW91, and augmented correlation-consistent basis sets for RI-MP2, DLPNO-CCSD(T0), and LNO-CCSD(T). The performance of ωB97X-D was evaluated across all basis sets, including a specialized variant for LNO-CCSD(T) that omits diffuse functions on hydrogen atoms to improve computational efficiency. The explicitly correlated (F12) technique was applied to CCSD(T*)-F12 and DLPNO-CCSD(T0)-F12 methods. The study also explored different truncation criteria for natural orbitals in DLPNO-CCSD(T0) and LNO-CCSD(T), impacting both accuracy and computational cost. QC calculations were performed using multiple software packages, with a notable reduction in computation time observed in the latest version of ORCA, despite negligible differences in binding energy calculations.
Results
The “Results” section of the research paper presents key findings derived from the conducted experiments or analyses. It highlights significant data trends, statistical outcomes, and any observed relationships between variables. The results are typically illustrated through tables, graphs, or figures, which provide visual representation and enhance the understanding of the data.
In this section, the authors may also discuss the implications of their findings in relation to existing literature, emphasizing how their results contribute to the broader field of study. Additionally, any limitations encountered during the research process are acknowledged, ensuring a comprehensive interpretation of the results. Overall, this section serves to validate the research hypothesis and lays the groundwork for subsequent discussions and conclusions.
Discussion
In this section, the authors discuss the benchmarking of quantum chemistry (QC) methods for predicting electronic binding energies, equilibrium geometries, and thermochemical properties of microhydrated clusters formed from atmospherically relevant acids and bases. A total of 395 unique clusters were generated, and various QC methods were evaluated against a reference method, CCSD(T*)-F12/cc-pVTZ-F12, which is considered the gold standard for energetics. The study highlights the performance of different methods, particularly the DLPNO NormalPNO-CCSD(T0)/aug-cc-pVTZ method, which demonstrated a low mean absolute error (MAE) of approximately 0.20 kcal mol⁻¹ and was selected for larger clusters due to its efficiency and accuracy.
The authors also examined the equilibrium geometries of the clusters, utilizing root-mean-square deviations (RMSDs) to assess the accuracy of various QC methods compared to a reliable reference. Additionally, they explored thermochemical properties, emphasizing the importance of vibrational frequencies in calculating Gibbs free energies of formation. The study found that corrections for anharmonicity and low vibrational frequencies are crucial for accurate thermochemical predictions, with the quasi-harmonic approximation (QHA) and vibrational perturbation theory (VPT2) being employed to enhance the accuracy of vibrational frequencies. Overall, the findings underscore the significance of selecting appropriate QC methods for studying the thermodynamics of hydrated molecular clusters, particularly in atmospheric chemistry contexts.