应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 这种形式的坐标该怎么得到？
111/2返回列表12下一页
查看: 2194  |  回复: 10
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

nedved33

金虫 (小有名气)

[求助] 这种形式的坐标该怎么得到？

想对过渡态做scan，好像要求的坐标不是笛卡尔坐标，在论坛里的帖子里看到这种形式的坐标，我想问下，类似下面这种坐标该怎么转换得到？可以自动生成吗？
%chk=C6H12O2.chk
# b3lyp/6-311+g(d,p)   scf=maxcycle=1000  scan

C6H6O2

0 1
C              
C                  1            B1
C                  2            B2    1            A1
C                  3            B3    2            A2    1            D1    0
C                  4            B4    3            A3    2            D2    0
C                  5            B5    4            A4    3            D3    0
O                  1            B6    2            A5    3            D4    0
O                  6            B7    5            A6    4            D5    0
H                  1            B8    7            A7    2            D6    0
H                  2            B9    1            A8    7            D7    0
H                  3           B10    2            A9    1            D8    0
H                  4           B11    3           A10    2            D9    0
H                  5           B12    4           A11    3           D10    0
H                  6           B13    5           A12    4           D11    0

   B1             1.47001481
   B2             1.35213504
   B3             1.44805746
   B4             1.34807607
   B5             1.46873534
   B6             1.21337327
   B7             1.21177979
   B8             1.10663924
   B9             1.08542210
   B10            1.08739136
   B11            1.08640585
   B12            1.08633080
   B13            1.11068087
   A1           125.94532223
   A2           127.11728249
   A3           123.17323344
   A4           121.42916035
   A5           122.70659721
   A6           124.00116584
   A7           119.98704577
   A8           114.47240470
   A9           117.57370885
   A10          118.75448281
   A11          121.87705051
   A12          115.15090368
   D1             0.01346178
   D2           179.92161275  S 180  1.0
   D3          -179.98264803
   D4          -179.99936667
   D5           179.98132399
   D6          -179.99685447
   D7             0.01790071
   D8          -179.97661468
   D9            -0.06752315
   D10            0.01014314
   D11           -0.01875418

另外想问一下，经常在一些帖子里看到震荡，有能量震荡，力震荡，位移震荡，请问这个是怎么回事，在输出文件里可以看到吗？具体该怎么看啊，是不是IRC扫描出来的这种图
回复此楼

» 收录本帖的淘帖专辑推荐

Gaussian

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2012-07-09 15:20:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hustzju

木虫 (著名写手)
关注下，帮顶
回复此楼
2楼2012-07-09 16:19:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gmy1990

荣誉版主 (著名写手)

优秀版主优秀版主

【答案】应助回帖

感谢参与，应助指数 +1
坐标用gaussview就可以得到，在保存的时候有选项选择
赞一下 回复此楼
3楼2012-07-09 20:04:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

daixing2011

铜虫 (初入文坛)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
gmy1990: 金币+2 2012-07-09 20:58:23
nedved33: 金币+5, ★★★★★最佳答案 2012-07-11 09:26:10
你是要得到那个震荡的图吗  ？？ 自己没做出过那个图是吧！！

那个是优化就可以看见的， 比如你优化一个结构，得到输出文件1.log
打开高斯view---file---open----找到1.log----注意，把下面方框里的钩打上----ok----这不就打开了吗能看到结构了吧----再点result----optimization----

好了  曲线出来了  够详细吗

这个总能量曲线说的就是 优化过程中  每一步的能量  如果震荡了 会造成优化不出来  就得适当调节了
赞一下 回复此楼
4楼2012-07-09 20:14:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mchen10

金虫 (正式写手)
zmatrix最好是手动建 软件给的多半是没有物理意义的。。
赞一下 回复此楼
更痛苦的是梦没醒路已经不能走了
5楼2012-07-10 05:06:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nedved33

金虫 (小有名气)
引用回帖:
4楼: Originally posted by daixing2011 at 2012-07-09 20:14:27
你是要得到那个震荡的图吗  ？？ 自己没做出过那个图是吧！！

那个是优化就可以看见的， 比如你优化一个结构，得到输出文件1.log
打开高斯view---file---open----找到1.log----注意，把下面方框里的钩打上----o ...

看到了，谢谢哈，但我的能量曲线是刚开始一直呈抛物线下降，后面一大段都是非常平滑的曲线，这样优化应该没问题吧。。
还有出来的是两个图，一个是总能量，还有一个RMS gradient，请问这个参数有什么含义？
赞一下 回复此楼
6楼2012-07-11 09:26:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nedved33

金虫 (小有名气)
引用回帖:
5楼: Originally posted by mchen10 at 2012-07-10 05:06:26
zmatrix最好是手动建 软件给的多半是没有物理意义的。。

请问该怎么建啊，自动的不行吗？不是专门搞计算的，实验室里就我一个人在学这个
赞一下 回复此楼
7楼2012-07-11 09:27:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nedved33

金虫 (小有名气)
引用回帖:
3楼: Originally posted by gmy1990 at 2012-07-09 20:04:46
坐标用gaussview就可以得到，在保存的时候有选项选择

