24小时热门版块排行榜    

查看: 2137  |  回复: 8

piii700

木虫 (小有名气)

[交流] 【求助】大家用DMOL3算的C60分子的LUMO-HOMO带隙是多少啊? 已有3人参与

大家用DMOL3算的C60分子的LUMO-HOMO带隙是多少啊?文献中说实验值是1.7
文献计算结果是1.683。为啥我怎么算,换了各种设置算都是1.8呢啊?
大家给指点一下啊
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

study163

木虫 (正式写手)

天生我材

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
zhang668(金币+1): 多谢交流 2011-01-18 15:53:39
我的计算结果是1.667 eV,你把你的INPUT文件贴出来
我的计算我做主
2楼2011-01-18 15:44:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

piii700

木虫 (小有名气)

# Task parameters
Calculate                     energy
Symmetry                      on
Max_memory                    2000

# Electronic parameters
Spin_polarization             unrestricted
Charge                        0
Basis                         dnp
Pseudopotential               none
Functional                    pbe
Aux_density                   octupole
Integration_grid              fine
Occupation                    thermal 0.0050
Cutoff_Global                 5.5000 angstrom
Scf_density_convergence       1.0000e-005
Scf_charge_mixing             0.2000
Scf_spin_mixing               0.5000
Scf_iterations                500
Scf_diis                      6 pulay

# Print options
Print                         eigval_last_it

# Calculated properties
Plot                          homo
Plot                          lumo
Plot                          fermi 5
Grid                          msbox  3 0.2500 0.2500 0.2500 3.0000
--------------------------------------------------------------------------------------------------------

# Task parameters
Calculate                     energy
Symmetry                      off
Max_memory                    2048

# Electronic parameters
Spin_polarization             restricted
Charge                        0
Basis                         dnp
Pseudopotential               none
Functional                    pwc
Harris                        off
Aux_density                   octupole
Integration_grid              fine
Occupation                    thermal 0.0050
Cutoff_Global                 3.7000 angstrom
Scf_density_convergence       1.0000e-006
Scf_charge_mixing             0.2000
Scf_iterations                50
Scf_diis                      6 pulay

# Print options
Print                         eigval_last_it

# Calculated properties
Plot                          homo
Plot                          lumo
Grid                          msbox  3 0.2500 0.2500 0.2500 3.0000
-----------------------------------------------------------------------------------

# Task parameters
Calculate                     energy
Symmetry                      off
Max_memory                    2048

# Electronic parameters
Spin_polarization             restricted
Charge                        0
Basis                         dnp
Pseudopotential               none
Functional                    pbe
Aux_density                   octupole
Integration_grid              fine
Occupation                    thermal 0.0050
Cutoff_Global                 3.7000 angstrom
Scf_density_convergence       1.0000e-006
Scf_charge_mixing             0.2000
Scf_iterations                50
Scf_diis                      6 pulay

# Print options
Print                         eigval_last_it

# Calculated properties
Plot                          homo
Plot                          lumo
Grid                          msbox  3 0.2500 0.2500 0.2500 3.0000
----------------------------------------------------------------
这3次都是1.8 啊
3楼2011-01-18 15:51:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zzy870720z

荣誉版主 (文坛精英)

优秀版主优秀版主优秀版主优秀版主

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
ellsaking(金币+1): 谢谢回帖交流 2011-01-18 18:47:46
个人的理解是和你的设置精度有关
计算前是否测试了cutoff等参数
还有结构是否进行过优化等等
博学、审问、慎思、明辨、笃学
4楼2011-01-18 18:42:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

study163

木虫 (正式写手)

天生我材

★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
zzy870720z(金币+2): 谢谢建议 2011-01-18 19:40:00
Occupation                    thermal 0.0050

这一项改成fermi,即不加smearing,改好后试一下
我的计算我做主
5楼2011-01-18 19:21:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

piii700

木虫 (小有名气)

谢谢楼上2位啊 我去试试看啊
6楼2011-01-18 19:34:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

piii700

木虫 (小有名气)

还是不行啊
7楼2011-01-18 20:12:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

study163

木虫 (正式写手)

天生我材

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
zzy870720z(金币+1): 谢谢指教 2011-01-19 12:33:31
你最开始计算的时候优化了吗

# Task parameters
Calculate                     energy
Symmetry                      off

# Electronic parameters
Spin_polarization             unrestricted
Charge                        0
Basis                         dnp
Pseudopotential               none
Functional                    gga(p91)
Aux_density                   octupole
Integration_grid              fine
Occupation                    thermal 0.0020
Cutoff_Global                 3.7000 angstrom
Scf_density_convergence       1.0000e-006
Scf_charge_mixing             0.2000
Scf_spin_mixing               0.5000
Scf_iterations                1000
Scf_diis                      6 pulay

# Print options
Print                         eigval_last_it
我的计算我做主
8楼2011-01-19 09:06:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

piii700

木虫 (小有名气)

谢谢 study163   
我再去试试
9楼2011-01-19 09:53:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 piii700 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见