关联成像

**量子成像**（英語：Quantum Imaging，QI），又稱**鬼成像**（英語：Ghost Imaging，GI）、**關聯成像**（英語：Correlated imaging，CI）等，是一種利用光場的二階或高階關聯獲得物體資訊的成像方法。量子成像屬於非定域成像，其概念起源於 20 世紀 50 年代的 HB-T 實驗。繼纏結光量子成像實驗之後，陸續有研究者提出了古典光量子成像、無透鏡量子成像、計算量子成像、差分量子成像等技術。量子成像技術在光刻、雷射雷達、生物組織造影、水下成像等領域都有應用。

## 歷史

### 基礎研究

20 世紀 50 年代，漢伯里・布朗（英語：R. Hanbury Brown）與特威斯（英語：R. Q. Twiss）進行了 HB-T 實驗，其主要內容是設計了一種干涉儀以解決在基線較長的情況下，無線電星角直徑的測量問題。這種干涉儀在無線電天文學中又被稱為強度干涉儀，其大致結構是：來自光源的光被半鍍銀鏡分為兩束，分別照射兩個光電倍增管的陰極，兩管的輸出電流經過放大以後在線性混頻器中相乘，如此持續運行大約一小時以獲得可觀測的結果。因此其在形態上類似邁克生干涉儀，但不同點在於，邁克生干涉儀是在檢測之前混合兩路信號，並且信號在混合之前始終擁有相位資訊。這種干涉儀則不然，進行混合的只是兩路信號的強度，相位資訊在經過光電倍增管之後已經被抹去了。但即便如此，兩路信號仍然表現出相干性，這意味著光場的相干性不僅在相位上體現出來，甚至僅僅在強度上也可以體現。文章發表之後，一些物理學家無法在 HB-T 實驗驗證中得出顯著的相干性，並聲稱「這種相干性的存在要求對量子力學中一些基本概念進行重大修訂」。漢伯里・布朗和特威斯進行了理論推導，並嘗試了對天狼星 A 的觀測，最終驗證了實驗結果的合理性。這個實驗將光子符合探測引入光學實驗，使得人們認識到光的二階干涉效應會造成光場的強度關聯（相關性）。

1970 年，美國國家宇航局電子研究中心的大衛・伯納姆（英語：David Burnham）與唐納德・溫伯格（英語：Donald Weinberg）在光子計數實驗中使用了名叫「自發參量下轉換」的非線性光學技術。兩人用氦 - 鎘雷射器產生 325nm 雷射束泵浦 25mm 長的磷酸二氫銨（ADP）晶體，測得入射光子轉化為相位匹配（入射光子動量和能量等於兩出射光子動量和與能量和）的雙光子的機率最大。這一技術使得纏結光子對可以很方便地獲得。1985 年以後，人們開始注意到自發參量下雙光子場的存在，並開始對其進行研究和應用。1988 年，蘇聯學者大衛・尼古拉耶維奇提出了一種用來驗證愛因斯坦 - 波多爾斯基 - 羅森實驗，並證明互補原理的方法，其中設計的一種裝置實際上可以作為量子成像設備使用。這種裝置進行符合檢測的方法是，用兩路探測器中的一路是否收到光子，來控制另一路探測的開閉。大衛・尼古拉耶維奇也由此被認為是量子成像方案的提出者。

### 纏結光量子成像

1994 年，大衛・尼古拉耶維奇在一篇嘗試同時使用古典表述和量子表述解釋光的量子性質的論文中，從量子力學的角度闡釋了雙光子非定域成像的理論基礎。1995 年，美國馬里蘭大學的史硯華小組以纏結雙光子作為光源，結合符合測量技術，實現了量子成像。史硯華小組用波長 351.1nm 的氬離子雷射器產生直徑 2mm 的光束作為泵浦源，泵浦被切割成為簡併 II 型相位匹配角的 ${\ce {\beta-BaB2O4}}$ （偏硼酸鋇，BBO）晶體，出射兩路正交偏振信號。分別將從 BBO 晶體的 e - 射線平面和 o - 射線平面出射的光束稱為信號光和閒置光，這兩路光的波長均為 702.2nm，即泵浦光波長的二倍。混合有未經自發參量下轉換的原泵浦光、信號光和閒置光的混合光束經過紫外級熔融石英色散稜鏡濾除剩餘的波長為 351.1nm 的原泵浦光，進入作為偏振分束器的湯普森稜鏡將信號光和閒置光分開。信號光束照射焦距為 400mm 的成像凸透鏡，透過的光束經過一個簡併波長為 702.2nm、帶寬為 83nm 的濾光片，再照射在物體（寫有「UMBC」字樣的光圈）上。經過物體之後，信號光束再通過一個焦距為 25mm 的收集凸透鏡，到達一個直徑為 0.8mm 的乾冰冷卻雪崩光電二極體 D1，這個雪崩光電二極體恰位於收集透鏡的焦點上。閒置光束經過另一個同樣的濾光片，不經過物體，直接照射在由 0.5mm 多模光纖尖端上，再到達另一個乾冰冷卻雪崩二極體 D2。兩個正交的編碼驅動檢測器在垂直於光束方向的平面上進行掃描，並將檢測器的輸出脈衝發送到接收窗口為 1.8ns 的符合計數器。設成像透鏡到收集透鏡距離為 S，成像透鏡到多模光纖尖端距離為 S'，成像透鏡焦距為 f，當三者滿足關係式

