计算化学公社

标题: 碳环激发态优化失败原因 [打印本页]

作者
Author:
summer@123    时间: 2026-3-24 10:01
标题: 碳环激发态优化失败原因
我在做碳环的基态和激发态(S1)优化计算时,基态优化能够正常结束并收敛,对于激发态优化,log文件显示:Gaussian 在读入结构后、正式计算前就被外部中断了,出现这种情况的原因是什么?应该如何解决?感谢各位老师解惑。
输入文件:
%chk=1C32S1.chk
%mem=40GB
%nprocshared=30
#p opt B3LYP/6-31G* TD(nstates=7,root=1) nosymm pop=full

1C32S1

0 1
C                  0.60060000    6.49268312    0.00000000
C                 -0.60060000    6.49268311    0.00000000
C                 -1.92976021    6.22829670    0.00000000
C                 -3.03952431    5.76861737    0.00000000
C                 -4.16633193    5.01570859    0.00000000
C                 -5.01570860    4.16633193    0.00000000
C                 -5.76861738    3.03952431    0.00000000
C                 -6.22829671    1.92976022    0.00000000
C                 -6.49268311    0.60060000    0.00000000
C                 -6.49268310   -0.60060000    0.00000000
C                 -6.22829670   -1.92976021    0.00000000
C                 -5.76861737   -3.03952431    0.00000000
C                 -5.01570860   -4.16633193    0.00000000
C                 -4.16633193   -5.01570860    0.00000000
C                 -3.03952431   -5.76861737    0.00000000
C                 -1.92976021   -6.22829671    0.00000000
C                 -0.60060000   -6.49268310    0.00000000
C                  0.60060000   -6.49268310    0.00000000
C                  1.92976021   -6.22829671    0.00000000
C                  3.03952431   -5.76861738    0.00000000
C                  4.16633193   -5.01570860    0.00000000
C                  5.01570859   -4.16633193    0.00000000
C                  5.76861737   -3.03952431    0.00000000
C                  6.22829670   -1.92976021    0.00000000
C                  6.49268311   -0.60060000    0.00000000
C                  6.49268311    0.60060000    0.00000000
C                  6.22829671    1.92976021    0.00000000
C                  5.76861737    3.03952431    0.00000000
C                  5.01570859    4.16633193    0.00000000
C                  4.16633192    5.01570859    0.00000000
C                  3.03952431    5.76861737    0.00000000
C                  1.92976021    6.22829671    0.00000000
输出文件:

Entering Gaussian System, Link 0=g16
Input=1C32S1.gjf
Output=1C32S1.log
Initial command:
/home/software/g16A03/g16/l1.exe "/tmp/29226.node00/Gau-93442.inp" -scrdir="/tmp/29226.node00/"
Entering Link 1 = /home/software/g16A03/g16/l1.exe PID=     93443.

Copyright (c) 1988,1990,1992,1993,1995,1998,2003,2009,2016,
            Gaussian, Inc.  All Rights Reserved.

This is part of the Gaussian(R) 16 program.  It is based on
the Gaussian(R) 09 system (copyright 2009, Gaussian, Inc.),
the Gaussian(R) 03 system (copyright 2003, Gaussian, Inc.),
the Gaussian(R) 98 system (copyright 1998, Gaussian, Inc.),
the Gaussian(R) 94 system (copyright 1995, Gaussian, Inc.),
the Gaussian 92(TM) system (copyright 1992, Gaussian, Inc.),
the Gaussian 90(TM) system (copyright 1990, Gaussian, Inc.),
the Gaussian 88(TM) system (copyright 1988, Gaussian, Inc.),
the Gaussian 86(TM) system (copyright 1986, Carnegie Mellon
University), and the Gaussian 82(TM) system (copyright 1983,
Carnegie Mellon University). Gaussian is a federally registered
trademark of Gaussian, Inc.

This software contains proprietary and confidential information,
including trade secrets, belonging to Gaussian, Inc.

