Hafnium oxide non-volatile memories have shown promise as an artificial synapse in neuromorphic computing architectures. However, there is still a need to fundamentally understand how to reliably control the analog resistance change induced by oxygen ions that partially rupture or re-form the conductive filament. In this work, the impact of measurement conditions (pulse amplitude and pulse width) and titanium dopants on the analog resistance change of atomic layer deposited hafnium oxide memristor synapses are studied. A lower pulse amplitude improves the linearity of resistance change as a function of the number of pulses but results in a smaller memory window. The addition of titanium dopants does not substantively change the analog resistance modulation of hafnium oxide. Density functional theory calculations show that titanium strongly impacts oxygen ion motion in the HfxTiyOz matrix but does not impact significantly in the HfTi metallic filament. This study demonstrates that the analog characteristic of HfxTiyOz artificial synapses is largely independent of the titanium doped bulk oxide since the resistance change is primarily controlled by the HfTi metallic conducting filament.

As traditional computing has been facing limits associated with the end of conventional device scaling,1,2 a new computing paradigm is desired. Neuromorphic computing3–13 has gained widespread attention to circumvent the challenges of traditional computing by mimicking information processing associated with the brain enabling computational tasks such as vector-matrix multiplication, pattern recognition, and adaptation in the real-world environment.14–17 Memristors7,18–22 are considered a leading candidate for artificial synapses in neuromorphic computing. Memristors based on transition metal oxides (e.g., HfO2,23–26 TaOx,27,28 TiO2,29,30 WO3,31–33 SrTiO3,33–35 ZnO,36,37 and SiO238,39) have gained widespread attention because of their multiple attractive features such as high scalability, multi-level capability, low energy consumption, and compatibility with large-scale integration. Among these, HfOx exhibits additional advantages such as CMOS compatibility, scalability, lower power consumption, large memory window (>100), high density, fast-switching speed (∼ns), and good endurance (≥1010 cycles).3,40–42 Beyond these features of digital Resistive Random Access Memory (RRAM), promising analog synaptic characteristics have also been reported for HfOx memristors. The change of the weight of biological synapses when pulsed with an action potential is associated with the motion of neurotransmitters.43–46 Similarly, the resistance change in filamentary oxide memristors is caused by the motion of oxygen vacancies (Vo) or ions across the conducting filament. A variety of biological synaptic responses such as short-term depression/potentiation (STD/STP),47 long-term depression/potentiation (LTD/LTP),48,49 and spike-time-dependent plasticity (STDP)48,50 have also been demonstrated for HfOx artificial synapses for analog applications.

Pattern recognition accuracy of analog neuromorphic architectures has been shown to be dependent on the control and linearity of the resistance change of the artificial memristor synapse as a function of applied voltage pulses.51 However, a linear change in resistance as a function of identical voltage pulses is not easily achieved for HfOx memristor synapses.52 Non-identical voltage pulses and transistors can be used for HfOx memristor synapses to improve the linearity. However, this requires complex peripheral circuits for pulse generation and increases the burden of latency, energy, and chip area.51,53–55

It is known that the motion of a small number of ions is involved in the analog resistive switching of HfOx memristors.56 Therefore, the incorporation of layers of other oxides with HfOx23,57 as well as the addition of dopants to HfOx58,59 have been considered to alter ion motion and improve resistive switching properties. For example, a previous study reported that adding TiOx layers to sputtered HfOx improves the uniformity in digital switching and provides multi-level capability.23,60 Another study showed that sputtered HfOx/TiOx multilayered devices consume less energy compared to HfOx devices and can be integrated in the large-scale neuromorphic visual system.40 Nickel doping in sputtered HfOx has been shown to decrease power consumption and improve multi-level resistance behavior with identical pulses.61 Synaptic behavior such as long-term potentiation and depression has been reported for a memristor with sputtered HfOx doped with manganese62 yet the switching window was small and an abrupt increase of the conductance with applied pulse voltage was observed. While sputtering can be useful for demonstration purposes, atomic layer deposition (ALD) is more compatible with large-scale integrated circuit fabrication.63 A previous study has reported using ALD deposited electro thermal modulation layer to improve the linearity of analog HfOx RRAM.64 However, a compliance transistor is used during measurement, and the impact of measurement conditions on the observed behavior was not studied. Aluminum doping of HfOx devices has been reported to improve analog switching; however, a complex pulsing scheme was used.58,65 Chakrabarti et al.66,67 and a recent study68 have explored the use of titanium doped HfOx deposited using ALD for multi-level memory; however, the impact of these dopants on the analog switching behavior was not evaluated.

This study investigates the analog switching properties of ALD HfOx and titanium doped HfOx memristors using identical voltage pulses. In addition, the impact of titanium doping on the analog resistance change is analyzed for different measurement conditions such as pulse amplitude and width. The impact of these measurement conditions on analog features such as off state resistance and linearity is assessed. While the measurement characteristics are observed to strongly impact the analog characteristics, the results are largely independent of titanium doping. Density functional theory (DFT) calculation of oxygen ion motion in these materials suggests that this is due to the resistance change being primarily controlled by the HfTi metallic conducting filament and not the oxide itself.

Metal-insulator-metal memristors with HfOx and HfxTiyOz amorphous oxides were fabricated as shown in Fig. 1(a). Silicon wafers with 310 nm of SiO2 were cleaned with acetone, methanol, and isopropanol. Bottom electrodes (BE) with size 10 × 10 μm2 were formed using mask-less ultraviolet photolithography followed by lift off. The BEs consisted of a ∼5 nm titanium adhesive layer and ∼70 nm of gold deposited using electron beam evaporation at a rate of 0.1 nm/s at a pressure 2.5 × 10−6 Torr. For the standard device, the ∼5 nm HfOx active layer was synthesized using thermal ALD at 250 °C with de-ionized water and tetrakis(dimethylamido) hafnium (TDMAHF) as precursors. For the titanium doped devices, HfxTiyOz was synthesized using thermal ALD by sequential deposition of three cycles of tetrakis(dimethylamido)hafnium (TDMAHF) and de-ionized water followed by one cycle of tetrakis(dimethylamido)titanium (TDMAT) and de-ionized water. This was repeated 12 times for 48 total cycles to achieve a thickness ∼5 nm HfxTiyOz oxide. X-ray photoelectron spectroscopy analysis indicates a Hf:Ti ratio of 3.14:1 (Fig. S2 in the supplementary material) which is close to the estimated 3:1 ratio from the synthesis conditions. Spectroscopic ellipsometry was utilized to verify the thicknesses of the oxides. The ∼5 nm titanium capping layer and ∼150 nm gold top electrode were deposited using electron beam evaporation and patterned using lift off. Before the deposition of the TE, BE, and oxide, an oxygen plasma descum process for 30 s (at a flow-rate of 50 SCCM, plasma power of 150 W RF, and pressure of 60 mT) was performed to remove the residual photoresist and surface hydrocarbons.

