Electro-Thermal Model of Threshold Switching in TaOx-Based Devices

Mar 15, 2017 - Pulsed and quasi-static current–voltage (I–V) characteristics of threshold switching in TiN/TaOx/TiN crossbar devices were measured...
0 downloads 11 Views 912KB Size
Research Article www.acsami.org

Electro-Thermal Model of Threshold Switching in TaOx‑Based Devices Jonathan M. Goodwill,*,† Abhishek A. Sharma,‡ Dasheng Li,† James A. Bain,‡ and Marek Skowronski*,† †

Department of Materials Science and Engineering, Carnegie Mellon University, Pittsburgh, Pennsylvania 15213, United States Department of Electrical and Computer Engineering, Carnegie Mellon University, Pittsburgh, Pennsylvania 15213, United States



S Supporting Information *

ABSTRACT: Pulsed and quasi-static current−voltage (I−V) characteristics of threshold switching in TiN/TaOx/TiN crossbar devices were measured as a function of stage temperature (200−495 K) and oxygen flow during the deposition of TaOx. A comparison of the pulsed and quasistatic characteristics in the high resistance part of the I−V revealed that Joule self-heating significantly affected the current and was a likely source of negative differential resistance (NDR) and thermal runaway. The experimental quasi-static I−V’s were simulated using a finite element electro-thermal model that coupled current and heat flow and incorporated an external circuit with an appropriate load resistor. The simulation reproduced the experimental I−V including the OFF-state at low currents and the volatile NDR region. In the NDR region, the simulation predicted spontaneous current constriction forming a small-diameter hot conducting filament with a radius of 250 nm in a 6 μm diameter device. KEYWORDS: tantalum oxide, Poole−Frenkel conduction, threshold switch, negative differential resistance, thermal runaway

1. INTRODUCTION The ever-increasing need for faster and more energy- and space-efficient ways to store data has driven the interest in passive crossbar memory arrays. This design offers excellent scalability with the smallest possible cell size of 4F2, where F is the minimum feature size, but suffers from the so-called sneak path problem where all cells along the word and bit lines corresponding to the selected cell are half-selected. The parasitic current flowing through these cells can cause read errors and severely limits the array size. One way to reduce the sneak path problem is to incorporate a highly nonlinear selection device in series with the memory element.1−3 Threshold switching devices are one of the candidates for this application.4,5 This report is focused on one class of such devices, namely, structures based on transition metal oxides with highly temperature dependent conductivity, and the underlying mechanism of their operation. The current−voltage (I−V) characteristics of a threshold switch resemble the letter “S” and consist of three branches. The branch with the lowest current is referred to as the high resistance OFF-state, which continues until the slope of the I− V becomes infinite at the threshold voltage (VTH). At higher currents, the slope of the I−V becomes negative, giving rise to the intermediate branch known as the negative differential resistance (NDR) region. The last region, the low resistance ON-state, is defined as the high-current region above the NDR region where the slope of the I−V becomes positive once again. © XXXX American Chemical Society

While the conductivity in OFF- and ON-states is thought to be mostly uniform, the NDR region exhibits formation of high current density filaments.6 A number of mechanisms have been suggested to give rise to NDR in different materials/devices. For example, in amorphous chalcogenides, the proposed mechanisms include double injection,7 field-induced filling of traps,8 energy gain by electrons hopping between defect sites,9 electric field-induced secondary phase nucleation,10 and even thermal runaway,11,12 among others. In addition to amorphous chalcogenides, threshold switching has been extensively studied in VO2 and NbO2 where it was generally interpreted as due to a field- or Joule heating-induced insulator-to-metal transition.13−18 There are only a few reports of threshold switching in dielectric oxides such as TiOx, TaOx, or HfOx, and the mechanism behind this phenomenon is not well established.19−21 Three recent reports on NbOx offer some insight into the proposed mechanism of threshold switching in dielectric oxides, namely, thermal runaway.22−24 The three publications report the quasi-static I−V’s of NbOx-based devices and fit the data with a Poole−Frenkel-type formula to obtain conductivity as a function of temperature and electric field. Simulations including self-heating showed that, at a certain voltage, the Received: December 23, 2016 Accepted: March 15, 2017 Published: March 15, 2017 A

DOI: 10.1021/acsami.6b16559 ACS Appl. Mater. Interfaces XXXX, XXX, XXX−XXX

Research Article

ACS Applied Materials & Interfaces

