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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » 新手求助pwscf的进行磁性运算和LDA+U运算
北京石油化工学院2026年研究生招生接收调剂公告
362/4返回列表上一页1234下一页
查看: 4231  |  回复: 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下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿郑州大学085600求调剂 +17 吃的不少 2026-04-05 20/1000 2026-04-06 00:32 by T可可西里T
[考研] 机械专硕274求调剂，不挑专业学校 +6 泛泛2333 2026-04-05 7/350 2026-04-05 23:20 by chyhaha
[考研] 求调剂 +6 chenxrlkx 2026-04-05 6/300 2026-04-05 22:05 by imissbao
[考研] 考研生物与医药调剂 +8 铁憨憨123425 2026-03-31 8/400 2026-04-05 20:28 by lys0704
[考研] 0703调剂 +11 拾玖壹 2026-04-04 12/600 2026-04-05 10:29 by 果冻大王
[考研] 求 +5 化工专硕323分 2026-04-04 5/250 2026-04-05 08:02 by 544594351
[考研] 266分，一志愿电气工程，本科材料，求材料专业调剂 +11 哇呼哼呼哼 2026-04-01 12/600 2026-04-04 23:17 by 永字号
[考研] 材料调剂 +15 一样YWY 2026-04-01 15/750 2026-04-04 22:23 by hemengdong
[考研] 材料科学与工程调剂 +19 深V宿舍吧 2026-03-30 20/1000 2026-04-04 22:13 by hemengdong
[考研] 283分求调剂 +7 小聂爱学习 2026-04-03 7/350 2026-04-04 21:51 by hemengdong
[考研] 环境科学与工程334分求调剂 +9 王一一依依 2026-03-30 12/600 2026-04-04 20:55 by dongzh2009
[考研] 一志愿华北电力大学（北京），材料科学与工程学硕265，求调剂 +11 yelck 2026-04-03 12/600 2026-04-04 19:52 by dongzh2009
[考研] 考研调剂 +4 美丽的youth_ 2026-04-04 5/250 2026-04-04 17:16 by imissbao
[考研] 考研调剂 +3 15615482637 2026-04-03 3/150 2026-04-03 22:50 by ms629
[考研] 求材料调剂 一志愿南昌大学 328分 +5 yyy..... 2026-04-03 5/250 2026-04-03 13:46 by 百灵童888
[考研] 农学考研求调剂 +3 dkdkxm 2026-04-01 3/150 2026-04-02 16:04 by wangjagri
[考研] 0710生物学求调剂 +9 manman511 2026-04-01 9/450 2026-04-02 10:00 by zxl830724
[考研] 一志愿北京科技，085601总分305求调剂 +9 半生瓜！ 2026-04-01 11/550 2026-04-02 08:28 by Wang200018
[考研] 求调剂 生物学 377分 +6 zzll03 2026-03-31 6/300 2026-03-31 17:33 by 唐沐儿
[考研] 080500-315分复试调剂 +9 上岸3821 2026-03-31 9/450 2026-03-31 17:29 by 唐沐儿
信息提示
关闭
请填处理意见
关闭