Control of electrical resistance and magnetoresistance by electric-field-driven oxygen ion migration in a single GdOx wire

Electric-field-driven ion migration can significantly modulate the electric and magnetic properties of solids, creating novel functionalities useful for advanced electromagnetic devices. Earlier works have used vertically stacked structures for this purpose, in which the redox process results from ion migration driven by a vertical electric field through the interfaces. However, the existence of the interfaces between the dissimilar layers causes the oxidation and reduction processes to have high and asymmetric energy barriers, which means that a large electric field is required to control the devices. Here, we show that in a partially oxidized single GdOx wire using a lateral electric field configuration, low and symmetric energy barriers for the oxidation and reduction processes can be achieved. We provide evidence that the redox process is the result of the lateral motion of oxygen ions by directly visualizing the electric-field-driven real-time ionic motion using an optical microscope. An electric field as low as 105 V/m was able to drive oxygen ions at room temperature, allowing controllable modulation of the electrical resistance using a lateral electric field. A large negative magnetoresistance was also observed in the GdOx wire, and its magnitude was significantly enhanced up to 20% at 9 T through oxygen ion control. Our results suggest that the electrical and magnetic properties of single GdOx can be efficiently controlled through oxygen ion motion driven by a lateral electric field, which paves the way for fully functional electromagnetic devices such as artificial synapses. A new approach for creating tiny wires that move oxygen ions instead of electrons can offer researchers efficient control over devices with brain-like conductivity. According to neuroscientists, synapses, the junctions where two neurons meet, store information by activating different types of conductive states. Soogil Lee and Kab-Jin Kim from the KAIST institute in Daejeon, South Korea, and colleagues report similar behavior in nanoscale-thin wires made from gadolinium oxide, a rare-earth material in which oxygen ions can move freely. By pushing ions in the wire laterally between a cathode and an anode using an electric field, the team could establish two different conductivity levels due to slight changes in redox chemistry. The simplicity of this design enables lower-energy ion manipulation than previous studies employing vertically stacked films of gadolinium oxide and other metals. Electric-field-driven ion migration is an effective way to control many physical properties such as electric and magnetic properties of materials. In this research, we demonstrated that remarkably small electric field as low as 105 V/m can effectively modulate resistance of a GdOx microwire due to the lateral motion of oxygen ions. This result enables mimicking a synaptic device with low power consumption like as a brain. Furthermore, huge MR modulation of the GdOx microwire with the electric field provides novel functionality as an electromagnetic device.


Introduction
Gadolinium (Gd) is one of the most interesting rareearth metals. It exhibits fascinating electrical and magnetic properties [1][2][3][4] and has practical value in applications including hydrogen storage 5 and CMOS technology as gate electrodes 6 and spintronics [7][8][9][10][11][12] . One distinctive property of Gd is that it is very easily oxidized, forming oxide compounds, GdO x . Since the oxygen ions are mobile in GdO x 13, 14 , it has been used as an oxygen source (or sink) in ferromagnet/GdO x heterostructures [14][15][16][17][18][19] . In these heterostructures, the electric and magnetic properties of the ferromagnet can be modulated by the redox process through electric-field-driven O 2− migration between the ferromagnet and GdO x . However, the vertical oxygen migration in the heterostructures has an inherent asymmetry, which leads to different energy barriers for oxidation and reduction 18 . Furthermore, the energy barrier for oxygen migration is generally large, so the use of a high electric field is required to control the device 16,17 . Therefore, for GdO x to be properly utilized in practical devices, it is highly necessary to achieve a symmetric and low energy barrier for oxygen ion motion.
In this work, we investigate the possibility of utilizing oxygen ions in a single GdO x device instead of heterostructures. Contrary to previous works [14][15][16][17][18][19] , we explore the lateral motion of oxygen ions by applying a lateral electric field. By employing optical and electrical measurements, we find that the oxygen ions are transported along the GdO x layer in a controllable way with almost the same energy barriers for oxidation and reduction. Furthermore, the energy barrier is found to be much smaller than in previous reports, which we ascribe to the low concentration of oxygen ions in our GdO x . The possibility of practical devices is presented by showing the modulation of the electrical resistance and magnetoresistance (MR) through oxygen ion control in a single GdO x device.

