|
|
各位好,我在计算一个简单化学反应的速率常数的时候,出现了 ground state energy/minimal potential energy has not been initialized 这个报错问题,大家能帮我看看这个输入文件哪里出错了吗?我的输入文件:TemperatureList[K] 298.15 400. 500. 600. 700. 800. 900. 1000. 1100. 1200. 1300. 1400. 1500. 1600. 1700. 1800. 1900. 2000.
PressureList[bar] 1.0
EnergyStepOverTemperature 0.5
ExcessEnergyOverTemperature 30
ModelEnergyLimit[kcal/mol] 400
WellCutoff 10
ChemicalEigenvalueMax 0.2
ChemicalEigenvalueMin 1.e-6
CalculationMethod direct
EigenvalueOutput eigenvalue.out
MicroRateOutput ke.out
MicroEnerMin[kcal/mol] 0.
MicroEnerMax[kcal/mol] 20.
MicroEnerStep[kcal/mol] 0.1
Model
EnergyRelaxation
Exponential
Factor[1/cm] 200
Power 0.85
ExponentCutoff 15
End
CollisionFrequency
LennardJones
Epsilons[1/cm] 240 68.0015
Sigmas[angstrom] 3.7 3.610
Masses[amu] 32 28
End
Well W1 # CH3+NH3
Species #C1H3
RRHO
Geometry[angstrom] 8
C -1.88170200 0.00005000 0.01115800
H -1.89962200 0.00215300 1.09358200
H -1.90152300 -0.93827000 -0.52851500
H -1.89959900 0.93632100 -0.53213000
H 0.72311000 -0.00032300 -0.11531500
N 1.74212500 -0.00005200 -0.11623900
H 2.03670300 -0.81620300 0.41464800
H 2.03626600 0.81638200 0.41445600
Core RigidRotor
SymmetryFactor 1.0
End
Frequencies[1/cm] 18
58.4281 62.0133 81.1457
121.3686 134.5527 182.3571
559.1751 1060.8615 1378.9159
1380.5449 1630.8986 1632.6334
3088.7186 3269.4477 3269.9595
3088.7186 3269.4477 3269.9595
ZeroEnergy[kcal/mol] -6.032
ElectronicLevels[1/cm] 1
0 2
End
End
Bimolecular P1 # CH4+NH2
Fragment C1H4
RRHO
Geometry[angstrom] 5
C 0.00000000 0.00000000 0.00000000
H 0.62967300 0.62967300 0.62967300
H -0.62967300 -0.62967300 0.62967300
H -0.62967300 0.62967300 -0.62967300
H 0.62967300 -0.62967300 -0.62967300
Core RigidRotor
SymmetryFactor 4.0
End
Frequencies[1/cm] 9
1339.5198 1339.5198 1339.5198
1557.6391 1557.6391 3026.2487
3131.5537 3131.5537 3131.5537
ZeroEnergy[kcal/mol] 0
ElectronicLevels[1/cm] 1
0 1
End
Fragment N1H2
RRHO
Geometry[angstrom] 3
N 0.00000000 0.00000000 0.14236700
H 0.00000000 0.80641200 -0.49828500
H 0.00000000 -0.80641200 -0.49828500
Core RigidRotor
SymmetryFactor 2.0
End
Frequencies[1/cm] 3
1507.3578 3349.4140 3440.8127
ZeroEnergy[kcal/mol] 0
ElectronicLevels[1/cm] 1
0 2
End
GroundEnergy[kcal/mol] 0.0
End
Barrier B1 P1 W1 # NH2 + CH4 = CH3+NH3
RRHO
Stoichiometry C1N1H6
SymmetryFactor 1.
End
Geometry[angstrom] 8
C -1.26894800 -0.00002500 0.01374300
H -1.47325700 -0.01893600 1.08146000
H -1.59542800 -0.89910100 -0.50514600
H -1.59240900 0.91833500 -0.47221200
H 0.04005700 -0.00023000 -0.16341200
N 1.31683700 -0.00002400 -0.13531700
H 1.50863200 -0.81042800 0.46217500
H 1.50823400 0.81068100 0.46189700
End
Frequencies[1/cm] 17
45.9112 382.8654 560.8094 783.7523 858.2216 1165.4347
1336.4621 1386.6097 1416.9911 1436.6192 1484.0563 1557.9383
3054.1643 3180.2561 3184.8097 3399.4036 3493.5757
ZeroEnergy[kcal/mol] 28.60
ElectronicLevels[1/cm] 1
0 2
End
End
Tunneling Eckart
ImaginaryFrequency[1/cm] 1650.3711
WellDepth[kcal/mol] 28.60
WellDepth[kcal/mol] 34.63
End
End
|
|