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

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 有什么方法可以得到影响自旋态转变的振动模式？或是分析自旋振动耦合的方法求助
71/1返回列表
查看: 1524  |  回复: 6
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
【悬赏金币】回答本帖问题，作者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
回复此楼

» 猜你喜欢
1楼2025-04-30 22:10:24
已阅   回复此楼   关注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 的主题更新
71/1返回列表
不应助 确定回帖应助 (注意：应助才可能被奖励，但不允许灌水，必须填写15个字符以上)
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