24小时热门版块排行榜    

查看: 1400  |  回复: 4
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

伽罗佤

银虫 (初入文坛)

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

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

» 猜你喜欢

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

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

已阅   回复此楼   关注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的回帖
查看全部 5 个回答

童生

银虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
伽罗佤: 金币+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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 k0dTPqJtl0jt 2026-08-14 5/250 2026-08-15 04:45 by 4wMiSEwB6436
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 k0dTPqJtl0jt 2026-08-14 4/200 2026-08-15 02:21 by 4wMiSEwB6436
[基金申请] 小木虫上这么多卖论文的,真有人买论文么?感觉没必要啊 +11 Tide man 2026-08-10 12/600 2026-08-15 02:12 by home3163
[论文投稿] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 HFw0lei2R37i 2026-08-14 5/250 2026-08-15 00:16 by 4wMiSEwB6436
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 22/1100 2026-08-14 23:45 by Noways
[基金申请] 是这周出结果还是下周出结果? +4 yuleib84 2026-08-11 4/200 2026-08-14 23:05 by lfy8008
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +5 Ldrop2023 2026-08-13 5/250 2026-08-14 18:48 by ssxclkj
[基金申请] 奇怪,两个人的filecode固定段从头到尾一模一样 +8 布布和一二 2026-08-10 11/550 2026-08-14 14:58 by Equinoxhua
[基金申请] 有时候,自然基金真的不能太认真 (我的申报经验) +6 majunge000 2026-08-11 8/400 2026-08-14 14:44 by 医学老男孩
[硕博家园] 读博的好处 +4 lnee 2026-08-11 4/200 2026-08-14 10:20 by ahsoarli
[硕博家园] 一作与独作在应聘高校教师时区别大吗 +3 mbygzh 2026-08-08 4/200 2026-08-13 19:31 by 龙-樱
[基金申请] 重要来源:本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 2019年青年基金涵评意见,大家看看几个A,几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
[基金申请] 综述论文作为代表作会不会影响评审专家的印象分? +11 yufeiwaner 2026-08-09 13/650 2026-08-12 08:17 by yufeiwaner
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[基金申请] 确定了,国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 这样的filecode谁见过 +11 布布和一二 2026-08-08 22/1100 2026-08-10 11:10 by wmfsnow
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
信息提示
请填处理意见