Sample preparation
For this study, a 10-nm-thick Gd film was deposited on a Si/SiN x (200 nm) substrate by direct current (DC) magnetron sputtering under a working pressure of 3 mTorr. To investigate the oxygen ion migration effect, we used three different capping layers, as shown in Fig. 1a. The first sample had 3-nm-thick SiN x as a capping layer. For this sample, we intentionally increased the radio frequency (RF) sputtering power up to 150 W to create a bumpy surface. This allowed oxygen to penetrate into the Gd layer and form GdO x when the sample was exposed to air. The second sample was capped by a 1.5nm-thick Ta layer using DC magnetron sputtering. When the sample was exposed to air, the Ta oxidized, and the oxygen ions in TaO x diffused into the Gd layer and formed GdO x . The third sample was prepared as a control sample with a much thicker SiN x capping layer. For this, 20-nm-thick SiN x was deposited by RF magnetron sputtering. To make a smooth surface, we reduced the RF sputtering power to as low as 50 W. We expected that in this case, Gd would not be oxidized at all and would remain in a metallic state (see Supplementary Information 1 for more details). For the electrical measurements, the films were patterned into microwires with a Hall-bar structure by conventional photolithography and Ar-ion milling techniques, as shown in Fig. 1b. The wire lengths were set to 55 μm with varying widths of 5, 10, and 20 μm. We used the four-probe measurement technique, which separates the contacts for the input current and output voltage (see Fig. 1b).

Results and discussion
Electrical and optical measurements of ion migration Figure 1c-e shows the temporal variation in the longitudinal resistivity (ρ xx ) of the three different samples. To monitor the differences in resistance, we fixed the current, whose magnitude corresponded to an electric field of 10 4 -10 5 V/m, depending on the resistance. The resistivity of partially oxidized Gd, that is, SiN x /GdO x and TaO x / GdO x (Fig. 1c, d), decreased drastically over time, while that of metallic Gd was almost constant over time (Fig. 1e).
This large decrease in resistivity indicates that a reduction process, from insulator-like GdO x to metal-like Gd, occurs when the oxygen ions are driven by an electric field. To provide more direct evidence, we measured the real-time color change of TaO x /GdO x using an optical microscope because the change in oxygen ion concentration is reflected by changes in the optical contrast 20 . We first applied a negative electric field for 1000 s to push the oxygen ions to the right edge of the wire and then reversed the electric field polarity. Figure 1f shows optical microscope images taken after reversing the polarity of the electric field. We set the initial image (t = 0) as a background image and subtracted it from all other images taken at t > 0. The black region indicates a decrease in oxygen ions, while the white region represents an increase in oxygen ions. The result shows that the accumulated oxygen ions move towards the anode under the electric field. Figure 1g shows the variations in the resistances of electrodes A-B and electrodes B-C, which were measured at the same time as the optical images. The results show that there is a time delay in the resistance variation between the physically separated regions, which suggests that the variation in resistivity is caused by the oxygen ion migration driven by the electric field. To further confirm the oxygen migration, we measured the time delay by changing the electric field. Figure 1h clearly shows that the time delay decreases with increasing electric field, which is consistent with the oxygen migration model. We observed similar oxygen ion motion in SiN x /GdO x , but for metallic Gd, we did not observe any color change (Supplementary Information 2). All these observations suggest that the variation in resistivity observed in GdO x originates from the electric-field-driven oxidation and reduction processes induced by lateral oxygen ion migration (further discussions are presented in Supplementary Information 3 and 4).