Figure 1. (a) Quasi-static I−V measured at 300 K for a device with a 7.5% O2 flow TaOx film. The inset shows the low voltage characteristics. (b) Low field conductivities of devices fabricated with three different flow rates of oxygen during sputtering. (c) A comparison between the pulsed I−V (triangles) and quasi-static I−V (continuous lines) at different stage temperatures for a 7.5% oxygen film and (d) pulsed I−V for the same film (triangles) versus the fit of the Poole−Frenkel formula (lines) at different stage temperatures. fraction is not the same as the percentage of oxygen in the TaOx layer, nor are these two proportional. The devices fabricated with these films will be referred to by the respective O2 percentages of the ambient in which they were sputtered. The lateral size of the devices was 5 × 5 μm, and they were patterned in a ground-signal-ground configuration allowing for short pulse measurements. The quasi-static I−V’s were obtained using a circuit consisting of the source and load resistor along with the device under test connected in series. The transient data were obtained using Time Domain Transmissometry.27 The details of the fabrication, device layout, and experimental setups can be found in our earlier publications.21,28−30 All data are plotted versus the voltage across the device rather than source voltage.

devices undergo thermal runaway. The estimated functional layer temperature increase right before the threshold was estimated to be 35 K24 and 150 K22 above the stage temperature of 300 K. Both of these are far below the insulator-to-metal transition in NbO2 at 1081 K.25,26 The devices discussed in these publications needed to be electroformed, a procedure that leads to the formation of a smalldiameter permanent filament consisting of a conductive reduced oxide. Neither the size nor the composition of the filament is known. Moreover, the inhomogeneous material does not allow for observation and/or simulation of the salient feature of a threshold switch, specifically, the spontaneous current constriction in a structurally uniform device. In this work, we present experimental data on I−V characteristics of TiN/TaOx/TiN devices as a function of stage temperature and deposition conditions of the TaOx film. These were simulated using a finite element model incorporating coupled heat and current flow. The model correctly describes the S-NDR characteristics and predicts the spontaneous current constriction to a small-diameter filament of elevated current density and temperature. These insights reveal some of the limitations of TaOx-based threshold switches as selectors.

3. RESULTS A typical experimental quasi-static I−V of a tantalum oxide device structure is shown in Figure 1a. These data were obtained at a stage temperature of 300 K on a 7.5% oxygen device with a series load resistor of 26.7 kΩ. The devices were not electroformed, and the current and dissipated power were intentionally limited to 1 mA and 3 mW to prevent permanent changes to the device. The initial characteristics remained unchanged throughout the entire series of experiments, which was confirmed by measuring the high resistance state characteristics before and after the transition to the NDR region. The I−V shows the OFF-state at low voltages/currents with gradually increasing slope until eventually becoming vertical at 10 V. The functional I(V) dependence in the OFFstate is discussed in detail below. At higher current, the slope of the I−V becomes negative as the device enters the NDR region. Between 9.5 and 3.5 V, the I−V shows a straight line corresponding to a transition between two points in the NDR region along the load line. Above the transition, the device current increases rapidly with small changes in voltage. The inset in Figure 1a shows the I−V at low voltages where the initial current dependence is linear. This slope was used to determine the low field conductivity for the simulations discussed in the section below. The slope of the linear part of the I−V as a function of the stage temperature and oxygen fraction is shown in Figure 1b.

2. MATERIALS AND METHODS Devices with a TiN/TaOx/TiN structure were fabricated on a silicon wafer with 1 μm of thermal oxide. The metal−insulator−metal structure consisted of 30 nm thick sputtered TiN electrodes and 150 nm of functional reactively sputtered TaOx layer. Since the functional layer was deposited on a nonplanar surface with a patterned bottom electrode, the relative variations in film thickness over the edges of the electrode are expected to be large. In order to minimize the effect of such variations, a relatively large thickness of the functional layer was chosen. The TiN electrodes were sputtered at room temperature with an argon flow rate of 60 sccm and a chamber pressure of 3 mTorr. The TaOx films were sputtered at room temperature with a chamber pressure of 3 mTorr, an oxygen flow of 3, 4.5, and 6 sccm, and a constant total flow of Ar and O2 at 60 sccm. This corresponds to 5, 7.5, and 10% of oxygen fraction in the sputtering ambient. This B

DOI: 10.1021/acsami.6b16559 ACS Appl. Mater. Interfaces XXXX, XXX, XXX−XXX

Research Article

