Atomistic study of finite temperature properties in ferroelectric BiAlO$_3$

Technische Zusammenfassung: Atomistische Untersuchung der temperaturabhängigen Eigenschaften von ferroelektrischem BiAlO3

Problemstellung
Bleifreie Perowskit-Ferroelektrika sind entscheidend für nachhaltige technologische Anwendungen, doch die Suche nach praktikablen Alternativen zu bleihaltigen Materialien bleibt eine Herausforderung. Bismut-Aluminat (BiAlO3) hat sich als ein vielversprechender Kandidat erwiesen, der mittels Hochdruck-Hochtemperatur-Techniken synthetisiert wird. Während experimentelle Studien seinen rhomboedrischen ($R3c$) Grundzustand und seine ferroelektrische Natur etabliert haben, gibt es eine signifikante Diskrepanz zwischen experimentellen Daten und theoretischen Vorhersagen. Insbesondere berichten Dichtefunktionaltheorie-Berechnungen (DFT), die auf Null Kelvin beschränkt sind, über spontane Polarisationswerte ($72,3–82,5 \, \mu\text{C/cm}^2$), die erheblich von experimentellen Messungen ($9–29 \, \mu\text{C/cm}^2$) abweichen. Darüber hinaus fehlt es bestehenden computergestützten Ansätzen an der Fähigkeit, das Verhalten von BiAlO3 bei endlichen Temperaturen zu modellieren, was Lücken im Verständnis der Temperaturabhängigkeit der spontanen Polarisation, der Dehnung, der Koerzitivfelder und der Kopplung zwischen ferroelektrischen (FE) und antiferrodistortiven (AFD) Moden hinterlässt.

Methodik
Um diese Einschränkungen zu adressieren, entwickelten die Autoren ein auf der ersten Prinzipien basierendes atomistisches Modell, um BiAlO3 bei endlichen Temperaturen zu simulieren. Die Methodik umfasste die folgenden Schritte:

1. Phononenanalyse: Mittels der Dichtefunktionalen Störungstheorie (DFPT) innerhalb des VASP-Pakets identifizierten die Autoren strukturelle Instabilitäten in der kubischen ($Pm\bar{3}m$) Phase. Dies offenbarte imaginäre Frequenzen, die dem $\Gamma$-Modus (ferroelektrische Verzerrung) sowie den $R$- und $M$-Punkten (antiferrodistortive Sauerstoffoktaeder-Kippungen) entsprechen.
2. Konstruktion eines effektiven Hamiltonians: Basierend auf den identifizierten Instabilitäten wurde ein effektiver Hamiltonian parametrisiert. Dieser Hamiltonian umfasst Terme für das ferroelektrische lokale Modus ($u$), das antiferrodistortive Modus ($\omega$, entsprechend den $R_{25}$-Verzerrungen) und die Dehnung ($\eta$). Entscheidend ist, dass er die Kopplung zwischen diesen Moden berücksichtigt (FE-Dehnung, AFD-Dehnung und FE-AFD). Die Parameter wurden aus DFT-Berechnungen abgeleitet und sind in atomaren Einheiten angegeben.
3. Monte-Carlo-Simulationen: Der parametrisierte Hamiltonian wurde in Single-Flip-Metropolis-Monte-Carlo-Simulationen (MC) verwendet. Eine Superzelle mit $20 \times 20 \times 20$ Einheitszellen unter Verwendung periodischer Randbedingungen wurde eingesetzt. Die Simulationen nutzten eine Simulated-Annealing-Technik, wobei von 1500 K auf 10 K abgekühlt wurde. Um die konkurrierenden Instabilitäten nahe dem Phasenübergang genau zu erfassen, wurde die Anzahl der MC-Schritte in der Nähe der Curie-Temperatur ($T_C$) erhöht.
4. Externe Perturbationen: Das Modell wurde erweitert, um die Auswirkungen externer elektrischer Felder (entlang [111], [110] und [100]), hydrostatischen Drucks sowie uniaxialer/biaxialer Spannungen zu untersuchen.

Wesentliche Ergebnisse

