エントロピー最大原理の非平衡への拡張
ー ミクロ力学からの構成と実験への予言
佐々 真一 （京大 理）（共同研究者は最後に紹介します）
熱場の量子論とその応用@ 京都

2019/09/03
1

PART I

変分原理の誘惑と罠

2

魅惑のポイント１「美しい」
例：最小作用の原理
力学：微分方程式による局所的な記述

微分方程式の解を大域的に特徴づける
見通しがよくなる、計算が簡単になる、構造（対称性など）が分かりやすくなる・・

3

魅惑のポイント２「新しい体系の創出」
例：
期待値
量子力学

最小作用の原理

古典力学

経路積分
変分構造が未知の体系を示唆

4

非平衡系における変分原理－ 色々 －
最小散逸原理最小発熱原理エントロピー生成最小原理
最大流速原理過剰エントロピー生成最小原理
大偏差原理

相加性原理

5

綺麗なバラにはとげがある6

綺麗なバラにはとげがある
標準的な批判「既知の結果を言い直しただけ。」「実験への予言がない。」「その世界で閉じていて、他と関係ない。」
・・・・「物理じゃない！」（無意味語・・）

形式的な議論で、「特定の専門家が」好きなのは分かるが、「科学全体の発展」には関係ないことも多い。

7

変分原理の提案（普通の）研究への要請「既知の理論で導出されていない結果を含む。」「実験への新しい予言を与える。」「他の理論との新しい関係が見える。」
・・・・「物理！」（無意味語・・）

8

現象例

界面温度は？
9現象例：簡単な考察
界面位置を指定
熱流が空間的に一定 ⇒ 温度場
圧力が空間的に一定 ⇒ 密度場
界面位置を決める原理？
例：化学ポテンシャルの連続性⇒ 界面温度＝平衡転移温度

LiquidGas

平衡定積相共存
化学ポテンシャル連続 ⇔ 圧力一定

10

講演のメッセージ
熱力学変分原理 - 復習
非平衡熱力学変分原理 - 現代的概説
オンサーガ原理、エントロピー生成最小原理

新しい変分原理 – 最近の研究から「修正エントロピー」最大原理
11

セミマクロからの導出と実験への予言（５分）

PART II

熱力学変分原理とその拡張
12平衡系の変分原理１（エンタルピー）

例1
エントロピー関数
変分関数

変分公式（エントロピー最大原理）

等圧断熱環境で平衡値 を決める

熱力学第２法則から証明

13

平衡系の変分原理2例２自由エネルギー関数

変分関数

変分公式（自由エネルギー最小原理）

等圧等温環境で平衡値 を決める

熱力学第２法則から証明

14変分原理からの発展１
例1変分公式（エントロピー最大原理）

ゆらぎ公式（Einstein）

等圧断熱環境でのゆらぎと応答の関係

ミクロカノニカル分布の示唆

15変分原理からの発展２
例2変分公式（自由エネルギー最小原理）

ゆらぎ公式（Einstein）

等温等圧環境でのゆらぎと応答の関係