FIG. 1.

Schematics of the fabricated (M–I–M) device structure having ∼5 nm active layer of HfOx or HfxTiyOz. Oxygen vacancy motion in the active layer oxide between the top and bottom electrodes in the synaptic device can imitate the biological synaptic dynamics.

FIG. 1.

Schematics of the fabricated (M–I–M) device structure having ∼5 nm active layer of HfOx or HfxTiyOz. Oxygen vacancy motion in the active layer oxide between the top and bottom electrodes in the synaptic device can imitate the biological synaptic dynamics.

Close modal

A Keithley 4200 SCS semiconductor parameter analyzer was used to conduct the electrical testing. Identical square pulses of a given pulse width and amplitude were applied. The resistance of the device was measured between each pulse by performing a current–voltage sweep from 0 to 0.1 V. These measurements were repeated using different pulse amplitudes (−0.1 to −0.8 V) and pulse widths (160 ns to 1 ms), on both HfOx and HfxTiyOz devices. Each measurement condition was repeated on multiple devices. Figure S1 in the supplementary material illustrates the conditions used for the measurements.

We employed the density functional theory (DFT) provided by VASP to calculate the diffusion barrier of neutral oxygen vacancy. For each case, we obtained the optimized structure of a supercell cell with one oxygen vacancy (HfxTiyOz) or one interstitial oxygen (HfTi metal) first. Then, the nudged elastic band method combined with the climbing image method69 was used to calculate the neutral oxygen vacancy diffusion barrier. For HfO2, the structure optimization was performed on a supercell of 2 × 2 × 2 with one oxygen vacancy. For the titanium doped case, the effect is studied by replacing Hf with Ti performing the DFT with neutral oxygen vacancies in both oxides and metals. For Hf1.5Ti0.5O2, we first performed the structure optimization in the unit cell with the ratio of Hf:Ti being 3:1, and then, we constructed the supercell of 2 × 2 × 2 using the unit cell obtained above. A similar method was used to obtain the structure of HfTi metal with a supercell of 4 × 4 × 3.

The memristor devices were formed and stabilized prior to performing pulse measurements. The conductive filament was formed by applying a positive voltage sweep at the top electrode from 0 to 3.5 V with a current compliance ∼0.1 mA (Fig. S3 in the supplementary material). The forming voltage is observed to be lower for the HfxTiyOz devices (∼2.5 V) compared to the HfOx devices (∼3.1 V). The forming step was followed by a gradual reset from 0 to a maximum negative reset voltage with an increment of ∼−0.1 V and several set–reset stabilization loops (Figs. S4 and S5 in the supplementary material). This stabilization is performed to ensure the devices are toggling between the HRS and the LRS and that the conducting filaments are stable. After forming and stabilization, the resistances of the devices were measured as shown in Fig. S6 in the supplementary material. It is observed that after stabilization the resistance states were nearly identical for both HfOx and HfxTiyOz devices, confirming that any difference in the analog temporal response is independent of the initial device state.

Figures 2(a) and 2(b) show the change in resistance under identical negative pulses for the HfOx and HfxTiyOz devices, respectively. As mentioned previously, the application of identical pulses is advantageous because it dramatically simplifies the pulse generation module design.40,41 In general, the resistance is observed to increase dramatically upon the first pulse and then eventually saturates. This highly nonlinear “first pulse effect” is typically observed in HfOx based devices,62 and it degrades the performance of analog neuromorphic architectures.70 The large resistance change with the first pulse is less dramatic with decreasing pulse amplitude and pulse width. The maximum temperature of the filament is known to be strongly dependent on the current and, hence, voltage. For the first negative voltage pulse, the temperature initially increases causing oxygen ions to move from the titanium capping layer into the filament. The oxidation of the filament increases the resistance, thus quickly reducing the current and the temperature. Subsequent pulses at the same negative voltage, therefore, result in a lower temperature and a slower change in resistance as compared to the first pulse. A smaller pulse width generally results in a smaller initial change in resistance and temperature (Fig. S7 in the supplementary material) because the time for oxygen ion diffusion at peak temperature is reduced.

FIG. 2.

Measured resistance as a function of pulse number with variable pulse amplitudes and widths for devices with (a) HfOx and (b) HfxTiyOz devices. The data for each condition are an average from ten separate device measurements. Open and closed symbols indicate 1 ms and 160 ns pulse widths, respectively. An abrupt increase of resistance at the first pulse is observed at −0.8 V for both HfOx and HfxTiyOz devices. At a smaller pulse amplitude, this first pulse effect is reduced. It is observed that a smaller pulse width (160 ns) also decreases the first pulse effect for both HfOx and HfxTiyOz devices.

FIG. 2.

Measured resistance as a function of pulse number with variable pulse amplitudes and widths for devices with (a) HfOx and (b) HfxTiyOz devices. The data for each condition are an average from ten separate device measurements. Open and closed symbols indicate 1 ms and 160 ns pulse widths, respectively. An abrupt increase of resistance at the first pulse is observed at −0.8 V for both HfOx and HfxTiyOz devices. At a smaller pulse amplitude, this first pulse effect is reduced. It is observed that a smaller pulse width (160 ns) also decreases the first pulse effect for both HfOx and HfxTiyOz devices.

Close modal

As shown previously, the resistance of the devices saturate after multiple pulses (∼200) to a value that we term as off-resistance (Roff), which depends on the pulsing conditions. Figures 3(a) and 3(b) show Roff as a function of pulse amplitude, in the range from −0.5 to −0.8 V with two different pulse widths (160 ns and 1 ms). It is observed that Roff is strongly dependent on the pulse voltage while being almost independent of the pulse width. This trend is observed in both standard HfOx and HfxTiyOz devices. The strong effect of negative pulse voltage on Roff can be attributed to the fact that voltage primarily controls the reset behavior of the bipolar memristors.71,72 The impact of pulse width on Roff is not significant compared to the pulse amplitude because total energy input has a second-order relationship with the voltage and a linear relationship with time.73 Moreover, there is a voltage below which no significant change in Roff is observed because the energy input to move oxygen ions is too low to cause any change in resistance. This voltage at which a resistance change is measured is similar (∼ − 0.6 V) for both the HfOx and HfxTiyOz devices. The observed results indicate that by adjusting the pulse amplitude, the saturated resistance of the device can be increased, which can be applied to neuromorphic applications.

FIG. 3.

