跳转到内容
← 返回知识库
核物理约 19 分钟阅读

核磁共振

核物理核磁共振NMRMRI核磁矩拉莫尔进动化学位移自旋弛豫医学成像拉比

关键词

核磁共振; NMR; MRI; 核磁矩; g因子; 拉莫尔进动; 射频脉冲; 化学位移; 自旋弛豫; T1; T2; 傅里叶变换; 拉比; 布洛赫; 珀塞尔

MRI扫描有辐射吗?"核"字引发的误解

标题:"MRI(磁共振成像)里有辐射吗?"——NMR用的是无害的射频波,不是X射线

医院里的核磁共振成像(MRI,magnetic resonance imaging)经常让患者感到不安——"磁共振""核"这个字让人联想到核辐射。实际上,核磁共振(NMR/MRI)不使用任何电离辐射:它用的是射频电磁波(频率约1—1000 MHz),能量约 10−710^{-7}—10−410^{-4} eV,远远低于可见光的能量(约1—3 eV),更远低于可以电离原子的X射线($> 10$ eV)。射频波不能破坏DNA,不会引起辐射损伤。

"核"磁共振中的"核"指的是原子核(如氢原子核,即质子),不是"核辐射"。NMR利用的是原子核的量子自旋——在外磁场中,带自旋的原子核像小陀螺一样进动,其进动频率(拉莫尔频率)与磁场强度成正比,可以被特定频率的射频波共振激发。

射频落在电磁波谱的无线电窗口里,和可见光、X 射线是同一类场,只是频率低了大约九到十二个数量级。一份射频光子按 E=hνE=h\nu 只有约 10−710^{-7} 到 10−410^{-4} 电子伏特,到不了打断化学键的门槛;它能做的,是推动已经在拉莫尔进动中的核磁矩翻一个身。把 MRI 的「核」读成核辐射,等于把核自旋共振和电离径迹叠在同一个恐惧词上。

核磁共振是现代分析化学(测定有机分子结构)和医学成像(MRI)最重要的工具之一,也是精密测量(质子磁矩、基本常数)和量子计算(利用核自旋作为量子比特)的平台。

从分子束到病房:四次诺贝尔奖背后的NMR发现史

标题:1938—1952年:拉比、布洛赫与珀塞尔的三次NMR诺贝尔奖

核磁共振的历史从精密的分子束实验开始:

1938年,伊西多·拉比(Isidor Rabi,1898—1988)在哥伦比亚大学用分子束磁共振方法,首次精确测量了氢核和氘核的核磁矩,并证明了核磁共振的存在——分子束通过非均匀磁场时,特定频率的射频场可以使核自旋翻转,产生可探测的信号变化。拉比于1944年获诺贝尔物理学奖。

1946年,菲利克斯·布洛赫(Felix Bloch)(斯坦福)和爱德华·珀塞尔(Edward Purcell)(哈佛)各自独立在液态和固态样品中实现了核磁共振的直接探测:布洛赫用感应线圈探测进动核的磁场;珀塞尔通过测量射频场在共振条件下的功率吸收来探测。两人于1952年共同获诺贝尔物理学奖。

拉比的分子束已经把核磁矩写成可测的进动频率,布洛赫和珀塞尔要做的是让整块液体或固体自己发声,而不是再拆成一束分子。感应线圈里的电压,和吸收曲线上的功率缺口,量的是同一件事实:核自旋在静磁场里按拉莫尔频率集体进动。没有这条进动,后面的化学位移和成像对比都没有频率可编码。

1966年,理查德·恩斯特(Richard Ernst)引入了傅里叶变换NMR:用宽频脉冲同时激发所有频率,再通过傅里叶变换获得全频谱,将NMR灵敏度提高了数百倍,使高分辨NMR成为有机化学的标准工具。恩斯特于1991年获诺贝尔化学奖。

恩斯特把「扫一个频率」换成「一记宽脉冲再做频谱」,灵敏度跳了两个数量级,物理对象却没有变。谱仪看见的仍是 ω0=γB0\omega_0=\gamma B_0 附近那一小簇频率;劳特伯后来叠上梯度,只是让不同位置唱不同的拉莫尔调。化学家不加梯度时,同一套进动给出化学位移,见核磁谱学结构解析;加上梯度如何变成图像,见核磁共振成像MRI原理。时域怎么换成频谱,约束见傅里叶分析。

1973年,保罗·劳特伯(Paul Lauterbur)提出用梯度磁场对NMR信号进行空间编码,实现了第一幅NMR图像(MRI的雏形)。彼得·曼斯菲尔德(Peter Mansfield)发展了回波平面成像(EPI),实现了快速MRI扫描。两人于2003年共同获诺贝尔生理学或医学奖。

