有機半導体の性能を大きく左右する結晶構造を、分子の設計段階から予測する新たな計算手法が開発された。東京大学と山形大学の共同研究チームは、分子結晶の対称性と層状構造に着目し、分子の配置を決める変数を段階的に最適化することで、実験で知られる4種類の結晶構造を再現することに成功した。従来のように大量の候補構造を生成して安定なものを探すのではなく、分子の傾きや回転、層の重なり方を順に調べることで、それぞれの構造が安定化する理由も明らかにした。

研究を行ったのは、東京大学大学院工学系研究科の大野亮太大学院生、都築誠二特任研究員、長谷川達生教授、山形大学の井上悟研究専任准教授らのグループだ。成果は2026年10月23日付の科学誌『Communications Chemistry』に掲載された（DOI: 10.1038/s42004-026-02252-8）。対象となったのは、ベンゼン環が2〜6個連なった5種類のポリアセン分子であり、既知の結晶構造を計算の出発点として与えずに、その分子配列を再現した。新しい結晶を合成した研究ではなく、結晶構造の予測と、その形成を支配する分子間相互作用の理解を進めた成果である。

AD

膨大な候補の網羅探索から、分子配置の段階的な最適化へ

有機半導体では、同じ分子であっても、結晶中での並び方によって電気的な性質が大きく変わる。例えば、有機トランジスタに使われるペンタセンは、分子が形成する結晶構造によって電荷の移動しやすさが異なる。そのため、分子を合成する前に安定な結晶構造を予測できれば、材料開発に必要な試行錯誤を減らせる可能性がある。しかし、多数の分子が取り得る位置や向きの組み合わせは膨大であり、安定な構造を計算によって特定することは容易ではない。

従来の結晶構造予測（Crystal Structure Prediction：CSP）では、大量の候補構造を生成し、それぞれのエネルギーを計算して安定な構造を探す方法が広く使われてきた。この方法は有効だが、最終的な構造が得られても、なぜその分子配列が安定するのかを説明するには追加の解析が必要になる。研究チームは、この問題に対して、結晶が持つ幾何学的な規則性を利用し、探索する変数を体系的に整理する方法を採用した。

鍵となるのは「映進対称性」と「層状性」だ。映進対称性とは、ある平面に対する鏡映と、その平面に沿った平行移動を組み合わせると、分子の配置が重なる幾何学的な関係を指す。一方、層状性は、棒状の分子が長軸をおおむね同じ方向に揃えながら、横方向に並んで分子層を形成する性質である。有機半導体には、このような規則性を持つ結晶が多く存在する。

研究チームは、ベンゼン環が直線状に結合したポリアセン系列を検証対象に選んだ。具体的には、環が2個のナフタレン、3個のアントラセン、4個のテトラセン、5個のペンタセン、6個のヘキサセンの5種類だ。これらは骨格が硬い棒状分子であり、層状構造と対称性を利用した計算モデルの検証に適している。

開発した手法では、まず単純な分子層の基本構造を求め、次に分子の長軸を傾け、さらに分子面の回転や層間の相対位置を最適化する。各段階で分子間相互作用のエネルギーを計算し、安定な配置がどのように現れるかを追跡する仕組みだ。これにより、最終的な結晶構造を求めるだけでなく、構造の変化に伴うエネルギーの増減を「ポテンシャルマップ」として可視化できる。

分子が約25度傾くと安定する「魚の骨」構造

最初の計算では、分子の長軸が互いに平行で、長軸方向への相対的なずれを持たない二次元の分子層を想定した。映進対称性を持つ分子対を基本単位とし、それを周期的に並べることで層全体の構造を表す。この配置は、分子面と映進面のなす角度 と、分子間の相対位置を表す格子定数などによって記述できる。

研究チームは、分散力補正を取り入れた密度汎関数理論（DFT）を用い、分子間の相互作用エネルギーを計算した。分散力は、電子分布の瞬間的な揺らぎによって生じる引力であり、分子結晶の安定性を左右する重要な要因だ。角度 を変えながら分子層のエネルギーを調べたところ、約25度と65度の付近で、対称性によって等価な安定構造が得られた。

このとき、隣接する分子面がなすヘリンボーン角は約51度となる。これは、実験で観測されているペンタセン結晶の約53度に近い値だ。研究チームは、この基本となる安定な分子配列を「R構造（Regular）」と名付けた。

R構造では、隣接する分子が互いに角度を持ち、魚の骨のような模様を形成する。これは「層状ヘリンボーン（Layered Herringbone：LHB）構造」と呼ばれ、有機半導体で広く見られる代表的な分子配列だ。計算によって、この構造の安定性が、分散力や静電力による引力と、分子が近づきすぎた際に生じる交換反発のバランスによって決まることが確認された。