時，裝置可產生清晰物像。其中符合計數器的原理是，兩探測器均有響應時才計數，如果只有一個有響應或都沒有響應則不計數。探測器 D1 沒有空間解析度，D2 所在的閒置光（參考光）光路中沒有待成像物體，所以兩者都不能單獨完成物體成像；但處於無物體參考光路的 D2 探測器經過逐點掃描再與 D1 結果關聯之後得到了物體的像，使得量子成像又被稱為「鬼」成像。1997 年，索托・里貝羅（英語：P. H. Souto Ribeiro）小組發現，當自發參量下轉換晶體足夠薄（在泵浦光傳播方向上足夠短）時，泵浦光束的角頻譜將轉化為自發參量下轉換過程生成的雙光子狀態。小組通過將自發參量下轉換光束的橫向特性與泵浦光束的橫向特性聯繫起來，提出了控制信號光子與閒置光子之間橫向空間相關性的方法，實現了對纏結光子關聯性質的調製，使得在不向信號光路和閒置光路添加任何光學器件的情況下，也可以控制符合檢測入射光束的橫向特性，而不對強度產生影響。

### 古典光量子成像

史硯華小組實驗驗證量子成像之後，古典光是否能實現量子成像在學術界產生了爭議。2001 年，波士頓大學的艾曼・阿布瓦迪（英語：Ayman F. Abouraddy）小組分別用古典光和纏結光進行了量子成像實驗，之後發表評論稱，纏結光可以成部分相干像或者甚至完全相干像，古典光只能成非相干像；量子成像是纏結光的特性，其他雙光子源並不能模仿。然而，2002 年，羅切斯特大學的瑞安・本寧克（英語：Ryan S. Bennink）團隊就使用隨機掃描雷射光源證明了了古典光源也可以進行量子成像。Bennink 團隊對連續雷射束進行斬波，通過隨機旋轉反射鏡產生不同方向的偏轉光，再通過分束器來產生古典光源；用桶探測器檢測信號光是否未被測試圖案遮擋，並記錄信號光未被遮擋時 CCD 拍攝參考光記錄下的幀，即用桶探測器門控 CCD。實驗結果是，未門控時 CCD 無法拍攝到圖案，經過門控之後 CCD 可以拍攝到測試圖案。這證明了，儘管古典光源沒有纏結光源所特有的一些全局性質，但具有一定相關性的古典光源仍然可以通過符合計數產生量子成像現象。

2004 年，義大利的路易吉・盧吉亞托小組提出了熱光源的量子成像理論方案。在這種方案中，來自熱光源的光束被分束器分為兩束，後續處理與纏結光源成像相同，小組通過類比纏結光源和熱光源的成像過程，推斷可以使用純粹的非相干光源（熱光源）實現量子成像。2005 年，史硯華小組用氦氖雷射器產生的雷射入射旋轉毛玻璃盤產生的贗熱光源為光源，發現了非纏結光源的雙光子成像高斯薄透鏡方程式，實現了量子成像。同年，中科院物理所的吳令安小組實現了真熱光的雙光子二階關聯，不久後用 20mA 直流供電、諧振波長 780nm 的空心陰極燈（銣燈）實現了真熱光源量子成像。

儘管已經證明無論古典光源還是纏結光源都可以進行量子成像，但費里（F. Ferri）小組於 2005 年發表文章稱，成像的解析度具有一個上限，這個上限只有使用纏結光源可以達到。這證明，纏結光源相對古典光源量子成像，具有資訊可見性和解析度更高的優勢，這種優勢在高精度測量和量子通信領域都有應用。另一方面，熱光源更容易生成和測量，但成像的解析度會更低，背景雜訊更強。

2008 年，美國麻省理工學院教授傑弗里・夏皮羅使用高斯態光模型理論對量子成像原理進行了統一的解釋。