拉莫尔进动、化学位移与T₁/T₂:成像对比度的物理根源

标题:拉莫尔进动、射频脉冲与化学位移

核磁共振的物理基础是原子核的量子自旋在外磁场中的行为:

核磁矩:带自旋(I≠0I \neq 0)的原子核具有核磁矩 μ=γNℏI\boldsymbol{\mu} = \gamma_N \hbar \mathbf{I},其中 γN\gamma_N 是旋磁比(gyromagnetic ratio,不同核不同,γH=2.675×108 rad/(s\cdotT)\gamma_H = 2.675 \times 10^8\ \text{rad/(s\cdot T)} 对氢核)。

氢核自旋 $1/2$、天然丰度接近百分之百、旋磁比又大,所以它是 NMR 和 MRI 的默认探针。氮-14 自旋是 1,有电四极矩,线更宽;碳-13 自旋 $1/2$ 但丰度只有约 1.1%,要靠富集或长时间累加。拉莫尔频率表其实是一张核素身份证,不是一份检查菜单。

拉莫尔进动:在外磁场 B0B_0(沿z轴)中,核磁矩以拉莫尔频率进动: ω0=γNB0\omega_0 = \gamma_N B_0 对氢核,B0=1 TB_0 = 1\ \text{T}(1特斯拉)对应 f0=ω0/(2π)≈42.6 MHzf_0 = \omega_0/(2\pi) \approx 42.6\ \text{MHz};临床MRI常用磁场为1.5 T(f0≈64 MHzf_0 \approx 64\ \text{MHz})、3 T(f0≈128 MHzf_0 \approx 128\ \text{MHz})或7 T(f0≈300 MHzf_0 \approx 300\ \text{MHz})。

拉莫尔公式在经典图像里是磁矩绕磁场进动的角速度,在量子图像里就是相邻磁能级的间隔 ΔE=ℏγB0\Delta E=\hbar\gamma B_0。氢核的旋磁比是核自己的属性,所以 1 特斯拉对应约 42.6 兆赫,不是厂家挑选的工作频率。临床常用 1.5 特斯拉和 3 特斯拉,只是把同一条比例尺放大,共振窗口仍落在调频广播附近的射频段。

化学位移(chemical shift):同一种核(如1^1H)在不同化学环境中,周围电子云屏蔽外磁场的程度不同,使实际感受的磁场略有偏移,导致共振频率偏移(通常用ppm表示,相对参考频率的百万分之一)。这使得NMR谱能区分分子中不同位置的氢原子(或13^{13}C、31^{31}P等),是有机分子结构测定的核心工具。

核实际感受到的不是裸露的 B0B_0,而是被电子云屏蔽后的 (1−σ)B0(1-\sigma)B_0,σ\sigma 随化学环境改变,于是共振频率变成 ω=γB0(1−σ)\omega=\gamma B_0(1-\sigma)。相对主频率的偏移只有百万分之几,谱图因此用 ppm 标度,并约定以四甲基硅烷一类参考物为零点。1951 年 Arnold、Dharmatti 与 Packard 在乙醇的核感应信号里看见三组峰,分别对应甲基、亚甲基和羟基上的氢。Ramsey 在 1950 年给出分子磁屏蔽的量子理论,化学位移从此可以计算,而不只是经验花样。

化学位移把峰的位置写成电子环境,自旋-自旋耦合(J 耦合)把峰的裂分写成键连关系。J 以赫兹计、不随 B0B_0 变,化学位移以 ppm 计、赫兹分裂随场强放大——高场谱仪更「散得开」,靠的就是这一条。一张一维氢谱能读骨架,常常是位置加裂分一起用,不是只看哪台机器场强更高。

弛豫时间:射频脉冲激发后,核磁化恢复平衡的过程由两个时间常数描述:

  • T1T_1(纵向弛豫时间,自旋-晶格弛豫):纵向磁化恢复时间,与分子运动/声子环境有关(典型MRI中,人体不同组织 T1∼0.3T_1 \sim 0.3—3 s)。
  • T2T_2(横向弛豫时间,自旋-自旋弛豫):横向磁化衰减时间,总是 T2≤T1T_2 \leq T_1(典型 T2∼0.03T_2 \sim 0.03—3 s)。

T1T_1 描述沿 B0B_0 的纵向磁化如何把能量交还给周围分子运动,T2T_2 描述横向磁化如何因核与核之间的磁场涨落而散相。Bloembergen、Purcell 与 Pound 在 1948 年把这两套速率写成分子相关时间的函数:运动太慢或太快,弛豫都变慢,中间有一个匹配拉莫尔频率的窗口。T2T_2 总是短于或等于 T1T_1,因为凡是让能量回晶格的过程都会拆掉横向相干,反过来却不必。组织之间的对比来自这两套时间常数本来就不一样,读出方式只是让已经存在的物理差更好认。

