-
Superconductivity is a macroscopic quantum phenomenon. Flux quantization and the Josephson effect are two physical phenomena which can best reflect the macroscopic quantum properties. Superconducting quantum interference device (SQUID) is one type of superconducting devices which uses these two characteristics. SQUID devices are widely used in the sensitive detection of magnetic signals. This paper briefly introduces the background and recent developments of low temperature superconductor and high temperature superconductor SQUID devices.
1. 引 言
根据BCS理论, 超导体中的载流子在临界温度Tc以下会形成库珀对, 并发生凝聚形成一种宏观有序态. 这种有序态可以用一个宏观波函数(序参量)描述. 事实上, 在BCS理论出现之前, London[1]就指出超导是一种宏观量子现象. 他还进一步预言一个闭合超导环所包含的磁通是量子化的, 是磁通量子Φ0的整数倍, 其中磁通量子是普朗克常数h除以载流子的有效电荷[1]. 这一预言后来得到了很好的实验验证, 实验结果显示Φ0 = h/2e ≈ 2.0678 × 10–15 Wb[2,3] (这里e为电子电荷). 分母为2e表明了超导体中的超导载流子是两个电子配对形成的库珀对. 高温超导材料发现以后, 实验上也很快验证了高温超导闭合环路中的磁通同样是量子化的[4], 最小单元也是Φ0. 磁通量子化是由超导宏观波函数的单值性决定的, 即绕环路一圈回到原点后波函数的值不变, 相位变化须是2π的整数倍.
1962年约瑟夫森预言[5], 被一层很薄绝缘体分开的两个超导体形成的超导-绝缘-超导(S-I-S)隧道结中, 两侧的超导波函数会穿过绝缘层彼此重叠, 使得库珀对可以相干地隧穿通过绝缘层, 也就是说超导电流可以无阻地流过隧道结. 超导电流Is和绝缘层两端超导体波函数相位差φ之间的关系为
Is=Icsinφ. (1) 如果结上有电压V, 则φ随时间的变化率为
∂φ∂t=2eℏV. (2) 约瑟夫森的理论预言很快为实验所证实[6], 这一效应也就称之为约瑟夫森效应. 进一步的研究显示约瑟夫森效应不但在S-I-S隧道结上可以观察到, 在其他形式的弱连接结上也可观察到.
根据(1)式和(2)式很容易计算得到一个约瑟夫森结的耦合能为
EJ=ℏIc2e(1−cosφ). (3) 另外, 也可看出约瑟夫森结可以看作一个非线性电感, 电感值为
LJ=ℏ2eIccosφ=Lc1cosφ. (4) 自上世纪60年代约瑟夫森效应和磁通量子化现象被实验所证实以来, 基于这两种现象的超导器件和应用不断出现, 在许多领域发挥着重要的作用. 本文介绍的超导量子干涉器件(superconducting quantum interference device, SQUID)是其中一个主要的方面. 关于SQUID的研究已有五十多年的历史, 是一个发展相对成熟的方向, 有很多的综述文献和书籍可供参考[7—14].
2. 基本原理
SQUID器件通常由包含一个或两个约瑟夫森结的超导环路构成. 一般称前者为射频SQUID(RF-SQUID), 后者为直流SQUID(DC-SQUID). 当穿过SQUID环路的磁通Φ发生变化时, SQUID器件的电感(射频SQUID)或临界电流(直流SQUID)会发生相应的变化, 因此可以用射频和直流方法分别来探测外界磁通变化的影响. 利用磁通耦合器, 任何可以转化为穿过SQUID环路磁通的被测量都可以利用SQUID传感器来测量. 例如, 用SQUID可以做成极其灵敏的磁强计或电流传感器等. 屏蔽电流或穿过SQUID超导环路的磁通使SQUID电压-磁通或电流-磁通曲线出现周期性振荡, 振荡的周期就是一个磁通量子Φ0.
图1是一个DC-SQUID的示意图. 假定流过SQUID的电流是I, 超导环路中的环流是J, 两个约瑟夫森结的临界电流都是Ic0, 相位差分别是φ1和φ2. 根据约瑟夫森方程, 流过两个结的电流满足如下关系:
12I+J=Ic0sinφ1, (5) 12I−J=Ic0sinφ2. (6) 根据超导环路量子化关系(假定超导环路电感很小, 可以忽略)
φ1−φ2=2πΦΦ0. (7) 定义SQUID的相位差为φ = (φ1 + φ2)/2, 根据(5)式—(7)式可以得到
I=2Ic0|cos(ΦΦ0π)|sinφ, (8) J=Ic0sin(ΦΦ0π)cosφ. (9) 从方程(8)可以看出, SQUID等效于一个临界电流可变的约瑟夫森结, 临界电流
Ic=2Ic0|cos(ΦΦ0π)| (10) 随磁通周期变化. 显然, 出现这种现象的原因是磁通产生的环流电流使得SQUID环路两侧约瑟夫森结上的相位差不同而引起干涉效应, 这也是这类器件称为量子干涉器件的原因. Jaklevic等[15]最早在实验上观察到这种量子干涉现象.
在大多数应用中, SQUID器件中超导环路的几何电感L不能忽略. 通常用屏蔽参数
βL=2LIc/Φ0 , 来描述环路几何电感的影响. 图2所示的是不同βL参数值对应的磁通对DC-SQUID器件临界电流调制的情况.由于磁通对SQUID临界电流的调制, 当SQUID电流偏置在略大于临界电流的状态时, SQUID两端的电压也会受磁场调制而在最大值和最小值间以Φ0为周期变化. 图3是根据DC-SQUID的等效电路给出的不同外加磁通Φa下的I-V曲线, 以及对应的结上电压随磁通变化的情况. 其中和约瑟夫森结并联的电容C是结超导电极和绝缘层形成的电容, 电阻R是外加的并联电阻或本征的并联电阻. 合适的R值保证约瑟夫森结的I-V曲线无回滞. 电压的变化幅度表明了磁通对SQUID电压的调制的大小, 被称为电压调制深度. 而电压变化与磁通变化之比ΔV/ΔΦ被称为SQUID的传递函数. 因此, SQUID本质上是一个磁通电压的转换器.
图 3 DC-SQUID (a)等效电路示意图; (b)磁通分别为整数个和半整数个Φ0时的I-V曲线; (c)电压-磁通曲线[7]Fig. 3. The DC-SQUID: (a) Schematic electric circuit; (b) current-voltage characteristics at integer and half-integer values of applied flux; the operation point is set by the bias current Ib; (c) voltage vs. flux Φa/Φ0 for constant bias current[7].SQUID可以测量远远小于一个磁通量子的磁通. 但是由于SQUID的环孔很小, 因此对磁场强度的测量并不灵敏. 实际应用中, 一般将SQUID与不同形式的探测线圈耦合起来构成SQUID磁强计或梯度计. 磁强计和梯度计可以测量很微弱的磁场或磁场梯度信号, 因此可以应用在很多领域.
除传统的直流和射频SQUID器件外, 另外一种称为超导量子干涉滤波器(superconducting quantum interference filter, SQIF)的器件近年也受到很多关注. 这种器件包含许多连接在一起的SQUID构成的阵列, 其中每个SQUID环路的大小不一, 这样使得各个SQUID对应的磁通电流干涉曲线在零磁通以外的区域相互抵消, 而只在零磁通的区域出现一个峰. 因此, 与传统SQUID相比, 电压-磁通曲线不再是周期变化的, 而是只出现一个峰. 这种器件可以用来进行磁场绝对值的测量[16].
由于SQUID电压-磁通曲线是非线性的, 通常需要采用反馈电路来得到线性化的输出, 从而实现磁通测量. 这样的基本电路称为磁通锁定反馈电路(flux-locked feedback loop或flux-locked loop, FLL). 图4给出了两种最常用的DC-SQUID读出电路的示意图. 其中左图是直接读出FLL读出电路, 右图是磁通调制的FFL读出电路. 虽然前一种读出电路更加简单, 但是在SQUID发展的历史上, 却是第二种磁通调制读出电路最早被广泛使用, 因为这种方式工作稳定, 噪声低. 磁通调制可以有效地减少DC-SQUID偏置电流噪声以及前置放大器(Preamplifier)低频噪声的影响. 而且置于SQUID和前置放大器之间的低温电压变压器可以放大SQUID上的对应磁通调制的交变电压信号, 减少了前置放大器本身噪声的影响. 磁通调制读出电路的缺点是电路相对复杂, 工作带宽受到磁通调制频率以及电压变压器的限制. 一般常规的磁通调制读出电路工作的工作频率在约100 kHz以内. 通过电路改进, 工作频率可以提高到几MHz. 随着多通道SQUID系统的发展, 以及半导体放大器性能的不断提高, 直接读出FLL电路由于其电路结构简单、工作带宽大, 而逐渐得到广泛应用.
针对实际应用, 人们还发展出了多种其他形式的读出电路, 可根据噪声、动态范围、频率响应等方面的要求选取合适的读出电路[7,14,17]. 目前利用常规的半导体读出电路, 可以检测10–6Φ0大小的磁通变化.
图5所示的是RF-SQUID和用于读出的谐振电路以及前置放大器的等效电路示意图. 前面已经提到约瑟夫森结可看成一个非线性电感. 当外加磁通变化时, RF-SQUID环路的电感发生变化导致与其耦合的谐振电路的谐振频率和相位等发生变化. 因此可以通过测量谐振电路的输出电压来测量RF-SQUID感受的磁通. 同样可以采用磁通锁定反馈电路来获得线性化的输出. 需要指出的是, 目前绝大多数应用都采用DC-SQUID作为探测器件. 低温超导RF-SQUID应用方面的报道较少, 高温超导RF-SQUID有可能在一些方面得到应用.
3. 噪声特性
SQUID器件本身的噪声主要分为三种: 分路电阻上的热涨落噪声(Johnson噪声)、散粒噪声和1/f噪声. 其中热涨落噪声是由于电子的热运动产生的在分路电阻两端的噪声电压. 散粒噪声是由于电流中运动电子数目的起伏导致的噪声电流.
散粒噪声和热涨落噪声都与频率无关, 被称为白噪声. 一般散粒噪声只有在电流很小时才能表现出来, 对于高温超导器件来说, 热涨落噪声远远大于散粒噪声的影响, 因此一般不考虑散粒噪声的影响.
1/f噪声的方均噪声电压频谱密度与频率的倒数成正比, 因此被称为1/f噪声. 1/f噪声起源还不完全清楚, 但有两个来源可以确定: 一是临界电流和结电阻的涨落, 二是源于薄膜中被钉扎的磁通的运动. 第一种来源的1/f噪声可以通过磁通调制和偏置电流反转技术来抑制[7,17]. 抑制第二种噪声需要采用高质量的超导薄膜和对SQUID进行特殊的设计来减少磁通钉扎. 一种常用的方法是将器件中超导薄膜部分制备成网格状, 使磁通局限在网格的孔中[18].
如果要用SQUID测量磁通, 热噪声的能量必须远小于磁通的能量, 否则被测量的磁通信号就会被热噪声掩盖. 这样就对SQUID的电感有限制, 不能太大. 通常在SQUID设计时, 会考虑工作温度和屏蔽参数βL~1等条件, 再根据约瑟夫森结无回滞条件βC < 1来确定SQUID电感、结的临界电流和正常态电阻等参数, 以使SQUID处于最佳的工作状态[7,19].
对低温SQUID磁强计,白噪声频率区域典型的磁场噪声水平为几个fT/Hz1/2. 高温超导SQUID磁强计的噪声水平变化比较大, 受工艺参数的影响明显. 高质量器件白噪声区间的噪声可以低于50 fT/Hz1/2. 此外, 高温SQUID磁强计的1/f噪声比低温器件显著大, 因此低频噪声增加较多. 采用偏置电流反转的方式可以减小由于结的临界电流涨落导致的1/f噪声.
原则上, 高温超导SQUID器件工作在77K, 相比于工作在液氦温度的低温SQUID器件, 热噪声大了约20倍, 因此高温超导SQUID器件呈现出高的噪声水平似乎也不奇怪. 但是根据SQUID噪声的理论分析, 磁通白噪声与SQUID电感L, 结电容C等有如下关系:
SΦ≈32kBTL(LCβC)1/2, (11) 这里βC = 2πICRn2C ≈ 0.4是约瑟夫森结的Stewart-McCumber参数. 高温超导结的βC参数基本不随结的大小变化, 而低温结中正常态电阻Rn通常由外加并联电阻决定, 因此SΦ与结电容C无关. 基于这些考虑, 由于高温超导约瑟夫森结的电容很小, 因此高温超导SQUID器件仍然可以实现低的噪声[20]. 需要注意的是, 上述讨论不适用于冻结磁通跳跃、结区载流子俘获释放过程以及器件线路中寄生的弱连接等引起的噪声, 而这些常常是影响高温超导SQUID噪声性能的主要因素. 为了获得性能良好的器件, 需要生长高质量的高温超导薄膜, 薄膜要均匀外延、单一取向、致密且孔洞少.
4. 低温超导和高温超导SQUID
根据所采用的超导材料, SQUID器件大致可分为低温超导SQUID器件和高温超导SQUID器件. 本节将从几个方面比较低温和高温超导SQUID器件.
材料 在具体的超导材料方面, 虽然在早期采用不同材料来制备SQUID器件, 但自上世纪80年代以后, 低温SQUID器件一般都用超导金属铌, 采用薄膜平面工艺制备. 铌是难熔金属, 机械性能好, 用铌薄膜制备的约瑟夫森结非常稳定. 高温超导SQUID器件一般采用YBCO/ReBCO超导薄膜制备.
工作温区 低温超导SQUID器件一般工作在4.2 K, 用液氦冷却. 高温超导SQUID器件则工作在77 K, 用液氮冷却. 人们也在研究利用基于小型制冷机的恒温器来冷却SQUID器件, 这样可以方便应用, 但面临如何有效抑制制冷机产生的电磁干扰的困难.
器件制备 SQUID器件中的关键是约瑟夫森结, 低温超导SQUID器件中的约瑟夫森结为Nb/AlOx/Nb三层结构. 其中铝在铌表面浸润性好, 可以很好地覆盖在铌膜的表面, 绝缘层的厚度可以通过控制铝氧化过程的氧分压和时间来很好地控制. 器件制备中可以通过选用合适的刻蚀气体对铌选择性刻蚀[21]. 由于工艺成熟, 器件的参数控制比较容易, 器件的参数一致性和良率都能够得到较好保证.
与之相比, 高温超导约瑟夫森结的制备就要困难得多. 这主要归咎于材料复杂的化学和结构特性、非常短的相干长度、对氧含量敏感等原因, 导致薄膜质量难以控制、平整度和均匀性不好, 很难在高温超导薄膜上制备非常薄的绝缘层. 因此尚没有利用高温超导薄膜制作三层结构的约瑟夫森结的方法. 早期高温超导约瑟夫森结主要有双晶结、台阶结、台阶SNS结(即超导-正常金属-超导结)以及斜边结等形式[22,23], 如图6所示.
利用不同取向的晶界形成的弱连接称为晶界结. 双晶结是将基片定向后, 将两片不同取向的基片对齐后重新生长在一起, 形成双晶基片. 在双晶基片上外延生长高温超导薄膜, 由于薄膜取向与基片取向一致, 便在晶界处形成双晶结. 双晶结的临界电流密度随双晶取向间的夹角增加而指数减小[24]. 通常采用24°, 45°和36°的双晶基片制作SQUID. 用双晶结制备的SQUID器件通常有较好的噪声特性, 但缺点是需要双晶基片, 成本高而且结的位置不能任意选择.
台阶结是在基片上合适位置刻蚀出一个较缓的台阶, 再在基片上外延生长超导薄膜, 在台阶处形成弱连接. 由于薄膜在台阶斜坡上的生长取向和基片上不同, 台阶处的弱连接实际上也是晶界结. 台阶结的优点是可以在基片的任何位置刻蚀台阶, 制备工艺相对简单. 但是台阶结的质量依赖于台阶的微观结构, 因此器件性能可重复性不是很好. 不过, 近十多年里, 随着工艺的不断改进, MgO基片上台阶结的质量一直在提高, 器件制备的一致性有较大提高[25]. 德国Juelich研究中心的研究人员利用MgO基片台阶结获得了白噪声水平达到约5 fT/Hz1/2的高温超导直流SQUID器件, 并利用该器件成功检测出人体脑磁信号[20] (图7).
图 7 MgO衬底上45°台阶上生长的YBCO薄膜的扫描电镜(SEM)图像(左图)和高分辨透射电镜HRTEM图像(中图). 右图: 一个16 mm大小、采用台阶结的高温超导DC-SQUID磁强计在超导屏蔽环境下测量的噪声谱[20]Fig. 7. SEM image (left) and HRTEM image (middle) of an YBCO film deposited on a double-layer-buffered 45° step on an MgO substrate. A 45° [100]-tilted GB is clearly shown. Right: Noise spectral density of a 16 mm high-Tc DC-SQUID magnetometer with step-edge junctions measured in a superconducting shield[20].台阶SNS结是利用光刻和Ar离子刻蚀的方法在基片上刻出一个很陡的台阶. 然后在基片上生长超导薄膜, 由于台阶边缘很陡, 在边缘处没有超导薄膜. 再镀上一层Ag或Au膜覆盖在台阶边缘,形成SNS结. 这种结的正常态电压IcRn在77 K可以做到1 mV, 然而其输运和噪声性质均不理想.
斜边结是先在基片上生长一层超导薄膜和绝缘层, 再刻蚀出一个斜边, 然后沉积一层阻挡层和顶电极. 优点是由于顶电极的屏蔽作用使得其抗干扰能力强, 具有热循环稳定性. 日本的研究人员在2000年前后, 投入了较大的力量研究斜边结, 采用SrSnO3/La-YBa2Cu3O7/SrSnO3/La-YBa2Cu3O7多层结构. 约瑟夫森结的临界电流偏差分布可以在8%以内[26].
除了上述的四种方法外, 人们还采用离子辐照的方式来制备高温超导约瑟夫森结, 并制备约瑟夫森结和SQUID器件[27]. 由于高温超导材料对氧的含量很敏感, 利用几十到几百keV的氦、氖、氮甚至氧等轻元素离子可以有效地在高温超导体中形成氧错位和空位等缺陷, 使得材料中受辐照区域变成非超导的. 最近的研究显示, 利用氩离子辐照在狭缝区域局部去除超导材料形成沟槽型弱连接, SQUID器件的噪声水平可以达到63 fT/Hz1/2, 与双晶结SQUID相当[28]. 利用有纳米尺度狭缝的金属掩模进行局部辐照的方式, 可以根据需要在指定位置制备出约瑟夫森结. 此外, 近年来随着聚焦离子束技术的发展, 利用聚焦氦离子直接局部轰击可以有效形成约瑟夫森结[29,30].
磁通变换器 为了获得高的磁场灵敏度, SQUID通常需要通过磁通变换器和一个更大面积的探测线圈耦合. 对低温超导SQUID, 可以用低温超导线绕制探测线圈, 也可以利用薄膜制成单圈或多圈平面探测线圈. 探测线圈可以是一阶或二阶梯度计的形式, 便于减少环境干扰. 磁通变换器一般做成平面多圈形式, 置于SQUID环孔上(图8).
对高温超导SQUID, 绕制探测线圈和利用多层薄膜技术制备具有良好超导连接的多圈磁通变换器很困难. 目前还没有柔性且易于形成超导连接的高温超导线, 因此很难绕制探测线圈. 另外, 高温超导薄膜材料需要在结晶性好的取向基底上才能很好地外延生长, 要制备类似图8所示的多圈磁通变换器比较困难, 工艺难度大. 国际上只有少数研究组能够制备出, 器件良率不易提高. 因此, 高温SQUID磁强计多采用探测线圈与SQUID回路直接耦合的方式, 或者采用倒扣的方式将探测线圈和SQUID耦合.
应用 有关SQUID器件应用更加详细的介绍将在下一节中给出. 这里只是简短地对比一下低温超导和高温超导SQUID器件在实际应用场景中的区别. 经过50多年的发展, SQUID器件的应用研究已经延伸到许多领域. 原则上在绝大部分应用中, 既可以使用低温超导SQUID器件, 也可以使用高温超导SQUID器件. 而且, 高温超导器件可以用廉价和容易获得的液氮冷却, 易于应用推广. 因此, 高温超导体发现以来, 高温超导SQUID器件的制备和应用研究一直是一个重要的研究方向, 经过30多年的努力, 器件的性能也一直在不断提升中, 报道的最好指标与商用低温超导SQUID器件的指标已经比较接近了. 但是高性能器件的批量化产出则较多地受到参数的一致性、良率和较高成本的影响. 目前, 在仅需要单个或少量SQUID器件的一些应用中, 已经使用高温超导SQUID器件. 而在需要大批量SQUID器件的系统中, 如脑磁图测量系统, 基本上还是使用低温超导SQUID器件.
5. SQUID应用
自上世纪60年代SQUID出现以来, SQUID的应用研究一直在不断地持续开展. SQUID器件具有以下特点: 1)接近量子极限的磁探测灵敏度, 可以探测低至10–15 T/Hz1/2的磁场信号; 2)具有很宽的频响特性, 由于超导体对磁场而非磁场的变化率产生响应, SQUID器件可以探测直流到GHz的信号, 通常测试系统的频响受读出电路的限制, 一般可达10—100 MHz量级; 3)很大的动态范围, SQUID器件的动态范围很容易达到130 dB水平, 磁信号幅度变化6—7个量级不会影响系统的工作稳定性; 4)很高的线性响应特性. 由于这些特点, 几十年来, SQUID器件的应用研究已经涉及到医学科学、科学仪器、地球物探、工业无损检测以及精密测量等许多领域. 关于已有的SQUID应用领域可以参考相关文献[8,10].
5.1 生物磁探测
生物体内血液和神经中离子流动会产生非常微弱的磁信号. 利用SQUID的高灵敏磁探测能力, 可以对人体或其他生物体内不同器官和组织产生的微弱磁信号进行无创无接触检测, 在疾病检测和医学研究中极具应用前景. 上世纪70年代起, 就已经开始相关应用研究, 首先是心脏图和脑磁图的测量, 后来也拓展到身体其他部位, 如肺、肠、胃、脊椎和胎儿等的磁信号[11]. 由于生物磁探测在医学健康领域的巨大应用前景, 相关研究是目前SQUID应用研究领域最大的一个分支.
5.1.1 脑磁信号探测
大脑神经活动产生的磁场信号可为研究大脑的行为、脑部疾病产生的原因和确定病灶等提供非常有用的信息, 在脑科学和脑疾病研究中越来越受到重视[31]. 脑磁图(magnetoencephalography, MEG)系统大多采用多路信号探测器, 收集大脑外表皮各处的磁信号. MEG与研究大脑行为常用的脑电图(EEG)和脑部功能磁共振成像系统相比, 各有优缺点, 作用具有互补性. 与脑电图相比, 脑磁图测量具有无接触无创的优势, 而且磁信号和电信号是很好的互补. 与磁共振成像相比, MEG虽然不能提供组织的成像, 但具有毫秒量级的时间分辨, 可以获得快速的时域响应信息.
人的脑磁信号幅度约在1—2 pT的范围, MEG测量需要fT的磁场分辨率. 因此, SQUID器件非常合适于MEG应用. 目前SQUID器件是MEG系统中采用最广泛的探测器件, 这一领域也是SQUID最大的应用领域[12]. 1993年, 第一台多通道全头型MEG的出现, 标志着MEG医学临床应用研究的开始. 据估计, 目前约有250台脑磁图商业系统分布在世界各地, 主要集中在发达国家. 这些系统包含总计约50000个SQUID器件.
利用MEG, 头部重要功能皮质区(eloquent cortex)产生的诱发场、脑震荡源和癫痫活动可以在一次检测中确定, 也可以研究大脑两个半球之间的相互作用[32]. 二十多年来MEG的研究内容也在不断演变和深入, 从感知、处理、分析大脑不同区域的功能连通性, 到诸如婴儿期到成年期的大脑发育相关问题研究等[31]. 目前, 癫痫患者的手术前评估, 脑肿瘤患者雄辩皮质的鉴定或血管畸形已被公认为临床MEG研究的适应症[33].
由于MEG具有很好的时间分辨性和一定的空间分辨率, 在研究大脑重要功能皮质层不同区域的功能和彼此之间的连接性方面能够发挥重要作用, 因此MEG在注意、记忆等认知功能和神经紊乱性疾病等脑科学和脑疾病的研究、诊断以及治疗疗效跟踪和评估等方面可发挥重要作用. 图9是SQUID脑磁测量装置以及研究视觉感知与脑磁信号关系的示意图和一些结果.
全球有数亿人患有癫痫、中风、阿尔兹海默及抑郁症等与大脑相关的疾病, 形成巨大的社会负担. 通过早期的准确检测和诊断方式改进, 可以早发现早治疗,从而有可能减少这些疾病造成的损害和负担. 因此, 随着SQUID脑磁系统的不断改善和技术的进步, SQUID的潜在需求有可能增加十倍甚至一百倍.
SQUID-MEG系统的进一步广泛推广应用需至少在以下两方面做大量工作. 首先是大量临床研究数据的积累和分析, 以期建立更多脑磁信号与特定脑功能或脑疾病之间的对应关系, 以便临床应用; 此外还需在提高测量脑磁信号的信噪比、降低制造和使用成本等方面做出不断的改进. 其中制造成本在一定程度上可以随着需求的增加而降低. 目前的SQUID-MEG系统采用工作在液氦温度的低温超导SQUID器件, 液氦的使用是使用成本高的重要原因. 近年来, 人们已经研发了利用小型制冷机重新液化挥发氦气的闭路恒温器[34], 这样可以减少液氦消耗和补充方面的成本.
使用液氦的另一个缺点是低温恒温器的存在使得SQUID探测器和信号源之间的距离(典型值约为30 mm)变大了, 既减少了信号强度又降低了空间分辨率. 近年来, 随着高温超导SQUID器件性能的不断提高, 利用高温SQUID构建MEG系统的研发工作受到更多关注[20,35]. 由于工作在液氮温区, 低温恒温器壁厚可大大减小到几毫米甚至1 mm. 这样, 高温超导SQUID器件可以更加接近大脑表皮, 提升空间分辨率, 同时又增加了信号强度.
SQUID器件在MEG中的应用也面临其他技术的竞争. 高灵敏原子磁强计可以具有与SQUID器件相当的磁场灵敏度, 基于这种器件的MEG系统正在发展中. 其最大的优势是不需要低温系统, 可望显著降低运行成本. 系统重量也可大大减少, 探测头盔可以做成穿戴式. 此外, 器件可以更加靠近大脑表皮, 获得更好的空间分辨率[36,37]. 但是原子磁强计工作在标量模式, 无方向分辨. 而且其频率响应范围小, 只能工作在约30 Hz以下的频率范围, 动态范围也不大. 由于这些原因, 未来的MEG可能是两种探测器并存, 根据应用需求的不同而决定采用原子磁强计还是SQUID器件.
5.1.2 心磁信号探测
心脏的电生理活动产生的磁场信号可以用SQUID器件进行探测. 采用多路探测系统, 可以获得胸部心磁信号的分布, 构成二维心磁图(magnetocardiography, MCG). 心磁信号的幅度约100 pT, 测量MCG一般需要0.1—1 pT的磁场分辨率.
由于心磁信号幅度比脑磁信号大约两个量级, 检测的技术难度相对较低. 人们预计心磁图可作为目前临床普遍使用的心电图的补充, 有望在一些致命性心脏疾病的早期诊断中发挥作用. 全球心脏疾病患者的数目非常庞大, 据统计, 我国冠心病患者约有1000万人, 而且每年有大量的新增病人. 现在临床广泛应用的心电图仍有局限性, 它不能反映所有心脏病变. 尤其是心电图对可能引起心肌梗死的冠心病的漏诊率约有40%, 因此迫切需要其他的快速方便的检测手段. 在这样的情况下, 心磁图有可能成为一种很好的预先诊断心脏猝死的补充手段.
目前国际上已有一些公司推出了商业化产品. 日本Tsukuba大学医院利用一台64通道的心磁系统, 在近十年的时间, 对10000多位病人进行了心磁图测量, 结果清楚表明心磁图可以给出比心电图更多的信息[38].
国内包括北京阜外医院、天津泰达国际心血管医院以及上海第六人民医院等也先后开展了相关临床研究. 北京大学以及中科院物理所分别研制了基于高温超导SQUID器件的心磁测量试验系统, 中科院上海微系统所研制了多通道的低温SQUID心磁系统.
目前限制SQUID心磁图系统广泛推广使用的主要原因一是建造和使用成本, 二是缺乏建立在大量数据系统分析基础上的具有明确医学意义的心磁图结果与各种心脏疾病的对应关系. 高温超导SQUID器件的灵敏度能较好地满足心磁测量, 主要的挑战是提高器件制备工艺的可靠性和良率, 以便用于多通道心磁测量系统.
5.2 其他生物医学医用
除上述的脑磁和心磁等生物磁测量应用外, SQUID弱磁检测在生物医学中还有一些其他应用. 例如, 用SQUID追踪检测包附有生物功能基团磁性纳米粒子的SQUID免疫测定[39], 和用SQUID作为探测器的超低磁场核磁共振和成像等[11]. 由于SQUID具有高灵敏度, SQUID免疫测定可检测到极少量的病原. 超低磁场磁共振成像中, 具有不同自旋晶格弛豫时间生物组织的成像对比度更大, 有可能便于区分病变组织. 因此, 这两方面的应用将有可能在诸如癌症和肿瘤组织的早期检测领域发挥作用.
包附有生物功能基团的磁性纳米粒子的追踪检测是生物医学上常用的方法之一. 免疫测定在生物和医学中广泛用于复杂试样如血液中的被分析物含量的确定. 它通过抗原与相应抗体之间的特定结合, 确定试样中是否存在某种特定的抗原及其含量. 随着纳米材料的发展, 利用磁性纳米粒子作为生物标记进行免疫测定的技术得到发展. 利用纳米磁性粒子作为标记并再利用SQUID所具有的非常高的磁测量灵敏度, 可以测定出已有方法不能测到的微量抗原. 磁性纳米粒子具有稳定性好、无毒、可以通过磁场进行操控的特点. 实验中通过测量磁性粒子磁矩的变化, 可以定量地确定待测试样中抗原的含量. 这种方法还避免了已有方法中需要洗去未结合抗体的步骤, 使得测定过程更加简单方便. 因此, 这种方法对于生物和医学研究以及更进一步的临床诊断都具有重要的意义.
检测的方式有弛豫时间、交流共振和剩余磁场检测等几种方式[39]. 有研究报道显示, 利用SQUID和有特定生物功能集团的磁性标记, 可以从血液中区分检测出早期奥兹海默症对应的Aβ1–42蛋白[40].
这方面的研究大多采用高温超导SQUID, 应用的推广需要研制针对各种适应检测的磁性标示物.
磁共振成像(MRI)是超导技术最广泛的商业应用, 已经成为临床医学检测和研究不可或缺的工具. 虽然目前磁共振成像技术发展的大趋势是更高磁场, 但在过去二十多年中, 人们也一直探索在超低磁场环境下的核磁共振谱和磁共振成像技术, 以期开发出新的应用. 与常规高场磁共振成像系统相比较, 超低磁场磁共振成像(ULF-MRI)可以获得更好的自旋晶格弛豫时间T1对比度, 这有可能在肿瘤以及其他检测中发挥作用[41—43]. 超低磁场核磁成像系统的磁场一般在100 μT的量级, 这样的系统简单轻便, 便于移动. 为了补偿由于磁场降低而导致的信号大幅降低, 一般采用施加预极化磁场和采用灵敏的SQUID探测器的方式来提升信噪比. 已有的研究显示, 在如此低的磁场下不但可以得到人体组织成像, 而且一系列肿瘤组织的体外测量也表明有可能区分肿瘤和正常组织[44].
对于实际应用而言, 最主要的挑战是如何减少系统中的各种噪声, 提升整体探测灵敏度. 预计, 基于SQUID探测器的超低磁场核磁共振和成像系统有可能在肿瘤组织、脑部受损组织、阿尔兹海默症的检测等方面获得应用.
由于ULF-MRI和MEG都采用SQUID器件作为探测器, 可以将两种检测技术集成起来. 这样, 将有可能集合MEG测量时间分辨性高和MRI空间分辨率好的优点, 有助于脑磁信号的定位和源的重构, 以及获得神经电流的分布图等[45]. 图10给出了在两个不同MEG-MRI集成实验系统上测量的结果. 其中图10(a)和图10(b)分别给出在一个MEG-MRI集成系统上和一个常规的先进的MEG上测量得到的脑磁信号结果, 数据是视觉诱发反应产生的等效偶极和磁场分布. 两组结果是在同样的视觉激励方式下得到的. 图10(c)和图10(d)则分别给出同一个测试样本在另一个MEG-MRI集成系统得到的96 μT磁场下的ULF-MRI切片图和一个常规3T的MRI切片图. 其中ULF-MRI系统也测量听觉反应产生的脑磁信号, 相应的偶极子的位置和方向在图10(c)中标出.
图 10 (a) 在MEG-MRI集成实验系统上测量的以及(b)同样刺激在最先进的MEG系统上得到的视觉诱发反应产生的等效偶极和磁场分布; (c)用另一个MEG-MRI集成实验系统在96 μT磁场下得到的超低磁场MRI切片图和记录的听觉反应脑磁信号对应的偶极子; (d)常规3 T磁场下同一个样本的切片图[11]Fig. 10. Equivalent dipoles and field patterns of the visually evoked responses using (a) the MEG–MRI system and (b) state-of-the-art MEG with the same stimulus protocol. MRI slices (c) at 96 μT, with the registered equivalent dipole of the auditory response overlaid, and (d) from an uncoregistered 3 T image acquired separately from the same subject[11].5.3 地球物理应用
SQUID器件的高灵敏性和很好的低频特征非常适合于地球物理方面的应用, 用以探测微弱的地球磁场变化或激励的响应信号. 自SQUID出现以来, 已经在多种场合开展过示范应用, 一些商业化设备也已经出现[46—49]. 这方面的应用大致可分为被动型和主动型两类. 在被动型探测中, SQUID器件磁强计或梯度计直接探测扫描区域内微弱的地磁场变化, 发现磁异常区域. 这些磁异常通常与各种金属矿或其他导电型矿有关. 这类应用中, 通常采用梯度计构型, 或者利用多个梯度计进行磁场分布的全张量测量.
在主动型探测中, 通过测量在地表发射的激励脉冲电磁信号的响应, 可以得到地表下不同深度的信息. SQUID器件的灵敏性和低频特性有助于获得常规感应式探测线圈所不能得到的深部信息.
当前我国经济快速发展, 对各种矿产资源的需求十分旺盛, 但是我国地下矿产探测主要局限在浅层(< 800 m), 远低于发达国家的水平(约2000 m), 深层资源的探测急需开展. 高灵敏度的SQUID磁探测系统对深层地质构造的磁探测能力独特, 远强于传统感应线圈式的探头, 具有较大优势, 将在探寻深层资源方面发挥重要作用. 国内相关单位近年来也加强了相关方面的研究, 在器件制备、系统集成和实地应用研究方面都取得了显著进展.
中国科学院上海微系统研究所成功研制出国际第二套航空低温超导全张量磁梯度探测系统[50], 并进行了多次野外飞行测量, 得到了航空超导全张量磁梯度分布图. 此外, 研制的用于瞬变电磁(TEM)探测超导SQUID系统, 噪声为7 fT, 动态范围达160 dB[51]. 在实际应用中, 多地实际反演深度大于2000 m. 先后成功探到铅锌矿, 并在云南、内蒙等常规技术疑难矿区获得突破, 探测结果得到钻井资料的验证. 图11是超导瞬变电磁探矿系统和在内蒙古四子王旗实地测量结果, 数据的反演结果中超过2000 m的低阻区域能够清楚地显示出来.
SQUID器件在地球物理探测中的应用技术还可应用于其他许多场合, 只要有局部的磁或磁响应异常, 包括军事上的潜艇探测、地下溶洞探测、发现考古遗址以及战争遗留的未爆炸武器等[52].
对于野外的现场测试, 使用液氮冷却的高温SQUID器件更加方便. 但是由于高温超导SQUID器件的可靠性差, 更容易受到环境干扰的影响, 比较成熟的系统中还是主要采用低温超导SQUID器件. 近年来, 受液氦供应短缺的影响, 低温SQUID地球物理探测方面的应用受到较大影响[53].
5.4 科学仪器中的应用
SQUID器件在各种科学仪器和科学研究领域中都有很多的应用. 目前广泛用于材料磁性测量和表征的美国Quantum Design公司的磁性质测量系统(magnetic properties measurement system, MPMS)产品是一种大家熟悉的使用SQUID的科学仪器. 该类系统中, 一个SQUID器件与超导线绕成的二阶梯度线圈构成的二阶梯度计用于测量样品的磁信号.
在各种大型天文望远镜中, 用于测量深空各种波段电磁辐射信号的探测器大多用到诸如超导转变边沿探测器(TES)、磁性微型量热器(MMC)、以及热电子热辐射器(HEB)等器件阵列, SQUID作为灵敏的探测器用于探测这些器件的微小电流变化[54]. 由于阵列中器件数目庞大, 需要大量SQUID读出电路, 人们发展了多种包括时分、码分和频分方式的多路复用方案[55]. 此外, 基于SQUID的低温射频放大器可具有接近量子极限的噪声水平, 在诸如暗物质探测等基础科学领域中也扮演着重要的角色, 可望将探测效率提升2—3个量级[56,57].
扫描SQUID显微镜是结合SQUID和扫描探测技术发展起来的高灵敏度微区磁信号探测仪器. 扫描SQUID主要应用于基础科学、纳米材料科学、工业无损检测和分子生物技术等研究领域. IBM实验室的研究人员利用扫描SQUID显微镜证实了高温超导体具有d波配对性, 大大推动了高温超导机理的研究[58]. 半导体芯片中纳米量级芯片导线的无损检测也是扫描SQUID显微镜的一个可能应用. 目前在市场上已有扫描SQUID显微镜的产品.
为了提高扫描SQUID的空间分辨率, 以色列Weizmann科学研究所的一个研究组发展了一种很有创意的称为SQUID-on-Tip的方案. 他们将石英管拉成直径约100 nm的空心细管, 在尖端做成针尖Nano-SQUID, 获得了低于100 nm的空间分辨率, 可测量单个电子自旋. 而且,他们利用SQUID临界电流随温度的变化, 实现了好于1 μK的局部温度测量. 他们用该系统研究了超导薄膜中磁通线的动力学行为(图12), 以及多种量子物态的局部耗散机制[59,60].
5.5 工业无损检测
无损检测在工业生产环节的质量控制、设备运行的可靠性保障、材料服役性能的评估等各个领域具有重要的地位. 随着人们对质量和可靠性要求的提高, 无损检测技术越来越受到重视, 新的检测方法和手段也在不断地研究和发展. 超导量子干涉器(SQUID)无损检测是研究比较多的一种检测技术. 虽然基于SQUID的无损检测系统价格较贵、操作不方便(附带低温冷却装置)会限制其使用范围, 但SQUID的高带宽、高磁场灵敏度特性, 在一些应用中会有很大的优势和发展潜力. 高温SQUID器件制备技术的发展和成熟将会促进SQUID在无损检测中的广泛应用.
SQUID无损检测实质上是探测被测物体附近的磁场分布从而非接触地获得与磁信号有关的信息(如缺陷大小、深度等等). 根据磁场产生的特点通常可分为无源法和有源法. 前者无需激励场, 主要用于探测非铁磁材料内的磁掺杂、检测铁磁材料的应力分布、探测铁磁体表面的磁通泄漏(MFL)等. 后者需要对被测材料施加电流或利用电磁激励在导体中产生涡流, 电流(或涡流)的分布与导体的导电特性有关, 如果导体中存在缺陷, 该缺陷将会严重影响导体的导电性能, 例如导体块中如果存在裂缝, 电流(或涡流)在导体中流动时会绕过裂缝而重新分布, 通过SQUID检测该电流(或涡流)产生的磁场, 就可得到有关缺陷的信息. 由于SQUID具有极高的探测灵敏度, 因此在检测材料的深层缺陷、细小缺陷等方面具有很好的应用前景[10]. 已有的应用示范包括: 航空器零部件检测、半导体掺杂均匀性检测、焊接工艺影响检测、碳纤维复合材料检测、金属材料的应力应变检测、磁泄漏(MFL)检测、金属材料腐蚀性检测、材料或食品中的磁性或金属粒子检测等许多方面[61-63].
6. 小结和展望
超导量子干涉器件作为一种最灵敏的磁探测器件, 自从出现以来, 一直受到广泛的关注. 在科学研究、生物医学、资源探测和工业检测等多个领域的应用研究已经取得显著进步, 在一些地方发挥着难以替代的作用. 需要低温环境, 尤其是液氦的使用, 是限制其应用推广的一个重要因素, 采用高温超导材料制备的器件可以在很大程度上减少这种限制的影响. 因此, 不断改进高温超导SQUID器件的制备技术, 提升器件性能及一致性是未来SQUID应用的一个关键. 同时, SQUID技术的发展也依赖于新的创新应用的开发. 进入二十一世纪以来, 对于脑科学的研究越来越受到重视, 基于SQUID的脑磁图测量有可能发挥更大的作用. 目前量子计算的热潮也促进和拓展了SQUID器件的应用, SQUID可以用于构建超导量子比特, 也可用于超导量子比特的直接读出, 还可以做成高灵敏度的单光子水平微波信号放大器来探测量子比特态的信息[64-67].
[1] London F 1950 Superfluids (New York: Wiley)
[2] Deaver B S, Fairbank W M 1961 Phys. Rev. Lett. 7 43
Google Scholar
[3] Doll R, Näbauer M 1961 Phys. Rev. Lett. 7 51
Google Scholar
[4] Gough C E, Colclough M S, Forgan E M, Jordan R G, Keene M, Muirhead C M, Rae A I M, Thomas N, Abell J S, Sutton S 1987 Nature 326 855
Google Scholar
[5] Josephson B D 1962 Phys. Lett. 1 251
Google Scholar
[6] Anderson P W, Rowell J M 1963 Phys. Rev. Lett. 10 230
Google Scholar
[7] Clarke J, Braginski A I 2004 The SQUID Handbook (Vol. 1) (New Jersey: John Wiley & Sons)
[8] Clarke J, Braginski A I 2004 The SQUID Handbook (Vol. 2) (New Jersey: John Wiley & Sons)
[9] Fagaly R 2006 Rev. Sci. Instrum 77 101101
Google Scholar
[10] Anders S, Blamire M G, Buchholz F Im, Crété D G, Cristianoe R, Febvref P, Fritzsch L, Herr A, Ilichev E, Kohlmann J, Kunert J, Meyer H G, Niemeyer J, Ortlepp T, Rogalla H, Schurig T, Siegel M, Stolz R,Tarte E, Brake H J M ter, Toepfer H, Villegier J C, Zagoskin A M, Zorin A B 2010 Physica C 470 2079
Google Scholar
[11] Körber R, Storm J H, Seton H 2016 Supercond. Sci. Technol. 29 113001
Google Scholar
[12] Clarke J, Lee Y H, Schneiderman J 2018 Supercond. Sci. Technol. 31 080201
Google Scholar
[13] Pulizzi F 2011 Nature Materials 10 262
Google Scholar
[14] Zhang Y, Dong H, Krause H J, Zhang G F, Xie X M 2020 SQUID Readout Electronics and Magnetometric Systems for Practical Applications (New York: Wiley)
[15] Jaklevic R C, Lambe J, Silver A H, Mercereau J E 1964 Phys. Rev. Lett. 12 159
[16] Caputo P, Oppenläder J, Hässler Ch, Tomes J, Friesch A, Träble T, Schopohl N 2004 Appl. Phys. Lett. 85 1389
Google Scholar
[17] Drung D 2003 Supercond. Sci. Technol. 16 1320
Google Scholar
[18] Dantsker E, Tanaka S, Clarke J 1997 Appl. Phys. Lett. 70 2037
Google Scholar
[19] Tesche C D, Clarke J 1977 J. Low Temp. Phys. 27 301
[20] Faley M I, Pratt K, Reineman R, Schurig D, Gott S, Atwood C G, Sarwinski R E, Paulson D N, Starr T N, Fagaly R L 2017 Supercond. Sci. Technol. 30 083001
Google Scholar
[21] Gurvitch M, Washington M, Higgins H 1983 Appl. Phys. Lett. 42 472
Google Scholar
[22] Koelle D, Kleiner R, Ludwig F, Dantsker E, Clarke J 1999 Rev. Mod. Phys. 71 631
Google Scholar
[23] Tafuri F, Kirtley J R 2005 Rep. Prog. Phys. 68 2573
Google Scholar
[24] Dimos D, Chaudhari P, Mannhart J, LeGoues F K 1988 Phys. Rev. Lett. 61 219
Google Scholar
[25] Du J, Lazar J Y, Lam S K H, Mitchell E E, Foley C P 2014 Supercond. Sci. Technol. 27 095005
Google Scholar
[26] Wakana H, Adachi S, Kamitani A, Nakayama K, Ishimaru Y, Tarutani Y, Tanabe K 2005 IEICE Trans. Electron. E88-C 208
Google Scholar
[27] Bergeal N, Lesueur J, Faini G, Aprili M, Contour J 2006 Appl. Phys. Lett. 89 112515
Google Scholar
[28] Trabaldo E, Ruffieux S, Andersson E, Arpaia R, Montemurro D, Schneiderman J F, Kalaboukhov A, Winkler D, Lombardi F, Bauch T 2020 Appl. Phys. Lett. 116 132601
Google Scholar
[29] Cybart S A, Cho E, Wong T, Wehlin B H, Ma M K, Huynh C, Dynes R 2015 Nat. Nanotechnol. 10 598
Google Scholar
[30] Cho E Y, Li H, LeFebvre J C, Zhou Y W, Dynes R, Cybart S A 2018 Appl. Phys. Lett. 113 162602
Google Scholar
[31] Hari R, Salmelin R 2012 NeuroImage 61 386
Google Scholar
[32] Hämäläinen M, Hari R, Ilmoniemi R J, Knuutila J, Lounasmaa O V 1993 Rev. Mod. Phys. 65 413
[33] Mäkelä J P, Forss N, Jääskeläinen J, Kirveskari E, Korvenoja A, Paetau R 2006 Neurosurgery 59 493
Google Scholar
[34] Lee Y H, Kwon H, Yu K K, Kim J M, Lee S K, Kim M Y, Kim K 2017 Supercond. Sci. Technol. 30 084003
Google Scholar
[35] Öisjöen F, Schneiderman J F, Figueras G A, Chukharkin M L, Kalabukhov A, Hedström A, Elam M, Winkler D 2012 Appl. Phys. Lett. 100 132601
Google Scholar
[36] Sander T H, Preusser J, Mhaskar R, Kitching J, Trahms L, Knappe S 2012 Opt. Express 3 981
Google Scholar
[37] Boto E, Holmes N, Leggett J, Roberts G, Shah V, Meyer S S, Muñoz L D, Mullinger K J, Tierney T M, Bestmann S, Barnes G R, Bowtell R, Brookes M J 2018 Nature 555 657
Google Scholar
[38] Inaba T, Nakazawa Y, Yoshida K, Kato Y, Hattori A, Kimura T, Hoshi T, Ishizu T, Seo Y, Sato A, Sekiguchi Y, Nogami A, Watanabe S, Horigome H, Kawakami Y, Aonuma K 2017 Supercond. Sci. Technol. 30 114003
Google Scholar
[39] Enpuku K, Tsujita Y, Nakamura K, Sasayama T, Yoshida T 2017 Supercond. Sci. Technol. 30 053002
Google Scholar
[40] Yang S Y, Chieh J J, Yang C C, Liao S H, Chen H H, Horng H E, Yang H C, Hong C Y, Chiu M J, Chen T F, Huang K W, Wu C C 2013 IEEE Trans. Appl. Supercond. 23 1600604
Google Scholar
[41] Clarke J, Hatridge M, Möβle M 2007 Annu. Rev. Biomed. Eng. 9 389
Google Scholar
[42] Dong H, Hwang S M, Wendland M, You L, Clarke J, Inglis B 2017 Magn. Reson. Med. 78 2342
Google Scholar
[43] Buckenmaier K, Pedersen A, SanGiorgio P, Scheffler K, Clarke J, Inglis B 2019 Neuroimage. 186 185
Google Scholar
[44] Busch S E, Hatridge M, Mößle M, Myers W, Wong T, Mück M, Chew K, Kuchinsky K, Simko J, Clarke J 2012 Magn. Reson. Med. 67 1138
Google Scholar
[45] Magnelind P E, Gomez J J, Matlashov A N, Owens T, Sandin J H, Volegov P L, Espy M A 2011 IEEE Trans. Appl. Supercond. 21 456
Google Scholar
[46] Leslie K E, Binks R A, Lam S K H, Sullivan P A, Tilbrook D L, Thorn R G, Foley C P 2008 Leading Edge 27 70
Google Scholar
[47] Hato T, Tsukamoto A, Adachi S, Oshikubo Y, Watanabe H, Ishikawa H, Sugisaki M, Arai E, Tanabe K 2013 Supercond. Sci. Technol. 26 115003
Google Scholar
[48] Chwala A, Stolz R, Schmelz M, Zakosarenko V, Meyer M, Meyer H G 2015 IEICE Trans. Electron. E98.C 167
Google Scholar
[49] Supracon A Ghttp://www.supracon.de/
[50] Song Z, Dai H, Rong L, et al. 2019 IEEE Trans. Appl. Supercond. 29 1600205
[51] Rong L, Bao S, Wu J, et al. 2019 IEEE Trans. Appl. Supercond. 29 1601704
[52] Stolz R, Zakosarenko V, Schulz M, Chwala A, Fritzsch L, Meyer H G, Kötlin E O 2006 Leading Edge 25 178
Google Scholar
[53] Halperin W P 2014 Nat. Phys. 10 467
Google Scholar
[54] Yang C, Si M T, You L X 2020 Sci. China Inf. Sci. 63 180502
Google Scholar
[55] Irwin K D, Cho H M, Doriese W B, Fowler J W, Hilton G C, Niemack M D, Reintsema C D, Schmidt D R, Ullom J N, Vale L R 2012 J. Low Temp. Phys. 167 588
Google Scholar
[56] Sloan J V, Hotz M, Boutan C, et al. 2016 Phys. Dark Universe 14 95
Google Scholar
[57] Asztalos S J, Carosi G, Hagmann C, et al. 2010 Phys. Rev. Lett. 104 041301
Google Scholar
[58] Tsuei C C, Kirtley J R 2000 Rev. Mod. Phys. 72 969
Google Scholar
[59] Embon L, Anahory Y, Jelić Ž L, Lachman E O, Myasoedov Y, Huber M E, Mikitik G P, Silhanek A V, Milošević M V, Gurevich A, Zeldov E 2017 Nat. Commun. 8 85
Google Scholar
[60] Halbertal D, Cuppens J, M. Ben Shalom Shalom M. B, Embon L, Shadmi N, Anahory Y, Naren H R, Sarkar J, Uri A, Ronen Y, Myasoedov Y, Levitov L S, Joselevich E, Geim A K, Zeldov E 2016 Nature 539 407
Google Scholar
[61] Hatsukade Y, Inaba T, Maruno Y, Tanaka S 2005 IEEE Trans. Appl. Supercond. 15 723
Google Scholar
[62] Krause H J, Michael M M, Tanaka S 2015 SQUlDs in Nondestructive Evaluation. In Applied Superconductivity: Handbook on Devices and Applications (Vol. 2) (Weinheim: Wiley)
[63] Faley M, Kostyurina E, Kalashnikov K, Maslennikov Y, Koshelets V, Dunin-Borkowski R 2017 Sensors 17 2798
Google Scholar
[64] Clarke J, Wilhelm F K 2008 Nature 453 1031
Google Scholar
[65] Krantz J P, Kjaergaard M, Yan F, Orlando T P, Gustavsson S, Oliver W D 2019 Appl. Phys. Rev. 6 021318
Google Scholar
[66] Kwon S, Tomonaga A, Gopika L B, Devitt S J, Tsai J S 2020 arXiv:2009.08021 [quant-ph]
[67] Kjaergaard M, Schwartz M E, Braumuller J, Krantz P, Wang J I J, Gustavsson S, Oliver W D 2020 Annu. Rev. Condens. Matter Phys. 11 369
Google Scholar
期刊类型引用(18)
1. 陈佳钰,赵凯明,张斌珍. SQUID传感器测量的微弱磁场异常信号检测方法. 微纳电子技术. 2025(02): 148-155 . 百度学术
2. 王学锋,刘崇泰,卢向东,李建军,邓意成,徐强锋. 基于CPT效应矢量磁场测量技术发展现状与展望. 激光与光电子学进展. 2025(01): 125-133 . 百度学术
3. 刘斌,谢明玲,员朝鑫,强进,王向谦. 隧穿磁电阻效应磁场传感器中的噪声分析及降低噪声的方法. 甘肃科技. 2024(03): 1-4+9 . 百度学术
4. Wenzhu Wu,Ming Xu,Changlin Han,Junquan Tang,Jia Xu,Lin Xu,Mingxin Qin. Design and Implementation of a High-Sensitivity Magnetic Sensing System Based on GMI Effect. Journal of Beijing Institute of Technology. 2024(03): 237-247 . 必应学术
5. 朱承峰,李万,陈建,孙潇莹,曾九孙,蔡晋辉,王振宇,李劲劲,王雪深. 柔性超导电路制备工艺技术研究. 计量学报. 2024(07): 929-934 . 百度学术
6. 王睿奇,封燮,张斌珍,汤苏晋. 基于改良Pan-Tompkins算法的心磁R波定位方法. 舰船电子工程. 2024(07): 29-34 . 百度学术
7. 耿佳琪,赵龙,仇如嘉,柯艳国,赵博文. 金刚石NV色心磁力仪测磁性能优化. 电子测量技术. 2024(17): 115-122 . 百度学术
8. 李青竹,李晶,李志宁,石志勇,文雪忠. 磁梯度张量系统目标探测技术研究进展. 兵工学报. 2024(12): 4205-4230 . 百度学术
9. 汤苏晋,封燮,王睿奇,邢彪,段俊萍,张斌珍. 基于迈斯纳效应的二阶梯度计不平衡度调节方法. 测试技术学报. 2023(04): 316-322 . 百度学术
10. 王睿奇,封燮,汤苏晋,陈铭远,张斌珍,段俊萍. 应用于无屏蔽环境的心磁测量系统研究与设计. 传感技术学报. 2023(07): 1153-1158 . 百度学术
11. 汤苏晋,封燮,王睿奇,陈铭远,张斌珍,段俊萍. 基于SQUID梯度计的心磁图仪标定及通道差异性研究. 传感技术学报. 2023(11): 1740-1745 . 百度学术
12. 尹渌,张朝晖,张天尧,米旺,帕夫洛·萨特科夫斯基. SQUID二阶磁梯度计标定系统的空间位置优化. 仪器仪表学报. 2023(09): 155-164 . 百度学术
13. 郑文,于扬. 超导量子计算核心器件. 物理. 2023(11): 731-743 . 百度学术
14. 秦晟昊,张昊晖,吴瑶,周晗,黄一聪,林建新. 基于双晶退火衬底的高温超导纳米量子干涉仪在液氦液氮温区下的输运特性. 低温物理学报. 2023(06): 330-336 . 百度学术
15. 丁立,鲁云峰,赵建亭,鄢志丹,贺青. 基于低温电流比较仪的微弱电流放大器设计. 仪器仪表学报. 2022(05): 86-92 . 百度学术
16. 周涛,陆颖,孔祥燕,史继强. 趋势观察:超导量子干涉器件相关研究发展态势与启示. 中国科学院院刊. 2022(11): 1666-1670 . 百度学术
17. 张博,贺青,杨欣达,欧阳鹏辉,王轶文,韦联福. 共面波导型超导微波功分器:设计、制备和测试. 物理学报. 2021(15): 325-331 . 百度学术
18. 韩金舸,欧阳鹏辉,李恩平,王轶文,韦联福. 超导约瑟夫森结物理参数的实验推算. 物理学报. 2021(17): 31-37 . 百度学术
其他类型引用(26)
-
图 3 DC-SQUID (a)等效电路示意图; (b)磁通分别为整数个和半整数个Φ0时的I-V曲线; (c)电压-磁通曲线[7]
Figure 3. The DC-SQUID: (a) Schematic electric circuit; (b) current-voltage characteristics at integer and half-integer values of applied flux; the operation point is set by the bias current Ib; (c) voltage vs. flux Φa/Φ0 for constant bias current[7].
图 7 MgO衬底上45°台阶上生长的YBCO薄膜的扫描电镜(SEM)图像(左图)和高分辨透射电镜HRTEM图像(中图). 右图: 一个16 mm大小、采用台阶结的高温超导DC-SQUID磁强计在超导屏蔽环境下测量的噪声谱[20]
Figure 7. SEM image (left) and HRTEM image (middle) of an YBCO film deposited on a double-layer-buffered 45° step on an MgO substrate. A 45° [100]-tilted GB is clearly shown. Right: Noise spectral density of a 16 mm high-Tc DC-SQUID magnetometer with step-edge junctions measured in a superconducting shield[20].
图 10 (a) 在MEG-MRI集成实验系统上测量的以及(b)同样刺激在最先进的MEG系统上得到的视觉诱发反应产生的等效偶极和磁场分布; (c)用另一个MEG-MRI集成实验系统在96 μT磁场下得到的超低磁场MRI切片图和记录的听觉反应脑磁信号对应的偶极子; (d)常规3 T磁场下同一个样本的切片图[11]
Figure 10. Equivalent dipoles and field patterns of the visually evoked responses using (a) the MEG–MRI system and (b) state-of-the-art MEG with the same stimulus protocol. MRI slices (c) at 96 μT, with the registered equivalent dipole of the auditory response overlaid, and (d) from an uncoregistered 3 T image acquired separately from the same subject[11].
-
[1] London F 1950 Superfluids (New York: Wiley)
[2] Deaver B S, Fairbank W M 1961 Phys. Rev. Lett. 7 43
Google Scholar
[3] Doll R, Näbauer M 1961 Phys. Rev. Lett. 7 51
Google Scholar
[4] Gough C E, Colclough M S, Forgan E M, Jordan R G, Keene M, Muirhead C M, Rae A I M, Thomas N, Abell J S, Sutton S 1987 Nature 326 855
Google Scholar
[5] Josephson B D 1962 Phys. Lett. 1 251
Google Scholar
[6] Anderson P W, Rowell J M 1963 Phys. Rev. Lett. 10 230
Google Scholar
[7] Clarke J, Braginski A I 2004 The SQUID Handbook (Vol. 1) (New Jersey: John Wiley & Sons)
[8] Clarke J, Braginski A I 2004 The SQUID Handbook (Vol. 2) (New Jersey: John Wiley & Sons)
[9] Fagaly R 2006 Rev. Sci. Instrum 77 101101
Google Scholar
[10] Anders S, Blamire M G, Buchholz F Im, Crété D G, Cristianoe R, Febvref P, Fritzsch L, Herr A, Ilichev E, Kohlmann J, Kunert J, Meyer H G, Niemeyer J, Ortlepp T, Rogalla H, Schurig T, Siegel M, Stolz R,Tarte E, Brake H J M ter, Toepfer H, Villegier J C, Zagoskin A M, Zorin A B 2010 Physica C 470 2079
Google Scholar
[11] Körber R, Storm J H, Seton H 2016 Supercond. Sci. Technol. 29 113001
Google Scholar
[12] Clarke J, Lee Y H, Schneiderman J 2018 Supercond. Sci. Technol. 31 080201
Google Scholar
[13] Pulizzi F 2011 Nature Materials 10 262
Google Scholar
[14] Zhang Y, Dong H, Krause H J, Zhang G F, Xie X M 2020 SQUID Readout Electronics and Magnetometric Systems for Practical Applications (New York: Wiley)
[15] Jaklevic R C, Lambe J, Silver A H, Mercereau J E 1964 Phys. Rev. Lett. 12 159
[16] Caputo P, Oppenläder J, Hässler Ch, Tomes J, Friesch A, Träble T, Schopohl N 2004 Appl. Phys. Lett. 85 1389
Google Scholar
[17] Drung D 2003 Supercond. Sci. Technol. 16 1320
Google Scholar
[18] Dantsker E, Tanaka S, Clarke J 1997 Appl. Phys. Lett. 70 2037
Google Scholar
[19] Tesche C D, Clarke J 1977 J. Low Temp. Phys. 27 301
[20] Faley M I, Pratt K, Reineman R, Schurig D, Gott S, Atwood C G, Sarwinski R E, Paulson D N, Starr T N, Fagaly R L 2017 Supercond. Sci. Technol. 30 083001
Google Scholar
[21] Gurvitch M, Washington M, Higgins H 1983 Appl. Phys. Lett. 42 472
Google Scholar
[22] Koelle D, Kleiner R, Ludwig F, Dantsker E, Clarke J 1999 Rev. Mod. Phys. 71 631
Google Scholar
[23] Tafuri F, Kirtley J R 2005 Rep. Prog. Phys. 68 2573
Google Scholar
[24] Dimos D, Chaudhari P, Mannhart J, LeGoues F K 1988 Phys. Rev. Lett. 61 219
Google Scholar
[25] Du J, Lazar J Y, Lam S K H, Mitchell E E, Foley C P 2014 Supercond. Sci. Technol. 27 095005
Google Scholar
[26] Wakana H, Adachi S, Kamitani A, Nakayama K, Ishimaru Y, Tarutani Y, Tanabe K 2005 IEICE Trans. Electron. E88-C 208
Google Scholar
[27] Bergeal N, Lesueur J, Faini G, Aprili M, Contour J 2006 Appl. Phys. Lett. 89 112515
Google Scholar
[28] Trabaldo E, Ruffieux S, Andersson E, Arpaia R, Montemurro D, Schneiderman J F, Kalaboukhov A, Winkler D, Lombardi F, Bauch T 2020 Appl. Phys. Lett. 116 132601
Google Scholar
[29] Cybart S A, Cho E, Wong T, Wehlin B H, Ma M K, Huynh C, Dynes R 2015 Nat. Nanotechnol. 10 598
Google Scholar
[30] Cho E Y, Li H, LeFebvre J C, Zhou Y W, Dynes R, Cybart S A 2018 Appl. Phys. Lett. 113 162602
Google Scholar
[31] Hari R, Salmelin R 2012 NeuroImage 61 386
Google Scholar
[32] Hämäläinen M, Hari R, Ilmoniemi R J, Knuutila J, Lounasmaa O V 1993 Rev. Mod. Phys. 65 413
[33] Mäkelä J P, Forss N, Jääskeläinen J, Kirveskari E, Korvenoja A, Paetau R 2006 Neurosurgery 59 493
Google Scholar
[34] Lee Y H, Kwon H, Yu K K, Kim J M, Lee S K, Kim M Y, Kim K 2017 Supercond. Sci. Technol. 30 084003
Google Scholar
[35] Öisjöen F, Schneiderman J F, Figueras G A, Chukharkin M L, Kalabukhov A, Hedström A, Elam M, Winkler D 2012 Appl. Phys. Lett. 100 132601
Google Scholar
[36] Sander T H, Preusser J, Mhaskar R, Kitching J, Trahms L, Knappe S 2012 Opt. Express 3 981
Google Scholar
[37] Boto E, Holmes N, Leggett J, Roberts G, Shah V, Meyer S S, Muñoz L D, Mullinger K J, Tierney T M, Bestmann S, Barnes G R, Bowtell R, Brookes M J 2018 Nature 555 657
Google Scholar
[38] Inaba T, Nakazawa Y, Yoshida K, Kato Y, Hattori A, Kimura T, Hoshi T, Ishizu T, Seo Y, Sato A, Sekiguchi Y, Nogami A, Watanabe S, Horigome H, Kawakami Y, Aonuma K 2017 Supercond. Sci. Technol. 30 114003
Google Scholar
[39] Enpuku K, Tsujita Y, Nakamura K, Sasayama T, Yoshida T 2017 Supercond. Sci. Technol. 30 053002
Google Scholar
[40] Yang S Y, Chieh J J, Yang C C, Liao S H, Chen H H, Horng H E, Yang H C, Hong C Y, Chiu M J, Chen T F, Huang K W, Wu C C 2013 IEEE Trans. Appl. Supercond. 23 1600604
Google Scholar
[41] Clarke J, Hatridge M, Möβle M 2007 Annu. Rev. Biomed. Eng. 9 389
Google Scholar
[42] Dong H, Hwang S M, Wendland M, You L, Clarke J, Inglis B 2017 Magn. Reson. Med. 78 2342
Google Scholar
[43] Buckenmaier K, Pedersen A, SanGiorgio P, Scheffler K, Clarke J, Inglis B 2019 Neuroimage. 186 185
Google Scholar
[44] Busch S E, Hatridge M, Mößle M, Myers W, Wong T, Mück M, Chew K, Kuchinsky K, Simko J, Clarke J 2012 Magn. Reson. Med. 67 1138
Google Scholar
[45] Magnelind P E, Gomez J J, Matlashov A N, Owens T, Sandin J H, Volegov P L, Espy M A 2011 IEEE Trans. Appl. Supercond. 21 456
Google Scholar
[46] Leslie K E, Binks R A, Lam S K H, Sullivan P A, Tilbrook D L, Thorn R G, Foley C P 2008 Leading Edge 27 70
Google Scholar
[47] Hato T, Tsukamoto A, Adachi S, Oshikubo Y, Watanabe H, Ishikawa H, Sugisaki M, Arai E, Tanabe K 2013 Supercond. Sci. Technol. 26 115003
Google Scholar
[48] Chwala A, Stolz R, Schmelz M, Zakosarenko V, Meyer M, Meyer H G 2015 IEICE Trans. Electron. E98.C 167
Google Scholar
[49] Supracon A Ghttp://www.supracon.de/
[50] Song Z, Dai H, Rong L, et al. 2019 IEEE Trans. Appl. Supercond. 29 1600205
[51] Rong L, Bao S, Wu J, et al. 2019 IEEE Trans. Appl. Supercond. 29 1601704
[52] Stolz R, Zakosarenko V, Schulz M, Chwala A, Fritzsch L, Meyer H G, Kötlin E O 2006 Leading Edge 25 178
Google Scholar
[53] Halperin W P 2014 Nat. Phys. 10 467
Google Scholar
[54] Yang C, Si M T, You L X 2020 Sci. China Inf. Sci. 63 180502
Google Scholar
[55] Irwin K D, Cho H M, Doriese W B, Fowler J W, Hilton G C, Niemack M D, Reintsema C D, Schmidt D R, Ullom J N, Vale L R 2012 J. Low Temp. Phys. 167 588
Google Scholar
[56] Sloan J V, Hotz M, Boutan C, et al. 2016 Phys. Dark Universe 14 95
Google Scholar
[57] Asztalos S J, Carosi G, Hagmann C, et al. 2010 Phys. Rev. Lett. 104 041301
Google Scholar
[58] Tsuei C C, Kirtley J R 2000 Rev. Mod. Phys. 72 969
Google Scholar
[59] Embon L, Anahory Y, Jelić Ž L, Lachman E O, Myasoedov Y, Huber M E, Mikitik G P, Silhanek A V, Milošević M V, Gurevich A, Zeldov E 2017 Nat. Commun. 8 85
Google Scholar
[60] Halbertal D, Cuppens J, M. Ben Shalom Shalom M. B, Embon L, Shadmi N, Anahory Y, Naren H R, Sarkar J, Uri A, Ronen Y, Myasoedov Y, Levitov L S, Joselevich E, Geim A K, Zeldov E 2016 Nature 539 407
Google Scholar
[61] Hatsukade Y, Inaba T, Maruno Y, Tanaka S 2005 IEEE Trans. Appl. Supercond. 15 723
Google Scholar
[62] Krause H J, Michael M M, Tanaka S 2015 SQUlDs in Nondestructive Evaluation. In Applied Superconductivity: Handbook on Devices and Applications (Vol. 2) (Weinheim: Wiley)
[63] Faley M, Kostyurina E, Kalashnikov K, Maslennikov Y, Koshelets V, Dunin-Borkowski R 2017 Sensors 17 2798
Google Scholar
[64] Clarke J, Wilhelm F K 2008 Nature 453 1031
Google Scholar
[65] Krantz J P, Kjaergaard M, Yan F, Orlando T P, Gustavsson S, Oliver W D 2019 Appl. Phys. Rev. 6 021318
Google Scholar
[66] Kwon S, Tomonaga A, Gopika L B, Devitt S J, Tsai J S 2020 arXiv:2009.08021 [quant-ph]
[67] Kjaergaard M, Schwartz M E, Braumuller J, Krantz P, Wang J I J, Gustavsson S, Oliver W D 2020 Annu. Rev. Condens. Matter Phys. 11 369
Google Scholar
期刊类型引用(18)
1. 陈佳钰,赵凯明,张斌珍. SQUID传感器测量的微弱磁场异常信号检测方法. 微纳电子技术. 2025(02): 148-155 . 百度学术
2. 王学锋,刘崇泰,卢向东,李建军,邓意成,徐强锋. 基于CPT效应矢量磁场测量技术发展现状与展望. 激光与光电子学进展. 2025(01): 125-133 . 百度学术
3. 刘斌,谢明玲,员朝鑫,强进,王向谦. 隧穿磁电阻效应磁场传感器中的噪声分析及降低噪声的方法. 甘肃科技. 2024(03): 1-4+9 . 百度学术
4. Wenzhu Wu,Ming Xu,Changlin Han,Junquan Tang,Jia Xu,Lin Xu,Mingxin Qin. Design and Implementation of a High-Sensitivity Magnetic Sensing System Based on GMI Effect. Journal of Beijing Institute of Technology. 2024(03): 237-247 . 必应学术
5. 朱承峰,李万,陈建,孙潇莹,曾九孙,蔡晋辉,王振宇,李劲劲,王雪深. 柔性超导电路制备工艺技术研究. 计量学报. 2024(07): 929-934 . 百度学术
6. 王睿奇,封燮,张斌珍,汤苏晋. 基于改良Pan-Tompkins算法的心磁R波定位方法. 舰船电子工程. 2024(07): 29-34 . 百度学术
7. 耿佳琪,赵龙,仇如嘉,柯艳国,赵博文. 金刚石NV色心磁力仪测磁性能优化. 电子测量技术. 2024(17): 115-122 . 百度学术
8. 李青竹,李晶,李志宁,石志勇,文雪忠. 磁梯度张量系统目标探测技术研究进展. 兵工学报. 2024(12): 4205-4230 . 百度学术
9. 汤苏晋,封燮,王睿奇,邢彪,段俊萍,张斌珍. 基于迈斯纳效应的二阶梯度计不平衡度调节方法. 测试技术学报. 2023(04): 316-322 . 百度学术
10. 王睿奇,封燮,汤苏晋,陈铭远,张斌珍,段俊萍. 应用于无屏蔽环境的心磁测量系统研究与设计. 传感技术学报. 2023(07): 1153-1158 . 百度学术
11. 汤苏晋,封燮,王睿奇,陈铭远,张斌珍,段俊萍. 基于SQUID梯度计的心磁图仪标定及通道差异性研究. 传感技术学报. 2023(11): 1740-1745 . 百度学术
12. 尹渌,张朝晖,张天尧,米旺,帕夫洛·萨特科夫斯基. SQUID二阶磁梯度计标定系统的空间位置优化. 仪器仪表学报. 2023(09): 155-164 . 百度学术
13. 郑文,于扬. 超导量子计算核心器件. 物理. 2023(11): 731-743 . 百度学术
14. 秦晟昊,张昊晖,吴瑶,周晗,黄一聪,林建新. 基于双晶退火衬底的高温超导纳米量子干涉仪在液氦液氮温区下的输运特性. 低温物理学报. 2023(06): 330-336 . 百度学术
15. 丁立,鲁云峰,赵建亭,鄢志丹,贺青. 基于低温电流比较仪的微弱电流放大器设计. 仪器仪表学报. 2022(05): 86-92 . 百度学术
16. 周涛,陆颖,孔祥燕,史继强. 趋势观察:超导量子干涉器件相关研究发展态势与启示. 中国科学院院刊. 2022(11): 1666-1670 . 百度学术
17. 张博,贺青,杨欣达,欧阳鹏辉,王轶文,韦联福. 共面波导型超导微波功分器:设计、制备和测试. 物理学报. 2021(15): 325-331 . 百度学术
18. 韩金舸,欧阳鹏辉,李恩平,王轶文,韦联福. 超导约瑟夫森结物理参数的实验推算. 物理学报. 2021(17): 31-37 . 百度学术
其他类型引用(26)
Catalog
Metrics
- Abstract views: 23344
- PDF Downloads: 1835
- Cited By: 44