Conduction Mechanism of Valence Change Resistive Switching Memory: a Survey

Resistive switching effect in transition metal oxide (TMO) based material is often associated with the valence change mechanism (VCM). Typical modeling of valence change resistive switching memory consists of three closely related phenomena, i.e., conductive filament (CF) geometry evolution, conduction mechanism and temperature dynamic evolution. It is widely agreed that the electrochemical reduction-oxidation (redox) process and oxygen vacancies migration plays an essential role in the CF forming and rupture process. However, the conduction mechanism of resistive switching memory varies considerably depending on the material used in the dielectric layer and selection of electrodes. Among the popular observations are the Poole-Frenkel emission, Schottky emission, space-charge-limited conduction (SCLC), trap-assisted tunneling (TAT) and hopping conduction. In this article, we will conduct a survey on several published valence change resistive switching memories with a particular interest in the I-V characteristic and the corresponding conduction mechanism.


Overview of Resistive Switching Memory
Typically, the resistive switching memory is a metal-insulator-metal (MIM) structure built with two metallic electrodes sandwiching an oxide layer.These bistable MIM cells switches back and forth between high resistance state (HRS) and low resistance state (LRS) when an appropriate voltage is applied across the electrodes.This resistance change phenomena is generally attributed to the formation and rupture process of the conductive filament in the oxide as a result of electrode redox reactions and nanoionic transport [66]. Figure 1 illustrates schematic diagram of the switching mechanism in VCM type resistive switching memory.There are two types of switching mode, namely, unipolar and bipolar.For unipolar, both CF formation and rupture process take place under the same electrical polarity, whereas bipolar requires opposite electrical polarities as shown in Figure 2a,b.Additionally, compliance current (ICC) is often applied to avoid excessive current flow and to prevent the irreversible hard breakdown.Devices with same materials can exhibit different switching mode depending on other factors such as electroforming, doping concentration, stoichiometry and physical parameters [67].ECM (cation-based) and VCM (anion-based) memories exhibit bipolar switching model, whereas TCM (both cation and anion-based) memory exhibits unipolar switching model [66].For high density implementation, resistive switching devices can be arranged in crossbar architecture which has minimum memory cell area of 4F 2 [68,69].Figure 3 shows three-dimensional modeling of a 2-by-2 resistive random access memory (ReRAM or RRAM) implemented in crossbar structure.

Typical Conduction Mechanism in Resistive Switching Memory
The electron transport in a MIM structure takes place via series of mechanisms.Some of these conduction mechanisms rely on the electrical property at the electrode-dielectric interface, e.g., energy barrier height of the interface and conduction carriers in dielectric films; the others depend on the properties of the dielectrics itself, e.g., trap level, trap spacing, trap density, carrier drift mobility, dielectric relaxation time and density of states in conduction [70].Among the electrode-limited conduction mechanisms are (i) Schottky emission; (ii) Fowler-Nordheim (F-N) tunneling; and (iii) direct tunneling.Bulk-limited conduction mechanism includes (i) Poole-Frenkel (P-F) emission; (ii) Ohmic conduction; (iii) SCLC; (iv) ionic conduction; (v) hopping conduction; and (vi) trap-assisted tunneling (TAT).

Schottky Emission
The Schottky or thermionic emission happens when thermally-activated electrons injected over the energy barrier into the conduction band of the oxide as shown in Figure 4 [55].This type of thermionic emission is one of the most often observed conduction mechanism in oxide especially in relative high temperature [71].
where m* is the electron effective mass in the oxide, k is Boltzmann's constant, T is the absolute temperature, h is Planck's constant, E is the electric field across the oxide, ΦB is the junction barrier height and ε is the permittivity of the oxide.In the case where the mean free path of the electrons is less than the Schottky barrier width, Simmons' modified Schottky equation (Equation ( 2)) is more appropriated [72,73].
where α is a constant, μ is bulk mobility and m0 is electron mass.