Energy barriers for ion migration
We next studied the oxygen ion migration in a more quantitative way. Figure 2a shows the time-dependent resistance variation in GdO x at temperatures from 75°C to 90°C at a 5°C interval. Here, the temperatures were set higher than room temperature to obtain a fast resistance change (i.e., rapid oxygen ion migration). We applied an electric field of 1 × 10 4 V/m. As we increase the temperature, the slope of the resistance variation increases for both the oxidation and reduction processes. This implies that the ionic motion occurs through a thermal activation process. To quantitatively determine the activation energy barrier, we plot the logarithm of the resistance variation rate (ΔR xx =Δt) with respect to the inverse of temperature in Fig. 2b. The linear dependence in Fig. 2b suggests that the ionic motion follows the Arrhenius law, which we can use to estimate the activation energy barrier.
The obtained energy barriers were found to be 0.49 ± 0.04 and 0.48 ± 0.01 eV for the oxidation and reduction processes, respectively. These values are much smaller than the binding energy of Gd-O in gadolinium oxide (Gd 2 O 3 ), but comparable to a previously reported oxygen ion diffusion energy barrier 13 . This suggests that the motion of oxygen ions in GdO x occurs via successive jumps from Gd 2 O 3 to Gd. In this process, the energy used to destroy the old bond and that used to create the new bond are almost the same, so the energy barrier to movement of the oxygen ions can be small 14 . It should be noted that the energy barriers obtained here are different from those in vertically stacked heterostructures. In heterostructures, the energy barriers are higher than ours and are asymmetric for oxidation and reduction (the energy barrier for the reduction process is higher than that for the oxidation process) 18 . Considering that the diffusion energy barrier depends on the oxygen ion concentration 21 , the smaller energy barrier in our device is possibly due to the low oxygen concentration of our GdO x , which can be supported by the low resistivity of our GdO x (~10 3 μΩ·cm) compared to that of GdO x in previous works (~10 11 μΩ·cm) 18 . The asymmetry in the energy barrier in the heterostructures arises because the oxidation and reduction occur as a result of oxygen ion motion across the interfaces between dissimilar layers. In our case, on the other hand, oxidation and reduction occur in a single GdO x wire layer through the lateral motion of oxygen ions, resulting in similar energy barriers for oxidation and reduction.
It should also be noted that despite having similar energy barriers, the oxidation and reduction processes are not exactly equivalent in our GdO x wires. This can be inferred from the slight offset shown in Fig. 2b. The offset indicates that the reduction process is slower than the oxidation process due to the concentration gradient of oxygen ions 22 . In the reduction process, the distributed oxygen ions are driven by an electric field and accumulate at the anode. In the oxidation process, on the other hand, the accumulated oxygen ions are spread out over the wire. As a result, the concentration-gradient-induced ion diffusion is additive to the electric-field-induced ion migration for the oxidation process, while it is subtractive for the reduction process. This generates the slightly different resistance variation rates for oxidation and reduction (for a more quantitative analysis, see Supplementary Information 5 and 6).
Since the oxygen ions are driven by a thermal activation process with finite energy barriers, ionic motion can be generally suppressed at low temperature. Figure 2c shows the time-dependent resistance variation in the GdO x wire at various temperatures. The results show that the resistance variation is gradually suppressed with decreasing temperature, and the resistance becomes stable below 200 K, suggesting that the thermal activation of oxygen ions is suppressed below 200 K.

Mimicking the synaptic function using ion migration
The electric-field-driven oxygen ion migration and resulting resistance variation enable control of the device resistance in a more active way. To verify this, we repeatedly switched the electric field polarity at the lower bound (R low xx = 3.00 kΩ) and upper bound (R high xx = 3.06 kΩ) as shown in Fig. 3a. The device resistance periodically changes between the two resistance states, indicating that oxygen ions move back and forth within the wire. This controllable quasilinear variation in device resistance can be used in artificial synaptic devices 22 . Figure 3a-c shows the variation in resistance for various electric fields (0.55-1.64 × 10 5 V/m). The higher-amplitude and longer-duration electric fields lead to a larger change in resistance. If we treat the device resistance as a synaptic weight, then this phenomenon closely mimics the transmission characteristics of a biological synapse. Furthermore, the sign of the resistance variation can also be controlled by the electric field polarity, which can be used to mimic the potentiation and depression functions of a biological synapse.

