Research ArticleSPINTRONICS

Anomalous spin-orbit torque switching due to field-like torque–assisted domain wall reflection

+ See all authors and affiliations

Science Advances  21 Apr 2017:
Vol. 3, no. 4, e1603099
DOI: 10.1126/sciadv.1603099

Abstract

Spin-orbit torques (SOTs) allow the electrical control of magnetic states. Current-induced SOT switching of the perpendicular magnetization is of particular technological importance. The SOT consists of damping-like and field-like torques, and understanding the combined effects of these two torque components is required for efficient SOT switching. Previous quasi-static measurements have reported an increased switching probability with the width of current pulses, as predicted considering the damping-like torque alone. We report a decreased switching probability at longer pulse widths, based on time-resolved measurements. Micromagnetic analysis reveals that this anomalous SOT switching results from domain wall reflections at sample edges. The domain wall reflection was found to strongly depend on the field-like torque and its relative sign to the damping-like torque. Our result demonstrates a key role of the field-like torque in deterministic SOT switching and the importance of the sign correlation of the two torque components, which may shed light on the SOT switching mechanism.

Keywords
  • Spintronics
  • spin-orbit torque
  • magnetization switching
  • magneto optical Kerr effect
  • domain wall

INTRODUCTION

Spin-orbit coupling can convert charge currents to spin currents (18). The ability to generate spin currents without the help of ferromagnets constitutes a core building block of an emerging research field, spin orbitronics, which uses spin-orbit coupling as a spin-current source in spintronic devices. When the spin current from spin-orbit coupling is absorbed by a ferromagnet, it exerts a spin-orbit torque (SOT) on the ferromagnet. The SOT can switch magnetization (9, 10) and induce fast domain wall motion (1113) in ferromagnet/heavy metal bilayers. The SOT switching of the perpendicular magnetization is of particular technological relevance because perpendicular magnetic random access memories have a better scalability than in-plane ones.

From the viewpoint of fundamental physics and its applications, it is of critical importance to understand the detailed SOT characteristics and consequent magnetization dynamics. The microscopic origin of SOT remains under debate (9, 10, 1426), but it is commonly decomposed into two mutually orthogonal vector components, the damping-like torque (DLT) and field-like torque (FLT). SOT-induced magnetization dynamics, including the DLT and FLT terms, is described by the Landau-Lifshitz-Gilbert equationEmbedded Image(1)where γ is the gyromagnetic ratio, Embedded Image is the unit vector along the magnetization, μ0Heff is the effective magnetic field, including the external, anisotropic, magnetostatic, and exchange fields, and α is the damping constant. τd (Embedded Image) and τf (Embedded Image) describe the magnitudes of DLT and FLT, respectively (in magnetic field units), where J is the current density, Ms is the saturation magnetization, d is the thickness of ferromagnet, c|| and c are the DLT and FLT efficiencies, respectively, and ŷ is the unit vector perpendicular to both the current direction and the inversion asymmetry direction (that is, thickness direction; see coordinate system in Fig. 1A). From Eq. 1, we found that the two torque components affect magnetization dynamics in distinctly different ways: The DLT directs the magnetization toward the y axis, whereas the FLT induces magnetization precessions around the y axis.

Fig. 1 TR-MOKE experimental setup and dc-induced magnetization switching in Ta/CoFeB/MgO.

(A) Schematic illustration of TR-MOKE measurements. μ0Hx is the external magnetic field. The dc or pulse current is applied along the x axis. The picosecond laser is shined as a probe beam. The patterned perpendicular anisotropy Ta (6 nm)/CoFeB (0.8 nm)/MgO (2 nm) square is connected to a ground (G)–signal (S)–ground (G) coplanar waveguide. (B) dc-induced magnetization switching with various μ0Hx. The data are shifted vertically for clarity. (C) dc switching current density Jc versus μ0Hx.

