计算化学公社

标题: 用dmol3应该怎么计算这种磷钨酸团簇 [打印本页]

作者
Author: zyj19831206 时间: 2017-1-5 19:10
标题: 用dmol3应该怎么计算这种磷钨酸团簇
初始参数为这些：
# Task parameters
Calculate                   optimize
Opt_energy_convergence       2.0000e-005
Opt_gradient_convergence    4.0000e-003 A
Opt_displacement_convergence  5.0000e-003 A
Opt_iterations             1000
Opt_max_displacement       0.3000 A
Initial_hessian             improved
Symmetry                   off
Max_memory                   6000
File_usage                   smart
Scf_density_convergence    1.000000e-005
Scf_charge_mixing          1.000000e-001
Scf_spin_mixing             2.000000e-001
Scf_diis                   10 pulay
Scf_iterations             1000
Preconditioner             4.000000

# Electronic parameters
Spin_polarization          unrestricted
Charge                      0
Basis                      dnd
basis_version                basfile_v4.4
Pseudopotential             dspp
Functional                   pbe
Aux_density                hexadecapole
Dftd                         TS
Integration_grid             medium
Occupation                   thermal 0.0050
Cutoff_Global                4.5000 angstrom

# Calculated properties
Print_eigval_window          -1.d9
Plot                         homo
Plot                         lumo
Plot                         fermi 10
Plot                         density
Plot                         deformation
Plot                         spindensity
Plot                         potential
Nuclear_EFG                on
Electrostatic_moments       on
Plot                         fukuip

根本不能收敛，是不是初始自旋需要修改，磷钨酸为H2PW12O40.

作者
Author: 柚子不开花 时间: 2017-1-6 15:28
是优化结构不能收敛吗？帖子上有回答增加迭代次数和调整Smearing值能让体系收敛。

作者
Author: zyj19831206 时间: 2017-1-6 18:56

柚子不开花发表于 2017-1-6 15:28
是优化结构不能收敛吗？帖子上有回答增加迭代次数和调整Smearing值能让体系收敛。

这两个都修改了，还是不行，我猜测是因为这个是阴离子需要加负电荷，结果还是不行。

欢迎光临计算化学公社 (http://bbs.keinsci.com/)