How can I force the Gaussian to carry optimization calculation in each states while using a checkpoint file in calculation. The Gaussian is taking a guess from the previous checkpoint but it is not optimizing and finding the new ground state symmetry for each excited state. The script I use is below.

%chk=benzene_avtz.chk

#P CAM-B3LYP/aug-cc-pVTZ TD(nstates=10)

C6H6

0 1

C 0.0026528488 0.0045977473 0.0

C 1.3879634464 0.0049736534 0.0

C 2.0802932009 1.2048757761 0.0

C 1.3873123579 2.4044019927 0.0

C 0.0020017603 2.4040260866 0.0

C -0.6903279943 1.2041239639 0.0

H -1.7714606137 1.2038305969 0.0

H -0.5388186127 3.3401677164 0.0

H 1.9276246044 3.3408369895 0.0

H 3.1614258204 1.2051691431 0.0

H 1.9287838194 -0.9311679764 0.0

H -0.5376593977 -0.9318372495 0.0

--Link1--

%chk=benzene_avtz.chk

#P geom=allcheck guess=read CAM-B3LYP/aug-cc-pVTZ density=scf out=wfx

benzene_wfx_s0.wfx

--Link1--

%chk=benzene_avtz.chk

#P geom=allcheck guess=read CAM-B3LYP/aug-cc-pVTZ TD(nstates=1) density=(cis=1) out=wfx

benzene_wfx_s1.wfx

--Link1--

%chk=benzene_avtz.chk

#P geom=allcheck guess=read CAM-B3LYP/aug-cc-pVTZ TD(nstates=1,add=1) density=(cis=2) out=wfx

benzene_wfx_s2.wfx

--Link1--

%chk=benzene_avtz.chk

#P geom=allcheck guess=read CAM-B3LYP/aug-cc-pVTZ TD(nstates=2,add=1) density=(cis=3) out=wfx

benzene_wfx_s3.wfx

--Link1--

%chk=benzene_avtz.chk

#P geom=allcheck guess=read CAM-B3LYP/aug-cc-pVTZ TD(nstates=3,add=1) density=(cis=4) out=wfx

benzene_wfx_s4.wfx

--Link1--

%chk=benzene_avtz.chk

#P geom=allcheck guess=read CAM-B3LYP/aug-cc-pVTZ TD(nstates=4,add=1) density=(cis=5) out=wfx

benzene_wfx_s5.wfx

--Link1--

%chk=benzene_avtz.chk

#P geom=allcheck guess=read CAM-B3LYP/aug-cc-pVTZ TD(nstates=5,add=1) density=(cis=6) out=wfx

benzene_wfx_s6.wfx

--Link1--

%chk=benzene_avtz.chk

benzene_wfx_s7.wfx

--Link1--

%chk=benzene_avtz.chk

#P geom=allcheck guess=read CAM-B3LYP/aug-cc-pVTZ TD(nstates=7,add=1) density=(cis=8) out=wfx

benzene_wfx_s8.wfx

--Link1--

%chk=benzene_avtz.chk

#P geom=allcheck guess=read CAM-B3LYP/aug-cc-pVTZ TD(nstates=8,add=1) density=(cis=9) out=wfx

benzene_wfx_s9.wfx

--Link1--

%chk=benzene_avtz.chk

#P geom=allcheck guess=read CAM-B3LYP/aug-cc-pVTZ TD(nstates=9,add=1) density=(cis=10) out=wfx

benzene_wfx_s10.wfx

In the attachment one could see that the guassian is taking a initial guess from checkpoint file. However, within single cycle it stopped the SCF calculation. I am afraid the program is not trying to find the actual optimized state. How can I force the program to find the correct optimized and actual symmetric state? Thanking you in advance.

More Munavvar Husain's questions See All
Similar questions and discussions