24小时热门版块排行榜    

查看: 6477  |  回复: 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.55.1.O.54,科目齐全,可伽急 +4 Z7qfCsdIiD1T 2026-11-10 5/250 2026-11-11 22:23 by 9fOTqj074tIg
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 Z7qfCsdIiD1T 2026-11-10 3/150 2026-11-11 22:20 by 9fOTqj074tIg
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 czAiolzR0IJC 2026-11-10 4/200 2026-11-11 22:18 by 9fOTqj074tIg
[基金申请] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 czAiolzR0IJC 2026-11-10 6/300 2026-11-11 22:12 by 9fOTqj074tIg
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 czAiolzR0IJC 2026-11-10 3/150 2026-11-11 22:09 by 9fOTqj074tIg
[考研] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 G2rhSvqvhNzl 2026-11-10 4/200 2026-11-11 21:58 by 9fOTqj074tIg
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 C6CRn66RYYXS 2026-11-10 4/200 2026-11-11 21:36 by 9fOTqj074tIg
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 YHiw1gou1Eda 2026-11-10 4/200 2026-11-11 21:14 by 9fOTqj074tIg
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +5 FChKX8meu1ds 2026-11-09 8/400 2026-11-11 20:49 by 9fOTqj074tIg
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 cqv1zKnzd7hr 2026-11-09 6/300 2026-11-11 20:41 by 9fOTqj074tIg
[考研] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 czAiolzR0IJC 2026-11-10 3/150 2026-11-11 18:31 by TgvkbZZl9tqB
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 czAiolzR0IJC 2026-11-10 4/200 2026-11-11 18:20 by TgvkbZZl9tqB
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 G2rhSvqvhNzl 2026-11-10 5/250 2026-11-11 17:58 by TgvkbZZl9tqB
[基金申请] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 C6CRn66RYYXS 2026-11-10 4/200 2026-11-11 17:53 by TgvkbZZl9tqB
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 FChKX8meu1ds 2026-11-09 6/300 2026-11-11 17:14 by TgvkbZZl9tqB
[考博] 27博士 +3 锦鲤送福至 2026-11-10 3/150 2026-11-11 15:53 by 可淡不可忘
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 0TEMgxdT1eIw 2026-11-09 7/350 2026-11-11 13:31 by cdqGD6ftHmbl
[文学芳草园] 厉娜 +3 张丹枫1991 2026-11-09 5/250 2026-11-09 16:53 by 张丹枫1991
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +5 5sJFDKp2kwuI 2026-11-06 16/800 2026-11-09 11:47 by GwrQVVR6Oenx
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +5 5sJFDKp2kwuI 2026-11-06 9/450 2026-11-08 22:46 by l48hFd8Nwae2
信息提示
请填处理意见