Saturated resistance (Roff) dependency on pulse amplitude and pulse width. (a) HfOx device and (b) HfxTiyOz devices. Roff indicates the maximum resistance that can be observed at a particular measurement condition. Below ∼−0.6 V, no change in Roff is observed for HfOx and HfxTiyOz devices. It is observed that the saturated resistance (Roff) is mostly dominated by the pulse amplitude. At any measurement condition, the maximum and minimum resistances that can be achieved are denoted by Rmax and Rmin, respectively, indicated as shaded (purple) regions.

FIG. 3.

Saturated resistance (Roff) dependency on pulse amplitude and pulse width. (a) HfOx device and (b) HfxTiyOz devices. Roff indicates the maximum resistance that can be observed at a particular measurement condition. Below ∼−0.6 V, no change in Roff is observed for HfOx and HfxTiyOz devices. It is observed that the saturated resistance (Roff) is mostly dominated by the pulse amplitude. At any measurement condition, the maximum and minimum resistances that can be achieved are denoted by Rmax and Rmin, respectively, indicated as shaded (purple) regions.

Close modal

Figures 4(a) and 4(b) illustrate the analog resistance change of HfOx and HfxTiyOz devices under pulse amplitudes of −0.7 and −0.8 V and pulse widths of 1 ms and 160 ns. It can be observed that for the −0.8 V and 160 ns pulsing condition, the titanium doped device (dark red curve-closed symbol) has a lower first pulse effect and improved linearity compared to the HfOx device (light red curve-closed symbol). A parameter α can be introduced here to quantify linearity, which is defined as the smaller the value of α higher the linearity is.74,75 Details of the model describing the α parameter can be found in Eq. (S1) in the supplementary material. While for the HfOx device, the α parameter value is ∼15, for titanium doped device, the value is ∼2, for the −0.8 V and 160 ns pulsing condition, shown in Figs. 4(c) and 4(d). Although a previous study has reported an alpha value of −0.63, the applied voltage (−1.5 V) was higher than in our case (−0.8 V). Therefore, the titanium doped HfOx device has a comparable alpha value and is a low-power device with a comparable alpha value.

FIG. 4.

Titanium dopant impact on the analog resistance change: (a) measured resistance as a function of pulse number for first 20 pulses and (b) measured resistance as a function of pulse number for 100 pulses. Here, the inset shows Rmax and Rmin distribution for 100 pulses, and saturated resistance (Roff) is well below that of Rmax. Open and closed symbols indicate 1 ms and 160 ns pulse widths, respectively. (c) Alpha parameter fitting for HfOx at −0.8 V and 160 ns and (d) alpha parameter fitting for HfxTiyOz at −0.8 V and 160 ns. A reduction of the first pulse effect and an improvement in the linearity are observed for the titanium doped HfOx device. For each experimental condition, eight to ten devices were measured. The standard deviation of the device data is shown in Fig. S9 in the supplementary material.

FIG. 4.

Titanium dopant impact on the analog resistance change: (a) measured resistance as a function of pulse number for first 20 pulses and (b) measured resistance as a function of pulse number for 100 pulses. Here, the inset shows Rmax and Rmin distribution for 100 pulses, and saturated resistance (Roff) is well below that of Rmax. Open and closed symbols indicate 1 ms and 160 ns pulse widths, respectively. (c) Alpha parameter fitting for HfOx at −0.8 V and 160 ns and (d) alpha parameter fitting for HfxTiyOz at −0.8 V and 160 ns. A reduction of the first pulse effect and an improvement in the linearity are observed for the titanium doped HfOx device. For each experimental condition, eight to ten devices were measured. The standard deviation of the device data is shown in Fig. S9 in the supplementary material.

Close modal

Therefore, although the saturated resistance (Roff) is larger for the HfOx device than the titanium doped device, due to the “first pulse effect,” it has a higher value of α and suffers from poor linearity. The addition of titanium dopant appears to lower the “first pulse effect,” improving the linearity. Moreover, the alpha fittings at different pulsing conditions for HfOx and HfxTiyOz devices are shown in Fig. S8 in the supplementary material. It is also that for a given device, the alpha value decreases for the lower pulse amplitude condition. This provides further evidence for the conclusion made earlier, which is lower pulse amplitude improves the linearity. The linearity improvement is a core requirement for the synaptic devices to achieve a high pattern recognition accuracy.41 

While there are subtle differences in the linearity of the HfOx and HfxTiyOz devices, the analog behavior is quite comparable. It has been reported that for the switching of oxide-based devices, oxygen ion inducing redox reactions and mixed electronic–ionic transport in the oxide region play a critical role. As the oxygen ions will diffuse through the pathways of the lowest barrier, the activation energy of oxygen ion migration has a direct impact on the analog switching. Theoretically, the migration barrier in bulk oxide is expected to be impacted significantly by adding titanium as the Gibbs free energy of formation HfOx is −1088.2 eV and TiOx is −888.8 eV. Therefore, a quantitative analysis of the oxygen vacancy migration barrier is needed once the conductive filament is formed for both HfOx and HfxTiyOz devices.

To analyze the effect of titanium doping on the oxygen ion migration, we have employed the density functional theory (DFT) with the nudged elastic band (NEB) method to theoretically calculate the neutral oxygen ion diffusion barrier in both HfOx and HfxTiyOz systems. A point should be noted that once the CF forms, most oxygen vacancies have aggregated in the CF, and the bulk region has a low concentration of oxygen vacancies. In the system of low concentration of oxygen vacancies, the hopping of oxygen ions is equivalent to oxygen vacancies moving forward. Therefore, the diffusion barrier of neutral oxygen ion (vacancy) is relevant in this case. Moreover, an interstitial oxygen ion has very high formation energy (1.69–2.76 eV) and charged oxygen vacancies are likely to be present near the interfaces. As a result, these are not likely true in this case. A detailed description of our calculations for the neutral oxygen ions can be found in Sec. II.

The DFT calculation results of the neutral oxygen ion diffusion barrier for HfOx and HfxTiyOz are illustrated in Fig. 5. It can be observed that the diffusion barrier is ∼four times lower for HfxTiyOz (∼0.6 eV) compared to HfOx (∼2.49 eV) oxide. The calculated barrier energy in HfOx oxide is consistent with the previously published results.76,77 The lower diffusion barrier for titanium doped HfOx oxide is expected since the neutral oxygen ion diffusion barrier in pure TiO2 is much lower than in pure HfO2. A ∼fourfold decrease in the energy barrier of titanium doped HfOx oxide compared to pure HfOx oxide is expected to cause a pronounced change in the analog responses in HfxTiyOz devices.

FIG. 5.

DFT calculation of oxygen vacancy diffusion barrier for HfOx and HfxTiyOz systems. The result shows the diffusion barrier in the HfxTiyOz system is almost ∼four times lower compared to the HfOx system.

FIG. 5.

DFT calculation of oxygen vacancy diffusion barrier for HfOx and HfxTiyOz systems. The result shows the diffusion barrier in the HfxTiyOz system is almost ∼four times lower compared to the HfOx system.