保存出来的格式（也是直角坐标）和上面的不一样哈
类似下面这种是不是得手动一个一个改哈
C                  3            B3    2            A2    1            D1    0
C                  4            B4    3            A3    2            D2    0
C                  5            B5    4            A4    3            D3    0
O                  1            B6    2            A5    3            D4    0
O                  6            B7    5            A6    4            D5    0
H                  1            B8    7            A7    2            D6    0
H                  2            B9    1            A8    7            D7    0
H                  3           B10    2            A9    1            D8    0
H                  4           B11    3           A10    2            D9    0
H                  5           B12    4           A11    3           D10    0
H                  6           B13    5           A12    4           D11    0
赞一下 回复此楼
8楼2012-07-11 09:29:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

weiyin1977

金虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与，应助指数 +1
nedved33: 金币+5, ★★★很有帮助, 3 2012-07-11 18:07:42
用gaussview保存时，把 write cartesians这选项点一下，再存就是internal coordination的形式。
赞一下 回复此楼
9楼2012-07-11 11:34:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nedved33

金虫 (小有名气)
引用回帖:
9楼: Originally posted by weiyin1977 at 2012-07-11 11:34:08
用gaussview保存时，把 write cartesians这选项点一下，再存就是internal coordination的形式。

C              
C           1    1.45231259
C           2    1.45177930    1  107.83199097
C           1    1.44952058    2  108.00481196    3    0.15749167    0
C           4    1.44914021    1  108.11579206    2   -0.03248183    0
C           4    3.14565269    1   87.72104764    2   89.73313775    0
C           6    1.44695371    4   89.61687775    1   -2.28995970    0
C           7    1.45164473    6  108.17532301    4  -87.85908816    0
C           6    1.44898258    4   87.84355427    1 -110.38036368    0
C           8    1.45382778    7  107.90399278    6    0.12519680    0
Fe         10    1.93005963    8   68.14146672    7   58.87656122    0
H           1    1.09338270    4  125.84107418    5  177.69719471    0
H           2    1.09276885    1  125.84371417    4  175.67152182    0
H           5    1.09388148    4  125.48632910    1 -175.83151876    0
H           7    1.09403207    6  125.80921570    4   93.79817582    0
H           8    1.09310117    7  125.83280674    6 -175.84614409    0
H           9    1.09236644    6  125.40835865    4  -93.28015433    0
H          10    1.09239462    8  125.76170978    7 -175.46035241    0
H           6    1.09429780    4   92.27175791    1  123.59630656    0
H           4    1.09391268    1  125.92745526    2 -178.93694302    0
H           3    S             2  125.90936355    1  175.74056549    0
S=1.09026300
按您说的已经做到这一步了，然后是不是要手动把这些数字都弄成变量
比如把最后H和3号碳原子的键长设为S，令S=1.09026300即可
最右边那一串0是什么意思？我看其他人的输入文件里好像没这串0
赞一下 回复此楼
10楼2012-07-11 15:37:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 nedved33 的主题更新
111/2返回列表12下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 基金系统什么内容也没有 30+3 winsaint 2026-08-27 6/300 2026-08-28 01:46 by alongwaytogo
[基金申请] 国自然面上复盘~欢迎讨论 (金币+5) +11 晴天加油 2026-08-26 12/600 2026-08-28 01:22 by li09465290@
[硕博家园] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 XLGfUIdOmEAO 2026-08-27 3/150 2026-08-28 00:49 by 3RAFyYXYSwGS
[教师之家] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 XLGfUIdOmEAO 2026-08-27 4/200 2026-08-27 22:26 by 3RAFyYXYSwGS
[教师之家] 导师吐槽：我怎么摊上了这么个极品研究生！ +7 苏东坡二世 2026-08-23 7/350 2026-08-27 18:14 by 瞬息宇宙
[基金申请] 申请删除本帖 +6 lyz123lyz 2026-08-27 7/350 2026-08-27 17:31 by 宁静致远sy
[基金申请] 基金不中，共勉 +10 eulota 2026-08-26 10/500 2026-08-27 17:23 by lqllinqiaoli
[基金申请] 基金未中，这种答复是模板吗？ +5 zhaosm1982 2026-08-27 6/300 2026-08-27 16:00 by lfy8008
[文学芳草园] 梦想 +3 myrtle 2026-08-26 3/150 2026-08-27 10:01 by angelyueyi
[教师之家] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 LIbGuocjEEYw 2026-08-26 4/200 2026-08-27 04:33 by Ie9AyIAvGbvs
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 LIbGuocjEEYw 2026-08-26 4/200 2026-08-27 02:01 by Ie9AyIAvGbvs
[基金申请] 我不理解！ +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +7 process2012 2026-08-23 10/500 2026-08-26 19:23 by hmhminy
[基金申请] 范进中举一文的中心思想 +9 炎黄贵胄 2026-08-22 10/500 2026-08-26 15:40 by semaglutide
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件，有关系吗？ +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 国合现在查不到了吗？ +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
[基金申请] 在坚冰还盖着北海的时候，我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
信息提示
关闭
请填处理意见
关闭