24小时热门版块排行榜    

查看: 1241  |  回复: 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的回帖
相关版块跳转 我要订阅楼主 伽罗佤 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 070303 总分349求调剂 +3 LJY9966 2026-03-15 5/250 2026-03-16 14:24 by xwxstudy
[考研] 0703化学调剂 ,六级已过,有科研经历 +5 曦熙兮 2026-03-15 5/250 2026-03-16 14:10 by 可淡不可忘
[考研] 0854控制工程 359求调剂 可跨专业 +3 626776879 2026-03-14 8/400 2026-03-16 13:42 by Dream007008
[考研] 本人考085602 化学工程 专硕 +9 不知道叫什么! 2026-03-15 10/500 2026-03-16 10:04 by houyaoxu
[考研] 化学调剂0703 +7 啊我我的 2026-03-11 7/350 2026-03-15 23:03 by 凌千颂111
[考研] 22408总分284求调剂 +3 InAspic 2026-03-13 3/150 2026-03-15 11:10 by zhq0425
[考研] 289求调剂 +5 步川酷紫123 2026-03-11 5/250 2026-03-15 00:45 by kruisytel
[考研] 中科大材料专硕319求调剂 +3 孟鑫材料 2026-03-13 3/150 2026-03-14 18:10 by houyaoxu
[考研] 【0703化学调剂】-一志愿华中师范大学-六级475 +5 Becho359 2026-03-11 5/250 2026-03-14 11:35 by 哦哦123
[考研] 0805,333求调剂 +3 112253525 2026-03-10 3/150 2026-03-13 23:42 by JourneyLucky
[考研] 341求调剂 +4 番茄头--- 2026-03-10 4/200 2026-03-13 23:12 by JourneyLucky
[考研] 337一志愿华南理工0805材料求调剂 +7 mysdl 2026-03-11 9/450 2026-03-13 22:43 by JourneyLucky
[考研] 0703化学调剂 +4 快乐的香蕉 2026-03-11 4/200 2026-03-13 22:41 by JourneyLucky
[考研] 26调剂/材料科学与工程/总分295/求收留 +9 2026调剂侠 2026-03-12 9/450 2026-03-13 20:46 by 18595523086
[考研] 332求调剂 +3 Zz版 2026-03-13 3/150 2026-03-13 20:36 by 18595523086
[考研] 293求调剂 +3 世界首富 2026-03-11 3/150 2026-03-13 16:27 by JourneyLucky
[考研] 277求调剂 +4 anchor17 2026-03-12 4/200 2026-03-13 11:15 by 白夜悠长
[基金申请] 提交后的基金本子,已让学校撤回了,可否换口子提交 +3 dut_pfx 2026-03-10 3/150 2026-03-11 08:38 by kudofaye
[考研] 调剂 +5 呵唔哦豁 2026-03-10 5/250 2026-03-10 22:00 by 28375m
[考博] 26申博求助 +3 跳跃饼干 2026-03-10 4/200 2026-03-10 21:15 by Tntcnn
信息提示
请填处理意见