Most of the previous SOT-switching studies have considered DLT as a main driving source but ignored the FLT. When considering DLT alone, the switching trajectory in the macrospin approximation is expected to be simple without magnetization precessions (27). This DLT-dominated switching leads to an increased switching probability with the current pulse width, consistent with previous quasi-static measurements (28, 29) and also in accordance with common belief. However, in some ferromagnet/heavy metal bilayers (for example, Ta-based bilayers), the FLT is significant (3032). For a sizable FLT, it has been theoretically predicted that the magnetization precession induced by the FLT complicates magnetization dynamics, especially for the case wherein the sign of FLT is the opposite to that of DLT (in our sign convention; see Eq. 1) (3335). In this respect, it is important to experimentally investigate the role of FLT in the SOT switching. Note that previous quasi-static measurements based on Hall bar detection (9, 10, 28, 29) or magneto-optical Kerr effect (MOKE) microscopy (36) would be unable to capture the core effect of FLT because this may induce fast dynamics (that is, magnetization precessions around the y axis). To overcome this limitation, it is essential to perform time-resolved (TR) measurements (37, 38), which provide an important step toward a better understanding of SOT-induced magnetization dynamics.

Here, we report SOT-induced magnetization dynamics in time domain by TR-MOKE measurements for Ta/CoFeB/MgO heterostructures with perpendicular magnetic anisotropy. This Ta-based structure has a large FLT (that is, |τfd| ≈ 4) whose sign is the opposite to that of DLT (32), such that it allows a detailed study of the FLT effect on SOT switching. The temporal evolution of the magnetization is detected by the stroboscopic pump-probe technique with an electrical pulse generator (pump) and a picosecond laser (probe) (see section S1), as shown schematically in Fig. 1A. We observe an anomalous SOT-induced switching behavior, in which the switching probability increases at short current pulses but decreases at longer pulses. On the basis of micromagnetic simulations, we interpret this anomalous switching behavior as a consequence of FLT-assisted domain wall reflection at sample edges.

RESULTS

TR measurements of perpendicular magnetization switching by SOT

We first perform a static polar MOKE measurement using an in-plane dc I with an external in-plane magnetic field μ0Hx along the x axis to examine the dc SOT switching characteristics of the Ta (6 nm)/CoFeB (0.8 nm)/MgO (2 nm) sample (device 1). The pattern size of device 1 is 3 μm × 3 μm, which is large enough to detect SOT-induced change in the MOKE signal. Figure 1B shows MOKE signals as a function of dc with various in-plane magnetic fields. As a polar MOKE signal probes the average z component of CoFeB magnetization (〈Mz〉) in a laser spot, where a high (low) signal corresponds to the 〈Mz〉-up (〈Mz〉-down) magnetic state, the hysteretic curves show SOT-induced deterministic magnetization switching. The switching polarity is determined by the direction of the current and Hx (9, 10): The up-to-down switching occurs in a positive (negative) current with a negative (positive) Hx. As shown in Fig. 1C, the dc switching current density Jc decreases with increasing Hx, in agreement with previous reports using Hall bar experiments (9, 20).

We next carry out TR measurements by injecting a current pulse with various pulse widths (that is, tpw ≤ 5 ns, J = 5.2 × 1011 A m−2, and μ0Hx = −168 mT) and by detecting the magnetic state in time domain through TR-MOKE signal (Fig. 2A). The current pulse is turned on at t = 0 (that is, current-on) and turned off at a time indicated by a red triangle in each curve (that is, current-off). The horizontal dashed lines serve as guides for the maximum change in MOKE signal (~7 μV), corresponding to full magnetization switching from the up to the down state.

Fig. 2 TR-MOKE measurements of SOT-induced perpendicular magnetization switching in Ta/CoFeB/MgO.

(A) Temporal evolutions of TR-MOKE signals corresponding to the average z component of magnetization (Mz) in an applied current density (J = 5.2 × 1011 A m−2) of various pulse widths (tpw) from 1 to 5 ns for μ0Hx = −168 mT. The data are shifted vertically for clarity. The current pulse starts at t = 0, and the end of the current pulse is indicated as a red triangle in each curve. The horizontal dashed lines serve as guides to the maximum change in MOKE signal (~7 μV), corresponding to the full magnetization switching from the up to the down state. (B) Time-varying MOKE signal for tpw = 1.8 and 5.0 ns. (C) Switching probability (Psw) as a function of tpw, extracted from Fig. 2A.

For short pulses (tpw ≤ 1.6 ns), the MOKE signal change is smaller than 7 μV, indicating that the current pulse is too short to switch the magnetization. For an increase of tpw to 1.8 ns, a complete switching is achieved, as evidenced by a signal change of ~7 μV (by defining the final magnetic state at t = 8 ns). However, when tpw > 2.5 ns, we observe an anomalous temporal change in the TR-MOKE signal: It decreases in the initial time stage (<2 ns) but increases back even before the current pulse is turned off.

