Numerical treatment of a transient, nonlinear, metal-ferroelectric insulator-metal system using a finite element/boundary element numerical code

The recent interest in capacitor structures, having as dielectric a series of compounds with distorted perovskite topologies of both decentration of ion B in oxygen octahedra and, with gathering of oxygen octahedra, has led to the need for a better quantitative macroscopic evaluation, not only for nonlinearity and anisotropy of such ferroelectric structures, but also for the transient behavior of their electric field distribution. This paper is the first in a series attempting to use advanced Finite Element Method (FEM) and Boundary Element Method (BEM) numerical codes, as well as the latest technology in massive numerical processing, to give a quantitative picture of the transient phenomena in ferroelectrics with distorted, perovskite topologies as dielectrics between layers of metal, starting, initially, only with the treatment of the nonlinearity in transient regimes and adding, later, the anisotropy and tensorial approach, and trying to reconcile the experimental and theoretical data with the numerical results