用dmol3应该怎么计算这种磷钨酸团簇

zyj19831206

初始参数为这些：
# 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.

zyj19831206

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

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

柚子不开花

是优化结构不能收敛吗？帖子上有回答增加迭代次数和调整Smearing值能让体系收敛。