Fowler-Nordheim (F-N) and Direct Tunneling
In the presence of high electric field, the Fowler-Nordheim (F-N) and direct tunneling across the oxide may take place.Generally, direct tunneling is more dominant in oxide thinner than 3 nm [74], while F-N tunneling is usually more dominant at thicker oxide [71,75] and has been well accepted as charge changing mechanism in floating-gate memory [74].Figure 5 and 6 show the energy band diagram for F-N and direct tunneling respectively.
The current density due to F-N mechanism can be represented with Equation (3) [76]: where m* is the electron effective mass in the oxide, h is Planck's constant, and ΦB is the junction barrier height; whereas, the current density due to direct tunneling can be approximated to Equation (4) [71]: where κ is the relative dielectric constant of the oxide layer and tox,eq is the equivalent oxide thickness (EOT).Note that the current density (JDT) for a given oxide thickness is dependent to the κ value of the dielectric and tunneling barrier height (ΦB) [77].

Poole-Frenkel (P-F) Emission
The P-F emission happens when trapped electrons get excited into the conduction band of the oxide.The electric field decreases the Coulombic potential barrier of the electrons and subsequently increases its probability for being thermally excited out from the traps [70].The exponential portion of P-F expression is very similar to Schottky emission, except the junction barrier height (ΦB) is replaced with the depth of traps' potential well (ΦT) and the barrier lowering effect in P-F is double of Schottky emission due to immobility of positive charge [78,79].The schematic energy band diagram of P-F emission is depicted in Figure 7.
The current density JPF can be expressed as Equation ( 5) [71]: where µ is the electronic drift mobility, NC is the density of states in conduction band, E is the applied electric field, ΦT is the depth of traps potential well, k is Boltzmann's constant and T is the absolute temperature.

Space-Charge-Limited-Conduction (SCLC)
A trap-controlled SCLC is constituted of three portions: The Ohmic region (I ∝ V), the Child's square law region (I ∝ V 2 ) and the steep increase in the high field region [80].This mechanism could be easily identified whenever Ohmic conduction is observed at low field, followed by power-law dependence (I ∝ V 2 ) of the current density is observed in high field [80][81][82].In low field regime, the conduction mechanism is dominated by the thermally generated free electrons in the oxide film.When the field exceeds the square-law onset voltage, the density of electrode-injected electrons gradually exceeds the equilibrium concentration and dominating the conduction [43].The I ∝ V 2 region is known as trap-unfilled SCLC region, while the third region is often referred as trap-filled SCLC region.Hence, the probability of SCLC conduction is usually higher if the electrode's contact is highly carrier injective [52].According to the SCLC theory [83], for the simple case of a single discrete trapping level, the current density at high field is given by Equation ( 6) [84]: where εi is the permittivity of the oxide, µ is the mobility, θ is the ratio of free and shallow trapped charge, and d is thickness of the oxide.

Ionic Conduction
Movement of ions under influence of an applied electric field attributes to ionic conduction which is similar to a diffusion process [78].The ionic current density is proportional to ion drift velocity and the general expression is given by Equation ( 7) [85,86]: where v is the Debye frequency, r is the jump distance, ∆G ≠ is the free activation enthalpy and E is the electric field.In the case of low electric fields (E ≪ kT/re), the current density is linearly related to electric field and can be approximated to Equation ( 8): For high field (E ≫ kT/re), the current density is exponentially related to electric field as shown in Equation ( 9):

Ohmic Conduction
There are a small number of mobile electrons that get generated due to thermal excitation, even though the energy band gap of the oxide is by definition large.These electrons attribute to the Ohmic conduction in which the current density is proportional to the electric field and is relatively obvious in HRS.Equation (10) represents the current density of the Ohmic conduction.
where σ is electrical conductivity, µ is the electron mobility, Nc is the effective density of states of the conduction band, EC is the conduction band and EF is the Fermi energy level.

Nearest Neighbor Hopping (NNH)
NNH also known as fixed range hopping (FRH).Oxide layer often contains defects known as traps.Electrons that trapped in one trap site can hop into other nearest trap sites through tunneling effect.The current density due to this mechanism can be represented by Equation (11) [87].
where σ0 is the electrical conductivity at temperature of T0.

