I am trying Geometry Optimization in solvent of the triplet state of a molecule. I am using Q-Chem 6.0.1 and chose PCM as the solvent method. However I get the error message “Q-Chem fatal error occurred in module libgen/rem_setup.C, line 2279: Post-HF gradients are not available with implicit solvent models.”
How should I set up solvent model and parameters to avoid this error? The job as submitted:
$molecule
0 3
C
N 1 B1
C 2 B2 1 A1
C 3 B3 2 A2 1 D1
C 4 B4 3 A3 2 D2
C 5 B5 4 A4 3 D3
C 6 B6 5 A5 4 D4
C 7 B7 6 A6 5 D5
C 8 B8 7 A7 6 D6
C 9 B9 8 A8 7 D7
C 10 B10 9 A9 8 D8
C 6 B11 7 A10 8 D9
H 12 B12 6 A11 7 D10
O 11 B13 12 A12 6 D11
C 14 B14 11 A13 12 D12
C 15 B15 14 A14 11 D13
C 16 B16 15 A15 14 D14
C 17 B17 16 A16 15 D15
C 18 B18 17 A17 16 D16
C 19 B19 18 A18 17 D17
H 20 B20 19 A19 18 D18
H 19 B21 18 A20 17 D19
H 18 B22 17 A21 16 D20
H 17 B23 16 A22 15 D21
H 16 B24 15 A23 14 D22
O 10 B25 11 A24 12 D23
H 9 B26 10 A25 11 D24
H 8 B27 9 A26 10 D25
C 3 B28 4 A27 5 D26
H 29 B29 3 A28 4 D27
H 5 B30 6 A29 7 D28
H 4 B31 5 A30 6 D29
C 2 B32 3 A31 4 D30
H 33 B33 2 A32 3 D31
H 33 B34 2 A33 3 D32
H 33 B35 2 A34 3 D33
H 1 B36 2 A35 3 D34
H 1 B37 2 A36 3 D35
H 1 B38 2 A37 3 D36
B1 1.466943116
B2 1.310854437
B3 1.425190437
B4 1.325760200
B5 1.431340399
B6 1.434384467
B7 1.487331641
B8 1.320021159
B9 1.445488267
B10 1.400706764
B11 1.380350822
B12 1.092994756
B13 1.359191289
B14 1.400704178
B15 1.341638153
B16 1.355071838
B17 1.337770598
B18 1.353011580
B19 1.341514521
B20 1.101204894
B21 1.095027136
B22 1.089127350
B23 1.090118803
B24 1.095181059
B25 1.307187707
B26 1.102303952
B27 1.082591373
B28 1.425388828
B29 1.066101539
B30 1.072263416
B31 1.059657311
B32 1.467092646
B33 1.070952145
B34 1.071402287
B35 1.055881131
B36 1.070551260
B37 1.055178974
B38 1.070045185
A1 120.5803077
A2 122.9843407
A3 120.7396504
A4 123.1962083
A5 116.1390038
A6 125.4770905
A7 124.2727072
A8 127.1102301
A9 121.8117897
A10 124.7943023
A11 117.6818906
A12 111.7085107
A13 118.0458564
A14 117.0446329
A15 119.9699833
A16 120.2317463
A17 118.5072833
A18 121.6928671
A19 120.9402104
A20 119.7426185
A21 121.0292674
A22 119.4340833
A23 118.4500107
A24 119.7454920
A25 114.3434018
A26 118.3119574
A27 115.2164051
A28 117.5476982
A29 117.3053100
A30 120.9163473
A31 120.6460809
A32 110.8801914
A33 112.5195594
A34 107.4597720
A35 111.3521716
A36 107.4534796
A37 111.6570261
D1 -2.113878843
D2 -173.8437533
D3 0.318973048
D4 -7.899020279
D5 -173.6552528
D6 -37.13114696
D7 10.94200123
D8 37.44348948
D9 -0.284269242
D10 -162.8574343
D11 162.8334471
D12 159.4048624
D13 -90.78483687
D14 177.5391546
D15 5.711056249
D16 1.943497778
D17 -7.267827696
D18 -175.1215746
D19 177.7195697
D20 179.7474774
D21 -174.9473100
D22 2.035796158
D23 150.9381382
D24 -145.6045537
D25 179.1323188
D26 6.954174044
D27 171.7983517
D28 175.0504389
D29 -177.5355722
D30 175.9775522
D31 -53.19999934
D32 69.29943227
D33 -171.2536656
D34 59.95544162
D35 178.7094021
D36 -62.55027287
$end
$rem
JOBTYPE OPT
METHOD RIMP2
UNRESTRICTED TRUE
BASIS cc-pVTZ
AUX_BASIS rimp2-cc-pVTZ
MEM_TOTAL 22000
MEM_STATIC 1000
MAX_SCF_CYCLES 800
GEOM_OPT_MAX_CYCLES 800
SOLVENT_METHOD pcm
$end
$solvent
DIELECTRIC 24.6
$end