24小时热门版块排行榜    

查看: 1870  |  回复: 11

[交流] 【求助完结】calculation = "relax" 改为 “vc-relax",

我优化一个六面体的结构。在不优化晶胞的情况下可以顺利进行,优化已经结束。但是想优化晶胞试试,结果就报错了,除了calculation = "relax"这个参数,其他均未修改,不知道错在何处。
&CONTROL
  calculation    = "vc-relax",
  prefix         = "h-BN-v",
  nstep          = 500,
  pseudo_dir   = "/home/PWscf/BN",
  outdir       = "/home/PWscf/outdir",
/
&SYSTEM
  ibrav                  = 4,
  celldm(1)              = 4.7318,
  celldm(3)              = 5.02308,
  nat                    = 4,
  ntyp                   = 2,
  ecutwfc                = 120.0D0,
  nspin                  = 1,
/
&ELECTRONS
  conv_thr    = 1.0D-8,
  mixing_beta = 0.3D0,
  electron_maxstep = 400,
/
&IONS
  ds                = 1.D-8,
/
ATOMIC_SPECIES
B  10.811       B.UPF
N  14.00674     N.UPF
ATOMIC_POSITIONS { angstrom }
B 0.00000000    0.00000000    1.66400000
N 1.25200000    0.72284254    1.66400000
N 0.00000000    0.00000000    4.99200000
B 1.25200000    0.72284254    4.99200000
K_POINTS {automatic}
10 10 10 0 0 0

[ Last edited by minmin_0082003 on 2011-1-6 at 18:28 ]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
什么错误总得给出来大家看下啊
2楼2011-01-06 17:17:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
引用回帖:
Originally posted by zzy870720z at 2011-01-06 17:17:32:
什么错误总得给出来大家看下啊

错误如下。不好意思
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
     from  read_namelists  : error #         1
      reading namelist cell
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

     stopping ...

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
     from  read_namelists  : error #         1
      reading namelist cell
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

     stopping ...

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
     from  read_namelists  : error #         1
      reading namelist cell
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

     stopping ...

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
     from  read_namelists  : error #         1
      reading namelist cell
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

     stopping ...
application called MPI_Abort(MPI_COMM_WORLD, 0) - process 3
application called MPI_Abort(MPI_COMM_WORLD, 0) - process 0
application called MPI_Abort(MPI_COMM_WORLD, 0) - process 2
3楼2011-01-06 17:18:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
★
youzhizhe(金币+1):谢谢交流。 2011-01-06 22:43:43
提示读取namelist错误
你加上
&CELL
/
试试
4楼2011-01-06 17:25:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
引用回帖:
Originally posted by zzy870720z at 2011-01-06 17:25:31:
提示读取namelist错误
你加上
&CELL
/
试试

试过,不行,错误相同
5楼2011-01-06 17:31:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
★
youzhizhe(金币+1):谢谢交流。 2011-01-06 22:43:56
引用回帖:
Originally posted by minmin_0082003 at 2011-01-06 17:31:20:

试过,不行,错误相同

不过感觉你的晶格参数有些问题
celldm(3)=5.xx*celldm(1)
好像很大啊
6楼2011-01-06 17:33:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
★ ★
minmin_0082003(金币+10):很奇怪,你修改其他的了吗,我修改了加上cell,还是不行。后面直接拷贝你给我的输入就运行了..... 2011-01-06 18:04:12
youzhizhe(金币+2):谢谢交流。 2011-01-06 22:44:20
引用回帖:
Originally posted by minmin_0082003 at 2011-01-06 17:31:20:

试过,不行,错误相同

我试过了,没有问题,都出几步了
CODE:
&CONTROL
  calculation    = "vc-relax",
  prefix         = "h-BN-v",
  nstep          = 500,
  pseudo_dir   = "./",
  outdir       = "./outdir",
/
&SYSTEM
  ibrav                  = 4,
  celldm(1)              = 4.7318,
  celldm(3)              = 5.02308,
  nat                    = 4,
  ntyp                   = 2,
  ecutwfc                = 120.0D0,
  nspin                  = 1,
/
&ELECTRONS
  conv_thr    = 1.0D-8,
  mixing_beta = 0.3D0,
  electron_maxstep = 400,
/
&IONS
  ds                = 1.D-8,
/
&CELL
/
ATOMIC_SPECIES
B  10.811       B.pw91-n-van_ak.UPF
N  14.00674   N.pw91-van_ak.UPF
ATOMIC_POSITIONS { angstrom }
B 0.00000000    0.00000000    1.66400000
N 1.25200000    0.72284254    1.66400000
N 0.00000000    0.00000000    4.99200000
B 1.25200000    0.72284254    4.99200000
K_POINTS {automatic}
10 10 10 0 0 0

7楼2011-01-06 17:55:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
引用回帖:
Originally posted by zzy870720z at 2011-01-06 17:33:08:


