计算化学公社

标题: Gaussian片段初猜报错 [打印本页]

作者
Author: snljty 时间: 2019-9-21 12:00
标题: Gaussian片段初猜报错
首先顺便请问一下，论坛首页那个谷歌镜像的链接是不是挂了，我这边今天用不了了。

我在想乙烷的碳碳键断裂成一个甲基碳正离子+一个甲基碳负离子的电子结构能否找到，就考虑了片段初猜，设置两个片段分别带1个正电和1个负电，然后降低点群从D3d到C3v，对碳碳键键长进行软扫描。
但是刚一开始算就在link 122报错，加了#p后由于出错部分信息没有变化，就不往上贴了。之前看到介绍l122报错可能是片段的电荷或自选多重度设错了，但是我检查了一下没有问题，是我的这个计算想法不可行么？谢谢！

输入文件部分如下：
# opt=modredundant b3lyp/gen guess=(fragment=2,always)
pop=(always,cm5) empiricaldispersion=gd3bj symmetry=pg=c3v

Title Card Required

0 1 1 1 -1 1
C(Fragment=1)    0.00000000 0.00000000 0.49954362
H(Fragment=1)    0.00000000 1.00905043 1.12707501
H(Fragment=1)    -0.87386330 -0.50452521 1.12707501
H(Fragment=1)    0.87386330 -0.50452521 1.12707501
C(Fragment=2)    -0.00000000 -0.00000000 -0.50045638
H(Fragment=2)    0.87344751 0.50428516 -0.85682724
H(Fragment=2)    -0.00000000 -1.00857031 -0.85682724
H(Fragment=2)    -0.87344751 0.50428516 -0.85682724

基组用的ma-SVP，基组部分没有报错就不贴了。

输出文件最后部分如下：
dumping /fiocom/, unit = 3 NFiles =    1 SizExt = 524288 WInBlk =    512
               defal = F LstWrd =    67072 FType=2 FMxFil=10000

Number             0
Base          40448
End          67072
End1          67072
Wr Pntr       40448
Rd Pntr       40448
Length       26624
Int.Num          1
FileIO operation on non-existent file.
FileIO: IOper= 2 IFilNo(1)=  -584 Len=       24 IPos=          0 Q=       236128864
Error termination in NtrErr:
NtrErr Called from FileIO.

完整输出文件在附件中

谢谢大家！

作者
Author: zjxitcc 时间: 2019-9-21 13:59
片段初猜只能用来构造初猜，以供随后算单点；无法直接与结构优化、扫描等任务在一个Link里同时实现。

有两种解决方案（1）使用高斯自带的Link分步任务。Link0采用片段初猜（C-C距离较远），Link1开始执行柔性扫描，从C-C远距离往近距离扫描；（2）自己写脚本实现，想算哪个结构算哪个，想批量就批量。
PS：其实没必要如此操作点群，直接写nosymm更省事，你那么写，高斯未必真那么做。

作者
Author: snljty 时间: 2019-9-21 14:19

zjxitcc 发表于 2019-9-21 13:59
片段初猜只能用来构造初猜，以供随后算单点；无法直接与结构优化、扫描等任务在一个Link里同时实现。

有 ...

好的，我明白了，谢谢您！我写脚本处理一下

作者
Author: sobereva 时间: 2019-9-21 17:07
搜索框没挂
暂时不能用就过阵子再试

用片段初猜明显得不到你期望的那种电子结构，那需要靠约束性DFT

作者
Author: snljty 时间: 2019-9-21 18:25

sobereva 发表于 2019-9-21 17:07
搜索框没挂
暂时不能用就过阵子再试

谢谢老师！那个尝试波函数确实收敛到D3d的样子了。

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