24小时热门版块排行榜    

查看: 1379  |  回复: 4

伽罗佤

银虫 (初入文坛)

[求助] 请问gamess如何做分子或晶体缺陷的频率计算(振动分析)

知道这个问题比较大,但希望热心的虫友们略微指点一二,不甚感激!
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

童生

银虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
伽罗佤: 金币+15, ★★★很有帮助 2013-10-14 20:06:07
$contrl scftyp=rhf runtyp=hessian mplevl=2 $end
$contrl ispher=1 icharg=0 $end
$basis gbasis=acct $end
$SCF DIIS=.T. SOSCF=.F. DIRSCF=.T. FDIFF=.F. NPUNCH=1 $END
$system mwords=2000 memddi=20000 $end
$data
Title: CF3H
C1
C           6.0   0.0000000003   0.0000004598   0.3400540053
H           1.0  -0.0000000001   0.0000000519   1.4256251686
F           9.0  -0.0000000004   1.2513342310  -0.1282022331
F           9.0  -1.0836865140  -0.6256678710  -0.1282024703
F           9.0   1.0836865143  -0.6256678717  -0.1282024704
$end

一个例子。
ANYONEWHOISNOTSHOCKEDBYQUANTUMTHEORYHASNOTUNDERSTOODIT.--NIELSBOHR(1885-1962)
2楼2013-10-14 17:01:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

伽罗佤

银虫 (初入文坛)

引用回帖:
2楼: Originally posted by 童生 at 2013-10-14 17:01:26
$contrl scftyp=rhf runtyp=hessian mplevl=2 $end
$contrl ispher=1 icharg=0 $end
$basis gbasis=acct $end
$SCF DIIS=.T. SOSCF=.F. DIRSCF=.T. FDIFF=.F. NPUNCH=1 $END
$system mwords=2000 memddi=20 ...

非常感谢提供的例子!不知gamess倒底能不能做晶体点缺陷结构的振动分析以及缺陷能级计算
3楼2013-10-14 20:04:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

童生

银虫 (小有名气)

FREQUENCIES IN CM**-1, IR INTENSITIES IN DEBYE**2/AMU-ANGSTROM**2,
     REDUCED MASSES IN AMU.

                          1           2           3           4           5
       FREQUENCY:         4.37        4.04        1.34        0.14        0.72
        SYMMETRY:         A           A           A           A           A
    REDUCED MASS:     12.84737    11.80408    13.59424    13.98000    12.34755
    IR INTENSITY:      0.04138     0.05424     0.01456     0.00010     0.00172

  1   C            X -0.02881784  0.00001063  0.06411300 -0.00016099 -0.00133521
                   Y -0.00000181 -0.03036473  0.00012450 -0.00590945  0.12671828
                   Z -0.00001892 -0.00422795 -0.00002886  0.11935559  0.00455212
  2   H            X -0.16214076 -0.00006507  0.14317033 -0.00015037 -0.00134836
                   Y  0.00006299 -0.18301398  0.00012541 -0.01235102  0.15382532
                   Z -0.00001893 -0.00422821 -0.00002886  0.11935422  0.00455189
  3   F            X -0.03669169  0.00004857 -0.07612483 -0.00027647 -0.00174627
                   Y -0.00002976  0.03547148  0.00012410 -0.00313177  0.11502216
                   Z -0.00009362  0.17171786 -0.00002991  0.12678239 -0.02670064
  4   F            X  0.06138046  0.00003521  0.08308093 -0.00011131 -0.00112040
                   Y -0.05666265  0.03548683 -0.09179071 -0.00322715  0.11466072
                   Z -0.13306212 -0.09227469  0.07888851  0.11565827  0.02016246
  5   F            X  0.06138047  0.00004607  0.08308092 -0.00010892 -0.00112185
                   Y  0.05660316  0.03547765  0.09203894 -0.00303504  0.11538262
                   Z  0.13309898 -0.09212357 -0.07894520  0.11563707  0.02018864

TRANS. SAYVETZ    X  1.12595443  0.00252878  2.62420580 -0.01151988 -0.09315707
                   Y -0.00165397  1.47329007  0.00869403 -0.26183114  8.23134093
                   Z -0.00132456 -0.29590405 -0.00202070  8.35545917  0.31854980
               TOTAL  1.12595642  1.50271396  2.62422098  8.35956854  8.23802924

   ROT. SAYVETZ    X -0.00552200 13.00628415 -0.00007704  0.54874389 -2.31013960
                   Y-11.35954815 -0.00644975  6.73619430  0.00090488 -0.00111833
                   Z  8.81265520 -0.00071367 14.30446400  0.01494802  0.05617098
               TOTAL 14.37714354 13.00628577 15.81119869  0.54894819  2.31082266


