24小时热门版块排行榜

>论坛更新日志 (4067)
>文献求助 (431)
>导师招生 (375)
>虫友互识 (369)
>休闲灌水 (248)
>考博 (168)
>硕博家园 (125)
>招聘信息布告栏 (118)
>博后之家 (96)
>论文投稿 (90)
>基金申请 (71)
>考研 (66)
>公派出国 (65)
>论文道贺祈福 (55)
>找工作 (45)
>教师之家 (40)

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

左边8399

捐助贵宾 (正式写手)

应助: 2 (幼儿园)
金币: 1290.2
散金: 16
红花: 1
帖子: 415
在线: 65.4小时
虫号: 792579
注册: 2009-06-12
性别: MM
专业: 理论和计算化学

[求助] MOLPRO优化计算为什么如此神奇？急啊！

熟悉MOLPRO的朋友们，我想请问一个问题。我测试了水分子用直角坐标计算的例子，输入文件是这样的
***,h2o
memory,500,M
geomtyp=xyz
geometry={
3

O          0.00000000 0.00000000 0.00000000
H          0.00000000 0.00000000 0.96000000
H          0.90509668 0.00000000 -0.32000001
}
basis=aug-cc-pvtz          !use Pople basis set
hf
CCSD                      !closed-shell scf
optg
freq,analytical

输出结果只计算到单点计算完毕不调用结构优化就结束了，如下：
RESULTS
=======

  Reference energy                   -76.056530205072
  CCSD singlet pair energy          -0.182312150042
  CCSD triplet pair energy          -0.084899250723
  CCSD correlation energy             -0.267211404782

!CCSD total energy                -76.323741609854

Timing summary (sec):

STEP                CPU(USER) SYS    CPU(TOT) WALL
Transformation       0.04    0.02    0.06    0.06
CCSD iterations       0.14    0.04    0.18    0.86

Program statistics:

Available memory in ccsd:             499999836
Min. memory needed in ccsd:                61793
Max. memory used in ccsd:                69544
Max. memory used in cckext:             110816 (11 integral passes)

**********************************************************************************************************************************
DATASETS  * FILE NREC LENGTH (MB) RECORD NAMES
            1    19       4.16    500    610    700    900    950    970    1000    129    960    1100
                                       VAR BASINP GEOM SYMINP ZMAT AOBASIS BASIS    P2S ABASIS    S
                                       1400    1410    1200    1210    1080    1600    1650    1700    1380
                                       T       V    H0    H01    AOSYM    SMH MOLCAS OPER    JKOP

            2    3       0.33    700    1000    2100
                                       GEOM    BASIS    RHF

PROGRAMS *       TOTAL    CCSD       HF    INT
CPU TIMES  *       0.65    0.27    0.03    0.14
REAL TIME  *       1.75 SEC
DISK USED  *       67.41 MB
SF USED *       2.80 MB
GA USED *       0.00 MB    (max)    0.00 MB    (current)
**********************************************************************************************************************************

不知道是什么原因，如果用Z-MAT坐标，就可以完整地计算，麻烦有谁能告知是什么原因造成的。因为我有些相对大一点的分子，按对称性转换成内坐标很麻烦，我就是不知道用直角坐标为什么不能算，我加入symm=nosymm也不行

回复此楼