ACS Applied Materials & Interfaces The three functional layers spanned a conductivity range of over 2 orders of magnitude at room temperature. The conductivity of each device shows thermally activated behavior with activation energies ranging from 0.26 ± 0.01 eV at low oxygen content to 0.44 ± 0.005 eV at high oxygen content. This suggests that the Fermi level in TaOx layers, which is located close to the middle of the band gap in a stoichiometric material, moves closer to the conduction band with decreasing oxygen content. The change is likely due to creation of traps associated with oxygen deficiency. The presence of traps in TaOx with energies of 0.2, 0.3, and 0.6 eV below the conduction band edge has been reported in thermally stimulated current measurements.31 Similar energies of traps have been extracted from fitting the leakage current in TaOxbased structures with the Poole−Frenkel formula.32,33 The effect of Joule self-heating on the I−V characteristics of a typical device can be assessed based on data shown in Figure 1c. All data in this figure have been obtained on the same 7.5% O2 device at three different stage temperatures. The continuous lines represent experimental quasi-static voltage sweeps at a ramp rate of 0.4 V/s, while the triangles denote pulsed I−V characteristics. The dwell time for each voltage step was orders of magnitude longer than the device thermal time constant τo (estimated in our earlier publication to be 3 μs),28 and the device was able to reach a steady state temperature at every point on the curve. In contrast, the pulsed I−V’s were measured with 30 ns long rectangular pulses with the readings taken at 10 ns. A single pulse was applied every few seconds. The RC time constant determined by the device capacitance and its resistance was about 3 ns, much shorter than the length of the pulse. The short pulse length Δt ≪ τo ensured the device did not heat up during the pulse as attested by current and voltage remaining constant during the pulse for most data points. Only at large applied voltages did the current increase slightly with a temperature rise calculated to be less than 10 K. As seen in Figure 1c, the quasi-static curves follow the pulsed ones at low voltages but then increase faster with voltage. The divergence is due to Joule heating and the conductivity increasing with temperature. Additionally, by overlaying pulsed I−V characteristics obtained at different stage temperatures onto a single quasi-static I−V curve, it is possible to estimate the increase of device temperature; the method was outlined by Gala et al.34 As an example, the quasi-static I−V curve at the stage temperature of 400 K diverges from the pulsed I−V at approximately 2.5 V. By the end of the quasi-static curve at 4 V, the current coincides with that obtained in a pulsed I−V corresponding to a stage temperature of 435 K (not shown). Thus, we can infer that the device temperature increased by 35 K due to Joule heating before the device reached VTH. It is also apparent that, even in the absence of Joule heating, the current increases superlinearly with voltage. The I(V) dependence of metal/TaOx/metal structures is typically interpreted as due to either a Schottky effect in thin interfacial barriers,35,36 field-enhanced emission from traps (Poole−Frenkel model),23,37−40 or polaronic conductivity.24,41 Most of the work to date indicates the Poole−Frenkel model as best describing the observations. We adopt this model here, taking the functional dependence of conductivity on field and temperature as23,42

σPF(F ,T ) =

σo(T ) 2 +

2 ⎛β F ⎞ ⎛ β F ⎞⎫ σo(T ) ⎛ kBT ⎞ ⎧ − 1⎟exp⎜ ⎟⎬ ⎟ ⎨1 + ⎜ ⎜ F ⎝ β ⎠⎩ ⎝ kBT ⎠ ⎝ kBT ⎠⎭ ⎪







(1)

where ⎛ Nd ⎞2 ⎛ Ed + Et ⎞ σo(T ) = qμNc ⎜ ⎟ exp⎜ − ⎟ ⎝ 2kBT ⎠ ⎝ Nt ⎠

⎛ q3 ⎞1/2 ⎟ β=⎜ ⎝ πεoεi ⎠ (2)

and constants Nd, Ed and Nt, Et correspond to the densities and ionization energies of donors and traps, respectively, Nc is the effective density of states in the conduction band, σo(T) is the low field conductivity, μ is the electron mobility, kB is the Boltzmann constant, q is the elementary charge, εo is the permittivity of free space, εi is the relative dielectric constant of the material, F is the electric field, and T is the stage temperature. The σo(T) in the formula was determined from the experimental data in Figure 1b (energies (Ed + Et)/2 in each sample are equal to the activation energies). An additional parameter, the dielectric constant, was used as an adjustable parameter, and the best fit was obtained for its value of 22. This corresponds to the value of the low frequency constant, in agreement with other reports.32,33,43,44 It should be noted, however, that many publications argue for use of the high frequency constant.35,38,39,45 The experimental data and fitting results are shown in Figure 1d for the 7.5% O2 TaOx film at three different stage temperatures. The fits are in good agreement with experiment for all three types of functional layers (not shown), stage temperatures, and fields. The finite element simulation assumed a 6 μm diameter cylindrical device as part of the circuit consisting of the voltage source and load resistor with the coupled current and heat flow. For electrical conductivity, we have used σPF(F,T) obtained from fitting pulsed I−V’s and extrapolated the experimental data range to higher temperatures as needed. The details of the simulation are given in the Supporting Information. The procedure generated I−V curves that reproduce the experimental quasi-static I−V, including the superlinear behavior in the OFF-state and the NDR region. The results are plotted in Figure 2 for three different stage temperatures and the devices with the functional layer deposited with 7.5% O2. The blue, red, and green lines represent the experimental data, while the black dotted lines correspond to the simulation. The same data are plotted in the inset on a logarithmic scale to show the I−V’s at low current values. The simulation reproduced the exper-

