计算化学公社

标题: 求助一个测试K点和encut收敛性的脚本 [打印本页]

---

作者
Author: wjh1122    时间: 2023-10-21 17:04
标题: 求助一个测试K点和encut收敛性的脚本
求助各位大佬一个测试K点和ENCUT收敛性的脚本

---

作者
Author: fux    时间: 2023-10-24 21:01


ENCUT:

1. #PBS -N test
   
2. #PBS -l nodes=1:ppn=8
   
3. #PBS -q high
   
4. #PBS -j oe
   
5. #PBS -l walltime=1000:00:00  
   
6. 
   
7. ulimit -s unlimited
   
8. cd $PBS_O_WORKDIR
   
9. NP=`cat $PBS_NODEFILE|wc -l`
   
10. 
    
11. source /public/software/profile.d/conda.sh
    
12. source /public/software/profile.d/hdf5.sh
    
13. source /public/software/profile.d/compiler_intel-compiler-2017.5.239.sh
    
14. VASP="/public/software/mpi/openmpi/2.1.3/bin/mpirun --mca btl_tcp_if_include enp24s0f0 --machinefile $PBS_NODEFILE -np $NP /public/software/apps/vasp/vasp-6/vasp_std >&vasp.log"
    
15. 
    
16. echo 'ENCUT/eV,energy/eV'>finish.csv
    
17. 
    
18. for i in 400 440 480 520 560 600 640 680
    
19. do
    
20. mkdir $i
    
21. cd $i
    
22. cp ../INCAR .
    
23. cp ../PO* .
    
24. cp ../KPOINTS .
    
25. sed -i "s/^ENCUT = /&$i/g" INCAR
    
26. $VASP
    
27. if [ -e OUTCAR ];then
    
28.    count=`grep "without" OUTCAR | tail -1 | awk -v FS=" " '{print $4}'`
    
29.    echo $i','$count >>../finish.csv
    
30. fi
    
31. cd ..
    
32. done

复制代码

有些细节可能要你自己动手改一下，k点测试也类似就是INCAR一致用不同的kmesh就行

---

欢迎光临 计算化学公社 (http://bbs.keinsci.com/)

Powered by Discuz! X3.3