|
|
本帖最后由 〇66 于 2024-11-19 20:45 编辑
各位老师好，我们通过调节F原子的取代位置和个数，以及醛基的个数，合成了一些图上的具有激发态分子内质子转移（ESIPT）特性的分子（为保护组内成果，中间不重要的片段打码处理）。标有绿色笑脸的是可以发生光致变色的分子，而灰色的不可以发生。这些分子都是固态。
我们推测，ESIPT过程中，激发态的能垒是影响这些分子能否发生光致变色的原因，所以我们想对其展开研究，看看能垒是否与其能否发生光致变色有直接的关系。若有关系，我们也进一步也想通过能垒的计算直接预测和判断出还未合成的新分子能否发生光致变色，以便于对未来的实验进行指导。
这是我设想的计算方法：
1. 用高斯直接建模，先对这几个分子在质子转移前（醇的形式）的基态结构和频率进行优化（OPT FREQ M062X/TZVP），得到S0（enol）
2. 对优化好的构型S0（enol）进行激发态计算（td=(nstates=5, root=1) M062X/TZVP），得到S1（enol）
3. 对这几个分子在质子转移后（酮的形式）的构型进行的基态结构和频率进行优化（OPT FREQ M062X/TZVP），得到S0（keto）
4. 对优化好的构型S0(keto)进行激发态计算（opt freq td=(nstates=5, root=1) M062X/TZVP）,得到S1（keto）
5. 寻找S1（enol）和S1(keto）之间的过渡态S1(TS)，（opt=(calcfc,ts,noeigen) freq td=(nstates=5, root=1) M062X/TZVP）
6. 用S1(TS)-S1(enol)的能量与分子能否发生光致变色的实验结果进行对照，看看有无一定的相关关系
目前我有以下几个简单的计算问题想请教各位老师：
1. 这种寻找过渡态的方式有问题吗？这个过程中H应该怎么摆放比较好？
2. 上述第2步，醇形式的分子在做激发态结构优化的时候，是不是会自动优化为酮的形式?所以我没有在这一步加opt和freq，但是这样算出来的醇形式的分子S1的能量就会和S0一样。那么如何计算醇形式分子在激发态S1的能量呢？
3. 对于后续的分子性能预测，除了能垒，是不是还可以关注斯托克斯位移等方面因素？
4. 是否有必要对电子结构的性质进行研究，如静电势、轨道、空间位阻等，以及进行ele-hole分析，研究其对光致变色性能的影响。
5. 对于计算方案的设定以及泛函基组等关键词的使用有没有什么建议和需要注意的地方？
谢谢各位老师解答
|
|
发表于 Post on 2024-11-19 20:40:12
|
只看该作者 Only view this author
|只看大图 View big picture
|倒序浏览 Reverse view
|阅读模式 Reading model
收藏 Add to favorites
楼主 Author