🔬 materials science

Relaxation approach to quantum-mechanical modeling of ferroelectric and antiferroelectric phase transitions

This paper introduces a novel quantum-mechanical relaxation framework that replaces traditional classical Arrhenius assumptions to enable efficient, first-principles modeling of ferroelectric and antiferroelectric phase transitions and their characteristic hysteresis loops.

Original authors: Nikhilesh Maity, Sergey Lisenkov, Arlies Valdespino, Milan Haddad, Lewys Jones, Amit Kumar, Nazanin Bassiri-Gharb, Inna Ponomareva

Published 2026-06-30

📖 4 min read☕ Coffee break read

CC BY 4.0

Original authors: Nikhilesh Maity, Sergey Lisenkov, Arlies Valdespino, Milan Haddad, Lewys Jones, Amit Kumar, Nazanin Bassiri-Gharb, Inna Ponomareva

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine a material that acts like a tiny, internal switch. When you apply an electric field, the atoms inside shift their positions, creating a new state. This is the world of ferroelectrics and antiferroelectrics. Think of them as the electrical cousins of magnets: just as magnets can be "on" or "off" (or point in different directions), these materials can be polarized in different ways.

For a long time, scientists tried to predict how these materials switch states using classical physics. They imagined the atoms as little balls sitting in a valley (a stable state) and needing a big push (heat or an electric field) to roll over a hill (an energy barrier) into a new valley. This is called the "Arrhenius" model.

The Problem with the Old Way
The paper explains that this classical "ball rolling over a hill" idea has a major flaw. When scientists used it to calculate how much electric field is needed to flip these materials, the numbers were wildly wrong. The model predicted you would need an electric field thousands of times stronger than what is actually observed in real life. It was like predicting you need a rocket ship to push a shopping cart, when in reality, a gentle nudge is enough.

The old model assumed the atoms were heavy, classical objects that had to wait for a lucky thermal bump to jump the barrier. But the paper argues this is wrong.

The New "Quantum Relaxation" Approach
The authors propose a completely different way to look at the problem. Instead of a ball waiting to jump a hill, they suggest the atoms behave like quantum waves that naturally "relax" or settle down into their new position.

Here is the analogy they use:

The Old View: Imagine a marble stuck in a small dip on a hill. It can't get out unless a random gust of wind (heat) blows hard enough to push it over the edge. If the hill is high, it might wait forever.
The New View: Imagine the marble is actually a ghostly wave. It doesn't need to "jump" over the wall. Instead, it slowly leaks through the wall or naturally settles into the lower valley because it is interacting with its surroundings (the environment). It's a process of relaxation, not a sudden jump.

How They Tested It
The researchers built a new computer model based on this "quantum relaxation" idea. They treated the atoms as quantum waves that evolve over time, slowly losing energy and settling into the most stable state.

They tested this on two famous materials:

PbTiO3 (a ferroelectric): The old model predicted it would need a massive electric field to switch. The new model predicted a field that matched real-world experiments perfectly.
PbZrO3 (an antiferroelectric): The old model predicted a weird, mixed behavior that never happens in reality. The new model correctly predicted the smooth, double-loop shape seen in experiments.

Why This Matters
The paper claims that by treating these phase transitions as quantum mechanical relaxation processes rather than classical jumps, they can finally simulate these materials accurately using "first-principles" (calculations based purely on the laws of physics, without guessing).

Key Takeaways from the Paper:

Classical physics fails here: The old "ball on a hill" model gives unrealistic results for these materials.
Quantum mechanics is essential: Even at room temperature, these materials behave like quantum systems, not classical ones.
Relaxation is the key: The switching happens because the system naturally relaxes toward equilibrium, not because it waits for a thermal jump.
Broad potential: The authors suggest this "relaxation" framework could eventually help understand other types of transitions (like magnetic or chemical changes), but for now, they have proven it works specifically for ferroelectric and antiferroelectric materials.

In short, the paper says: "Stop trying to push a quantum ball over a hill. Instead, let the quantum wave flow naturally into its new home, and the math finally works."

Technical Summary: Relaxation Approach to Quantum-Mechanical Modeling of Ferroelectric and Antiferroelectric Phase Transitions

Problem Statement
Ferroelectric and antiferroelectric materials exhibit characteristic hysteresis loops driven by temperature- and electric-field-induced phase transitions. Traditional modeling of these phenomena relies on classical frameworks, typically treating transitions as Arrhenius-type activated processes occurring well within the classical regime. These models assume that a system in a metastable state must overcome an energy barrier via thermal activation, with transition rates proportional to $e^{-E_A/k_BT}$ .

However, the authors identify critical failures in this traditional approach when applied to first-principles simulations. Specifically, classical models predict unphysically high coercive and switching fields (e.g., predicting ~3018 kV/cm for PbTiO3 versus an experimental ~175 kV/cm) and fail to reproduce experimental hysteresis loops at zero Kelvin. The discrepancy arises because, in classical transition state theory, the system remains trapped in a metastable state until the energy barrier vanishes (infinitely long lifetime), whereas real systems switch at finite fields. Furthermore, standard Density Functional Theory (DFT) simulations are limited by small supercell sizes that restrict phase space (e.g., preventing domain formation) and often assume pristine materials, ignoring defects that lower barriers in real samples. While previous studies have suggested quantum effects (such as zero-point energy) play a role in materials like SrTiO3 and BaTiO3, it remained unestablished whether quantum mechanical treatments could resolve the specific controversies regarding hysteresis loops and switching fields.

