Signals of entanglement and nonlocality are quantitatively evaluated at zero and finite temperature in an analogue black hole realized in the flow of a quasi one-dimensional Bose-Einstein condensate. The violation of Lorentz invariance inherent to this analog system opens the prospect to observe 3-mode quantum correlations and we study the corresponding violation of bipartite and tripartite Bell inequalities. It is shown that the long wavelength modes of the system are maximally entangled, in the sense that they realize a superposition of continuous variable versions of Greenberger-Horne-Zeilinger states whose entanglement resists partial tracing.
论文ID : 2404.16497标题 : Violation of Bell inequalities in an analogue black hole作者 : Giorgio Ciliberto, Stephanie Emig, Nicolas Pavloff, Mathieu Isoard分类 : quant-ph, cond-mat.quant-gas发表时间 : 2024年4月 (arXiv预印本,2025年10月更新)论文链接 : https://arxiv.org/abs/2404.16497 本文在准一维玻色-爱因斯坦凝聚体(BEC)流动中实现的模拟黑洞系统中,定量评估了零温和有限温度下的纠缠和非局域性信号。该模拟系统固有的洛伦兹不变性破缺为观察三模量子关联提供了可能性,研究了相应的双粒子和三粒子贝尔不等式违背。结果表明,系统的长波长模式是最大纠缠的,它们实现了连续变量版本的Greenberger-Horne-Zeilinger(GHZ)态的叠加,其纠缠性质能够抵抗部分迹操作。
模拟引力物理 :模拟引力领域旨在为无法在传统引力环境中直接观测的广义相对论现象提供实验室模型,特别是黑洞超辐射和霍金辐射。量子纠缠与非局域性 :在模拟黑洞系统中,霍金辐射的量子特性表现为相关准粒子对的产生,这为研究量子纠缠和非局域性提供了独特平台。技术挑战 :现有研究主要关注双粒子纠缠,而三粒子或多粒子量子关联的研究相对较少,特别是在连续变量系统中的贝尔不等式违背。理论完善 :系统研究模拟黑洞中的量子纠缠性质,特别是三模式系统的非局域性。实验指导 :为未来的BEC模拟黑洞实验提供理论预测和量化评估。基础物理 :探索洛伦兹不变性破缺如何影响量子纠缠的性质。系统性理论分析 :首次系统性地研究了BEC模拟黑洞中双粒子和三粒子贝尔不等式的违背。GHZ态结构发现 :证明了长波长模式具有连续变量GHZ态的结构,且在部分迹操作后仍保持纠缠。温度效应量化 :详细分析了有限温度对量子纠缠和非局域性的影响。多种黑洞配置比较 :研究了瀑布型、δ峰型和平坦型三种不同的模拟黑洞配置。光学等效模型 :建立了参数下转换和分束器组成的光学等效系统。研究准一维BEC流动中实现的模拟黑洞系统的量子纠缠和非局域性,具体包括:
输入:BEC背景流的参数(马赫数、温度等) 输出:CHSH参数、Svetlichny参数、Mermin参数等非局域性量度 约束:一维平均场近似的有效性条件 系统由量子场算符 Ψ ^ ( x , t ) \hat{\Psi}(x,t) Ψ ^ ( x , t ) 描述,满足:
i ℏ ∂ t Ψ ^ = − ℏ 2 2 m ∂ x 2 Ψ ^ + [ U ( x ) + g n ^ − μ ] Ψ ^ i\hbar\partial_t\hat{\Psi} = -\frac{\hbar^2}{2m}\partial_x^2\hat{\Psi} + [U(x) + g\hat{n} - \mu]\hat{\Psi} i ℏ ∂ t Ψ ^ = − 2 m ℏ 2 ∂ x 2 Ψ ^ + [ U ( x ) + g n ^ − μ ] Ψ ^
采用Bogoliubov近似:Ψ ^ ( x , t ) = Φ ( x ) + ψ ^ ( x , t ) \hat{\Psi}(x,t) = \Phi(x) + \hat{\psi}(x,t) Ψ ^ ( x , t ) = Φ ( x ) + ψ ^ ( x , t )
出射模式与入射模式通过散射矩阵S(ω)关联:
( c ^ 0 c ^ 1 c ^ 2 † ) = ( S 00 S 01 S 02 S 10 S 11 S 12 S 20 S 21 S 22 ) ( b ^ 0 b ^ 1 b ^ 2 † ) \begin{pmatrix} \hat{c}_0 \\ \hat{c}_1 \\ \hat{c}_2^\dagger \end{pmatrix} = \begin{pmatrix} S_{00} & S_{01} & S_{02} \\ S_{10} & S_{11} & S_{12} \\ S_{20} & S_{21} & S_{22} \end{pmatrix} \begin{pmatrix} \hat{b}_0 \\ \hat{b}_1 \\ \hat{b}_2^\dagger \end{pmatrix} c ^ 0 c ^ 1 c ^ 2 † = S 00 S 10 S 20 S 01 S 11 S 21 S 02 S 12 S 22 b ^ 0 b ^ 1 b ^ 2 †
其中模式2具有负范数,需要反向量子化。
采用GKMR伪自旋算符进行离散化:
Π ^ x ( j ) = ∫ 0 + ∞ d q ( ∣ q ⟩ j j ⟨ q ∣ − ∣ − q ⟩ j j ⟨ − q ∣ ) \hat{\Pi}^{(j)}_x = \int_0^{+\infty} dq (|q\rangle_j {}_j\langle q| - |-q\rangle_j {}_j\langle -q|) Π ^ x ( j ) = ∫ 0 + ∞ d q ( ∣ q ⟩ j j ⟨ q ∣ − ∣ − q ⟩ j j ⟨ − q ∣ ) Π ^ y ( j ) = i ∫ 0 + ∞ d q ( ∣ q ⟩ j j ⟨ − q ∣ − ∣ − q ⟩ j j ⟨ q ∣ ) \hat{\Pi}^{(j)}_y = i\int_0^{+\infty} dq (|q\rangle_j {}_j\langle -q| - |-q\rangle_j {}_j\langle q|) Π ^ y ( j ) = i ∫ 0 + ∞ d q ( ∣ q ⟩ j j ⟨ − q ∣ − ∣ − q ⟩ j j ⟨ q ∣ ) Π ^ z ( j ) = ∫ − ∞ + ∞ d q ∣ q ⟩ j j ⟨ − q ∣ \hat{\Pi}^{(j)}_z = \int_{-\infty}^{+\infty} dq |q\rangle_j {}_j\langle -q| Π ^ z ( j ) = ∫ − ∞ + ∞ d q ∣ q ⟩ j j ⟨ − q ∣
连续变量贝尔不等式 :将离散变量的贝尔不等式扩展到连续变量系统,通过伪自旋算符实现离散化。三模式分析 :突破传统双模式分析限制,系统研究三模式量子关联。Wigner函数方法 :利用Wigner变换技术计算高斯态上伪自旋算符的期望值。遗传算法优化 :采用遗传算法数值优化多参数贝尔算符的期望值。瀑布配置 :下游马赫数 m d = 2.9 m_d = 2.9 m d = 2.9 (对应Technion 2019实验)温度范围 :T = 0 T = 0 T = 0 到 T = 0.2 g n u T = 0.2 g n_u T = 0.2 g n u 频率范围 :ω ∈ [ 0 , Ω ] \omega \in [0, \Omega] ω ∈ [ 0 , Ω ] ,其中Ω \Omega Ω 为阈值频率PPT纠缠度量 :Λ ( i ∣ j ) ∈ ] − ∞ , 1 ] \Lambda^{(i|j)} \in ]-\infty, 1] Λ ( i ∣ j ) ∈ ] − ∞ , 1 ] CHSH参数 :B ( i ∣ j ) = max a , a ′ , b , b ′ ⟨ B ^ ( i ∣ j ) ⟩ B^{(i|j)} = \max_{a,a',b,b'} \langle\hat{B}^{(i|j)}\rangle B ( i ∣ j ) = max a , a ′ , b , b ′ ⟨ B ^ ( i ∣ j ) ⟩ Svetlichny参数 :S ( 0 ∣ 1 ∣ 2 ) = max ⟨ S ^ ( 0 ∣ 1 ∣ 2 ) ⟩ S^{(0|1|2)} = \max \langle\hat{S}^{(0|1|2)}\rangle S ( 0∣1∣2 ) = max ⟨ S ^ ( 0∣1∣2 ) ⟩ Mermin参数 :M ( 0 ∣ 1 ∣ 2 ) = max ∣ ⟨ M ^ ( 0 ∣ 1 ∣ 2 ) ⟩ ∣ M^{(0|1|2)} = \max |\langle\hat{M}^{(0|1|2)}\rangle| M ( 0∣1∣2 ) = max ∣ ⟨ M ^ ( 0∣1∣2 ) ⟩ ∣ 瀑布配置 (waterfall) δ峰配置 (delta peak) 平坦配置 (flat profile) 零温下,CHSH参数在瀑布配置(m d = 2.9 m_d = 2.9 m d = 2.9 )中达到最大值约2.25 模式(0|2)和(1|2)都表现出纠缠,但模式(0|1)始终可分离 贝尔不等式违背对温度比纠缠更敏感 零温长波长极限下,Svetlichny参数达到上界2 2 2\sqrt{2} 2 2 Mermin参数在ω = 0 , T = 0 \omega = 0, T = 0 ω = 0 , T = 0 时达到最大值4 系统表现出GHZ悖论的特征 有限温度显著抑制非局域性信号 T = 0.1 g n u T = 0.1 g n_u T = 0.1 g n u 时,三粒子非局域性基本消失纠缠度量对温度的抗性优于贝尔参数 通过分析不同马赫数配置发现:
m u > 0.6 m_u > 0.6 m u > 0.6 :纠缠主要集中在模式(0|2)m u < 0.2 m_u < 0.2 m u < 0.2 :纠缠主要集中在模式(1|2)极限情况m u → 0 , 1 m_u \to 0,1 m u → 0 , 1 时系统退化为EPR态 零温长波长极限下,入射真空态可表示为:
∣ 0 ω = 0 ⟩ i n = ∑ l , m , n C l , m , n ( ∣ x l + , x m − , x n − ⟩ + ∣ x l − , x m + , x n + ⟩ ) |0_{\omega=0}\rangle_{in} = \sum_{l,m,n} C_{l,m,n} (|x_l^+, x_m^-, x_n^-\rangle + |x_l^-, x_m^+, x_n^+\rangle) ∣ 0 ω = 0 ⟩ in = ∑ l , m , n C l , m , n ( ∣ x l + , x m − , x n − ⟩ + ∣ x l − , x m + , x n + ⟩)
这是无穷多个简并GHZ态的叠加,解释了为什么部分迹后仍保持纠缠。
表面重力波 :Weinfurtner等人的水槽实验光学系统 :Philbin等人的光纤模拟事件视界BEC实验 :Steinhauer等人在Technion的开创性工作伪自旋方法 :Banaszek-Wodkiewicz和GKMR伪自旋贝尔不等式 :从CHSH到Svetlichny不等式的扩展纠缠度量 :PPT判据、对数负性、纠缠熵等BEC模拟黑洞能够产生显著的双粒子和三粒子量子纠缠 长波长模式具有连续变量GHZ态的独特结构 洛伦兹不变性破缺使得三模式纠缠成为可能 温度是观测量子非局域性的主要限制因素 温度敏感性 :三粒子非局域性对温度极其敏感,实验观测困难理论近似 :基于高斯态假设,未考虑非线性回馈效应测量挑战 :某些量子关联函数的实验测量技术复杂非线性效应 :研究非高斯修正对量子关联的影响其他平台 :扩展到激子-极化子等其他模拟系统实验实现 :开发更精确的量子关联测量技术理论完整性 :系统性地研究了模拟黑洞的量子信息性质数学严谨性 :采用Wigner函数等成熟的量子光学工具物理洞察 :发现了连续变量GHZ态的新结构实验相关性 :参数选择贴近实际BEC实验条件实验可行性 :三粒子非局域性的观测在技术上极具挑战性近似限制 :一维平均场近似可能低估某些量子效应普适性 :结果主要针对特定的BEC系统,推广性有限学术价值 :为模拟引力和量子信息的交叉研究提供了重要理论基础实用价值 :为未来BEC量子模拟实验提供了明确的理论预测可复现性 :数学推导详细,数值方法清晰,便于验证和扩展基础研究 :量子引力、模拟引力物理量子技术 :连续变量量子信息处理实验物理 :超冷原子、BEC量子模拟论文引用了116篇相关文献,涵盖了模拟引力、量子信息、BEC物理等多个领域的重要工作,特别是:
Hawking (1975) 的黑洞辐射理论 Steinhauer等人的BEC模拟黑洞实验 Bell不等式和量子纠缠的理论基础 连续变量量子信息的数学工具