Mott Variable Range Hopping (VRH)
In addition to NNH, Mott VRH is another widely studied hopping conduction mechanism.Mott VRH describes the electron hopping behavior in strongly disordered systems with localized charge-carrier states [87], and has been suggested to be the main conduction mechanism in the semiconducting CFs [22].In contrast with NNH, where the trap electrons are expected to hop into the nearest neighboring trap, the electrons under the VRH scheme may hop into traps that are further away but have lower trap energy.In term of the I-V characteristic, the Mott VRH current is similar to Ohmic conduction where it is proportional to applied electric field.However, the conductance (σ) of Mott VRH demonstrates temperature dependency (ln(σ) ∝ T −1/4 ) as shown in Equation ( 12) [87]: where σ0 is the electrical conductivity at temperature of T0.

Trap-Assisted Tunneling (TAT)
TAT is attributed to the tunneling currents assisted by the defects in the dielectric [88,89].In contrasts to single step tunneling process like F-N tunneling or direct tunneling, the traps (e.g., oxygen vacancies) in oxide layer assist the electrons to tunnel from cathode to anode through a 2-steps process [90].The electrons are first captured from the cathode into the traps and subsequently emitted to anode.For inelastic TAT, the electron will relax to the trap's energy level by emitting one or more phonons with the energy ℏω [91].The generalized expression of the TAT current is quite similar to the F-N tunneling except for the E 2 factor, and can be written as in Equation ( 13) [92]: where A is a constant and ΦT is the energy of the electron traps with respect to the conduction band edge of the oxide.

Summary
Table 1 summarizes the basic conduction expression in oxide layer and highlights the electric field and temperature dependency.

Conduction Mechanism of Published Resistive Switching Memories
In this section, we will review selected published works where the dominant conduction mechanism are explicitly studied and reported.These works are categorized and listed in Tables 2-6 based on corresponding dominant conduction mechanism.Please note that the works reported are by far not exhaustive.The tables here only cover a subset of the huge collection of published resistive switching devices.Others transport models that requires iterative or numerical calculation, such as Menzel model [93], domain-based model [94,95], drift-diffusion model [96], ion-transport-recombination model [97][98][99] and hourglass model [100,101] are deliberately omitted in this article.

