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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » VASP对于极性薄膜的计算，薄膜不带电，但电荷分布不均匀，如何设置DIPOLE的相关参数
91/1返回列表
查看: 2244  |  回复: 8
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

linda19918

银虫 (小有名气)

[求助] VASP对于极性薄膜的计算，薄膜不带电，但电荷分布不均匀，如何设置DIPOLE的相关参数 已有1人参与

slab体系整体是不带电的，但是因为电荷分布不均匀，使模型有偶极矩，应该也是要偶极校正的吧？偶极修正的目的是什么，是修正所建立体系的电偶极矩，还是由于周期性，体系与体系之间的电偶极矩？
回复此楼

» 收录本帖的淘帖专辑推荐

第一性原理相关文档 VASP

» 猜你喜欢
1楼 2015-08-14 14:29:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

枫寒

铁虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与，应助指数 +1
linda19918: 金币+5, ★★★很有帮助 2015-08-15 14:51:44
IDIPOL = 3         
LDIPOL = .TRUE.
赞一下 回复此楼
2楼2015-08-14 14:52:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

linda19918

银虫 (小有名气)
引用回帖:
2楼: Originally posted by 枫寒 at 2015-08-14 14:52:35
IDIPOL = 3         
LDIPOL = .TRUE.

能解释一下它的作用是什么吗？说明书上说得不是很清楚，它说要pre-converging the orbital without LDIPOL-tag,还有就是celll 要足够大是指厚度么？
赞一下 回复此楼
3楼2015-08-14 15:00:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

linda19918

银虫 (小有名气)
引用回帖:
2楼: Originally posted by 枫寒 at 2015-08-14 14:52:35
IDIPOL = 3         
LDIPOL = .TRUE.

没有用诶，加了这个参数能量还是发散的都不收敛
赞一下 回复此楼
4楼2015-08-16 15:00:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

枫寒

铁虫 (小有名气)
引用回帖:
3楼: Originally posted by linda19918 at 2015-08-14 15:00:50
能解释一下它的作用是什么吗？说明书上说得不是很清楚，它说要pre-converging the orbital without LDIPOL-tag,还有就是celll 要足够大是指厚度么？...

加极化校正的原因是：以前大家建模型是中间固定，两面对称（就比如正面吸附一个分子，反面同等位置也要吸附相同的分子），现在建模型都是一边吸附，这就带来了部分偶极距，所以要加极化校正。 IDIPOL=3意思是沿着z轴方向加校正。具体你为什么能量不收敛我不知道。可能跟极化校正没关系。不知道你算的什么体系。
赞一下 回复此楼
5楼2015-08-16 17:10:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur
因为采用Slab模型，但采用的计算方法还是计算固体（3D）的方法，边界上是周期化的（电势计算过程很自然引入Ewald求和，如果不是这样的边界条件这个问题自然不存在，当然也就没有办法使用平面波来进行处理）。这样，如果把Slab模型沿着Z方向建立为超胞，当Slab体系上下不对称的时候，靠上一层的底面和靠下一层的顶面发生作用（形成偶极）。但实际这种作用并不存在，所以需要通过抵消这种作用进行校正，其中一种方法就是自洽过程中加上一个与这种电势差等大的电势来抵消。

如果SCF考虑偶极校正的话还难以收敛，应当考虑磁矩、结构等其他是否设置正确。
赞一下(1人) 回复此楼
不一定挂在论坛，计算问题问题欢迎留言。
6楼2015-08-16 17:34:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

linda19918

银虫 (小有名气)
引用回帖:
6楼: Originally posted by 卡开发发 at 2015-08-16 17:34:54
因为采用Slab模型，但采用的计算方法还是计算固体（3D）的方法，边界上是周期化的（电势计算过程很自然引入Ewald求和，如果不是这样的边界条件这个问题自然不存在，当然也就没有办法使用平面波来进行处理）。这样， ...

结构肯定是对的，计算的能量是这样的