不过感觉你的晶格参数有些问题
celldm(3)=5.xx*celldm(1)
好像很大啊

谢谢,cell3没有单位了,晕乎乎了,想成是Ang直接又除以0.5....转化成波尔了,灰常感谢
8楼2011-01-06 18:07:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
[quote]Originally posted by zzy870720z at 2011-01-06 17:55:01:



我想再废话一句
难道Cell加的位置不一样会导致出错。我把cell加在&ELECTRONS后来着
9楼2011-01-06 18:09:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
youzhizhe(金币+1):谢谢交流。 2011-01-06 22:44:39
引用回帖:
Originally posted by minmin_0082003 at 2011-01-06 18:09:08:
[quote]Originally posted by zzy870720z at 2011-01-06 17:55:01:



我想再废话一句
难道Cell加的位置不一样会导致出错。我把cell加在&ELECTRONS后来着

我就添加了cell,其他的修改了一下目录
就可以跑了
你在windows下建立的吗?
dos2linux一下文件试试
10楼2011-01-06 22:20:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
引用回帖:
Originally posted by zzy870720z at 2011-01-06 22:20:17:


我就添加了cell,其他的修改了一下目录
就可以跑了
你在windows下建立的吗?
dos2linux一下文件试试

建立输入吗?我都是手写的,在linux下,话说在windows下怎么弄?是不是安装哪个pwgui?我装过,出现错误信息,就罢了。呵呵
我试过了就是因为我放的cell的位置不对导致的
11楼2011-01-07 17:32:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
gzqdyouxia(金币+1):谢谢回复 2011-01-07 22:22:23
引用回帖:
Originally posted by minmin_0082003 at 2011-01-07 17:32:00:

建立输入吗?我都是手写的,在linux下,话说在windows下怎么弄?是不是安装哪个pwgui?我装过,出现错误信息,就罢了。呵呵
我试过了就是因为我放的cell的位置不对导致的

呵呵
我还以为你是windows输入的,可以用pwgui试试
那个建输入挺好的,恭喜你解决问题
12楼2011-01-07 22:02:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 minmin_0082003 的主题更新
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 2rKfVQg2jQww 2026-11-06 15/750 2026-11-09 05:37 by l48hFd8Nwae2
[论文投稿] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 5sJFDKp2kwuI 2026-11-06 12/600 2026-11-09 05:26 by l48hFd8Nwae2
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 5sJFDKp2kwuI 2026-11-06 13/650 2026-11-09 05:22 by l48hFd8Nwae2
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 2rKfVQg2jQww 2026-11-06 11/550 2026-11-09 00:26 by l48hFd8Nwae2
[公派出国] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 5sJFDKp2kwuI 2026-11-06 8/400 2026-11-09 00:14 by l48hFd8Nwae2
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +8 ncTVBM6AOJHF 2026-11-04 23/1150 2026-11-07 06:49 by Equinoxhua
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +10 ncTVBM6AOJHF 2026-11-04 22/1100 2026-11-06 17:18 by 5sJFDKp2kwuI
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +6 ELBZ6WwEaRj6 2026-11-04 11/550 2026-11-06 16:32 by 5sJFDKp2kwuI
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +6 ncTVBM6AOJHF 2026-11-04 15/750 2026-11-06 13:49 by BkA9Bi5cMhkD
[考研] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +7 ncTVBM6AOJHF 2026-11-04 14/700 2026-11-06 13:47 by BkA9Bi5cMhkD
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +8 ncTVBM6AOJHF 2026-11-04 12/600 2026-11-06 13:47 by BkA9Bi5cMhkD
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +8 ncTVBM6AOJHF 2026-11-04 18/900 2026-11-06 13:36 by BkA9Bi5cMhkD
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 ELBZ6WwEaRj6 2026-11-05 5/250 2026-11-06 09:01 by 2lRNluKrJpwL
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 ELBZ6WwEaRj6 2026-11-04 12/600 2026-11-06 05:32 by qL1cU46YIdK0
[考博] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 kvewwGtKktqB 2026-11-04 11/550 2026-11-06 05:10 by qL1cU46YIdK0
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 kvewwGtKktqB 2026-11-04 14/700 2026-11-06 05:10 by qL1cU46YIdK0
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 ELBZ6WwEaRj6 2026-11-04 9/450 2026-11-06 02:15 by qL1cU46YIdK0
[考研] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 ttH3f1oA6ace 2026-11-04 10/500 2026-11-06 02:04 by qL1cU46YIdK0
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 ttH3f1oA6ace 2026-11-04 9/450 2026-11-06 02:02 by qL1cU46YIdK0
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +6 ncTVBM6AOJHF 2026-11-04 10/500 2026-11-06 01:51 by qL1cU46YIdK0
信息提示
请填处理意见