آموزش کد محاسباتی سیستا (چگونه یک فایل شبه پتانسیل تولید کنیم)

در این آموزش قصد داریم شبه پتانسیل مورد نظر خودمون برای کد سیستا رو با استفاده از کد ATOM تولید کنیم. برای اینکار نیاز داریم که یک فایل ورودی به صورت زیر آماده کنیم. 

بعد از آماده کردن فایل ورودی دستورات زیر رو اجرا می کنیم:

 

دستور اول اسکریپتی از قبل آماده شده رو اجرا می کنه که دستورات مورد نیاز از قبل در اون قرار داده شدن. فایل های شبه پتانسیل در همون دایرکتوری مادر با پسوند های زیر قرار داده میشن:

.vps (بدون قالب بندی)

.psd (قالب بندی شده)

.xml (در قالب XML

فایل های خروجی مختلف در یک دایرکتوری جدید قرار داده میشن (مشابه نام همون فایل ورودی بدون پسوند .inp )

حال میخوایم نمودار های چگالی بار شبه پتانسیل و تمام الکترونی رو رسم کنیم و اون ها رو با هم مقایسه کنیم:

برای رسم نمودار می تونیم از دستور زیر استفاده کنیم:

# gnuplot –persist charge.gplot 

می تونیم از خروجی های بدست اومده نمودار های دیگه ای هم رسم کنیم و اطلاعات مفیدی بدست بیاریم برای این کار می تونیم از دستور زیر استفاده کنیم:

# gnuplot –persist pseudo.gplot

می تونیم می تونیم این نمودار ها رو برای عدد کوانتومی بالاتر هم رسم کنیم. برای رسم شبه پتانسیل ها به ازای عدد کوانتومی مختلف در فضای حقیقی می تونیم از دستور زیر استفاده کنیم:

# gnuplot –persist pots.gplot

همچنین می تونیم با استفاده از دستور زیر شبه پتانسیل های استتار شده و استتار نشده رو رسم کنیم:

# gnuplot –persist scrpots.gplot

 

 


بنده دانشجوی دکترای فیزیک ماده چگال از دانشگاه تربیت مدرس تهران هستم. حوزه مورد علاقه من فیزیک محاسباتی (به طور خاص نظریه تابعیت چگالی) و همچنین سیستم های توپولوژیک است.


پاسخ دهید

نشانی ایمیل شما منتشر نخواهد شد. بخش‌های موردنیاز علامت‌گذاری شده‌اند *