Satoshi Yukawa

Satoshi Yukawa

Department of Earth and Space Science, Osaka University

非平衡現象論

理学研究科宇宙地球科学専攻大学院科目, 大阪大学, 理学部, 2021

大学院科目、金曜日2限、理F202教室

最終課題

以下の二つの課題を提出せよ。
締め切り: 2月10日(厳守)
提出媒体:紙媒体もしくは電子媒体
紙媒体はF517号室前レポート入れ、電子媒体は yukawa あっと ess.sci.osaka-u.ac.jp まで。

  • Quizの中から、4問を選択し解答せよ
  • ゆらぎ、緩和、輸送、応答、ゆらぎの定理、情報熱力学など講義で出てきたキーワードを含む論文を読んでレポートするか、そのキーワードを含む自分の研究内容をレポートせよ。

スケジュール

10月1日
非平衡現象とは
平衡近傍でのゆらぎ: Boltzmann-Einsteinの原理。孤立系、開放系での表現

10月8日
平衡近傍でのゆらぎ: 最小仕事による表現、具体例
ゆらぎの動力学: 確率過程

10月15日
ゆらぎの動力学: 現象論的方程式、Onsagerの相反関係

10月22日
ゆらぎの動力学: Onsagerの最小散逸の原理、Onsager係数の決定、Onsager-Machlup過程(Langevin)

10月29日
ゆらぎの動力学: Onsager-Machlup過程(Fokker-Planck)、例: Brown運動

11月5日
(学祭のため休講)

11月12日
線形応答理論: 静的な場合、動的な場合(Liouville方程式、adjoint表現)

11月19日
線形応答理論: 動的な場合(摂動計算)、静的な場合との比較

11月26日
線形応答理論: 揺動散逸定理
ゆらぎの定理: 動機

12月3日
ゆらぎの定理: 能勢・Hoover動力学

12月10日
(休講)

12月17日
ゆらぎの定理: ゆらぎの定理、Jarzynski等式

12月24日
ゆらぎの定理: 剪断系における具体例

1月7日
情報と熱力学: Maxwellの悪魔、Szilard機関

1月14日
(大学共通テストのため休講)

1月21日
情報と熱力学: 情報量、KLダイバージェンス

1月28日
情報と熱力学: フィードバック系のゆらぎの定理

Quiz

Quiz 1(10月1日出題)
マクロ量 \( \{ a_i \} \)の平衡値 \( \{ \overline{a_i} \} \)からのずれは エントロピーで支配され \[ P(\{a_i\}) = (2\pi)^{-n/2} \left[ \det \left( -\dfrac{1}{k_\mathrm{B}} \dfrac{\partial^2 S}{\partial a_i \partial a_j } \right) \right]^{1/2} \exp \left( \dfrac{1}{2 k_\mathrm{B}} \sum_{ij} \dfrac{\partial^2 S}{\partial a_i \partial a_j} (a_i - \overline{a_i}) (a_j - \overline{a_j}) \right) \] と書ける。この多次元Gauss積分を実行し、規格化定数がこれで良いことを確認せよ。

Quiz 2(10月8日出題)
マクロ量 \( \{ x_i \} \)の分布 \[ P(\{x_i\}) \propto \exp \left( -\dfrac{1}{2} \sum_{ij} \beta_{ij} x_i x_j \right) \] に対し、期待値 \( \left\langle x_i x_j \right\rangle = \beta_{ij}^{-1}\)となることを示せ。また \(X_i = \sum_j \beta_{ij} x_j\)と置いたとき、 \( \left\langle X_i x_j \right\rangle = \delta_{ij} \)、 \( \left\langle X_i X_j \right\rangle = \beta_{ij} \)となることを示せ。

Quiz 3(10月15日出題)
Hall効果をDrudeモデルにより考察する。 電荷\(q\)をもつ質量\(m\)の荷電粒子が運動方程式 \[ m \dfrac{d\mathbf{v}}{dt} = q \mathbf{E} + q \mathbf{v} \times \mathbf{B}- \gamma \mathbf{v} \] に従うとする。ここで\(\mathbf{E}, \mathbf{B}\)はそれぞれ電場、磁場、\(\gamma\)は散逸を表す。 電場、磁場が\(\mathbf{E} = (E_x, E_y, 0) , \mathbf{B} = (0,0,B)\)のとき、 定常速度\( (v_x, v_y) \)を求め、電気伝導率テンソル\( \sigma_{xx}, \sigma_{xy},\sigma_{yx},\sigma_{yy} \)を \[ v_x = \sigma_{xx} E_x + \sigma_{xy} E_y \] \[ v_y = \sigma_{yx} E_x + \sigma_{yy} E_y \] で定義したとき、電気伝導率テンソル\( \sigma_{xy} , \sigma_{yx} \)に関してOnsagerの相反定理が成立していることを確認せよ。

