24小时热门版块排行榜    

查看: 591  |  回复: 7
当前主题已经存档。

meijuancao

铁虫 (初入文坛)

[交流] 【求助】Dmol3问题

我算的231个原子的钯团簇,dmol3出现如下错误信息:
Problem with geometry-independent mesh, switching to adaptive mesh

Integration points and checksum:    1075509    4158.205293  302

Error: poor numerical integration accuracy, use finer integration
但我已经用了fine的integration grid了

请热心人指教。另外,dmol3能不能算这样大的过渡金属团簇。
input 参数如下:
INPUT_DMOL keywords (for archive):
______________________________________________________________________>8

                                                                                <--
# Task parameters                                                               <--
Calculate                     optimize                                          <--
Opt_energy_convergence        1.0000e-005                                       <--
Opt_gradient_convergence      2.0000e-003 A                                     <--
Opt_displacement_convergence  5.0000e-003 A                                     <--
Opt_iterations                50                                                <--
Opt_max_displacement          0.3000 A                                          <--
Symmetry                      on                                                <--
                                                                                <--
# Electronic parameters                                                         <--
Spin_polarization             restricted                                        <--
Charge                        0                                                 <--
Basis                         dnp                                               <--
Pseudopotential               ecp                                               <--
Functional                    vwn                                               <--
Harris                        on                                                <--
Aux_density                   monopole                                          <--
Integration_grid              fine                                              <--
Occupation                    thermal 0.0050                                    <--
Cutoff_Global                 4.5000 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                                              <--
Mulliken_analysis             charge                                            <--
Hirshfeld_analysis            charge                                            <--
Frequency_analysis            on                                                <--
Grid                          msbox  3 0.2500 0.2500 0.2500 3.0000              <--
______________________________________________________________________>8


Density functional:
Vosko Wilk Nusair local correlation


ECP takes precedence over using symmetrized basis




Calculation is Spin_restricted

Molecule has been put into center of mass coordinate system

Symmetry group of the molecule: oh

[ Last edited by freshgirl on 2009-6-25 at 19:14 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhengjinde

木虫 (正式写手)

似乎你的自洽场轮数设置太少了
Opt_iterations                50         
Scf_iterations                50
231个金属原子我想50轮之内是无法结束计算的。建议改大些再看看

[ Last edited by zhengjinde on 2008-10-14 at 16:44 ]
2楼2008-10-14 16:39:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

meijuancao

铁虫 (初入文坛)

我试过了,还是不行,根本就没有迭代就退出来了
我现在越来越怀疑dmol3能不能算动这个体系了
3楼2008-10-14 17:09:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

solarman

新虫 (正式写手)

多半是你机器配置不够,231个过渡金属,terribly large!
4楼2008-10-15 00:20:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

canbill

银虫 (小有名气)

确实太大了,估计在软件的迭代上限内无法收敛。
5楼2008-10-15 10:01:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhenben

铜虫 (初入文坛)

估计算不了吧……
6楼2008-10-15 10:44:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

donkeypku

至尊木虫 (著名写手)

两院院士


spur(金币+1,VIP+0):3Q!
应该是体系太大,dmol这个程序变得很有问题,比如对内存利用很不好,对大体系会莫名其妙的无法运算
博观约取,厚积薄发
7楼2008-10-15 23:50:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

沙风

铁杆木虫 (正式写手)

有为青年

把Integration_grid   改成       xfine
8楼2008-10-16 10:56:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 meijuancao 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见