用途:QuTiP 是 Python 量子工具箱,用于模拟和分析量子力学系统,覆盖封闭(幺正)与开放(耗散)系统,提供针对不同场景优化的多种求解器。
主要能力:
- 量子对象与态:创建 Fock 态 basis(N,n)、相干态 coherent(N,alpha)、热密度矩阵 thermal_dm,以及湮灭算符 destroy、粒子数算符 num、泡利矩阵等,并可用 tensor 构建复合系统。
- 时间演化:sesolve 处理纯态幺正演化;mesolve 处理混合态与耗散;mcsolve 做量子轨迹(蒙特卡洛);另有 brmesolve(弱系统-浴耦合)、fmmesolve(周期哈密顿量/Floquet)。
- 分析与测量:expect 求期望值,entropy_vn 计算冯诺依曼熵,concurrence 度量纠缠,fidelity、tracedist 比较态,correlation_2op_1t 与 spectrum_correlation_fft 计算关联函数与谱,steadystate 求稳态。
- 可视化:Bloch 布洛赫球、wigner 维格纳函数、plot_fock_distribution、hinton 与 matrix_histogram 矩阵可视化。
- 高级方法:Floquet 理论、HEOM 非马尔可夫强耦合、置换不变性(Dicke 态、集体算符)等。
使用方式:安装 uv pip install qutip;量子信息处理(线路、门)需 uv pip install qutip-qip;量子轨迹查看器需 uv pip install qutip-qtrl。典型流程为定义哈密顿量与塌缩算符、设定初态与时间列表、调用求解器并绘图。
适用场景:量子光学(如阻尼谐振子、Jaynes-Cummings 模型)、量子比特纠缠动力学、开放量子系统耗散研究等。
限制与提示:需截断希尔伯特空间;纯态优先用 sesolve 更快;时间依赖用字符串格式最快;用 e_ops 只存所需数据;通过 Options 调整容差;mcsolve 可并行。内存不足可减小维度或用 Krylov 方法;数值不稳定可减小步长或提高容差。文档见 https://qutip.readthedocs.io/,教程见 https://qutip.org/qutip-tutorials/,源码见 https://github.com/qutip/qutip。