24小时热门版块排行榜    

查看: 5870  |  回复: 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个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[博后之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 IDs3scOF0tjC 2026-09-28 4/200 2026-10-01 21:59 by Equinoxhua
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +7 IDs3scOF0tjC 2026-09-28 8/400 2026-10-01 19:08 by Equinoxhua
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 CfXuS1rDhLYN 2026-09-28 4/200 2026-10-01 09:23 by Equinoxhua
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 cqwQDCxMcL3I 2026-09-29 4/200 2026-09-30 16:33 by Z9YWQ5EAO3qp
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 cqwQDCxMcL3I 2026-09-29 6/300 2026-09-30 16:26 by Z9YWQ5EAO3qp
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 cqwQDCxMcL3I 2026-09-29 4/200 2026-09-30 13:35 by C3vABopYbZJS
[教师之家] 某top大学教授说“能够在市场中兑现的能力才是真能力”无比同意! +10 zju2000 2026-09-26 13/650 2026-09-30 10:26 by yexuqing
[有机交流] 同一个分子,一条来自文献,一条来自AI——不告诉你答案,你会选哪条? +4 tianxiaoxian 2026-09-28 7/350 2026-09-30 08:34 by 88817753
[博后之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 IDs3scOF0tjC 2026-09-28 4/200 2026-09-29 19:27 by cqwQDCxMcL3I
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 ZYXdhzDAy9ZX 2026-09-28 4/200 2026-09-29 16:47 by yCO1Ll7aHtsw
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 16:03 by etmYJ6d2rquH
[博后之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 15:04 by etmYJ6d2rquH
[找工作] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +6 IDs3scOF0tjC 2026-09-28 6/300 2026-09-29 14:13 by etmYJ6d2rquH
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 IDs3scOF0tjC 2026-09-28 5/250 2026-09-29 13:50 by etmYJ6d2rquH
[教师之家] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 GDBe8tDZqE8z 2026-09-28 4/200 2026-09-29 09:49 by JzYBbHIXWSrW
[找工作] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 06:28 by ZvyPK8n6Nfki
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 GDBe8tDZqE8z 2026-09-28 3/150 2026-09-28 23:12 by b8fCuTHqckEt
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 CfXuS1rDhLYN 2026-09-28 4/200 2026-09-28 22:55 by ez6fGg9abYaj
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 CfXuS1rDhLYN 2026-09-28 3/150 2026-09-28 19:08 by 9lS3ad5oOymn
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 4/200 2026-09-28 18:51 by 9lS3ad5oOymn
信息提示
请填处理意见