Poole-Frenkel (P-F) Emission
Devices with conduction dominated by P-F emission should exhibits linear fitting when the I-V characteristic is plotted on ln(J/E) vs. (E 1/2 /T) graph.This conduction mechanism has been widely observed in HRS of several resistive switching devices that based on materials such as SnOx [24], ZnO [25,26], AlOx [27], CeOx [28], WOx [29], LaHoO3 [30] and etc. as listed in Table 2. Figure 9a,b show the I-V characteristic of unipolar SnO2-based device with three different materials as top electrodes (Pt, Au and Ti) as reported by Nagashima et al. [24].As shown in Figure 9d, the I-V curve in HRS is well fitted to P-F emission model (ln(J/E) ∝ E 1/2 ).This indicates that the conduction mechanism of this 40 nm thick crystalline SnO2 is dominated by the electrons that hopping from the trap states into conduction band under high electric field.As reported by Lee et al. [25], comparable behavior is observed in the flexible ZnO-based device with Au and stainless steel as electrodes.In LRS, the J vs. E plot shows strong linear relationship with slope close to one, indicating that conduction mechanism in LRS is governed by Ohmic conduction.For HRS regime, it is observed that the conduction mechanism continue to be governed by Ohmic conduction at low field, while the nonlinear effect is getting dominance at high voltages (E > 0.1 MV/cm).This non-linear region has good agreement with the P-F emission model (ln(J/E) ∝ E 1/2 ).In a separate work, Chen et al. [26] also observed P-F emission on self-rectifying memory device based on same material (ZnO).The HRS current of the device, which using Al and Si as electrode, can be fit to the P-F emission with a slope of ~1.56.
Kim et al. [27] obtained excellent fitting for both P-F and Schottky emission in the nitrogen-doped forming-free AlOx device, because both mechanisms have similar linear dependency between current and the square root of voltage when measured under fixed temperature.Although Schottky emission has strong temperature dependency (J ∝ T 2 ) and could be used to differentiate it from P-F emission, however Kim et al. has to omit this option as infeasibility to increase the temperature of the test structure by few orders of magnitude.In turn, another approach was proposed, where the trap energy barriers (ΦT) of P-F emission model and Schottky barrier (ΦB) of Schottky emission model are extracted from respective slope in the fitted curve and compared with the expected values.It was found that the extracted barrier (ΦB) from Schottky emission plot is unreasonably low and the authors conclude that P-F emission is more appropriate description for the conduction mechanism in N-doped AlOx device.
Sharma et al. reported similar conduction mechanism in rare-earth based oxide memories, such as Pt/SmGdO3/Pt [102] and Pt/LaHoO3/Pt [30].In both cases, a linear relationship between ln(J/E) vs. E 1/2 confirmed the domination of Poole-Frenkel emission in OFF state.In the SmGdO3 sample, Ohmic conduction was observed at low voltage regime (<0.6 V) of the ON state.Whereas minor deviation from initial linearity was observed in high voltage regime (>0.6 V).The deviation in resistance is probably due to Joule-heating.
Similarly, the I-V characteristic of Cerium oxide-based memory reported by Ismail et al. [28] resembles the findings in [27], which have good linear fit to both P-F emission and Schottky emission.Anyhow, the observed current does not demonstrate strong correlation in the plot of ln(I/T 2 ) against 1/T, which is expected in Schottky emission.Therefore, the authors suggest that the conduction in this 15 nm thick CeOx device is dominated by P-F emission.
Overall, linearity relation in the ln(J/E) vs. E 1/2 plot is a well-accepted approach in identifying the P-F emission.Additional checks such as temperature dependency (ln(J/T 2 ) ∝ 1/T) and I-V symmetric check are usually employed to distinguish P-F emission from the electrode-limited Schottky emission.

Schottky Emission
As tabulated in the Table 1, current density attributed to Schottky emission has dependency on both temperature (T) and electric field (E).Hence, the existence of Schottky emission can be ascertained by fitting to two functions, which are: (i) ln(J/T 2 ) ∝ 1/T under fixed electric field; and (ii) ln(J) ∝ E 1/2 under fixed temperature.There are numerous published resistive switching devices that suggested Schottky emission as the dominant conduction mechanism, such as [33][34][35][36][37][38][39][40][41][42], etc.The combination of electrodes and materials of these works is listed in Table 3. Cao et al. [33] compared I-V characteristic of the unipolar switching Pt/TiO2/Pt device with several non-linear conduction mechanisms including P-F emission, SCLC and Schottky emission, and found that the experimental data fit perfectly to Schottky emission model with slope of 3.553 in the double-logarithmic plot of I-V.
Huang et al. [34] observed comparable conduction mechanism in the resistive switching Ti/TiO2/Pt diode which prepared using a simple evaporation process.The asymmetry of J-V curves as shown in Figure 10 is the result of unequal Schottky barriers at the Ti/TiO2 (0.13 eV) and the TiO2/Pt (0.73 eV) interfaces.Since Schottky emission is the electrode-limited conduction mechanism, the I-V characteristic of devices with heterogeneous electrode material tends to demonstrate asymmetric property, which can be used to differentiate Schottky emission from other bulk-limited conduction mechanism such as P-F emission.Similar asymmetrical property of I-V relation is also witnessed in Lee et al.'s experiment on TiN/Ti/TiOx/TiN structure and the double-layered TiN/Ti/TiOx/HfOx/TiN structure [35].It is clearly seen that the single layer structure demonstrates stronger symmetrical property in I-V characteristic compared to double layered structure.In the single layer structure, Schottky barrier at Ti/TiOx and TiOx/TiN interface are 0.17 eV and 0.3 eV respectively.Meanwhile, the Schottky barrier of additional HfOx/TiN interface in the double-layered structure is relatively higher at 2.36 eV.
Zhu et al. [36] observed that the anodic-AlOx device with homogenous Al electrode is fit to both P-F and Schottky emission expression.Since symmetrical check on I-V relation is not possible due to homogenous electrodes, Zhu et al. resorted to compare the activation energies (Ea) extracted from the respective Arrhenius plots, but found that the extracted values are close (0.41-0.43 eV) and indecisive.Thus, this is a rare case where no clear evidence on whether P-F or Schottky emission is governing the conduction.
Sun Kuo et al. employed Schottky emission to explain the increasing HRS resistance during cycling operation on the Al/VOx/Cu [38].Noted that Al electrode has lower work function (4.28 eV) thus can form a Schottky barrier with p-type VOx that has higher work function of 4.7 eV.A good linear relation of ln(J) and V 1/2 , indicating the Schottky emission is the dominant conduction mechanism in HRS.For LRS, the I-V fitting and conductive atomic force microscopy (CAFM) image suggest that Ohmic conduction is the dominant current.
Hsieh et al. [41] compared the I-V characteristic of double-layered ZnO/HfO2 memory with the P-F and Schottky emission fittings.The plots are found to have better fit with Schottky current model.
By and large, fitting of ln(J)-E 1/2 function with fixed temperature is broadly adopted in classifying Schottky emission.The additional temperature dependency check (ln(J/T 2 ) ∝ 1/T) is sometime omitted due to infeasibility in adjusting the temperature.In some cases where the top and bottom electrodes have relatively big difference in work function, asymmetry I-V characteristic could be observed.

