Optimal preparation and reachable-state constraints in the Mpemba effect

技術的要約：ムペンバ効果における最適準備と到達可能状態の制約

問題提起
ムペンバ効果は、初期に高温の状態にある系が、低温の状態にある系よりも共通の最終状態に向かってより速く緩和するという現象であり、通常、非定常な初期条件を事前に（a priori）選択すると仮定した速度論的枠組みの中で分析される。しかし、このアプローチは、これらの状態を選択するために使用された特定の加齢プロトコル（aging protocol）によって課される物理的な制約をしばしば無視している。速度論的な定式化において、この効果は、平衡から逸脱した内部変数に符号化された、系の熱履歴の記憶に依存している。近年の研究では、準備時間（Pontus–Mpembaなど）を考慮したプロトコルが導入されているが、根本的な問いは依然として残っている。すなわち、「どの非平衡状態が、許容可能な駆動メカニズムを通じて物理的に到達可能なのか」であり、この「到達可能性の制約」がムペンバ効果の大きさをどのように制限するのかという点である。本論文は、確率論的なサーモスタットによって課される到達可能集合を明示的に考慮した上で、速度論的枠組みにおける、ムペンバ効果を最大化するための最適な加齢プロトコルを決定するという問題に取り組む。対象とするのは、一様に加熱された粒状流体である。

手法
著者らは、白色雑音による確率論的サーモスタットによって駆動される、非弾性衝突（反発係数 $\alpha$）を持つ硬球の $d$ 次元粒状ガスをモデル化している。系のメゾスコピックな状態は、粒状温度 $T$ と、速度分布関数（VDF）の非ガウス性を定量化する過剰尖度 $a_2$ によって記述される。解析は、第1ソネ近似（first Sonine approximation）内で行われ、ここでは $T$ と $a_2$ の進化は、結合された微分方程式の閉じた集合によって支配される。

準備段階は、以下の最適制御問題として定式化される：

1. 状態変数： 粒状温度 $T$ および過剰尖度 $a_2$。
2. 制御変数： 範囲 $[\chi_{\min}, \chi_{\max}]$ に制限された、確率論的サーモスタットの強度 $\chi(t)$。
3. 目的関数： 最終時刻 $t_f$ における過剰尖度 $a_2$ を極値化する制御プロトコル $\chi(t)$ を見つけること。固定された温度における冷却率が $a_2$ に直接比例するため、$a_2$ を極値化することは、その後の緩和率を最大化または最小化することに対応する。
4. 解法テクニック： 著者らは、最適性の必要条件を導出するために**ポントリャーギンの最大原理（PMP）**を適用する。彼らは、ハミルトニアン系、余剰コスト方程式、および横断条件を分析し、最適制御の構造を決定する。

主要な貢献と結果

• 最適プロトコルの構造： PMPの適用により、最適な準備プロトコルは常に**「ワンバング（one-bang）」プロトコル**であることが明らかになった。これは、制御 $\chi(t)$ が、全準備区間 $(0, t_f)$ を通じて、最大値 ($\chi_{\max}$) または最小値 ($\chi_{\min}$) のいずれかに一定であることを意味する。具体的な選択は、均質冷却状態（HCS）における過剰尖度の符号に依存し、その符号は反発係数 $\alpha$ によって決定される。

  • 「大きな」非弾性（$\alpha < \alpha_c \equiv 1/\sqrt{2}$）の場合、$a_2$ を最小化するための最適制御は $\chi_{\max}$ である。
  • 「小さな」非弾性（$\alpha > \alpha_c$）の場合、最適制御は $\chi_{\min}$ である。
  • $a_2$ を最大化する場合は、この論理が逆転する。

• 到達可能状態の制約： 中心的な発見は、非自明な到達可能状態の制約が存在することである。$(T, a_2)$ 平面におけるアクセス可能な非平衡状態の集合は、有界である。具体的には、アクセス可能な過剰尖度の値は、ガウス状態 ($a_2 = 0$) と、均質冷却状態における過剰尖度の値 ($a_2^{\text{HCS}}$) の間に限定される。

  • $\chi_{\min} = 0$（自由冷却）かつ $\chi_{\max} \to \infty$（瞬時かつ強力な加熱）という理想的な極限において、過剰尖度の極値は正確に $0$と$a_2^{\text{HCS}}$ となる。
  • これは、確率論的サーモスタットが、任意に大きな非ガウス性を生成することはできないことを意味する。すなわち、非ガウス性の大きさは、$a_2^{\text{HCS}}$ の値によって厳密に制限されている。

• 数値的検証： 理論的な予測は、直接シミュレーションモンテカルロ（DSMC）法を用いて検証された。シミュレーションは、制御の境界が緩和されるにつれて、過剰尖度の極値が理論的な境界（$0$および$a_2^{\text{HCS}}$）に収束することを確認した。最適制御理論と数値データの間の合致は、全反発係数の範囲において極めて良好であった。

• ムペンバ効果への示唆： 本研究は、速度論的なムペンバ効果の強さが、基礎となる非平衡状態へのアクセスのしやすさと密接に関連していることを示している。$a_2$ の到達可能集合は有界であるため、高温の試料と低温の試料の間で達成できる緩和率の非対称性には、根本的な限界が存在する。以前の重要な研究で仮説とされていた「大きな」非ガウス性は、確率論的サーモスタットのみでは生成できない。

意義と主張
本論文は、異常緩和、最適制御、および状態のアクセシビリティの間の厳密な関係を確立したと主張している。ムペンバ効果を、動的な到達可能性の制約下での最適状態準備の問題として再構成することで、著者らは、この効果を準備プロトコルの物理的な限界から独立して論じることはできないと主張している。

本研究は、過剰尖度が、コバックス効果における役割と同様に、系の記憶を符号化する重要な内部変数として機能することを強調している。結果は、最大に達しうるムペンバ効果は、材料固有の性質だけでなく、駆動メカニズムによって定義される到達可能な非平衡状態空間の幾何学的形状によって根本的に制約されることを示唆している。著者らは、確率論的サーモスタットがこれらの特定の境界を課しているものの、代替の駆動メカニズムは異なる到達可能集合をもたらす可能性があり、それがより強力な緩和異常を可能にする可能性があると述べているが、これは今後の調査課題である。結論として、緩和経路のエンジニアリングには、アクセス可能な状態空間の特定の幾何学をナビゲートすることが必要である。