Electrically programmable magnetoresistance in \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\text{AlO}_{x}$$\end{document}AlOx-based magnetic tunnel junctions

We report spin-dependent transport properties and I–V hysteresis characteristics in an \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\text{AlO}_{x}$$\end{document}AlOx-based magnetic tunnel junction (MTJ). The bipolar resistive switching and the magnetoresistances measured at high resistance state (HRS) and low resistance state (LRS) yield four distinctive resistive states in a single device. The temperature dependence of resistance at LRS suggests that the resistive switching is not triggered by the metal filaments within the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\text{AlO}_{x}$$\end{document}AlOx layer. The role played by oxygen vacancies in \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\text{AlO}_{x}$$\end{document}AlOx is the key to determine the resistive state. Our study reveals the possibility of controlling the multiple resistive states in a single \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\text{AlO}_{x}$$\end{document}AlOx-based MTJ by the interplay of both electric and magnetic fields, thus providing potential applications for future multi-bit memory devices.

Manipulation and control of transport properties of spin-polarized electrons in nano-structured magnetic multilayers, such as Giant Magnetoresistance (GMR) in spin valves and Tunneling Magnetoresistance (TMR) in magnetic tunnel junctions (MTJs) [1][2][3][4] , offer versatile possibilities for developing novel spintronic applications in information storage and memory devices. However, to meet the increasing demand for faster, smaller and lower energy consumption in the relevant technologies, the development of devices equipped with multi-functionalities and nonvolatile electronic properties is of intensive interest. Recently, there have been a lot of efforts aiming at the electric-field control of magnetism in MTJs to manipulate their magneto-electronic transport properties [5][6][7] . Such a strategy could lead to desirable ultra-low energy consumption, a demand increasingly sought after in nano-scaled devices. Furthermore, by introducing a capacitive charge accumulation at the multilayer interfaces, the density of states at the Fermi level could be altered by an external electric field, which in turn could lead to the changes of magnetism and magnetic anisotropy [8][9][10] .
In contrast, another type of electric-field controlled, nonvolatile resistive switching (RS) effect based on the reversible resistive changes between the high and the low resistance states (HRS and LRS) of a metal/insulator/ metal structure (identical to an MTJ, except the use of a non-ferromagnetic metal as the electrodes) is regarded one of the most promising new memories with the advantages of high integration, low power consumption, high read-write speed, non-volatility and compatibility with the existing CMOS (Complementary Metal-Oxide-Semiconductor) technology [11][12][13] . It would be an important step towards multi-functionality if both resistive switching and magnetoresistance can be integrated into a single device and its multiple states can be electrically/ magnetically controlled at the same time. Previous studies have investigated the RS effect in AlO x and MgObased MTJs [14][15][16] ; nevertheless, after the electroforming of the resistive memory operation process, it is difficult to simultaneously obtain both observable on-off and MR ratios. In most cases either the TMR is absent or the on-off ratio is suppressed 17,18 . In this study, we demonstrate that both resistive switching and magnetoresistance in an AlO x -based MTJ can be achieved via electric/magnetic control of the resistance states. The bipolar resistive switching ratio between HRS and LRS is up to 202% and the device shows a simultaneous tunneling MR ratio of 10% and 4 % at HRS and LRS, respectively, promising a multi-bit memory of four distinctive memory states.

Experimental section
The structure of the AlO x -based MTJ is stacked in the sequence of NiFe (15 nm)/CoFe (10 nm)/AlO x (1.5 nm)/ CoFe (35 nm) and patterned in a crossed-bar configuration, with the bottom NiFe/CoFe layer serves as the soft ferromagnetic (FM) electrode and the top CoFe as the hard FM electrode, as illustrated in the lower-right inset of Fig. 1 www.nature.com/scientificreports/ mbar. The AlO x layer was prepared by nature oxidation of Al followed by an O 2 plasma treatment at 3 mbar for 20 seconds. Details of the sputtering and the patterning processes can be retrieved in our previous work 19 . The current-voltage (I-V) characteristics were measured at various voltages and the magnetoresistance was measured in the Current-Perpendicular-to-the-Plane (CPP) configuration using the four-point-probe method with the bias voltage applied on the top electrode while the bottom electrode was grounded, whereas the temperature dependence of the resistance was measured in a Quantum Design Physical Property Measurement System (PPMS).

