تحسين ذرات هيرشفيلد الدورية Periodic Hirshfeld Atom Refinement

شارك:
المجلة: The Journal of Physical Chemistry Letters، المجلد: 17، العدد: 11
DOI: https://doi.org/10.1021/acs.jpclett.5c03918
PMID: https://pubmed.ncbi.nlm.nih.gov/41758498
تاريخ النشر: 2026-02-27
المؤلف: Kanghyun Chu وآخرون
الموضوع الرئيسي: الهياكل البلورية للمركبات الكيميائية

نظرة عامة

تناقش هذه القسم تطوير نوع جديد من تحسين ذرات هيرشفيلد (HAR)، الذي يُطلق عليه اسم HAR الدوري (pHAR)، والذي يوسع تطبيق HAR ليشمل أي هيكل شبكة دورية. تم تحقيق هذا التقدم من خلال استخدام مدارات غاوسيّة مركزية حول الذرات ضمن صيغة موجة بلوتش، مع الحفاظ على التوافق مع منهجيات HAR التقليدية.

تم تقييم فعالية pHAR باستخدام بيانات حيود الأشعة السينية عالية الجودة من البلورات الأحادية للبورانات والبورات، التي تحتوي على روابط N-H و B-H في بيئات كيميائية متنوعة. تشير النتائج إلى أن pHAR يوفر أطوال روابط X-H تتماشى بشكل وثيق مع بيانات المرجع المستمدة من تجارب حيود النيوترونات، مما يظهر دقة محسّنة. ومن الجدير بالذكر أن هذه الدراسة توسع بشكل كبير البيانات الهيكلية التجريبية المتاحة حول روابط B-H، حيث تضاعف تقريبًا مجموعة البيانات الحالية.

طرق

تحدد هذه القسم الطرق التجريبية والحسابية المستخدمة في الدراسة، معتمدة بشكل أساسي على حزم البرمجيات Crystal23 و Tonto لحسابات الكيمياء الكمومية وتحليل كثافة الإلكترون. تقوم Crystal23 بإجراء الحسابات تحت ظروف حدود دورية باستخدام مدارات ذرية من نوع غاوسي، بينما يركز Tonto على تقسيم كثافة الإلكترون وتحسين المعلمات الهيكلية. تم تعيين عتبة تقارب المجال الذاتي (SCF) إلى $10^{-7}$ هارتري، وتم أخذ عينات من منطقة بريلوان باستخدام طريقة باك-مونكهورست، مع تحديد نقاط أخذ عينات محددة لمختلف أنظمة البلورات. كان نصف قطر الخلية الفائقة حوالي 20 Å لالتقاط ذيول كثافة الإلكترون بشكل كافٍ، وتم تسهيل تبادل البيانات بين حزم البرمجيات من خلال مصفوفات الكثافة.

لزيادة الكفاءة الحسابية، يتخطى Tonto حساب حاصل دوال الأساس تحت ظروف معينة، وتم تأكيد صلاحية كثافة الإلكترون المعاد إنتاجها من خلال المقارنات مع خرائط الكثافة من Crystal23. كما توضح الدراسة عمليات تحسين ذرات الهيدروجين واستخدام مجموعة الأساس pob-TZVP-rev2 جنبًا إلى جنب مع الوظيفة الهجينة B3LYP. ومن الجدير بالذكر أن تنفيذ روتين فرعي لتوحيد عامل الشكل في Tonto يسمح بتحسين التعامل مع الأنظمة ذات التناظر العالي ويعزز الاستقرار العددي لمصفوفات المربعات الصغرى. يمكّن هذا التقدم المنهجي من إجراء مقارنات فعالة بين تحسين ذرات الهيدروجين التقليدي (HAR) ومنهج pHAR المقترح، مما يساهم في دقة وموثوقية التحليلات الهيكلية التي تم إجراؤها في الدراسة.

نقاش

تسلط قسم النقاش في ورقة البحث الضوء على التقدم في تقنيات حيود الأشعة السينية لتحديد الهياكل الذرية بدقة، مع التركيز بشكل خاص على قيود نموذج الذرة المستقلة (IAM) وفوائد الأساليب الحديثة مثل نموذج المضاعفات (MM) وتحسين ذرات هيرشفيلد (HAR). على الرغم من استخدام IAM على نطاق واسع، إلا أنه يفشل في حساب التباين الناتج عن الروابط، مما يؤثر بشكل خاص على موضع ذرات الهيدروجين. بالمقابل، تتضمن MM و HAR نماذج أكثر تعقيدًا لكثافة الإلكترون، مما يحسن دقة معلمات ذرات الهيدروجين وأطوال الروابط. يسمح الأسلوب الجديد لتحسين ذرات هيرشفيلد الدوري (pHAR) بمعالجة المركبات الشبكية الدورية، مما يعزز دقة قياسات أطوال الروابط، خاصة لروابط B-H و N-H، مقارنة ببيانات حيود النيوترونات.