This software is provided under written license and may be
used, copied, transmitted, or stored only in accord with that
written license.

The following legend is applicable only to US Government
contracts under FAR:

                    RESTRICTED RIGHTS LEGEND

Use, reproduction and disclosure by the US Government is
subject to restrictions as set forth in subparagraphs (a)
and (c) of the Commercial Computer Software - Restricted
Rights clause in FAR 52.227-19.

Gaussian, Inc.
340 Quinnipiac St., Bldg. 40, Wallingford CT 06492


---------------------------------------------------------------
Warning -- This program may not be used in any manner that
competes with the business of Gaussian, Inc. or will provide
assistance to any competitor of Gaussian, Inc.  The licensee
of this program is prohibited from giving any competitor of
Gaussian, Inc. access to this program.  By using this program,
the user acknowledges that Gaussian, Inc. is engaged in the
business of creating and licensing software in the field of
computational chemistry and represents and warrants to the
licensee that it is not a competitor of Gaussian, Inc. and that
it will not use this program in any manner prohibited above.
---------------------------------------------------------------


Cite this work as:
Gaussian 16, Revision A.03,
M. J. Frisch, G. W. Trucks, H. B. Schlegel, G. E. Scuseria,
M. A. Robb, J. R. Cheeseman, G. Scalmani, V. Barone,
G. A. Petersson, H. Nakatsuji, X. Li, M. Caricato, A. V. Marenich,
J. Bloino, B. G. Janesko, R. Gomperts, B. Mennucci, H. P. Hratchian,
J. V. Ortiz, A. F. Izmaylov, J. L. Sonnenberg, D. Williams-Young,
F. Ding, F. Lipparini, F. Egidi, J. Goings, B. Peng, A. Petrone,
T. Henderson, D. Ranasinghe, V. G. Zakrzewski, J. Gao, N. Rega,
G. Zheng, W. Liang, M. Hada, M. Ehara, K. Toyota, R. Fukuda,
J. Hasegawa, M. Ishida, T. Nakajima, Y. Honda, O. Kitao, H. Nakai,
T. Vreven, K. Throssell, J. A. Montgomery, Jr., J. E. Peralta,
F. Ogliaro, M. J. Bearpark, J. J. Heyd, E. N. Brothers, K. N. Kudin,
V. N. Staroverov, T. A. Keith, R. Kobayashi, J. Normand,
K. Raghavachari, A. P. Rendell, J. C. Burant, S. S. Iyengar,
J. Tomasi, M. Cossi, J. M. Millam, M. Klene, C. Adamo, R. Cammi,
J. W. Ochterski, R. L. Martin, K. Morokuma, O. Farkas,
J. B. Foresman, and D. J. Fox, Gaussian, Inc., Wallingford CT, 2016.

******************************************
Gaussian 16:  ES64L-G16RevA.03 25-Dec-2016
                23-Jan-2026