Figure 2. Experimental I−V characteristics of TaOx based devices with a film deposited with 7.5% O2 in the chamber (thick color lines) and finite element simulation of the same device (dotted lines) as a function of voltage and stage temperature. C

DOI: 10.1021/acsami.6b16559 ACS Appl. Mater. Interfaces XXXX, XXX, XXX−XXX

Research Article

ACS Applied Materials & Interfaces

resistance to 0.5 K/μW. Also worth noticing is that the current density at the lateral edge is an order of magnitude lower compared to the current density at the same location at VTH. This drop is due to the lower voltage across the device at point “2” than at VTH (as apparent in Figure 2). The temperature at the center of the filament at this point is 1200 K. With a further increase of total current, the calculated temperature is increasing linearly with current with the fwhm increasing only slightly. The simulation is not reliable much beyond point “2”. The temperature is over 700 K past the highest stage temperature used in the calibration shown in Figure 1, and the Poole−Frenkel formula likely deviates from the real values of electrical conductivity. Also, the thermal properties, which are assumed independent of temperature, could have changed significantly. At this point, we would like to stress that the observed current constriction to form the conducting filament is purely caused by local Joule heating of the functional layer. The device characteristics have not changed as the result of testing, indicating that no physical changes occurred in the structure. If changes such as crystallization of the initially amorphous material have occurred,46 they have not affected the electrical characteristics. This is different than in similar memory-type cells which undergo ion redistribution.47,48 Permanent physical changes do occur in these structures at current densities and time-at-voltage just slightly above those used in this work. These were intentionally avoided during testing.

imental data quite well: the curve in the OFF-state at all temperatures overlaps with and/or is parallel to the experimental curve, the threshold voltage is within 5% of the experimental value, and the NDR regions coincide. Similar agreement was obtained for devices with different oxygen flow rates. More insight into the characteristics of threshold switching can be obtained from evolution of temperature and current density with voltage. The 3-D temperature distribution maps for the two points marked “1” and “2” in Figure 2 for a device with a 7.5% O2 layer at a stage temperature of 300 K are shown in Figure 3a. Figure 3b shows the corresponding line profiles of

4. DISCUSSION For all biasing conditions and stage temperatures, the electrical conductivity of the functional layer is much lower than the conductivity of the metallic electrodes, resulting in a laterally uniform voltage across the electrodes. This results in a uniform current flow and uniform heating at low dissipated power. Upon approaching VTH, the temperature in the center of the device increases and the functional layer becomes a little warmer and more conductive than the lateral boundaries of the device. The current and temperature form a positive feedback loop where increases in current cause increases in temperature, and this, in turn, increases current further. The conditions for the runaway can be estimated in a simple model in which the electrical resistivity as a function of temperature is approximated by the linear term in the Taylor series

Figure 3. (a) Temperature and (b) current density distribution in the 7.5% O2 device corresponding to points “1” (at VTH, continuous blue line) and “2” (within NDR region, dotted red line) in Figure 2.

current density along the radius of the device. Point “1” corresponds to the knee of current at VTH. The average temperature of the functional layer increased by about 50 K above the stage temperature of 300 K with the center of the device being 30 K higher than the lateral edges. The corresponding current density ratio is a factor of 2. The variations across the device at this point are due to lateral heat flow. At point “2” in the NDR region, the current density and temperature distribution are highly nonuniform even though the voltage across the device is the same everywhere. The current density in the center of the device is about 4 orders of magnitude higher than the current density at the lateral edge with the full width at half-maximum (fwhm) of 250 nm. It is interesting to note that the thermal resistances of the OFF-state and NDR regions are different from each other. In the OFFstate, the current and power dissipation is more or less uniform across the diameter of the device, corresponding to a thermal resistance of 0.05 K/μW. In the NDR region, power dissipation is confined to a much smaller area, increasing the thermal

