XCP entdeckt 44.030 mögliche 2D-Materialien in nicht geschichteten Kristallen

Am 29. Oktober 2026 richtet sich der Blick an der TU Dresden auf Kristalle, die mehr können sollen, als ihre dreidimensionale Form zunächst vermuten lässt. Ein Team des Helmholtz-Zentrums Dresden-Rossendorf (HZDR) und der TU Dresden hat ein Rechenverfahren entwickelt, das spaltbare Kristallflächen und mögliche zweidimensionale Materialien systematisch aufspürt.

Seit der Herstellung von Graphen im Jahr 2004 ist das Interesse an 2D-Materialien groß. Sie können elektronische, magnetische und chemische Eigenschaften entwickeln, die in ihren dreidimensionalen Ausgangsmaterialien nicht auftreten. Bisher suchten Forschende vor allem in Van-der-Waals-Kristallen: Deren schwach gebundene Schichten lassen sich vergleichsweise leicht ablösen. Bei nicht geschichteten Kristallen sitzen die Atome dagegen in einem dreidimensionalen Netzwerk, verbunden durch kovalente, ionische oder metallische Bindungen. Genau dort wurde die Suche bisher knifflig.

Ein Rechenmodell für schwierige Kristalle

Das neue Verfahren heißt „eXfoliation and Cleavage Potential“, kurz XCP. Es schätzt anhand von Atomsorten und interatomaren Abständen die Stärke atomarer Bindungen. Dafür kombiniert das Modell Morse- und Yukawa-Potentiale: Der Morse-Anteil beschreibt vor allem kovalente und metallische Bindungen, der Yukawa-Anteil ionische Wechselwirkungen.

Die fünf Potentialparameter wurden für 1540 Atompaar-Kombinationen mithilfe von DFT-Daten aus der AFLOW-ICSD-Datenbank angepasst. Zusammen umfasst XCP 7700 Parameter. Weil die Referenzdaten auf PBE(+U)-DFT beruhen, übernimmt das Modell auch deren systematische Eigenschaften, etwa die tendenzielle Unterschätzung dispersiver Wechselwirkungen.

Zwei Algorithmen gehen anschließend auf unterschiedliche Weise vor. BONDDEL entfernt schrittweise die schwächsten Bindungen, bis verbundene zweidimensionale Teilstrukturen übrig bleiben. Dabei erkennt das Verfahren auch atomar gekrümmte oder nicht-planare Bruchflächen. HKLSEARCH bewertet dagegen die Bindungsenergiedichte verschiedener Miller-Ebenen und sucht nach Flächen mit niedriger Oberflächenenergie. So lassen sich auch atomar unebene Trennflächen erfassen.

Ein Exfoliierbarkeitsverhältnis R vergleicht die Energie, die nötig ist, um ein 2D-Blatt aus dem 3D-Kristall herauszulösen, mit der Energie, die für dessen weiteres Zerlegen benötigt würde. Hohe Werte sprechen für eine leichte Abtrennung bei zugleich stabiler zweidimensionaler Struktur. Der Rechenaufwand bleibt überschaubar: Große Materialdatenbanken können zunächst mit XCP durchsucht werden; aufwendige elektronische Strukturrechnungen folgen erst für aussichtsreiche Kandidaten.

24.566 Kristalle, 44.030 mögliche Flächen

Für die Studie, veröffentlicht am 25. August 2026 in Nature Communications, untersuchten Tom Barnowsky, Carsten Timm und Rico Friedrich experimentell bekannte, nicht geschichtete Bulk-Strukturen aus der AFLOW-ICSD-Datenbank. Aus 24.566 dreidimensionalen Kristallen ermittelte XCP 44.030 spaltbare Flächen beziehungsweise mögliche 2D-Strukturen.

Nach dem Ausschluss von Systemen, deren chemische Zusammensetzung sich beim Ablösen verändert, blieben 17.897 Kandidaten. Für 2531 compositionserhaltende Systeme mit R > 1,5 folgten detaillierte DFT-Rechnungen. Diese Auswahl ist bewusst konservativ: Auch Systeme mit niedrigeren R-Werten können grundsätzlich exfolierbar sein.

Bei den Kandidaten bleiben chemische Zusammensetzung und berechnete Bindungsverhältnisse für die Bildung zweidimensionaler Strukturen günstig. Für den Großteil wurden vergleichsweise geringe Ablösungsenergien berechnet. Sie reichen von nahezu null bis etwa 0,5 eV/Ų; 85,9 Prozent der relaxierten Kandidaten liegen unter 140 meV/Ų, dem für Hematene angegebenen Vergleichswert.