Space Charge Limited Conduction (SCLC)
SCLC is another commonly observed conduction mechanism in published resistive memories.Signature of SCLC is relatively observant with two regions of distinct slope when the I-V characteristic is plotted in double-logarithmic scale.The I-V relation follows Ohmic law (I ∝ V) in low field regime, and Child's Law (I ∝ V 2 ) in high field regime.In HRS, the third region could be observed in which the current will increase rapidly due to the formation of CF.Table 4 lists some published works in which the SCLC characteristic is observed.
Yu et al. [43] reported that both LRS and HRS of the Al/TiOx/Al device is governed by the SCLC conduction mechanism.The Ohmic conduction behavior (I ∝ V) is observed in low voltage region where V < 0.1 V, and followed by a quadratic term (I ∝ V 2 ) as defined by Child's Law in V > 0.1 V.
The SCLC mechanism is also identified to be the dominant conduction mechanism in Au/Cr/Zr + -implanted-ZrO2/Si device as reported by Liu et al. [44].The fitting of each SCLC regions for both positive and negative voltage range are clearly shown SCLC behavior.By leveraging the expression for trap-controlled SCLC as listed in Equation ( 5), the role of implanted Zr + in resistive switching phenomenon is explained in the work.Pt/ZrO 2 /TiO 2 /Pt Ohmic SCLC Kim et al. [47] reported the observation of the SCLC conduction in Al/TiOx/TiO2/Al device.Three regions with distinct slope differences are noted in the double-logarithm I-V plots.The author also includes the P-F emission and thermionic (Schottky) emission curve in same plot for comparison.Fu et al. [51] exhibited that the conduction mechanism in La0.5Sr0.5CoO3(LSCO)-baseddevice is also agree with SCLC description.As shown in Figure 11, both positive and negative voltage range demonstrates the clear signature of SCLC mechanism, which has three regions with noticeable slope difference.The conduction characteristics of this perovskites-based device can be explained by the injected hole carrier and is depend on Ag/LSCO and Pt/LSCO interfaces traps that are exponentially distributed in energy.
Among the reviewed literatures, whenever SCLC behavior is observed in HRS, the same mechanism tends to dominate in LRS as well.However, there are few exceptions, such as [45,46,50,52,53] where only Ohmic conduction is observed in LRS.
In Cu/AlN/Pt device developed by Chen et al. [50], the square-law onset voltage observed in HRS is about 0.125 V.However, the current in LRS is already concealed at compliance current (Icomp) around 0.11 V.
In [45], both the ZnMn2O4 and ZnMnO3 samples have the square-law onset voltage around 2 V in HRS.The Vreset of ZnMnO3 sample is lower than this square-law onset voltage and the current shows a strong linearity relation with the voltage over the whole LRS range.In contrast, some minor nonlinearity is observed during the LRS to HRS transition in the ZnMn2O4 sample, which has slightly higher Vreset as depicted in Figure 12.Therefore, it is believed that the CF is ruptured and reset to HRS before the ZnMn2O4 device moved into trap-unfilled SCLC region (I-V 2 ).Wang et al. [46] reported that the SCLC is the dominant conduction mechanism in HRS and negative-based LRS of the Ti/CuxO/Pt device.However, there is no I-V 2 region observed in positive-biased LRS as shown in Figure 13, which could be capped by the current compliance imposed on the device.Similarly, for the Zr/CeOx/Pt memory device [52], the Ohmic conduction in LRS can be explained by the current compliance limit in negative voltage region and low Vreset in positive voltage region as shown in Figure 14.Generally, observation of three apparent regions where (i) I ∝ V, (ii) I ∝ V 2 and (iii) steep increase of I, signifies the dominant of SCLC mechanism.To summarize, SCLC conduction has been observed in wide range of materials and it is commonly accepted that the defect-related traps, e.g., oxygen vacancies, are contributing to the SCLC conduction.