⎛ ⎞ 1 ∂σ V = IRo⎜1 − IVRTH + ···⎟ σ(To) ∂T ⎝ ⎠

(3)

where Ro is the resistance at the stage temperature, σ is the electrical conductivity, and RTH is the thermal resistance between the functional layer and the heat sink. The runaway will occur when ∂V/∂I is zero, leading to IV =

σ(To) ⎛ ∂σ ⎞−1 ⎜ ⎟ 2RTH ⎝ ∂T ⎠

(4)

This gives the value of resistance at runaway R(T) = Ro/2. Thus, runaway will occur when the Joule heating causes the resistance to drop by a factor of 2 and can occur at very different temperatures of the filament depending on the stage temperature. It is easy to notice that the quasi-static current at voltages approaching VTH is, in fact, about twice that of the current in the pulsed I−V’s (Figure 1c). D

DOI: 10.1021/acsami.6b16559 ACS Appl. Mater. Interfaces XXXX, XXX, XXX−XXX

Research Article

ACS Applied Materials & Interfaces

and load resistor, and heat and current flow. The current density as a function of electric field and temperature was obtained from short pulse I−V measurements. It was shown that the threshold switching in TaOx could be well accounted for by a thermal runaway process leading to the formation of a small-diameter hot filament. The thermal runaway is shown to be a direct consequence of the strong dependence of conductivity on temperature.

There are two possible scenarios for the end of a runaway. If the circuit includes a small load resistor, the voltage across the device will remain approximately constant during the runaway. The filament will expand and its temperature will keep increasing until the device is destroyed. On the other hand, if the load resistance is significantly larger than the resistance of the filament, the voltage drop across the device will decrease, limiting the total dissipated power. The device will come to the steady state with a small-diameter hot filament. Figure 4a compares the experimental and simulated values of the threshold voltage for three different functional layers as a



ASSOCIATED CONTENT

S Supporting Information *

The Supporting Information is available free of charge on the ACS Publications website at DOI: 10.1021/acsami.6b16559. Device structure and material parameters and justification used in COMSOL simulation of experimental I−V curves (PDF)



AUTHOR INFORMATION

Corresponding Authors

*E-mail: [email protected] (J.M.G.). *E-mail: [email protected] (M.S.). ORCID

Jonathan M. Goodwill: 0000-0002-3466-3350 Marek Skowronski: 0000-0002-2087-0068 Author Contributions

The manuscript was written through contributions of all authors. All authors have given approval to the final version of the manuscript. Notes

The authors declare no competing financial interest.



Figure 4. Comparison between the experimentally measured (solid lines) and simulated (dashed lines) threshold voltage (a) and threshold power (b) for devices with three different O2 flows.

ACKNOWLEDGMENTS This work was supported in part by FAME, one of six centers of STARnet, a Semiconductor Research Corporation program sponsored by MARCO and DARPA and by NSF Grant DMR 1409068.

function of stage temperature. There are two apparent trends: VTH is decreasing with increasing oxygen deficiency and stage temperature. It should be kept in mind that the conductivities of all devices are increasing with stage temperature and the current at VTH is increasing even though VTH is decreasing. The power at the beginning of the runaway (Figure 4b) is increasing with the stage temperature for both experiment and simulation. The decrease of VTH at lower O2 content in the ambient is associated with the decreasing activation energy of the conductivity. It results in a lower ∂σ/∂T, which means that a greater change of temperature is necessary to halve the initial resistance and, hence, higher dissipated power is required for the runaway. At low stage temperatures and oxygen deficiencies, the only way to increase dissipated power is to increase the VTH. A similar decrease of VTH with an increase of stage temperature was reported in our earlier work on TaOx34 and other oxides.22,23,49 Lastly, it is important to notice that the predicted size of the hot filament, while small, is much bigger than the size of the access devices to be used in memory arrays. The 10 nm-size devices retain the general shape of their I−V characteristics (simulated), i.e., exhibit NDR and threshold switching. The scaling behavior is somewhat complex and will be described in a separate publication.



REFERENCES

