24小时热门版块排行榜    

查看: 2879  |  回复: 14
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

peng85722

木虫 (小有名气)

[交流] 具有范德华力的分子晶体结构优化

各位虫子好,
最近在用vasp优化分子晶体,考虑了范德华力修正,优化用的ISIF=3,但是出来的原子位置不对,晶胞参数也与实验值发生很大偏差。我还尝试用实验晶格常数,采用ISIF=2,原子位置也是优化乱了。请问有经验的各位能否给些建议,是需要采用其他计算方法,比如杂化泛函还是用其他的计算软件。
谢谢!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

peng85722

木虫 (小有名气)

引用回帖:
9楼: Originally posted by piaoma at 2016-05-04 10:49:30
估计是你哪设错了。我用vasp和QE都算过分子晶体,加了vdw修正,晶格参数都挺好的

我算的p3ht晶体,
INCAR
SYSTEM = p3ht
ISTART = 0   
ICHARG = 2  
LWAVE = .FALSE.   
LCHARG = .FALSE.
PREC = A
ENCUT = 500   

ALGO = Fast  
NELM = 60

IBRION = 2   ! -1: ions are not moved; 2: conjugate-gradient algorithm
NSW = 160   ! ions are not moved, NSW=0
ISIF = 3
POTIM = 0.2   ! if IBRION=1,2,3 (relaxation), default 0.5
EDIFF = 1.0E-5
EDIFFG = -0.02

#vdW-DF functional
LVDW= .TRUE.
VDW_RADIUS = 30.0
VDW_SCALING = 0.75
VDW_D = 20.0
VDW_C6 =0.14 1.75 5.57
VDW_R0 =1.001 1.452 1.683

ISMEAR = 0   
SIGMA = 0.2
NEDOS = 1000  
LORBIT = 11  

警惕的对称性第一点POSCAR
1.00000000000000     
    16.0000000000000000    0.0000000000000000    0.0000000000000000
     0.4897658570000000    7.7846086851000003    0.0000000000000000
     0.0000000000000000    0.0000000000000000    7.8000001906999996
....
帮忙看一下INCAR 有什么问题?
11楼2016-05-04 14:34:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 15 个回答

jmlv

金虫 (著名写手)

★ ★ ★ ★ ★ ★
小木虫: 金币+0.5, 给个红包,谢谢回帖
peng85722: 金币+5 2016-04-30 16:28:20
我计算最多的是金属晶体,加范德华力修正,计算没问题的。
分子晶体是指哪类晶体?
2楼2016-04-29 14:32:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

klskyqh

铜虫 (小有名气)

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
peng85722: 金币+1 2016-04-30 16:27:20
peng85722: 金币+4 2016-04-30 16:27:32
peng85722: 金币+5 2016-04-30 16:27:43
hse算试试?

发自小木虫IOS客户端
3楼2016-04-30 12:56:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

klskyqh

铜虫 (小有名气)


小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
3楼: Originally posted by klskyqh at 2016-04-30 12:56:20
hse算试试?

还有加u可能也有影响

发自小木虫IOS客户端
4楼2016-04-30 12:58:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见