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.