Effect of ion migration on the MR
We finally investigated the effect of oxygen ions on the MR of GdO x . To investigate the MR of GdO x , we first controlled the device resistivity by applying an electric field at room temperature, as shown in Fig. 1c, d, and then cooled it down to a low temperature to stabilize the Fig. 2 Determination of energy barriers for ion migration. a Temporal dependence of longitudinal resistance R xx at various temperatures in the SiN x (3nm)/GdO x (10nm) device. Here, R xx was measured by repeatedly switching the electric field polarity at the lower bound (R low xx =4.99kΩ) and upper bound (R high xx =5.01kΩ), as shown in Fig. 3. An electric field of 1×10 4 V/m was applied. A positive slope indicates the oxidation process, and a negative slope corresponds to the reduction process. Here, ΔTime=Time-Time at 5.00kΩ. b Arrhenius plot obtained from a. The detailed fitting function is discussed in Supplementary Information 5 and 6. c Variation in R xx produced by the reduction process in the GdO x device at various temperatures. Here, all the measurements were performed in a vacuum chamber under a 10mbar pressure and 0.7% humidity.
resistance. By applying an appropriate electric field at room temperature, we achieved two different oxygen ion concentration states with resistivity, ρ xx , values of 2200 and 1700 μΩ·cm. Figure 4a shows the temperature-dependent resistivity for two oxygen-concentration-controlled samples (red and blue symbols) and metallic Gd (green symbol). The temperature coefficient of resistivity, dρ xx dT , was found to be positive for metallic Gd (green symbols), while it was negative for the oxygen-concentration-controlled samples (red and blue symbols). This indicates that oxygen ion inclusions make the system exhibit a more insulator-like state.
The MR was then measured at 10 K, which is far lower than the Curie temperature of GdO x (for a discussion about the Curie temperature, see Supplementary Information 7). Figure 4b shows that all three samples exhibit a negative MR, that is, the resistance decreases as the magnetic field increases 23 . In metallic ferromagnets, a negative MR is usually attributed to suppression of magnons (or spin fluctuations) by an external magnetic field [23][24][25] . However, here, we found that the slope of the negative MR could be remarkably enhanced by increasing the oxygen ion concentration. The variation in MR reached 20% at 9 T for the ρ xx = 2200 μΩ·cm sample.
The large enhancement in MR observed in Fig. 4b can be explained by the tunneling giant MR (tunneling GMR) that generally appears in inhomogeneous granular films [26][27][28] . In magnetically inhomogeneous media with nonaligned ferromagnetic entities, it has been found that the resistivity can be gradually decreased by applying a magnetic field since the magnetic field aligns the ferromagnetic entities parallel to each other. This leads to a low resistance state in tunneling transport.
The tunneling model can be further supported by percolation theory. According to percolation theory, metallic conduction and hopping (or tunneling) conduction Longitudinal resistance (R xx ) with repeated switching of the current (I) direction under a current of a 1mA (E=0.55×10 5 V/m), b 2mA (E=1.1×10 5 V/m), or c 3mA (E=1.6×10 5 V/m) in the SiN x (3nm)/GdO x (10nm) device. When R xx <R LOWER =3.00kΩ or R xx >R HIGH =3.06kΩ was met, the current direction (electric field polarity) was switched. Fig. 4 Effect of ion migration on the magnetoresistance. a Temperature-dependent resistivity for oxygen-concentrationcontrolled GdO x (red and blue symbols) and metallic Gd (green symbol). b Magnetic-field-dependent resistivity of the three samples at 10K. The colors correspond to those in a. c Relationship between 1/ρ xx and T −1/2 (the data are replotted from a). The solid lines represent fitting curves based on Eq. (1). coexist in partially oxidized systems. Then, the temperature dependence of the resistivity can be given by [28][29][30] where c is the activation energy for electron tunneling, ρ t is the resistivity of the tunneling channel, and ρ c is the resistivity of the metallic channel. Figure 4c shows ρ À1 xx versus T −1/2 for two insulator-like samples (replotted from Fig. 4a). The curve shows the best fitting based on Eq. (1). The results indicate that ρ t = 2900 μΩ·cm and ρ c = 4200 μΩ·cm for the red curve and that ρ t, = 2700 μΩ·cm and ρ c = 2900 μΩ·cm for the blue curve. This implies that the contribution of tunneling conduction is comparable to or even more dominant than that of metallic conduction in our GdO x , which verifies the tunneling-induced MR. Last, we would like to emphasize that unlike previous reports where the tunneling GMR was controlled by the material composition [26][27][28] , in our GdO x device, the MR can be easily controlled by tuning the oxygen ion migration using an electric field.

Conclusions
In summary, we investigated lateral-electric-field-driven oxygen ion migration in a partially oxidized Gd device. Using an optical microscope, we directly observed the directional motion of oxygen ions driven by an electric field. The resistance was found to be closely linked to the oxygen ion concentration, allowing the electrical resistance to be greatly modulated by an electric field. We also investigated the magnetic properties of Gd and found that oxygen ions greatly enhanced the magnetoresistance of Gd. Our results therefore suggest that the electrical resistance and magnetoresistance of Gd can be significantly modulated by electric-field-driven oxygen ion migration, which may enable functional electromagnetic devices such as artificial synapses.