To clarify the anomalous SOT switching behavior, we compare two cases, tpw = 1.8 and 5.0 ns, in Fig. 2B. For tpw = 1.8 ns, which corresponds to normal switching, the TR-MOKE signal decreases monotonously when the current is on, indicating that 〈Mz〉 changes monotonously from the up to the down state by SOT. This decreased signal is maintained even after the current is turned off so that the up-to-down switching is completed. On the other hand, for tpw = 5.0 ns, which corresponds to anomalous switching, the signal decreases until t = 2 ns is reached and then increases back. Note that this anomalous increase in the signal is present even before the current is off. It implies that the SOT is responsible not only for the initial decrease in the MOKE signal (that is, switching) but also for the anomalous increase in the signal (that is, switching back). The anomalous switching back phenomena are observed for all cases with tpw > 2.0 ns, as summarized in Fig. 2C, which shows the switching probability (PSW) as a function of tpw, where PSW (%) = [1 – VMOKE (t = 8 ns)/7 μV] × 100.

We next show that the anomalous switching back phenomena become more pronounced as the current amplitude increases. We perform TR-MOKE measurements at various current densities and pulse widths for a structure of Ta (3 nm)/CoFeB (1.2 nm)/MgO (2 nm) with a pattern size of 3 μm × 6 μm (device 2). As an example, Fig. 3A shows TR-MOKE signals normalized by the maximum signal change as a function of time at various current densities for tpw = 30 ns. Focusing on the normalized MOKE signal at t = 45 ns (that is, after the current is turned off), we find that a higher current density causes more switching back. Figure 3B summarizes PSW versus tpw at various current densities. It shows that the switching back phenomena become more noticeable at higher current densities.

Fig. 3 Current-dependence of SOT-induced magnetization switching.

The sample structure is Ta (3 nm)/CoFeB (1.2 nm)/MgO (2 nm)/SiO2 (3 nm). (A) Temporal evolutions of TR-MOKE signals, normalized by the maximum signal change, for various current densities for tpw = 30 ns and μ0Hx = −90 mT. (B) Switching probability (PSW) as a function of tpw at various current densities.

Anomalous SOT switching due to domain wall reflections at edges: Micromagnetic analysis

To understand the anomalous switching back phenomena, we perform micromagnetic simulations at 0 K (Fig. 4; see Materials and Methods). We use parameters of magnetic properties and SOTs [DLT (c|| = −0.07) and FLT (c = +0.28)], deduced from magnetometer (see section S2) and harmonic Hall measurements for Ta/CoFeB/MgO samples (32). Figure 4A shows temporal evolutions of normalized 〈mz〉 at various current pulse widths for J = 15 × 1011 A m−2 and μ0Hx = −200 mT. The current-off time is depicted as a vertical dashed line for each case. For tpw = 1.5 and 1.6 ns, the current pulse is too short to switch the magnetization, whereas for tpw = 1.7 ns, a full switching is achieved. The switching back is observed for a longer pulse (tpw = 1.8 ns): 〈mz〉 returns back to the initial state (〈mz〉 = +1) after the current is turned off. Figure 4B shows a switching parameter P as a function of tpw, where P is defined as “1” (“0”) for the switching (no-switching) event. We find that the switching back (equivalently no-switching event) is not unique for the case of tpw = 1.8 ns but appears in a somewhat oscillatory manner for tpw > 1.8 ns. Note that the oscillation is not periodic. This nonperiodic oscillatory switching obtained from the zero-temperature calculation (Fig. 4B) can explain the decreased switching probability at longer pulses observed in the room temperature measurement (Fig. 2C) because the thermal effect randomizes the oscillatory switching dynamics. In section S3, we show micromagnetic simulation results for detailed oscillatory switching dynamics.

Fig. 4 Micromagnetic simulation results for SOT switching at 0 K.

(A) Temporal evolutions of average Mz/Ms (〈mz〉) at various current pulse widths. (B) Switching parameter (P; 0 = no-switching and 1 = switching) as a function of tpw for J = 15 × 1011 A m−2 and μ0Hx = −200 mT. (C) Temporal evolutions of 〈mz〉 at various current densities for tpw = 1.2 ns and μ0Hx = −200 mT. Snapshots of magnetization configuration (mz = Mz/Ms) at time t for tpw = 1.7 ns (D) and 1.8 ns (E). Yellow arrows show the direction of domain wall motion. (F) Switching parameter P as a function of tpw for various current densities (J = 16 × 1011, 18 × 1011, and 20 × 1011 A m−2).

