DOI: https://doi.org/10.1103/physreva.111.032818
تاريخ النشر: 2025-03-26
المؤلف: K. Nordlund وآخرون
الموضوع الرئيسي: الجيولوجيا الفيزيائية عالية الضغط والمواد
نظرة عامة
في هذه الدراسة، يحقق المؤلفون في الإمكانيات بين الذرات الطاردة التي تعتبر حاسمة لفهم قوة التوقف النووي وتشتت الذرات خلال أحداث الإشعاع keV و MeV. يستخدمون نموذج Ziegler-Biersack-Littmark (ZBL) ضمن نظرية الكثافة الوظيفية بدون مدارات (OF-DFT) ويقارنون توقعاته مع طرق الكيمياء الكمومية الأكثر تقدمًا، وتحديدًا نظرية الاضطراب من الدرجة الثانية لمولر-بليست (MP2) ونظرية الكثافة الوظيفية (DFT) باستخدام المدارات الذرية العددية (NAOs). تشير النتائج إلى أن الإمكانيات الطاردة المستمدة من حسابات MP2 و DFT تتماشى بشكل وثيق، حيث تختلف بنسبة تقارب 1% للطاقة فوق 30 eV، بينما يظهر نموذج ZBL اختلافات أكبر، مع اختلافات تبلغ حوالي 5-10% عند الطاقات الأعلى.
كما يقدم المؤلفون ملاءمات جديدة محددة للزوج لوظائف التصفية، والتي تمثل كمجموع من ثلاثة أسية، والتي تظهر مستوى عالٍ من الاتفاق (ضمن ∼ 2%) مع بيانات الكيمياء الكمومية. يتم استخدام هذه الإمكانيات المنقحة لمحاكاة ملفات عمق زراعة الأيونات في السيليكون البلوري الأحادي، مما ينتج عنه نتائج تتطابق بشكل وثيق مع الملاحظات التجريبية. ومع ذلك، تشير الدراسة إلى أنه تحت ظروف القناة، يمكن أن تؤثر المكونات الجذابة للإمكانات بشكل كبير على ملفات العمق، مما يشير إلى أن النماذج الطاردة البحتة قد لا تكون كافية لتوقعات دقيقة لزراعة الأيونات. تتوفر مجموعات البيانات الشاملة والملاءمات التحليلية كموارد مفتوحة الوصول.
مقدمة
تسلط مقدمة هذه الورقة البحثية الضوء على أهمية تأثيرات الإشعاع عبر مجالات مختلفة، بما في ذلك التكنولوجيا النووية وعلوم المواد، خاصة في سياق زراعة الأيونات لتصنيع الدوائر المتكاملة. تؤكد على أن الضرر الناتج عن الإشعاع يتأثر بشكل أساسي بالسلاسل التصادمية التي تبدأها الجسيمات النشيطة، والتي تؤدي إلى تعديلات في المادة بسبب طاقة الارتداد المنقولة إلى النوى داخل العينة. يتم التحكم في ديناميات التفاعل بواسطة طاقة الإمكانات الطاردة بين النوى، والتي تعتمد بشكل أساسي على المسافات بين الذرات بدلاً من عوامل أخرى مثل زوايا الروابط أو أرقام التنسيق.
يحدد المؤلفون نهجهم للتحقيق في الإمكانيات بين الذرات من خلال حسابات ab initio على ثلاثة مستويات نظرية: نظرية الاضطراب من الدرجة الثانية لمولر-بليست للعناصر الأخف، ونظرية الكثافة الوظيفية مع التأثيرات النسبية للعناصر الأثقل، ونظرية الكثافة الوظيفية بدون مدارات (OF-DFT) لنموذج إمكانيات عالمي. يخططون لمقارنة منهجية للنتائج من هذه الطرق، خاصة ضد الإمكانية المستخدمة على نطاق واسع Ziegler-Biersack-Littmark (ZBL)، وتحليل الآثار المترتبة على نتائجهم على الكميات القابلة للقياس تجريبيًا. تهدف الورقة إلى تقديم فهم شامل للإمكانات الطاردة بين الذرات وأهميتها لتأثيرات الإشعاع في المواد.
طرق
في هذا القسم، يوضح المؤلفون المنهجيات المستخدمة في بحثهم، مع التركيز بشكل خاص على دمج مختلف الإمكانيات بين الذرات لمحاكاة زراعة الأيونات في السيليكون. يقدمون الوصول إلى مجموعات بيانات ZBL المحددة للزوج، DMol، و MP2، جنبًا إلى جنب مع معاملات إمكانيات NLH، من خلال أرشيف مفتوح الوصول. يتم تقديم تحليل مقارن للمدى المتوسط لمختلف الأيونات وطاقات الزراعة، مع تسليط الضوء على دقة النتائج المحاكية مقابل البيانات التجريبية، مع انحرافات متوسطة تتراوح من 0.23% إلى 4.9%.
يؤكد المؤلفون على أهمية دمج التفاعلات متعددة الجسيمات القريبة من التوازن مع الإمكانيات الطاردة عالية الطاقة دون إدخال معلمات حرة، لتقليل الشكوك في نتائج القناة. يقترحون طريقتين للتداخل لضمان انتقال سلس بين الإمكانيات الطاردة وإمكانات التوازن، مع الحفاظ على الاستمرارية في الطاقة ومشتقاتها. تتناول المناقشة أيضًا تطبيق الإمكانيات المدربة بواسطة الآلة، والتي يتم تدريبها بشكل منهجي ضد مجموعة شاملة من بيانات DFT، مما يسمح بدمج متسق للإمكانات قصيرة المدى وإمكانات التوازن. تهدف هذه الطريقة إلى تعزيز دقة المحاكاة المتعلقة بتأثيرات الإشعاع في المواد، خاصة في سياق كل من تقريب التصادم الثنائي (BCA) ومحاكاة الديناميكا الجزيئية (MD).
نتائج
يقدم قسم النتائج تحليلًا مقارنًا لمحاكاة زراعة الأيونات لـ As-Si و B-Si و P-Si، كما هو موضح في الأشكال 18 و 19 و 20، على التوالي. تقارن الدراسة هذه النتائج مع تلك التي تم الحصول عليها من كود SRIM، الذي لا يأخذ في الاعتبار تأثيرات القناة ويظهر اختلافات كبيرة مع البيانات التجريبية. يتم تقديم ملخص عددي للنتائج في الجدول VI، الذي يستبعد مدى SRIM المتوسط بسبب عدم توافقه مع النتائج التجريبية. تظهر نتائج محاكاة MDRANGE ارتباطًا قويًا مع البيانات التجريبية، يُعزى إلى استخدام نموذج كثافة إلكترونية ثلاثي الأبعاد لقوة التوقف الإلكترونية.
تظهر الدراسة أيضًا أن إمكانيات MP2 و DMol و NLH تعطي نتائج مشابهة، تتماشى مع أسسها النظرية، بينما تتماشى أيضًا الإمكانيات المحددة للزوج ZBL بشكل وثيق مع النتائج التجريبية، حيث تظهر اختلافات متوسطة تبلغ فقط 2.2% و 2.8% و 3.3% لـ B-Si و P-Si و As-Si، على التوالي، فوق 100 eV. بالمقابل، تظهر الإمكانية العالمية ZBL انحرافًا أكبر، خاصة في ذيل القناة لـ B-Si و P-Si. على الرغم من عدم وجود إمكانيات واحدة تقدم توافقًا مثاليًا مع التجارب، فإن الأساليب من المبادئ الأولى (MP2 و DMol) تعطي عمومًا أقرب النتائج. تشير التحليلات إلى أن الانحرافات قد تنشأ من قيود نموذج التوقف الإلكتروني وتبرز حساسية ملفات النطاق لتباين الشعاع، والذي لم يتم الإبلاغ عنه سابقًا في الدراسات التجريبية. بشكل عام، تؤكد النتائج على ضرورة دمج بنية البلورة في محاكاة ملفات نطاق زراعة الأيونات لتعزيز الدقة.
مناقشة
في هذا القسم، يناقش المؤلفون الإطار النظري لحساب الإمكانيات الطاردة بين الذرات، مع التركيز بشكل خاص على صيغة وظيفة التصفية. تعتبر الإمكانية الكولومبية، \( V_{\text{Coul}}(r) = \frac{1}{4\pi \varepsilon_0} \frac{Z_1 Z_2 e^2}{r} \)، أساسًا لفهم التفاعلات بين النوى العارية، بينما يسمح تقريب بورن-أوبنهايمر بفصل الحركات النووية والإلكترونية. يتم اشتقاق الإمكانية بين الذرات \( V(r) \) من الطاقات الكلية لنظام ثنائي الذرة وذراته المكونة، مما يؤدي إلى التعبير \( V(r) = V_{\text{Coul}}(r) + E_{\text{el}}^{1+2}(r) – E_{\text{el}}^1 – E_{\text{el}}^2 \). يؤكد المؤلفون على أهمية وظيفة التصفية \( \phi(r) = \frac{V(r)}{V_{\text{Coul}}(r)} \)، التي تقترب من 1 عندما \( r \to 0 \) و 0 عندما \( r \to \infty \)، مما يسهل المحاكاة العددية الأكثر سلاسة عبر المسافات النووية المتغيرة.
تقدم المناقشة أيضًا نهج الذرة الموحدة، الذي يفترض أنه مع اقتراب ذرتين من بعضهما البعض، تشبه كثافتهما الإلكترونية المجمعة تلك الخاصة بذرة موحدة برقم ذري \( Z_1 + Z_2 \). يصحح هذا النهج قيود نموذج Ziegler-Biersack-Littmark (ZBL)، الذي يفترض تراكب الكثافات الذرية ويفشل في حساب التعديلات الكبيرة في دالة الموجة الإلكترونية عند المسافات الصغيرة. يبرز المؤلفون التحديات في نمذجة الهياكل الإلكترونية بدقة عند المسافات المتوسطة، حيث تحدث استرخاءات مدارية، ويلاحظون أن تغييرات كبيرة في الطاقة يمكن أن تنشأ من هذه التفاعلات. يختتم القسم بالتأكيد على الحاجة إلى منهجيات متقدمة لالتقاط تعقيدات الإمكانيات الطاردة، خاصة في الأنظمة التي تشمل المعادن الانتقالية أو اللانثانيدات، حيث قد لا توفر الكيمياء الكمومية التقليدية رؤى كافية.
DOI: https://doi.org/10.1103/physreva.111.032818
Publication Date: 2025-03-26
Author(s): K. Nordlund et al.
Primary Topic: High-pressure geophysics and materials
Overview
In this study, the authors investigate the repulsive interatomic potentials crucial for understanding nuclear stopping power and atom scattering during keV and MeV radiation events. They employ the Ziegler-Biersack-Littmark (ZBL) model within orbital-free density functional theory (OF-DFT) and compare its predictions with more advanced quantum chemical methods, specifically second-order Møller-Plesset perturbation theory (MP2) and density functional theory (DFT) using numerical atomic orbitals (NAOs). The results indicate that the repulsive potentials derived from MP2 and DFT calculations align closely, differing by approximately 1% for energies above 30 eV, while the ZBL model shows greater discrepancies, with differences of about 5-10% at higher energies.
The authors also introduce new pair-specific fits for the screening functions, represented as a sum of three exponentials, which demonstrate a high level of agreement (within ∼ 2%) with the quantum chemical data. These refined potentials are utilized to simulate ion implantation depth profiles in single crystalline silicon, yielding results that closely match experimental observations. However, the study notes that under channeling conditions, the attractive components of the potential can significantly influence the depth profiles, suggesting that purely repulsive models may not suffice for accurate ion implantation predictions. The comprehensive datasets and analytical fits are made available as open access resources.
Introduction
The introduction of this research paper highlights the significance of radiation effects across various domains, including nuclear technology and material science, particularly in the context of ion implantation for integrated circuit manufacturing. It emphasizes that the damage caused by radiation is primarily influenced by the collision cascades initiated by energetic particles, which lead to modifications in the material due to the recoil energy transferred to nuclei within the sample. The interaction dynamics are governed by the repulsive potential energy between nuclei, which is predominantly dependent on interatomic separations rather than other factors like bond angles or coordination numbers.
The authors outline their approach to investigating interatomic potentials through ab initio calculations at three theoretical levels: second-order Møller-Plesset perturbation theory for lighter elements, density-functional theory with relativistic effects for heavier elements, and orbital-free density functional theory (OF-DFT) for a universal potential model. They plan to systematically compare the results from these methods, particularly against the widely used Ziegler-Biersack-Littmark (ZBL) potential, and analyze the implications of their findings on experimentally measurable quantities. The paper aims to provide a comprehensive understanding of repulsive interatomic potentials and their relevance to radiation effects in materials.
Methods
In this section, the authors detail the methodologies employed in their research, particularly focusing on the integration of various interatomic potentials for simulating ion implantation in silicon. They provide access to the ZBL pair-specific, DMol, and MP2 datasets, along with NLH potential coefficients, through an open-access archive. A comparative analysis of mean ranges for different ions and implantation energies is presented, highlighting the accuracy of simulated results against experimental data, with average deviations ranging from 0.23% to 4.9%.
The authors emphasize the importance of systematically combining near-equilibrium many-body interactions with high-energy repulsive potentials without introducing free parameters, to minimize uncertainties in channeling results. They propose two interpolation methods to ensure a smooth transition between the repulsive and equilibrium potentials, maintaining continuity in energy and its derivatives. The discussion also touches on the application of machine-learned potentials, which are systematically trained against a comprehensive set of DFT data, allowing for a consistent integration of short-range and equilibrium potentials. This approach aims to enhance the accuracy of simulations related to radiation effects in materials, particularly in the context of both binary collision approximation (BCA) and molecular dynamics (MD) simulations.
Results
The results section presents a comparative analysis of ion implantation simulations for As-Si, B-Si, and P-Si, as illustrated in figures 18, 19, and 20, respectively. The study contrasts these results with those obtained from the SRIM code, which does not account for channeling effects and shows significant discrepancies with experimental data. A numerical summary of the findings is provided in table VI, which excludes SRIM mean ranges due to their lack of agreement with experimental results. The MDRANGE simulation results demonstrate a strong correlation with experimental data, attributed to the use of a 3D electron density model for electronic stopping power.
The study further reveals that the MP2, DMol, and NLH potentials yield similar results, consistent with their theoretical underpinnings, while the ZBL pair-specific potentials also align closely with experimental outcomes, exhibiting average differences of only 2.2%, 2.8%, and 3.3% for B-Si, P-Si, and As-Si, respectively, above 100 eV. In contrast, the ZBL universal potential shows greater deviation, particularly in the channeling tail for B-Si and P-Si. Although no single potential consistently provides the best agreement with experiments, the first-principles approaches (MP2 and DMol) generally yield the closest results. The analysis suggests that deviations may stem from limitations in the electronic stopping model and highlights the sensitivity of range profiles to beam divergence, which was not previously reported in experimental studies. Overall, the findings underscore the necessity of incorporating crystal structure in simulations of ion implantation range profiles to enhance accuracy.
Discussion
In this section, the authors discuss the theoretical framework for calculating repulsive interatomic potentials, particularly focusing on the screening function formalism. The Coulomb potential, \( V_{\text{Coul}}(r) = \frac{1}{4\pi \varepsilon_0} \frac{Z_1 Z_2 e^2}{r} \), serves as the foundation for understanding interactions between bare nuclei, while the Born-Oppenheimer approximation allows for the separation of nuclear and electronic motions. The interatomic potential \( V(r) \) is derived from the total energies of a diatomic system and its constituent atoms, leading to the expression \( V(r) = V_{\text{Coul}}(r) + E_{\text{el}}^{1+2}(r) – E_{\text{el}}^1 – E_{\text{el}}^2 \). The authors emphasize the significance of the screening function \( \phi(r) = \frac{V(r)}{V_{\text{Coul}}(r)} \), which approaches 1 as \( r \to 0 \) and 0 as \( r \to \infty \), facilitating smoother numerical simulations across varying internuclear distances.
The discussion also introduces the united atom approach, which posits that as two atoms approach each other, their combined electron density resembles that of a united atom with atomic number \( Z_1 + Z_2 \). This approach corrects the limitations of the Ziegler-Biersack-Littmark (ZBL) model, which assumes superposition of atomic densities and fails to account for significant electronic wave function modifications at small separations. The authors highlight the challenges in accurately modeling electronic structures at intermediate distances, where orbital relaxation occurs, and note that substantial energy changes can arise from these interactions. The section concludes by underscoring the need for advanced methodologies to capture the complexities of repulsive potentials, particularly in systems involving transition metals or lanthanides, where conventional quantum chemistry may not provide sufficient insights.