特に重要なのが、分子の端にある水素原子と、隣接するベンゼン環の中心付近との接触である。分子がT字型に近い配置を取ると、水素原子が隣の分子の環中心付近に近づき、分散力などによる安定化が強まる。一方、分子同士が過度に接近すると、電子雲の重なりに伴う交換反発が急激に増加する。約25度という角度は、こうした引力と斥力の釣り合いから現れた安定配置だと説明できる。

続いて研究チームは、R構造の分子層に対し、分子の長軸を一様に傾斜させる操作を導入した。傾斜方向を表す方位角 \phi_} と、傾斜角 \theta_} を変化させると、隣接分子が長軸方向に沿って相対的にずれ、ベンゼン環約1個分の変位に対応する位置で新たなエネルギーの極小点が現れた。

この操作によって得られたのが、映進対称性を保った「G構造（Glide）」と、映進対称性が破れた「N構造（Non-glide）」という2種類の準安定構造だ。つまり、基本となるR構造から分子を傾斜させるだけで、異なる結晶構造につながる2つの安定な配置が導かれたことになる。

AD

ベンゼン環の数によって変わる4種類の結晶構造

R、G、Nという3種類の分子層を得た研究チームは、さらに分子面に垂直な軸の周りの回転と、層同士の積み重なり方を計算に取り入れた。これらの配置を段階的に最適化した結果、ポリアセン系列で実験的に知られている4種類の結晶構造タイプ（I〜IV）を再現できた。

結晶構造 主な幾何学的特徴 代表的な実験例 主な安定化要因
多型I R構造を基本とする、分子長軸がほぼ平行な層状構造 テトラセン、ペンタセンの薄膜相 層内のT字型接触と分散力・静電力・交換反発の均衡
多型II 映進対称性が破れたN構造を基礎とする積層 ペンタセン単結晶、ヘキサセン 分子の長軸方向のずれと、層間の密な接触
多型III G構造から分子面を回転させ、映進対称性を維持 ナフタレン、アントラセン 分子の回転による層間の空隙充填と新たな分子間接触
多型IV N構造の異なる層間配置と、分子長軸のわずかな不均一傾斜 テトラセンのバルク相、ペンタセンの別の結晶相 層同士の噛み合わせと層間相互作用の増大

環数の少ないナフタレンとアントラセンでは、G構造から分子を回転させると、実験構造とよく一致する多型IIIが得られた。分子面を回転させることで層内に隙間が生まれ、そこへ隣接する層の分子が入り込む。さらに、水素原子とベンゼン環の新たな接触も形成され、層間相互作用が強まる。この安定化は分子長の短いポリアセンで特に有効であり、長い分子では同じ回転によるエネルギー上の利点が小さくなる。

一方、ペンタセンやヘキサセンのような長い分子では、N構造を基礎とする積層が有利になる。N構造の分子層を水平方向にずらすと、層表面の凹凸が互いに噛み合う配置で格子体積が小さくなり、2つの極小点が現れる。これらの配置は、実験で知られる多型IIと多型IVに対応した。

さらに研究チームは、分子長軸のわずかな不均一傾斜も考慮し、両構造の安定性を比較した。その結果、環数が4個のテトラセンでは多型IVが有利になり、環数が5個のペンタセンと6個のヘキサセンでは多型IIが優勢になることが分かった。分子が長くなるほど、長軸方向にずれた配置で分子同士の接触を維持しやすくなり、N構造を基礎とする密な積層が安定するためと考えられる。

こうした計算結果から、ポリアセン系列では、環数の増加に伴って最も安定な結晶構造が多型IIIからIV、IIへと変化する傾向が説明できた。これは実験で知られる構造の違いともよく一致している。

また、研究チームは4種類の多型について、層内でT字型に接する分子間相互作用と、平行に並ぶ分子間相互作用を比較した。その結果、両者の強度比は、いずれの多型でもおおむね2対1に保たれていた。分子の傾きや層の重なり方が異なっても、ヘリンボーン構造を支える基本的な分子間相互作用のバランスは共通していることを示す結果だ。

結晶構造の再現に成功、それでも残る基板と熱の影響

今回の手法は、既知の結晶構造を計算の初期配置として使わず、5種類のポリアセンに見られる結晶構造タイプを再現した。ただし、これを未知の化合物の結晶構造を完全に予測できるようになった成果と捉えることはできない。対象としたのは、すでに実験によって構造が明らかになっている分子であり、層状性や高い対称性を仮定したモデルの有効性を検証した研究だからだ。

