sobereva 发表于 2019-12-20 12:58
靠电子能量比较就够了

好的，谢谢社长！

Frank 发表于 2019-12-20 08:27
谢谢社长回复， 那这样的话，为了保持一致，我如果比较基态和激发态质子转移反应的能垒就只能都用电子能 ...

靠电子能量比较就够了

sobereva 发表于 2019-12-19 18:47
直接优化激发态极小点、过渡态，程序直接就会告诉你激发态电子能量。一般只需要计算激发态电子能量，没太 ...

谢谢社长回复， 那这样的话，为了保持一致，我如果比较基态和激发态质子转移反应的能垒就只能都用电子能量吗？还是说用比如M062X/def2-TZVP 甚至是double hybrid DFT 来直接优化基态和激发态的结构然后用Gaussian输出的总能量(而不去把热力学校正、高精度单点、溶剂化自由能分开去算)？

Frank 发表于 2019-12-19 05:59
请教楼主一个问题，激发态质子转移的反应物、过渡态和产物的能量应该怎么计算？是用基态反应物、过渡态、产 ...

直接优化激发态极小点、过渡态，程序直接就会告诉你激发态电子能量。一般只需要计算激发态电子能量，没太大必要计算激发态自由能

请教楼主一个问题，激发态质子转移的反应物、过渡态和产物的能量应该怎么计算？是用基态反应物、过渡态、产物的总能量（热力学校正+高精度单点+溶剂化自由能）+绝热激发能吗？

sobereva 发表于 2018-2-12 15:55
多电子波函数之间的重叠没有一个可以唯一衡量的方式，不像轨道波函数那样int(psi1(r)psi2(r))就完了
氢 ...

谢谢，我们组做超快光谱的，所以wigner这样的热力学公式基本不适用，不过知道多电子波函数重叠积分很难处理这个信息就非常感谢了！

wswrpd 发表于 2018-2-12 11:16
sob老师，我想追加一下楼主的这个问题。

楼主的第二张图里，S1上有两个local minimum，我取这个两个lo ...

多电子波函数之间的重叠没有一个可以唯一衡量的方式，不像轨道波函数那样int(psi1(r)psi2(r))就完了
氢迁移的隧道效应可直接通过wigner、eckart等公式来计算，取决于氢转移过程的势能面特征（比如wigner需要提供过渡态虚频，eckart需要提供TS和极小点的能量信息）

sobereva 发表于 2017-1-16 07:16
是很明显，但只要用额外方式考虑隧道效应即可，诸如wigner、eckart，或者polyrate支持的更高级的方式。 ...

sob老师，我想追加一下楼主的这个问题。

楼主的第二张图里，S1上有两个local minimum，我取这个两个local minimum的电子波函数为psi_1和psi_2，能不能计算这两个电子波函数的重叠积分？这个问题类似于博文sobereva.com/163，区别是这是分子内的两个电子波函数重叠（由于氢原子转移），这个重叠积分是不是直接关联氢原子隧穿？

另外由于骨架（除氢原子以外的原子们）的波函数可能变化不大，所以算出来的重叠积分是不是也有可能不直接反映氢原子的随穿概率？

这个问题的另一种表述就是，（局域的）电子波函数为psi_1和psi_2之间的耦合(coupling)在ESIPT这个过程里如何理论计算（我们先不考虑溶剂影响）？

0.063eV才6.1kJ/mol，应该能在发荧光前很大比例翻过去。但对于这类问题有没有什么rule of thumb我就不知道了

ggdh 发表于 2017-1-16 09:58
一般认为基态反应的能垒20Kcal /mol就能在室温下反应了。
虽说激发态寿命短，没有太多时间给他反应。但 ...

ggdh老师好，谢谢老师指点，这个20Kcal/mol的标准有没有什么标准，那这样的话。这个0.230eV在激发态上，也可以回到normal了？
我这里的中间，是找的过渡态，这样对吗？
谢谢老师指点。

小范范1989 发表于 2017-1-16 08:37
谢谢sob老师和ggdh老师的指点。我想咨询一下sob老师，如果我就这样计算，能发文章吗？
我的疑惑是别人会 ...

一般认为基态反应的能垒20Kcal /mol就能在室温下反应了。
虽说激发态寿命短，没有太多时间给他反应。但是你这里还有隧道效应
另外你这个扫描是柔性还是刚性扫描？

sobereva 发表于 2017-1-16 07:16
是很明显，但只要用额外方式考虑隧道效应即可，诸如wigner、eckart，或者polyrate支持的更高级的方式。 ...

谢谢sob老师和ggdh老师的指点。我想咨询一下sob老师，如果我就这样计算，能发文章吗？
我的疑惑是别人会不会问，为什么这个0.063eV和0.004eV能够跨过去，而0.230eV和0.382eV不能呢？
一般这个能垒有一个什么标准吗？谢谢老师指点。

ggdh 发表于 2017-1-16 00:05
我有个问题请教：
H的迁移是不是隧道效应比较明显。
在这种情况下。基于BO近似的势能面是不是不再那么可 ...

偶，这样子啊，第一次接触，不是很理解。谢谢你的指点，非常感谢顺带问一句：我的疑惑是别人会不会问，为什么这个0.063eV和0.004eV能够跨过去，而0.230eV和0.382eV不能呢？
谢谢ggdh老师

ggdh 发表于 2017-1-16 00:05
我有个问题请教：
H的迁移是不是隧道效应比较明显。
在这种情况下。基于BO近似的势能面是不是不再那么可 ...

是很明显，但只要用额外方式考虑隧道效应即可，诸如wigner、eckart，或者polyrate支持的更高级的方式。BO势能面还是照用不误，那些考虑隧道效应的方法也都是基于势能面上的信息的，不用BO近似的话也没法定义势能面。