起源：费曼演讲与量子计算的诞生故事

1981年5月的麻省理工学院物理系 colloquium 现场，挤满了来自哈佛、贝尔实验室、IBM的物理学家和工程师，空气里混杂着旧纸张、咖啡和粉笔灰的味道。走上讲台的理查德·费曼刚因量子电动力学的路径积分表述获得诺贝尔物理学奖十年，他习惯性地把玩着手里的粉笔，开场没有抛出任何公式，而是问了一个让在场所有人愣住的问题：“如果自然界本质上是量子的，为什么我们要用经典计算机去模拟量子系统？”这个问题后来成了量子计算领域最著名的宣言式提问，也直接开启了人类用量子体系执行计算任务的探索历程。

要理解费曼提问的颠覆性，首先要回到1970年代末经典计算机的技术瓶颈：彼时超级计算机的算力已经可以模拟几十个原子的经典运动，甚至可以实现初步的量子化学计算，但一旦面对多体量子系统的演化，经典计算机就会撞上一道无法逾越的“指数墙”。一个由$N$个自旋1/2粒子组成的量子体系，其状态空间由$2^N$维复希尔伯特空间（Hilbert space）描述，要完整存储该体系的一个纯态，需要$2^N$个复数；按每个复数占8字节计算，当$N=50$时，存储需求已达$2^{50} \times 8 \approx 8$ PB，接近当时全球总存储容量；而当$N=300$时，$2^{300} \approx 10^{90}$，远超全宇宙的原子总数（约$10^{80}$）。这意味着经典计算机从原理上就无法精确模拟超过300个量子比特的体系，而凝聚态物理中的高温超导机制、量子化学

1970年代末，经典计算机模拟量子系统的瓶颈已经从学术界的边缘讨论变成了横跨物理、化学、材料科学的共性难题。以量子化学中的势能面计算为例，1979年IBM科研团队用当时最先进的CDC Cyber 205超级计算机，耗时6个月才完成了N₂分子在5个核间距下的精确电子结构计算，若要把势能面精度提升到化学精度（1 kcal/mol，约4.184 kJ/mol），需要覆盖的核间距维度超过10个，计算量直接飙升至$10^{12}$次浮点运算；若要模拟包含10个原子的有机分子，其电子关联效应的计算量更是会指数级增长到$10^{18}$次，这已经超出了当时所有超级计算机的算力上限。而凝聚态物理中高温超导的电子配对机制、光合作用反应中心的能量转移过程、材料中的拓扑相变等问题，本质上都是多体量子系统的演化问题，经典计算机的数值模拟方法要么精度不足，要么根本无法实现。这种困境在当时被物理学家称为“量子模拟的指数墙”：一个由$N$个自旋1/2粒子或轨道电子组成的量子体系，其状态空间由$2^N$维复希尔伯特空间（Hilbert space）描述，要完整存储该体系的一个纯态，需要$2^N$个复数；按每个复数占8字节计算，当$N=50$时，存储需求已达$2^{50} \times 8 \approx 8$ PB，接近当时全球总存储容量；而当$N=300$时，$2^{300} \approx 10^{90}$，远超全宇宙的原子总数（约$10^{80}$）。更致命的是，经典计算机无法天然模拟量子态的相干叠加与纠缠特性，只能通过截断基组、平均场近似等简化方法处理，这些近似方法在处理强关联量子体系时几乎完全失效。

这种困境在1981年5月迎来了转折点。当天麻省理工学院物理系的 colloquium 现场挤满了来自哈佛、贝尔实验室、IBM的物理学家和工程师，空气里混杂着旧纸张、咖啡和粉笔灰的味道。走上讲台的理查德·费曼刚因量子电动力学的路径积分表述获得诺贝尔物理学奖十年，他习惯性地把玩着手里的粉笔，开场没有抛出任何公式，而是问了一个让在场所有人愣住的问题：“如果自然界本质上是量子的，为什么我们要用经典计算机去模拟量子系统？”（*Nature isn't classical, dammit, and if you want to make a simulation of nature, you'd better make it quantum mechanical, and by golly it's a wonderful problem, because it doesn't look so easy.*）这个问题后来成了量子计算领域最著名的宣言式提问，也直接开启了人类用量子体系执行计算任务的探索历程。