計算結果には、相互作用をどこまで考慮するかによる違いも見られた。近接する分子間の相互作用だけでは、多型間のわずかな安定性の差を十分に説明できない場合があったが、離れた分子からの寄与も取り入れることで、実験で知られる安定性の傾向に近づいた。局所的な分子配列は主に近くの分子との相互作用によって決まる一方、結晶全体としてどの多型が安定するかには、より広い範囲に及ぶ相互作用も重要となる。

また、ペンタセン薄膜に見られる多型Iでは、計算上の理想的なR構造と実験構造の間にわずかな違いが確認された。実際の薄膜では、基板との相互作用や界面の影響によって分子配列が変化する可能性があるためだ。研究チームは、計算結果と実験構造の差にこうした基板・界面効果が関わっていると考えている。分子間相互作用だけを扱うモデルでは、薄膜成長時の条件まで完全に再現できるわけではない。

熱力学的な制約も残る。今回の計算は静的な分子配置の相互作用エネルギーを主に評価したものであり、有限温度での相安定性を厳密に判断するには、エントロピーを含むギブズ自由エネルギー $G=H-TS$ を考慮する必要がある。研究チームは、計算された安定性と実験結果がよく一致することから、今回のポリアセン系列ではエントロピーの影響が比較的小さい可能性を指摘している。しかし、柔軟な分子や複雑な結晶にまで、この前提をそのまま適用できるとは限らない。

さらに、今回の研究では、結晶構造を求めるための幾何学的な変数を段階的に導入しているが、一部の変数の選択には既知の実験構造との比較から得られた知見も利用している。そのため、実験構造を初期値として使用していないことと、実験情報を一切利用せずに未知の構造を予測したことは区別する必要がある。未知の分子に対しても同じ手順で正しい結晶構造を導けるかどうかは、今後の検証課題となる。

AD

分子のわずかな位置の違いが電荷輸送特性を大きく変える

研究チームは、再現した結晶構造について、半導体の電荷輸送に関係する「分子間移動積分（Transfer Integral）」も計算した。移動積分は、隣接する分子の電子軌道がどの程度結び付いているかを表す指標であり、電荷が分子間を移動する際の電子的な結合の強さに関係する。

計算の結果、4種類の結晶構造では、分子間相互作用の基本的なバランスがほぼ共通しているにもかかわらず、移動積分の大きさには顕著な違いがあった。特にペンタセンでは、多型Iの移動積分が大きく、多型IIやIVではそれより小さくなる傾向が確認された。

この違いを生む主な要因は、隣接分子の長軸方向への相対的なずれだ。分子がずれると、電荷輸送に関わる分子軌道の重なり方が変化する。今回の計算では、そのずれがベンゼン環の約半分に相当する位置で移動積分が小さくなり、約1個分に達すると再び大きくなる挙動が見られた。分子軌道の波動関数には正負の位相を持つ部分があり、分子の位置関係によって軌道間の結合が強まったり、互いに打ち消し合ったりするためだ。

これは、結晶構造の微小な違いが、有機半導体の電子物性に大きく影響することを示している。ただし、移動積分が大きいからといって、実際のトランジスタの電荷移動度が必ず高くなるわけではない。デバイスの性能は、結晶粒界や欠陥、薄膜の形態、電極との接触抵抗、格子振動による散乱などにも左右される。今回明らかになったのは、結晶構造と電荷輸送に関係する電子的な結合の理論的な関係であり、実際のデバイス性能を測定した結果ではない。

柔軟な側鎖を持つ有機半導体への拡張が次の課題

東京大学の研究発表によると、研究チームは今後、アルキル鎖などの柔軟な側鎖を持つ有機半導体材料へ手法を拡張する方針だ。実用的な有機半導体には、溶液から薄膜を形成しやすくするため、分子骨格に長い側鎖を導入したものが多い。こうした材料では、硬い主骨格の配置だけでなく、側鎖の形状や運動も結晶構造に影響するため、探索すべき自由度が増える。

今回の研究で有効性が確認されたのは、対称性が高く、剛直な棒状分子からなるポリアセン系列である。柔軟な側鎖を持つ分子や、より複雑な結晶に対しても、同じように配置変数を整理し、安定な構造を導けるかはまだ分かっていない。さらに、基板との相互作用や温度による自由エネルギーの変化も取り込めれば、実際の製造条件に近い環境で結晶構造を予測するための重要な前進となる。

本研究の意義は、結晶構造を計算で再現したことだけでなく、分子の形状や配置の変化が、どのような相互作用を通じて結晶を安定化させるのかを体系的に説明した点にある。従来の網羅的な結晶構造予測を置き換えるというよりも、そこに物理的な理解を加える新たな手法と位置付けられる。未知の分子に対する予測精度が検証されれば、望ましい結晶構造や電荷輸送特性を持つ有機半導体を、分子設計の段階から合理的に開発するための有力な手段となる可能性がある。