ریاضیات

تحلیل ریاضی برد میانگین ذرات آلفا

در این نوشته تعاریف، روابط پایه و روش محاسبهٔ برد یک ذرهٔ آلفا و نیز نحوهٔ به‌دست‌آوردن برد میانگین برای یک توزیع انرژی ارائه شده است — به زبان ساده و همراه با مثال عددی.

۱) تعریف برد

برد \(R\) برای یک ذرهٔ آلفا با انرژی اولیه \(E_0\) معمولاً در تقریب توقف پیوسته (CSDA) به‌صورت زیر تعریف می‌شود:

R(E_0) = ∫₀^{E₀} dE / |dE/dx (E)|

که در آن \(|dE/dx|\) توان توقف خطی (میزان افت انرژی به ازای واحد طول) است.

۲) تقریب‌های رایج برای توان توقف و رابطهٔ تجربی برد

برای ذرات سنگین باردار مانند آلفا، رابطهٔ بتّه (تقریباً) نشان می‌دهد که توان توقف با توان معکوس سرعت وابسته است. در عمل برای ساده‌سازی اغلب از تقریب تجربی «قانون توان» بین برد و انرژی استفاده می‌شود:

R(E) ≈ a · Eⁿ

که پارامترهای تجربی \(a\) و \(n\) به مادهٔ هدف و واحدهای انتخابی بستگی دارند. در مواد معمول، مقدار \(n\) برای آلفا معمولاً بین ۱.۳ تا ۱.۸ قرار دارد.

۳) برد میانگین برای توزیع انرژی

اگر ذرات آلفا دارای توزیع انرژی با چگالی احتمال \(f(E)\) باشند (نرمالیزه شده)، برد میانگین برابر است با:

R̄ = ∫₀^{∞} R(E) · f(E) dE

با استفاده از تقریب \(R(E)=aE^{n}\) این رابطه به شکل زیر ساده می‌شود:

R̄ = a · ∫₀^{∞} E^{n} f(E) dE

یعنی برد میانگین برابرِ پارامتر a ضرب در ممان (moment) مرتبه nِ توزیع انرژی است.

۴) چند مثال تحلیلی ساده

• توزیع δ (تک‌انرژی): اگر همهٔ ذرات انرژی ثابت \(E_0\) داشته باشند (f(E)=δ(E−E₀)) آنگاه R̄ = a·E₀ⁿ.
• توزیع نمایی: برای \(f(E)=λ e^{−λE}\) (E≥0) و با \(R(E)=aE^{n}\) داریم:

  R̄ = a · λ^{-n} · Γ(n+1)

  که Γ تابع گاما است.

• توزیع تقریب گاوسی: اگر توزیع تقریباً گاوسی با میانگین μ و واریانس σ² باشد و σ≪μ، می‌توان از بسط تیلور استفاده کرد:

  R̄ ≈ a[ μⁿ + (n(n−1)/2) μ^{n−2} σ² + … ]

مثال عددی (نمونه):

فرض کنید پارامترها به‌صورت زیر باشند: \(a=0.01\)، \(n=1.5\)، و توزیع انرژی نمایی با λ=0.2.

با استفاده از رابطهٔ نمایی می‌نویسیم:

R̄ = a · λ^{-n} · Γ(n+1)

که به‌صورت عددی می‌شود (حدود):

R̄ ≈ 0.1486 (یکا وابسته به a و E)

توجه کنید که این مقدار نمونه‌ای است و مقادیر حقیقی بسته به مادهٔ هدف، چگالی و واحدها تغییر می‌کند.

۵) نکات عملی

• در گزارش‌های تجربی معمولاً برد را به صورت جرم بر سطح (g/cm²) نیز می‌دهند تا اثر چگالی ماده جدا شود: \(R_{\text{(g/cm²)}} = ρ · R_{\text{(cm)}}\).
• برای محاسبات دقیق از داده‌های CSDA یا شبیه‌سازهای اختصاصی (مثل SRIM) استفاده می‌شود؛ در انرژی‌های پایین نیز تصحیحات اضافی لازم است (تعاملات شیمیایی، تبادل بار، و غیره).
• اگر طیف انرژی منبع را دارید (اندازه‌گیری یا مدل)، بهترین روش عددی این است که تابع R(E) معتبر برای مادهٔ هدف را بردارید و انتگرال \(∫ R(E) f(E) dE\) را با روش عددی محاسبه کنید.

تذکر: در صورتی که بخواهید من برایتان یک نسخهٔ آمادهٔ HTML با فرمت‌بندی دقیق‌تر (مثلاً شامل تصاویری از نمودارهای R(E) یا جدول‌های عددی) بسازم، بگویید — من همین الآن آن را تهیه می‌کنم.