Methodology
The authors propose a fundamentally different framework: modeling phase transitions as relaxational processes within a quantum mechanical regime, rather than as classical activated events. They develop and compare three specific models:

Ground-State Relaxation (GSR) Model:
- The system is described by a time-propagating quantum state $|\psi(t)\rangle$ .
- To simulate relaxation from an excited (metastable) state to the ground state, the energy eigenvalues $E_n$ are modified to include an imaginary component: $E'_n = E_n - i\hbar\gamma_n/2$ , where $\gamma_n$ is the relaxation rate.
- The relaxation rate is defined as proportional to the sum of Bohr frequencies associated with all lower-energy states: $\gamma_n = k_{scale} \sum_{E_m < E_n} (E_n - E_m)/\hbar$ .
- The polarization is calculated as the expectation value of the polarization operator. This model is limited to zero Kelvin.
Density Matrix Based (DMB) Model:
- To incorporate finite temperatures, the authors utilize the density matrix formalism $\rho(t)$ .
- The time evolution follows a modified Liouville equation that includes a relaxation term driving the system toward the equilibrium density matrix $\rho_{eq}$ (canonical ensemble) rather than just the ground state.
- The equation accounts for both population decay and gain, normalized by the Boltzmann factor $e^{-(E_m - E_n)/k_BT}$ .
- A trace-conserving term is added to ensure physical validity.
Lindblad Operators (LO) Model:
- This approach implements the relaxation within the Lindblad formalism, using master equations with Lindblad operators $L_{kj}$ to describe interactions with a thermal bath.
- This method explicitly preserves the trace of the density matrix and satisfies the condition of detailed balance, ensuring that the ratio of transition rates between states matches the Boltzmann distribution at equilibrium.

Key Inputs and Parameters

Energy Landscapes: Computed using DFT (VASP) for materials including PbTiO3, PbZrO3, PbHfO3, CsGeBr3, HfO2, and AlN. The energy profiles $U(P)$ are derived from linear interpolation between phases, representing homogeneous polarization reversal (nucleation of a nanodomain).
Mass and Charge: The effective mass of the polar mode is derived from DFT-calculated harmonic frequencies. The Born effective charge ( $Z^*$ ) converts polarization to supercell distortion.
Relaxation Rates: While analytical expressions exist (Fermi's Golden Rule), the rates are fitted to experimental data (IR/Raman linewidths) via a scaling factor $k_{scale}$ .
Dynamics: Intrinsic dynamics (Bohr frequencies) are largely turned off in the primary simulations, as the authors argue their effects average out on the relaxation timescale.

Results

Antiferroelectrics (PbZrO3): The GSR model successfully reproduces the antiferroelectric hysteresis loop and predicts switching fields comparable to experimental values (315 kV/cm) at 0 K, a result unachievable with classical models. The DMB and LO models extend this to finite temperatures (10–1000 K).
- The DMB model predicts step-like polarization evolution and square hysteresis loops at low temperatures.
- The LO model predicts exponential evolution and slanted hysteresis loops, showing better agreement with the specific experimental shape of PbZrO3 loops.
Ferroelectrics (PbTiO3): The models reproduce ferroelectric hysteresis loops and correctly predict the trend of significantly higher coercive fields for non-perovskites (AlN, HfO2) compared to perovskites.
Temperature Dependence: The models capture the slimming of hysteresis loops as temperature increases, transitioning from antiferroelectric to paraelectric behavior.
Coexistence Field Correction: The authors note that DFT underestimates the coexistence field (the field where polar and antipolar phases have equal energy). To align simulations with experimental loop widths, they adjusted the AC field frequency used in the simulations for antiferroelectrics.

Significance and Claims
The paper claims that treating (anti)ferroelectric phase transitions as quantum mechanical relaxation processes overcomes the limitations of traditional classical, Arrhenius-based models. The primary significance lies in the ability to perform efficient, first-principles-based simulations that yield switching fields and hysteresis shapes in good agreement with experiments, even at zero Kelvin where classical thermal activation is impossible.

The authors assert that this framework demonstrates the essential role of quantum mechanics in transitions long regarded as purely classical. They suggest the approach is applicable to a broad range of phase transitions (magnetic, elastic, multiferroic, electronic) and other phenomena involving quantum tunneling and chemical reaction rates. The work provides a new perspective for understanding the underlying quantum-mechanical features of phase transitions, specifically highlighting that the "unconventional" assumption of relaxation-driven quantum evolution is necessary to resolve discrepancies in ferroic material modeling.

Technical Summary: Relaxation Approach to Quantum-Mechanical Modeling of Ferroelectric and Antiferroelectric Phase Transitions

More like this