如何进行ESR顺磁谱图的模拟？

青花

wzkchem5 发表于 2024-12-4 18:39
对，调节lwpp直至和实验谱图大致相符，然后明确说明这里lwpp是调出来的，不是算出来的，因此峰宽和实验相 ...

老师好，他们实验峰宽只有零点几，按照实验调小了lwpp之后，反而峰变多变杂了，这种情况要怎么办呢

wzkchem5

青花 发表于 2024-12-5 10:59
老师好，他们实验峰宽只有零点几，按照实验调小了lwpp之后，反而峰变多变杂了，这种情况要怎么办呢

检查是不是有构象平均效应，在计算中没有考虑到。例如甲基的三个氢，因为甲基快速旋转，实验上只能观察到一个A值，但是计算默认考虑不到转动平均效应，可能给出三个不同的A值，导致峰数比实验更多

青花

wzkchem5 发表于 2024-12-5 21:25
检查是不是有构象平均效应，在计算中没有考虑到。例如甲基的三个氢，因为甲基快速旋转，实验上只能观察到 ...

这个要怎么考虑呢，是要手动跑多个构象取平均吗

wzkchem5

青花 发表于 2024-12-5 15:05
这个要怎么考虑呢，是要手动跑多个构象取平均吗

如果不同的构象化学不等价，那么需要做构象搜索，确保把相关的构象都搜全了，然后分别跑，取Boltzmann平均。
如果不同的构象化学等价，那么把等价的氢的A值做加权平均就行了。
举例：乙基自由基的三个甲基氢的A值手动平均就行了，不需要额外做计算；对甲基环己基自由基因为甲基可能是平键也可能是直键，这两个构象要分别算，取Boltzmann平均

青花

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

谢谢您。您看这个是我之前的输入文件，是不是很多东西没考虑到啊，

zjxitcc

青花 发表于 2024-12-5 22:05
这个要怎么考虑呢，是要手动跑多个构象取平均吗

wzkchem5老师已经解释过了，如果你想看具体细节，可以参考我写的这篇《NMR化学位移和耦合常数计算简介》

青花

zjxitcc 发表于 2024-12-6 11:22
wzkchem5老师已经解释过了，如果你想看具体细节，可以参考我写的这篇《NMR化学位移和耦合常数计算简介》
...

好的好的谢谢您，感谢两位老师