Tính mật độ điện tích

Để chạy được tập lệnh tính mật độ điện tích, ta cần chạy chương trình pw.x để tính toán hàm sóng tự hợp trước. Lệnh tính hàm sóng:

pw.x < si.scf.in > si.scf.out

cũng như cấu trúc file lệnh si.scf.in đã mô tả trong bài viết về tính toán hàm sóng. Một khi thư mục dữ liệu outdir=‘../tmp/’ đã xuất hiện, ta đã có thể sử dụng file si.pp.in để tính mật độ điện tích theo cú pháp:

pp.x < si.pp.in > si.pp.out

Nội dung điển hình của tệp chứa lệnh tính mật độ điện tích (.pp.in):

&INPUTPP
outdir='../tmp/',
prefix='si',
plot_num=0,
/
&PLOT
iflag=3,
output_format=6,
fileout='si_rho.cube',
nx=64,ny=64,nz=64,
/

Sau khi tính xong sẽ xuất hiện file .pp.out báo tin về quá trình tính toán và file kết quả .cube, .pp. Ta có thể quan sát kết quả trực quan qua file .cube bằng phần mềm VESTA.

Quan sát mật độ điện tích qua phần mềm Vesta.