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

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

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

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

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

.xml (در قالب XML

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

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

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

# gnuplot –persist charge.gplot 

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

# gnuplot –persist pseudo.gplot

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

# gnuplot –persist pots.gplot

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

# gnuplot –persist scrpots.gplot

-----------------------------------

منبع : http://www.comphys.ir