python使用 psi4 计算分子 NMR

要使用 psi4 计算分子 NMR，需要安装 psi4 和相关的 Python 库。

以下是一个简单的示例代码，演示如何使用 psi4 计算甲烷的 NMR 谱：

import psi4

# 设置 psi4 计算选项
psi4.set_options({'basis': '6-31G*', 'reference': 'RHF'})

# 定义甲烷分子
mol = psi4.geometry("""
0 1
C  0.000000  0.000000  0.000000
H  0.628319  0.628319  0.628319
H -0.628319 -0.628319  0.628319
H  0.628319 -0.628319 -0.628319
H -0.628319  0.628319 -0.628319
""")

# 运行 NMR 计算
nmr = psi4.prop('nmr(h2o)', mol, return_wfn=True)

# 提取化学位移和各个核的张量
shieldings = nmr["nuclear_shielding"].np
tensors = nmr["nuclear_spin_tensor"].np

# 输出结果
print("Chemical shift (ppm): ", shieldings)
print("Nuclear spin tensors: ", tensors)

在这个例子中，我们使用 psi4.prop() 函数运行 NMR 计算，并将结果存储在变量 nmr 中。我们可以使用 nmr["nuclear_shielding"] 和 nmr["nuclear_spin_tensor"] 提取化学位移和各个核的张量，并将它们存储在变量 shieldings 和 tensors 中。

最后，我们使用 print() 函数输出化学位移和张量的结果。