شبیه سازی سطح (100) شبکه fcc پلاتین و طلا به روش مونت-کارلو
نویسندگان
1
2
doi
چکیده
با استفاده از روش مونت ـ کارلو سطح (100) شبکه مکعبی مرکز سطحی فلزات طلا و پلاتین شبیه سازی شده است. معادله پتانسیلی که برای توصیف رفتار اتمهای تشکیل دهنده این سطوح مورد استفاده قرار گرفته, پتانسیل ساتن و چن است. این تابع پتانسیل اختصاصاً برای سیستمهای متشکل از هستههای غوطهور در دریای الکترونی طراحی شده است که در آن انرژی پتانسیل هر اتم به چگالی الکترونی محیطی بستگی دارد که در آن قرار گرفته است. همین بستگی, این معادله را در زمره پتانسیلهای بس ذرهای قرار میدهد. با افزایش دما, اتمهای سازنده سطح به موقعیتهای خالی روی سطح منتقل میشوند. این جابهجاییها منجر به افزایش تعداد ناراستیها, افزایش سطح شبکه و تغییر عدد کوئوردینانسیون میشود. در هر انتقال, فقط یک گام, یعنی انتقال به نزدیکترین موقعیت خالی مجاور مجاز است و همین باعث میشود که سیستم در یک حالت جدید انرژی قرار بگیرد. حرکتهای اولیه منجر به افت و خیزهای شدیدی در انرژی سیستم میشود ولی با گذشت مدت زمان نسبتاً طولانی, یعنی هنگامی که تعداد گامها از مرتبه چند صدهزار شود, تغییرات انرژی سیستم به ازای یک گام خیلی جزئی خواهد شد و با تقریب بسیار خوبی میتوان از افت و خیز در انرژی چشم پوشید و آرایش حاصل را به عنوان حالت پایدار سیستم در دمای مورد بررسی در نظر گرفت. در این محاسبات, برهم کنش هر ذره با اولین و دومین همسایه نزدیک آن در نظر گرفته میشود و انتقالهایی پذیرفته میشوند که انرژی آنها با توجه به ضریب بولتزمن مناسب باشد. در نهایت اندازه سطح را به صورت تابعی از دما به دست آوردیم. نتایج محاسبات نشان میدهند که وقتی قسمتی از سطح از اتم خالی میشود, سایر انتقالها عمدتاً طوری صورت میگیرند که اندازه این حفره بزرگ شود. جهت حصول اطمینان از صحت عملیات شبیه سازی, از معادله شوتکی استفاده شد؛ به این ترتیب که خطی بودن لگاریتم طبیعی نسبت تعداد نقصهای ایجاد شده به کل تعداد اتمها بر حسب معکوس دما مورد تأیید قرار گرفت.