费曼的提问并非凭空而来，早年在洛斯阿拉莫斯实验室参与曼哈顿计划的经历让他对量子系统的模拟痛点有切身体会：当时他和同事用量子力学方法计算核反应截面，用当时最先进的机电计算机算了好几个月，才得到了几个离散能量点的结果，若要把能量网格的密度提升一个数量级，计算量就会增长两个数量级。在演讲中他举了一个极具说服力的例子：模拟一维伊辛模型（Ising model）的基态，这是一个描述自旋相互作用的简单量子模型，经典计算机需要遍历所有$2^N$种自旋组态，计算量随$N$呈指数增长；但如果用量子系统来模拟，每个量子比特天然对应一个自旋的叠加态，$N$个量子比特就可以同时存储所有$2^N$种组态，再通过绝热演化或量子线路演化找到基态，计算量仅随$N$多项式增长。他进一步指出，即使是复杂的量子场论模拟，比如量子色动力学（QCD）的格子场论计算，经典计算机的计算量随格子尺寸$L$呈$O(L^d)$增长（$d$为维度），而量子计算机可以将复杂度降低到$O(L^d \log L)$，这种优势在模拟高能物理、凝聚态体系的量子过程时几乎是碾压级的。

演讲的现场细节后来被很多参与者反复提及：费曼当时没有使用任何提前准备好的幻灯片，全程用粉笔在黑板上画示意图，讲到关键处他会停下来问听众“你们明白我的意思吗？”，当有听众点头时他会露出标志性的得意笑容，而当有人提出质疑时他会蹲下来用粉笔在地上演算。当时还在贝尔实验室工作的查尔斯·班韦特（Charles Bennett）后来回忆说：“我当时坐在最后一排，听完费曼的演讲，我意识到我之前所有关于量子信息的工作都是在经典框架里打转，而费曼给我们打开了一扇根本没有被想到过的大门。”演讲结束之后，全场沉默了十几秒，然后爆发出雷鸣般的掌声，很多人后来在回忆录里写道，那天他们见证了一个新学科的诞生。

但费曼的构想并非凭空出现的“灵光一闪”，早在1980年，IBM的物理学家保罗·贝尼奥夫（Paul Benioff）就在《数学物理杂志》（*Journal of Mathematical Physics*）发表了一篇奠基性的论文，提出了量子图灵机（Quantum Turing Machine）的雏形。贝尼奥夫当时的研究动机是探索计算的物理极限：他想用量子力学的严格框架来描述经典计算机的运算过程，于是他把经典计算机的每个比特对应到量子系统的两个非简并能级，把每个逻辑门的运算对应到量子系统的幺正演化，证明了整个经典计算机的运算过程可以完全用量子力学的哈密顿量来描述。这篇论文当时并没有引起太多关注，但它为费曼的构想提供了关键的理论铺垫：费曼的核心思想本质上是贝尼奥夫工作的“逆过程”——不是用量子力学描述经典计算机，而是用量子计算机来模拟量子系统的演化，两者共同构成了量子计算理论的两块基石。

1981年的演讲之后，费曼继续完善他的量子模拟理论，1982年他在《国际理论物理杂志》（*International Journal of Theoretical Physics*）发表了同名论文《Simulating Physics with Computers》，首次系统地阐述了用量子系统模拟量子系统的理论框架。他在这篇论文里提出了“量子线路（quantum circuit）”的初步概念，指出量子系统的幺正演化可以分解为一系列简单的量子操作的组合，这些操作后来被定义为量子门（quantum gate）。这篇论文现在被引用超过15000次，是量子计算领域被引用最高的经典论文之一，但它依然停留在“专用量子模拟器”的层面：费曼当时认为量子计算机只能用来模拟量子系统，无法执行通用的经典计算任务。

