24小时热门版块排行榜    

查看: 5025  |  回复: 6
【悬赏金币】回答本帖问题,作者18366将赠送您 20 个金币

18366

新虫 (小有名气)

[求助] 有什么方法可以得到影响自旋态转变的振动模式?或是分析自旋振动耦合的方法求助 已有1人参与

在Gaussian计算高低自旋态能量和结构优化,有没有什么方法可以计算出影响自旋态转变的关键振动模式?
就像图片里(文章出处:Cite this: J. Chem. Theory Comput. 2022, 18, 3, 1329–1339)文章1332页码里,https://doi.org/10.1021/acs.jctc.1c01184   ,  中的TSH方法“Figure 2a shows the dynamical normal mode activity for the TSH simulation, clearly identifying the dominance of three modes, ν13, ν14, and ν15.”,说明这几个模式影响显著。

不知道有没有什么方法能更简单的得到这个关键振动模式,所以提出问题。

类似的,通过什么方法可以得到自旋振动耦合数据,进而可以分析不同振动模式影响自旋变化的程度。@jiewei
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

paramecium86

版主 (著名写手)

【答案】应助回帖

感谢参与,应助指数 +1
现在阶段Gaussian16并没有surface hopping相关的任何功能。据说25版是要加入相关功能。那等出来之后才能具体看看支持到什么程度。
至于振动自旋耦合部分 也是没办法算的。
2楼2025-05-02 10:05:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

18366

新虫 (小有名气)

引用回帖:
2楼: Originally posted by paramecium86 at 2025-05-02 10:05:49
现在阶段Gaussian16并没有surface hopping相关的任何功能。据说25版是要加入相关功能。那等出来之后才能具体看看支持到什么程度。
至于振动自旋耦合部分 也是没办法算的。

我最近查了下资料,也询问了老师,可以做类似重组能这样的数据。我用的g09,然后分子不太大50个原子,但是振动过程有部分有一点旋转,用了dushin文件运行结果出了一个C1就结束了,就是没运行完,应该是不知道哪里没满足要求,所以结束了。还看到了在g16里计算黄坤的python脚本程序说是也可以算重组能,但是我把g16前面给出的%oldchk=final.chk
#p geom=allcheck freq(readfc,fcht,readfcht)

initial=source=chk final=source=chk spectroscopy=onephotonemission
print=(huangrhys,matrix=JK)

initial.chk
final.chk放在我的g09文件里,gview直接不显示了。因为我不太懂这样的命令,所以我直接粘进g09里试了试,结果就是不行的。到这里,我又没辙了。
发自小木虫手机客户端
3楼2025-05-07 13:27:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

18366

新虫 (小有名气)

引用回帖:
3楼: Originally posted by 18366 at 2025-05-07 13:27:02
我最近查了下资料,也询问了老师,可以做类似重组能这样的数据。我用的g09,然后分子不太大50个原子,但是振动过程有部分有一点旋转,用了dushin文件运行结果出了一个C1就结束了,就是没运行完,应该是不知道哪里没 ...

我用g16打开也是没有显示
发自小木虫手机客户端
4楼2025-05-07 15:09:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

paramecium86

版主 (著名写手)

【答案】应助回帖

引用回帖:
3楼: Originally posted by 18366 at 2025-05-07 00:27:02
我最近查了下资料,也询问了老师,可以做类似重组能这样的数据。我用的g09,然后分子不太大50个原子,但是振动过程有部分有一点旋转,用了dushin文件运行结果出了一个C1就结束了,就是没运行完,应该是不知道哪里没 ...

黄昆引子可以算出来。不过算出来也是在输出文件里。 直接文本编辑器打开输出文件找到对应部分看黄昆引子就行了。这里用不到chk文件和gview。
5楼2025-05-07 16:18:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

18366

新虫 (小有名气)