。。。。
   free energy    TOTEN  =    -27013.50693412 eV
  free energy    TOTEN  =    -26816.92597603 eV
  free energy    TOTEN  =    -27176.35064060 eV
  free energy    TOTEN  =    -91456.19775036 eV
  free energy    TOTEN  =   -103242.14612303 eV
  free energy    TOTEN  =    -67400.54656992 eV
  free energy    TOTEN  =    -77627.93362477 eV
  free energy    TOTEN  =    -38063.32893445 eV
  free energy    TOTEN  =    -18288.43756682 eV
  free energy    TOTEN  =     -1932.74943813 eV
  free energy    TOTEN  =      1674.36195917 eV
  free energy    TOTEN  =    -46267.68872523 eV
  free energy    TOTEN  =    -21604.74777661 eV
  free energy    TOTEN  =    -31074.81302005 eV
  free energy    TOTEN  =    -30733.55777787 eV
  free energy    TOTEN  =    -82355.51252433 eV
  free energy    TOTEN  =    -39593.83933052 eV
  free energy    TOTEN  =    -42439.71113435 eV
  free energy    TOTEN  =    -36012.97209928 eV
  free energy    TOTEN  =    -40510.61434933 eV
  free energy    TOTEN  =    -34462.10154152 eV
  free energy    TOTEN  =    -39536.96984058 eV
  free energy    TOTEN  =    -33941.05344567 eV
  free  energy   TOTEN  =    -33941.053446 eV
  free energy    TOTEN  =    -68146.30470837 eV
  free energy    TOTEN  =    -36445.97061862 eV
  free energy    TOTEN  =    -37476.67348261 eV
  free energy    TOTEN  =    -31269.49501946 eV
  free energy    TOTEN  =    -36599.81124818 eV
  free energy    TOTEN  =    -50290.46448839 eV
  free energy    TOTEN  =    -75464.73280329 eV
  free energy    TOTEN  =   -108921.95507254 eV
  free energy    TOTEN  =   -100236.23979020 eV
  free energy    TOTEN  =   -119279.50945865 eV
  free energy    TOTEN  =   -123489.56841725 eV
  free energy    TOTEN  =   -117313.64138222 eV
  free energy    TOTEN  =    -57923.87098952 eV
  free energy    TOTEN  =    -80013.88046309 eV
  free energy    TOTEN  =    -94022.49486397 eV
  free energy    TOTEN  =    -83561.13141342 eV
  free energy    TOTEN  =    -79314.15060362 eV
  free energy    TOTEN  =     17584.00683691 eV
  free energy    TOTEN  =     34600.25250465 eV
  free energy    TOTEN  =     -8480.36214930 eV
  free energy    TOTEN  =     -5345.37903098 eV
  free energy    TOTEN  =    -28521.54315908 eV
  free energy    TOTEN  =    -43513.84362561 eV
  free energy    TOTEN  =    -58726.43416466 eV
  free energy    TOTEN  =    -92729.45742849 eV
  free energy    TOTEN  =    -19732.03384804 eV
  free energy    TOTEN  =    -32282.94544490 eV
  free energy    TOTEN  =    -60925.49871546 eV
  free energy    TOTEN  =   -109028.11307290 eV
  free energy    TOTEN  =    -34501.28002429 eV
  free energy    TOTEN  =    -45529.94231479 eV
  free energy    TOTEN  =    -46789.14053627 eV
  free energy    TOTEN  =    -45965.76236228 eV
  free energy    TOTEN  =    -39486.01567803 eV
  free energy    TOTEN  =    -32105.09912003 eV
  free energy    TOTEN  =    -21713.68330276 eV
  free energy    TOTEN  =    -46719.46336639 eV
  free energy    TOTEN  =    -16142.36303126 eV
  free energy    TOTEN  =    -35066.01543978 eV
  free energy    TOTEN  =    -52655.65474271 eV
  free energy    TOTEN  =    -24243.33046634 eV
  free energy    TOTEN  =    -22616.18034647 eV
  free energy    TOTEN  =    -16399.97625587 eV
  free energy    TOTEN  =   -146709.81100757 eV
  free energy    TOTEN  =    -17713.33309677 eV
  free energy    TOTEN  =    -40319.87901096 eV
  free energy    TOTEN  =   -150898.41812263 eV
  free energy    TOTEN  =    -32163.35804586 eV
  free energy    TOTEN  =    -55897.41764037 eV
  free energy    TOTEN  =   -113212.24789360 eV
  free energy    TOTEN  =    -30110.33716414 eV
  free energy    TOTEN  =   -129426.91210824 eV
  free energy    TOTEN  =    -77056.21588786 eV
  free energy    TOTEN  =    -74995.98011414 eV
  free energy    TOTEN  =    -72263.82877745 eV
  free energy    TOTEN  =    -74654.16280007 eV
  free energy    TOTEN  =    -91255.17300212 eV
  free energy    TOTEN  =     -7440.36451562 eV
  free energy    TOTEN  =    -25141.52530481 eV
  free energy    TOTEN  =    -22389.22167742 eV
  free  energy   TOTEN  =    -22389.221677 eV
  free energy    TOTEN  =    -33972.15094436 eV
  free energy    TOTEN  =    -13496.74830179 eV
  free energy    TOTEN  =     -5672.46025371 eV
  free energy    TOTEN  =     -5521.55767918 eV
  free energy    TOTEN  =     -1000.22551228 eV
  free energy    TOTEN  =     -2126.72772074 eV
  free energy    TOTEN  =     -3139.15831115 eV
  free energy    TOTEN  =     -5452.76892024 eV
  free energy    TOTEN  =     -1308.20870506 eV
  free energy    TOTEN  =      -753.49938920 eV
  free energy    TOTEN  =     -2002.97059421 eV
  free energy    TOTEN  =      -601.23336399 eV
  free energy    TOTEN  =     -1485.90940310 eV
  free energy    TOTEN  =      -834.01932567 eV
  free energy    TOTEN  =     -1397.17227226 eV
  free energy    TOTEN  =      -872.56400073 eV
  free energy    TOTEN  =      -327.23084247 eV
  free energy    TOTEN  =       433.06806795 eV
  free energy    TOTEN  =        92.14699098 eV
  free energy    TOTEN  =      -160.43244370 eV
  free energy    TOTEN  =       763.30842463 eV
  free energy    TOTEN  =      -144.66191026 eV
  free energy    TOTEN  =         3.30223627 eV
  free energy    TOTEN  =       -56.60487105 eV
  free energy    TOTEN  =        12.68521173 eV
  free energy    TOTEN  =      -970.72539298 eV
  free energy    TOTEN  =     -1223.31109998 eV
  free energy    TOTEN  =     -2156.09992354 eV
  free energy    TOTEN  =     -3009.96249140 eV
  free energy    TOTEN  =     -2633.25019754 eV
  free energy    TOTEN  =     -3061.92619150 eV
  free energy    TOTEN  =     -2761.18720117 eV
  free energy    TOTEN  =     -3043.60098772 eV
  free energy    TOTEN  =     -3466.95935152 eV
  free energy    TOTEN  =     -2786.03236398 eV
  free energy    TOTEN  =     -5567.11174658 eV
  free energy    TOTEN  =     -5066.43682904 eV
  free energy    TOTEN  =     -7213.00362756 eV
  free energy    TOTEN  =     -5173.77382366 eV
  free energy    TOTEN  =     -5221.74440492 eV
  free energy    TOTEN  =     -3561.76645405 eV
  free energy    TOTEN  =     -3924.05505322 eV
  free energy    TOTEN  =     -8515.67358844 eV
  free energy    TOTEN  =     -5608.96696042 eV
  free energy    TOTEN  =     -3075.71799628 eV
  free energy    TOTEN  =     -5568.66517526 eV
  free energy    TOTEN  =     -3491.03489852 eV
  free energy    TOTEN  =     -2709.90508221 eV
  free energy    TOTEN  =     -2917.52502232 eV
  free energy    TOTEN  =     -3237.91438610 eV
  free energy    TOTEN  =     -3410.08969989 eV
  free energy    TOTEN  =       118.50923221 eV
  free energy    TOTEN  =      -658.38364023 eV
  free energy    TOTEN  =       537.61600346 eV
  free energy    TOTEN  =       619.43139370 eV
  free energy    TOTEN  =       -18.95985060 eV
  free energy    TOTEN  =      -148.51472777 eV
  free energy    TOTEN  =     -1518.78836312 eV
  free energy    TOTEN  =    -32506.63201196 eV
  free energy    TOTEN  =    -14047.35412927 eV
  free  energy   TOTEN  =    -14047.354129 eV
  free energy    TOTEN  =     -3290.89385440 eV
  free energy    TOTEN  =     -4726.53036556 eV
  free energy    TOTEN  =     -4052.21090007 eV
  free energy    TOTEN  =     -3086.47560383 eV
  free energy    TOTEN  =     -1975.38694601 eV
  free energy    TOTEN  =     -1281.10207172 eV
  free energy    TOTEN  =     -3434.54860951 eV
  free energy    TOTEN  =     -3889.49035408 eV
  free energy    TOTEN  =     -3405.56621420 eV
  free energy    TOTEN  =     -5766.14029034 eV
  free energy    TOTEN  =     -1964.42158500 eV
  free energy    TOTEN  =     -7461.34816926 eV
  free energy    TOTEN  =    -22096.06846693 eV
  free energy    TOTEN  =    -15139.90888562 eV
  free energy    TOTEN  =    -22982.58728208 eV
  free energy    TOTEN  =    -23182.56810246 eV
  free energy    TOTEN  =    -24105.27768967 eV
  free energy    TOTEN  =    -24004.90369629 eV
  free energy    TOTEN  =    -22114.37862868 eV
  free energy    TOTEN  =     -5384.43436703 eV
  free energy    TOTEN  =     16597.51180044 eV
  free energy    TOTEN  =       308.02532860 eV
  free energy    TOTEN  =      -289.37660321 eV
  free energy    TOTEN  =      -797.15129365 eV
  free energy    TOTEN  =     -1165.06810491 eV
  free energy    TOTEN  =     -1549.71576650 eV

