查看: 1407  |  回复: 2

topten-gg

银虫 (正式写手)


[交流] 【求助】DMol3怎么使用点电荷模型?

Dmol3怎么用Point_Charge关键词添加点电荷
我用的是MS4.4,在help文件中是这样描述这个关键字的:
Point_ChargePoint_ChargePurposeThe Point_Charge keyword specifies that the calculation is to include point-charges.
SyntaxPoint_Charge value
[table][tr][td]command keyword [/td][td]values [/td][td]default [/td][td]meaning [/td][/tr][tr][td]Point_Charge [/td][td]on/off [/td][td]off [/td][td]Add point charges to the calculation if value is "on". The locations of point charges are read from the runname.pchg file. [/td][/tr][/table]DescriptionPoint charges can be used to simulate the chemical environment due to a large molecule or crystal. The validity of such a calculation depends, of course, on the validity of the point-charge model and the choice of charges used. Point charges can be used in combination with point group symmetry. You need to make sure, however, that the symmetry used in the calculation is a subgroup of the symmetry of both the central molecule and the point charge system.
The point charge file can be in the form of a .car file or in a generic .pchg format. The .pchg format is:
@PCHARGESX Y Z CHARGEX Y Z CHARGE.....$endWhere X, Y and Z are cartesian coordinates of the point charge, CHARGE is its charge and all fields are in floating-point format.

我的input文件是这样写的:
# Task parameters
Calculate                     energy
Symmetry                      on
Max_memory                    2048

# Electronic parameters
Spin_polarization             restricted
Charge                        0
Basis                         dnp
Pseudopotential               none
Functional                    pbe
Aux_density                   hexadecapole
Integration_grid              medium
Occupation                    fermi
Cutoff_Global                 5.5000 angstrom
Scf_density_convergence       5.0000e-002
Scf_charge_mixing             0.2000
Scf_iterations                500
Scf_diis                      off
Point_Charge               on

# Print options
Print                         eigval_last_it

# Calculated properties
Mulliken_analysis             off
Hirshfeld_analysis            off
Esp_fit                       on


点电荷的信息写在一个pointcharge.pchg文件中:
@PCHARGES
   2.605   4.980   0.255     1.265   
   5.175   4.898   1.964     1.485   
   5.133   5.586   3.447    -0.730   
   6.571   4.276   1.865    -0.657   
   4.880   5.990   0.918    -0.688   
......
$end
好像在做计算的时候程序根本就没有调用到pointcharge.pchg这个文件,把
Point_Charge               on
改成
Point_Charge               pointcharge.pchg
好像也没有调用到点电荷的信息,谁知道这个应该怎么设置的么?


有谁知道这里应该怎么设置吗?

谢谢了
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

流浪猫小琪

银虫 (正式写手)


★
小木虫: 金币+0.5, 给个红包,谢谢回帖
我也在用这个关键词,我用的是7.0,命令是Point_Charges on,其他同楼主一样,运行过程是调用了这个命令的,但是显示的是Error: point charges calculation requested but no PCHG file present。我觉得和pchg文件名字有关,导致软件没有读到这个文件。
请问楼主解决这个问题了吗
2楼2018-05-18 18:53:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

流浪猫小琪

银虫 (正式写手)


★
小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
2楼: Originally posted by 流浪猫小琪 at 2018-05-18 18:53:39
我也在用这个关键词,我用的是7.0,命令是Point_Charges on,其他同楼主一样,运行过程是调用了这个命令的,但是显示的是Error: point charges calculation requested but no PCHG file present。我觉得和pchg文件名 ...

自己回复下自己,记得在on后面加上1,表示从pchg文件读取。
3楼2018-05-18 20:04:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 topten-gg 的主题更新
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[有机交流] OpenAI又攻克数学难题了,离AI做有机合成还有多远? +8 wangyl_123 2026-10-10 12/600 2026-10-13 20:40 by Merryya
[硕博家园] 极限拉扯,如何抉择 +9 Isak0320 2026-10-12 13/650 2026-10-13 11:57 by Isak0320
[教师之家] 课题组会帮你通过考核么? +4 好运来~9527 2026-10-10 7/350 2026-10-13 08:41 by xli1984
[基金申请] 入职后第一次管课题组的专利缴费,对账才发现官费比别人多掏了一倍 +6 13108017953 2026-10-08 6/300 2026-10-12 20:07 by w3w3y3
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:34 by XQNyDCa1Tk5n
[找工作] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 22:26 by XQNyDCa1Tk5n
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:26 by XQNyDCa1Tk5n
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:15 by XQNyDCa1Tk5n
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 22:11 by XQNyDCa1Tk5n
[教师之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:06 by XQNyDCa1Tk5n
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 14:25 by VOMexwnDzOFZ
[教师之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 (EPI+-1)(金币-50) +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 14:05 by VOMexwnDzOFZ
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:51 by VOMexwnDzOFZ
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:45 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:54 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:34 by VOMexwnDzOFZ
[考博] 招收博士生一名 +3 mandolin 2026-10-10 3/150 2026-10-11 10:21 by 不是山谷!
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 qR2AdCWTVtkw 2026-10-09 4/200 2026-10-11 03:14 by M0zrYVSAnxpE
[有机交流] 各位大佬,路线设计求助 +6 zapen 2026-10-08 8/400 2026-10-10 15:59 by zapen
[有机交流] 我应该选哪条合成路线让实验成功率高一些? +5 jing816y 2026-10-08 9/450 2026-10-10 09:05 by zapen
信息提示
请填处理意见