如何在Qiskit中构造任意密度矩阵？以ρ=[7/9,0;0,2/9]为例

为什么幺正变换无法构造目标混合态 #

你提到的目标密度矩阵 $\rho=\begin{bmatrix}7/9 & 0 \ 0 & 2/9\end{bmatrix}$ 是混合态，它的秩为2（存在两个非零特征值）。而幺正变换作用在纯态 $|0\rangle\langle0|$ 上得到的仍然是纯态（秩为1的密度矩阵），因此确实不存在这样的幺正矩阵 $U$ 实现转换。

Qiskit中构造任意密度矩阵的简单方法 #

方法1：直接使用DensityMatrix类构造 #

这是最直接的方式，直接传入密度矩阵的数值即可：

from qiskit.quantum_info import DensityMatrix
import numpy as np

# 定义目标密度矩阵
rho = np.array([[7/9, 0], [0, 2/9]])
# 构造DensityMatrix对象
dm = DensityMatrix(rho)

# 验证结果
print(dm)

你可以通过dm.data查看底层的numpy数组，确认构造结果是否符合预期。

方法2：通过混合纯态的方式构造 #

混合态本质是纯态的概率加权和，比如 $\rho = p|0\rangle\langle0| + (1-p)|1\rangle\langle1|$，这里 $p=7/9$。你可以先构造对应的纯态，再按概率混合：

from qiskit.quantum_info import Statevector, DensityMatrix
import numpy as np

# 定义纯态和对应概率
state0 = Statevector([1, 0])
state1 = Statevector([0, 1])
probabilities = [7/9, 2/9]

# 按概率混合得到密度矩阵
dm = sum(p * DensityMatrix(state) for state, p in zip([state0, state1], probabilities))

print(dm)

方法3：在量子电路中模拟混合态 #

如果需要在量子电路上下文使用这个混合态，可以将构造好的密度矩阵作为初始状态传入模拟器：

from qiskit import QuantumCircuit
from qiskit.quantum_info import DensityMatrix
from qiskit_aer import AerSimulator

# 构造目标密度矩阵
rho = np.array([[7/9, 0], [0, 2/9]])
dm = DensityMatrix(rho)

# 创建量子电路（这里以空电路为例，你可以添加任意门操作）
qc = QuantumCircuit(1)

# 用模拟器运行电路，传入初始密度矩阵
simulator = AerSimulator()
result = simulator.run(qc, initial_state=dm).result()
final_dm = result.data()['density_matrix']

print(final_dm)

总结 #

对于任意密度矩阵，方法1是最高效直接的选择；如果需要基于纯态的物理意义构造，方法2更直观；如果要结合量子电路的演化，方法3可以将密度矩阵作为初始状态接入模拟器。

内容的提问来源于stack exchange，提问作者윤성원