这个结果肯定不对，
INCAR设置如下：

ISTART = 0
ICHARG=2
PREC = Accurate
ALGO = F
ENCUT = 500
EDIFF = 1E-5
IBRION = 2
NSW = 100
POTIM = 0.2
ISIF = 0
LREAL = A
ISPIN = 2
ISMEAR = 0
SIGMA = 0.05
LCHARG = F
LWAVE = T
SYMPREC = 1E-4
#LDAU = TRUE
#LDAUTYPE = 2
#LDAUL = -1 2
#LDAUU = 0  5.4
#LDAUJ = 0   0
IDIPOL=3
LDIPOL=.TRUE.

体系的结构参数肯定是对的，以及经过了多次验证
在log.vasp中看到了这样的提示：

  WARNING: DENTET: can't reach specified precision
  Number of Electrons is NELECT =   96.0004255627368
赞一下 回复此楼
7楼2015-08-16 20:43:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur
引用回帖:
7楼: Originally posted by linda19918 at 2015-08-16 20:43:38
结构肯定是对的，计算的能量是这样的

。。。。
   free energy    TOTEN  =    -27013.50693412 eV
  free energy    TOTEN  =    -26816.92597603 eV
  free energy    TOTEN  =    -27176.35064060 eV
   ...

OK，这样就不用多说了，问题一下子就清晰了，截取了几个数据

  free energy    TOTEN  =    -83561.13141342 eV
  free energy    TOTEN  =    -79314.15060362 eV
  free energy    TOTEN  =     17584.00683691 eV
  free energy    TOTEN  =     34600.25250465 eV
  free energy    TOTEN  =     -8480.36214930 eV

总结下：

1、首先电荷数没办法正确
2、能量时而正时而负，数量级变化过快

结论：程序编译有问题，暂与其他因素无关。
建议：可以换台服务器或电脑试试，如果确诊，重新编译一下。
赞一下 回复此楼
不一定挂在论坛，计算问题问题欢迎留言。
8楼2015-08-16 22:51:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kangkangkang

新虫 (小有名气)
引用回帖:
4楼: Originally posted by linda19918 at 2015-08-16 15:00:41
没有用诶，加了这个参数能量还是发散的都不收敛...

降低截断能
发自小木虫手机客户端
回复此楼
9楼2025-03-19 10:45:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 linda19918 的主题更新
91/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