Gaussian片段初猜报错

snljty

首先顺便请问一下，论坛首页那个谷歌镜像的链接是不是挂了，我这边今天用不了了。

我在想乙烷的碳碳键断裂成一个甲基碳正离子+一个甲基碳负离子的电子结构能否找到，就考虑了片段初猜，设置两个片段分别带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.

完整输出文件在附件中

谢谢大家！

zjxitcc

片段初猜只能用来构造初猜，以供随后算单点；无法直接与结构优化、扫描等任务在一个Link里同时实现。

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

snljty

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

有 ...

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

sobereva

搜索框没挂
暂时不能用就过阵子再试

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

snljty

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

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