24小时热门版块排行榜    

查看: 2265  |  回复: 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个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0703化学调剂 +6 妮妮ninicgb 2026-03-15 9/450 2026-03-16 16:40 by houyaoxu
[考研] 311求调剂 +5 26研0 2026-03-15 5/250 2026-03-16 16:21 by a不易
[考研] 材料与化工求调剂 +3 为学666 2026-03-16 3/150 2026-03-16 15:09 by 加号+
[考研] 285求调剂 +6 ytter 2026-03-12 6/300 2026-03-16 15:05 by njzyff
[考研] 344求调剂 +3 knight344 2026-03-16 3/150 2026-03-16 09:42 by 无际的草原
[考研] 材料与化工 323 英一+数二+物化,一志愿:哈工大 本人本科双一流 +4 自由的_飞翔 2026-03-13 5/250 2026-03-14 19:39 by hmn_wj
[考研] 330求调剂 +3 ?酱给调剂跪了 2026-03-13 3/150 2026-03-14 10:13 by JourneyLucky
[考研] 0703求调剂 +7 jtyq001 2026-03-10 7/350 2026-03-14 01:06 by JourneyLucky
[基金申请] 有必要更换申报口吗 20+3 fannyamoy 2026-03-11 3/150 2026-03-14 00:52 by zhanghaozhu
[考研] 341求调剂 +3 番茄头--- 2026-03-10 3/150 2026-03-13 23:07 by JourneyLucky
[考研] 0856材料与化工301求调剂 +5 奕束光 2026-03-13 5/250 2026-03-13 22:00 by 星空星月
[考研] 290求调剂 +9 ADT 2026-03-11 9/450 2026-03-13 21:55 by JourneyLucky
[考研] 0703化学一志愿211 总分320求调剂 +5 玛卡巴卡啊哈 2026-03-11 5/250 2026-03-13 21:40 by JourneyLucky
[考研] 四川大学085601材料工程专硕 初试294求调剂 +4 祝我们好在冬天 2026-03-11 4/200 2026-03-13 21:39 by peike
[考研] 315求调剂 +9 小羊小羊_ 2026-03-11 10/500 2026-03-13 21:13 by SXNU李老师
[考研] 301求调剂 +6 Liyouyumairs 2026-03-11 6/300 2026-03-13 20:11 by JourneyLucky
[考研] 308求调剂 +3 是Lupa啊 2026-03-12 3/150 2026-03-13 14:30 by 求调剂zz
[考研] 321求调剂(食品/专硕) +3 xc321 2026-03-12 6/300 2026-03-13 08:45 by xc321
[基金申请] 提交后的基金本子,已让学校撤回了,可否换口子提交 +3 dut_pfx 2026-03-10 3/150 2026-03-11 08:38 by kudofaye
[考博] 26申博求助 +3 跳跃饼干 2026-03-10 4/200 2026-03-10 21:15 by Tntcnn
信息提示
请填处理意见