### 無透鏡量子成像

2004 年，韓申生和程靜提出可以通過適當選擇成像幾何形狀，用非相干光源實現無透鏡傅立葉變換成像。這為非可見光量子成像，如 X 射線繞射成像，提供了理論上的可行性。2006 年，史硯華小組實現了無透鏡或其他等效成像系統的贗熱光源量子成像，這種方案適合任何波長的輻射作為光源，並且形成圖像的過程中不需要任何成像透鏡，因此這種方案對於 X 射線、伽瑪射線和其他波長光源的成像應用幫助很大。

### 計算量子成像

夏皮羅也提出了計算量子成像的理論。與傳統量子成像不同，計算量子成像方案僅保留了包含待成像物體的測量光路和桶探測器，通過雷射照射空間光調製器產生可調強度、相位等參數的空間調製光場（又稱為主動式光源），再根據繞射理論計算得到原參考光路在無透鏡量子成像中可以得到的特定位置的光強分布，與測量光場進行符合關聯得到圖像。這種方案所用的裝置可以生成沒有背景雜訊的圖像，其解析度和成像區域可以通過調整空間光調製器的參數來控制。2009 年，以色列科學家 B. Sun 等人在 3D 成像實驗中驗證了計算量子成像的可行性。

### 差分量子成像

2010 年，義大利學者路易吉・盧賈托等提出差分量子成像（差分鬼成像，英語：Differential Ghost Imaging，DGI）的方案，該方案提高了量子成像信噪比的數量級，大幅提高了量子成像的成像質量，在一些強幹擾和背景雜訊較大的環境下表現較為良好，但需要大量的測量數據和更複雜的計算。2012 年，羅開紅及其同事提出一種名為「對應成像」（英語：Correspondence Imaging，CI）的技術，可概述為選擇桶探測器在正向或負向的強度波動（正信號與負信號），對對應的參考光路所得的數據進行條件平均，而不是用桶探測器獲得的光強與參考光路所得數據直接相乘。2013 年，中科院物理研究所吳令安小組提出時間對應差分量子成像，將差分量子成像和對應成像的優點結合起來，降低了數據處理難度，縮短了成像所需的計算時間。

## 原理

以下嘗試分析史硯華小組 1995 年所做的纏結光源量子成像實驗，簡述量子成像的原理。

自發參量下轉換過程之後，泵浦光子分裂成為一對信號光子 - 閒置光子，其能量滿足

其中 $\omega _{s}$ 是信號光子的能量， $\omega _{i}$ 是閒置光子的能量， $\omega _{p}$ 是泵浦光子的能量，「 $\cong $ 」代表近似等於。雙光子的的橫向位矢和動量 EPR $\delta $ 函數表達式分別為

和其中， ${\overrightarrow {\rho _{s}}}$ 和 ${\overrightarrow {\rho _{i}}}$ 分別是信號光子和閒置光子的位矢， ${\overrightarrow {k_{s}}}$ 和 ${\overrightarrow {k_{i}}}$ 分別是信號光子和閒置光子的動量。雖然光子可能從非線性晶體表面的任何一點出射，但 $\delta ({\overrightarrow {\rho _{s}}}-{\overrightarrow {\rho _{i}}})$ 表明，如果在某一空間位置找到信號 - 閒置光子對中的一個光子，則必在對應位置找到另一個光子，即光子對必從同一點出射。 $\delta ({\overrightarrow {k_{s}}}+{\overrightarrow {k_{i}}})$ 定義了光子對的角度相關性： ${\overrightarrow {k_{s}}}=-{\overrightarrow {k_{i}}}$ ，即雙光子必在相對出射點大致相等但相反的角度存在。

定義 $G^{(2)}({\overrightarrow {\rho _{o}}},{\overrightarrow {\rho _{i}}})$ ； ${\overrightarrow {\rho _{o}}}$ 和 ${\overrightarrow {\rho _{i}}}$ 是物平面和像平面上點的橫向位矢。以下將證明，物平面和像平面之間存在一個 $\delta $ 函數，即存在一種「點對點」的對應關係，若在物平面上 ${\overrightarrow {\rho _{o}}}$ 位置找到信號光子，則必在滿足 $\delta ({\overrightarrow {\rho _{s}}}-{\overrightarrow {\rho _{i}}}/m)$ 關係的位置找到閒置光子。 $m=-{\frac {s_{i}}{s_{o}}}$ ，是像對物體的放大倍數，其中 $s_{o}$ 是物平面與透鏡的光學距離， $s_{i}$ 是像平面與透鏡的光學距離。