Close modal

However, the observed differences in the experimental results of analog responses for both HfOx and HfxTiyOz devices are not large, suggesting that the bulk oxide does not play the dominant role in the analog response suggesting that the diffusion of oxygen ions within the conducting filament should be explored.

Since there are much fewer oxygen ions inside the CF, the CF primarily consists of either a Hf or HfTi alloys system for HfOx and HfxTiyOz devices, respectively.69,70 As a result, we employ the DFT with NEB to calculate the diffusion barrier of the oxygen ion at interstitial sites in Hf and HfTi alloy metals, and our results are summarized in Figs. 6 and 7. Previous DFT study on the oxygen diffusion in HCP Hf metal has shown that stable oxygen interstitial sites are either octahedral O 2a(0,0,0)) or tetrahedral T 4f(23,13,18) with a site energy difference of ETEO=0.91eV. The diffusion barriers were found to be 1.97 eV in the path from O to T and 1.06 eV in the reverse direction.78 Same definition was followed to perform the NEB for HfTi alloy with Hf:Ti being 3:1 in this study. Four stable oxygen interstitial sites in HfTi alloy metal were found, and diffusion barriers in paths between them were calculated. In each path, three images were calculated using the NEB method to obtain the diffusion barrier height, the difference between the peak in energy and the endpoints (stable interstitial sites). Since the diffusion barrier is defined as the energy difference between peak and the endpoints along the path, site C, the site with lowest site energy from our calculations was chosen as the reference point. Tables I and II summarize the simulation results. It is observed that the minimal diffusion barrier occurs in the path from A to D with a value of 2.02 and 0.84 eV in the reverse direction, which is close to the values found in the pure Hf HCP metal. Our results suggest that Hf and HfTi metallic filaments have similar oxygen ion diffusion barriers.

FIG. 6.

DFT calculation for Hf metal and HfTi alloy metal systems. (a) The oxygen ion diffusion barrier in the Hf system is similar to that of the HfTi alloy system. (b) DFT with the nudged elastic band method (NEB) was used to calculate the diffusion barrier of oxygen ion in sites A, B, C, and D.

FIG. 6.

DFT calculation for Hf metal and HfTi alloy metal systems. (a) The oxygen ion diffusion barrier in the Hf system is similar to that of the HfTi alloy system. (b) DFT with the nudged elastic band method (NEB) was used to calculate the diffusion barrier of oxygen ion in sites A, B, C, and D.

Close modal
FIG. 7.

Schematics showing oxygen ion motion during analog pulse application on (a) HfOx and (b) HfxTiyOz devices. Although there is a difference in oxygen ion migration barriers in bulk oxides, migration barriers are similar both in Hf and HfTi metal-rich conducting filaments. Application of pulses in milliseconds or nanoseconds range primarily moves oxygen ions within the conducting filaments. Therefore, the resulting analog responses are not significantly different.

FIG. 7.

Schematics showing oxygen ion motion during analog pulse application on (a) HfOx and (b) HfxTiyOz devices. Although there is a difference in oxygen ion migration barriers in bulk oxides, migration barriers are similar both in Hf and HfTi metal-rich conducting filaments. Application of pulses in milliseconds or nanoseconds range primarily moves oxygen ions within the conducting filaments. Therefore, the resulting analog responses are not significantly different.

Close modal
TABLE I.

DFT calculation results for diffusion barriers in sites A, B, C, and D for the HfTi alloy metal system.

SiteWyckoff positionΔE (site energy relative to site C) (eV)
Octahedral, 2a(0, 0, 0) (no Ti around) 0.0303  
Hexahedral, 2d(23,13,14) (under Hf) 0.9408 
Octahedral, 2a(0, 0, 0) (with Ti around) 
Tetrahedral, 4f(23,13,18)(under Ti) 1.2156 
SiteWyckoff positionΔE (site energy relative to site C) (eV)
Octahedral, 2a(0, 0, 0) (no Ti around) 0.0303  
Hexahedral, 2d(23,13,14) (under Hf) 0.9408 
Octahedral, 2a(0, 0, 0) (with Ti around) 
Tetrahedral, 4f(23,13,18)(under Ti) 1.2156 
TABLE II.

DFT calculation results for diffusion barriers of possible paths for the HfTi alloy metal system.

PathDiffusion barrier (eV)Diffusion barrier (reverse) (eV)
A to B 2.0459  1.1353 
C to B 2.0761 1.1353 
C to A 2.5174 2.4871 
A to D 2.0253 0.8400 
Octahedral to tetrahedral (pure Hf) 2.0530 1.0610 
PathDiffusion barrier (eV)Diffusion barrier (reverse) (eV)
A to B 2.0459  1.1353 
C to B 2.0761 1.1353 
C to A 2.5174 2.4871 
A to D 2.0253 0.8400 
Octahedral to tetrahedral (pure Hf) 2.0530 1.0610 

Similar diffusion barriers of oxygen ion in both Hf metal and HfTi alloy metal systems indicate that the oxygen ion diffusion dynamics are same in both metallic systems as shown in Fig 7. Moreover, it has been reported that in HfOx memristors, digital switching involves movement of a few atoms within a fast-switching time scale. During analog switching, pulses are applied in even shorter timescales such as pulse widths of milliseconds or nanoseconds range. Therefore, very few numbers of oxygen ions are involved during the analog temporal response. If the oxygen ions primarily move through the Hf or HfTi metal-rich CFs having similar diffusion barriers, only then the analog responses will be comparable.

In summary, HfOx based memristors having a ∼5 nm HfOx and HfxTiyOz amorphous oxide were fabricated. The analog switching of the fabricated devices were analyzed under identical pulsing conditions. The impact of measurement conditions on the key analog features such as linearity and off-resistance were assessed. The results suggest that although a small pulse width helps decrease the abrupt resistance increase at the first pulse, pulse amplitude has a pronounced impact on the resistance change. A smaller pulse amplitude decreases the first pulse effect and improves the linearity. This is likely because a small pulse amplitude causes a smaller change in the initial filament temperature. However, a larger pulse amplitude is also necessary to achieve a higher saturated resistance (Roff) and a larger memory window. Therefore, there is a trade-off between the saturated resistance and linearity. Moreover, although an improvement in the linearity is observed for the titanium doped HfOx device, analog responses are comparable to the HfOx device. The DFT calculation verifies that the titanium addition in HfOx decreases the oxygen ion migration barrier in the bulk oxide by almost ∼four times; however, the migration barriers are similar in the HfTi and Hf metallic systems. Therefore, the Hf and HfTi conducting filaments, not the bulk oxide, play the dominating role during the analog temporal response of HfOx and HfxTiyOz devices.

