### modified anatolia 

修改anatolia的程序，进行全谱线分析。这个程序分为三部分：


* 第一步计算哈密顿方程的特征值和特征函数求谱峰的峰值高度和出峰位置；
* 第二步利用lorenztion函数计算谱图
* 第三步用bobyqa对化学位移和耦合常数，谱峰高度和线性宽度LB进行拟合。

ANATOLIA: NMR software for spectral analysis of total lineshape

https://analyticalsciencejournals.onlinelibrary.wiley.com/doi/full/10.1002/mrc.4689

### 给定化学位移和耦合常数计算谱线

为啥计算化学位移和薛定谔方程的关系，然后又为啥薛定谔方程与哈密顿矩阵的关系？


如何构建哈密顿矩阵，求解矩阵的特征值和特征函数，如何通过特征函数求解谱峰高度，还有如何筛选出有效的跃迁频率。学会这些就会明白谱线的出峰情况，包括一阶谱峰和高阶谱峰的分析。

参考文献：Harald Gunther, NMR Spectroscopy basic principles, concepts and applications in chemistry, Chapter 6


### determining the Number of signals in NMR Spectrometer

#### 氢谱的核磁信号个数

主要就是判断化学等价（Chemical shift equivalence）

- reflection through a mirror plane in the molecule
- rotation of the molecule itself
- fast bond rotation （困难，C-C的旋转）

后面才判断两个氢同接在一个碳上面的情况，分为Homotopic，enantiotopic, diastereotopic。


有相关的文章探讨了分子的对称性，[Simple Algorithms for Determining the Molecular Symmetry](https://pubs.acs.org/doi/10.1021/ci990322q), [Identification of Symmetries in Molecules and Complexes](https://pubs.acs.org/doi/full/10.1021/ci049966a),[New computer program to calculate the symmetry of molecules](https://www.degruyter.com/document/doi/10.2478/BF02475193/html)这篇文章提供了MATLAB算法。

接下来我用SYVA，因为他提供了代码。[SYVA: A program to analyze symmetry of molecules based on vector algebra](https://www.sciencedirect.com/science/article/pii/S0010465517300310)。刚才试了一下，发现这个分子对称和化学等价还有区别。比如苯环的一个位点取代，我们一般认为苯环上其他对称的氢是化学位移相同的，但是如果分子不对称，这样子的条件就不满足SYVA的归类了。

Program：

- 去掉氢原子的结构进行对称性判断，对称的原子上面的氢原子都是化学等价的
- 判断同一个碳原子上面的氢，是否是非对映质子


回答标题：首先通过对称性判断化学等价性，然后考察是否磁等价；这之外需要考察快速交换导致的化学等价性质。

**if a plane of symmetry or an axis of symmetry makes two or more nuclei equivalent by symmetry, then they are said to be chemically equivalent.** 这个出 Rotational axis or Symmetry plane (or by a rapid exchange mechanism)

困难在判断磁等价，以及快速交换的化学等价。
![](https://1.bp.blogspot.com/-H-hOyemf4HY/X3Mqgm_3LSI/AAAAAAAALoI/fvKVUli0j3EB8gMteSLeopVNBaTsnD1XACLcBGAsYHQ/w640-h360/pic%2B4.tif)

![](https://1.bp.blogspot.com/-jl-v1bt-9KM/X3MqpZUT_xI/AAAAAAAALoM/kYv6JR7BEdU8EODIteHupzK9eNc_bpLzACLcBGAsYHQ/w640-h264/Pic%2B8.jpg)

![](https://1.bp.blogspot.com/-Sj9ofM95cGQ/X3iRQP_3n5I/AAAAAAAALp8/L3dJvHpc1iYRyXYdAynBILxsQuTckogewCLcBGAsYHQ/w640-h360/pic%2B7.jpg)

![](https://1.bp.blogspot.com/-Oy8CsSp74EA/X3iRoN6aURI/AAAAAAAALqE/pPIPSDJjtm8Zo2swcTd1KiiJjZ0nB-q3gCLcBGAsYHQ/w640-h176/Pic%2B14.tif)

#### 13C谱的个数

#### 判断价键的旋转（rotateble bond）

The original Veber papers (J. Med. Chem. 2002, 45, 2615) define what they consider to be a rotatable bond: Rotatable bonds were defined as any single bond, not in a ring, bound to a nonterminal heavy (i.e., non-hydrogen) atom. Excluded from the count were amide C-N bonds because of their high rotational energy barrier.

使用rdkit的RotatableBondSmarts可以来判断哪些键是否是可以转动的。https://www.rdkit.org/docs/Cookbook.html