Micromagnetic simulations reveal that the switching back phenomena originate from the domain wall reflection at sample edges in the presence of current (thus SOT), as discussed below. We show temporal evolutions of mz for tpw = 1.7 and 1.8 ns, corresponding to the switching (Fig. 4D) and switching back (Fig. 4E), respectively. For both cases, in the initial time stage (t < 1.7 ns), a reversed domain is nucleated at a corner and then expands isotropically. For tpw = 1.7 ns, the current is turned off at the moment when the domain wall arrives at sample edges. In this case, the domain wall remains moving in the same direction as before because of inertia (3941), even after the current is turned off. As a result, a full switching is achieved. On the other hand, for tpw = 1.8 ns, the current is still turned on when the domain wall arrives at sample edges. In this case, the domain wall is reflected from the edges and moves in the opposite direction, leading to switching back. Therefore, the domain wall reflection at sample edges in the presence of SOT is key to explaining the anomalous switching back phenomena.

The domain wall reflection at sample edges also explains the current-dependent switching back behavior shown in Fig. 3. Figure 4C shows temporal evolutions of 〈mz〉 at various current densities when tpw = 1.2 ns and μ0Hx = −200 mT. At a low current density, a reversed domain is not nucleated (J = 14 × 1011 A m−2), or the domain wall is unable to reach the sample edges (J = 16 × 1011 A m−2). At an increased current density (J = 18 × 1011 A m−2), a full switching occurs. However, at a higher current density (J = 20 × 1011 A m−2), 〈mz〉 initially decreases but returns to the initial state because of the domain wall reflection (see section S3). Figure 4F shows the switching parameter P as a function of tpw at various current densities. Similar to the results shown in Fig. 4B, the switching back appears in a somewhat oscillatory manner. The no-switching event (P = 0) becomes more frequent for a higher current density when tpw exceeds a threshold to enable the switching. This behavior is consistent with experimental observations (Fig. 3), in which the switching back phenomena become more pronounced as the current amplitude increases.

Effect of FLT on the domain wall reflection: Collective coordinate analysis

To understand the effect of FLT on the domain wall reflection, we investigate domain wall dynamics based on a semi–one-dimensional micromagnetic model. We first show the domain wall moving along the bottom edge in a two-dimensional sample (Fig. 4, D and E), where a Néel-type domain wall is stabilized by Hx at this edge. Figure 5A shows temporal evolutions of the domain wall position at various FLT to DLT ratios (that is, τfd = c/c||) for J = 6 × 1011 A m−2 and μ0Hx = −200 mT. Here, we fix the DLT efficiency c|| as −0.07 and vary the FLT efficiency c. For all cases, the domain wall is reflected at the edge (located at 2 μm), exhibits a backward motion for a while, and then moves back to the edge again. After several reflections, the domain wall eventually annihilates at the edge, and the switching is completed. This domain wall reflection is understood by the reflection of a transverse wave at a fixed end. It is known that the phase of a transverse wave changes by π when the wave is reflected at a fixed end (Fig. 5B). Because a domain wall can be decomposed into transverse spin waves, the phase change corresponds to a change in the domain wall angle φ upon reflection (see Fig. 5C for schematics). The most important feature in Fig. 5A is that the distance Δq (defined as a positive value) for the backward domain wall motion strongly depends on the magnitude and sign of τfd. For a positive τfd, Δq is small and the domain wall annihilates soon after it reaches the sample edge. On the other hand, for a negative τfd, Δq is large, which in turn causes a noticeable switching back behavior, as shown in Fig. 4E. For a domain wall moving along the left edge (parallel to the y axis), where the Bloch-type domain wall is stabilized by Hx, it experiences a similar reflection process due to the same symmetry of domain wall angle and SOT (see section S4).

Fig. 5 Domain wall reflection in one-dimensional model.