Trap-Assisted Tunneling (TAT)
I-V characteristic of TAT conduction could be fit linearly in the plots of ln(J) versus 1/E.Besides this, microscopic models such as multi-phonon TAT (MTAT) model are often used.The TAT conduction modeling is found to be consistent with some experimental data based on HfOx [54][55][56], TaOx [57] and SiOx [58,59] as depicted in Table 5. Bersuker et al. [54] analyzed the conduction mechanism in TiN/HfO2/TiN device by using the statistical multi-phonon trap-assisted tunneling (MTAT) model.The model is a quantitative physical model that takes into account of the lattice relaxation due to electron-phonon coupling as well as the microscopic structural changes associated with the concentration of the defects [104,105].The simulated result has good agreement with experimental data when voltage >0.8 V and temperature between 300 K and 400 K.
Yu et al. [55] used process of elimination to pinpoint the dominant conduction mechanism in the TiN/HfOx/Pt device which fabricated using atomic layer deposition (ALD) method [106].The P-F emission option was dismissed due to nonlinear curve on ln(I/V) vs. V 1/2 plot as well as unrealistic extracted trap energy (ET) and dielectric constant (εr).Schottky emission and F-N tunneling were also ruled out due to weak dependency on temperature and electric field respectively.Direct tunneling was also excluded as it should be negligible in the 10 nm thick HfOx film.Along these lines, the remaining probable conduction mechanism is TAT.The hypothesis was then evaluated by using current modeling that takes into account of the electron transition rate and transmission probability that is calculated with Wigner-Kramers-Brillouin (WKB) approximation.Several follow-up models from the same group, i.e., stochastic model [107] and Kinetic Monte Carlo (KMC) model [108], further corroborated the TAT as dominant conduction mechanism.
Mehonic et al. [58] reported a comparable conduction mechanism in 35 nm thick SiO2-based device but utilizing a much compact approach which the current density is fitted to an expression with respect to electric field and trap's energy barrier (ΦT) [92] as shown in Equation 10.The negative linear relation in the ln(J) vs. 1/E plots clearly exemplifies the TAT conduction in both LRS and HRS.The extracted trap energy barrier (ΦT) of HRS and LRS are 1.06 eV and 0.41-0.63eV respectively [59].
Long et al. [56] affirmed the role of Ohmic conduction during LRS in 6 nm thick HfO2-based device.However, the conduction in HRS is found to be governed by combination of both P-F emission and TAT process.As shown in Figure 15, the experimental data is close to TAT fitting on 125-160 K, while P-F emission seen dominating at 265-365 K.
Similarly, TAT conduction is observed in addition to P-F emission in single layer Pt/TaOx/Pt device as reported by Jeon et al. [57].However, the energy of electron traps ΦT that extracted from plain P-F emission fitting is not consistent with value observed using spectroscopic ellipsometry (SE).Thus, the most appropriate conduction mechanism in this Pt/TaOx/Pt device is probably TAT.
When comparing to other conduction mechanisms, the usage of expression fitting is not as popular in identifying TAT.Most researchers relate TAT conduction via process of elimination and microscopic modeling.Generally, it is accepted that oxygen vacancy defects near the cathode is playing a critical role in conduction mechanism [55,109].Hence, TAT has broadly employed to explain the electron transport mechanism in oxide thin film with high defect density.