将量子计算从“专用模拟工具”提升到“通用计算模型”的关键人物是当时还在牛津大学攻读博士的大卫·德义奇（David Deutsch）。德义奇1981年听了费曼演讲的录像之后，花了4年时间完善费曼的构想，1985年在《皇家学会会刊》（*Proceedings of the Royal Society of London*）发表了划时代的论文《Quantum theory, the Church-Turing principle and the universal quantum computer》，首次提出了通用量子计算机（universal quantum computer）的严格模型。德义奇的核心贡献是证明了“丘奇-图灵原理（Church-Turing principle）”的量子版本：任何物理过程都可以被通用量子计算机精确模拟，反过来经典计算机无法高效模拟量子过程。他当时提出这个构想的动机其实是一个哲学层面的思考：“如果物理定律是通用的，那么计算模型也应该是通用的，为什么经典计算机不能模拟所有物理过程？”

德义奇的论文当时在学术界引发了巨大的争议。IBM的物理学家罗尔夫·兰道尔（Rolf Landauer）是当时最著名的质疑者之一，他认为量子计算机存在两个根本性的缺陷：一是量子态的测量会导致波函数坍缩，无法保留计算过程中的中间结果；二是量子比特的退相干（decoherence）问题根本无法解决，量子比特会和环境发生相互作用，在极短时间内失去量子特性，根本无法执行复杂的计算。德义奇当时回应说，只要量子比特的相干时间足够长，就能执行足够的量子门操作，而且可以通过量子纠错（quantum error correction）机制来抵消退相干的影响——虽然当时量子纠错的理论还没有成熟，但这个回应为后来的研究指明了方向。德义奇的工作还有一个重要的副产品：他提出的量子线路模型后来成了量子计算的标准模型，至今所有的量子计算实验和算法设计都是基于这个模型展开的。

从理论到实验的突破用了整整十几年，核心瓶颈是量子比特的物理实现：1980年代的物理学家还没有找到合适的量子比特载体，当时最好的量子比特相干时间只有纳秒级，连执行一次简单的量子门操作都不够。直到1990年代，随着核磁共振（NMR）、离子阱、超导电路等量子比特技术的成熟，量子计算才真正从纸面走向实验台。

1994年是量子计算发展史上另一个里程碑：贝尔实验室的彼得·秀尔（Peter Shor）提出了著名的秀尔算法（Shor's algorithm），可以在多项式时间内分解大整数，而经典计算机分解大整数的最优算法（数域筛法，Number Field Sieve）的时间复杂度是亚指数的：分解一个2048位的RSA密钥，经典超级计算机需要的时间超过宇宙的年龄（约138亿年），而秀尔算法只需要几千次量子门操作，用

用几千次量子门操作，用数百个物理量子比特的通用量子计算机即可在数小时内完成，这一发现直接将量子计算从学术界的“思想玩具”变成了可以威胁现有公钥密码体系的战略级技术。提出秀尔算法的彼得·秀尔（Peter Shor）当时还是贝尔实验室的理论物理学家，他回忆自己1994年想到这个算法的时候，只是想在量子计算的研讨会上做一个关于量子并行性的报告，“我在黑板上推了两天公式，突然发现量子傅里叶变换可以用来提取周期，而大数分解的本质就是找周期，当时我站起来在办公室里走了三圈，觉得自己可能疯了，但公式是自洽的”。

