P-E Hysteresis Loop Simulator: Ferroelectric Polarization

simulator intermediate ~10 min
Loading simulation...
BaTiO3-class loop — P_s = 26 uC/cm2, P_r = 22.4 uC/cm2, E_c = 15 kV/cm

With default parameters mimicking barium titanate, the hysteresis loop shows a remnant polarization of 22.4 uC/cm2 and a coercive field of 15 kV/cm, typical of a soft ferroelectric ceramic.

Formula

P(E) = P_s * tanh((E +/- E_c) * sq / E_c) (phenomenological hysteresis model)
W_loss = oint P dE (energy dissipated per hysteresis cycle)
P_r = P_s * tanh(sq) (remnant polarization from squareness)

The Signature of Ferroelectricity

The polarization-electric field (P-E) hysteresis loop is the definitive fingerprint of a ferroelectric material. When you cycle an electric field from positive to negative and back, the polarization traces a characteristic loop rather than retracing a straight line. This irreversibility — hysteresis — arises from the energy cost of switching electric dipoles between two or more stable orientations within the crystal lattice.

Anatomy of the Loop

Three parameters define the loop: saturation polarization P_s (the maximum polarization when all domains are aligned), remnant polarization P_r (the polarization remaining at zero field), and coercive field E_c (the field needed to reduce polarization to zero). Materials like PZT have large P_s and moderate E_c, making them excellent actuators. Materials like BiFeO3 can have enormous P_s but also high E_c, making switching energy-intensive.

Minor Loops and Sub-Switching

When the applied field is less than the coercive field, the material traces a minor loop — a thinner, lenticular path inside the full hysteresis curve. Minor loops represent partial domain switching and are important in fatigue studies, where repeated sub-coercive cycling can gradually degrade the switchable polarization through defect accumulation at domain walls.

Energy and Applications

The area enclosed by the hysteresis loop equals the energy dissipated per cycle. For memory applications (FeRAM), you want a square loop with minimal loss. For energy harvesting, you want a large loop area to maximize the electrical energy extracted from mechanical cycling. This tension between loss minimization and energy extraction drives much of ferroelectric materials engineering.

FAQ

What is a P-E hysteresis loop?

A P-E hysteresis loop plots the electric polarization P of a ferroelectric material as a function of the applied electric field E. The loop's shape reveals key material properties: remnant polarization (P at zero field), coercive field (field needed to depolarize), and saturation polarization (maximum achievable P).

Why does the loop have hysteresis?

Hysteresis arises because domain switching is an irreversible process — energy is dissipated as domain walls move through the crystal, overcoming pinning sites and defects. The area enclosed by the loop equals the energy lost per cycle.

What is remnant polarization?

Remnant polarization P_r is the polarization remaining after the external field is removed. It represents the material's memory of its polarization state and is the basis for ferroelectric random-access memory (FeRAM).

How does loop squareness affect applications?

A more square (rectangular) loop means sharper switching between polarization states, which is desirable for digital memory applications. A slanted loop indicates gradual switching, often seen in relaxor ferroelectrics used for actuators and capacitors.

Sources

Embed

<iframe src="https://homo-deus.com/lab/ferroelectrics/pe-hysteresis/embed" width="100%" height="400" frameborder="0"></iframe>
View source on GitHub