Hopping Conduction
Nearest Neighbor Hopping (NNH) and Mott Variable-range Hopping (VRH) are the hopping conductions, which commonly observed in resistive switching memories.The devices dominated by these hopping conductions exhibit temperature dependency, i.e., ln(σ) ∝ T −1 for NNH and ln(σ) ∝ T −1/4 for VRH.Table 6 lists some of the published works which hopping conduction is reported to be the dominant conduction mechanism.Ho et al. [60] reported that two different characteristic of I-V behaviors are observed in the WOx-based device, i.e., linear in LRS and non-linear in HRS.The LRS conduction is believed to be governed by metallic conduction because of the Ohmic I-V characteristic and finite residual resistance at low temperature.In HRS, the I-V curve is well fitted to a straight line in ln(J) ∝ T −1/4 plot, suggesting that the current is primarily due to Mott VRH.A good hyperbolic-sine fitting of the I-V curve further agreed the role of Mott VRH in HRS current conduction.
Both Mott VRH and Fixed-range Hopping (FRH) are observed in the Ir/Ta2O5−σ/TaOx/TaN device as reported in [61].The I-V plot fitted the Mott VRH criteria in temperature lower than 200 K, and FRH criteria in the higher temperature.From the plots of ln(σ) versus 1/T, an obvious slope-switching points is observed around 200 K, implying the dominating current mechanism is switching toward Mott VRH from FRH when the temperature decreases beyond 200 K.
Phan et al. [62] verified that I-V characteristic of Cr-doped SrTiO3 thin firm is matching the Mott VRH model in positive biased regime and trap-controlled SCLC model in negative biased regime.Within the Mott VRH dominating region, conductance shows temperature dependency (ln(I) ∝ T −1/4 ) when bias voltage is below +1 V and electrical field dependency (ln(I) ∝ E −1/4 ) when bias voltage is above +1 V.
Similar to [61], Zhang et al. [64] also observed temperature dependency of the hopping conduction in TaOy/Ta2O5−x device.Two distinct slopes can be seen in the plot of ln(σ) against 1000/T.There is good fit to NNH expression (slope ~1) at temperature above 50 K, and VRH expression (slope ~1/4) at temperature below 50 K.
Fang et al. [65] reported that the resistance in both LRS and HRS of Pt/HfOx/TiN device increases linearly with respect to exp(1/kT), and tends to saturate at temperature below 77 K.The slope of this fitted curve is also being used to extract device parameters such as hopping barrier Ea based on the hopping conduction expression as in Equation ( 14) [71,110]: where a is the mean spacing between trap sites or the hopping distance, n is the electron concentration in the conduction band of the dielectric, v is the frequency of thermal vibration of electrons at trap sites and Ea is the barrier height between two traps.This indicates that the slope of ln(J) versus 1/T plots can be related to the difference of field-driven barrier lowering (qaE) and hopping barrier between traps (Ea).
Generally, resistive switching materials tend to have high density of traps, hence it allows trap electrons to hop from one to other traps.The conductance due to these hopping electrons demonstrates temperature dependency, i.e., ln(σ) ∝ T −1/4 for Mott VRH and ln(σ) ∝ T −1 for NNH or FRH.

