zjxitcc 发表于 2024-12-6 11:22 好的好的谢谢您,感谢两位老师 |
青花 发表于 2024-12-5 22:05 wzkchem5老师已经解释过了,如果你想看具体细节,可以参考我写的这篇《NMR化学位移和耦合常数计算简介》 |
wzkchem5 发表于 2024-12-5 21:25 ! UKS PBE0 6-31G* * xyz 0 3 C 1.72419400 3.83488800 3.25260600 C 1.80773000 2.46836800 3.00652400 C 3.02635200 1.87611400 2.69509100 ..... H 0.02016500 3.91925600 5.76839900 H 0.45154200 5.67215700 5.29456100 * %eprnmr gtensor 1 Nuclei = all H {aiso, adip, aorb, fgrad, rho} end 谢谢您。您看这个是我之前的输入文件,是不是很多东西没考虑到啊, |
青花 发表于 2024-12-5 15:05 如果不同的构象化学不等价,那么需要做构象搜索,确保把相关的构象都搜全了,然后分别跑,取Boltzmann平均。 如果不同的构象化学等价,那么把等价的氢的A值做加权平均就行了。 举例:乙基自由基的三个甲基氢的A值手动平均就行了,不需要额外做计算;对甲基环己基自由基因为甲基可能是平键也可能是直键,这两个构象要分别算,取Boltzmann平均 |
wzkchem5 发表于 2024-12-5 21:25 这个要怎么考虑呢,是要手动跑多个构象取平均吗 |
青花 发表于 2024-12-5 10:59 检查是不是有构象平均效应,在计算中没有考虑到。例如甲基的三个氢,因为甲基快速旋转,实验上只能观察到一个A值,但是计算默认考虑不到转动平均效应,可能给出三个不同的A值,导致峰数比实验更多 |
wzkchem5 发表于 2024-12-4 18:39 老师好,他们实验峰宽只有零点几,按照实验调小了lwpp之后,反而峰变多变杂了,这种情况要怎么办呢 |
wzkchem5 发表于 2024-12-4 20:34 太感谢了,明白了,谢谢您 |
青花 发表于 2024-12-4 12:44 不是。EPR的纵轴是EPR响应信号对磁场的导数(dX/dB,X其实是希腊字母chi),仪器测量的时候通过让B小幅度振荡,测X振荡的幅度,这样推出dX/dB,Modulation Amplitude就是B小幅度振荡的幅度。根据导数的定义知,Modulation Amplitude越小,结果的系统误差越大,但实验噪音的影响也越大(也就是随机误差越大)。Modulation Amplitude越大,lwpp越大,但Modulation Amplitude趋于零的时候lwpp不趋于0,因为自然展宽、振动展宽等等也会导致非零lwpp,Modulation Amplitude只是总展宽里面仪器展宽的那部分,甚至可能都不是仪器展宽的全部。所以lwpp没有任何理由应该接近或等于Modulation Amplitude |
wzkchem5 发表于 2024-12-4 18:39 感谢,再麻烦您一下。仪器上设置的Modulation Amplitude是我们用的lwpp对吗,然后实验谱图出来之后展宽与Modulation Amplitude设定值会有一些出入,我们再以实际为准再来计算,是这样?麻烦了 |
青花 发表于 2024-12-4 11:36 对,调节lwpp直至和实验谱图大致相符,然后明确说明这里lwpp是调出来的,不是算出来的,因此峰宽和实验相符不说明算得准 |
wzkchem5 发表于 2024-12-4 18:04 谢谢老师,成功导出了。我看之前楼里之前的回复,可不可以理解为H数量过大时,允许增大lwpp(比如lwpp=1 → lwpp=2)用于抵消一部分未计算的较小A值的峰。 |
青花 发表于 2024-12-4 10:20 可以,matlab有一个函数叫做dlmwrite,可以把matlab里的任意向量或矩阵导出成csv格式,具体在matlab帮助窗口里搜 |
请问一下老师们,用orca+easyspin模拟EPR曲线,可以把数据导出用origin作图吗 |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2025-8-18 03:11 , Processed in 0.180098 second(s), 25 queries , Gzip On.