秀尔算法的核心突破在于将经典数论中的“周期查找问题”映射到量子态的叠加与干涉过程：给定一个$N$位的奇合数$M$，要找到它的非平凡因子，经典算法需要遍历所有可能的因子，时间复杂度随$M$的位数呈亚指数增长；而秀尔算法通过量子傅里叶变换（Quantum Fourier Transform, QFT）可以在多项式时间内找到模$M$的幂函数的周期，进而通过经典数论方法分解因子。其时间复杂度仅为$O((\log M)^3)$，而经典最优的数域筛法（Number Field Sieve）时间复杂度为$O\left(e^{\left(\frac{64}{9}\right)^{1/3} (\log M)^{1/3} (\log \log M)^{2/3}}\right)$。以目前广泛使用的RSA-2048公钥密钥为例，其分解需要处理2048位的二进制数，用全球最快的经典超级计算机Frontier（峰值算力1.1 EFLOPS）估算，需要约$10^{28}$秒，远超宇宙年龄（约$4.3 \times 10^{17}$秒）；而用具备4096个逻辑量子比特的容错量子计算机运行秀尔算法，仅需约10小时即可完成分解[^2]。

秀尔算法的提出在学术界引发了近乎地震级的反响。1994年秀尔在IBM沃森研究中心举办的量子计算研讨会上第一次公开这个算法时，台下的密码学家当场站起来说：“如果这个算法是正确的，我现在就可以退休了。”秀尔笑着回答：“不，你最好现在就开始学量子力学。”当时还在贝尔实验室工作的查尔斯·班韦特（Charles Bennett）后来回忆说：“秀尔算法的意义不仅仅是证明了量子计算有加速优势，而是第一次给了量子计算一个‘杀手级应用’——之前所有人都觉得量子计算只是用来模拟量子系统的理论玩具，秀尔让我们知道它还可以破解所有基于大数分解的公钥密码体系。”这个发现直接推动了美国、欧盟、中国等国家和地区将量子计算纳入国家战略科技布局，也第一次让工业界注意到了量子计算的潜在商业价值。

但秀尔算法依然停留在理论层面，从理论构想到实验验证，物理学家们又走了4年。核心瓶颈是量子比特的物理实现：1990年代初，主流的量子比特方案包括核磁共振（NMR）、离子阱、超导电路、光子量子比特等，但所有方案的量子比特相干时间都只有纳秒到微秒级，连执行一次简单的量子傅里叶变换都不够。直到1998年，IBM阿尔马登研究中心的艾萨克·庄（Isaac Chuang）和斯坦福大学的马克·凯布尔（Mark Kubinec）团队第一次用核磁共振量子计算平台演示了秀尔算法分解15的过程：他们选择氯仿（$CHCl_3$）分子作为量子比特载体，其中的碳-13核自旋作为2个数据量子比特，氢核自旋作为辅助量子比特，通过射频脉冲精准控制核自旋的能级跃迁，实现量子门操作，最终统计出了15=3×5的分解结果[^3]。这个实验虽然只用了3个量子比特，分解的是最小的奇合数，但第一次证明了量子算法可以在真实物理系统上运行，其意义相当于莱特兄弟的第一次飞行——虽然只飞了12秒，距离实用化还有很远，但证明了原理是可行的。

就在秀尔算法提出前后，量子计算的理论框架也在快速完善。1992年，大卫·德义奇和理查德·乔萨（Richard Jozsa）提出了第一个证明量子计算有明确加速优势的算法——德义奇-乔萨算法（Deutsch-Jozsa algorithm）：给定一个黑箱函数$f: \{0,1\}^n \to \{0,1\}$，判断它是常函数（所有输入输出相同）还是平衡函数（一半输入输出0，一半输入1），经典计算机在最坏情况下需要查询$2^{n-1}+1$次才能确定结果，而量子计算机仅需1次查询即可得到确定性答案[^4]。这个算法虽然没有什么实用价值，但它第一次从理论上证明了量子计算相对于经典计算的“量子优越性（quantum advantage）”：量子比特的叠加特性可以让量子计算机同时处理所有可能的输入，再通过量子干涉放大正确答案的概率，这是经典计算机根本无法实现的。德义奇后来回忆说：“我当时做这个算法的时候，只是想找一个足够简单的例子，让所有人都能明白量子计算到底快在哪里，没想到这个算法成了量子优越性的第一个‘标准测试题’。”

