Regulating oxygen ion transport at the nanoscale to enable highly cyclable magneto-ionic control of magnetism Zhengwei Tan,† Zheng Ma,† Laura Fuentes,‡, # Maciej Oskar. Liedke,¥ Maik Butterling,¥ Ahmed G. Attallah,¥ Eric Hirschmann,¥ Andreas Wagner,¥ Llibertat Abad,‡,# Nieves Casañ-Pastor,‡ Aitor F. Lopeandia,† Enric Menéndez,*, †, Jordi Sort*,†,⁋ † Departament de Física, Universitat Autònoma de Barcelona, 08193 Cerdanyola del Vallès, Spain ‡ Institut de Ciència de Materials de Barcelona, CSIC, Campus UAB, 08193 Bellaterra, Barcelona, Spain # Centre Nacional de Microelectrònica, Institut de Microelectrònica de Barcelona-CSIC, Campus UAB, 08193 Bellaterra, Barcelona, Spain ¥ Institute of Radiation Physics, Helmholtz-Zentrum Dresden - Rossendorf, Dresden 01328, Germany ⁋ Institució Catalana de Recerca i Estudis Avançats (ICREA), Pg. Lluís Companys 23, E-08010 Barcelona, Spain 1 KEYWORDS: magnetoelectricity, voltage control of magnetism, magneto-ionics, transition metal oxide, ion diffusion ABSTRACT Magneto-ionics refers to the control of magnetic properties of materials through voltage-driven ion motion. To generate effective electric fields, either solid or liquid electrolytes are utilized, which also serve as ion reservoirs. Thin solid electrolytes have difficulties to (i) withstand high electric fields without electric pinholes and (ii) maintain stable ion transport during long-term actuation. In turn, the use of liquid electrolytes can result in poor cyclability, thus limiting their applicability. Here we propose a nanoscale-engineered magneto-ionic architecture (comprising a thin solid electrolyte in contact with a liquid electrolyte), that drastically enhances cyclability while preserving sufficiently high electric fields to trigger ion motion. Specifically, we show that the insertion of a highly nanostructured (amorphous-like) Ta layer (with suitable thickness and electric resistivity) between a magneto-ionic target material (i.e., Co3O4) and the liquid electrolyte, increases magneto-ionic cyclability from < 30 cycles (when no Ta is inserted) to more than 800 cycles. Transmission electron microscopy together with variable energy positron annihilation spectroscopy reveal the crucial role of the generated TaOx interlayer as a solid-electrolyte (i.e., ionic conductor) that improves magneto-ionic endurance by proper tuning of the types of voltage- driven structural defects. The Ta layer is very effective in trapping oxygen and hindering O2– ions from moving into the liquid electrolyte, thus keeping O2– motion mainly restricted between Co3O4 and Ta when voltage of alternating polarity is applied. We demonstrate that this approach provides a suitable strategy to boost magneto-ionics by combining the benefits of solid and liquid electrolytes in a synergetic manner. 2 With the advent of global phenomena such as the Internet of Things, artificial intelligence, machine learning or Big Data, the demand for highly functional and energy-efficient miniaturized microelectronic components is growing exponentially.1,2 Spintronic systems3,4 operated using electric currents through spin-transfer torque5,6 or spin-orbit torque7,8 effects are becoming key elements for next-generation nanoelectronics with enhanced memory and information processing capabilities. However, magnetization switching using electric current involves undesirable Joule heating effect, which is detrimental to energy efficiency.9,10 An interesting alternative is to modulate magnetic properties by applying electric fields instead of electric current, thus minimizing power dissipation. This has rapidly evolved into a whole area of research referred to as voltage control of magnetism (VCM). Magneto-ionics refers to a particular mechanism for VCM in which voltage-driven ion transport (of e.g., O2−,11–15 Li+,16 F−,17 H+,18–21 or N3− 22–26 species) leads to a large and controllable modulation of magnetism without the need of strain transfer. This is different from voltage- controlled strain-mediated multiferroic heterostructures, which are less convenient for spintronics because repeated voltage actuation can lead to mechanical fatigue and eventual device failure. Archetypical magneto-ionic structures comprise a gate electrolyte (either solid or liquid)27,28 in contact with a ferromagnetic (FM) target material. Voltage is applied across the electrolyte, using the FM film as working electrode and causing a voltage polarity-dependent insertion/removal of ions into/from the target material. In this way, magnetic properties, such as saturation magnetization,11,22 magnetic anisotropy and coercivity,12,18 exchange bias field,13,29 or skyrmion generation/suppression,20 amongst others, can be reversibly controlled, making magneto-ionic materials highly promising for ultralow power magnetic devices. An extreme case is when voltage induces a complete reversible transition between FM and non-magnetic states, leading to voltage- 3 driven ON-OFF switching of ferromagnetism. This has been reported in a number of magneto- ionic systems, such as Co3O4,11,14,30 Li+ intercalated -Fe2O3,31 F− intercalated La2−2xSr1+2xMn2O7,17 CoN26,28, CoMnN25 or -Co(OH)2.32 ON-OFF switching of the ferromagnetic state has been also induced by electrostatic surface charging in FeS2.33 In most magneto-ionic systems, the source of the moving ions are electrolytes (e.g., GdOx, HfOx, propylene carbonate with dissolved LiPF6 or KI, etc.) that are in direct contact with pristine FM or ferrimagnetic layers (e.g., Co, Fe, Fe2O3), whose properties are manipulated with voltage.27,34– 36 Magneto-ionic systems based on O2− insertion from an external electrolyte to a target FM material often suffer from slow dynamics at room temperature and irreversible compositional/ structural changes in the FM phase, eventually leading to degradation and limited cyclability.12 In addition, for relatively thick films, it is challenging to achieve a fully OFF magnetic state with voltage due to the limited penetration of the external ions towards the interior of the FM target layers.37 Recently, smaller ions (e.g., H+, F− or Li+) have been introduced to achieve faster and more cyclable voltage-induced manipulation of magnetic properties. Fast, reversible and cyclable tuning of perpendicular magnetic anisotropy,18 Dzyaloshinskii-Moriya interaction38,39 or ferrimagnetic spin textures40 has been achieved through chemisorbed O2− or H+ ion species. However, systems relying on H+ are either sensitive to environmental conditions (e.g., humidity) or restricted by their incompatibility with traditional complementary metal-oxide semiconductor (CMOS)-based devices (i.e., standard fabrication processes for semiconductor devices, such as metal oxide field-effect transistors).28 An alternative approach is to use structural oxygen or nitrogen, self-contained in the target materials, as the source of the moving ions. Examples of such materials are Co3O4,11 CoN22,26 or CoMnN25 films. These target materials, which are CMOS compatible, exhibit an initially fully OFF (i.e., paramagnetic) state and provide “ready-prepared” 4 lattice sites for ion diffusion, allowing for net magneto-ionic generation of ferromagnetism by voltage-triggered O2− or N3− ion motion from the films towards a neighboring electrolyte. Unfortunately, achieving a high magneto-ionic cyclability in these materials (i.e., removing and reinserting the O2−/N3− ions many times by switching voltage polarity) remains a challenge. To induce magneto-ionics, either solid or liquid electrolytes can be utilized. Solid electrolytes (with ultra-thin dielectric layers) are preferred for solid-state spintronics. Ultra-thin solid dielectric layers are needed to induce sufficiently large electric fields under moderate applied voltages. However, at such small thicknesses, difficulties arise to withstand high electric fields without electric pinholes. Moreover, thin solid electrolytes offer a limited ion buffering capability (i.e., they easily become saturated with ions and cannot sustain stable ion transport during long-term operation, especially when the magneto-ionic layer is thicker than the solid electrolyte layer). Liquid electrolytes are convenient for other magnetoelectric applications, such as magnetophoresis/microfluidics41 or to emulate neuromorphic functionalities (the brain operates in a liquid environment).15,26 Owing to the formation of the “electric double layer”, whose thickness is < 1 nm, liquid electrolytes are able to generate ultra-large electric fields (hundreds of MV cm–1) at the interface between the liquid and the target magneto-ionic layer, without electric pinholes.27 Liquid electrolytes may also be good ion reservoirs.11,14,22 However, when voltage is applied, the mobile ions released from the target layer into the liquid can travel long distances towards the counter electrode, making their reintroduction into the magneto-ionic layer (with voltage of opposite polarity) difficult. In addition, at the counter electrode, if sufficiently high voltage is applied, the dissolved ions (e.g., O2−, H+ or N3−) can transform into the corresponding oxygen, hydrogen or nitrogen gases and be released to the atmosphere in the form of bubbles. These effects are difficult to control and are highly detrimental for magneto-ionic reversibility and endurance. 5 Restricting ion transport within the magneto-ionic layers could avoid these problems and is expected to improve the magneto-ionic cyclability. Here, we propose an improved nanoscale-engineered magneto-ionic structure that results from inserting an amorphous Ta layer (which gets spontaneously passivated in air) between the Co3O4 film (magneto-ionic target material) and propylene carbonate, PC (liquid electrolyte). Upon application of negative voltage, O2– ions migrate from Co3O4 to Ta, promoting the formation of TaOx (Ta is a good oxygen getter), which acts as a thin solid electrolyte with good ionic conductivity. This architecture preserves sufficiently high electric fields to trigger ion motion while allowing to repeatedly induce ON-OFF switching of ferromagnetism at room temperature. By optimizing the thickness of Ta layer, a drastic increase of cyclability is achieved in Co3O4/Ta (25 nm)/PC (> 800 cycles) compared to Co3O4/PC with no Ta insertion (< 30 cycles). The Ta layer hinders oxygen ions from entering the liquid electrolyte and allows for O2– redistribution inside the Co3O4 and Ta layers, as assessed by high-angle annular dark-field scanning transmission electron microscopy (HAADF-STEM) and electron energy loss spectroscopy (EELS). Positron annihilation spectroscopy is used to precisely characterize the defect structure in Ta, which strongly contributes to the O2– ion transport. An increase of the defect size during biasing is observed. The reduced electric conductivity of the passivated amorphous Ta layer is important to allow the penetration of electric field inside the Co3O4 layer and the concomitant O2– ion transport. Since O2– diffusion becomes mainly restricted within the Co3O4/TaOx structure (instead of O2– being released to PC), this greatly improves the efficiency of ion transport in the material, resulting in largely enhanced magneto-ionic cyclability. 6 RESULTS AND DISCUSSION The basic building block of the investigated magneto-ionic system is a 40 nm thick Co3O4 film grown by DC reactive sputtering at room temperature onto Cu (60 nm)/Ti (20 nm)/[100]Si (725 μm) substrates. The Co3O4 films were coated with sputtered amorphous Ta protective layers of variable thickness (from 5 to 50 nm) which were left to passivate in air (see Methods). Uncoated Co3O4 films were also grown as a reference. Figure 1(a) shows a low magnification high angle annular dark field - scanning transmission electron microscopy (HAADF-STEM) image of the as- prepared sample with a 25 nm-thick Ta protective capping layer. Clear interfaces between the various layers can be seen, and their thicknesses are in good agreement with the nominal ones. Note that the topmost Pt capping layer was grown only during the TEM lamella preparation (i.e., it was not present during magneto-ionic experiments). To investigate the Co, O and Ta distribution, Co, O and Ta electron energy loss spectroscopy (EELS) mappings from the area marked with an orange rectangle in Figure 1(a) were acquired. As shown in Figure 1(b), Co and O elements are homogeneously distributed within the Co3O4 films. Additionally, some O signal is also detected inside the Ta layer, particularly on its upper part. This is ascribed to surface self-passivation of Ta in air, which is known to result in the formation of a stable TaOx top layer.42,43 Figure 1(c) shows the θ/2θ X-ray diffraction (XRD) pattern of an as-grown film. The XRD peak observed at 2θ ≈ 19.1° is consistent with (111) Co3O4 (PDF® 00-009-0418), whereas that at 2θ ≈ 36.7° could be ascribed to either (111) Co3O4 or (111) CoO (PDF® 00-001-1227). Peaks corresponding to the Cu buffer layer are also detected. Conversely, no traces of Ta are evidenced, suggesting that Ta is highly nanostructure or even amorphous-like. Growth of amorphous Ta by physical deposition methods has been previously reported44,45 and it is promoted by the existence of ions and particles with high kinetic energy when sputtering at sufficiently high gun powers. Amorphous Ta exhibits 7 higher electric resistivity than its crystalline counterparts.44 This is beneficial to keep sufficient electric field strength inside the actuated films while voltage is applied, thereby enabling magneto- ionics. To further investigate the structure of the samples, high-resolution transmission electron microscopy (HRTEM) images of the cross-sections of as-prepared Co3O4 and Ta layers were taken, as shown in Figure 1(d) and 1(e), respectively. The areas marked with red squares were chosen for fast Fourier transform (FFT) analyses, as shown in the insets. In Figure 1(d), “1” spot corresponds to an interplanar distance of 0.248 nm and is consistent with (311) Co3O4 (PDF® 00-009-0418) or (111) CoO (PDF® 00-001-1227) interplanar distances (0.244 and 0.245 nm, respectively), while spot “2” corresponds to an interplanar distance of 0.225 nm and is unambiguously ascribed to (200) CoO (PDF® 00-001-1227) interplanar distance (0.212 nm). “3” spot corresponds to an interplanar distance of 0.448 nm, which is consistent with (311) Co3O4 (0.467 nm, PDF® 00-009-0418). This indicates that a mixture of CoO and Co3O4 phases is plausible in the Co oxide film, in agreement with previously reported results.15 For simplicity, in spite of the presence of CoO in particular close to the interface, the Co oxide film is labelled Co3O4 throughout the manuscript. Remarkably, no spots are observed for Ta, in agreement with its amorphous nature, as also evidenced by XRD (Figure 1(c)). 8 Figure 1. Structural and compositional characterization of as-deposited films. (a) HAADF-STEM micrograph of an as-grown sample. (b) EELS Co (in red), O (in blue) and Ta (in yellow) elemental mappings of the area marked with an orange rectangle in (a). (c) θ/2θ XRD patterns of the as- prepared samples. (d), (e) High-resolution TEM images of the cross section of the as-deposited Co3O4 and Ta films, respectively. The inset shows the fast Fourier transform of the areas marked with red squares. In panel e, “1” spot corresponds to an interplanar distance of 0.248 nm and is consistent with (311) Co3O4 (PDF® 00-009-0418) or (111) CoO (PDF® 00-001-1227) interplanar distances (0.244 and 0.245 nm, respectively), while spot “2” corresponds to an interplanar distance of 0.225 nm and is unambiguously ascribed to (200) CoO (PDF® 00-001-1227) interplanar distance (0.212 nm). “3” spot corresponds to an interplanar distance of 0.448 nm, which is consistent with (311) Co3O4 (0.467 nm, PDF® 00-009-0418). To induce magneto-ionics, electrolyte-gating was performed in a capacitor-like configuration (Figure 2(a))14 using a platinum wire as counter electrode and an aprotic, anhydrous polar liquid electrolyte composed of propylene carbonate with Na+ and OH– solvated species.27,46,47 When 9 voltage is applied, a sub-nm-thick electric double layer forms at the electrolyte side of the electrolyte/Ta interface, allowing for the generation of a high electric field.47 This electric field is ultimately responsible for driving oxygen ions from the Co3O4 to the Ta/TaOx layers (as illustrated in Figure 2(b)). Voltage treatments were performed in-situ, while hysteresis loops were recorded at room temperature by vibrating sample magnetometry (VSM), with an in-plane applied magnetic field. The total measured magnetic moment of the samples (in emu) was normalized to the area of the sample and the nominal thickness of Co3O4 (to obtain emu cm–3). As seen in Figure 3, irrespective of the thickness of the Ta layer (i.e., 0, 5, 10, 25 or 50 nm), all films in the as-grown state show very little ferromagnetic response (< 33 emu cm–3, which is equivalent to ≈ 2.3% the magnetization of pure FCC-Co48). This small ferromagnetic signal is common in sputtered Co3O4 15 and is likely due to either a small fraction of residual Co clusters that do not become fully oxidized during the sputtering process or substrate contamination. The deposition of Ta onto Co3O4 has a negligible effect on the initial ferromagnetic signal of the Co3O4 films. Upon negative biasing at –25 V for 30 min or 1 h, clear hysteresis loops are observed. The generated saturation magnetization (MS) is maximum for the sample without Ta (around 400 emu cm–3, suggesting that 27% of the volume of Co3O4 gets reduced to Co) and it slightly decreases when a thin Ta interlayer (5 – 10 nm) is inserted between the Co3O4 film and the liquid electrolyte. MS further decreases for sufficiently thick Ta. Remarkably, magneto-ionic effects when Pt (instead of Ta) is grown onto Co3O4 are strongly reduced: a small hysteresis loop is observed only for 5 nm Pt (compare Figure 3(f) with 3(b)), and no ferromagnetic response is obtained for thicker Pt. In all cases, the initial virtually non-magnetic state can be recovered by applying +25 V for 30 min. 10 Figure 2. Magneto-ionic characterization of Co3O4 thin films and Co3O4/Ta bilayer films under electrolyte gating. (a) Schematic of the designed structure for electrolyte actuation. (b) Left: Sketch of the formation of electric double layer during voltage actuation; Right: photograph of the home- 11 made electrolytic cell used to apply voltage to our system, which consists of an Eppendorf filled with propylene carbonate, our sample as working electrode and a Pt wire counter electrode. (c) Time evolution of the saturation magnetization, MS, for Co3O4 thin films and Co3O4/Ta bilayer films (with Ta thickness of 5, 10, 25 and 50 nm) under applied voltage of –25 V. (d) Estimation of the onset voltage required to trigger magneto-ionics in Co3O4 thin films and Co3O4/Ta (25 nm) bilayer films, as well as the recovery process. (e) Calculated absolute value of the voltage at the Co3O4/Ta interface, |Vint|, as a function of the TaOx conductivity, TaOx, when externally applying |V| = 25 V for different thickness of Ta layer, δTa, and a fixed Co3O4 conductivity, Co3O4 = 20 S m–1. (f) Calculated interface voltage as a function of TaOx, for different values of conductivity of Co3O4 and a fixed Ta thickness (δTa= 50 nm). The regions highlighted with a cyan rectangle in (e) and (f) are those where |Vint| would be below the magneto-ionic threshold voltage of the system. Figure 3 shows that, for all samples, most of the ferromagnetic response is induced in less than 30 min, since MS does not further increase when applying voltage for 1 h. The time evolution of MS for all samples under –25 V is shown in Figure 2(c). An external magnetic field of 10 kOe (i.e., above the anisotropy field of the generated ferromagnetic counterpart) was applied during these measurements to ensure magnetic saturation. In all films, an immediate increase of MS is observed in response to the applied V, evidencing a quick onset of the oxygen ionic motion, which leaves metallic ferromagnetic Co behind.14,15 The obtained steady-state value of MS decreases with the Ta thickness, from 401 emu cm–3 to 19 emu cm–3 (in agreement with the magnetic hysteresis loops shown in Figure 3(a)–(e)). One important parameter of large significance for device applications is the minimum threshold voltage required to trigger magneto-ionic effects.22,25 The onset voltage for Co3O4 films without and with a 25 nm–thick Ta capping layer was evaluated by subjecting the 12 films to increasing negative voltage steps of –2 V, until MS started to increase, as shown in Figure 2(d). The results reveal that the onset voltage is approximately –6 V for both samples, whereas a voltage of +8 V leads to complete recovery in both cases, which agrees with previous works on similar systems.11 Figure 3. In-plane VSM hysteresis loops (each lasting 30 min) of the as-prepared Co3O4/Ta films (black), the films biased under –25 V (for 30 min in red, and for 1h in blue), and subsequently 13 recovered after applying +25 V for 30 min. The different panels correspond to different Ta thicknesses: (a) 0 nm, (b) 5 nm, (c) 10 nm, (d) 25 nm, (e) 50 nm. Panel (f) corresponds to Co3O4/Pt without Ta interlayer (with Pt thickness of 5 nm). Figures 2(c) and 2(d) also show that the rate at which MS increases under voltage application is larger for smaller Ta layer thickness. This, together with the reduction of the steady-state MS for thicker Ta, indicates that the effective electric field acting on the Co3O4 layer becomes lower for larger Ta thickness. This suggests that, for thicker Ta, there is a more severe dissipation of space charge density while O2– is migrating from Co3O4 to form TaOx.49 In a first approximation, this effect can be modelled considering the system as a simple voltage divider. The voltage drop at the Co3O4/Ta interface can be estimated assuming that, during the magneto-ionic process, the Co3O4 and the newly formed TaOx layers act as two resistances connected in series. In this case, the voltage at the Co3O4/TaOx interface is: 𝑅𝐶𝑜3𝑂4 ∆𝑉𝑖𝑛𝑡 = ∆𝑉. Considering that R =  δ/A (where  is the resistivity of each material, 𝑅𝐶𝑜3𝑂4 +𝑅𝑇𝑎𝑂𝑥 δ is the film thickness and A is the lateral area of the films), it is straightforward to estimate ∆𝑉𝑖𝑛𝑡 as a function of the layers’ thicknesses and the electric conductivities of Co3O4 and TaOx. The electrical conductivity of Co3O4, measured by the Van der Pauw method, is approximately Co3O4 = 20 S m–1.14 The conductivity of TaOx strongly depends on the oxygen content and it varies by several orders of magnitude, from 105 S m–1 in amorphous metallic Ta44 to 10–6 S m–1 for highly oxidized Ta.50 The presence of oxygen in the naturally passivated Ta (Figure 1(b)) can easily bring  to values in the range 10–1 – 102 S m–1 before any voltage is applied.51 Figure 2(e) shows the evolution of the calculated interface voltage as a function of TaOx for the different Ta film thicknesses. For highly conductive Ta (or TaOx) spacer layers, ∆𝑉𝑖𝑛𝑡 becomes independent of the Ta thickness, which means that any eventual magneto-ionic effect would be independent of the Ta 14 thickness (and ∆𝑉𝑖𝑛𝑡 would be always equal to –25 V). However, in this case, electric field would be highly screened at the interface and would not penetrate inside Ta (since, in metals, electric field is confined within the Thomas-Fermi screening length, which is typically < 0.5 nm).27 Thus, for highly metallic interlayers (e.g., noble-metal like Pt, as in Figure 3(f)), little magneto-ionic effects are expected since there is no driving force (no electric field) to induce ion motion inside the metal.24 This is opposite to what happens in semiconductors (e.g., TaOx), where electric field will penetrate deeper into the layer. Interestingly, for low electrically conductive TaOx (i.e., TaOx < 100 S m–1), ∆𝑉𝑖𝑛𝑡 becomes clearly lower as the Ta layer thickness increases (Figure 2(e)). For δTa = 50 nm, if TaOx is lower than 5 S m–1, ∆𝑉𝑖𝑛𝑡 eventually drops below the threshold voltage (region highlighted in cyan), meaning that no magneto-ionic effects will be induced in this case (while, for thinner Ta layers, for the same TaOx value, ∆𝑉𝑖 will still be above the threshold). Note that during magneto-ionic motion of O2– from Co3O4 to TaOx, the values of TaOx will rapidly decrease, thus reducing the interface voltage and slowing down the magneto-ionic process. Another interesting effect is that when Co3O4 transforms to metallic Co, the conductivity of the magneto-ionic layer (Co3O4) increases. This, in turn, has an effect on ∆𝑉𝑖𝑛𝑡 , which is shown in Figure 2 (f). In bulk metallic Co, Co can reach values >107 S m–1.52 As plotted in Figure 2(f), when Co3O4 transforms to Co and Co3O4 increases, the interface voltage drastically drops, and ∆𝑉𝑖𝑚𝑡 easily falls below the threshold voltage. Then any magneto-ionic response will tend to stop. As an example, for the particular case of δTa = 50 nm, any eventual O2– ion motion triggered while TaOx is sufficiently large, will tend to stop as soon as Co3O4 increases above 103 S m–1 (see Figure 2(f)) or TaOx decreases below 5 S m–1 (Figure 2(e)). It is noteworthy that, although this intuitive picture is a simplified representation of the reality, it already provides some basic understanding of the role of the TaOx layer thickness and resistivity on the induced magneto-ionic effects. 15 Notwithstanding the decrease of the steady-state MS with the deposition of a Ta capping layer (by approximately 25% for δTa = 25 nm), the formation of TaOx (which acts as a solid electrolyte) has one very beneficial effect: it drastically enhances magneto-ionic cyclability. This is evidenced in Figure 4, which shows the cyclability results (i.e., repeated increase/decrease of M) in Co3O4 thin film and Co3O4/Ta (25 nm) bilayer films upon application of –25V/+25V voltage pulses with a duration of 80 s. As shown in Figs. 4(a) and 4(b), a very stable and reversible behavior is observed during the first 10 cycles for both samples. However, the Co3O4/Ta (25 nm) bilayer sample maintains a very stable cyclability even after 800 cycles (Figure 4(a)), whereas no sign of magneto- ionic effect is detected for Co3O4 films after 30 cycles (Figure 4(b)). This demonstrates a significant enhancement of the endurance of the system by the deposition of a Ta overlayer with appropriate thickness. Figure 4. Magneto-ionic cyclability of (a) Co3O4/Ta (25 nm) bilayer films and (b) the Co3O4 thin films subjected to –25 V / +25 V voltage pulses applied with a periodicity of 80 s. The improved cyclability in the Co3O4/Ta (25 nm) bilayer is ascribed to the role of Ta in allowing diffusion of oxygen ions. In this sample, when negative voltage is applied, the O2– ions exiting 16 Co3O4 are captured by Ta, forming a TaOx solid electrolyte, instead of being directly released to the liquid electrolyte. Ta is a well-known oxygen getter53 and TaOx is a good ionic conductor.54 In unprotected Co3O4 films, O2– ions are directly dissolved in the propylene carbonate and once solvated with the PC chains, they can travel long distances towards the positively charged counter electrode (i.e., the Pt wire), where they eventually form O2 gas bubbles and are released to the atmosphere at sufficiently high voltages. This long-distance transport of ions in the liquid is the main reason for the poor cyclability of uncapped Co3O4 films. However, when a Ta interlayer is grown adjacent to Co3O4, O2– ions move back and forth along relatively short distances (for alternating applied voltages of opposite polarities), and the process is confined between Co3O4 and the TaOx solid electrolyte, thus drastically improving cyclability. 17 Figure 5. Compositional characterization of the Co3O4/Ta (25 nm) films gated at −25 V for 1 h. (a) HAADF-STEM and (b–c) elemental EELS mappings corresponding to the area marked with an orange rectangle in the HAADF-STEM image. Cu and Pt layers serve as working electrode and protective capping layer during TEM lamellae preparation, respectively. Co, O and Ta are represented by red, blue and yellow colors in the EELS elemental mappings. (d) Depth profile of Co, O and Ta elements along the dark pink arrow drawn in (a). (e) High-resolution TEM image of the sample treated at −25 V for 1 h. To further understand the effect of the top Ta layer on the induced magnetic properties from the perspective of ion transport mechanism, cross-sectional lamellae of a bilayer sample (40 nm–thick Co3O4 plus 25 nm–thick Ta overlayer) electrolyte-gated at –25 V for 1 h were studied by HAADF- 18 STEM and EELS. As shown in Figure 5(a) – (c), an inhomogeneous microstructure inside the Co3O4 layer is generated upon voltage application, with Co-rich and O-rich regions, leading to void-like morphologies, in agreement with previous studies on this kind of magneto-ionic material.11 In turn, an enrichment in O is observed in the upper part of the Ta layer. This is clearly evidenced by quantitative EELS analysis (Figure 5(d)) and high-resolution TEM (Figure 5(e)). This indicates that, besides the O2– ions transport from Co3O4 towards and across the Ta layer, O2– ions also locally migrate and redistribute inside the Co3O4 layer, eventually forming metallic Co.11 Interestingly, the mixing of O and Ta elements on the top part of the Ta layer demonstrates the oxidation of Ta. Considering that the liquid electrolyte used in this work provides a non-aqueous environment, the observed O signal can only originate from oxygen transport from Co3O4. The accumulation of O2– on the top of the TaOx layer suggests that ion transport remains rather restricted within the Co3O4/TaOx bilayers (i.e., short-distance ion diffusion). In other words, the presence of the thin TaOx solid electrolyte, limits the amount of O2– released to the liquid, thereby enhancing cyclability. Furthermore, during this magneto-ionic experiment, Co3O4 tends to become more amorphous. This is corroborated by θ/2θ X-ray diffraction (see Figure 6(a)), where the peak from Co3O4 (111) planes disappears upon voltage treatment. In addition, the fast Fourier transform (FFT) spots obtained from high-resolution TEM prove the existence of metallic Co in the treated sample, as can be seen in Figure 6(b). The formation of metallic Co is responsible for the MS increase after voltage treatment (Figure 2(c) and Figure 3). To assess the electrochemical behavior of the system, cyclic voltammetry (CV) curves were recorded for all the surfaces involved in the main magneto-ionic setup (see Fig. S1 in the Supporting Information). When the Cu grown on Ti is exposed to the electrolyte, a cathodic wide peak is observed, centred at –0.9 V vs. Pt reference electrode, while further anodic run evidences 19 a flat peak. This is consistent, in principle, with the reduction of Cu2O (native oxide layer formed on Cu when exposed to air) to Cu (first wave) and consequent oxidation of Cu to Cu2O. The CV curve of Co3O4 grown onto the Ti/Cu does not show the Cu waves, confirming that the cobalt oxide layer does cover completely the Cu phase and prevents its oxidation/reduction. Once Ta is grown on Co3O4, a tiny wide oxidation wave is observed, centered at 0.6 V vs. Pt, in agreement with the expected oxidation of Ta. Since both coatings (Co3O4 and TaOx) are less conductive than Cu, the redox processes are hindered compared to the case of Cu. Figure 6. Structural characterization by X-ray diffraction (XRD) and high-resolution transmission electron microscopy (HRTEM). (a) θ/2θ XRD diffraction patterns of Co3O4/Ta (25 nm) films in the as-prepared state and after being treated at –25 V for 1 h. (b) HRTEM images of the cross section of Co3O4/Ta (25 nm) film after applying –25 V for 1 h. The inset shows the fast Fourier transform (FFT) of the area marked with a red rectangle, and it shows spots with an average interplanar distance of 0.19 nm, which corresponds to (101) hexagonal-closed-packed (HCP) metallic Co. For phase identification, the cards no. PDF® 00-005-0727 and PDF® 00-009-0418 were taken for Co and Co3O4, respectively. 20 To shed further light on the microstructure of the films upon magneto-ionic actuation, Doppler broadening variable energy positron annihilation spectroscopy (DB-VEPAS)55 and variable energy positron annihilation lifetime spectroscopy (VEPALS)56 experiments were performed (see Methods). PAS is a sensitive probe to open volume defects on atomic scales and across depth due to positrons preferential localization and subsequent annihilation with electrons in crystal empty spaces, i.e., vacancies and their agglomerates. The energy distribution and time necessary for annihilation depends on the local electron density.57 As shown in Figure 7(a), after negative biasing at –25 V, the low electron momentum fraction, S (directly proportional to the size and concentration of defects),58 increases in the Ta region but slightly drops in the Co3O4 layer, suggesting a substantial raise of defect density in Ta and a slight drop of defect concentration in Co3O4. The increase of defect density in Ta is probably a consequence of the O2– interdiffusion and the formation of TaOx. Additionally, the average defect size (which is proportional to the average positron lifetime, τav) strongly increases in both TaOx and Co3O4 layers after voltage actuation (Figure 7(b)). The expected average size is in the range of large vacancy agglomerates (8–10 vacancies).59 Figures 7(b) and 7(c) show the existence of three discrete lifetime components (τ1, τ2 and τ3),59 which correspond to three different average defect sizes.11 They were obtained by deconvolution of PALS spectra using a nonlinear least-squared fitting method (package PALSfit software).60 The corresponding relative intensities (I1, I2 and I3) relate to the concentration of each defect type (size). The larger τi is, the larger the defect size is since it takes longer for positrons to be annihilated with electrons.57 In the as-grown sample, only τ1 and τ2 lifetime components are detected, meaning an absence of large void-like structures (no τ3) in the pre-biasing state. The most abundant defect size (τ1) contributes as I1 = 62 – 77% of the overall signal and is in the range of 21 bi-vacancy clusters in the upper sub-surface Ta region (i.e., passivation layer) and single vacancy in the underneath Ta sub-layer.59 τ1 ≈ 213 ps for Co3O4 represents a defect size involving 3 – 4 mixed (Co and O) vacancies within a complex, based on our previous works/publications.11 After biasing at –25V, τ1 strongly increases both in Ta and Co3O4 layers. Open volume is generated in the size range of 4-6 vacancy agglomerations for Ta and > 7 – 8 for Co3O4, which result in an increase of intensity to I1 = 75 – 95%. The second lifetime component, τ2, for the as-grown sample is in the range of large vacancy agglomerations (>10 vacancies), typical for grain boundaries, with corresponding I2 = 28 – 37% of positrons being annihilated at these defect states. Negative biasing increases τ2 to the range of voids (the threshold is usually about 500 ps) with an average diameter in the range of 0.28 – 0.37 nm,61 while their density is quite small, up to I2 ≈ 10% and they do not reach deeper than the upper 10 nm of Ta. Clearly, these defects are associated with the formation of the upper TaOx (Figure 5). Finally, a larger pores population (0.7–0.8 nm in diameter) was found after biasing (τ3), with low but not negligible intensity (I3 ≈ 1.4%), which extends into the Co3O4 layer. The existence of pore related components (especially τ3) is a fingerprint of increasing amorphization of the TaOx layer whereas the Co3O4 layer remains less affected and nanocrystalline. The largest free volume is found in the direct vicinity of the surface, i.e., the most amorphized region. The increase of positron lifetimes in the Ta layer after biasing evidences the O2– transport generated expansion of available open volume channels, which enables large cyclability of the system. 22 Figure 7. Positron annihilation spectroscopy characterization of the as-grown and voltage-treated (–25 V for 1 h) Co3O4/Ta (25 nm) sample. (a) Dependence of the low electron momentum fraction (S) on the positron implantation energy (Ep). (b) Dependence of the average positron lifetime, τav, on Ep. (c) Dependence of the positron lifetime components (τi) on Ep. (d) Dependence of the relative defect type intensity (Ii) on Ep. CONCLUSIONS In summary, this work demonstrates the beneficial effect of adding a thin solid ionic conductor (in this case, a naturally passivated amorphous Ta interlayer) on the magneto-ionic response of Co3O4 films under the action of electric field when immersed in liquid electrolytes. A priori, adding 23 a capping layer between the Co3O4 film and the liquid electrolyte could be thought of as simply hindering magneto-ionics and reducing (or eventually suppressing) any oxygen ion migration triggered by the externally applied electrical voltage. This is indeed observed for thin metallic capping layers (e.g., 5 nm Pt) or thick highly resistive Ta layers (with thickness larger than 50 nm). However, for thinner Ta interlayers, in spite of a moderate reduction of the steady-state Ms (by a factor 25% for 25 nm–thick Ta), a drastic enhancement of magneto-ionic cyclability is observed, from less than 30 cycles in uncoated Co3O4 to more than 800 cycles for Ta thickness of 25 nm. Such enhancement of endurance is ascribed to the key role of the generated TaOx layer in preventing O2– from being released to the liquid electrolyte, i.e., limiting oxygen ion transport within the Co3O4/TaOx layers. This is confirmed by compositional/structural characterization using HAADF-STEM and EELS as well as positron spectroscopy experiments. Beyond magneto- ionics, the reported strategy to enhance cyclability can be easily extrapolated to other systems relying on ion transport mechanisms, such as iontronics, sensors or neuromorphic computing. EXPERIMENTAL SECTION Sample Fabrication: 40 nm-thick Co3O4 thin films were grown at room temperature by reactive sputtering in a high-vacuum chamber (with a base pressure < 8 ×10–8 Torr) on non-doped 〈100〉– oriented Si wafers previously coated with a 20 nm–thick titanium adhesion layer and 60 nm-thick copper seed layer. Prior to growing Co3O4, the Cu seed layers were partly masked to leave enough space for the electric contact (i.e., to later serve as a working electrode). The Co3O4 films were subsequently coated with Ta protective layers of variable thickness, ranging from 5 to 50 nm, which were left unprotected to self-passivate in air. Two reference films were also prepared to serve as references: uncoated Co3O4, and Co3O4 coated with Pt with the same range of thicknesses 24 as Ta. The growth of Co3O4 was carried out in a mixed Ar and O2 atmosphere using an oxygen partial pressure of 50% and a total working pressure of 3×10–3 Torr. The distance between the substrate and targets was around 10 cm, and the growth rate was approximately 0.6 Å s–1. Ti, Cu, Ta and Pt layers were grown under 3×10–3 Torr Ar. The gun power to grow Ta was 100 W. Magneto-ionic Characterization and cyclic voltammetry curves: Room-temperature magnetoelectric measurements were performed using a commercial vibrating sample magnetometer (VSM) from Micro Sense (LOT, Quantum Design), with a maximum applied in- plane magnetic field of 2 T. The samples were electrolyte-gated using an external Agilent B2902A power supply, applying voltage between the counter electrode (a Pt wire) and the working electrode (e.g., the investigated Si/Ti/Cu/Co3O4/Ta thin films) in a home-made electrolytic cell (see Fig. 2(b)). The electrolyte consisted of anhydrous propylene carbonate with Na+ and OH– solvated species (10 – 25 ppm), formed by immersing small pieces of metallic sodium that were able to react with any possible traces of water.11 Negative voltages in this work indicate the accumulation of negative charges at the working electrode (and vice versa for positive voltages). The magnetization (M) was obtained by normalizing the magnetic moment to the sample volume exposed to the electrolyte. Note that the linear slopes in the hysteresis loops at high fields (arising from diamagnetic or paramagnetic contributions) were subtracted by correcting the background signal (i.e., at fields always significantly larger than the saturation fields). Cyclic voltammetry (CV) curves for all the surfaces involved in the main magneto-ionic setup were also recorded. For these experiments, a scan speed of 10 mV/s, and a potential sweep range from 0 V to –2 V to +1.4 V and to 0 V vs. Pt were selected. The exposed surface was 0.5 cm2 and the distance between the reference and working electrode was 4 mm, thus mimicking the conditions utilized during the magneto-ionic experiments. Here, however, we used 3 electrodes 25 (including a reference electrode) whereas magneto-ionic experiments we performed in a 2 electrodes configuration. Thus, all voltages in the CV curves are given with respect to a reference electrode, in this case Pt (99.99 % Goodfellow), known to act well in organic solvents. The global cell (two-electrode) potential has been measured during these CV experiments to be in the order of 4 to 6 V depending on the system. Further potential could not be applied during the CV curves due to the limitations of the utilized potentiostat (VSP Biologic). However, the curves obtained in this potential range are already representative of the electrochemical behavior of the system. Structural and compositional measurements: θ/2θ X-ray diffraction (XRD) patterns were collected on a Materials Research Diffractometer (MRD) from Malvern PANalytical company, equipped with a PIXcel1D detector, using Cu Kα radiation. High-resolution transmission electron microscopy (HRTEM), high-angle annular dark-field scanning transmission electron microscopy (HAADF-STEM), and electron energy loss spectroscopy (EELS) were carried out on a TECNAI F20 HRTEM/STEM microscope operated at 200 kV. Cross-sectional lamellae were prepared by focused ion beam, placed onto a copper transmission electron microscopy grid, and topped with a protective platinum layer. Doppler broadening variable energy positron annihilation spectroscopy (DB-VEPAS) and Variable energy positron annihilation lifetime spectroscopy (VEPALS): DB-VEPAS measurements were conducted at the setup for in–situ defect analysis (AIDA)55 of the slow positron beamline (SPONSOR).62 Positrons were accelerated and mono energetically implanted into samples in the range between Ep = 0.05 – 35 keV, which allows for depth sensitive analysis. The mean positron implantation depth was approximated using a simple material density (ρ) dependent formula: =36/ρ·Ep1.62. 63 Since at the annihilation site thermalized positrons have very small momentum compared to the electrons, a broadening of the 511 keV line is observed 26 mostly due to momentum of the electrons, which is measured with a high-purity Ge detector (overall energy resolution of 1.09 ± 0.01 keV at 511 keV). This broadening is characterized by parameters S defined as a fraction of the annihilation distribution in the middle (511±0.93 keV). The S-parameter is a fraction of positrons annihilating with low momentum valence electrons and represents vacancy type defects and their concentration.58 VEPALS measurements were conducted at the mono-energetic positron source (MePS) beamline, which is an end station of the radiation source ELBE (Electron Linac for beams with high Brilliance and low Emittance) at Helmholtz- Zentrum Dresden-Rossendorf (Germany). 56 A digital lifetime CeBr3 scintillator detector was used, with a homemade software employing a SPDevices ADQ14DC-2X with 14 bit vertical resolution and 2GS s–1 (GigaSamples per second) horizontal resolution and with a time resolution function down to about 0.230 ns.64 The resolution function required for spectrum analysis uses two Gaussian functions with distinct intensities depending on the positron implantation energy, Ep, and appropriate relative shifts. All spectra contained at least 1 × 107 counts. The spectra were deconvoluted using the non-linearly least-squared-based package PALSfit fitting software into discrete lifetime components, which directly confirm different defect types (i.e., sizes).60 Typical lifetime spectrum N(t) is described by N(t) =  (1/i) Ii exp(-t/i), where i and Ii are the positron lifetime and intensity of the i-th component, respectively (Ii = 1). The corresponding relative intensities (Ii) reflect to a large extent the concentration of each defect type and positron lifetimes (τi) are directly proportional to defect size (i.e., the larger the open volume, the lower the probability and the longer it takes for positrons to be annihilated with electrons). The positron lifetime and its intensity were probed as a function of positron implantation energy Ep or, in other words, implantation depth (thickness). The average positron lifetime τav is defined as τav = ∑i τi·Ii, which is the weighted average of the defect size. The shortest 27 lifetime component (τ1 < 0.32 ns) represents positron annihilation inside vacancy clusters (likely within grains) and/or at the grain boundaries, depending on the films microstructure. The intermediate lifetime (0.35 < τ2 < 0.90 ns) accounts for annihilation at larger vacancy clusters (linked to grain boundaries and their intersections), surface states, and small voids/pores (0.28 − 0.37 nm in diameter, calculated based on the shape-free model for pore-size estimation of Wada et al.59 The longest lifetime component (2.3 < τ3 < 3.3 ns) indicates contributions of larger voids (0.58 – 0.74 nm in diameter). ASSOCIATED CONTENT Supporting Information. Supporting Information is provided. This includes cyclic voltammetry curves of the basic investigated layer structures. AUTHOR INFORMATION Corresponding Authors Enric Menéndez, Departament de Física, Universitat Autònoma de Barcelona, E-08193 Cerdanyola del Vallès, Spain; Email: enric.menendez@uab.cat Jordi Sort, Departament de Física, Universitat Autònoma de Barcelona, E-08193 Cerdanyola del Vallès, Spain; Institució Catalana de Recerca i Estudis Avançats (ICREA), Pg. Lluís Companys 23, E-08010 Barcelona, Spain; Email: jordi.sort@uab.cat Author Contributions 28 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. ACKNOWLEDGMENTS Financial support by the European Union’s Horizon 2020 Research and Innovation Programme (‘BeMAGIC’ European Training Network, ETN/ITN Marie Skłodowska–Curie grant Nº 861145), the European Research Council (2021-ERC-Advanced ‘REMINDS’ Grant Nº 101054687), the Spanish Government (CEX2019-000917-S y PID2021-123276OB-I00, PID2020-116844RB-C21 and PDC2021-121276-C31) and the Generalitat de Catalunya (2021-SGR-00651) is acknowledged. J. S. thanks the Spanish “Fábrica Nacional de Moneda y Timbre” (FNMT) for fruitful discussions. E. M. is a Serra Húnter Fellow. Parts of this research were carried out at ELBE at the Helmholtz-Zentrum Dresden - Rossendorf e. V., a member of the Helmholtz Association. We would like to thank the facility staff for assistance. This work was partially supported by the Impulse-und Net-working fund of the Helmholtz Association (FKZ VH-VI-442 Memriox), and the Helmholtz Energy Materials Characterization Platform (03ET7015). REFERENCES (1) Shalf, J. M.; Leland, R. Computing Beyond Moore's Law. Computer 2015, 48, 14–23. (2) Ramesh, R.; Manipatruni, S. Electric Field Control of Magnetism. Proc. R. Soc. A. 2021, 477, 20200942. 29 (3) Žutić, I.; Fabian, J.; Das Sarma, S. Spintronics: Fundamentals and Applications. Rev. Mod. Phys. 2004, 76, 323–410. (4) Wolf, S. A.; Awschalom, D. D.; Buhrman, R. A.; Daughton, J. M.; von Molnar, S.; Roukes, M. L.; Chtchelkanova, A. Y.; Treger, D. M. Spintronics: A Spin-Based Electronics Vision for the Future. Science 2001, 294, 1488−1495. (5) Kawahara T.; Ito K.; Takemura R.; Ohno H. Spin-Transfer Torque RAM Technology: Review and Prospect. Microelectron. Reliab. 2012, 52, 613−627. (6) Brataas, A.; Kent, A. D.; Ohno, H. Current-Induced Torques in Magnetic Materials. Nat. Mater. 2012, 11, 372–381. (7) Han, X.; Wang, X.; Wan, C.; Yu, G.; Lv, X. Spin-Orbit Torques: Materials, Physics, and Devices. Appl. Phys. Lett. 2021, 118, 120502. (8) Gambardella, P.; Miron, I. M. Current-Induced Spin-orbit Torques. Phil. Trans. R. Soc. A. 2011, 369, 3175–3197. (9) Shiota, Y.; Nozaki, T.; Bonell, F.; Murakami, S.; Shinjo, T.; Suzuki, Y. Induction of Coherent Magnetization Switching in a Few Atomic Layers of FeCo Using Voltage Pulses. Nat. Mater. 2012, 11, 39–43. (10)Song, C.; Cui, B.; Li, F.; Zhou, X.; Pan, F. Recent Progress in Voltage Control of Magnetism: Materials, Mechanisms, and Performance. Prog. Mater. Sci. 2017, 87, 33–82. (11)Quintana, A.; Menéndez, E.; Liedke, M. O.; Butterling, M.; Wagner, A.; Sireus, V.; Torruella, P.; Estradé, S.; Peiró, F.; Dendooven, J.; Detavernier, C.; Murray, P. D.; Gilbert, D. A.; Liu, K.; Pellicer, E.; Nogués, J.; Sort, J. Voltage-Controlled ON–OFF Ferromagnetism at Room Temperature in a Single Metal Oxide Film. ACS Nano 2018, 12, 10291–10300. 30 (12)Bauer, U.; Yao, L.; Tan, A. J.; Agrawal, P.; Emori, S.; Tuller, H. L.; Van Dijken, S.; Beach, G. S. D. Magneto-ionic Control of Interfacial Magnetism. Nat. Mater. 2015, 14, 174–181. (13)Gilbert, D. A.; Grutter, A. J.; Arenholz, E.; Liu, K.; Kirby, B. J.; Borchers, J. A.; Maranville, B. B. Structural and Magnetic Depth Profiles of Magneto-ionic Heterostructures Beyond the Interface Limit. Nat. Commun. 2016, 7, 12264. (14)Rojas, J.; Quintana, A.; Lopeandia, A.; Salguero, J.; Costa-Krämer, J. L.; Abad, L.; Liedke, M. O.; Butterling, M.; Wagner, A.; Henderick, L.; Dendooven, J.; Detavernier, C.; Sort, J.; Menéndez, E. Boosting Room-Temperature Magneto-ionics in a Non-Magnetic Oxide Semiconductor. Adv. Funct. Mater. 2020, 30, 2003704. (15)Martins, S.; De Rojas, J.; Tan, Z.; Cialone, M.; Lopeandia, A.; Herrero-Martín, J.; Costa- Krämer, J. L.; Menéndez, E.; Sort, J. Dynamic Electric-field-induced Magnetic Effects in Cobalt Oxide Thin Films: Towards Magneto-ionic Synapses. Nanoscale 2022, 14, 842–852. (16)Dasgupta, S.; Das, B.; Knapp, M.; Brand, R. A.; Ehrenberg, H.; Kruk, R.; Hahn, H. Intercalation-driven Reversible Control of Magnetism in Bulk Ferromagnets. Adv. Mater. 2014, 26, 4639–4644. (17)Vasala, S.; Jakob, A.; Wissel, K.; Waidha, A. I.; Alff, L.; Clemens, O. Reversible Tuning of Magnetization in a Ferromagnetic Ruddlesde-Popper-Type Manganite by Electrochemical Fluoride Ion Intercalation. Adv. Electron. Mater. 2020, 6, 1900974. (18)Tan, A. J.; Huang, M.; Avci, C. O.; Büttner, F.; Mann, M.; Hu, W.; Mazzoli, C.; Wilkins, S.; Tuller, H. L.; Beach, G. S. D. Magneto-Ionic Control of Magnetism Using a Solid-State Proton Pump. Nat. Mater. 2019, 18, 35–41. (19)Ye, X.; Singh, H. K.; Zhang, H.; Geßwein, H.; Chellali, M. R.; Witte, R.; Molinari, A.; Skokov, K.; Gutfleisch, O.; Hahn, H.; Kruk, R. Giant Voltage-Induced Modification of 31 Magnetism in Micron-Scale Ferromagnetic Metals by Hydrogen Charging. Nat. Commun. 2020, 11, 4849. (20)Chen, G.; Ophus, C.; Quintana, A.; Kwon, H.; Won, C.; Ding, H.; Wu, Y.; Schmid, A. K.; Liu, K. Reversible Writing/deleting of Magnetic Skyrmions Through Hydrogen Adsorption/desorption. Nat. Commun. 2022, 13, 1350. (21)Lu, N.; Zhang, Z.; Wang, Y.; Li, H.-B.; Qiao, S.; Zhao, B.; He, Q.; Lu, S.; Li, C.; Wu, Y.; Zhu, M.; Lyu, X.; Chen, X.; Li, Z.; Wang, M.; Zhang, J.; Tsang, S. C.; Guo, J.; Yang, S. Zhang, J.; Deng, K.; Zhang, D.; Ma, J.; Ren, J.; Wu, Y.; Zhu, J. Zhou, S.; Tokura,Y.; Nan, C.-W.; Wu, J.; Yu, P. Enhanced Low-Temperature Proton Conductivity in Hydrogen- Intercalated Brownmillerite Oxide. Nat. Energy 2022, 7, 1208−1216. (22)De Rojas, J.; Quintana, A.; Lopeandía, A.; Salguero, J.; Muñiz, B.; Ibrahim, F.; Chshiev, M.; Nicolenco, A.; Liedke, M. O.; Butterling, M.; Wagner, A.; Sireus, V.; Abad, L.; Jensen, C. J.; Liu, K.; Nogués, J.; Costa-Krämer, J. L.; Menéndez, E.; Sort, J. Voltage-Driven Motion of Nitrogen Ions: A New Paradigm for Magneto-Ionics. Nat. Commun. 2020, 11, 5871. (23)De Rojas, J.; Salguero, J.; Ibrahim, F.; Chshiev, M.; Quintana, A.; Lopeandia, A.; Liedke, M. O.; Butterling, M.; Hirschmann, E.; Wagner, A.; Abad, L.; Costa-Krämer, J. L.; Menéndez, E.; Sort, J. Magneto-ionics in Single-layer Transition Metal Nitrides. ACS Appl. Mater. Interfaces 2021, 13, 30826–30834. (24)De Rojas, J.; Salguero, J.; Quintana, A.; Lopeandia, A.; Liedke, M. O.; Butterling, M.; Attallah, A. G.; Hirschman, E.; Wagner, A.; Abad, L.; Costa-Krämer, J. L.; Sort, J.; Menéndez, E. Critical Role of Electrical Resistivity in Magnetoionics. Phys. Rev. Appl. 2021, 16, 034042. 32 (25)Tan, Z.; Martins, S.; Escobar, M.; De Rojas, J.; Ibrahim, F.; Chshiev, M.; Quintana, A.; Lopeandia, A.; Costa-Krämer, J. L.; Menéndez, E. From Binary to Ternary Transition-metal Nitrides: A Boost Toward Nitrogen Magneto-ionics. ACS Appl. Mater. Interfaces 2022, 14, 44581–44590. (26)Tan, Z.; De Rojas, J.; Martins, S.; Lopeandia, A.; Quintana, A.; Cialone, M.; Herrero-Martín, J.; Meersschaut, J.; Vantomme, A.; Costa-Krämer, J. L.; Sort, J.; Menéndez, E. Frequency- dependent Stimulated and Post-stimulated Voltage Control of Magnetism in Transition Metal Nitrides: Towards Brain-Inspired Magneto-ionics. Mater. Horiz. 2023, 10, 88–96. (27)Navarro-Senent, C.; Quintana, A.; Menéndez, E.; Pellicer, E.; Sort, J. Electrolyte-gated Magnetoelectric Actuation: Phenomenology, Materials, Mechanisms, and Prospective Applications. APL Mater. 2019, 7, 030701. (28)De Rojas, J.; Quintana, A.; Rius, G.; Stefani, C.; Domingo, N.; Costa-Krämer, J. L.; Menéndez, E.; Sort, J. Voltage Control of Magnetism with Magneto-Ionic Approaches: Beyond Voltage-Driven Oxygen Ion Migration. Appl. Phys. Lett. 2022, 120, 070501. (29)Zehner, J.; Wolf, D.; Hasan, M. U.; Huang, M.; Bono, D.; Nielsch, K.; Leistner, K.; Beach, G. S. D. Magnetoionic Control of Perpendicular Exchange Bias. Phys. Rev. Mater. 2021, 5, L061401. (30)Martins, S.; Ma, Z.; Solans-Monfort, X.; Sodupe, M.; Rodriguez-Santiago, L.; Menéndez, E.; Pellicer, E.; Sort, J. Enhancing Magneto-ionic Effects in Cobalt Oxide Films by Electrolyte Engineering. Nanoscale Horiz. 2023, 8, 118–126. (31)Zhang, Q.; Luo, X.; Wang, L.; Zhang, L.; Khalid, B.; Gong, J.; Wu, H. Lithium-Ion Battery Cycling for Magnetism Control. Nano Lett. 2016, 16, 583–587. 33 (32)Quintana, A.; Firme, A. A.; Jensen, C. J.; Zheng, D.; Liu, C.; Zhang, X.; Liu, K. Hydroxide- Based Magneto-ionics: Electric-Field Control of Reversible Paramagnetic-to- Ferromagnetic Switch in α-Co(OH)2 films. J. Mater. Chem. C 2022, 10, 17145. (33)Walter, J.; Voigt, B.; Day-Roberts, E., Heltemes, K.; Fernandes, R. M.; Birol, T.; Leighton, C. Voltage-induced Ferromagnetism in a Diamagnet. Sci. Adv. 2020, 6, eabb7721. (34)Molinari, A.; Hahn, H.; Kruk, R. Voltage-control of Magnetism in All-solid-state and Solid/liquid Magnetoelectric Composites. Adv. Mater. 2019, 31, 1806662. (35)Pachat, R.; Ourdani, D.; Syskaki, M.; Lamperti, A.; Roy, S.; Chen, S.; Pietro, A. D.; Largeau, L.; Juge, R.; Massouras, M.; Balan, C.; Van Der Jagt, J. W.; Agnus, G.; Roussigné, Y.; Gabor, M.; Chérif, S. M.; Durin, G.; Ono, S.; Langer, J.; Querlioz, D.; Ravelosona, D.; Belmeguenai, M.; Herrera Diez, L. Magneto-Ionics in Annealed W/CoFeB/HfO2 Thin Films. Adv. Mater. Interfaces 2022, 9, 2200690. (36)Zhou, X.; Yan, Y.; Jiang, M.; Cui, B.; Pan, F.; Song, C. Role of Oxygen Ion Migration in the Electrical Control of Magnetism in Pt/Co/Ni/HfO2 Films. J. Phys. Chem. C 2016, 120, 1633–1639. (37)Herrera Diez, L.; Kruk, R.; Leistner, K.; Sort, J. Magnetoelectric Materials, Phenomena, and Devices. APL Mater. 2021, 9, 050401. (38)Chen, G.; Mascaraque, A.; Jia, H.; Zimmermann, B.; Robertson, M.; Lo Conte, R.; Hoffmann, M.; González Barrio, M. A.; Ding, H.; Wiesendanger, R.; Michel, E. G.; Blügel, S.; Schmid, A. K.; Liu, K. Large Dzyaloshinskii-Moriya Interaction Induced by Chemisorbed Oxygen on a Ferromagnet Surface. Sci. Adv. 2020, 6, eaba4924. (39)Chen, G.; Robertson, M.; Hoffmann, M.; Ophus, C.; Fernandes Cauduro, A. L.; Lo Conte, R.; Ding, H.; Wiesendanger, R.; Blügel, S.; Schmid, A. K.; Liu, K. Observation of 34 Hydrogen-Induced Dzyaloshinskii-Moriya Interaction and Reversible Switching of Magnetic Chirality. Phys. Rev. X 2021, 11, 021015. (40)Huang, M.; Usama Hasan, M.; Klyukin, K.; Zhang, D.; Lyu, D.; Gargiani, P.; Valvidares, M.; Sheffels, S.; Churikova, A.; Büttner, F.; Zehner, J.; Caretta, L.; Lee, K.-Y.; Chang, J.; Wang, J.-P.; Leistner, K.; Yildiz, B.; Beach, G. S. D. Voltage Control of Ferrimagnetic Order and Voltage-assisted Writing of Ferrimagnetic Spin Textures. Nat. Nanotechnol. 2021, 16, 981−988. (41)Zehner, J.; Huhnstock, R.; Oswald, S.; Wolff, U.; Soldatov, I.; Ehresmann, A.; Nielsch, K.; Holzinger, D.; Leistner, K. Nonvolatile Electric Control of Exchange Bias by a Redox Transformation of the Ferromagnetic Layer. Adv. Electron. Mater. 2019, 5, 1900296. (42)Lecuyer, S.; Quemerais, A.; Jezequel, G. Composition of Natural Oxide Films on Polycrystalline Tantalum Using XPS Electron Take-Off Angle Experiments. Surf. Interface Anal. 1992, 18, 257–261. (43)Levine, B. R.; Sporer, S.; Poggie, R. A.; Della Valle, C. J.; Jacobs, J. J. Experimental and Clinical Performance of Porous Tantalum in Orthopedic Surgery. Biomaterials, 2006, 27, 4671–4681. (44)Kevin S.; Damian B.; Steffen F.; Oliver P.; Detlef D. Preparation and Properties of Thin Amorphous Tantalum Films Formed by Small e-beam Evaporators. J. Phys. D: Appl. Phys. 2009, 42, 135417. (45)Dalla Torre, J.; Gilmer, G. H.; Windt, D. L.; Kalyanaraman, R.; Baumann, F. H.; O’Sullivan, P. L.; Sapjeta, J.; Dı́Az De La Rubia, T.; Djafari Rouhani, M. Microstructure of Thin Tantalum Films Sputtered onto Inclined Substrates: Experiments and Atomistic Simulations. J. Appl. Phys. 2003, 94, 263–271. 35 (46)Liu, Q.; Wu, F.; Mu, D.; Wu, B. A Theoretical Study on Na+ Solvation in Carbonate Ester and Ether Solvents for Sodium-ion Batteries. Phys. Chem. Chem. Phys. 2020, 22, 2164– 2175. (47)Weisheit, M.; Fähler, S.; Marty, A.; Souche, Y.; Poinsignon, C.; Givord, D. Electric Field- Induced Modification of Magnetism in Thin-Film Ferromagnets. Science 2007, 315, 349– 351. (48)Silva, C.; Vovk, A.; Da Silva, R.C.; Strichovanec, P.; Algarabel, P.A.; Gonçalves, A.P.; Borges, R.P.; Godinho, M.; Cruz, M. M. Magnetic Properties of Co–N Thin Films Deposited by Reactive Sputtering. Thin Solid Films 2014, 556, 125–127. (49)Li, Z. L.; Du, B. X.; Yang, Z. R.; Li, J. Effects of Crystal Morphology on Space Charge Transportation and Dissipation of Sic/Silicone Rubber Composites. IEEE Trans. Dielectr. Electr. Insul. 2017, 24, 2616–2625. (50)Krikorian, E.; Sneed, R. J. Deposition of Tantalum, Tantalum Oxide, and Tantalum Nitride with Controlled Electrical Characteristics. J. Appl. Phys. 1996, 37, 3674–3681. (51)Gerstenberg, D.; Calbick, C. J. Effects of Nitrogen, Methane, and Oxygen on Structure and Electrical Properties of Thin Tantalum Films. J. Appl. Phys. 1964, 35, 402–407. (52)De Vries, J. W. C. Temperature and Thickness Dependence of the Resistivity of Thin Polycrystalline Aluminium, Cobalt, Nickel, Palladium, Silver and Gold Films. Thin Solid Films, 1988, 167, 25–32. (53)Zubler, E. G. The Gettering Properties of Tantalum. J. Electrochem. Soc. Jpn. 1963, 110, 1072. (54)Bean, C. P.; Fisher, J. C.; Vermilyea, D. A. Ionic Conductivity of Tantalum Oxide at Very High Fields. Phys. Rev. 1956, 101, 551. 36 (55)Liedke, M. O.; Anwand, W.; Bali, R.; Cornelius, S.; Butterling, M.; Trinh, T. T.; Wagner, A.; Salamon, S.; Walecki, D.; Smekhova, A. Open Volume Defects and Magnetic Phase Transition in Fe60Al40 Transition Metal Aluminide. J. Appl. Phys. 2015, 117, 163908. (56)Wagner, A.; Butterling, M.; Liedke, M. O.; Potzger, K.; Krause-Rehberg, R. Positron Annihilation Lifetime and Doppler Broadening Spectroscopy at the ELBE Facility. AIP Conference Proceedings, 2018, 1970, 040003. (57)Tuomisto, F.; Makkonen, I. Defect Identification in Semiconductors with Positron Annihilation: Experiment and Theory, Rev. Mod. Phys. 2013, 85, 1583–1631. (58)Clement, M.; De Nijs, J. M. M.; Balk, P., Schut, H.; Van Veen, A. Analysis of Positron Beam Data by the Combined Use of the Shape-and Wing-Parameters. J. Appl. Phys. 1996, 79, 9029–9036. (59)Xu, Q.; Popov, E.; Troev, T.; Zhang, J.; Dai, Y. Positron Lifetime Calculation of Vacancy Clusters in Tantalum Containing Hydrogen and Helium. J. Nucl. Mater. 2018, 506, 71–75. (60)Olsen, J. V.; Kirkegaard, P.; Pedersen, N.J.; Eldrup, M. PALSfit: A New Program for the Evaluation of Positron Lifetime Spectra, Phys. Status Solidi. 2007, 4, 4004–4006. (61)Wada K.; Hyodo T. A Simple Shape-Free Model for Pore-Size Estimation with Positron Annihilation Lifetime Spectroscopy, J. Phys. Conf. Ser. 2013, 443, 012003. (62)Anwand, W.; Brauer, G.; Butterling, M.; Kissener, H. R.; Wagner, A. Design and Construction of a Slow Positron Beam for Solid and Surface Investigations, Defect Diffus. Forum. 2012, 331, 25–40. (63)Dryzek, J.; Horodek, P. GEANT4 Simulation of Slow Positron Beam Implantation Profiles. Nucl. Instruments Methods Phys. Res. Sect. B Beam Interact. with Mater. Atoms. 2008, 266, 4000–4009. 37 (64)Hirschmann, E.; Butterling, M.; Hernandez Acosta, U.; Liedke, M.O.; Attallah, A.G.; Petring, P.; Görler, M.; Krause-Rehberg, R.; Wagner, A., A New System for Real-time Data Acquisition and Pulse Parameterization for Digital Positron Annihilation Lifetime Spectrometers With High Repetition rates, J. Instrum. 2021, 16, 08001. TABLE OF CONTENTS 38