• Grundzustand und Phasenübergang: Die Simulationen bestätigen einen strukturellen Phasenübergang von der paraelektrischen kubischen Phase zur ferroelektrischen rhomboedrischen ($R3c$) Phase. Die vorhergesagte Curie-Temperatur ist $1160 \, \text{K}$, und die spontane Polarisation bei niedrigen Temperaturen beträgt $81 \, \mu\text{C/cm}^2$ entlang der [111]-Richtung. Die Autoren merken an, dass $T_C$ höher ist als der experimentelle Wert ($>793 \, \text{K}$), was wahrscheinlich auf eine Überschätzung der Bindungsenergie in der Parametrisierung zurückzuführen ist, während die Polarisation gut mit den DFT-Grundzustandswerten übereinstimmt.
• Modenkopplung: Es wird eine starke Kopplung zwischen ferroelektrischen und antiferrodistortiven Moden über ein breites Spektrum von Temperaturen und elektrischen Feldern beobachtet. Der AFD-Ordnungsparameter folgt dem Trend des FE-Modus, wobei beide entlang der [111]-pseudokubischen Richtung orientiert sind.
• Reaktion auf elektrische Felder: Unter einem elektrischen Feld entlang [111] zeigt das Material eine einzelne Hystereseschleife, was die $R3c$-Phase bestätigt. Die Sättigungspolarisation ($P_{Sat}$) variiert leicht mit der Feldrichtung ($57,8–60,2 \, \mu\text{C/cm}^2$), während das Koerzitivfeld ($E_C$) eine signifikante Anisotropie aufweist ($2,4–3,8 \, \text{MV/cm}$). Die AFD- und Dehnungsordnungsparameter zeigen "Schmetterlings"-Hysteresekurven, wobei die Dehnung eine hohe Sättigung ($0,021$) und eine signifikante Remanenzdehnung ($0,018–0,019$) aufweist.
• Temperaturabhängigkeit: Mit steigender Temperatur nehmen die Fläche der Hystereseschleife, die Sättigungspolarisation, die Remanenzpolarisation und das Koerzitivfeld ab. Dieser Trend ist konsistent mit anderen ferroelektrischen Materialien.
• Druck- und Spanneffekte:
  • Hydrostatischer Druck: Positiver hydrostatischer Druck unterdrückt sowohl die spontane Polarisation als auch die Curie-Temperatur (Rate von $-45 \, \text{K/GPa}$), während negativer Druck den gegenteiligen Effekt hat.
  • Uniaxiale Spannung: Uniaxiale Zugspannung induziert eine Sequenz von zwei Phasenübergängen (kubisch $\to$ tetragonal $\to$ rhomboedrisch), während Druckspannung zu einem einzigen Übergang führt (kubisch $\to$ rhomboedrisch).
  • Biaxiale Spannung: Der Trend ist umgekehrt; biaxiale Zugspannung führt zu einem einzigen Übergang, während biaxiale Druckspannung zwei Übergänge induziert.

Bedeutung und Ansprüche
Das Paper beansprucht, eine kritische Lücke in der Literatur zu schließen, indem es den ersten computergestützten Ansatz liefert, um DFT-Simulationen auf endliche Temperaturen für BiAlO3 zu erweitern. Durch die Entwicklung eines parametrisierten effektiven Hamiltonians bietet die Studie atomistische Einblicke in das Verhalten des Materials unter experimentell relevanten Bedingungen, einschließlich der Rolle der FE-AFD-Kopplung und der Reaktion auf mechanische Lasten.

Die Autoren heben hervor, dass die durch ihre Simulationen vorhergesagten hohen Sättigungs- und Remanenzdehnungen darauf hindeuten, dass BiAlO3 ein vielversprechender Kandidat für Formgedächtnis-Anwendungen sein könnte. Zudem zeigt die Studie, dass die Curie-Temperatur und die Phasenstabilität von BiAlO3 durch hydrostatischen Druck und mechanische Spannung gesteuert werden können, was potenzielle Wege zur Optimierung des Materials für spezifische Geräteanwendungen eröffnet. Die Arbeit unterstreicht die Notwendigkeit, über Null-Kelvin-DFT-Berechnungen hinauszugehen, um die funktionellen Eigenschaften bleifreier Ferroelektrika vollständig zu verstehen.