確定 $\delta $ 函數對應的相關性後，以下展示物函數 $A({\overrightarrow {\rho _{o}}})$ 如何轉到 $A({\overrightarrow {\rho _{i}}}/m)$ 。上圖中的一條「光線」實際上表示一對信號 - 閒置光子。 $(r_{1},t_{1})$ 和 $(r_{2},t_{2})$ 在物平面和像平面上的對應關係實際上是雙光子振幅疊加的結果。

建立格林函數 $g({\overrightarrow {k_{s}}},\omega _{s},{\overrightarrow {\rho _{o}}},z_{o})$ 和 $g({\overrightarrow {k_{i}}},\omega _{i},{\overrightarrow {\rho _{2}}},z_{2})$ 。信號光路和閒置光路的測量方法非常不同，需要分開考慮。

### 信號光路

在信號光路，信號光在從源輸出到成像透鏡的 $d_{1}$ 的距離上自由傳播，通過 $s_{o}$ 距離的物光路，到達收集透鏡，再聚焦到 $D_{1}$ 上，光學傳遞函數為：

 ${\overrightarrow {\rho _{s}}}$ 和 ${\overrightarrow {\rho _{l}}}$ 分別是光源輸出平面和成像透鏡平面上定義的橫向位矢。上式中，前後兩個大括號分別代表光子從輸出平面到成像透鏡平面，和成像透鏡平面到物平面的自由空間傳播過程。 $e^{i{\frac {\omega _{s}}{2cd_{1}}}\left\vert {\overrightarrow {\rho _{s}}}-{\overrightarrow {\rho _{l}}}\right\vert ^{2}}$ 和 $e^{i{\frac {\omega _{s}}{2cs_{o}}}\left\vert {\overrightarrow {\rho _{l}}}-{\overrightarrow {\rho _{o}}}\right\vert ^{2}}$ 是菲涅爾相位。把成像透鏡視為薄透鏡，則其變換函數近似為一個高斯函數 $l(\left\vert {\overrightarrow {\rho _{l}}}\right\vert ,f)\cong e^{-i{\frac {\omega }{2cf}}\left\vert {\overrightarrow {\rho _{l}}}\right\vert ^{2}}$ 。

### 閒置光路