理论框架的另一块重要拼图是量子纠错。1994年秀尔提出算法的同时，也发现了量子比特的退相干（decoherence）问题：量子比特会和周围环境发生相互作用，在极短时间内失去叠加和纠缠特性，导致计算结果出错。当时IBM的物理学家罗尔夫·兰道尔（Rolf Landauer）和牛津大学的罗杰·彭罗斯（Roger Penrose）是量子计算最著名的质疑者，他们认为退相干问题根本无法解决，量子计算机永远不可能实现实用化。但1995年秀尔提出了第一个量子纠错码，证明了只要量子比特的错误率低于某个阈值，就可以通过冗余编码和集体测量抵消环境噪声的影响；1996年班韦特、丹尼尔·戈特斯曼（Daniel Gottesman）和伊曼纽尔·克尼尔（Emmanuel Knill）进一步提出了稳定子量子纠错码（stabilizer code）的通用框架，证明了容错量子计算（fault-tolerant quantum computing）在理论上是可行的[^5]。这个“量子纠错阈值定理”直接回应了所有质疑者的核心论点：退相干不是量子计算的“死穴”，只要物理量子比特的质量足够好，就可以通过编码实现逻辑量子比特的容错运算。

1990年代也是量子计算物理路线分化的关键时期。核磁共振（NMR）量子计算是第一个实现量子算法演示的路线，它的优势是量子比特的相干时间长（可以达到秒级），操控精度高，但缺点是量子比特数量很难扩展，因为NMR的磁场是全局的，无法精准操控单个分子的核自旋，所以NMR路线最多只演示到了7个量子比特，很快就被其他路线超越。离子阱路线是另一个早期的主流路线：1990年代，奥地利因斯布鲁克大学的莱纳·布拉特（Rainer Blatt）团队和美国国家标准与技术研究院（NIST）的大卫·维纳兰德（David Wineland）团队分别实现了囚禁离子的量子比特操控，他们用激光束精准操控单个钙离子或铍离子的内态自旋，相干时间可以达到秒级，2000年左右已经实现了10个离子量子比特的纠缠，这个路线现在的优势是量子比特的相干时间长、操控精度高，缺点是离子囚禁的系统扩展难度大，很难做到百比特以上的规模。超导电路路线是现在最接近实用化的路线：1999年，日本NTT的中村康徳（Yasunobu Nakamura）和米歇尔·德沃雷（Michel Devoret）团队首次提出了电荷量子比特的方案，用超导电路中的约瑟夫森结的非线性电感实现量子比特的能级，2002年，耶鲁大学的罗伯特·舍尔科普夫（Robert Schoelkopf）和德沃雷团队提出了传输子（transmon）量子比特的方案，通过增加超导电容抑制电荷噪声，将量子比特的相干时间从纳秒级提升到了微秒级，这个方案后来成了现在所有超导量子计算平台的基础，谷歌、IBM、阿里巴巴、本源量子的超导量子计算机都是基于传输子量子比特的。

2001年，IBM的团队再次用NMR平台演示了秀尔算法分解15的实验，这次的演示比1998年的版本更稳定，统计结果的信噪比更高，第一次让学界直观感受到了量子算法的加速优势。同一年，斯坦福大学的科学家用光子量子比特实现了德义奇-乔萨算法的演示，证明了光子也可以作为量子比特的载体。到2001年底，全球已经有超过50个大学和工业实验室建立了量子计算研究组，量子计算已经从少数物理学家的“个人爱好”变成了横跨物理、计算机、数学、密码学、材料科学的交叉学科。

费曼本人没有看到量子计算从理论到实验的突破：1998年，他在加州理工学院因癌症去世，享年79岁。他在1996年的一次访谈中说：“我当时在MI