molpro计算效率问题

诸葛壹次心

我有一个算例是用molpro跑的，输入文件见下，C4的体系，用了8个核，目前已经算了15个小时，但是仍然在计算UHF/cc-pvtz的freq一步上，感觉好奇怪啊，简直慢得不可思议，想了很久实在找不到原因。难道是在molpro中用z-matrix的控制某个键长不变会导致方法不是解析的？

1. ***,RO2_230
   
2. memory,1000,M
   
3. print,basis,orbitals
   
4. angstrom
   
5. geometry={                !define the nuclear coordinates
   
6. C;
   
7. H,1,B1;
   
8. H,1,B2,2,A1;
   
9. H,1,B3,3,A2,2,D1;
   
10. C,1,B4,3,A3,2,D2;
    
11. H,5,B5,1,A4,3,D3;
    
12. C,5,B6,1,A5,3,D4;
    
13. H,7,B7,5,A6,1,D5;
    
14. H,7,B8,5,A7,1,D6;
    
15. C,5,B9,1,A8,7,D7;
    
16. H,10,B10,5,A9,1,D8;
    
17. H,10,B11,5,A10,1,D9;
    
18. O,7,B12,5,A11,1,D10;
    
19. H,13,B13,7,A12,5,D11;
    
20. O,10,2.3,5,A13,1,D12;
    
21. O,15,B15,10,A14,5,D13;
    
22. }
    
23. B1=1.09331724
    
24. B2=1.09103427
    
25. B3=1.09270253
    
26. B4=1.53980793
    
27. B5=1.11826342
    
28. B6=1.52925655
    
29. B7=1.09913708
    
30. B8=1.09856152
    
31. B9=1.47120666
    
32. B10=1.08273318
    
33. B11=1.08216781
    
34. B12=1.42672228
    
35. B13=0.96129605
    
36. B15=1.23755496
    
37. A1=108.21528842
    
38. A2=108.77218725
    
39. A3=110.66888452
    
40. A4=107.21925057
    
41. A5=112.66540001
    
42. A6=109.24334854
    
43. A7=109.22114223
    
44. A8=113.3742421
    
45. A9=120.59855293
    
46. A10=121.04870225
    
47. A11=108.2641961
    
48. A12=108.34602049
    
49. A13=97.91817973
    
50. A14=108.49530908
    
51. D1=-116.86692823
    
52. D2=121.17143969
    
53. D3=-176.06339647
    
54. D4=-60.50298586
    
55. D5=-174.73599909
    
56. D6=-56.82232519
    
57. D7=128.87846502
    
58. D8=22.40409801
    
59. D9=-142.47354454
    
60. D10=64.28766706
    
61. D11=-176.65152912
    
62. D12=119.65744006
    
63. D13=-1.35903404
    
64. 
    
65. basis=vtz
    
66. uhf
    
67. optg
    
68. frequencies
    
69. {casscf
    
70. closed,25
    
71. occ,30
    
72. wf,57,1,1
    
73. !start,2200.2
    
74. }
    
75. optg
    
76. frequencies
    
77. rs2
    
78. ---
    

复制代码

诸葛壹次心

molpro计算效率问题