تشير النتائج إلى أنه بينما يظهر pHAR تحسينات كبيرة على IAM، فإن التحسينات على HAR التقليدي هامشية. تؤكد الدراسة على أهمية اختيار مجموعة الأساس، خاصة استخدام الدوال المشتتة، في تحقيق توافق هيكلي مثالي مع بيانات النيوترونات. ومن الجدير بالذكر أن pHAR يقلل من الانحرافات المطلقة المتوسطة في أطوال الروابط، خاصة للروابط الأكثر قطبية، مما يشير إلى أن الوصف الدوري الصحيح للبيئة أمر حاسم لتحديد الهيكل بدقة. بشكل عام، يمثل pHAR خطوة كبيرة إلى الأمام في تحسين الهياكل الجزيئية، مما يمكّن من اكتشاف الميزات الكيميائية الدقيقة التي كانت صعبة الحل سابقًا.

Journal: The Journal of Physical Chemistry Letters, Volume: 17, Issue: 11
DOI: https://doi.org/10.1021/acs.jpclett.5c03918
PMID: https://pubmed.ncbi.nlm.nih.gov/41758498
Publication Date: 2026-02-27
Author(s): Kanghyun Chu et al.
Primary Topic: Crystal structures of chemical compounds

Overview

The section discusses the development of a new variant of Hirshfeld Atom Refinement (HAR), termed periodic HAR (pHAR), which extends the applicability of HAR beyond molecular crystals to any periodic-network structure. This advancement is achieved by employing atom-centered Gaussian orbitals within a Bloch wave formalism, maintaining compatibility with conventional HAR methodologies.

The efficacy of pHAR was evaluated using high-quality single-crystal X-ray diffraction data from boranes and borates, which contain N-H and B-H bonds in various chemical environments. The findings indicate that pHAR provides X-H bond lengths that closely align with reference data obtained from neutron-diffraction experiments, showcasing enhanced precision. Notably, this study significantly expands the available experimental structural data on B-H bonds, nearly doubling the existing dataset.

Methods

The section outlines the experimental and computational methods employed in the study, primarily utilizing the software packages Crystal23 and Tonto for quantum-chemical calculations and electron-density analysis. Crystal23 conducts calculations under periodic boundary conditions using Gaussian-type atomic orbitals, while Tonto focuses on electron-density partitioning and structural parameter refinement. The self-consistent field (SCF) convergence threshold was set to $10^{-7}$ Hartree, and the Brillouin zone was sampled using the Pack-Monkhorst method, with specific sampling points determined for various crystal systems. The supercell radius was approximately 20 Å to adequately capture electron-density tails, and data exchange between the two software packages was facilitated through density matrices.

To enhance computational efficiency, Tonto skips the calculation of the product of basis functions under certain conditions, and the validity of the reproduced electron density was confirmed through comparisons with density maps from Crystal23. The study also details the refinement processes for hydrogen atoms and the use of the pob-TZVP-rev2 basis set along with the B3LYP hybrid functional. Notably, the implementation of a subroutine for form-factor symmetrization in Tonto allows for improved handling of high-symmetry systems and enhances the numerical stability of least-squares matrices. This methodological advancement enables effective comparisons between conventional hydrogen atom refinement (HAR) and the proposed pHAR approach, ultimately contributing to the accuracy and reliability of the structural analyses performed in the study.

Discussion

The discussion section of the research paper highlights advancements in X-ray diffraction techniques for accurately determining atomic structures, particularly focusing on the limitations of the Independent Atom Model (IAM) and the benefits of modern approaches like the multipole model (MM) and Hirshfeld Atom Refinement (HAR). The IAM, while widely used, fails to account for bonding-induced anisotropy, particularly affecting hydrogen atom positioning. In contrast, MM and HAR incorporate more complex electron density models, improving the accuracy of hydrogen atom parameters and bond lengths. The newly developed periodic Hirshfeld Atom Refinement (pHAR) method allows for the treatment of periodic network compounds, enhancing the precision of bond length measurements, especially for B-H and N-H bonds, compared to neutron diffraction data.

The findings indicate that while pHAR shows significant improvements over IAM, the enhancements over conventional HAR are marginal. The study emphasizes the importance of basis set selection, particularly the use of diffuse functions, in achieving optimal structural agreement with neutron data. Notably, pHAR reduces average absolute deviations in bond lengths, particularly for more polar bonds, suggesting that the correct periodic description of the environment is crucial for accurate structural determination. Overall, pHAR represents a significant step forward in the refinement of molecular structures, enabling the detection of subtle chemical features that were previously challenging to resolve.