محاسبه چگالی اتمی و کمیت‌های مورد نیاز در MCNP

قلب تپنده محاسبات در کد MCNP، بر پایه “احتمال برخورد” استوار است. این احتمال مستقیماً به تعداد هسته‌های موجود در واحد حجم بستگی دارد. بنابراین، محاسبه دقیق چگالی اتمی (Atomic Density) و درک تفاوت‌های آن با چگالی جرمی، اولین قدم برای اطمینان از صحت نتایج شبیه‌سازی است. هرگونه خطا در این محاسبات اولیه، منجر به انحراف در مسیر آزاد میانگین ذرات و در نتیجه، بروز خطای سیستماتیک در تمامی تالی‌ها خواهد شد.

ضرورت محاسبه چگالی اتمی در شبیه‌سازی‌های مونت‌کارلو

در فیزیک هسته‌ای، ذرات با هسته‌های اتمی برهم‌کنش می‌کنند، نه با جرم کلی ماده. کد MCNP برای تعیین اینکه یک ذره پس از طی چه مسافتی با یک هسته برخورد می‌کند، از مفهومی به نام “چگالی اتمی” استفاده می‌کند که بیانگر تعداد اتم‌ها در هر سانتی‌متر مکعب از ماده است. این کمیت به همراه “سطح مقطع میکروسکوپی” عناصر، احتمال وقوع واکنش را در هر نقطه از هندسه تعیین می‌کند. بدون داشتن این عدد، کد قادر به تشخیص میزان “فشردگی” هسته‌ها در مسیر حرکت ذره نخواهد بود.

مبانی ریاضی و فیزیکی محاسبه چگالی اتمی ($N$N)

برای محاسبه چگالی اتمی، از رابطه‌ای استفاده می‌شود که چگالی جرمی ماده ($\rho$ρ)، عدد آووگادرو ($N_A$NA​) و جرم اتمی عنصر ($M$M) را به هم پیوند می‌دهد. حاصل این فرمول، تعداد اتم‌ها در واحد حجم است. در محیط‌های پیچیده که از چندین عنصر تشکیل شده‌اند، کاربر باید چگالی اتمی هر عنصر را به صورت مجزا بر اساس کسر جرمی یا اتمی آن در ترکیب کل محاسبه کند. این مقدار معمولاً عددی بسیار بزرگ با توان‌های مثبت است که نشان‌دهنده تراکم عظیم اتم‌ها در مقیاس میکروسکوپی است.

تفاوت بین چگالی جرمی و چگالی اتمی در ورودی‌های MCNP

کد MCNP به کاربر اجازه می‌دهد تا چگالی ماده را در کارت‌های سلول به دو روش تعریف کند. روش اول، استفاده از چگالی جرمی (معمولاً با واحد گرم بر سانتی‌متر مکعب) است که با علامت منفی در فایل ورودی مشخص می‌شود. روش دوم، استفاده از چگالی اتمی است که با علامت مثبت وارد می‌شود. درک این تفاوت بسیار حیاتی است؛ چرا که اگر کاربری چگالی جرمی را بدون علامت منفی وارد کند، کد آن را به عنوان چگالی اتمی (تعداد اتم در واحد حجم) تفسیر کرده و محاسبات با خطایی در ابعاد چندین مرتبه بزرگی (مرتبه توان ۱۰) انجام خواهد شد.

محاسبه چگالی برای ترکیبات شیمیایی و آلیاژها

چالش اصلی زمانی بروز می‌کند که با موادی مانند بتن، فولاد یا سوخت‌های هسته‌ای مواجه هستیم که ترکیبی از چندین عنصر مختلف هستند. در این موارد، ابتدا باید “وزن مولکولی موثر” یا “جرم اتمی میانگین” ترکیب محاسبه شود. سپس با استفاده از درصدهای وزنی (Weight Percent) که معمولاً در استانداردهایی مانند ASTM یافت می‌شوند، چگالی هر جزء به طور دقیق استخراج می‌گردد. دقت در این مرحله بسیار مهم است، زیرا حتی مقادیر ناچیز از عناصری با قدرت جذب بالا (مانند بور یا کادمیوم) در یک آلیاژ، می‌تواند رفتار کل سیستم را تغییر دهد.

واحد “بارن-سانتی‌متر” و اهمیت آن در محاسبات هسته‌ای

در خروجی‌ها و ورودی‌های MCNP، اغلب با واحدی به نام “اتم بر بارن-سانتی‌متر” ($atoms/b-cm$atoms/b−cm) مواجه می‌شویم. از آنجا که واحد سطح مقطع هسته‌ای “بارن” (${10}^{-24}$10−24 سانتی‌متر مربع) است، استفاده از این واحد ترکیبی باعث می‌شود تا اعداد محاسباتی از حالت توان‌دار بسیار بزرگ خارج شده و به اعدادی قابل‌فهم و ساده‌تر (در محدوده ۰.۰۱ تا ۰.۱) تبدیل شوند. تسلط بر تبدیل واحدهای رایج آزمایشگاهی به این واحد استاندارد هسته‌ای، مهارتی است که هر کاربر حرفه‌ای MCNP باید به آن مسلط باشد تا از صحت توصیف مادی محیط اطمینان حاصل کند.