(A) Temporal evolutions of the domain wall position q at various FLT to DLT ratios (τfd = c/c||) for c|| = −0.07, J = 6 × 1011 A m−2, and μ0Hx = −200 mT. Δq is the distance for the backward motion of a reflected domain wall. Schematic illustrations of the reflection of a transverse wall at a fixed end (B), and the reflection of a domain wall at an edge (C). (D) Temporal evolutions of domain wall angle φ and domain wall position q for τfd = −2. Domain wall angles φstd, φM, φref, and φ0 are defined in (C) and (D) (see the text for details). (E) φref, φ0, and φstd as a function of τfd. (F) Δq as a function of τfd.

We adopt the collective coordinate approach for the domain wall position q and domain wall angle φ (see section S5) (4246) to explain the dependence of Δq on the magnitude and sign of τfd. We define three domain wall angles, φref, φ0, and φstd (see Fig. 5, C and D): φref is the angle just after the reflection, φ0 is the angle at which the backward domain wall motion reverses to the forward motion, and φstd is the angle for the steady-state motion before the reflection. Note that φref = 2φM − φstd, where φM(≡π − tan−1f/Hx)) describes the domain tilting in the film plane (see Fig. 5C). Because the edge acts as a fixed end, the domain wall component transverse to the azimuthal angle φM of the magnetization is reversed upon reflection. That is, φstd, the angle of incoming domain wall, can be rewritten as φstd = φM + (φstd − φM), where the first (second) term is longitudinal (transverse) to the domain angle φM. Upon reflection, only a transverse component is reversed (that is, changes its sign), whereas a longitudinal component is conserved. Therefore, φref = φM − (φstd − φM) = 2φM − φstd.

From the collective coordinate approach, we obtain for −Hx > τf > 0 (see section S5)Embedded Image(2)Embedded Image(3)Embedded Image(4)where Embedded Image, Embedded Image, Embedded Image, and HK,eff is the effective anisotropy field. We find that Eqs. 2 to 4 describe general tendencies of the numerically obtained three domain wall angles with respect to τfd (Fig. 5E). Some disagreement for φref can be attributed to the dynamically distorted domain wall profile just after the reflection.

From the collective coordinate approach, assuming a small damping, an approximate Δq is given as (see section S5 for details)Embedded Image(5)where λ is the domain wall width. Figure 5F shows that Eq. 5 qualitatively describes the numerical results of the dependence of Δq on τfd. Note that the first term of Eq. 5 dominates over the second term (Fig. 5F). Therefore, the dependence of Δq on τfd is mostly governed by cos φ0/cos φref. Also note that φref changes more rapidly with τfd than φ0 (Fig. 5E). Therefore, the FLT dependence of φref is key to understanding a large backward domain wall motion for a negative τfd. From Eq. 3, the FLT-dependence of φref is described by −2 tan−1f/Hx) (thus, φM) in a small damping approximation, which means that the FLT affects the backward domain wall motion through its effect on φM.

DISCUSSION

We demonstrate that the FLT has a crucial role in the SOT switching and causes the anomalous switching back phenomena when it is large and its sign is the opposite to that of DLT. Our result suggests that not only DLT but also FLT should be carefully examined to achieve deterministic SOT switching. Furthermore, our result raises a question on the microscopic origin of SOT: What determines the sign correlation between DLT and FLT? To our knowledge, the sign product of DLT and FLT in all previous experiments has been negative (in our sign convention) except for a data point by Kim et al. (30), of which the spin-orbit effective field is too small to unambiguously determine the sign product. For instance, in our previous work (26), the signs of DLT and FLT vary with the oxidation of a ferromagnet, but the sign product is always negative. If this fixed sign product is true, it indicates that a single dominant mechanism is responsible for both DLT and FLT. It is in contrast to the currently widely accepted argument that the DLT (FLT) originates from the bulk spin Hall effect in a heavy metal layer (the Rashba effect at the ferromagnet/heavy metal interface). Notably, in theories for the bulk spin Hall mechanism, the sign product is positive (20) for a positive imaginary part of the spin-mixing conductance (47). On the other hand, in theories for the Rashba mechanism, the sign product is negative (17, 25).

We have demonstrated the anomalous switching back phenomena in rather large samples (a few square micrometers), which is required to obtain detectable MOKE signals. However, note that these results could be applicable to nanosized samples (a few tens of nanometers), as long as the sample size is larger than the domain wall width because the domain wall dynamics is key to the anomalous switching back phenomena.