MRI成像利用不同组织(肌肉、脂肪、脑灰质/白质、肿瘤等)T1T_1、T2T_2值的差异产生对比度,可以在不同采集参数下(T1加权、T2加权、质子密度加权)突出不同组织结构。

把对比说成 T1 加权或 T2 加权,容易让人以为那是两套互不相干的机器模式。物理上它们是同一条磁化矢量的两个投影:一个沿静磁场恢复,一个在垂直平面里衰减。脂肪和水的氢核处在不同的相关时间窗口里,同样的拉莫尔激发会留下不同的恢复曲线;灰度是弛豫的读出,不是电离辐射的读出。与 X 射线成像如何分家,见CT与磁共振。梯度线圈把位置写成频率差的同时,也会在体内感应涡流与外周神经刺激;切换越快,编码越密,却也越接近安全上限。序列设计因此总是在分辨率、扫描时间与病人耐受之间折中,而不是把场强数字单独往上堆。

场强越高越好吗?7T乃至10.5T的代价与争议

标题:高场MRI的优势与代价——7 T对人体的影响

MRI的信号-噪声比(SNR)与磁场强度 B0B_0 大致成正比(精确地,与 B07/4B_0^{7/4} 成正比),因此提高场强是提升图像质量的直接手段。临床上目前标准为1.5 T和3 T;研究型系统已达7 T(已获FDA批准用于特定临床诊断),甚至10.5 T(明尼苏达大学研究型扫描仪)。

高场MRI的代价包括:

  • 射频功率沉积(SAR,specific absorption rate):人体吸收的射频功率与 B02B_0^2 成正比,7 T下的SAR限制使某些成像序列无法使用。
  • 化学位移伪影:磁场越强,不同组织间的化学位移越大,边界伪影越明显。
  • 磁化率效应:高场下不同组织间磁化率差异引起的磁场不均匀性增加,使几何畸变和信号丢失更严重。
  • 安全顾虑:高场下金属植入物(心脏起搏器、金属骨钉)的安全约束更严格;极高场(>8 T)对人体的生物效应尚在研究中。

场强升高时,化学位移以赫兹计的分裂按 B0B_0 放大,ppm 标度本身不变。这就是高场谱仪更「散得开」的原因,也是成像里脂肪与水边界伪影变明显的同一条账。信噪比随场强上升,是因为玻尔兹曼极化 μB0/kT\mu B_0/kT 变大,不是因为射频变得更「硬」。

对于是否应将MRI场强继续推向14 T甚至更高,学界有争议:支持者认为超高场将使脑功能成像、代谢成像分辨率大幅提高(神经科学研究的核心工具);反对者认为工程挑战和安全风险尚未解决,成本效益有待评估。超导主磁体要靠液氦维持低温,失超保护与匀场垫片的工程量随孔径与场强一起上涨;临床机房的屏蔽、急停与植入物筛查清单,往往比谱仪本身更能决定一台机器能不能真正开诊。

从蛋白质结构解析到量子比特:NMR跨越半个世纪的影响

标题:从蛋白质结构到量子计算——NMR的跨领域影响

NMR的影响早已超出核物理和医学:

有机化学:高分辨液态NMR(500 MHz以上,现在已达GHz级别)是测定有机分子和蛋白质三维结构的标准手段之一。库尔特·维特里希(Kurt Wüthrich)开发了蛋白质溶液NMR结构测定方法,于2002年获诺贝尔化学奖。目前PDB数据库中约12%的结构由NMR解析。

溶液里蛋白质的谱线宽度受 T2T_2 限制:分子越大、转动越慢,横向磁化散得越快,峰就越宽。维特里希的结构测定能走通,靠的是核欧沃豪斯效应一类交叉弛豫,而不是把场强无限制往上推。核自旋后来被拿来演示量子算法,也是因为室温下 T1T_1、T2T_2 可以到秒量级,远长过电子自旋的相干时间。

MRI的医学革命:全球目前有约40000台MRI扫描仪,每年进行超过4亿次扫描。MRI不仅用于解剖成像(脑、心脏、关节、腹部),还用于:功能MRI(fMRI,通过血氧水平依赖信号探测神经活动)、弥散加权成像(DWI,早期脑梗死诊断的金标准)、MR血管成像(MRA)、MR波谱(MRS,无损测量体内代谢物浓度)。

量子信息:原子核自旋是最早实现量子信息处理的系统之一。2001年,IBM使用13^{13}C和1^{1}H核自旋的NMR系统,在7个量子比特上演示了肖尔算法分解15=3×5——尽管NMR量子计算因可扩展性问题后来被其他平台超越,但它为量子计算的早期实验探索提供了宝贵平台。

