Computational Determination of Optimal Growth Protocols for Metastable Polymorphs

技術要約：メタステーブル（準安定）ポリモーフの最適成長プロトコルの計算論的決定

問題提起
有機–無機界面における特定のターゲット構造の確実な成長は、特に目的の構造が熱力学的にメタステーブル（準安定）である場合、中心的な課題として残っている。熱力学は系を最も低いギブス自由エネルギー状態へと駆動するが、速度論的効果によって系がメタステーブルな構成にトラップされることがよくある。このようなメタステーブルな構造は、しばしば優れた特性（例：高いバイオアベイラビリティや電荷キャリア移動度）を示すが、それらを選択的に成長させるための実験的プロトコルを特定することは困難である。温度を低減させる標準的な戦略は、競合する安定相の形成が同時に起こるため、許容できないほど長い実験時間や低い収率を招くことになる。さらに、複数の競合するメタステーブル構造が存在する場合があり、それがターゲット物質の収率をさらに低下させる。実験的に最適な時間依存の温度および圧力プロトコルを特定することは、時間とコストがかかる。

手法
著者らは、キネティック・モンテカルロ（kMC）シミュレーションと、パラメータ化された核生成・成長モデルを組み合わせ、最適制御理論（OCT）によって最適化された計算フレームワークを提案している。

1. モデル系： 本研究では、金属表面上の$\pi$共役有機分子の典型的なモデルを用い、これを4つの潜在的なポリモーフに抽象化している。
   • LHB： ヘリンボーン構造における分子の平伏状態（寄与は無視できる）。
   • LBW： ブリックウォール構造における分子の平伏状態（高温・低圧において熱力学的に安定）。
   • SHB： ヘリンボーン構造における分子の直立状態（低温・高圧において安定）。
   • SBW： ブリックウォール構造における分子の直立状態（メタステーブルであり、ターゲットとなる構造）。
2. 速度論的モデリング： 相転移速度を記述するために、パラメータ化された有効核生成・成長モデルを定式化している。表面占有率$c_i$の変化率は、温度（$T$）および分圧（$p$）に依存する核生成速度（$n_{ji}$）および成長速度（$k_{ji}$）によって支配される。これらの速度は、自由エネルギー障壁と試行頻度を含むアレニウス型の式に従う。
3. パラメータ化： 素過程の速度（拡散、再配向、吸着）は、$(T, p)$の細かいグリッド上で実行されたkMCシミュレーションから導出される。有効モデルのパラメータは、kMCデータに対する最小二乗法によるフィッティングを通じて得られる。
4. 最適制御の定式化： 成長プロトコルは最適制御問題として扱われる。目的は、固定された期間（$t_f = 1800$ s）内におけるメタステーブル相SBWの最終収率（$c_{SBW}(t_f)$）を最大化することである。制御変数は、時間依存の温度$T(t)$および対数圧力$q(t) = \log_{10} p(t)$である。
   • 制約条件： 現実的な実験的境界が課されており、温度変化は$1 \text{ K s}^{-1}$、圧力変化は100秒あたり1桁の範囲に制限されている。温度範囲は$300\text{--}525$ K、圧力は$10^{-8}\text{--}10^{-6}$ barに制限されている。開始点と終了点において時間微分がゼロになることを要求することで、滑らかな遷移が強制される。
   • ソルバー： 問題は、Pythonフレームワークyapssを用いた直接擬スペクトル法を用いて解かれ、Ipoptソルバーが使用される。

主な結果

• ベースライン性能： 手動で最適化された3ステップのプロトコル（「成長点」への急速な遷移、一定条件での保持、急速な復帰）がベースラインとして確立された。このアプローチでは、SBW相の最大収率は**73%**であり、残りの27%は望ましくない安定相であるSHBであった。
• 最適化されたプロトコル： OCTから導出されたプロトコルは、同じ総持続時間において、収率を**97%**へと大幅に向上させた。最適化された軌跡は、手動のアプローチとは根本的に異なっている。
  • 初期段階： SBWに有利な低温・高圧領域へ即座に移動するのではなく、プロトコルは最初に圧力を下げる。これにより、系をLBW安定領域に留めつつ、低圧下では速度論的に抑制されるSHBの蓄積を最小限に抑える。
  • 成長段階： 温度が十分に下がり、LBWからSBWへの転換速度がLBWからSHBへの転換速度を上回る（$k_{LBW\to SBW} > k_{LBW\to SHB}$）領域に入ると、プロトコルは圧力を最大許容値（$10^{-6}$ bar）まで急速に上昇させ、転換を駆動する。
  • 最終段階： プロトコルが終了する前に完全な転換を確実にするため、終盤に向けて温度を徐々に上昇させる。その後、最終的な実験条件へと急速に戻る。
• 検証： 最適化されたプロトコルはkMCシミュレーションによって検証された。有効モデルは97%の収率を予測したが、kMCシミュレーションでは最終段階においてわずかな偏差が見られ、最終的な収率はわずかに低くなった。著者らは、この原因を、有効モデルがトレーニングデータには十分に表現されていない領域（高温、SBWの被覆率がほぼ最大の状態）へと外挿されているためであり、モデルがSHBの核生成を過小評価しているためであるとしている。それでもなお、kMCの結果は、最適化された軌跡が高収率の成長を実現する条件を特定することに成功したことを裏付けている。

意義および主張
本論文は、時間依存の温度および圧力プロトコルを計算論的に最適化することで、特定のメタステーブルな表面ポリモーフを高収率で選択的に成長させられることを示している。成長プロセスを最適制御問題として定式化することにより、高収率の成長は、単に一つの有利な成長点へと可能な限り速く移動することによって達成されるのではなく、特定の速度論的領域を利用し、熱力学的に安定な構造が優先される領域を回避するように系をフェーズ空間内で誘導することによって達成されることを著者らは示している。

著者らは、このアプローチがメタステーブルな構造形成を制御するための合理的な手法を提供すると主張している。彼らは、マイクロスコピックな構造がデバイス性能を決定するという、有機エレクトロニクスにおける本研究の関連性を強調している。望ましい界面ポリモーフを選択的に成長させる能力は、高性能な有機電子デバイスの製造における現在の限界を克服するための経路を提供する。本研究は、この手法が物理的に解釈可能な速度論的モデルと最適制御理論の組み合わせに依存していることを強調しており、試行錯誤的な実験的最適化に代わる、より効率的な選択肢を提供している。