3
2

Delete article

Deleted articles cannot be recovered.

Draft of this article would be also deleted.

Are you sure you want to delete this article?

More than 1 year has passed since last update.

理物Advent Calendar 2022

Day 2

一発芸としての第二法則

Last updated at Posted at 2022-12-01

飲み会で「え,〇〇物理学科なの?なんか見せてよ」と無茶ぶりされたときに,あたふたしたくないですよね.そんなとき一瞬で示せるとカッコいい定理を紹介します.皆さんもぜひ自己紹介のタイミングで示してみてください.
(かなり粗い理解をしている部分が多々あります.ご指摘などありましたらご連絡いただけると幸いです.)

セットアップ

熱浴$B$と系$S$が接している状況を考えます.全系のハミルトニアンは,

\hat{H}:=\hat{H}_S+\hat{H}_I+\hat{H}_B

とかけます.ただし,$\hat{H}_S,~\hat{H}_B$はそれぞれ系,熱浴のハミルトニアンで,$\hat{H}_I$は系と熱浴の相互作用を表すハミルトニアンです.
 初期条件は,$S$,$B$の間で相関がないような状態でかけます.

\hat{\rho}(0)=\hat{\rho}_s(0)\otimes\hat{\rho}_B^{can}

$S$系は任意の状態を持ってきてよくて,$B$の初期状態はカノニカル分布$\hat{\rho}_B^{can}=\frac{\exp(-\beta \hat{H}_B)}{Z}$とします.$Z$は分配関数.系全体はユニタリ発展するから,$\hat{\rho}(t)=\hat{U}\hat{\rho}(0)\hat{U}^{\dagger}$と書けます.このときの系,熱浴の状態は$\rho(t)$の部分トレースでかけます.つまり,

\hat{\rho}_B=\mathrm{Tr}_S\hat{\rho}(t),~~\hat{\rho}_S=\mathrm{Tr}_B\hat{\rho}(t)

となります.

やりたいこと

熱力学第二法則は,大学1年生で習う熱力学では,定理というよりかは原理として導入されるものです.その表現は様々で,例えばトムソンの原理やクラウジウスの原理などが挙げられますね.(なんかの論文に第二法則の同値な言い換えが21個くらい載ってた気がするけどどれだか忘れちゃったし,ホントに載ってたかも忘れちゃった......)
 要は,平衡熱力学じゃない上のセットアップから第二法則を示すのが今からやりたいことです.実際,これ以降導入するvon Neumannエントロピー

S( \hat{\rho}_S)=-\mathrm{Tr}[\hat{\rho}_S \ln \hat{\rho}_S]

は,密度行列の固有値に注目すれば形式的にShannonエントロピーと一致します.

定理と証明

系のvon Neumannエントロピー変化を$\Delta S_S=S(\hat{\rho}_S(t))-S(\hat{\rho}_S(0))$, 系の吸熱を$Q=-\mathrm{Tr}_B[\hat{H}_B(\hat{\rho}_B(t)-\hat{\rho}_B(0))]$とします.そもそも,熱浴は仕事をしないと思えば,(系の受け取るエネルギー)$=$(熱浴の出すエネルギー)$=$(熱浴のエネルギー減少)となることから理解できると思います.このとき,以下が成立します.

\Delta S_S-\beta Q\geq 0

これが熱力学第二法則に相当します.(要は系と熱浴のエントロピー変化の和が非負ってコト.)
以下,証明.初期条件でのエントロピーは,

\begin{align}
S(\hat{\rho}(0))&=-\mathrm{Tr_S}[\hat{\rho}_S(0) \ln \hat{\rho}_S(0)]-\mathrm{Tr_B}[\hat{\rho}_B^{can} \ln \hat{\rho}_B^{can}]\\
&=S_S(0)-\mathrm{Tr_B}[\hat{\rho}_B^{can}\ln\frac{\exp(-\beta \hat{H}_B)}{Z}]\\
&=S_S(0)+\beta\mathrm{Tr_B}[\hat{H}_B\hat{\rho}_B(0)]+\ln Z
\end{align}

とかけます.熱浴$B$の初期状態はギブス分布である($\hat{\rho}_B(0)=\hat{\rho}_B^{can}$)ことに注意してください.時刻$t$におけるエントロピーは,相対エントロピー

D(\hat{\rho}||\hat{\sigma})=\mathrm{Tr}[\hat{\rho}(\ln \hat{\rho}-\ln\hat{\sigma})]

の非負性を用いれば,

\begin{align}
S(\hat{\rho}(t))&=-\mathrm{Tr}[\hat{\rho}(t) \ln \hat{\rho}(t)]\\
&\leq-\mathrm{Tr}[\hat{\rho}(t) \ln (\hat{\rho}_S(t)\otimes\hat{\rho}_B^{can})]\\
&=-\mathrm{Tr}[\hat{\rho}(t)\ln(\hat{\rho}_S(t)\otimes\mathbb{I}_B)+\hat{\rho}(t)\ln(\mathbb{I}_S\otimes\hat{\rho}_B^{can})]\\
&=S_S(t)+\beta\mathrm{Tr}_B[\hat{H}_B\hat{\rho}_B(t)]+\ln Z
\end{align}

とできます.ここで,全系についてはユニタリー発展をしているので,

\begin{align}
S(\hat{\rho}(t))&=-\mathrm{Tr}[\hat{\rho}(t) \ln \hat{\rho}(t)]\\
&=-\mathrm{Tr}[U(\hat{\rho}(0) \ln \hat{\rho}(0))U^\dagger]\\
&=-\mathrm{Tr}[\hat{\rho}(0) \ln \hat{\rho}(0)]=S(\hat{\rho}(0))
\end{align}

となります.以上の結果を整理すると,求めたかった不等式が示せました.

注意事項

 証明を見たときかなりコンパクトだったので「一発芸になるやん!」と思ってこの記事を書きましたが,この記事の執筆をしているうちに説明に思ったより時間がかかることに気が付いたので,これは飲み会の一発芸には向かないかもしれません.また,以上はIyoda,Kaneko,and Sagawa(2017)より引用した証明になります(究極の出典はあるかもだけど). このお三方が同じ飲み会にいらしている場合は控えることを強くお勧めします.

参考文献

Iyoda, E., Kaneko, K., & Sagawa, T., Phys. Rev. Lett. 119, 100601(2017)

3
2
0

Register as a new user and use Qiita more conveniently

  1. You get articles that match your needs
  2. You can efficiently read back useful information
  3. You can use dark theme
What you can do with signing up
3
2

Delete article

Deleted articles cannot be recovered.

Draft of this article would be also deleted.

Are you sure you want to delete this article?