请问各位，可以利用计算化学的方法来研究三重态分子对其他分子的氧化作用吗？

sobereva

RickyE 发表于 2017-1-17 06:28
sob老师您好，在我的势能面扫描中出现了您所说的单电子都在A 上然后跑到B上的情况，但是分子A的结构出现 ...

没具体信息不好说，结构的跳变往往是越过过渡态，比如在过渡态左边优化会跑到i，过渡态右边优化会跑到j，因此会看到扫描到这个点附近有个突跃。往往这个突跃点被拿来作为过渡态初猜。但也不一定对应于你的情况。
DFT我认为是够用的。在DFT基础上把具体过程、原因能解释通就完了。
你要想把转移过程理想地表达清楚，可考虑先让结构变化连续。比如在跳变前后的结构之间试图插一些结构算算是不是正好电子转移情况也正好介于中间。或者如果怀疑有过渡态走个IRC，这样结构是连续平滑变化的。

由于最近一些天正在讲gromacs培训班，没时间更具体看了。

RickyE

sobereva 发表于 2017-1-17 11:18
没具体信息不好说，结构的跳变往往是越过过渡态，比如在过渡态左边优化会跑到i，过渡态右边优化会跑到j， ...

Sob 老师，我尝试了寻找过渡态，但是过渡态并没有连接PES上的相关点。突跃点前后确实伴随着结构的突变。我的反应物是triplet A + singlet B，突变过程中体系变成了singlet A + triplet B，或者干脆自旋布局离域到两个分子上了，同时伴随着A分子结构的变化。我现在的尝试是直接从singlet A+triplet B 扫描势能面，这样应该能得到连续的势能面，但是也就和实验事实不符合了。现在确实没有什么好的思路了。我把概况了相关工作的poster放到附件了，尽管粗糙但基本也说明了问题。希望能得到您的进一步指教。谢谢Sob 老师了！

sobereva

RickyE 发表于 2017-1-30 04:42
Sob 老师，我尝试了寻找过渡态，但是过渡态并没有连接PES上的相关点。突跃点前后确实伴随着结构的突变。 ...

最近一阵子在香港旅游途中，实在没时间仔细看，你先尝试再搞搞

RickyE

sobereva 发表于 2017-1-31 00:12
最近一阵子在香港旅游途中，实在没时间仔细看，你先尝试再搞搞

哈哈，那我先搞搞，Sob老师旅途愉快！

horald

电子的转移会产生radicals (doublet ), 你说的triplet 交换是什么意思呢？对于你这个反应，参考 triplet ketone 的H－abstraction 反应的计算。

RickyE

horald 发表于 2017-2-1 02:56
电子的转移会产生radicals (doublet ), 你说的triplet 交换是什么意思呢？对于你这个反应，参考 triplet ke ...

反应物是triplet A + singlet B, 但是柔性扫描过程中变成了singlet A + triplet B，同时伴随着能量和结构的跳变。所以可能是我的路径选择问题，或者是势能面的交叉dft无法处理？
horald老师可以分享一下triplet ketone 的H－abstraction 反应的文献吗，谢谢您！

horald

给你找了两篇。
不知道你扫描的coordinate 是什么。尽然影响到反应物的能级 （S0, T1)。理想的情况反应物（amide) 不应有太低的triplet。如果有，反应的机理就不一样了。你再检查下单电子占有轨道的分布。

RickyE

horald 发表于 2017-2-2 02:58
给你找了两篇。
不知道你扫描的coordinate 是什么。尽然影响到反应物的能级 （S0, T1)。理想的情况反应物 ...

谢谢horald老师的文献，我仔细学习一下！
扫描的是三重态下N-H距离逐渐缩短的过程，H 为singlet B上的，N为triplet A 上的，过程中发现体系变成singlet A + triplet B（和实验报道底物相反）。您说到理想反应物不应有太低的triplet，让我想到我优化得到triplet A + singlet B（实验报道底物）和singlet A + triplet B两种体系时，后者能量明显低于前者，请问我是应该直接研究前者，还是假设triplet A 和singlet B 发生能量转移产生singlet A + triplet B，研究后者的H 转移呢？
您说检查单电子占有的轨道分布是检查spin density 和哪些轨道比较重合吗？
再次谢谢您！

horald

1. 假设你用open shell 来计算。对于triplet，你应该看到两个singly occupied orbitals, 如都为alpha 轨道 （自旋相同，triplet）。 你看轨道的电子云密度主要分布在哪个分子上。如果是triplet A + singlet B, 那singly occupied orbitals 的电子云密度应集中在A上。当然A和B 的轨道有degenerate 的情况，那得仔细点。
2. 不知你的triplet A 是不是nitrene, 若果是，去查nitrene 的计算。
3. 比较 energy ［triplet A （T1） － Singlet A （S0， ground state）］ 和 energy ［triplet B （T1） － Singlet B （S0， ground state）］，如果B的triplet （T1）低， energy transfer 可发生在H 转移之前。
4. 你说的triplet B 是失去H后的 吗？
5. 看triplet ketone的 反应应该有帮助。

RickyE

horald 发表于 2017-2-3 05:46
1. 假设你用open shell 来计算。对于triplet，你应该看到两个singly occupied orbitals, 如都为alpha 轨道  ...

1. 谢谢您的指导，我回去查看相应轨道，另外请问这个轨道查询是不是可以和open-shell spin density的结果相对应呢？
2. triplet A 为氰基（如附图所示），所以是不是和ketone的反应较为类似呢？
3. 果然如您所说，B的triplet能量更低，所以energy transfer可能发生在H 转移之前。
4. triplet B 并不是失去H后的，triplet B 中 H 明显仍然属于B分子，所以是先发生了energy transfer吗？
5. 看了您提供的文献，自己也做了一些检索，果然帮助很大，再次谢谢您！

horald

我很怀疑H abstraction 是通过 －CN group，产物你能确定吗。不知你怎么得到 triplet A （dicyanobenzene). 有没有可能A在激发后，产生自由基 （去Cl.), 再发生后面的反应。

RickyE

horald 发表于 2017-2-4 01:37
我很怀疑H abstraction 是通过 －CN group，产物你能确定吗。不知你怎么得到 triplet A （dicyanobenzene). ...

具体产物无法确定，但是单体A受激发在三重态下去Cl.确实是其主要的photochemical fate。您的意思是有可能先发生去Cl.的反应，接着自由基再成为H acceptor吗， 接受H的位点依然是-CN group 吗？另外如果机理确实如此，像这种三重态分子失去自由基的反应，得到的产物多重度是怎样的呢？

horald

实际上A 失去Cl. 是在excited singlet states or triplets 应该确定下。你可以看到很多bond breaking 发生在pi (or n) --> sigma* (Rydberg) 激发态 (singlet)上.
我倾向-CN group 不太可能参与H－abstraction. 失去Cl. 的位置比较可能。

你如果要算在triplet上的反应，整个体系应该保持在Triplet上。但是要注意spin contamination (bond breaking 都要注意下）.

RickyE

horald 发表于 2017-2-7 03:08
实际上A 失去Cl. 是在excited singlet states or triplets 应该确定下。你可以看到很多bond breaking 发生 ...

谢谢horald老师您的建议，我再去多尝试计算一下！