Quiz 4(10月22日出題)
講義ではマクロ物理量 \( \{ a_i \} \)の時間発展を条件付き確率 \( W_{\Delta t} (a_i+\Delta a_i | a_i) \)で表した。この条件付き確率が詳細つりあい \[ W_{\Delta t} (a_i+\Delta a_i | a_i) P_{eq}(a_i) = W_{\Delta t} (a_i | a_i+\Delta a_i) P_{eq}(a_i+\Delta a_i) \] を満たすことを示せ。

Quiz ５(10月29日出題)
ブラウン運動のLangevin表現、 \[ \dfrac{dp_i}{dt} = - \dfrac{f}{T_R M} p_i + \xi_i(t) \] ただし\( \langle \xi_i(t) \rangle_{eq} = 0 \enspace , \quad \langle \xi_i(t) \xi_j (t’) \rangle_{eq} = 2 k_\mathrm{B} f \delta_{ij} \delta (t-t’) \)に対して、 初期条件\( p_i(0) = p_i^0 \) のもと、\( p_i(t) \)の解を求めよ。 これを使い、 \[ \langle p_i (t+\tau) p_j (t) \rangle_{p_i(0)=p_i^0, p_j(0)=p_j^0} \] を計算せよ。また、初期値 \( p_i^0, p_j^0 \)を熱平衡にとり十分に時間が経った後の、 \[ \langle p_i (t+\tau) p_j (t) \rangle_{eq} - \langle p_i (t+\tau) \rangle_{eq} \langle p_j (t) \rangle_{eq} \] を計算せよ。

Quiz 6(11月12日出題)
講義ではハミルトニアンが無摂動状態\( H_0\)から外場\(E\)のもと\( H_0- M E \)と変化する時の物理量\(A \)の応答を見た。 逆温度で摂動を\( E = 1- \beta’ / \beta \)と与えたとき、\(A=M=H_0 \)として熱容量を求めよ。

Quiz 7(11月19日出題)
ハミルトニアン\(H\)とのPoisson括弧で定義される\( \mathcal{L}(\cdot) =\{ \cdot, H\}\) という線形微分演算子に対して、物理量\( M \)に対し \[ \int_0^t ds e^{s \mathcal{L}} \mathcal{L} M = e^{t \mathcal{L}} M - M \] となることを示せ。

Quiz 8(11月26日出題)
複素アドミッタンス \[ \chi_{AM}(\omega) = \lim_{\epsilon \downarrow 0}\int_{0}^{\infty} ds \phi_{AM}(s) e^{i\omega s} e^{-\epsilon s} \] に対し、Kramers-Kronig関係式 \[ \Re \chi_{AM}(\omega) = \pv\int_{-\infty}^{\infty} \dfrac{d\omega’}{\pi} \Im \dfrac{\chi_{AM}(\omega’)}{\omega’-\omega} \] \[ \Im \chi_{AM}(\omega) = - \pv\int_{-\infty}^{\infty} \dfrac{d\omega’}{\pi} \Re \dfrac{\chi_{AM}(\omega’)}{\omega’-\omega} \] が成立することを示せ。

