# Demo 2：Gibbs 态制备

> 第 6 章 量子科学计算 · [本章演示](index.md) · [代码演示总索引](../index.md)

**可视化**：单子图（4 柱对比） —— 4 β 扫描的 P(|00⟩)=1/Z 反推

**对应教程**：[第 6 章 量子科学计算 中相关小节](https://tutorial.chenzhaoyun.com/quantum-tutorial/ch06-scientific-computing/) · 依赖安装：`pip install "unified-quantum[visualization,simulation]"`

## 实现要点

给定 H = diag(0, 1, 2, 3) 和逆温 beta，用 uniqc.rotation_prepare 直接将
Boltzmann 分布 |psi> = sum_j sqrt(e^{-beta E_j}/Z) |j> 装到 2 个 qubit 上。
扫描 beta 从低温到高温，展示不同温度下的概率分布与 H 矩阵本征结构。

## 完整源码

```{literalinclude} 02_gibbs_state.py
:language: python
:linenos:
:caption: demos/ch06-scientific-computing/02_gibbs_state.py
```

## 运行结果

![Demo 2：Gibbs 态制备](02_gibbs_state.png)

## 运行方式

```bash
# 仓库根目录执行；输出 PNG 到 demos/ch06-scientific-computing/02_gibbs_state.png，关键统计到 stdout
.venv/bin/python3 demos/ch06-scientific-computing/02_gibbs_state.py
```

---

[← 本章演示](index.md) · [↑ 代码演示总索引](../index.md) · [⌂ 教程门户](https://tutorial.chenzhaoyun.com/)