| 查看: 1941 | 回复: 4 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
chenxiaoqian新虫 (小有名气)
|
[求助]
我在计算热力学，算热膨胀系数的时候遇到问题
|
||
|
我在用VASP计算热膨胀系数，输入./script-qha.sh -r 提交工作给GridEngine工作系统，可是计算特别慢，并且一直处于如图的状态，想知道是什么原因，因为是第一次计算这个，不知道这种情况正常不，求助各位大神~~ PS：我用的是学校的计算中心，一般计算是提交任务，不知道会不会有影响。  V8R)M}K(2YDACJ7T[{0`R8B.png |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有9人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有6人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有6人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有5人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有6人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有8人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有6人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有9人回复
-
国自然面上复盘~欢迎讨论
已经有15人回复
-
中青基了要发朋友圈吗？
已经有3人回复
chenxiaoqian
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 659.4
- 帖子: 152
- 在线: 76.5小时
- 虫号: 2952102
- 注册: 2014-01-31
- 专业: 原子和分子物理
|
这是script-qha.sh 文件 [chenq@su01 qha]$ cat script-qha.sh #!/bin/bash dim="1 1 1" pa="0 1/2 1/2 1/2 0 1/2 1/2 1/2 0" mp="3 3 3" nunit=32 while getopts ':mrc' opt; do case $opt in m) for i in `/bin/ls QHAPOSCAR*`;do qhaname=`echo $i|sed s/POSCAR//` mkdir $qhaname mv $i $qhaname/POSCAR cd $qhaname pwd phonopy -d --dim="$dim" for j in `/bin/ls POSCAR-*`;do dispname=`echo $j|sed s/POSCAR/disp/` mkdir $dispname mv $j $dispname/POSCAR cp ../{INCAR,KPOINTS,POTCAR} $dispname done mkdir perfect cp ../{INCAR,KPOINTS,POTCAR} perfect cp SPOSCAR perfect/POSCAR cd .. done ;; r) for i in `/bin/ls -d QHA-*`;do cd $i pwd vol=`echo $i|sed s/QHA-//` # For displacements for j in `/bin/ls -d disp-*`;do cd $j pwd num=`echo $j|sed s/disp-//` sed s/num/$num/ ../../calc.sh|sed s/vol/$vol/|mpirun -np 20 vasp cd .. done # For perfect cd perfect sed s/num/per/ ../../calc.sh|sed s/vol/$vol/|mpirun -np 20 vasp cd .. cd .. done ;; c) echo "# cell volume energy of cell other than phonon" > e-v.dat for i in `/bin/ls -d QHA-*`;do cd $i phonopy -f disp-*/vasprun.xml phonopy --mp="$mp" -t --dim="$dim" --pa="$pa" --tmax=2004 --tstep=2 volume=`grep volume perfect/vasprun.xml|tail -n 1|awk -F'<|>' '{printf("%20.13f", $3)}'` energy=`grep e_wo_entrp perfect/vasprun.xml|tail -n 1|awk -F'<|>' '{printf("%20.13f", $3)}'` cd .. echo `echo "$volume/$nunit"|bc -ls` `echo "$energy/$nunit"|bc -ls` >> e-v.dat done ;; esac done 这是calc.sh文件 #$ -S /bin/zsh #$ -cwd #$ -N Sivol-num #$ -pe mpi* 8 #$ -e err.log #$ -o std.log |
3楼2017-09-21 09:34:01
zzggbb
铁杆木虫 (著名写手)
- 应助: 2 (幼儿园)
- 金币: 13790.4
- 帖子: 1022
- 在线: 169.2小时
- 虫号: 484955
- 注册: 2007-12-27
- 专业: 凝聚态物性 II ：电子结构
2楼2017-09-21 06:23:14
小虫虫LG
木虫 (正式写手)
- 应助: 114 (高中生)
- 金币: 1775.9
- 散金: 1871
- 红花: 25
- 帖子: 975
- 在线: 678.5小时
- 虫号: 2128430
- 注册: 2012-11-15
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
4楼2017-09-21 14:53:59
卡开发发
专家顾问 (著名写手)
Ab Initio Amateur
-
专家经验: +224 - 1ST强帖: 2
- 应助: 47 (小学生)
- 金币: 3567.8
- 散金: 14
- 红花: 223
- 帖子: 2477
- 在线: 1254小时
- 虫号: 1369423
- 注册: 2011-08-16
- 性别: GG
- 专业: 金属材料的磨损与磨蚀
- 管辖: 第一性原理
5楼2017-09-24 08:23:11
10