Quiz 9(12月3日出題) 位相空間中のある領域\( \mathcal{D}(t) \)の体積のLagrange的な時間変化が、 \[ V(\Gamma(t+\Delta t)) = \int_{\mathcal{D}(t)} d\Gamma(t) (1+\Lambda (\Gamma(t)\Delta t)+O(\Delta t^2) \] となることを示せ。ただし \( \Lambda (\Gamma(t)) = \sum_i \left(\dfrac{\partial \dot{q}_i }{\partial q_i } + \dfrac{\partial \dot{p}_i} {\partial p_i} \right) \) である。

Quiz 10(12月17日出題)
運動量\( p\)、座標\( q\)に対し \[ \dfrac{dp}{dt} = 0 \enspace, \quad \dfrac{dq}{dt}= p \] という運動を考える。

  • この運動が可逆であることを示せ。
  • ある初期状態から時間\( t\)だけ時間発展させ、状態を時間反転し、おなじ時間\( t \)だけ時間発展し、時間反転するともとの初期状態に戻ることを示せ。 また \[ \dfrac{dp}{dt} = -p \enspace, \quad \dfrac{dq}{dt}= p \] という運動に対し、
  • この運動が可逆でないことを示せ。
  • ある初期状態から時間\( t\)だけ時間発展させ、状態を時間反転し、おなじ時間\( t \)だけ時間発展し、時間反転してももとの初期状態に戻らないことを示せ。

Quiz 11(12月24日出題)
\[ \sigma_{xy} = \dfrac{1}{V}\left( \sum_i \dfrac{p_x^i p_y^i}{m} + \sum_i q_y^i f_x^i \right) \] が応力テンソルの\( xy \)成分の表現である事を示せ。

Quiz 12(1月7日出題)
Szilardエンジンにおいて、系を\( \alpha : 1-\alpha \)で分割した場合、得られる仕事の期待値が \[ W= k_\mathrm{B} T \left( \alpha \log \dfrac{1}{\alpha} + (1-\alpha) \log \dfrac{1}{1-\alpha} \right) \] となることを示せ。

Quiz 13(1月21日出題)
系X、系Yの同時分布\(P_{X\otimes Y} \)に対する情報量 \[ S(X,Y) = -\sum_{x,y} P_{X\otimes Y} (x,y) \log P_{X\otimes Y} (x,y) \] が \[ S(X,Y)= S(Y)+S(X|Y) \] と書けることを示せ。ここで\( S(X|Y) \)は条件付き情報量 \[ S(X|Y)=- \sum_y P_Y(y) \left( \sum_x P_X(x|y) \log P_X(x|y) \right ) \] である。

評価

期末レポート、およびQuizの解答状況(数問選んで、最後に提出)をみて総合的に判断する。 詳細はシラバスを参照してください。

参考書、参考文献

一般的なこと

  • S. R. de Groot and P. Mazur “Non-equilibrium Thermodynamics”, Dover Publications
  • 一柳正和「不可逆過程の物理」日本評論社
  • 川崎恭治「非平衡と相転移-メソスケールの統計力学-」朝倉書店
  • 早川尚男「非平衡統計力学」サイエンス社
  • 北原和夫「非平衡系の統計力学」岩波書店
  • 鈴木増雄「統計力学、岩波講座現代の物理学」岩波書店
  • 戸田盛和、久保亮五、斎藤信彦、橋爪夏樹「統計物理学、現代物理学の基礎第二版」岩波書店
  • 関本謙「ゆらぎのエネルギー論」岩波書店
  • D. N. Zubarev “Nonequilibrium Statistical Thermodynamics”, Consultants Bureau
  • R. Zwanzig “Nonequilibrium Statistical Mechanics”, Oxford
  • M. Ichiyanagi, “Conceptual developments of non-equilibrium statistical mechanics in the early days of Japan”, Phys. Rep. 262 (1995) 227.
  • R. Kubo, M. Toda, and N. Hashitsume, “Statistical Physics II, Nonequilibirum Statistical Mechanics” Spinger-Verlag
  • J. A. McLennan, “Introduction to Non-equilibrium Statistical Mechanics” Prentice-Hall
  • 佐々真一「非平衡現象論」講義ノート

Einsteinのゆらぎの理論

  • L. Landau, E. Lifshitz 「統計物理学 (下)」岩波書店, L. Landau, E. Lifshitz, “Statistical Physics (course of theoretical physics volume 5) Butterworth-Heinemann”
  • A. Einstein, “Theorie der Opaleszenz von homogenen Flüssigkeiten und Flüssigkeitsgemischen in der Nähe des Kritischen Zustandes”, Ann. der Phys. 33 (1910), pp. 1275-1298. 邦訳 アインシュタイン選集1、監修 湯川秀樹 翻訳 井上健、谷川安孝、中村誠太郎 共立出版(1971)

確率過程

  • N. G. van Kampen, “Stochastic Processes in Physics and Chemistry”, Elsevier

ブラウン運動

  • A. Einstein “Investigations on the theory of the Brownian movement”, Dover

Onsagerの相反定理、Onsager-Machlup過程

線形応答

  • 戸田盛和、久保亮五、斎藤信彦、橋爪夏樹「統計物理学、現代物理学の基礎第二版」岩波書店
  • R. Kubo, M. Toda, and N. Hashitsume, “Statistical Physics II, Nonequilibirum Statistical Mechanics” Spinger-Verlag

ゆらぎの定理、Jarzynski等式

情報と熱力学

  • H. Leff and A. F. Rex ed. “Maxwell’s Demon 2 Entropy, Classical and Quantum Information, Computing”, Institute of Physics Publishing

大偏差性