24小时热门版块排行榜    

查看: 831  |  回复: 5
当前主题已经存档。
【有奖交流】积极回复本帖子,参与交流,就有机会分得作者 weimei212 的 2 个金币
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

weimei212

金虫 (初入文坛)

[交流] 【求助】Dmol中的电荷分析

我用DMOL3算了几个简单的取代苯,算了population analysis 中的MULLIKEN atom charge,我把总和加起来也不是0啊。为什么呢?
Modify中有个charges的calculate,这个计算结果与计算性质出来的原子电荷不一致啊?lable出来的电荷和计算的out文件中的不一样。
哪个准确呢?
请大家指点哈。
一下是Input文件。谢谢。。。
# Task parameters
Calculate                     energy
Symmetry                      on

# Electronic parameters
Spin_polarization             restricted
Charge                        0
Basis                         dnp
Pseudopotential               none
Functional                    pwc
Harris                        off
Aux_density                   octupole
Integration_grid              medium
Occupation                    fermi
Cutoff_Global                 4.3000 angstrom
Scf_density_convergence       1.0000e-005
Scf_charge_mixing             0.2000
Scf_iterations                500
Scf_diis                      6 pulay

# Print options
Print                         eigval_last_it

# Calculated properties
Plot                          homo
Plot                          lumo
Plot                          density
Plot                          deformation
Plot                          potential
Nuclear_EFG                   on
Electrostatic_moments         on
Mulliken_analysis             charge
Hirshfeld_analysis            off
Grid                          msbox  3 0.2500 0.2500 0.2500 3.0000

[ Last edited by weimei212 on 2009-6-9 at 20:59 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dongdong3881

木虫 (正式写手)

我也遇到这样的问题了  关注ing  顶一个
好好自己好好学习
3楼2009-06-09 21:36:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 6 个回答

loovfnd

至尊木虫 (著名写手)

★ ★
weimei212(金币+2,VIP+0):已经加了输入文件,帮忙看下吧。。谢谢~~ 6-9 20:59
最好把输入文件也列出来,这样解决起来容易。
莫强求!
2楼2009-06-09 20:44:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

study163

木虫 (正式写手)

天生我材

★ ★
weimei212(金币+1,VIP+0):谢谢,但是我把symmetry后有个charge选项设成0了啊。还请指点。 6-10 11:25
wuchenwf(金币+1,VIP+0):谢谢 6-10 21:25
MULLIKEN atom charge 指的是电荷转移的情况,相加起来应该是零,如果不是零,你在计算中是不是保持电中性了。
我的计算我做主
4楼2009-06-10 09:18:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

weimei212

金虫 (初入文坛)

引用回帖:
Originally posted by study163 at 2009/03/06 09:18:
MULLIKEN atom charge 指的是电荷转移的情况,相加起来应该是零,如果不是零,你在计算中是不是保持电中性了。

但是我把symmetry后有个charge选项设成0了啊
5楼2009-06-10 11:25:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 3n8v2C8RimXI 2026-09-13 4/200 2026-09-14 15:48 by Aj1rhIDL5ixY
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 3n8v2C8RimXI 2026-09-13 4/200 2026-09-14 15:45 by Aj1rhIDL5ixY
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 3n8v2C8RimXI 2026-09-13 4/200 2026-09-14 15:37 by Aj1rhIDL5ixY
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +5 QUjhNVAcOSff 2026-09-13 6/300 2026-09-14 15:21 by Aj1rhIDL5ixY
[考博] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 QUjhNVAcOSff 2026-09-13 4/200 2026-09-14 15:13 by Aj1rhIDL5ixY
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 QUjhNVAcOSff 2026-09-13 3/150 2026-09-14 15:09 by Aj1rhIDL5ixY
[博后之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 QUjhNVAcOSff 2026-09-13 4/200 2026-09-14 15:01 by Aj1rhIDL5ixY
[公派出国] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 QUjhNVAcOSff 2026-09-13 6/300 2026-09-14 15:00 by Aj1rhIDL5ixY
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 QUjhNVAcOSff 2026-09-13 3/150 2026-09-14 14:57 by Aj1rhIDL5ixY
[考博] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-14 3/150 2026-09-14 12:35 by 40ms4Wlo4umh
[基金申请] 求教各位大神:2026教育部人文社科青年基金项目何时公示呀? +5 云宴山人 2026-09-10 7/350 2026-09-14 12:11 by ecnu2013
[找工作] 浙江师范大学是怎么坑我的 +6 recruit123 2026-09-09 6/300 2026-09-14 12:00 by 萨仁其其格
[考研] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 LwdutQ8HoqWP 2026-09-13 3/150 2026-09-14 11:13 by 23jxep3nCNZb
[博后之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 3n8v2C8RimXI 2026-09-13 3/150 2026-09-14 10:56 by 23jxep3nCNZb
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 QUjhNVAcOSff 2026-09-13 3/150 2026-09-14 10:11 by 23jxep3nCNZb
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 QUjhNVAcOSff 2026-09-13 4/200 2026-09-14 10:01 by 23jxep3nCNZb
[考博] 西北工业大学材料学院电化学传感与催化课题组招收2027年推免研究生、博士生、博士后 +3 马方园 2026-09-10 4/200 2026-09-14 09:29 by Chem张zz
[考博] 27秋季材料博士申请求助 10+3 Cahal1 2026-09-11 4/200 2026-09-14 08:50 by 北京莱茵编辑
[教师之家] 浙江师范大学是怎么坑我的 +10 recruit123 2026-09-09 15/750 2026-09-13 10:45 by seaskyy
[基金申请] 面上有专家说实验设备不是我们单位的 +6 zj_muchong 2026-09-08 8/400 2026-09-11 06:56 by suifengdao
信息提示
请填处理意见