Conclusions
We described and analyzed some commonly observed conduction mechanisms of resistive switching memories, namely, Poole-Frenkel (P-F) conduction, space charge limited conduction (SCLC), Schottky emission, trap-assisted tunneling (TAT) and hopping conduction.We found that the resistive switching memories do not demonstrate the same I-V characteristic even when built with the same materials.This could be due to other attributes such as electrode material, oxide properties (thickness, bandgap, trap energy level), fabrication process (annealing temperature and ambient), operating temperature and interfacial barrier height between metal and oxide.Generally, fitting to expression of typical conduction mechanisms is the most popular approach to identify the dominant conduction contributor.In some cases, additional analysis is carried out on the symmetrical property of the I-V characteristic and extracted physic parameters.However, the drawback of these techniques is that they rely on the assumption that electron transport is only governed by a single mechanism.This ignores the fact that multiple mechanisms can collectively contribute to the conduction depending on the temperature, applied electric field and voltage polarity.Thus, a generalized conduction modeling which incorporating concurrent contribution from major conduction mechanisms should be leveraged to enable accurate I-V characteristic modeling over various materials and operating scenarios.

Figure 1 .
Figure 1.Schematic of the switching mechanism of bipolar metal oxide resistive switching memory (a) Switching layer in pristine state before electroformed; (b) A low resistance filament is formed when forming voltage is applied to top electrode; (c) The filament is ruptured when a negative voltage is applied and device turns into high resistant state (HRS).

Figure 2 .
Figure 2. (a) Unipolar; and (b) bipolar switching.ICC is the compliance current to prevent permanent breakdown.

Figure 3 .
Figure 3. Resistive random access memory (RRAM) with resistive switching devices arranged in crossbar architecture.

Figure 6 .
Figure 6.Energy band diagram for direct tunneling.

Figure 7 .
Figure 7. Schematic energy band diagram of Poole-Frenkel emission in MIM structures.

Figure 8
illustrates the energy band diagram of NNH and VRH mechanism.

Figure 8 .
Figure 8. Schematic energy band diagram of Nearest Neighbor Hopping (NNH) and Variable-range Hopping (VRH) model in an MIM structure.

a
range Hopping (VRH) High field only, low field follows Ohm's Law; b A & B are constants.

Figure 9 .
Figure 9. J-E plots of SnO2 in (a) LRS and (b) HRS when using Pt, Au and Ti as top electrodes.(c,d) show the comparisons with Schottky emission model and Poole-Frenkel emission model, respectively.Reprinted with permission from [24].Copyright © 2009, American Institute of Physics.

Figure 10 .
Figure 10.J-V characteristics of a typical Pt/TiO2/Ti MIM diode.The thickness of TiO2 is 8 nm.Insets show the band structures of the MIM under different biases and the Schottky emission fitting of ln(J) vs. E 1/2 from −0.4 V to −4 V, respectively.Reprinted with permission from [34].Copyright © 2010, American Institute of Physics.

Figure 11 .
Figure 11.(a) Relation of log(I) vs. log(V) for the positive applied voltage range from 0 V to 2 V; (b) Relation of log(I) vs. log(V) for the negative applied voltage range form −2 V to 0 V. Reprinted with permission from [51].Copyright © 2014, American Institute of Physics.

Figure 13 .
Figure 13.I-V characteristics of both the positive and negative bias regions of the Ti/CuxO/Pt device plotted in a double-logarithmic scale.Both ON-state and OFF-state are included in this plot.Reprinted with permission from [46].Copyright © 2010, American Institute of Physics.

Figure 14 .
Figure 14.I-V curves for Zr/CeOx/Pt memory device are displayed in double-logarithmic scale.The linear fitting results in both ON state and OFF state: (a) positive-voltage region; and (b) negative-voltage region.Adapted from [52].

Figure 15 .
Figure 15.HRS fitting of HfO2 based device for F-P emission between 190 K and 365 K and TAT fitting between 80 K and 97.6 K. Reprinted with permission from [56].Copyright © 2012, American Institute of Physics.

Table 1 .
Expression and field/temperature dependency of typical conduction mechanism in dielectric film.

Table 2 .
Published works with P-F emission observed.

Table 3 .
Published works with Schottky emission observed.

Table 4 .
Published works in which space charge limited conduction (SCLC) is observed.

Table 6 .
Published works with TAT observed.