应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » 新手求助pwscf的进行磁性运算和LDA+U运算
362/4返回列表上一页1234下一页
查看: 3999  |  回复: 35
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

huazhorg

铁杆木虫 (正式写手)
引用回帖:
10楼: Originally posted by 刘仕晨 at 2013-03-25 22:21:40
谢谢，谢谢，受益匪浅，对我这个新手来说收获很多，我现在怀疑我的K点也设的不够，所以我现在在测试K点，画出K点和能量的关系，然后拿那个K点再来加大ecut，谢谢
那在加U同时开启spin，那这些参数呢Hubbard_alpha ...

一般默认就可以。你也可以参照一下example25, 这就是一个spin+U的例子。
赞一下 回复此楼
11楼2013-03-26 07:51:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

刘仕晨

木虫 (正式写手)
引用回帖:
11楼: Originally posted by huazhorg at 2013-03-26 07:51:03
一般默认就可以。你也可以参照一下example25, 这就是一个spin+U的例子。...

谢谢，谢谢，我现在想看下ecut加到多少合适，正在以10往上加，然后每个ecut加密K点再计算，现在多算算，好多学点，先算到您建议的90。那个模守恒时，我看手册里写ecutrho是4倍的ecutwfc，超软是8-12倍，是这样的吧，那PAW呢，PAW是多少倍，因为我们用vasp计算的时候正常ecut是600，所以没想到要这么大
赞一下 回复此楼
12楼2013-03-26 14:28:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huazhorg

铁杆木虫 (正式写手)
引用回帖:
12楼: Originally posted by 刘仕晨 at 2013-03-26 14:28:59
谢谢，谢谢，我现在想看下ecut加到多少合适，正在以10往上加，然后每个ecut加密K点再计算，现在多算算，好多学点，先算到您建议的90。那个模守恒时，我看手册里写ecutrho是4倍的ecutwfc，超软是8-12倍，是这样的吧 ...

你按手册上来就没错。从文献上看，PAW的截断不一定比USPP的高。另外，espresso里的PAW赝势似乎比较少。
赞一下 回复此楼
13楼2013-03-26 19:09:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

刘仕晨

木虫 (正式写手)
引用回帖:
13楼: Originally posted by huazhorg at 2013-03-26 19:09:58
你按手册上来就没错。从文献上看，PAW的截断不一定比USPP的高。另外，espresso里的PAW赝势似乎比较少。...

您好，我在测试ecut和K点的时候遇到了新的问题，我是这样做的，分别设ecut=30 40 50 60 70 80 90 100 110 120,然后在每个ecut值下面慢慢增加K点，但是奇怪的是，在ecut=50 60和70的时候自洽运算会出现对称性操作丢失的警告：
    warning: symmetry operation #  3 not allowed.   fractional translation:
       0.5000000 -0.5000000  0.5000000  in crystal coordinates
     warning: symmetry operation #  4 not allowed.   fractional translation:
      -0.5000000  0.5000000  0.5000000  in crystal coordinates
     warning: symmetry operation #  7 not allowed.   fractional translation:
      -0.5000000  0.5000000 -0.5000000  in crystal coordinates
     warning: symmetry operation #  8 not allowed.   fractional translation:
       0.5000000 -0.5000000 -0.5000000  in crystal coordinates
     warning: symmetry operation # 11 not allowed.   fractional translation:
      -0.5000000  0.5000000 -0.5000000  in crystal coordinates
     warning: symmetry operation # 12 not allowed.   fractional translation:
       0.5000000 -0.5000000 -0.5000000  in crystal coordinates
     warning: symmetry operation # 15 not allowed.   fractional translation:
       0.5000000 -0.5000000  0.5000000  in crystal coordinates
     warning: symmetry operation # 16 not allowed.   fractional translation:
      -0.5000000  0.5000000  0.5000000  in crystal coordinates
我看到小木虫上有说通过增加ecut可以解决，可是为什么30 和40又能跑呢？这个问题怎么解决，我同时贴出我的IN文件：
&CONTROL
                 calculation = 'scf' ,
                restart_mode = 'from_scratch' ,
                      outdir = '/home/liushichen/work/pwscf/tmp/' ,
                  pseudo_dir = '/home/liushichen/tools/codes/espresso-4.3.2/pseudo/' ,
                      prefix = 'CrO2' ,
/
&SYSTEM
                       ibrav = 0,
                         nat = 6,
                        ntyp = 2,
                     ecutwfc = 50 ,
                     ecutrho = 200 ,
                 occupations = 'smearing' ,
                     degauss = 0.02 ,
                    smearing = 'gaussian' ,
                       nspin = 2 ,
   starting_magnetization(1) = 0.5,
                  lda_plus_u = .true. ,
                Hubbard_U(1) = 3,
/
&ELECTRONS
                    conv_thr = 1.0e-10 ,
                 mixing_beta = 0.7 ,