(1) Zhou, J.; Kim, K.-H.; Lu, W. Crossbar RRAM Arrays: Selector Device Requirements During Read Operation. IEEE Trans. Electron Devices 2014, 61 (5), 1369−1376. (2) Zhang, L.; Cosemans, S.; Wouters, D. J.; Groeseneken, G.; Jurczak, M.; Govoreanu, B. Selector Design Considerations and Requirements for 1S1R RRAM Crossbar Array. In IEEE 6th International Memory Workshop; IEEE Computer Society, 2014. (3) Kim, S.; Zhou, J.; Lu, W. D. Crossbar RRAM Arrays: Selector Device Requirements During Write Operation. IEEE Trans. Electron Devices 2014, 61 (8), 2820−2826. (4) Jo, S. H.; Kumar, T.; Narayanan, S.; Lu, W. D.; Nazarian, H. 3DStackable Crossbar Resistive Memory Based on Field Assisted Superlinear Threshold (FAST) Selector. In IEEE International Electron Devices Meeting; IEEE, 2014. (5) Song, J.; Woo, J.; Prakash, A.; Lee, D.; Hwang, H. Threshold Selector With High Selectivity and Steep Slope for Cross-Point Memory Array. IEEE Electron Device Lett. 2015, 36 (7), 681−683. (6) Ridley, B. K. Specific Negative Resistance in Solids. Proc. Phys. Soc., London 1963, 82 (6), 954−966. (7) Adler, D.; Henisch, H. K.; Mott, N. The Mechanism of Threshold Switching in Amorphous Alloys. Rev. Mod. Phys. 1978, 50 (2), 209− 220. (8) Adler, D.; Shur, M. S.; Silver, M.; Ovshinsky, S. R. Threshold Switching in Chalcogenide-Glass Thin Films. J. Appl. Phys. 1980, 51 (6), 3289−3309.

5. CONCLUSIONS The quasi-static I−V characteristics of different compositions of TaOx were experimentally measured and simulated by a finite element model that included power distribution between device E

DOI: 10.1021/acsami.6b16559 ACS Appl. Mater. Interfaces XXXX, XXX, XXX−XXX

Research Article

ACS Applied Materials & Interfaces

