Hypothesis and Theory ARTICLE
Programmable ferroelectric tunnel memristor
- Max Planck Institute of Microstructure Physics, Halle (Saale), Germany
We report a programmable analog memristor based on genuine electronic resistive switching combining ferroelectric switching and electron tunneling. The tunnel current through an 8 unit cell thick epitaxial Pb(Zr0.2Ti0.8)O3 film sandwiched between La0.7Sr0.3MnO3 and cobalt electrodes obeys the Kolmogorov-Avrami-Ishibashi model for bidimensional growth with a characteristic switching time in the order of 10−7 s. The analytical description of switching kinetics allows us to develop a characteristic transfer function that has only one parameter viz. the characteristic switching time and fully predicts the resistive states of this type of memristor.
In analogy to the resistor, inductor, and capacitor the memristor (memory resistor) is the passive circuit element, which was until recently considered to be “the missing circuit element” [1, 2]. Memristance, the memristor's main characteristic, which has the dimension of a resistance, should be tunable and history-dependent. Almost 40 years after the theoretical description, memristive behavior was firstly demonstrated in nanoscale systems possessing ionic transport properties . Very recently genuine electronic effects, such as the direct tunnel effect in ultrathin ferroelectric films, so-called ferroelectric tunnel junctions (FTJ),  have been proposed as memristors [5, 6]. It has been shown that the resistance of a ferroelectric tunnel junction is indeed tunable and history-dependent. Nevertheless, programming the resistance of such a device in a designed way remains a major goal to be achieved.
In this letter, we investigate memristor characteristics and switching kinetics of a memristor based on a ferroelectric tunnel junction that has been shown to obey proper quantum tunneling [7, 8] and not thermionic injection, as it has been shown on ferroelectric semiconductors . We demonstrate that the memristive states are univocally related to the polarization switching and that every memristor state can be achieved by design. Based on the analytical description given by the Kolmogorov-Avrami-Ishibashi [KAI] model [10, 11] we are able to explain the memristive behavior as being a result of coexisting, parallel-connected ferroelectric domains of opposite polarity inside 3600 μm2 capacitor devices. This said, in principle every ferroelectric tunnel junction with a bigger capacitor area than the nuclei size of its ferroelectric domains is intrinsically a memristor, irrespective of the injection mechanism.
Materials and Methods
According to the definition given by Chua  the time-dependent characteristics of a charge-controlled memristor can be described by:
where v(t) is the output voltage characteristics, M the memristance (with the dimension of a resistance) and the time-dependency given through the charge q(t) and the current i(t). In the very specific situation where M = constant, the equation describes the behavior of a simple resistor.
In the present case of a ferroelectric tunnel junction, the electronic tunneling probability through the ferroelectric layer and thus the resistance is controlled by the ferroelectric polarization, which is equivalent to a surface charge at the electrodes . Therefore, we can state that a FTJ is a charge-controlled memristor and Equation (1) applies. The charge q(t) can now simply be replaced by the time-dependent polarization P(t):
The memristive state of a FTJ is therefore dependent on the ferroelectric polarization and its time-dependent evolution is univocally related to the time dependence of the polarization and it can in principle be predicted. This is based on the working principle of FTJ, i.e., the electric current through a metal-ferroelectric-metal heterostructure is polarization dependent through the band alignment at the metal-ferroelectric interface [13, 14]. On the other hand, the time-dependent function P(t) is the result of the ferroelectric switching process. The polarization switching kinetics are usually described in terms of nucleation and growth of ferroelectric domains. The most common theoretical description of ferroelectric switching kinetics is based on the growth of nuclei that are randomly distributed across the ferroelectric capacitor  by ferroelectric domain wall propagation [16, 17]. After an initial nucleation process that takes place in the sub-picosecond range, ferroelectric domains grow laterally until they coalesce and the entire polarization will be in the opposite direction . The time dependence of the switching polarization P is given by the KAI equation: 
with (t) = P(t)/Ps being the normalized switched polarization, where Ps is the saturated value of polarization, τ is the characteristic switching time and n the domain switching dimensionality, which describes the domain wall movement. The exponent n can have ideally only integer values of 1, 2, and 3 representing one-, bi-, and three-dimensional domain growth during the switching process. From Equations (2) and (3) the memristance characteristics should be given by the following function:
In other words the memristance would be directly proportional to the amount of polarization switched by an applied external pulse of time interval t.
We have applied the above concept to ferroelectric tunnel junctions based on epitaxial oxide layers. The active layer of the structure is an ultrathin lead zirconate titanate Pb(Zr0.2Ti0.8)O3 (PZT) film epitaxially grown on a La0.7Sr0.3MnO3 (LSMO)/ SrTiO3 (STO) (100) substrate. A high-resolution TEM image of part of the structure is shown in Figure 1. Rather large-area capacitor-like active devices (60 × 60 μm2) were defined by thermal evaporation of cobalt through a shadow mask.
Figure 1. High-resolution transmission electron micrograph of the Co/PZT/LSMO tunnel junction used in the present study. The thickness of PZT is 8 unit cells (about 3.2 nm).
A typical current-voltage characteristic of a FTJ device is shown in Figure 2. This figure shows for the first time a hysteretic behavior similar to the i–v data of a TiO2-based memristor discussed by Strukov et al.  The switching voltages, viz. the voltages at which the resistance switches from the OFF to the ON state and vice-versa, correspond to the coercive fields of the ferroelectric thin film,  while the above dynamic current-voltage characteristic (see Figure 2) is not necessarily identical to the memristor characteristics. The latter should reveal the memristance dependence on time as for example given by Equation (4). Assuming that the memristor states in the FTJ are univocally defined by the polarization switching we may presume that switching kinetics of a FTJ is identical with the switching kinetics of polarization in ultrathin ferroelectric films.
Figure 2. Current-voltage characteristic of the Co/PZT/LSMO structure shown in Figure 1. The switching from RON to ROFF states and vice-versa is correlated to ferroelectric switching .
In order to experimentally study this phenomenon we have determined the memristor characteristics of a FTJ by measuring the resistive states induced by voltage pulses with variable width and amplitude. As the memristance is a history-dependent quantity, it is desirable to set the memristor device into a defined state before setting any other intermediate state. In the present case we used a negative voltage pulse with constant amplitude and width as a reset pulse, which sets the FTJ into the ON state defined by the lowest achievable tunnel resistance RON. From this reference state the FTJ can be driven into any intermediate state between RON and ROFF using a writing pulse, i.e., a positive voltage pulse with variable amplitude and width. This can be done repeatedly until the polarization is fully switching into the opposite direction driving the junction into the OFF state defined by the highest achievable ROFF resistance. Figure 3A shows the memristor characteristics determined using pulses of constant amplitude and variable width, following the scheme in Figure 3B. Between the RON and ROFF states there is in principle an infinite number of memristance states. The time scale of the writing pulses to access these intermediate states is in the sub-μ s range, which is comparable with the switching time of ferroelectric polarization that has been shown in previous studies . While the absolute values of the saturated (RON and ROFF) tunnel resistances might slightly vary, mostly due to additional conducting paths represented by defects,  the RON/ROFF ratio remains constant irrespective of the write amplitude and/or width.
Figure 3. Memristor characteristics of a Co/PZT/LSMO tunnel junction with PZT thickness of 8 u.c. (A) Tunneling resistance as function of the writing pulse width. The black data points were measured after each reset pulse, the gray data points after each writing pulse. (B) The applied voltage train of the measuring setup: in this case the amplitude and width of the reset pulse were −6 V and 200 ms, respectively, the writing pulse had the same amplitude as the reset pulse. The measurement of the tunnel current is performed with an applied voltage of 0.1 V. (C) Normalized ON/OFF ratios for different amplitudes of the write pulse. The lines are fits with Equation (9).
To normalize the measured data, we define RN as:
where RN takes values between 0 and 1. If the ROFF/RON ratio or the tunneling electro resistance value is very large, i.e., larger than 100, Equation (5) can be simplified to RN = (R − RON).
In Figure 3C, RN is plotted vs. the applied pulse time for different amplitude values of the applied pulses. For higher pulse voltages the system switches faster into the saturated ROFF state, which we contribute to the field-dependent domain wall velocity , leading to a faster unipolar formation of the domains in the ferroelectric tunnel barrier.
We now postulate that a memristive state can be described by a resistance resulting from a parallel connection of individual resistances given by coexisting ferroelectric domains of opposite polarity inside one capacitor similar to the case of a multilevel ferroelectric memory  (Figure 4). In other words, the effective resistance R of the tunnel junction after applying a certain pulse is given by the equivalent resistance of a circuit of resistors connected in parallel, which are defined by the area of domains with polarization up and down, R↑ and R↓ respectively. The areas which define R↑ and R↓ are complementary and the growth of one is on the expense of the other and happens according to the KAI dependence Equation (3). We can now write the reciprocal value of the resistance as a function of the normalized polarization P(t)/Ps = (t), which represents the time evolution of the switched domains by an applied pulse:
Figure 4. (A) A schematic cross section of the ferroelectric layer with visible domains. The red areas inside the sketch visualize the domains with pointing down polarization, marked by green arrows. The bright yellow area is the “unswitched” area and defined by polarization pointing up (blue arrows). The black arrows indicate the growth direction of the ferroelectric domains with opposite polarity. (B) Simplifying (A) to a model of parallel connected resistances, where R↑ is the summarized resistance of areas with polarization up and R↓ with polarization down accordingly.
From Equations (5) and (6) the time dependence of the normalized resistance can be written as (see also supplementary online material):
In the case of very large TER, Equation (7) simplifies to:
Equations (7) or (8) together with Equation (3) represent the characteristic transfer function of the ferroelectric memristor:
The experimental data can be fitted with Equation (9) in order to obtain the chatacteristic value of the switching time τ. As we know that in ultrathin epitaxial ferroelectric films the polarization switching is mostly bi-dimensional, the value of the parameter n in the above relation can be assumed to be 2, and thus one can constrain the fitting only to one parameter τ. As it can be easily seen in Figure 3C, the fit of the experimental data is very good which we interpret as the univocal result of ferroelectric domain wall kinetics. The value of the characteristic switching time is in the range of sub-μ s wich is, as expected, similar to the characteristic switching time in thin ferroelectric epitaxial PZT films of the same composition .
The analytical expression (9) can be used to predict or—better to say—to program any desired memristor state, independent from the memristor history. To show this experimentally, we analyzed the response of a ferroelectric memristor by successively applying programming pulses of constant amplitude, starting from a well-defined reference state ON, using different pulse width tA = 800 ns and tB = 1600 ns. In Figure 5, we present the time evolution of the resistance of a junction, where it can be easily seen that memristor states marked by horizontal lines can be achieved either with a long pulse tB or with a larger number of shorter pulses tA demonstrating the ability of programming any intermediate state between ROFF and RON.
Figure 5. Memristor programming. The data points represent the resistance after switching from a pre-defined ON-state with a pulse width of tA = 8 · 10−7 s (red) and tB = 16 · 10−7 s (blue). The vertical lines between the data points are the indicators for pulses that were applied. The horizontal lines are the programmed memristor states. The long time relaxation after applying the switching pulse is related to back-switching and relaxation of polarization .
According to the model, applying one programming pulse with the total width tt should always drive the FTJ to the same resistance as a sequence of pulses with width tf < tt of which the total time is tt = ∑n tf where n = tt/tf.
In summary, we have shown that a ferroelectric tunnel junction memristor based on epitaxial ferroelectric oxide ultra-thin films can be programmed in a designed way. The characteristic transfer function of this device is directly related to the polarization switching mechanism, viz. to the Kolmogorov-Avrami-Ishibashi model with two-dimensional domain growth. The characteristic transfer function of the FTJ memristor is dependent only on the characteristic switching time of the ferroelectric polarization which is in the present case of the order of 10−7 s. The physical model of MTJ-based memristors proposed above satisfies the mathematical equations describing the general features of memristors.
In comparison to other memristor devices,  the present device is based on a genuine electronic effect, respectively on quantum electronic transport and does not rely on ionic transport or phase change. The writing process, i.e., the partial switching of ferroelectric polarization, is achieved with rather large voltage approaching the coercive voltage of the ferroelectric layer, but still sufficiently low to comply to the nowadays microelectronic requirements. The current-based nondestructive readout, although performed at low voltages (<0.2 V) provides sufficient current density to allow downscaling and high integration density. However, further investigations have to be performed in order to identify and understand issues related to reliably, polarization back-switching, fatigue phenomena, etc.
Andy Quindeau, PhD student at the Max Planck Institute of Microstructure Physics. Dietrich Hesse, Professor at the Max Planck Institute of Microstructure Physics. Marin Alexe, Professor at Warwick University.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
This research was financed by the German Research Foundation (DFG) via SFB762. The authors acknowledge D. Pantel and S. Götze for the TER sample.
8. Pantel D, Lu H, Goetze S, Werner P, Kim DJ, Gruverman A, et al. Tunnel electroresistance in junctions with ultrathin ferroelectric Pb(Zr0.2Ti0.8)O3 barriers. Appl Phys Lett. (2012) 100:23. doi: 10.1063/1.4726120
15. Gruverman A, Wu D, Scott JF. Piezoresponse force microscopy studies of switching behavior of ferroelectric capacitors on a 100-ns time scale. Phys Rev Lett. (2008) 100:097601. doi: 10.1103/PhysRevLett.100.097601
16. McQuaid RGP, McMillen M, Chang L-W, Gruverman A, Gregg JM. Domain wall propagation in meso- and nanoscale ferroelectrics. J Phys Condens Matter. (2012) 24:024204. doi: 10.1088/0953-8984/24/2/024204
17. Yang SM, Jo JY, Kim DJ, Sung H, Noh TW, Lee HN, et al. Domain wall motion in epitaxial Pb (Zr, Ti) O 3 capacitors investigated by modified piezoresponse force microscopy. APL (2008) 92:252901. doi: 10.1063/1.2949078
Keywords: memristor, KAI, ferroelectric switching, ferroelectric tunnel junction, multiferroic oxides
Citation: Quindeau A, Hesse D and Alexe M (2014) Programmable ferroelectric tunnel memristor. Front. Physics 2:7. doi: 10.3389/fphy.2014.00007
Received: 23 October 2013; Accepted: 03 February 2014;
Published online: 25 February 2014.
Edited by:Nicholas X. Fang, Massachusetts Institute of Technology, USA
Reviewed by:Byungchan Han, Daegu Gyeongbuk Institute of Science and Technology, South Korea
Pavel Lukashev, University of Nebraska - Lincoln, USA
Copyright © 2014 Quindeau, Hesse and Alexe. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Marin Alexe, Max Planck Institute of Microstructure Physics, Weinberg 2, D-06120 Halle (Saale), Germany e-mail: email@example.com