/
CELL_PARAMETERS cubic
     8.350300000    0.000000000    0.000000000
     0.000000000    8.350300000    0.000000000
     0.000000000    0.000000000    5.502600000
ATOMIC_SPECIES
   Cr   52.00000  Cr.pz-hgh.UPF
    O   16.00000  O.pz-hgh.UPF
ATOMIC_POSITIONS crystal
   Cr      0.000000000    0.000000000    0.000000000   
   Cr      0.500000000    0.500000000    0.500000000   
    O      0.303000000    0.303000000    0.000000000   
    O      0.697000000    0.697000000    0.000000000   
    O      0.803000000    0.197000000    0.500000000   
    O      0.197000000    0.803000000    0.500000000   
K_POINTS automatic
  6 6 9   0 0 0
2，第二个问题是：我发现计算结果的out文件中的原子位置的Z方向有移动，但是我没有relax阿，我只是跑的自洽阿，这是为什么，贴出那部分文件：
Cartesian axes

     site n.     atom                  positions (alat units)
         1           Cr  tau(   1) = (   0.0000000   0.0000000   0.0000000  )
         2           Cr  tau(   2) = (   0.5000000   0.5000000   0.3294852  )
         3           O   tau(   3) = (   0.3030000   0.3030000   0.0000000  )
         4           O   tau(   4) = (   0.6970000   0.6970000   0.0000000  )
         5           O   tau(   5) = (   0.8030000   0.1970000   0.3294852  )
         6           O   tau(   6) = (   0.1970000   0.8030000   0.3294852  )
麻烦帮我看下，谢谢！
赞一下 回复此楼
14楼2013-03-27 15:39:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huazhorg

铁杆木虫 (正式写手)
NCPP要求高截断，30，40即使能跑也没什么意义，50，60，70测试结果无论出不出问题，都没什么意义。因为这是NCPP。
赞一下 回复此楼
15楼2013-03-27 21:15:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

刘仕晨

木虫 (正式写手)
引用回帖:
15楼: Originally posted by huazhorg at 2013-03-27 21:15:43
NCPP要求高截断，30，40即使能跑也没什么意义，50，60，70测试结果无论出不出问题，都没什么意义。因为这是NCPP。

那我的第二个问题，您知道么，这是什么情况
赞一下 回复此楼
16楼2013-03-27 22:21:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huazhorg

铁杆木虫 (正式写手)
引用回帖:
16楼: Originally posted by 刘仕晨 at 2013-03-27 22:21:49
那我的第二个问题，您知道么，这是什么情况...

结构有没有变形，看这几行数字还不能判别，因为这几行数字的含义还不知道。用xcrysden看，就会一目了然了。
赞一下 回复此楼
17楼2013-03-28 08:55:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gemucai

木虫 (正式写手)
引用回帖:
14楼: Originally posted by 刘仕晨 at 2013-03-27 15:39:27
您好，我在测试ecut和K点的时候遇到了新的问题，我是这样做的，分别设ecut=30 40 50 60 70 80 90 100 110 120,然后在每个ecut值下面慢慢增加K点，但是奇怪的是，在ecut=50 60和70的时候自洽运算会出现对称性操作丢 ...

你用8.3503乘以0.3294852，看看跟5.5026乘以0.5是不是一样的，就是给你换了坐标而已。
赞一下(1人) 回复此楼
18楼2013-03-28 10:23:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

刘仕晨

木虫 (正式写手)
引用回帖:
18楼: Originally posted by gemucai at 2013-03-28 10:23:11
你用8.3503乘以0.3294852，看看跟5.5026乘以0.5是不是一样的，就是给你换了坐标而已。...

谢谢，这么一说我就清楚了，我看了手册的，上面说以晶格常数a为单位，但是我一看就Z方向变了，表示很想不通，不过现在这么一说我突然想通了，因为a,b常数本来就是一样的，都是8.3503，根本不用变的，谢谢，太感谢你了！
赞一下 回复此楼
19楼2013-03-28 13:55:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

刘仕晨

木虫 (正式写手)
引用回帖:
18楼: Originally posted by gemucai at 2013-03-28 10:23:11
你用8.3503乘以0.3294852，看看跟5.5026乘以0.5是不是一样的，就是给你换了坐标而已。...

您好，我想问下，pwscf计算band时KPOINTS的高对称点的选择，一定要一个一个把所有的点都写出来么，能不能像vasp里面的那样，只给出高对称点和告诉它之间产生的点数，自动生成呢，比如像vasp里这样的：
k-points along high symmetry lines
10 ! 10 intersections
Line-mode
rec
0 0 0         ! gamma
0.5 0.5 0   ! X
0.5 0.5 0  ! X
0.5 0.75 0.25 ! W
0.5 0.75 0.25 ! W
0 0 0         ! gamma
谢谢了，麻烦你 了。
赞一下 回复此楼
20楼2013-04-22 13:27:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 刘仕晨 的主题更新
362/4返回列表上一页1234下一页
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