在閒置光路，閒置光從源輸出到成像透鏡，整個 $d_{2}$ 距離上都是自由傳播的，光學傳遞函數為：

 ${\overrightarrow {\rho '_{s}}}$ 和 ${\overrightarrow {\rho _{2}}}$ 分別是源平面和 $\delta $ 光電探測器平面上的橫向位矢。

### 全光路（物平面 - 像平面）

為簡化計算，假設 $\omega _{s}=\omega _{i}=\omega $ ，SPDC 雙光子共線，忽略所有比例常數（因此將 $=$ 替換為 $\propto $ ），將 $g({\overrightarrow {k_{s}}},\omega _{s},{\overrightarrow {\rho _{o}}},z_{o})$ 和 $g({\overrightarrow {k_{i}}},\omega _{i},{\overrightarrow {\rho _{2}}},z_{2})$ 代入像平面有效波函數（英語：effective wave function）

可求出雙光子有效波函數表達式

完成二重積分

後，雙光子有效波函數表達式變形為

完成對 $d{\overrightarrow {\rho _{s}}}$ 的積分後，可變形為

用 $s_{i}$ 代替 $d_{1}+d_{2}$ ， $d{\overrightarrow {\rho _{l}}}$ 上的積分在物平面和像平面之間產生的點對點對應關係，由高斯薄透鏡方程式定義：

用 ${\overrightarrow {\rho _{2}}}$ 等效代替 ${\overrightarrow {\rho _{i}}}$ ，則上式中的積分值趨向無窮大，應用高斯薄透鏡方程式

將 $m={\frac {s_{i}}{s_{o}}}$ 定義為成像系統放大倍數， $\delta ({\overrightarrow {\rho _{o}}}+{\frac {\overrightarrow {\rho _{i}}}{m}})$ 表示物平面上的每一個 ${\overrightarrow {\rho _{o}}}$ 與像平面上的每一個 ${\overrightarrow {\rho _{i}}}$ 唯一對應，兩向量方向相反，大小比為 $m={\frac {\left\vert {\overrightarrow {\rho _{i}}}\right\vert }{\left\vert {\overrightarrow {\rho _{o}}}\right\vert }}$ 。
考慮到成像透鏡尺寸有限（設半徑為 R），則積分會產生點擴散函數 $somb(x)$ ：

其中 $somb(x)=2J_{1}(x)$ ， $J_{1}(x)$ 是一階貝索函數，將物 - 像平面的點對點關係轉化為點對像素關係，限制了成像的空間解析度。在代入高斯薄透鏡方程式後，橫向雙光子有效波函數近似為一個 $\delta $ 函數：

這表現了物 - 像平面之間的點對點相關性，即若在物平面 ${\overrightarrow {\rho _{o}}}$ 位置找到信號光子，則必然在像平面 ${\overrightarrow {\rho _{i}}}$ 位置找到閒置光子， ${\overrightarrow {\rho _{o}}}$ 和 ${\overrightarrow {\rho _{i}}}$ 滿足 ${\overrightarrow {\rho _{o}}}+{\frac {\overrightarrow {\rho _{i}}}{m}}=0$ ， $m={\frac {s_{i}}{s_{o}}}$ 。
在信號光路的光學傳遞函數中再加入作為物體的光圈、收集透鏡和光子計數檢測器，其中後兩者整體可視為一個桶探測器。桶探測器對通過物體（光圈） $A({\overrightarrow {\rho _{o}}})$ ，觸發光電聯合檢測的光子進行積分：

同時，閒置光路中的探測器 $D_{2}$ 再次掃描像平面， ${\displaystyle {\overrightarrow {\rho _{2}}}}={\displaystyle {\overrightarrow {\rho _{i}}}}$ 。
綜上所述，位置 EPR 相關性是雙光子振幅相干疊加的結果：

原則上，一對信號 - 閒置光子對包含用於在物平面和像平面間產生點對點對應關係的雙光子振幅，因此量子成像又稱為「雙光子關聯成像」。

## 相關領域與應用

### 光刻技術

2000 年，博托（英語：A. N. Boto）小組提出了 N 光子吸收光刻技術，這種技術使用纏結光子流取代古典光的非相干性光子流，降低光刻技術中的最小可分辨特徵尺寸。高度纏結的光子可以使得光刻技術的最小可分辨特徵尺寸突破瑞立繞射極限規定的最小值 $\lambda /2$ ，達到 $\lambda /2N$ （N 是一個泵浦光子分裂成為的一組纏結光子的數量，這些光子最後都會被光刻膠吸收）。其原理可以簡述為：在古典光中，N 個光子到達某一特定空間區域的機率是單個光子到達該範圍的 N 次方，但纏結光中只要確定其中一個光子到達的區域，其他 N-1 個光子會到達的區域是確定的，如果光學系統對的足夠准，則 N 個纏結光子到達某一特定區域的機率就只需要計算一次，這使得 N 光子光刻不需要光焦度達到不切實際的程度，僅使用與古典器件相同的功率級別，即可使光刻最小可分辨特徵尺寸降低 N 倍。

### 雷射雷達技術

傳統雷射雷達分為兩種類型：掃描成像雷射雷達和非掃描成像雷射雷達。掃描成像雷射雷達通過用脈衝雷射逐點掃描目標區域來獲得目標的真實空間圖像，這種雷達難以對高速運動物體進行成像；非掃描成像雷射雷達用脈衝閃光雷射源和高解析度成像系統進行成像，一次曝光即可獲得目標的真實空間圖像，但目標反射的光強是由 CCD 相機許多小像素接收到的，因此檢測靈敏度較低，其檢測距離受到成像系統光路和整個成像平面的信噪比的影響。相較之下，量子成像雷射雷達具有遙感距離長、成像速度快和成像解析度高等優點。

2009 年，以色列科學家亞倫・希爾伯格（英語：Yaron Silberberg）等進行了計算量子成像的驗證實驗，提出了用於贗熱量子成像的圖像重建高級算法，並提出計算量子成像可以用於雷射雷達。2011 年，美國麻省理工學院學者提出計算量子成像用於遙感成像的方案，並分析了這種方案的性能。2012 年，上海光機所韓申生團隊實現了基於稀疏約束的關聯成像雷達，並在約 1km 的距離上進行了高解析度成像。2013 年，韓申生團隊實現了 3D 量子成像雷射雷達。2015 年，韓申生小組提出了結構化圖像重構方法，這種方法可以更準確地恢復具有各種稀疏比率的場景切片，並且具有明顯的增強解析度的效果，這適用於測量次數較少的量子成像過程，大幅提高了三維量子成像雷射雷達的成像質量。2016 年，清華大學戴瓊海小組提出了內容自適應計算量子成像方法，用以在測量次數較少的情況下重建較高質量的圖像，從而完成對動態目標的成像。

## 參見

* 量子力學
* 量子力學入門
* 量子光學

## 參考文獻