计算化学公社

标题: orca单点能得到的gap显著大于MS及文献值 [打印本页]

作者
Author: Apethie 时间: 2022-6-28 16:27
标题: orca单点能得到的gap显著大于MS及文献值
用ORCA对一个量子点模型进行单点能计算，而后通过Multiwfn进行分析，发现得到的HOMO-LUMO gap高达11eV。我用MS（Material Studio）没有加U得到的值为2~3eV。虽然不加U会导致gap变小，但是11和2的差值太大了，文献中使用MS加U计算得到gap也只是3点多eV。
有些疑惑，请问这是为什么，或者给个方向……

作者
Author: sobereva 时间: 2022-6-28 16:44
先把体系特征说清楚，俩程序都用什么具体设置计算的，一切能影响结果的因素全都说清楚，否则根本没回答

仔细把此贴看了
在网上求助计算化学问题的时候必须把问题描述得详细、具体、准确、清楚
http://sobereva.com/620（http://bbs.keinsci.com/thread-25787-1-1.html）

作者
Author: wzkchem5 时间: 2022-6-28 17:13
仔细看http://sobereva.com/543，了解为什么不同计算方法的HOMO-LUMO gap差别很大，以及为什么这一点是完全正常的、没有任何影响

作者
Author: Apethie 时间: 2022-6-28 22:06

wzkchem5 发表于 2022-6-28 17:13
仔细看http://sobereva.com/543，了解为什么不同计算方法的HOMO-LUMO gap差别很大，以及为什么这一点是完全 ...

嗯嗯，谢谢您，我看一下

作者
Author: Apethie 时间: 2022-6-28 22:35

sobereva 发表于 2022-6-28 16:44
先把体系特征说清楚，俩程序都用什么具体设置计算的，一切能影响结果的因素全都说清楚，否则根本没回答

...

抱歉模型是一个1.2nm的ZnO团簇（162个原子），表面悬挂键H原子钝化，且内部一个Zn被替换为Mn原子
orca设置如下
! D4 def2-TZVPP def2/J def2-TZVPP/C RIJCOSX tightSCF noautostart miniprint nopop
%method
Exchange X_PBE
Correlation C_P86
ScalHFX 0.69
ScalDFX 0.31
ScalGGAC 0.4210
ScalLDAC 0.4210
ScalMP2C 1.0
LDAOpt C_VWN5
D3S6 0.5132
D3S8 0.0
D3A1 0.44
D3A2 3.60
end
%mp2
DoSCS true
PS 0.5922
PT 0.0636
end

MS参数设置如下：
comment : CASTEP calculation from Materials Studio
xc_functional : PBE
spin_polarized : true
spin :       5
opt_strategy : Default
page_wvfns :       0
cut_off_energy :    700.000000000000000
grid_scale :       1.500000000000000
fine_grid_scale :       1.500000000000000
finite_basis_corr :       2
finite_basis_npoints :       3
elec_energy_tol : 2.000000000000000e-006
max_scf_cycles :    100
fix_occupancy : false
metals_method : EDFT
num_occ_cycles : 6
nextra_bands :    89
smearing_width :       0.100000000000000
spin_fix :       6
num_dump_cycles : 0
geom_energy_tol : 1.000000000000000e-003
geom_force_tol :       0.100000000000000
geom_stress_tol :       0.200000000000000
geom_disp_tol :       0.010000000000000
geom_max_iter :    200
geom_method : BFGS
fixed_npw : false
geom_modulus_est :    200.000000000000000  GPa

作者
Author: sobereva 时间: 2022-6-29 01:33

Apethie 发表于 2022-6-28 22:35
抱歉模型是一个1.2nm的ZnO团簇（162个原子），表面悬挂键H原子钝化，且内部一个Zn被替换为Mn原子
orca设 ...

连泛函都不统一，对比毫无意义
明明CASTEP用的PBE，干嘛ORCA不用PBE并且记得CASTEP计算时要留出足够的真空区。

作者
Author: Apethie 时间: 2022-6-29 12:12

sobereva 发表于 2022-6-29 01:33
连泛函都不统一，对比毫无意义
明明CASTEP用的PBE，干嘛ORCA不用PBE并且记得CASTEP计算时要留出足够的真 ...

嗯嗯，明白了，谢谢您，我再试试

作者
Author: Apethie 时间: 2022-7-1 11:44

sobereva 发表于 2022-6-29 01:33
连泛函都不统一，对比毫无意义
明明CASTEP用的PBE，干嘛ORCA不用PBE并且记得CASTEP计算时要留出足够的真 ...

请问社长我要用Multiwfn生成上面MS那样的PBE的ORCA输入文件，是不是不能直接生成，要手动再修改。我在Multiwfn里没有看到对应的选项，不知道是不是因为我没搞明白

作者
Author: sobereva 时间: 2022-7-2 05:59

Apethie 发表于 2022-7-1 11:44
请问社长我要用Multiwfn生成上面MS那样的PBE的ORCA输入文件，是不是不能直接生成，要手动再修改。我在Mul ...

量化程序算分子体系用PBE的场合很少，所以没留出相应选项

欢迎光临计算化学公社 (http://bbs.keinsci.com/)