カノニカル分布(T-p分布）の示唆

16

Einstein の地図
熱力学第２法則
導出

統計力学
「導出」

変分原理
示唆

ゆらぎの公式

新しい現象の予言

示唆

確率分布

現象の新しい理解
17

変分原理の拡張１（エンタルピー）

例1
エントロピー関数
変分関数

定常場 for変分公式（エントロピー最大原理）

等圧断熱環境で平衡値 を決める

線形応答領域：現象論で導出
18

変分原理の拡張2例２変分関数と大域温度

定常場 for変分公式（自由エネルギー最小原理）

線形応答領域：現象論から導出！

2温度環境で平衡値 を決める

Section 3 in [Nakagawa and Sasa,
Global thermodynamics for heat
Conduction systems,
https://arxiv.org/abs/1906.07022]19非自明な現象へ

熱力学の拡張大域熱力学 ！
自由エネルギー最小原理
Nakagawa and Sasa,
Global thermodynamics for heat
Conduction systems,
https://arxiv.org/abs/1906.07022

メソスケールダイナミクスの解析拡張された変分原理の導出「修正エントロピー」最大原理

Sasa, Nakagawa, Itami, Nakayama,
Stochastic order parameter dynamics
for phase coexistence in heat conduction,
https://arxiv.org/abs/1908.03029

20

本講演のターゲット

Order-disorder 転移（外場０）
平衡では１次転移境界で熱流制御
空間３次元
21（１２分）PART III
非平衡熱力学の視点からのメソスケールモデル
オンサーガ―原理とエントロピー生成最小原理

22

メソスケールモデルの作り方
1. 変数の選択 － slow variables の完全な組
2. Onsager 原理 － 第２法則および可逆性との整合性
3. 非平衡性 － 境界条件のみ（仮定）4. 非平衡定常状態 － エントロピー生成最小原理

23

もっとも簡単な例（＝定数）slow variable
熱力学法則からの制約(決定論）:

B >0 :オンサーガ―定数
可逆性からの制約（軌道のゆらぎ）:24

注釈
多変数系に一般化
オンサーガ係数が正値対称行列
オンサーガの相反定理
熱力学的力を固定して得られる定常状態
最小散逸原理などなど・・・
25

応用例：１次元単純熱伝導
Slow variable

Internal energy density field

entropy density as a function of

Onsager principle

非平衡境界条件：case I
熱流制御境界条件

非平衡境界条件：case II
温度制御境界条件

単純熱伝導＝内部エネルギー密度のみが遅い変数になっている熱伝導26定常状態：１次元単純熱伝導
エントロピー生成率
for

決定論力学
温度制御境界条件

エントロピー生成最小原理
熱流制御境界条件では不成立
Onsager係数が物理量に依存すると不成立

色々な提案・・
27

２１世紀の理解：前提
確率分布

Local detailed balance(Onsager principle)
28

２１世紀の理解：定理
線形応答領域

定理１エントロピー生成最小原理

定理２ ズバーレフ＝マクレナン表現
過剰エントロピー生成
29

以上を踏まえて、熱伝導下相共存のメソスケールモデル
30

完全な熱力学変数の組
order parameter density field

momentum density fieldenergy density fieldinternal energy density field

オーダーパラメータの運動エネルギー

界面熱力学
gradient term

Entropy functional
31
局所熱力学
entropy density

界面熱力学
gradient term

１次転移のLandau 理論

forfor32

エントロピー密度定数33Onsager 原理
可逆項34

超絶微妙な話

ミクロ長さ

界面厚み
系の大きさ
（粗視化スケール）

III C 節参照：
Sasa, Nakagawa, Itami, Nakayama,
Stochastic order parameter dynamics
for phase coexistence in heat conduction,
https://arxiv.org/abs/1908.03029

35関連する先行研究のモデル
Penrose and Fife, Physica D, 1990;
over-damped version, without noise, andHalperin, Hohenberg, and Ma, Phys. Rev. B, 1974;
over-damped version, near the critical, andBedeaux, Johannessen, and R{¥o}sjorde, Physica A, 2003:
hydrodynamics with gradient terms, without noise
Onuki, Phys. Rev. E, 2007;
hydrodynamics with gradient terms, without noise
Fujitani, J. Phys. Soc. Jpn., 2010;
Analysis of hydrodynamics with gradient terms, without noise36熱伝導下相共存のメソモデル

境界条件
37（エネルギー保存）

スケーリング、無次元化: microscopic length ～: システムサイズ ～: coarse graining length

Make dimensionless quantities with appropriate estimationeg.

38
weak noise and thin interface limitもっとも大事なこと
マクロな時定数で無次元化したときのノイズ強度はマクロとミクロの体積比 に比例する － 弱ノイズ極限
無次元化した界面の厚みが であることから、界面運動の無次元化した時定数はになる - 薄界面極限

通常の流体力学極限の処方箋（時空間スケーリング領域）では、「界面運動」が記述されない。

ミクロとマクロのスケール分離比をパラメータにいれて、弱ノイズ極限と薄界面極限を「同時に」とる必要がある。
39

ズバーレフ＝マクレナン表現

修正エントロピー

過剰エントロピー生成

40
界面の運動の時定数（３０分）

PART IV

解析と結果
修正エントロピー最大原理の導出と実験への予言
41

平衡系 （断熱環境）

The equilibrium profile:

エントロピー最大原理
maximizes
withorderdisorder42熱伝導系 （熱流制御）
The steady state profile: 「修正エントロピー」最大原理
maximizes
with
界面位置

界面温度
43

この寄与を計算！

単一界面配置「界面の位置」 がXにある配置

界面領域

44

単一界面配置に対する過剰エントロピー生成を評価する

45の評価界面位置given
熱力学量緩和の時定数

界面運動の時定数

決定論的力学の時空積分で評価

空間x 時間界面内ノイズの影響を無視できない

46

の評価

秩序相での値（定数）47過剰エントロピー生成
バルクの寄与のまとめ

48
界面の位置X に対するquasi-equilibrium 状態で決まる

の評価

各バルクの端点での温度差（決定論ダイナミクス）で決まる
時間領域の積分しか効かない

疑似平衡状態を保ったまま界面がゆっくり運動するときの界面でのの温度ギャップがエントロピー生成を決める

49

参考：実験！

50

温度ギャップの評価
疑似平衡状態における各バルク領域での温度プロファイルの界面運動による摂動的補正
全エネルギーが一定
潜熱の界面領域からバルクへの配分がオンサーガ―原理で決定

51

温度ギャップの結果
界面が端にあるときに温度ギャップが消える条件

52

過剰エントロピー生成
結果53
修正エントロピー最大原理を使って、物理量を計算できる！

界面温度の計算結果
orderdisordersuper-heated ordered statefor

準安定状態の熱伝導による安定化という現象が生じる！
54

PART V

展望
55

数値実験

56

さて・・

熱力学の拡張大域熱力学 ！
自由エネルギー最小原理
Nakagawa and Sasa,
Global thermodynamics for heat
Conduction systems,
https://arxiv.org/abs/1906.07022

メソスケールダイナミクスの解析拡張された変分原理の導出「修正エントロピー」最大原理

Sasa, Nakagawa, Itami, Nakayama,
Stochastic order parameter dynamics
for phase coexistence in heat conduction,
https://arxiv.org/abs/1908.03029

57

今日の講演の共同研究者
大域熱力学、熱伝導下相共存全般、ゆらぎの解析:中川尚子（茨城大）

オンサーガ―の標準形、ノイズの特異性、ゆらぎの解析:
伊丹将人（京大）

メソ記述のモデル、ゆらぎの解析:
中山洋平（中央大）熱伝導下相共存の大規模数値実験：
小林未知数（京大）

58