|
|
谢谢回复,我展示我的计算分子给你看看,请帮助看看为什么第五步开始就计算一直很慢了,前4步都是1天完成,第五步开始就一周左右(目前只计算到第6步):
第一步:
%nprocshared=20
%mem=55GB
%chk=2NF1.chk
# Opt Freq b3lyp/6-31G** scrf=(solvent=water,pcm)
2NF1 ground state
2 1
C -8.35184600 -4.37086300 -2.04627200
N -9.69511900 -4.09121800 -1.53104700
C -9.95150800 -2.64930000 -1.52529000
C -8.94078600 -1.94531700 -0.61971000
N -7.57146100 -2.24232700 -1.06487500
C -7.31522600 -3.69885900 -1.15142700
C -6.54551500 -1.40991700 -0.64025100
C -6.71510900 -0.45090700 0.37305700
C -5.68173700 0.42249700 0.74937700
C -4.41526000 0.34187700 0.12484500
C -4.23671300 -0.61190600 -0.89787000
C -5.26542000 -1.43579800 -1.26374900
N -5.87674000 1.37800100 1.76419000
C -4.87448900 2.18817000 2.11448700
C -3.60569300 2.16861800 1.54208900
C -3.30831300 1.21612700 0.49572800
C -7.21037600 1.46406200 2.43010700
C -7.34226500 2.53819200 3.50019900
C -2.60562700 3.11925500 2.02701500
O -1.44188800 3.21253500 1.62592900
O -2.19189700 1.11085400 -0.07938800
F -5.06302900 -2.28251700 -2.31191700
H -8.18044900 -5.45160600 -2.03783700
H -8.21671800 -4.01100500 -3.08183100
H -10.39192400 -4.56475000 -2.09981800
H -10.95861300 -2.46294800 -1.14021900
H -9.88634500 -2.20817400 -2.53585800
H -9.10629100 -0.86703200 -0.67843000
H -9.09504300 -2.26704100 0.42329600
H -7.36910300 -4.13176100 -0.13982600
H -6.31721900 -3.87538500 -1.54604800
H -7.66470500 -0.41092800 0.88169900
H -3.28548200 -0.66665900 -1.41378300
H -5.08693700 2.89754900 2.89910500
H -7.94571900 1.63944600 1.64061300
H -7.40940100 0.48287000 2.86895200
H -8.36393100 2.48236700 3.88565500
H -6.66475600 2.37706400 4.34353000
H -7.19586200 3.54651800 3.10307400
O -3.05097200 3.94060100 2.99553200
H -2.32113000 4.53000700 3.25398100
Ca 0.00000700 2.07663300 0.00000000
C 8.35182500 -4.37093100 2.04622200
N 9.69512000 -4.09123100 1.53108100
C 9.95146000 -2.64930400 1.52536400
C 8.94076800 -1.94533900 0.61973700
N 7.57142800 -2.24240200 1.06481900
C 7.31523400 -3.69894500 1.15132900
C 6.54547800 -1.41001400 0.64016400
C 6.71509500 -0.45096500 -0.37310400
C 5.68172400 0.42243700 -0.74943300
C 4.41522800 0.34178300 -0.12494100
C 4.23665400 -0.61204900 0.89772400
C 5.26535700 -1.43594300 1.26360600
N 5.87675500 1.37800900 -1.76417800
C 4.87452400 2.18821900 -2.11443500
C 3.60571600 2.16864300 -1.54206400
C 3.30829600 1.21606100 -0.49579700
C 7.21040900 1.46413000 -2.43004800
C 7.34235800 2.53846200 -3.49992800
C 2.60568200 3.11935600 -2.02691100
O 1.44194300 3.21263600 -1.62582400
O 2.19187700 1.11077800 0.07930900
F 5.06294300 -2.28269500 2.31174400
H 8.18046500 -5.45168000 2.03775500
H 8.21662400 -4.01109800 3.08177900
H 10.39190700 -4.56475000 2.09988400
H 10.95858100 -2.46291200 1.14035600
H 9.88622400 -2.20819800 2.53593600
H 9.10623400 -0.86705000 0.67848700
H 9.09509600 -2.26703900 -0.42326500
H 7.36918200 -4.13182600 0.13972400
H 6.31721100 -3.87551000 1.54589200
H 7.66471800 -0.41093800 -0.88169000
H 3.28540800 -0.66682800 1.41360600
H 5.08700400 2.89766100 -2.89898900
H 7.94575400 1.63931600 -1.64051100
H 7.40938000 0.48301600 -2.86908800
H 8.36401000 2.48263100 -3.88542100
H 6.66480900 2.37755500 -4.34326900
H 7.19605000 3.54672100 -3.10259800
O 3.05105500 3.94077000 -2.99535700
H 2.32123200 4.53021900 -3.25375900
第二步:
%oldchk=2NFCa1
%chk=2NFCa2
%nprocshared=20
%mem=55GB
# B3LYP/6-31G(d,p) TD=NStates=6 SCRF=(Solvent=water) Geom=Check Guess=Read
NF: linear response vertical excited states
2 1
第三步:
%mem=55GB
%nprocshared=20
%oldchk=2NFCa1
%chk=2NFCa3
# B3LYP/6-31G(d,p) SCRF=(Solvent=Water,Read) Geom=Check Guess=Read
2NFCa: prepare for state-specific non-eq solvation
by saving the solvent reaction field from the ground state
2 1
NonEq=write
--link1--
%mem=55GB
%nprocshared=20
%chk=2NFCa3
# B3LYP/6-31G(d,p) TD(NStates=6,Root=1) SCRF=(Solvent=Water,ExternalIteration,Read) Geom=Check Guess=Read
2NFCa: read non-eq solvation from ground state and compute energy of the first excited with the state-specific method
2 1
NonEq=read
第四步:
%mem=55GB
%nprocshared=20
%oldchk=2NFCa2
%chk=2NFCa4
# B3LYP/6-31G(d,p) TD=(Read,NStates=6,Root=1) SCRF=(Solvent=water)
Geom=Modify Guess=Read Opt=ReadFC
2NFCa: excited state opt
Modify geometry to break Cs symmetry
since first excited state is A"
2 1
13 9 10 11 -160.0
53 52 51 50 -188.0
第五步:
%mem=55GB
%nprocshared=20
%oldchk=2NFCa4
%chk=2NFCa5
# B3LYP/6-31G(d,p) TD=(Read,NStates=6,Root=1) Freq
SCRF=(Solvent=water) Geom=Check Guess=Read
2NFCa excited state freq
2 1
第六步:
%mem=55GB
%nprocshared=20
%oldchk=2NFCa5
%chk=2NFCa6
# B3LYP/6-31G** TD=(read,NStates=6,Root=1) SCRF=(Solvent=Water,ExternalIteration,read) Geom=Check Guess=Read
2NFCa emission state-specific solvation at first excited state optimized geometry
2 1
NonEq=write
|
|