用这个语句可以算频率,再算出来的 out 文件中搜“FREQUENCIES”就可以找到如上的结果(当然只是一部分),也可以用 wxMacMolPlt 直接查看 out 文件就和 gaussview 查看 log 文件的频率一样。
至于“缺陷能级计算”,我没做过,所以不清楚。你可以查看手册看是否可以。

» 本帖已获得的红花(最新10朵)

ANYONEWHOISNOTSHOCKEDBYQUANTUMTHEORYHASNOTUNDERSTOODIT.--NIELSBOHR(1885-1962)
4楼2013-10-15 09:15:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

伽罗佤

银虫 (初入文坛)

送红花一朵
引用回帖:
4楼: Originally posted by 童生 at 2013-10-15 09:15:05
FREQUENCIES IN CM**-1, IR INTENSITIES IN DEBYE**2/AMU-ANGSTROM**2,
     REDUCED MASSES IN AMU.

                          1           2           3           4           5
       FREQUENCY:       ...

谢谢耐心的答复,
5楼2013-10-15 12:57:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 伽罗佤 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[找工作] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 Vi50GxzrFcSG 2026-08-07 4/200 2026-08-08 10:39 by 3OOjAIS77qg2
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 HEQlVqMTIA7d 2026-08-07 5/250 2026-08-08 08:59 by 6vVgjDL4CnGu
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 HEQlVqMTIA7d 2026-08-07 4/200 2026-08-08 08:27 by 6vVgjDL4CnGu
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 HEQlVqMTIA7d 2026-08-07 3/150 2026-08-08 03:47 by 6vVgjDL4CnGu
[基金申请] 娱乐 +6 Tide man 2026-08-03 6/300 2026-08-07 22:40 by 铁帽子农民
[基金申请] 化学口download_prp&fileCode的固定段好像这几天一直没变,有变的大神么? +3 Tide man 2026-08-07 4/200 2026-08-07 22:39 by Tide man
[基金申请] 固定端突然变了,今天 +6 archvillain 2026-08-06 10/500 2026-08-07 16:03 by 医学老男孩
[基金申请] 听说今天filecode变了 +24 布布和一二 2026-08-06 47/2350 2026-08-07 16:02 by zhiyanjiang
[基金申请] filecode与中标关系的预测 +3 布布和一二 2026-08-07 3/150 2026-08-07 15:09 by gltch
[基金申请] filecode变化情况 +6 布布和一二 2026-08-07 22/1100 2026-08-07 14:45 by 且听虎啸
[基金申请] 大家散了吧,后缀研究没有意义,别浪费时间了,过好目前的每一天,不要焦虑 +5 Tide man 2026-08-06 7/350 2026-08-07 13:11 by 医学老男孩
[基金申请] 关于豆爷回答的JTJC与%2F数量 +4 yang182083 2026-08-06 6/300 2026-08-07 12:43 by 虫友是什么虫
[基金申请] filecode +14 等待解的谜 2026-08-06 19/950 2026-08-07 12:20 by wlwhappy
[基金申请] filecode +8 布布和一二 2026-08-06 11/550 2026-08-06 20:41 by tangpu318
[基金申请] 影响面上的因素 +8 布布和一二 2026-08-05 11/550 2026-08-06 10:41 by 宝贝虫子
[基金申请] 8月时间戳变的,举个手。玩一下,释放压力 +9 archvillain 2026-08-04 11/550 2026-08-05 20:06 by wlwhappy
[基金申请] 好消息?这个有何含义??? +8 Tide man 2026-08-05 10/500 2026-08-05 16:14 by xmuxiaoyu
[考博] 【2027博士申请】纳米药物递送方向 20+3 13586093586 2026-08-03 4/200 2026-08-05 09:59 by lfy8008
[基金申请] 有没有H口的?有收到消息的吗? +3 超级海虾 2026-08-04 3/150 2026-08-04 17:26 by 学教育滴
[基金申请] 纯娱乐,不喜欢勿喷 +7 Tide man 2026-08-04 10/500 2026-08-04 15:10 by loufangrui
信息提示
请填处理意见