Process in TiO2 Based Resistive Switching Devices. Appl. Phys. Lett. 2013, 102 (2), 023507. (29) Noman, M.; Sharma, A. A.; Lu, Y. M.; Kamaladasa, R.; Skowronski, M.; Salvador, P. A.; Bain, J. A. Mechanism of Localized Electrical Conduction at the Onset of Electroforming in TiO2 Based Resistive Switching Devices. Appl. Phys. Lett. 2014, 104 (11), 113510. (30) Sharma, A. A.; Noman, M.; Skowronski, M.; Bain, J. A. HighSpeed In-Situ Pulsed Thermometry in Oxide RRAMs. In Proceedings of the Technical Program - International Symposium on VLSI Technology, Systems and Application (VLSI-TSA); IEEE, 2014. (31) Lau, W. S.; Tan, T. S.; Sandler, N. P.; Page, B. S. Characterization of Defect States Responsible for Leakage Current in Tantalum Pentoxide Films for Very-High-Density Dynamic Random Access Memory (DRAM) Applications. Jpn. J. Appl. Phys. 1995, 34 (Pt. 1, No. 2B), 757−761. (32) Wu, X. M.; Soss, S. R.; Rymaszewski, E. J.; Lu, T.-M. Dielectric Constant Dependence of Poole-Frenkel Potential in Tantalum Oxide Thin Films. Mater. Chem. Phys. 1994, 38 (3), 297−300. (33) Ezhilvalavan, S.; Tseng, T.-Y. Conduction Mechanisms in Amorphous and Crystalline Ta2O5 Thin Films. J. Appl. Phys. 1998, 83 (9), 4797−4801. (34) Gala, D. K.; Sharma, A. A.; Li, D.; Goodwill, J. M.; Bain, J. A.; Skowronski, M. Low Temperature Electroformation of TaOx-Based Resistive Switching Devices. APL Mater. 2016, 4 (1), 016101. (35) Blonkowski, S.; Regache, M.; Halimaoui, A. Investigation and Modeling of the Electrical Properties of Metal-Oxide-Metal Structures Formed from Chemical Vapor Deposited Ta2O5 Films. J. Appl. Phys. 2001, 90 (3), 1501−1508. (36) Zeng, W.; Eisenbraun, E.; Frisch, H.; Sullivan, J. J.; Kaloyeros, A. E.; Margalit, J.; Beck, K. CVD of Tantalum Oxide Dielectric Thin Films for Nanoscale Device Applications. J. Electrochem. Soc. 2004, 151 (8), F172−F177. (37) Young, P. L. DC Electrical Conduction in Thin Ta2O5 Films. I. Bulk-Limited Conduction. J. Appl. Phys. 1976, 47 (1), 235−241. (38) Oehrlein, G. S. Oxidation Temperature Dependence of the Dc Electrical Conduction Characteristics and Dielectric Strength of Thin Ta2O5 Films on Silicon. J. Appl. Phys. 1986, 59 (5), 1587−1595. (39) Choi, W. K.; Ling, C. H. Analysis of the Variation in the FieldDependent Behavior of Thermally Oxidized Tantalum Oxide Films. J. Appl. Phys. 1994, 75 (8), 3987−3990. (40) Spassov, D.; Atanassova, E.; Virovska, D. Electrical Characteristics of Ta2O5 Based Capacitors with Different Gate Electrodes. Appl. Phys. A: Mater. Sci. Process. 2006, 82 (1), 55−62. (41) Bryksin, V. V.; Goltsev, A. V.; Khanin, S. D.; Novotelnova, A. V.; Vasilev, A. N. Nonlinear Current-Voltage Characteristics of Ta2O5 and Nb2O5 Amorphous Oxides. Phys. Status Solidi B 1990, 161 (2), 777− 781. (42) Hartke, J. L. The Three-Dimensional Poole-Frenkel Effect. J. Appl. Phys. 1968, 39 (10), 4871−4873. (43) Hughes, D. M.; Jones, M. W. Electrical Conduction in Reactively Sputtered Tantalum Oxide Thin Films. J. Phys. D: Appl. Phys. 1974, 7, 2081−2096. (44) Chiu, F.; Wang, J.; Lee, J. Y.; Wu, S. C. Leakage Currents in Amorphous Ta2O5 Thin Films. J. Appl. Phys. 1997, 81 (10), 6911− 6915. (45) Oehrlein, G. S.; Reisman, A. Electrical Properties of Amorphous Tantalum Pentoxide Thin Films on Silicon. J. Appl. Phys. 1983, 54 (11), 6502−6508. (46) Kwon, J.; Sharma, A. A.; Chen, C. Y.; Fantini, A.; Jurczak, M.; Herzing, A. A.; Bain, J. A.; Picard, Y. N.; Skowronski, M. Transient Thermometry and High-Resolution Transmission Electron Microscopy Analysis of Filamentary Resistive Switches. ACS Appl. Mater. Interfaces 2016, 8 (31), 20176−20184. (47) Wong, H.-S. P.; Lee, H.-Y.; Yu, S.; Chen, Y.-S.; Wu, Y.; Chen, P.-S.; Lee, B.; Chen, F. T.; Tsai, M.-J. Metal-Oxide RRAM. Proc. IEEE 2012, 100, 1951−1970. (48) Kim, K. M.; Jeong, D. S.; Hwang, C. S. Nanofilamentary Resistive Switching in Binary Oxide System; a Review on the Present Status and Outlook. Nanotechnology 2011, 22 (25), 254002.

