深入探索:参数平移法则 (PSR) 的证明与适用范围
在正文中,我们介绍了参数平移法则(Parameter-Shift Rule, PSR):只需在两个”平移后”的电路上测量期望值,就能得到关于门参数的解析梯度(analytical gradient)。本词条将给出这一法则的完整证明,并严格刻画它的适用范围。整个推导只用到矩阵欧拉公式(其推导见《含参量子门》)和几个基本的三角恒等式,我们对每一个等号都注明所用的依据,读者可以逐步核验。
1. 问题的设定:把一般电路化简为单个参数门
考虑一个含参数的量子电路,其中只有一个门 依赖于参数 ,其余的门与测量都与 无关。为讨论方便,我们把 无关的门按位置归并到电路的前段和后段,即把整个电路写成 ,其中 和 都是与 无关的幺正算符。设输入态为 、可观测量(observable)为 ,则期望值为
其中第二个等号只是把 吸收进初态、把 吸收进观测量。具体地,代入 并展开 ,得
因此只需取 、并把 重新记作 (它仍是厄米算符),就得到了最简形式。这一化简表明:无论参数门位于电路的什么位置、周围还有多少别的门,我们都只需要研究 这一种情形。
参数门本身具有如下形式:
其中 是厄米的生成元(generator)。对于旋转门 ,生成元分别是对应的泡利矩阵 。这些生成元满足一个关键性质:(单位矩阵)。
2. 核心引理:期望值是 的正弦型函数
引理(正弦响应) 设 与 都是厄米算符且 , 为任意归一化初态。记 , 为对易子(commutator)。定义 ,其中 。则存在与 无关的实常数 ,使得
其中
证明 分三步进行。
第一步:把 展开成 与 的线性组合。 由 和矩阵欧拉公式(含参量子门)可得
对其取伴随(利用 以及 、 取实数值)得
为书写简洁,记 ,。
第二步:展开 并对初态取期望。 将上面两式代入 ,逐项相乘可得四项:
四项的来源分别是:;;;,最后一步用了 。把中间两项合并为对易子,得到
于是
第三步:用半角恒等式把 的函数改写为 的函数。 所需的三个恒等式为
它们都可由倍角公式 与 直接得到。先处理含 、 的两项:
再处理交叉项:
把三部分按常数项、 项、 项归拢,即得
其中 、、。
最后我们说明这三个系数都是实数,且与 无关。先看实数性: 是厄米的,而 的伴随为 ,故它也是厄米的,因此 与 都是实数,从而 。再看 :对易子满足
即 是反厄米的,因此 是厄米算符,其期望值 是实数,故 。至于与 无关这一点: 的构造只涉及固定初态 上固定算符 、、 的期望值,不含任何 。证毕。
推论(振幅–相位形式) 实系数的正弦与余弦的线性组合总可以合并成一个带相位的正弦波。取
则由加法公式 展开可得
这正是正文中所说”期望值必然是正弦函数”的准确含义:振幅 、相位 和垂直偏移 完全由初态 、生成元 和观测量 决定,与 无关;随 变化的只有 这一个因子。特别地,当 (即 )时 是常数,梯度恒为零,参数无法被训练——这与”平坦损失景观”的现象一致。
3. 参数平移法则的证明
定理(参数平移法则) 在引理的条件下,
证明 对 逐项求导,得
再计算平移后的两个函数值。由诱导公式 、、、,得
两式相减时,常数项 互相抵消:
证毕。用振幅–相位形式来读,这个证明说的是:,而 ,两者恰好吻合。PSR 之所以给出的是解析梯度而非数值近似,正是因为它利用了”期望值是频率为 的正弦波”这一先验结构:对任意正弦型函数,其在某点的导数都等于它向前、向后各平移 后的函数值之差的一半。
4. 数值验证:一个完整的算例
下面我们构造一个可以手工算透的例子,把引理中的常数逐一算出,并在具体的一点上验证 PSR 两边相等。取生成元 、观测量 ,初态为
该态满足归一化条件,因为 。
先计算引理所需的初态期望值。对形如 的一般单比特态,直接把泡利矩阵与态向量相乘可以验证:
本例中 、,于是
再利用两个泡利恒等式 与 (两者都可以用矩阵直接相乘验证),代入引理的系数公式:
因此期望值函数为
其中振幅与相位为
作为一条独立的交叉验证,我们绕开引理、直接化简海森堡绘景中的算符 :由第 2 节第二步的展开式,
(其中 ),于是 ,与引理给出的结果一致。
现在取 验证 PSR。左边是解析导数:
右边是两次平移后的期望值之差的一半。利用 、、、,以及 、,得
两边在保留的四位小数内完全相等,PSR 在该点得到验证。顺带一提,,也可以用振幅–相位形式 复核,说明推论的形式同样正确。
5. PSR 的适用范围:哪些门可以用?
适用的条件:生成元恰有两个不同本征值
回顾第 2、3 节的证明,我们只在两处用到了生成元的性质:其一,把 展开为 需要 ,等价于 的谱(本征值集合)为 ;其二,半角换元之后函数的频率恰好为 。把这两点一并推广,就得到最一般的适用条件。
定理(两本征值生成元的广义参数平移) 设 ,其中 是厄米算符且恰有两个不同本征值 。记谱隙(spectral gap),并设电路的其余部分与 无关。则
证明 记 、,并定义 。由于 厄米, 也厄米; 的两个本征值为 ,由谱分解可知 。又因单位矩阵与一切算符对易,指数可以拆开:
式中的标量相位 在 中与自身的复共轭相乘而抵消(,其中 ),所以 只通过组合 依赖于 。把第 2 节的引理应用于生成元 、参数 ,得
其中 与 无关。对上式求导,并对任意平移量 用和差化积公式 、,得
只要 ,两式相除便得 。取 ,即 ,系数化为 。证毕。
推论(谱为 的情形) 若 的本征值为 (),则 ,平移量为 、系数为 :
泡利生成元对应 ,此时平移量为 、系数为 ,即标准的 PSR。常见的适用例子包括:
- 单比特旋转门 、、,生成元谱均为 ,直接用标准 PSR。
- 两比特相互作用门 :生成元 的本征值为 ,标准 PSR 直接适用。
- 加权单比特生成元 ( 为不全为零的实数):生成元是厄米无迹的 矩阵,其本征值为 ,仍属两本征值情形,取 即可用上式求梯度。
注意参数化约定。 平移量的具体数值依赖于我们把门写成 还是 。若采用后一种约定且谱为 ,则有效参数为 ,同样的推导给出平移量 、系数 。两种写法在文献中都很常见,对照不同资料时务必先统一参数化约定;本词条统一采用 。
无法直接应用 PSR 的情况
以下类型的参数门不能直接套用标准(两项)PSR,需要额外处理。
-
生成元有三个及以上不同本征值。 例如三比特门 :把 对 枚举取值,可得生成元的谱为 ,其两两之差有 和 两个值,因此 同时含有频率 和频率 的正弦成分,任何单一平移量的两项公式都无法同时匹配两个频率。再如受控旋转 ,其生成元为 ,谱为 ,同样不满足两本征值条件。对这类门,标准做法是先做电路分解(例如 可以分解为两个 CNOT 与若干个 ,再对每个 分别使用 PSR),或使用下面提到的推广形式。
-
参数化操作不是幺正的。 本词条的推导前提是参数门为幺正演化。需要说明的是,若噪声信道以与 无关的方式作用在参数门之外,则可以把它并入等效初态或等效观测量(利用信道的伴随 ),PSR 依然成立;但若噪声进入了参数门内部(例如把参数门整体替换为含参数的噪声信道),幺正性被破坏,期望值一般不再是有限个正弦分量之和,需要另行分析。
-
参数以非线性方式进入。 若门的形式为 ,其中 是非线性函数(例如 ),则 ,其中 是引理给出的正弦型函数。此时不能对 直接套用 PSR,但只要 可微且 已知,链式法则给出
即先对中间变量 使用 PSR,再乘上 。
-
同一参数出现在多个门中(共享参数)。 设 同时出现在两个门中。把 看成二元函数 ,其中 只进入第一个门、 只进入第二个门,则链式法则给出 (在 处取值)。固定另一个变元时,每个偏导数都对应单个被夹在中间的参数门,可分别用两项 PSR 计算。结论是:梯度等于”每次只平移一个出现位置、各做一次两项差分、再求和”。
总结表
| 门类型 | 是否可直接用 PSR | 原因 |
|---|---|---|
Rx(θ), Ry(θ), Rz(θ) | ✅ | 生成元为泡利矩阵,谱为 ,平移 ,系数 |
| ✅ | 谱为 ,,平移 ,系数 | |
| ✅ | 生成元谱为 | |
CRz(θ)(受控旋转) | ⚠️ | 生成元谱为 ;需分解为 CNOT+Rz 或用推广公式 |
| ❌ | 生成元谱 ,谱隙不唯一,期望值含多个频率 | |
| ❌ | 参数非线性进入;可用链式法则补救 | |
| 含噪声通道的参数门 | ⚠️ | 视噪声位置而定:门外的 无关噪声不影响 PSR,门内噪声破坏推导前提 |
概括地说,无法直接应用参数平移法则的参数门,是那些生成元不恰有两个本征值(例如谱为 或 )、参数以非线性方式进入门的定义、或者参数化操作本身非幺正的门;其余情形(包括一切谱为 的生成元)都有精确的两项平移公式可用。
好消息有两点。第一,绝大多数量子机器学习模型使用的参数化门(如 )都属于可直接使用 PSR 的类型。第二,对于更复杂的情况,学术界已经发展出推广的参数平移法则(generalized parameter-shift rules),它们可能需要更多次平移或不同的平移量,但基本思想与本词条相同;像 PennyLane 这样的框架会自动识别生成元的谱结构,并调用合适的梯度计算方法。
总结:实践者的启示
- PSR 是解析梯度,不是数值近似:它的精确性源于量子门 的数学结构——期望值关于 是频率为由生成元谱隙决定的正弦波,这保证了梯度计算的可靠性。
- 检查你的门:当你设计自己的含参数电路时,应确认所用参数门的生成元恰有两个本征值(绝大多数常见门都满足),否则需要分解电路或改用推广公式。
- 信任你的框架:在使用 PennyLane 等工具时,通常不需要手动实现 PSR。框架会自动处理梯度计算,你只需要专注于模型的设计与训练。
通过本词条的推导,我们不仅知道了如何使用 PSR,更理解了它为何精确、以及它的边界在哪里。
参考资料
- 本词条整理自《量子计算算法教程》原文:ch02-quantum-nn/parameter-shift-appendix
- 内容遵循 CC BY-NC-SA 4.0 许可协议