******************************************
%chk=1C32S1.chk
%mem=40GB
%nprocshared=30
Will use up to   30 processors via shared memory.
--------------------------------------------------------
#p opt B3LYP/6-31G* TD(nstates=7,root=1) nosymm pop=full
--------------------------------------------------------
1/18=20,19=15,26=3,38=1/1,3;
2/9=110,12=2,15=1,17=6,18=5,40=1/2;
3/5=1,6=6,7=1,11=9,25=1,30=1,71=1,74=-5/1,2,3,8;
4//1;
5/5=2,38=5/2;
8/6=1,10=2,107=1,108=7/1;
9/8=1,15=2,41=7,42=1,70=2/14;
10/5=4,28=1,31=1/2;
6/7=3/1;
7/12=6,13=5,14=1,30=1/1,2,3,16;
1/18=20,19=15,26=3/3(2);
2/9=110,15=1/2;
99//99;
2/9=110,15=1/2;
3/5=1,6=6,7=1,11=9,25=1,30=1,71=1,74=-5/1,2,3,8;
4/5=5,16=3,69=1/1;
5/5=2,38=5/2;
8/6=1,10=2,107=1,108=7/1;
9/8=1,15=2,41=7,42=1,49=4,70=2/14;
10/5=4,28=1,31=1/2;
7/12=6,13=5,14=1,30=1/1,2,3,16;
1/18=20,19=15,26=3/3(-8);
2/9=110,15=1/2;
6/7=3/1;
99//99;
Leave Link    1 at Fri Jan 23 09:19:20 2026, MaxMem=  5368709120 cpu:               1.0 elap:               0.0
(Enter /home/software/g16A03/g16/l101.exe)
------
1C32S1
------
Symbolic Z-matrix:
Charge =  0 Multiplicity = 1
C                     0.6006    6.49268   0.
C                    -0.6006    6.49268   0.
C                    -1.92976   6.2283    0.
C                    -3.03952   5.76862   0.
C                    -4.16633   5.01571   0.
C                    -5.01571   4.16633   0.
C                    -5.76862   3.03952   0.
C                    -6.2283    1.92976   0.
C                    -6.49268   0.6006    0.
C                    -6.49268  -0.6006    0.
C                    -6.2283   -1.92976   0.
C                    -5.76862  -3.03952   0.
C                    -5.01571  -4.16633   0.
C                    -4.16633  -5.01571   0.
C                    -3.03952  -5.76862   0.
C                    -1.92976  -6.2283    0.
C                    -0.6006   -6.49268   0.
C                     0.6006   -6.49268   0.
C                     1.92976  -6.2283    0.
C                     3.03952  -5.76862   0.
C                     4.16633  -5.01571   0.
C                     5.01571  -4.16633   0.
C                     5.76862  -3.03952   0.
C                     6.2283   -1.92976   0.
C                     6.49268  -0.6006    0.
C                     6.49268   0.6006    0.
C                     6.2283    1.92976   0.
C                     5.76862   3.03952   0.
C                     5.01571   4.16633   0.
C                     4.16633   5.01571   0.
C                     3.03952   5.76862   0.
C                     1.92976   6.2283    0.





作者
Author:
wzkchem5    时间: 2026-3-24 10:14
先试试在同一台机器上能否优化更简单的分子的激发态
作者
Author:
wal    时间: 2026-3-24 12:13
http://bbs.keinsci.com/forum.php ... 3&fromuid=63020
你看看是不是跟这贴一样的问题
作者
Author:
summer@123    时间: 2026-3-24 16:00
wzkchem5 发表于 2026-3-24 10:14
先试试在同一台机器上能否优化更简单的分子的激发态

老师您好,根据您提供的解决方法,我以简单的苯环分子作为测试结构,进行激发态优化,优化命令与碳环分子一致,结构显示计算正常结束。以下是苯环激发态优化的输入文件。
输入文件:
%chk=test.chk
%mem=10GB
%nprocshared=10
#p opt B3LYP/6-31G* TD(nstates=7,root=1) nosymm pop=full

test

0 1
C                  1.21364800    0.70070000   -0.00000000
C                 -0.00000000    1.40140000   -0.00000000
C                 -1.21364800    0.70070000   -0.00000000
C                 -1.21364800   -0.70070000   -0.00000000
C                  0.00000000   -1.40140000   -0.00000000
C                  1.21364800   -0.70070000   -0.00000000
H                  2.14029518    1.23570000   -0.00000000
H                  0.00000000    2.47140000   -0.00000000
H                 -2.14029518    1.23570000   -0.00000000
H                 -2.14029518   -1.23570000   -0.00000000
H                 -0.00000000   -2.47140000   -0.00000000
H                  2.14029518   -1.23570000   -0.00000000


作者
Author:
sobereva    时间: 2026-3-25 02:40
用opt=cartesian
作者
Author:
summer@123    时间: 2026-3-25 08:03
sobereva 发表于 2026-3-25 02:40
用opt=cartesian

问题已解决,感谢sobereva老师!




欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3