事实卡

  1. 卡1

    核磁共振由拉比(1938年,分子束)、布洛赫、珀塞尔(1946年,液/固态)分别发现,三人先后获诺贝尔物理学奖(1944年、1952年)。

  2. 卡2

    MRI不使用任何电离辐射,依赖射频电磁波(约1—300 MHz)激发原子核磁共振,利用T1T_1/T2T_2弛豫时间差异区分不同组织。

  3. 卡3

    全球约40000台MRI扫描仪,每年进行超过4亿次扫描;MRI临床常用场强为1.5 T和3 T,最高研究型扫描仪已达10.5 T。

  4. 卡4

    恩斯特因傅里叶变换NMR获1991年诺贝尔化学奖;劳特伯和曼斯菲尔德因MRI技术获2003年诺贝尔生理学或医学奖。

关键洞察

核磁共振把原子核的量子自旋变成了无损透视人体的工具。

它探测的是原子核(主要是人体内丰富的氢核)在强磁场中按拉莫尔频率进动、并对射频脉冲产生共振吸收的量子效应;不同组织中氢核的弛豫时间(T1、T2)不同,于是磁共振成像(MRI)能在不使用电离辐射、不切开身体的情况下分辨软组织。从 1938 年拉比的分子束实验,到 1946 年布洛赫与珀塞尔在液/固态中直接探测,再到 1973 年劳特伯实现成像——一个最初纯属基础物理的现象,最终成了现代医学最重要的诊断手段之一。

跨域连接

  • 磁共振成像:给均匀磁场叠加线性梯度,不同位置的核就有不同的共振频率——频率变成位置的编码,NMR 从一台谱仪变成一台相机。推论:空间分辨率由梯度强度与带宽决定,成像的一切参数都是这条编码关系的化身。
  • 核磁谱学结构解析:同类核在不同化学环境里被电子云屏蔽的程度不同,共振频率偏移百万分之几——这个化学位移让谱峰位置变成官能团的指纹。推论:一张一维氢谱就能读出分子骨架,结构化学因此不必每次都依赖晶体。
  • 傅里叶分析:用一记宽频脉冲同时激发所有共振,再把自由感应衰减信号做傅里叶变换还原成谱——灵敏度因此提高数百倍。推论:谱的分辨率由采样时长决定、谱宽由采样间隔决定,时域与频域的这笔交换是硬性约束。
  • 量子计算理论:核自旋是天然的双态系统,七个自旋的 NMR 系统曾演示肖尔算法分解 15。推论:弛豫时间 T₁/T₂ 就是这类量子比特的寿命上限——NMR 平台后来因扩展性被超越,但它最早把量子算法跑在了实物上。
  • CT与磁共振:X 射线成像看的是密度差异,代价是电离辐射;磁共振看的是 T₁/T₂ 弛豫差异,用的是能量远低于可见光的射频波。推论:软组织的对比度来自弛豫时间的物理差异而非吸收系数——这是两种成像原理分工的根源,也解释了磁共振为何没有剂量限制。

参考文献

  1. Rabi, Isidor I., et al. "A New Method of Measuring Nuclear Magnetic Moment." Physical Review, 53(4): 318, 1938.
  2. Bloch, Felix. "Nuclear Induction." Physical Review, 70(7–8): 460–474, 1946. DOI: 10.1103/PhysRev.70.460
  3. Purcell, Edward M., Henry C. Torrey, and Robert V. Pound. "Resonance Absorption by Nuclear Magnetic Moments in a Solid." Physical Review, 69(1–2): 37–38, 1946.
  4. Lauterbur, Paul C. "Image Formation by Induced Local Interactions: Examples Employing Nuclear Magnetic Resonance." Nature, 242(5394): 190–191, 1973. DOI: 10.1038/242190a0
  5. Bloembergen, Nicolaas, Edward M. Purcell, and Robert V. Pound. "Relaxation Effects in Nuclear Magnetic Resonance Absorption." Physical Review, 73(7): 679–712, 1948. DOI: 10.1103/PhysRev.73.679
  6. Ramsey, Norman F. "Magnetic Shielding of Nuclei in Molecules." Physical Review, 78(6): 699–703, 1950. DOI: 10.1103/PhysRev.78.699
  7. Arnold, J. T., S. S. Dharmatti, and M. E. Packard. "Chemical Effects on Nuclear Induction Signals from Organic Compounds." The Journal of Chemical Physics, 19(4): 507, 1951.

延伸阅读

  1. Levitt, Malcolm H. Spin Dynamics: Basics of Nuclear Magnetic Resonance (2nd ed.). Wiley, 2008.
  2. 赵南明, 周海梦. 《生物物理学》. 高等教育出版社, 2000. (包含NMR生物医学应用章节)