Pressure-Induced Structural Phase Transitions in MBANP: A First-Principles Study of Elastic, Phonon, and Optical Properties

技術要約：MBANPにおける圧力誘起構造相転移

問題提起
有機非線形光学（NLO）結晶、特に2-(α-メチルベンジルアミノ)-5-ニトロピリジン（MBANP）は、その高い二次高調波発生（SHG）効率と構造的柔軟性により、次世代のフォトニクスおよびオプトエレクトロニクス技術において極めて重要である。MBANPは温度依存的な挙動や分子間相互作用に関しては広く研究されているが、静水圧に対する応答については全く未解明である。この応答を理解することは、非共有結合相互作用がいかにして結晶の安定性とNLO活性を支配しているかを解明する上で不可欠であり、また、応力を受ける小型化デバイスにおけるこれらの材料の機械的安定性を評価する上で実用的に不可欠である。

手法
本研究では、CASTEPコードを用いた密度汎関数理論（DFT）に基づく第一原理計算手法を用いる。計算には、弱い分子間相互作用を正確にモデル化するために、Perdew–Burke–Ernzerhof（PBE）汎関数にTkatchenko–Scheffler（TS）ファンデルワールス補正を加えた一般化勾配近似（GGA）を用いる。ノルム保存型擬ポテンシャルおよび遮断エネルギー830 eVの平面波基底関数を使用した。本研究では、0から20 GPaの範囲の静水圧下におけるMBANP単結晶を系統的に調査した。主な解析項目は以下の通りである：

• 構造進化： 格子定数および分子コンフォメーションのモニタリング。
• 弾性的性質： Voigt-Reuss-Hill平均法による弾性定数（$C_{ij}$）の算出および機械的弾性率（体積、剛性、ヤング、ポアソン比）の導出、ならびにBorn安定性基準によるチェック。
• 電子構造： バンド構造、状態密度（PDOS）、電荷密度、電子局在化関数（ELF）、および電荷移動メカニズムの解析。
• 分子間相互作用： Hirshfeld表面解析および2Dフィンガープリントプロットを用いた接触タイプの定量化。
• 振動および光学特性： ラマンスペクトルのシミュレーション、および周波数依存的な光学指標（誘電関数、屈折率、吸収、反射率）の算出。

主な貢献および結果

1. ~14 GPaにおける構造相転移：
   本研究は、14 GPa付近での圧力誘起等構造相転移を特定した。この閾値以下では、結晶は顕著な異方性圧縮性（$a$軸 > $c$軸 > $b$軸）を示す。14 GPaにおいて、$a$軸の異常な収縮に伴い、$b$軸および$c$軸の不連続な膨張と$\beta$角の跳躍が起こる。これは、協同的な分子再配列またはせん断様の歪みを意味しており、結果として準等方的な機械的特性を持つ高圧相をもたらす。これらの変化にもかかわらず、空間群の対称性（$P2_1$）は保持される。

2. 機械的安定性と硬化：
   MBANPは、修正Born安定性基準を満たし、フォノン分散スペクトルに虚数周波数が存在しないことから、0–20 GPaの全範囲にわたって機械的および動的に安定していることが確認された。材料は漸進的な硬化を起こし、体積弾性率は10.7 GPaから105.2 GPaへ、剛性率は3.76 GPaから27.2 GPaへと増加する。ポアソン比（$\nu$）は0–12 GPaの範囲で0.343から0.4へと単調に増加した後、0.38付近で安定し、延性的な性質を示している。$B/G$比は0–12 GPaの間で2.84から4.7へと大きく変化し、14 GPaで構造転移と相関する変曲点を経て、20 GPaでは3.86へと減少する。

3. 電子構造および電荷移動の増強：
   圧力は、ドナー・アクセプター軸（メチルベンジルアミノからニトロ基）に沿った分子内電荷移動（CT）を著しく増強させる。バンドギャップは、価電子帯（ピリジン $\pi$軌道）の広がりと伝導帯（ニトロ $\pi^*$ 状態）の広がりに起因して、0 GPaでの2.226 eVから20 GPaでの1.056 eVへと減少する。電荷密度解析により、ドナー窒素における電子減少とアクセプター酸素における電子蓄積の増加が明らかになり、N–H⋯O水素結合の長さの短縮とともに、「プッシュ・プル」効果の強化が確認された。

4. 分子間相互作用の再編成：
   Hirshfeld表面解析は、圧力駆動による相互作用の優位性のシフトを定量化している。弱く等方的なファンデルワールス接触（H⋯HおよびC⋯H）の寄与が減少する一方で、強く方向性のある相互作用が増加する。特に、$\pi$–$\pi$ スタッキングを示すC⋯C接触が急増し、特定のC–H⋯O接触は14 GPa以上で真正の水素結合（距離 < 2.6 Å）へと変化する。この再編成が、極限圧縮下での安定性を維持するための鍵となるメカニズムとして特定された。

5. 振動および光学応答：
   ラマンスペクトルは、結合の圧縮と力定数の増加に起因して、すべてのモードにおいて一様なブルーシフトを示す。14 GPaにおいて、ピーク強度と広がりにおける明確な不連続性が観察され、これが相転移を裏付けている。光学的には、静的誘電率と屈折率が圧力とともに増加する（それぞれ約3.56から4.4へ、1.86から2.1へ）。しかし、吸収端と主要な誘電関数ピークは高エネルギー側へシフト（ブルーシフト）しており、これは吸収端が高エネルギー側にシフトするという結論と一致している。その結果、可視領域における材料の吸収は弱まり、透過率の低下（100 nmの試料で87%から~73%へ）と反射率の増加を招く。

意義
本研究は、MBANPの高圧下における初の包括的な理論的研究を提供し、その構造・物性相関における重要な空白を埋めるものである。本研究は、静水圧が対称性を変えることなく、有機NLO結晶の機能的特性を調整するための強力なツールとして機能することを示している。機械的硬化、電荷移動増強、および分子間再編成のメカニズムを解明することで、本知見は小型化されたオプトエレクトロニクスデバイスへのMBANPの安定性を裏付けるとともに、調整可能な電気機械的および光学的特性を持つ次世代分子結晶の設計への指針を与えるものである。