QuTiP(Quantum Toolbox in Python)是用于模拟与分析量子力学系统的 Python 库,同时处理封闭(幺正)与开放(耗散)量子系统,提供针对不同场景优化的多个求解器。安装:uv pip install qutip;可选包 qutip-qip(量子信息处理、线路与门)、qutip-qtrl(量子轨迹查看器)。
用途:研究主方程、Lindblad 动力学、退相干、量子光学、腔 QED 等。
主要能力:一、量子对象与态,如 basis、coherent、thermal_dm、destroy、num、Pauli 矩阵、tensor 张量积;二、时间演化求解器,sesolve(纯态幺正)、mesolve(混合态与耗散)、mcsolve(量子轨迹/蒙特卡洛)、brmesolve(弱系统-浴耦合)、fmmesolve(周期哈密顿量/Floquet);三、分析测量,包括 expect、entropy_vn、concurrence、fidelity、tracedist、关联函数与 steadystate;四、可视化,Bloch 球、Wigner 函数、Fock 分布、Hinton 图与矩阵直方图;五、高级方法,Floquet 理论、HEOM(非马尔可夫强耦合)、置换不变性(Dicke 态)、随机求解器、超算符与性能优化。
使用方式:导入 qutip 后构建哈密顿量与塌缩算符,选择合适求解器并传入初态与时间列表,用 e_ops 提取期望值或分析函数。例如 sesolve(H, psi0, tlist, e_ops=[num(N)])、mesolve(H, psi0, tlist, c_ops, e_ops=[num(N)])、mcsolve(..., ntraj=500)。
适用场景:阻尼谐振子、双量子比特纠缠动力学、Jaynes-Cummings 模型等开放量子系统仿真。
限制与提示:需截断希尔伯特空间;纯态优先用 sesolve 更快;时间依赖用字符串格式最快;仅存储所需数据用 e_ops;通过 Options 调整容差;内存不足时减小空间维度或使用 Krylov 方法;收敛性需变化 ntraj、空间尺寸与容差检查;量子门需要 qutip-qip。参考文档 https://qutip.readthedocs.io/ ,教程 https://qutip.org/qutip-tutorials/ ,GitHub https://github.com/qutip/qutip 。