(9) Ielmini, D. Threshold Switching Mechanism by High-Field Energy Gain in the Hopping Transport of Chalcogenide Glasses. Phys. Rev. B: Condens. Matter Mater. Phys. 2008, 78 (3), 035308. (10) Karpov, V. G.; Kryukov, Y. A.; Karpov, I. V.; Mitra, M. FieldInduced Nucleation in Phase Change Memory. Phys. Rev. B: Condens. Matter Mater. Phys. 2008, 78 (5), 052201. (11) Tsendin, K. Electro-Thermal Theory of the Switching and Memory Effects in Chalcogenide Glassy Semiconductors. Phys. Status Solidi B 2009, 246 (8), 1831−1836. (12) Le Gallo, M.; Athmanathan, A.; Krebs, D.; Sebastian, A. Evidence for Thermally Assisted Threshold Switching Behavior in Nanoscale Phase-Change Memory Cells. J. Appl. Phys. 2016, 119 (2), 025704. (13) Stefanovich, G.; Pergament, A.; Stefanovich, D. Electrical Switching and Mott Transition in VO2. J. Phys.: Condens. Matter 2000, 12 (41), 8837−8845. (14) Zimmers, A.; Aigouy, L.; Mortier, M.; Sharoni, A.; Wang, S.; West, K. G.; Ramirez, J. G.; Schuller, I. K. Role of Thermal Heating on the Voltage Induced Insulator-Metal Transition in VO2. Phys. Rev. Lett. 2013, 110 (5), 056601. (15) Pickett, M. D.; Stanley Williams, R. Sub-100 fJ and SubNanosecond Thermally Driven Threshold Switching in Niobium Oxide Crosspoint Nanodevices. Nanotechnology 2012, 23 (21), 215202. (16) Zhou, Y.; Chen, X.; Ko, C.; Yang, Z.; Mouli, C.; Ramanathan, S. Voltage-Triggered Ultrafast Phase Transition in Vanadium Dioxide Switches. IEEE Electron Device Lett. 2013, 34 (2), 220−222. (17) Li, D.; Sharma, A. A.; Gala, D. K.; Shukla, N.; Paik, H.; Datta, S.; Schlom, D. G.; Bain, J. A.; Skowronski, M. Joule Heating-Induced Metal − Insulator Transition in Epitaxial VO2/TiO2 Devices. ACS Appl. Mater. Interfaces 2016, 8 (20), 12908. (18) Brockman, J. S.; Gao, L.; Hughes, B.; Rettner, C. T.; Samant, M. G.; Roche, K. P.; Parkin, S. S. P. Subnanosecond Incubation Times for Electric-Field-Induced Metallization of a Correlated Electron Oxide. Nat. Nanotechnol. 2014, 9 (6), 453−458. (19) Shin, J.; Kim, I.; Biju, K. P.; Jo, M.; Park, J.; Lee, J.; Jung, S.; Lee, W.; Kim, S.; Park, S.; Hwang, H. TiO2-Based Metal-Insulator-Metal Selection Device for Bipolar Resistive Random Access Memory CrossPoint Application. J. Appl. Phys. 2011, 109 (3), 033712. (20) Sharma, A. A.; Noman, M.; Abdelmoula, M.; Skowronski, M.; Bain, J. A. Electronic Instabilities Leading to Electroformation of Binary Metal Oxide-Based Resistive Switches. Adv. Funct. Mater. 2014, 24 (35), 5522−5529. (21) Sharma, A. A.; Karpov, I. V.; Kotlyar, R.; Kwon, J.; Skowronski, M.; Bain, J. A. Dynamics of Electroforming in Binary Metal OxideBased Resistive Switching Memory. J. Appl. Phys. 2015, 118 (11), 114903. (22) Slesazeck, S.; Mähne, H.; Wylezich, H.; Wachowiak, A.; Radhakrishnan, J.; Ascoli, A.; Tetzlaff, R.; Mikolajick, T. Physical Model of Threshold Switching in NbO2 Based Memristors. RSC Adv. 2015, 5 (124), 102318−102322. (23) Gibson, G. A.; Musunuru, S.; Zhang, J.; Vandenberghe, K.; Lee, J.; Hsieh, C.-C.; Jackson, W.; Jeon, Y.; Henze, D.; Li, Z.; Stanley Williams, R. An Accurate Locally Active Memristor Model for S-Type Negative Differential Resistance in NbOx. Appl. Phys. Lett. 2016, 108 (2), 023505. (24) Funck, C.; Menzel, S.; Aslam, N.; Zhang, H.; Hardtdegen, A.; Waser, R.; Hoffmann-Eifert, S. Multidimensional Simulation of Threshold Switching in NbO2 Based on an Electric Field Triggered Thermal Runaway Model. Adv. Electron. Mater. 2016, 2 (7), 1600169. (25) Janninck, R. F.; Whitmore, D. H. ELECTRICAL CONDUCTIVITY AND THERMOELETRIC POWER OF NIOBIUM DIOXIDE. Solid State Commun. 1966, 4 (3), xxv−xxvi. (26) Sakata, T.; Sakata, K.; Nishida, I. Study of Phase Transition in NbO2. Phys. Status Solidi B 1967, 20 (2), K155−K157. (27) Strickland, J. A. Time-Domain Reflectometry Measurements; Tektronix: Beaverton, OR, 1973. (28) Noman, M.; Sharma, A. A.; Lu, Y. M.; Skowronski, M.; Salvador, P. A.; Bain, J. A. Transient Characterization of the Electroforming F

DOI: 10.1021/acsami.6b16559 ACS Appl. Mater. Interfaces XXXX, XXX, XXX−XXX

Research Article

ACS Applied Materials & Interfaces (49) Chudnovskii, F. a.; Odynets, L. L.; Pergament, A. L.; Stefanovich, G. B. Electroforming and Switching in Oxides of Transition Metals: The Role of Metal−Insulator Transition in the Switching Mechanism. J. Solid State Chem. 1996, 122 (1), 95−99.

G

DOI: 10.1021/acsami.6b16559 ACS Appl. Mater. Interfaces XXXX, XXX, XXX−XXX