Results and discussion
The current-voltage (I-V) curves of the device after electroforming are shown in Fig. 1, which exhibit reversible, nonvolatile bipolar resistive switching characteristic. The room temperature MR effect can be observed at both HRS and LRS (see the discussion of Fig. 2 below). The newly fabricated device shows a high resistance ( ∼ 1.3 k ) in its pristine state. After the electroforming process at ∼ 1 volt (see the upper-left inset of Fig.1), the device switches from the pristine state to the low resistance state (LRS). Here a compliance current (ICC) of 9.4 mA is applied to avoid excessive current flow and to prevent the irreversible hard breakdown. The sudden increase of current level during the electroforming process causes localized heating and the formation of oxygen vacancies in the AlO x barrier, producing a new remnant insulating gap between the two ferromagnet electrodes along with an irreversible decrease of resistance from the pristine state 20 . Following the electroforming process, the voltage is swept between + 1.4 and − 1.4 V in a cyclic manner with the arrows in Fig. 1 indicating the directions of the voltage sweeps. When increasing the voltage from 0 to + 1.4 V (step 1 in Fig. 1) during the sweep, the current increases abruptly at ∼ 1.1 V which transits the device to low resistance state (LRS). From 1.1 to 1.4 V, despite an applied ICC of 9.4 mA, a flattened slope is present which can be attributed to the self-rectifying characteristic of AlO x 21 . On the negative bias side, the same behavior is observed at ∼ −1.0 V which switches the device from LRS to HRS (step 3 in Fig. 1). Here we designate LRS as the "ON" or "'1" state, and HRS as the "OFF" or "0" state. The on-off ratio of the device is 202% at 0.9 V. To further demonstrate the non-volatility of the resistive switching, a test was performed in such a way that the applied voltage was suddenly turned off at 1.4 V. An I-V curve was taken immediately after between + 0.5 and − 0.5 V (represented by the green open circles in the figure). As shown in Fig. 1, the curve follows closely the footprint of the major loop. The same procedure applied on the negative side results in a similar behavior (see the orange open circles), showing that the resistance states are reproducible. The I-V and the corresponding MR measurements have been performed for more than 100 cycles, during which the MR variation is within 10% and 29% for HRS and LRS, respectively. The overall decrease of MR after 100 cycles is about 3 % for HRS and 11% for LRS. These systematic results suggest a robust nonvolatile memory characteristic of the device. Figure 2 shows the magneto-transport characterization of the AlO x -based MTJ, in which the cycling of junction resistance is plotted against the applied magnetic field. The resistance curves display the typical pseudo spin-valve type characteristics at room temperature (RT) with a magneto-resistive (MR) ratio defined as where R P and R AP depict the resistances when the magnetizations of the FM electrodes are in the parallel and the antiparallel configurations, respectively. Electrical control of magnetoresistance has been regarded an alternative  Fig. 2b). Repeating the procedure indicates that the effect is reversible. Figure 2c gives a resistance of 473 (HRS) and an MR ratio of 10.7% upon setting the bias voltage back to -1.4 V. Figure 2d-f show that by further switching the bias between − 1.4 and + 1.4 V, the dual resistance states and MR states are both preserved. The same behavior was observed in other MTJs made in the same run under identical conditions. In a set of 19 simultaneously-made samples, 13 samples exhibit a uniform distribution of junction resistances with variation within 10% in the pristine state, while the other 6 samples show larger discrepancies due to the edge effect of the deposition rate. By carefully controlling the compliance current during the electroforming process, the resistance of either HRS or LRS can be controlled within 20% from sample to sample. This variation is likely due to the differences of oxygen vacancies in AlO x . Figure 2 presents strong evidences for multifunctional operation in these MTJ devices. Both junction resistance and magnetoresistance can be manipulated reversibly by the interplay of electric/magnetic fields in a nonvolatile manner. The results lead directly to the possibility of integrating multiple-state memories into a single device, which goes beyond the binary (0-1) limit.
To further investigate the mechanism of resistive switching in our device, we re-plot the positive sweep region of Fig. 1 in double-logarithmic scales (see Fig. 3a), in which the fitting results are indicated by black lines. In the low-voltage regime (0-0.5 V) of HRS, the curve exhibits a linear relation between current and voltage with a slope of ≈ 1.054. When V ≥ 0.66 V, the curve becomes nonlinear, and a linear regression results in an enhanced slope of ≈ 2.352. By further increasing the voltage, the current saturates to a constant value. Upon passing the threshold voltage ( ∼1.1 V), the device switches from HRS to LRS. When reversing and decreasing the voltage, the I-V curve goes through a nonlinear region (above 0.48 V, with a regression slop of 1.950) and a linear region (from 0.48 to 0 V, with a slope of 1.043). The power law fittings ( I ∝ V m ) to both the HRS and the LRS curves yield two distinctive regimes, i.e. the low bias-voltage regime with m ≈ 1 which exhibits Ohmic characteristics (I ∝ V), and the high bias regime with m ≈ 2 which obeys the Child's law ( I ∝ V 2 ) 23 . The conduction mechanism in the low bias regime is dominated by the thermally generated free electrons, such as the Schottky or the thermionic emission, which is one of the most often observed conduction mechanisms in oxides and semiconductors 24 . After which, the field exceeds the square-law onset voltage and the electron density gradually surpasses the equilibrium concentration and dominates the conduction in the higher bias regime. Such a behavior is common in oxides consisting of a discrete trapping level with both trap-unfilled and trap-filled transport characteristics 23 . Figure 3b shows the fitting results (blue curve for HRS and red curve for LRS) of the low-voltage I-V curves using the Simmons equation I = aV + bV 3 , where a and b are two barrier shape-related coefficients 25 . Assuming a rectangular type barrier, both the HRS and the LRS curves can be well fitted, indicating a direct tunneling www.nature.com/scientificreports/ transport mechanism 26 . Typically, for metal/insulator/metal junctions, when V > φ (where φ is the barrier height), a change occurs in the barrier shape from rectangular to triangular, leading to a transition of transport mechanism from direct tunneling (DT) to field emission, also known as the Fowler-Nordheim tunneling (FN) 27 .
In this high-voltage regime, the current is described by where W is the barrier width, m e is the effective electron mass in the oxide, is the reduced Planck's constant, φ is the junction barrier height and q is the electronic charge. The transition from DT to FN has been reported in RRAM (resistive random-access memory) devices 28 . A close fit of Eq. (2) is often considered an indication of trap-assisted tunneling (TAT), which has been verified a predominant current transport mechanism in RRAM devices. In Fig. 3c, we plot ln(I/V 2 ) as a function of 1/V at room temperature, which shows clear transitions from DT to FN at V T ∼ 0.83 V (HRS) and 0.64 V (LRS). In the present study, a control sample with fully oxidized Al 2 O 3 barrier exhibits an MR ratio of 30% at RT as prepared (see Fig. 3d), and the I-V curve shows no resistive switching. Compared to the I-V curve of the MTJ with partially oxidized barrier (Fig. 1), we conclude that the resistive switching is likely resulted from the oxygen deficiencies created during the preparation process. More details of the transport mechanism can be obtained by analyzing the temperature dependent characteristics. Figure 4a gives the temperature dependence of the MTJ junction resistance. The measurements were performed at a constant current mode at 10 μA, which corresponds to the low voltage regime (from ∼ 4.7 to 13 mV for HRS, and from ∼ 1.67 to 5 mV for LRS). Both the HRS and the LRS resistances increase with decreasing temperature between 300 K and 10 K, indicating a semiconducting behavior in this temperature range for both resistance states. In the inset of Fig. 4a, we have plotted the temperature dependence of conductance of the MTJ by applying the three-dimensional hopping model: where G, T, and A denote the conductance, temperature, and a constant related to the bulk-related inverse localization radius ( l loc ), respectively 29 . If the three-dimensional hopping model is obeyed, the relation between ln(△G ) and T −1/4 should be linear. Linear fits are obtained between RT and 25 K for both LRS and HRS. This linear behavior is a typical characteristic of inter-particle hopping in metal/semiconductor/metal thin films 30,31 . A close examination of the inset of Fig. 4a shows that the LRS curve agrees to the hopping model better than the www.nature.com/scientificreports/ HRS. A more complicated analysis combining the inter-particle tunneling and the higher-order inelastic hopping terms is described by Eq. (4), and the fitting results are shown in Fig. 4b for both HRS and LRS.
here G tun is the tunneling conductance, G hop is the spin-independent higher-order inelastic hopping conductance, G 0 and C are free parameters, △ = 4E/k B , where E is the tunneling activation energy, k B is the Boltzmann constant, and γ = N − [2/(N + 1)] with N indicating the number of localized states in the barrier 32 .
The temperature dependence of conductivity at either LRS or HRS can be described by Eq. (4), which implies that the transport mechanism comprises partially the thermal-assisted inelastic hopping through inter-particles. As seen in Fig. 4b, our experimental data are well fitted by Eq. (4) with γ = 1.33 ( N = 2 ) at LRS, suggesting that the conduction is contributed mainly by the tunneling channels and the second-order hopping. As for HRS, the best fit contains not only the second-order hopping ( γ = 1.33 ) but also the higher-order terms (i.e. third-order with γ = 2.5 , and fourth-order with γ = 3.6 ). In metal oxide based devices, the conduction rises primarily through the conduction channels formed in the regions with high concentration of oxygen vacancies ( V O ). Our results show that in addition to the tunneling channels, higher-order inelastic hopping due to the presence of localized states in the AlO x barrier also contributes significantly to the charge transport process. Such transport process is intriguing in MTJs and nanoparticle networks exhibiting sequential and co-tunneling characteristics, especially at low temperature 33,34 . Further studies at low temperature are on demand to gain insights into the temperature and bias dependences of the sequential and co-tunneling regime.
By correlating the results of Figs. 3 and 4, we conclude that the high MR ratio at HRS corresponds to the insulating regime of the device where MR arises from the tunneling effect, whereas the low MR ratio at LRS is due to the high level of localized/trapping states. In such a case, the magnetoresistance results mainly from the hopping of spin-polarized electrons through the localized states.

Conclusion
In summary, we have demonstrated that in an AlO x -based MTJ device, both resistive switching and magnetoresistance are observed and are electrically tunable. The room-temperature bipolar resistive switching ratio is ∼ 202% and the MR ratio is 10% and 4 %for HRS and LRS, respectively. The temperature dependence of resistance suggests that the switching mechanism is not the typical filament-type switching. The fitting of conductance to a model combining the inter-particle tunneling and the higher order hopping channels reveals that by electrically altering the oxygen vacancies in the AlO x barrier, the magnetoresistance can be tuned simultaneously. The results lead to a promising possibility of integrating MTJ into the future multi-bit spintronic memory devices.