引用回帖:
5楼: Originally posted by paramecium86 at 2025-05-07 16:18:33
黄昆引子可以算出来。不过算出来也是在输出文件里。 直接文本编辑器打开输出文件找到对应部分看黄昆引子就行了。这里用不到chk文件和gview。...

g09或者是g16计算需不需要%oldchk=final.chk
#p geom=allcheck freq(readfc,fcht,readfcht)

initial=source=chk final=source=chk spectroscopy=onephotonemission
print=(huangrhys,matrix=JK)

initial.chk
final.chk这一段添加在输入文件的前面?我只看到帖子说g16计算黄坤,g09也是可以的吗?
发自小木虫手机客户端
6楼2025-05-07 23:10:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

paramecium86

版主 (著名写手)

引用回帖:
6楼: Originally posted by 18366 at 2025-05-07 10:10:30
g09或者是g16计算需不需要%oldchk=final.chk
#p geom=allcheck freq(readfc,fcht,readfcht)
initial=source=chk final=source=chk spectroscopy=onephotonemission
print=(huangrhys,matrix=JK)
initial.chk
final ...

g09未必可以。  太久远了。  g16是官方给的教程里带有黄昆引子的计算。 oldchk 就是把final.chk里的内容复制一份到新的计算里来读,这样新任务的计算就不会覆盖原始的final.chk  这个计算肯定是要读取final.chk里的内容的。所以是必要的。
7楼2025-05-08 09:15:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 18366 的主题更新
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 科研孤儿太难了 +14 我4大白菜 2026-08-20 15/750 2026-08-20 22:33 by 袁向阳007
[基金申请] 让我中一个面上吧! +8 大萍1987 2026-08-20 10/500 2026-08-20 21:48 by 大萍1987
[基金申请] 时间戳又变了 +12 wuchongjun 2026-08-20 17/850 2026-08-20 20:05 by spzhao
[基金申请] 今天系统多次维护,明天很可能放榜! +10 zju2000 2026-08-16 11/550 2026-08-20 20:02 by cl479861084
[基金申请] 只有每年这种时候来逛逛小木虫 +21 yaoyewhu2008 2026-08-20 22/1100 2026-08-20 17:03 by yaoyewhu2008
[基金申请] 放榜前的不淡定 40+3 snowwithsea 2026-08-19 12/600 2026-08-20 15:56 by lfy8008
[基金申请] 重要消息,中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
[基金申请] filecode,4个jtjc了 +11 ziyangfang 2026-08-19 13/650 2026-08-19 21:28 by aasahr
[教师之家] 为什么余额宝的年化利率越来越低?主要原因有哪些? +5 瞬息宇宙 2026-08-15 5/250 2026-08-19 20:23 by super2002521
[基金申请] 今天放榜没戏了吧 +4 yuleib84 2026-08-19 5/250 2026-08-19 19:52 by hhs666
[基金申请] 今天放榜吗? +14 布布和一二 2026-08-19 15/750 2026-08-19 18:07 by gltch
[基金申请] filecode=后面第一个是大写字母 +10 wangze12014 2026-08-14 12/600 2026-08-19 16:56 by 苦难博士
[基金申请] 朋友圈看到的 +6 wangzilk 2026-08-18 8/400 2026-08-19 10:55 by Haru815
[基金申请] 快农历七夕节了,轻松一下,男人悄悄话,女施主请不要进来。 +6 Tide man 2026-08-14 7/350 2026-08-19 09:56 by ZJTJZ
[基金申请] 什么时候开奖? +6 CrisMessi 2026-08-18 6/300 2026-08-19 08:06 by Equinoxhua
[论文投稿] 投稿咨询 +4 wwm09 2026-08-17 6/300 2026-08-18 15:36 by wwm09
[基金申请] 时间戳又变了8-15 +14 archvillain 2026-08-15 26/1300 2026-08-18 13:37 by phantomgost
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
[精细化工] 招聘 金属平磨液,抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
信息提示
请填处理意见