24小时热门版块排行榜

返回列表

emma82

金虫 (小有名气)

应助: 0 (幼儿园)
金币: 950.6
帖子: 153
在线: 56.6小时
虫号: 787504
注册: 2009-06-05
专业: 催化化学

[交流] 【求助】如何在dmol中对不同的元素使用不同的基组已有3人参与

请教高手们一个问题，如何在dmol中对不同的元素使用不同的基组，比如说对H，O使用全电子，对铅使用ECP.
是在electronic中直接对core treatment选择还是在input文件中写指令？

回复此楼

1楼 2010-05-20 16:27:32

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

天棚元帅

木虫 (著名写手)

应助: 100 (初中生)
金币: 3141.3
散金: 625
红花: 3
帖子: 1186
在线: 283.9小时
虫号: 609050
注册: 2008-09-21
性别: GG
专业: 理论和计算化学

★
emma82(金币+5): 2010-05-21 09:33:36
gavinliu7390(金币+1):谢谢交流！ 2010-05-21 10:20:01

选择ECP，默认的就是轻的元素全电子，钪（Sc）之后是ECP

赞一下(1人)

回复此楼

A man gets what he earns when he earns it.

2楼2010-05-20 17:15:51

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

ym23

铁杆木虫 (著名写手)

应助: 4 (幼儿园)
贵宾: 0.186
金币: 9769.9
散金: 690
红花: 10
帖子: 1462
在线: 258.8小时
虫号: 837382
注册: 2009-08-31
性别: GG
专业: 凝聚态物性 II ：电子结构

dmol中对不同的元素使用不同的基组

★ ★
zzy870720z(金币+2):xiexie 2010-05-20 21:58:11
emma82(金币+5):谢谢了,很详细! 2010-05-21 09:34:01

◆ DMOL3里面可以采用不同的basis对同一个体系计算，如H2O，H原子采用DND，而O采用DNP，这时候需要关键词为：
Basis  value_keyword
                     min
                     dn
                     dnp
                     dnd
                     mixed
                     extended
                     all
                  user
mixed, extended, all and user. The customer can use their own basis by the ohther atomic programes.
DND=6-31G+polarized orbitals;    DNP=6-31G+double polarized orbitals;
example:
For CH3OH
Basis    mixed
6       dnp
8       dnd
1          extended
如此定义basis后在对CH3OH计算的时候C，O和H采用三种不同的basis函数了，注意keywords里面只列出原子序数即可。

◆ Dmol3中plot脚本编辑主要用于计算结果输出内容的设置，如电子密度分布，差分密度分布(dmol3中称为deformation density), spin density, spin, electrostatic potential, homo, lumo,orbitals and e't'c
File format:
Plot                   value keyword
                        density
                        deformation
                        alpha_density
                        beta_density
                        spin
                        difference
                        potential
                        homo
                        lumo
                        orbital ###
                        fermi ###
                        off
example:
plot density    potential
plot orbital 7  8  9
plot lumo
plot homo
plot fermi    1
这表示在计算结果中输出电子密度，静电势，分子轨道7，8和9，lumo，homo以及Fermi面形状。

赞一下(1人)

回复此楼

博观约取、厚积薄发---壁立千仞、无欲则刚

3楼2010-05-20 18:54:05

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