计算化学公社

标题: 加溶剂的激发态优化求助 [打印本页]

作者
Author: caocao321 时间: 2014-12-15 09:26
标题: 加溶剂的激发态优化求助
我想做C14N4H12分子在溶剂水下的激发态优化，给出的输入文件，
$contrl scftyp=rhf  dfttyp=blyp  tddft=excite runtyp=optimize$end
$pcm  solvnt=water $end
$tddft  nstate=1 mult=1 iroot=1 $end
输出文件总是给出下面信息：
! ! ! ! WARNING ! ! ! !
PCM ITERATIVE PROCEDURE NOT CONVERGED IN                50  ITERATIONS
GAMESS WILL PROCEED DESPITE THIS ...

ITER# 4 ERROR =  1.926E-05
ITER# 5 ERROR =  4.972E-06
ITER# 6 ERROR =  7.385E-07
ITER# 7 ERROR =  1.138E-07
ITER# 8 ERROR =  2.085E-08

! ! ! ! WARNING ! ! ! !
PCM ITERATIVE PROCEDURE NOT CONVERGED IN                50  ITERATIONS
GAMESS WILL PROCEED DESPITE THIS ...

ITER# 9 ERROR =  5.616E-09
ITER#10 ERROR =  9.495E-10
ITER#11 ERROR =  1.831E-10

这样做下来的结果还有意义么？PCM 到底需不需要收敛？

作者
Author: sobereva 时间: 2014-12-15 18:07
原则上应该收敛，否则溶剂化效应描述不准确。如果优化到最后的几步PCM都能收敛那就肯定没关系。
如果虽然没达到收敛限，但是已经收敛到不错精度了，离收敛限不远，那也可以接受。
也可以用Gaussian等程序算算对照验证一下。

作者
Author: caocao321 时间: 2014-12-15 20:55

sobereva 发表于 2014-12-15 18:07
原则上应该收敛，否则溶剂化效应描述不准确。如果优化到最后的几步PCM都能收敛那就肯定没关系。
如果虽然 ...

谢谢Sob
这个是我最后一步优化给出关于溶剂化部分的信息，

      TD-DFT/PCM GRADIENT CODED BY YALI WANG AND HUI LI
               UNIVERSITY OF NEBRASKA-LINCOLN
-----------------------------------------------------------------------

-------------------------
   START Z-VECTOR LOOP
-------------------------
INITIAL ERROR =  1.153E-01, MUST CONVERGE TO 1.000E-10
ITER# 1 ERROR =  2.623E-03
ITER# 2 ERROR =  4.105E-04

! ! ! ! WARNING ! ! ! !
PCM ITERATIVE PROCEDURE NOT CONVERGED IN                50  ITERATIONS
GAMESS WILL PROCEED DESPITE THIS ...

ITER# 3 ERROR =  5.655E-05
ITER# 4 ERROR =  1.042E-05
ITER# 5 ERROR =  1.563E-06
ITER# 6 ERROR =  2.814E-07
ITER# 7 ERROR =  1.542E-08
ITER# 8 ERROR =  2.519E-09
ITER# 9 ERROR =  5.032E-10
ITER#10 ERROR =  1.081E-10
ITER#11 ERROR =  1.696E-11
-------------------------
      Z-VECTOR CONVERGED
-------------------------
MAXIMUM MEMORY DURING TDDFT RESPONSES= 1652301 WORDS.
..... DONE WITH TD-DFT EXCITED STATE RESPONSE AND DENSITY MATRIX .....
我不太清楚这个溶剂化计算的流程，这个算是收敛了么？

作者
Author: sobereva 时间: 2014-12-15 21:26
这个应该是没问题的。这应该是ITER# 2这个Z-vector迭代内部的迭代没收敛，但是后来的Z-vector迭代一直也没再出现这个没收敛提示，并且顺利收敛，说明是可靠的。

作者
Author: caocao321 时间: 2014-12-16 09:12

sobereva 发表于 2014-12-15 21:26
这个应该是没问题的。这应该是ITER# 2这个Z-vector迭代内部的迭代没收敛，但是后来的Z-vector迭代一直也没 ...

好的。谢谢你了的解答。我很想知道这个是Z-vector最后收敛了，这个和PCM迭代有关系么？我对溶剂化一点都不了解。

作者
Author: sobereva 时间: 2014-12-16 17:10
应该是Z-vector迭代的每一步里面都要都要做PCM迭代至收敛

作者
Author: caocao321 时间: 2014-12-16 18:54

sobereva 发表于 2014-12-16 17:10
应该是Z-vector迭代的每一步里面都要都要做PCM迭代至收敛

恩好的。谢谢了

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