MATERIALS AND METHODS

Sample fabrication and TR-MOKE measurements

The magnetic films of Ta (6 nm)/Co40Fe40B20 (0.8 nm)/MgO (2 nm)/SiO2 (3 nm) and Ta (3 nm)/Co40Fe40B20 (1.2 nm)/MgO (2 nm)/SiO2 (3 nm) were deposited on the thermally oxidized silicon substrates by magnetron sputtering with a base pressure of <2 × 10−9 torr at room temperature and patterned into a square of 3 μm × 3 μm or 3 μm × 6 μm. For a short pulse excitation, the ground-signal-ground (GSG) coplanar waveguide was patterned with electron-beam lithography and deposited with Ta (3 nm)/Cu (75 nm)/Ta (4.5 nm). In the stroboscopic pump-probe experiments, the pulse generator (pump) and the picosecond laser (probe) controller were synchronized by the pattern generator with a triggering frequency of 100 kHz, and each data point of TR-MOKE signal corresponded to an average of 60,000 events. The spot diameter of the laser beam was 2 μm. The reflected laser beam was measured by a balanced photodetector to obtain the MOKE signal. Measurements were carried out on three devices and showed similar results.

Micromagnetic simulations

Micromagnetic simulations were carried out by numerically solving Eq. 1 at zero temperature. The following parameters were used: Ms = 1.0 × 106 A m−1, exchange stiffness constant Aex = 1.0 × 10−11 J m−1, perpendicular anisotropy constant K = 0.9 × 106 J m−3, α = 0.02, c|| = −0.07, and c = +0.28. The sample dimension for Fig. 4 was 200 nm × 200 nm × 1 nm, and the unit cell size was 2 nm × 2 nm × 1 nm. The sample dimension for Fig. 5 was 2000 nm × 50 nm × 1 nm, and the unit cell size was 2 × 50 × 1 nm3. For the current pulse, both rise and fall times were 100 ps. In our sign convention, a negative DLT efficiency (c|| < 0) induced an up-to-down switching for J > 0 and Hx < 0. For two-dimensional micromagnetic simulations (Fig. 4), we introduced an artificial defect at the bottom-left corner and considered local demagnetization fields to mimic a domain wall nucleation at room temperature.

SUPPLEMENTARY MATERIALS

Supplementary material for this article is available at http://advances.sciencemag.org/cgi/content/full/3/4/e1603099/DC1

section S1. Stroboscopic pump-probe MOKE experiments using a picosecond laser

section S2. Characterization of magnetic films

section S3. Oscillatory SOT-induced magnetization switching: Micromagnetic simulations

section S4. Domain wall moving along the left edge

section S5. Backward motion of a domain wall reflected at an edge

fig. S1. Stroboscopic pump-probe MOKE setup.

fig. S2. Vibrating sample magnetometry and MOKE measurements.

fig. S3. Time-varying z component of the magnetization and its configurations.

fig. S4. Domain wall types formed in the two-dimensional sample.

fig. S5. Temporal evolutions of domain wall position q for the Bloch-type domain wall.

References (4250)

This is an open-access article distributed under the terms of the Creative Commons Attribution-NonCommercial license, which permits use, distribution, and reproduction in any medium, so long as the resultant use is not for commercial advantage and provided the original work is properly cited.

REFERENCES AND NOTES

Funding: This work was supported by the National Research Foundation (NRF), Prime Minister’s Office, Singapore, under its Competitive Research Programme (CRP award no. NRFCRP12-2013-01). K.-J.L. was supported by the NRF of Korea grant funded by the Korea government (MSIP) (2011–0027905, 2015M3D1A1070465). Author contributions: J.Y., J.H.K., and H.Y. initiated this work. J.Y. and X.Q. deposited the films. J.Y., J.H.K., and J.S. fabricated devices. J.Y. performed the switching measurements. J.M.L. characterized the films. S.-W.L., J.Y., and K.-J.L. performed theoretical and numerical studies. All authors discussed the results. J.Y., S.-W.L., K.-J.L., and H.Y. wrote the manuscript. H.Y. supervised and led the project. Competing interests: The authors declare that they have no competing interests. Data and materials availability: All data needed to evaluate the conclusions in the paper are present in the paper and/or the Supplementary Materials. Additional data related to this paper may be requested from the authors.
View Abstract

Related Content

Navigate This Article