See the supplementary material for the detailed fabrication method, characterization of the HfOx and HfxTiyOz films, and specifics of the alpha parameter used to measure the linearity of the devices. Digital measurement results (forming, gradual reset, and stabilization), simulated temperature distribution during the analog temporal response, and standard distribution of the measured data can also be found in the supplementary material.

This work was supported by the Air Force Office of Scientific Research MURI entitled, “Cross-disciplinary Electronic-ionic Research Enabling Biologically Realistic Autonomous Learning (CEREBRAL)” under Award No. FA9550-18-1-0024. This work was performed in part at the Georgia Tech Institute for Electronics and Nanotechnology, a member of the National Nanotechnology Coordinated Infrastructure (NNCI), which is supported by the National Science Foundation (NSF) (No. ECCS-2025462). This material is based upon the work supported by the Georgia Tech ECE Fellowship and National Science Foundation Graduate Research Fellowship under Grant No. DGE-1650044.

The authors have no conflicts to disclose.

The data that support the findings of this study are available from the corresponding author upon reasonable request.

1.
B.
De Salvo
,
C.
Gerardi
,
S.
Lombardo
,
T.
Baron
,
L.
Perniola
,
D.
Mariolle
,
P.
Mur
,
A.
Toffoli
,
M.
Gely
, and
M.
Semeria
, “
How far will silicon nanocrystals push the scaling limits of NVMs technologies
,” in
IEEE International Electron Devices Meeting 2003
(
IEEE
,
2003
), pp.
26.1.1
26.1.4
.
2.
L.
Wang
,
C.
Yang
,
J.
Wen
, and
S.
Gai
, “
Emerging nonvolatile memories to go beyond scaling limits of conventional CMOS nanodevices
,”
J. Nanomater.
2014
,
927696
(
2014
).
3.
J. J.
Yang
,
M. D.
Pickett
,
X.
Li
,
D. A. A.
Ohlberg
,
D. R.
Stewart
, and
R. S.
Williams
, “
Memristive switching mechanism for metal/oxide/metal nanodevices
,”
Nat. Nanotechnol.
3
(
7
),
429
433
(
2008
).
4.
J. J.
Yang
,
D. B.
Strukov
, and
D. R.
Stewart
, “
Memristive devices for computing
,”
Nat. Nanotechnol.
8
(
1
),
13
24
(
2013
).
5.
A.
Sebastian
,
T.
Tuma
,
N.
Papandreou
,
M.
Le Gallo
,
L.
Kull
,
T.
Parnell
, and
E.
Eleftheriou
, “
Temporal correlation detection using computational phase-change memory
,”
Nat. Commun.
8
(
1
),
1
10
(
2017
).
6.
A.
Sebastian
,
M.
Le Gallo
,
R.
Khaddam-Aljameh
, and
E.
Eleftheriou
, “
Memory devices and applications for in-memory computing
,”
Nat. Nanotechnol.
15
(
7
),
529
544
(
2020
).
7.
Y.
Zhang
,
Z.
Wang
,
J.
Zhu
,
Y.
Yang
,
M.
Rao
,
W.
Song
,
Y.
Zhuo
,
X.
Zhang
,
M.
Cui
, and
L.
Shen
, “
Brain-inspired computing with memristors: Challenges in devices, circuits, and systems
,”
Appl. Phys. Rev.
7
(
1
),
011308
(
2020
).
8.
A.
Mehonic
,
A.
Sebastian
,
B.
Rajendran
,
O.
Simeone
,
E.
Vasilaki
, and
A. J.
Kenyon
, “
Memristors—from in-memory computing, deep learning acceleration, and spiking neural networks to the future of neuromorphic and bio-inspired computing
,”
Adv. Intell. Syst.
2
(
11
),
2000085
(
2020
).
9.
D.
Ielmini
and
S.
Ambrogio
, “
Emerging neuromorphic devices
,”
Nanotechnology
31
(
9
),
092001
(
2020
).
10.
D.
Ielmini
, “
Brain-inspired computing with resistive switching memory (RRAM): Devices, synapses and neural networks
,”
Microelectron. Eng.
190
,
44
53
(
2018
).
11.
P. A.
Merolla
,
J. V.
Arthur
,
R.
Alvarez-Icaza
,
A. S.
Cassidy
,
J.
Sawada
,
F.
Akopyan
,
B. L.
Jackson
,
N.
Imam
,
C.
Guo
, and
Y.
Nakamura
, “
A million spiking-neuron integrated circuit with a scalable communication network and interface
,”
Science
345
(
6197
),
668
673
(
2014
).
12.
G.
Indiveri
and
S.-C.
Liu
, “
Memory and information processing in neuromorphic systems
,”
Proc. IEEE
103
(
8
),
1379
1397
(
2015
).
13.
D.
Ielmini
and
H.-S. P.
Wong
, “
In-memory computing with resistive switching devices
,”
Nat. Electron.
1
(
6
),
333
343
(
2018
).
14.
J.
Zhu
,
T.
Zhang
,
Y.
Yang
, and
R.
Huang
, “
A comprehensive review on emerging artificial neuromorphic devices
,”
Appl. Phys. Rev.
7
(
1
),
011312
(
2020
).
15.
N. K.
Upadhyay
,
H.
Jiang
,
Z.
Wang
,
S.
Asapu
,
Q.
Xia
, and
J.
Joshua Yang
, “
Emerging memory devices for neuromorphic computing
,”
Adv. Mater. Technol.
4
(
4
),
1800589
(
2019
).
16.
S.
Lawrence
,
A.
Yandapalli
, and
S.
Rao
, “
Matrix multiplication by neuromorphic computing
,”
Neurocomputing
431
,
179
187
(
2021
).
17.
H.
Abbas
,
Y.
Abbas
,
S. N.
Truong
,
K.-S.
Min
,
M. R.
Park
,
J.
Cho
,
T.-S.
Yoon
, and
C. J.
Kang
, “
A memristor crossbar array of titanium oxide for non-volatile memory and neuromorphic applications
,”
Semicond. Sci. Technol.
32
(
6
),
065014
(
2017
).
18.
L.
Chua
,
G. C.
Sirakoulis
, and
A.
Adamatzky
,
Handbook of Memristor Networks
(
Springer Nature
,
2019
).
19.
Z.
Wang
,
S.
Joshi
,
S. E.
Savel’ev
,
H.
Jiang
,
R.
Midya
,
P.
Lin
,
M.
Hu
,
N.
Ge
,
J. P.
Strachan
, and
Z.
Li
, “
Memristors with diffusive dynamics as synaptic emulators for neuromorphic computing
,”
Nat. Mater.
16
(
1
),
101
108
(
2017
).
20.
T.
Fu
,
X.
Liu
,
H.
Gao
,
J. E.
Ward
,
X.
Liu
,
B.
Yin
,
Z.
Wang
,
Y.
Zhuo
,
D. J.
Walker
, and
J. J.
Yang
, “
Bioinspired bio-voltage memristors
,”
Nat. Commun.
11
(
1
),
1
10
(
2020
).
21.
P.
Mazumder
,
S. M.
Kang
, and
R.
Waser
, “
Memristors devices, models, and applications
,”
Proc. IEEE
100
(
6
),
1911
1919
(
2012
).
22.
M. A.
Zidan
,
J. P.
Strachan
, and
W. D.
Lu
, “
The future of electronics based on memristive systems
,”
Nat. Electron.
1
(
1
),
22
29
(
2018
).
23.
S.
Yu
,
Y.
Wu
,
R.
Jeyasingh
,
D.
Kuzum
, and
H.-S. P.
Wong
, “
An electronic synapse device based on metal oxide resistive switching memory for neuromorphic computation
,”
IEEE Trans. Electron Devices
58
(
8
),
2729
2737
(
2011
).
24.
S.
Dirkmann
,
J.
Kaiser
,
C.
Wenger
, and
T.
Mussenbrock
, “
Filament growth and resistive switching in hafnium oxide memristive devices
,”
ACS Appl. Mater. Interfaces
10
(
17
),
14857
14868
(
2018
).
25.
S.
Kumar
,
Z.
Wang
,
X.
Huang
,
N.
Kumari
,
N.
Davila
,
J. P.
Strachan
,
D.
Vine
,
A. L. D.
Kilcoyne
,
Y.
Nishi
, and
R. S.
Williams
, “
Oxygen migration during resistance switching and failure of hafnium oxide memristors
,”
Appl. Phys. Lett.
110
(
10
),
103503
(
2017
).
26.
F.
Cüppers
,
S.
Menzel
,
C.
Bengel
,
A.
Hardtdegen
,
M.
Von Witzleben
,
U.
Böttger
,
R.
Waser
, and
S.
Hoffmann-Eifert
, “
Exploiting the switching dynamics of HfO2-based ReRAM devices for reliable analog memristive behavior
,”
APL Mater.
7
(
9
),
091105
(
2019
).
27.
M.
Yin
,
Y.
Yang
,
Z.
Wang
,
T.
Zhang
,
Y.
Fang
,
X.
Yang
,
Y.
Cai
, and
R.
Huang
, “
TaOx based memristors with recessed bottom electrodes and built-in ion concentration gradient as electronic synapses
,” in
2016 13th IEEE International Conference on Solid-State and Integrated Circuit Technology (ICSICT)
(
IEEE
,
2016
), pp.
1113
1115
.
28.
S. H.
Lee
,
J.
Moon
,
Y.
Jeong
,
J.
Lee
,
X.
Li
,
H.
Wu
, and
W. D.
Lu
, “
Quantitative, dynamic TaOx memristor/resistive random access memory model
,”
ACS Appl. Electron. Mater.
2
(
3
),
701
709
(
2020
).
29.
K.
Seo
,
I.
Kim
,
S.
Jung
,
M.
Jo
,
S.
Park
,
J.
Park
,
J.
Shin
,
K. P.
Biju
,
J.
Kong
, and
K.
Lee
, “
Analog memory and spike-timing-dependent plasticity characteristics of a nanoscale titanium oxide bilayer resistive switching device
,”
Nanotechnology
22
(
25
),
254023
(
2011
).
30.
Q.
Xia
,
W.
Robinett
,
M. W.
Cumbie
,
N.
Banerjee
,
T. J.
Cardinali
,
J. J.
Yang
,
W.
Wu
,
X.
Li
,
W. M.
Tong
, and
D. B.
Strukov
, “
Memristor—CMOS hybrid integrated circuits for reconfigurable logic
,”
Nano Lett.
9
(
10
),
3640
3645
(
2009
).
31.
T.
Chang
,
S.-H.
Jo
, and
W.
Lu
,
ACS Nano
5
,
7669
7676
(
2011
);
[PubMed]
T.
Ohno
,
T.
Hasegawa
,
T.
Tsuruoka
,
K.
Terabe
,
J. K.
Gimzewski
, and
M.
Aono
,
Nat. Mater.
10
,
591
(
2011
).
[PubMed]
32.
R.
Yang
,
K.
Terabe
,
Y.
Yao
,
T.
Tsuruoka
,
T.
Hasegawa
,
J. K.
Gimzewski
, and
M.
Aono
, “
Synaptic plasticity and memory functions achieved in a WO3−x-based nanoionics device by using the principle of atomic switch operation
,”
Nanotechnology
24
(
38
),
384003
(
2013
).
33.
T. D.
Dongale
,
S. V.
Mohite
,
A. A.
Bagade
,
R. K.
Kamat
, and
K. Y.
Rajpure
, “
Bio-mimicking the synaptic weights, analog memory, and forgetting effect using spray deposited WO3 memristor device
,”
Microelectron. Eng.
183–184
,
12
18
(
2017
).
34.
H.
Nili
,
T.
Ahmed
,
S.
Walia
,
R.
Ramanathan
,
A. E.
Kandjani
,
S.
Rubanov
,
J.
Kim
,
O.
Kavehei
,
V.
Bansal
, and
M.
Bhaskaran
, “
Microstructure and dynamics of vacancy-induced nanofilamentary switching network in donor doped SrTiO3−x memristors
,”
Nanotechnology
27
(
50
),
505210
(
2016
).
35.
C.
Funck
,
C.
Bäumer
,
S.
Wiefels
,
T.
Hennen
,
R.
Waser
,
S.
Hoffmann-Eifert
,
R.
Dittmann
, and
S.
Menzel
, “
Comprehensive model for the electronic transport in Pt/SrTiO3 analog memristive devices
,”
Phys. Rev. B
102
(
3
),
035307
(
2020
).
36.
J. T.
Jang
,
D.
Kim
,
W. S.
Choi
,
S.-J.
Choi
,
D. M.
Kim
,
Y.
Kim
, and
D. H.
Kim
, “
One transistor–two memristor based on amorphous indium–gallium–zinc-oxide for neuromorphic synaptic devices
,”
ACS Appl. Electron. Mater.
2
(
9
),
2837
2844
(
2020
).
37.
Y.
Wang
,
Y.
Gong
,
L.
Yang
,
Z.
Xiong
,
Z.
Lv
,
X.
Xing
,
Y.
Zhou
,
B.
Zhang
,
C.
Su
, and
Q.
Liao
, “
MXene-ZnO memristor for multimodal In-sensor computing
,”
Adv. Funct. Mater.
31
(
21
),
2100144
(
2021
).
38.
X.
Lian
,
X.
Shen
,
M.
Zhang
,
J.
Xu
,
F.
Gao
,
X.
Wan
,
E.
Hu
,
Y.
Guo
,
J.
Zhao
, and
Y.
Tong
, “
Resistance switching characteristics and mechanisms of MXene/SiO2 structure-based memristor
,”
Appl. Phys. Lett.
115
(
6
),
063501
(
2019
).
39.
C.
He
,
J.
Li
,
X.
Wu
,
P.
Chen
,
J.
Zhao
,
K.
Yin
,
M.
Cheng
,
W.
Yang
,
G.
Xie
, and
D.
Wang
, “
Tunable electroluminescence in planar graphene/SiO2 memristors
,”
Adv. Mater.
25
(
39
),
5593
5598
(
2013
).
40.
S.
Yu
,
B.
Gao
,
Z.
Fang
,
H.
Yu
,
J.
Kang
, and
H.-S. P.
Wong
, “
A low energy oxide-based electronic synaptic device for neuromorphic visual systems with tolerance to device variation
,”
Adv. Mater.
25
(
12
),
1774
1779
(
2013
).
41.
H.-S. P.
Wong
,
H.-Y.
Lee
,
S.
Yu
,
Y.-S.
Chen
,
Y.
Wu
,
P.-S.
Chen
,
B.
Lee
,
F. T.
Chen
, and
M.-J.
Tsai
, “
Metal-oxide RRAM
,”
Proc. IEEE
100
(
6
),
1951
1970
(
2012
).
42.
P.
Basnet
,
D. G.
Pahinkar
,
M. P.
West
,
C. J.
Perini
,
S.
Graham
, and
E. M.
Vogel
, “
Substrate dependent resistive switching in amorphous-HfOx memristors: An experimental and computational investigation
,”
J. Mater. Chem. C
8
(
15
),
5092
5101
(
2020
).
43.
E.
Başar
and
B.
Güntekin
, “
A review of brain oscillations in cognitive disorders and the role of neurotransmitters
,”
Brain Res.
1235
,
172
193
(
2008
).
44.
K.
Kuriyama
and
S.
Ohkuma
, “
Role of nitric oxide in central synaptic transmission: Effects on neurotransmitter release
,”
Jpn. J. Pharmacol.
69
(
1
),
1
8
(
1995
).
45.
A. E.
Abdel Ghafar
,
S.
E. Elkowrany
,
S.
A. Salem
,
A.
A. Menaisy
,
W.
A. Fadel
, and
W.
M. Awara
, “
Effect of some parasitic infection on neurotransmitters in the brain of experimentally infected mice before and after treatment
,”
J. Egypt. Soc. Parasitol.
26
(
2
),
497
508
(
1996
).
46.
L. F.
Abbott
and
W. G.
Regehr
, “
Synaptic computation
,”
Nature
431
(
7010
),
796
803
(
2004
).
47.
H.
Ryu
and
S.
Kim
, “
Self-rectifying resistive switching and short-term memory characteristics in Pt/HfO2/TaOx/TiN artificial synaptic device
,”
Nanomaterials
10
(
11
),
2159
(
2020
).
48.
Q.
Wang
,
G.
Niu
,
S.
Roy
,
Y.
Wang
,
Y.
Zhang
,
H.
Wu
,
S.
Zhai
,
W.
Bai
,
P.
Shi
, and
S.
Song
, “
Interface-engineered reliable HfO2-based RRAM for synaptic simulation
,”
J. Mater. Chem. C
7
(
40
),
12682
12687
(
2019
).
49.
S.
Kim
,
J.
Chen
,
Y.-C.
Chen
,
M.-H.
Kim
,
H.
Kim
,
M.-W.
Kwon
,
S.
Hwang
,
M.
Ismail
,
Y.
Li
, and
X.-S.
Miao
, “
Neuronal dynamics in HfOx/AlOy-based homeothermic synaptic memristors with low-power and homogeneous resistive switching
,”
Nanoscale
11
(
1
),
237
245
(
2019
).
50.
Y.
Matveyev
,
K.
Egorov
,
A.
Markeev
, and
A.
Zenkevich
, “
Resistive switching and synaptic properties of fully atomic layer deposition grown TiN/HfO2/TiN devices
,”
J. Appl. Phys.
117
(
4
),
044901
(
2015
).
51.
K.
Moon
,
S.
Lim
,
J.
Park
,
C.
Sung
,
S.
Oh
,
J.
Woo
,
J.
Lee
, and
H.
Hwang
, “
RRAM-based synapse devices for neuromorphic systems
,”
Faraday Discuss.
213
,
421
451
(
2019
).
52.
C.
Baeumer
,
R.
Valenta
,
C.
Schmitz
,
A.
Locatelli
,
T. O.
Mentes
,
S. P.
Rogers
,
A.
Sala
,
N.
Raab
,
S.
Nemsak
, and
M.
Shim
, “
Subfilamentary networks cause cycle-to-cycle variability in memristive devices
,”
ACS Nano
11
(
7
),
6921
6929
(
2017
).
53.
D.
Kuzum
,
R. G. D.
Jeyasingh
,
B.
Lee
, and
H.-S. P.
Wong
, “
Nanoelectronic programmable synapses based on phase change materials for brain-inspired computing
,”
Nano Lett.
12
(
5
),
2179
2186
(
2012
).
54.
G.
Indiveri
,
B.
Linares-Barranco
,
R.
Legenstein
,
G.
Deligeorgis
, and
T.
Prodromakis
, “
Integration of nanoscale memristor synapses in neuromorphic computing architectures
,”
Nanotechnology
24
(
38
),
384010
(
2013
).
55.
Y.
Li
and
K.-W.
Ang
, “
Hardware implementation of neuromorphic computing using large-scale memristor crossbar arrays
,”
Adv. Intell. Syst.
3
(
1
),
2000137
(
2021
).
56.
Y.-E.
Syu
,
T.-C.
Chang
,
J.-H.
Lou
,
T.-M.
Tsai
,
K.-C.
Chang
,
M.-J.
Tsai
,
Y.-L.
Wang
,
M.
Liu
, and
S. M.
Sze
, “
Atomic-level quantized reaction of HfOx memristor
,”
Appl. Phys. Lett.
102
(
17
),
172903
(
2013
).
57.
X.
Ding
,
D.
Zhu
,
Y.
Feng
,
L.
Cai
,
P.
Huang
,
L.
Liu
, and
J.
Kang
, “
Low-power, multilevel and analog characteristics in the multi-layer-oxide based RRAM devices compatible with CMOS technology
,” in
2018 14th IEEE International Conference on Solid-State and Integrated Circuit Technology (ICSICT)
(
IEEE
,
2018
), pp.
1
3
.
58.
S.
Roy
,
G.
Niu
,
Q.
Wang
,
Y.
Wang
,
Y.
Zhang
,
H.
Wu
,
S.
Zhai
,
P.
Shi
,
S.
Song
, and
Z.
Song
, “
Toward a reliable synaptic simulation using Al-doped HfO2 RRAM
,”
ACS Appl. Mater. Interfaces
12
(
9
),
10648
10656
(
2020
).
59.
N.-N.
Ge
,
C.-H.
Gong
,
X.-C.
Yuan
,
H.-Z.
Zeng
, and
X.-H.
Wei
, “
Effect of Mn doping on electroforming and threshold voltages of bipolar resistive switching in Al/Mn:NiO/ITO
,”
RSC Adv.
8
(
52
),
29499
29504
(
2018
).
60.
Z.
Fang
,
H. Y.
Yu
,
X.
Li
,
N.
Singh
,
G. Q.
Lo
, and
D. L.
Kwong
, “
HfOx/TiOx/HfOx/TiOx multilayer-based forming-free RRAM devices with excellent uniformity
,”
IEEE Electron Device Lett.
32
(
4
),
566
568
(
2011
).
61.
T.
Tan
,
Y.
Du
,
A.
Cao
,
Y.
Sun
,
G.
Zha
,
H.
Lei
, and
X.
Zheng
, “
The resistive switching characteristics of Ni-doped HfOx film and its application as a synapse
,”
J. Alloys Compd.
766
,
918
924
(
2018
).
62.
S.
Mandal
,
A.
El-Amin
,
K.
Alexander
,
B.
Rajendran
, and
R.
Jha
, “
Novel synaptic memory device for neuromorphic computing
,”
Sci. Rep.
4
(
1
),
1
10
(
2014
).
63.
J.
Ahn
,
J. Y.
Lee
,
J.
Kim
,
J.
Yoo
, and
C.
Ryu
, “
Comparison study from sputtering, sol-gel, and ALD processes developing embedded thin film capacitors
,” in
2006 8th Electronics Packaging Technology Conference
(
IEEE
,
2006
), pp.
10
14
.
64.
W.
Wu
,
H.
Wu
,
B.
Gao
,
P.
Yao
,
X.
Zhang
,
X.
Peng
,
S.
Yu
, and
H.
Qian
, “
A methodology to improve linearity of analog RRAM for neuromorphic computing
,” in
2018 IEEE Symposium on VLSI Technology
(
IEEE
,
2018
), pp.
103
104
.
65.
S.
Chandrasekaran
,
F. M.
Simanjuntak
,
R.
Saminathan
,
D.
Panda
, and
T.-Y.
Tseng
, “
Improving linearity by introducing Al in HfO2 as a memristor synapse device
,”
Nanotechnology
30
(
44
),
445205
(
2019
).
66.
B.
Chakrabarti
,
Investigation of Metal Oxide Dielectrics for Non-Volatile Floating Gate and Resistance Switching Memory Applications
(
The University of Texas at Dallas
,
2014
).
67.
B.
Chakrabarti
and
E. M.
Vogel
, “
Effect of Ti doping and annealing on multi-level forming-free resistive random access memories with atomic layer deposited HfTiOx nanolaminate
,”
Microelectron. Eng.
109
,
193
196
(
2013
).
68.
F. F.
Athena
,
M. P.
West
,
J.
Hah
,
R.
Hanus
,
S.
Graham
, and
E. M.
Vogel
, “
Towards a better understanding of the forming and resistive switching behavior of Ti-doped HfOx RRAM
,”
J. Mater. Chem. C
10
,
5896
5904
(
2022
).
69.
G.
Henkelman
,
B. P.
Uberuaga
, and
H.
Jónsson
, “
A climbing image nudged elastic band method for finding saddle points and minimum energy paths
,”
J. Chem. Phys.
113
(
22
),
9901
9904
(
2000
).
70.
E. A.
Cartier
,
W.
Kim
,
N.
Gong
,
T.
Gokmen
,
M. M.
Frank
,
D. M.
Bishop
,
Y.
Kim
,
S.
Kim
,
T.
Ando
, and
E. Y.
Wu
, “
Reliability challenges with materials for analog computing
,” in
2019 IEEE International Reliability Physics Symposium (IRPS)
(
IEEE
,
2019
), pp.
1
10
.
71.
J.
Lee
,
J.-H.
Ryu
,
B.
Kim
,
F.
Hussain
,
C.
Mahata
,
E.
Sim
,
M.
Ismail
,
Y.
Abbas
,
H.
Abbas
, and
D. K.
Lee
, “
Synaptic characteristics of amorphous boron nitride-based memristors on a highly doped silicon substrate for neuromorphic engineering
,”
ACS Appl. Mater. Interfaces
12
(
30
),
33908
33916
(
2020
).
72.
Z.
Li
,
B.
Tian
,
K.-H.
Xue
,
B.
Wang
,
M.
Xu
,
H.
Lu
,
H.
Sun
, and
X.
Miao
, “
Coexistence of digital and analog resistive switching with low operation voltage in oxygen-gradient HfOx memristors
,”
IEEE Electron Device Lett.
40
(
7
),
1068
1071
(
2019
).
73.
C.
Bloch
, “
Power
,”
Tikkun
23
(
3
),
49
(
2008
).
74.
S.
Agarwal
,
S. J.
Plimpton
,
D. R.
Hughart
,
A. H.
Hsia
,
I.
Richter
,
J. A.
Cox
,
C. D.
James
, and
M. J.
Marinella
, “
Resistive memory device requirements for a neural algorithm accelerator
,” in
2016 International Joint Conference on Neural Networks (IJCNN)
(
IEEE
,
2016
), pp.
929
938
.
75.
J.
Chen
,
W.-Q.
Pan
,
Y.
Li
,
R.
Kuang
,
Y.-H.
He
,
C.-Y.
Lin
,
N.
Duan
,
G.-R.
Feng
,
H.-X.
Zheng
, and
T.-C.
Chang
, “
High-precision symmetric weight update of memristor by gate voltage ramping method for convolutional neural network accelerator
,”
IEEE Electron Device Lett.
41
(
3
),
353
356
(
2020
).
76.
S. R.
Bradley
,
A. L.
Shluger
, and
G.
Bersuker
, “
Electron-injection-assisted generation of oxygen vacancies in monoclinic HfO2
,”
Phys. Rev. Appl.
4
(
6
),
064008
(
2015
).
77.
G.
Bersuker
,
D.
Gilmer
, and
D.
Veksler
, “
Metal-oxide resistive random access memory (RRAM) technology: Material and operation details and ramifications
,” in
Advances in Non-Volatile Memory and Storage Technology
(
Elsevier
,
2019
), pp.
35
102
.
78.
H. H.
Wu
,
P.
Wisesa
, and
D. R.
Trinkle
, “
Oxygen diffusion in hcp metals from first principles
,”
Phys. Rev. B
94
(
1
),
014307
(
2016
).

Supplementary Material