Die chemische Landschaft ist breit. Rund 30 Prozent der Oberflächen sind durch Sauerstoff terminiert. Häufig folgen H⁺ mit 4,6 Prozent, Cs⁺ mit 4,3 Prozent, F⁻ mit 4,2 Prozent und K⁺ mit 4,1 Prozent. Bei der tatsächlichen Herstellung können Lösungsmittel oder Passivierung diese Oberflächenchemie allerdings verändern.

Am häufigsten kommen bei den dreidimensionalen Elternstrukturen die Raumgruppen P̅1 (10,9 Prozent), P2₁/c (10,5 Prozent), I4/mmm (6,2 Prozent), Pnma (5,9 Prozent) und Cmcm (5,6 Prozent) vor. Die gefundenen Schichten decken alle fünf möglichen 2D-Bravais-Gitter ab: rechteckig (43,3 Prozent), schiefwinklig (18,7 Prozent), quadratisch (17,4 Prozent), hexagonal (13,4 Prozent) und rechteckig-zentriert (7,2 Prozent). Genannt werden außerdem 504 rechteckige, 295 schiefwinklige, 121 quadratische, 112 hexagonale und 73 rechteckig-zentrierte oxidische Schichten.

Stabil auf dem Papier – und nach der Ablösung?

Die Rechnung macht nicht bei der möglichen Spaltfläche halt. Nach dem Kriterium ΔE_hull < 100 meV/Atom gelten 49,0 Prozent der Monolagen als thermodynamisch stabil. Bei zweilagigen Systemen steigt der Anteil auf 70,0 Prozent. Strukturrelaxationen können die Exfoliationsenergien deutlich senken, in einzelnen Fällen bis um den Faktor vier.

Auch die Dynamik wurde geprüft. 64,4 Prozent der untersuchten Kandidaten gelten anhand von Phononrechnungen als dynamisch stabil; als Instabilitätsschwelle gilt eine maximale imaginäre Frequenz von mehr als 0,1 THz. Eine Phononstatistik ließ sich für 65,3 Prozent der Kandidaten erstellen, weil die verwendeten MLIP-relaxierten Strukturen ausreichend mit den DFT-Strukturen übereinstimmten.

Beim Heraustrennen kann sich eine Schicht kräftig verformen. Fläche und Dicke ändern sich häufig in entgegengesetzte Richtungen. Beim (010)-MgSiO₃ wächst die Fläche um 15 Prozent, während die Dicke um 30 Prozent sinkt. Bei (111)-K₂O schrumpft die Fläche um 18 Prozent, die Dicke nimmt um 23 Prozent zu. (001)-SiO₂ legt durch eine weitreichende Rekonstruktion sogar 40 Prozent Fläche zu. Und bei (001)-OsP₄ gehen 18 Prozent Fläche verloren, die Dicke steigt um 23 Prozent; zugleich werden die Phosphoratome nahezu vollständig neu geordnet.

Die Methode wurde an bereits experimentell hergestellten Materialien getestet. Bei nicht-vdW-2D-Materialien erkannte sie in etwa 80 Prozent der Fälle die richtigen Exfoliationsebenen – 23 von 29 Strukturen. Ein mehrdeutiger Pyrit-Fall wurde wegen widersprüchlicher experimenteller Angaben nicht in die Statistik aufgenommen. Traditionelle Van-der-Waals-Materialien erkennt XCP wegen ihrer schwachen Zwischenlagenbindungen ebenfalls zuverlässig. Bei Silizium fand das Verfahren die experimentell bevorzugten {111}-Spaltebenen.

Vom Halbleiter zum Metall

Nicht-vdW-2D-Materialien entstehen aus nicht geschichteten Kristallen, wobei starke kovalente oder ionische Bindungen gebrochen werden. Bekannte Beispiele sind MXene, 2D-WO₃, Hematene – also 2D-Fe₂O₃ – und Ilmenene, 2D-FeTiO₃. Frühere Rechenverfahren konzentrierten sich überwiegend auf Van-der-Waals-geschichtete Materialien und konnten dreidimensionale Bindungsnetzwerke nicht zuverlässig bewerten.

Gerade beim Übergang in die Fläche kann sich die elektronische Struktur drehen, als hätte jemand im Material einen Schalter umgelegt. B₆O wird vom direkten Bulk-Halbleiter zum 2D-Metall. Ba₅Ta₄O₁₅, Ca₂Al₂SiO₇ und Cs₂SiP₂ entwickeln einen direkten Bandabstand, obwohl ihre dreidimensionalen Phasen indirekte Bandlücken besitzen – interessant für optoelektronische Bauelemente.

Sr₃Fe₂O₇ bleibt als quadratische Monolage ferromagnetisch und halbmetallisch. Rb₂TlAu₃ behält als rechteckige Schicht seinen metallischen Charakter. Viele Kandidaten enthalten magnetische Übergangsmetalle und