Skip to main content

Thank you for visiting You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

Wigner solids of domain wall skyrmions


Detection and characterization of a different type of topological excitations, namely the domain wall (DW) skyrmion, has received increasing attention because the DW is ubiquitous from condensed matter to particle physics and cosmology. Here we present experimental evidence for the DW skyrmion as the ground state stabilized by long-range Coulomb interactions in a quantum Hall ferromagnet. We develop an alternative approach using nonlocal resistance measurements together with a local NMR probe to measure the effect of low current-induced dynamic nuclear polarization and thus to characterize the DW under equilibrium conditions. The dependence of nuclear spin relaxation in the DW on temperature, filling factor, quasiparticle localization, and effective magnetic fields allows us to interpret this ground state and its possible phase transitions in terms of Wigner solids of the DW skyrmion. These results demonstrate the importance of studying the intrinsic properties of quantum states that has been largely overlooked.


The skyrmion, a topological soliton solution for the description of hadrons in 1960s1, has attracted renewed interest in recent years due to its potential applications in information and communication technology demonstrated in magnetic materials2,3. The magnetic skyrmion is classified by a topological charge \(Q=\frac{1}{4\pi }\int {{{{{\bf{n}}}}}}\cdot \left({\partial }_{x}{{{{{\bf{n}}}}}}\times {\partial }_{y}{{{{{\bf{n}}}}}}\right){{{{{{\mathrm{d}}}}}}x{{{{{\mathrm{d}}}}}}y}\) counting the number of times that a unit vector \({{{{{\bf{n}}}}}}\) of the local spin wraps around the unit sphere in a spherically symmetric or combed hedgehog configuration. The symmetric hedgehog skyrmion has also been found in an isotopic (fully spin-polarized) quantum Hall ferromagnet (QHF) at filling factor ν = 1 of a two-dimensional electron gas (2DEG) in GaAs (called QH skyrmions)4,5: it is energetically cheaper to constitute this spin texture than to have single spin flips for the GaAs 2DEG with sufficiently low bare Zeeman energy determined by a small effective g factor of −0.44. There is an alternative possibility for the realization of skyrmions once low-lying excited states of the domain wall (DW) as a one-dimensional (1D) object have a nontrivial topological charge (referred to as DW skyrmions)6,7,8,9: the xy (z) component of a unit vector of spins gains a 2π (π) phase within the length (width) of this excitation along (across) the DW, defining a mapping topologically equivalent to the combed hedgehog sphere. The ubiquitous nature of the DW may hold promise for the creation and application of skyrmions in nonrelativistic and relativistic systems, which is of considerable importance from both theoretical and practical points of view.

More recently, the DW skyrmion has been demonstrated in chiral magnets with a strong Dzyaloshinskii–Moriya interaction induced by the spin–orbit coupling (SOC)10,11,12. Unlike the magnetic skyrmion that can form a lattice as the ground state, the magnetic DW skyrmion is shown to be low-lying spin excitations. In this study, we present data demonstrating the ground state of DW skyrmions in an Ising (easy-axis) QHF (IQHF). In contrast to the magnetic DW skyrmion, the DW skyrmion in the IQHF is expected to carry an electrical charge due to the identity between topological and electrical charge densities and to dominate the electronic properties of the IQHF as similar to the QH skyrmion6,13. Therefore, we combine nonlocal resistance measurements with a resistively detected NMR (RDNMR) technique (hereafter called NRDNMR) as a highly sensitive and minimally invasive approach to determine the DW skyrmion in the simplest IQHF at ν = 2 of an InSb 2DEG14,15. The nuclear spin relaxation (NSR) measurement, a generally accepted methodology for determining the low-frequency spin textures in the QH system (e.g., QH skyrmions)16, shows an increase in magnetic fluctuations in the DW coupling to the nuclei with decreasing temperature and increasing localization. This strongly suggests that the DW skyrmion may condense into a 1D Wigner crystal stabilized by long-range Coulomb repulsions, which displays both continuous and discontinuous phase transitions depending on the way that the effective magnetic field17 between two approaching Landau levels (LLs) with opposite spins changes.

Results and discussion

Nonlocal resistance measurements

The nonlocal resistance measurement (see “Methods” section and Supplementary Note 1) is known to play a unique role in the study of edge or surface states18,19. In this work, we demonstrate that the nonlocal measurement can be used to determine the intrinsic properties of bulk states. The bulk state under study here is the ν = 2 IQHF that is formed when the tilt angle θ (Fig. 1a) is tuned to bring the pseudospin-down [(n,σ) = (0,↓)] (where n and σ are the orbital and spin indices, respectively) and pseudospin-up (1,↑) LLs into degeneracy (Fig. 1b). It is shown in Fig. 1c that a resistance spike in longitudinal resistance RSD,12 occurs between two separate LL peaks (one locates at ~11 T and the other above the maximum field of 15 T in our measurements), signaling the presence of the IQHF20,21. Note that, for a simple level crossing, the two approaching LL peaks will merge into one14. The IQHF was also identified by the RDNMR measurement of the Knight shift at ν = 215 and by the calculation based on the Hartree–Fock theory (see below). Large Zeeman and cyclotron splittings22 effectively decouple the edge and bulk channels in the InSb 2DEG with relatively low mobility, making it possible to perform the nonlocal measurement. Although the ν = 2 IQHF is known to originate from dissipative transport through the DW channels6,21, its nonlocal transport properties have not been described yet. Figure 1c shows that the resistance of this IQHF consists of fine structures (FSs) in the nonlocal measurement, in contrast to that (RLL) of the LL peak outside the IQHF region. The FSs of two nonlocal configurations R23,14 and R14,23 are found to be different, however, they obey the reciprocity theorem \({R}_{{{{{{\rm{kl}}}}}},{{{{{\rm{mn}}}}}}}^{+B}={R}_{{{{{{\rm{mn}}}}}},{{{{{\rm{kl}}}}}}}^{-B}\) (Supplementary Fig. 3) derived from the Onsager relation23. Furthermore, the response of all FSs to temperature (T) is monotonic and in the opposite direction to that of RLL (Fig. 1d). The T dependence of the FSs is dominated by the variable range hopping (VRH) mechanism (Supplementary Fig. 4) as similar to that of the spike22, while the decrease in RLL at high temperatures originates from the phonon-assisted inter-edge-bulk scattering24. From the VRH formula, the B dependence of localization length ξ of the ν = 2 IQHF for different measurement configurations is plotted in Fig. 1e. Note that the edge states corresponding to the two intersecting LLs become part of an array of the (bulk) domains in the ν = 2 IQHF15, and in this case, the current is mainly carried by the edge channel related to the lowest LL with negligible scattering to the bulk state. Therefore, the FSs are believed to be determined by the bulk properties of the IQHF. Moreover, these FSs are different from nonlocal resistance fluctuations observed in the LL peak of quantum wires that are induced by the resonant tunneling between edge states through the bulk localized state25. Based on the bulk-edge model with an independent edge- and bulk-components (see ref. 18 and Supplementary Notes 1, 2), the four-terminal resistance is only parametrized by the longitudinal resistance of the bulk state and in this sense, all these resistances (e.g., R23,14, R14,23, and RSD,12) are essentially identical. However, the bulk current \({I}_{i}^{N}\) and the corresponding Hall electric-field εH in each segment i depend on the measurement configurations, accounting for the difference between the spike and its nonlocal counterpart. The presence (absence) of the FSs in the B-ξ dependence of R23,14 (RSD,12) with εH ~ 0.1 V m−1 (~14 V m−1) (Supplementary Note 2) in Fig. 1e indicates an electric field-dependent hopping transport; all FSs are broadened at large εH and thus merged into a single peak when the graded percolation problem is considered26. This conclusion is further supported by the data in Fig. 1f, where the FSs for R23,14 become a single peak as I23 increases up to 3.16 μA with εH ~ 10 V m1 and \({I}^{N}\) ~ 10 nA (where \({I}^{N}\) is the sum of \({I}_{i}^{N}\)). In this case, the hopping transport is still in equilibrium (Supplementary Fig. 5), and the segments are electrically linear as evidenced by the fact that RLL obeys the reciprocity theorem (Supplementary Fig. 6). However, the nonlocal resistance in the IQHF region does not satisfy the reciprocity theorem (Supplementary Fig. 6) because the Overhauser shift induced by dynamic nuclear polarization(DNP) breaks the Onsager relation locally. The effect of the DNP on nonlocal resistances was further investigated by the NRDNMR measurement (see “Methods” section).

Fig. 1: Nonlocal resistance measurements of the Ising quantum Hall ferromagnet (IQHF) at filling factor ν = 2.
figure 1

a Schematics of a two-dimensional electron gas (2DEG, dashed line) Hall bar with six terminals (1–4, S, and D) surrounded by a coil that produces a continuous-wave radiofrequency (RF) field. The tilt angle θ is defined between the magnetic field B and its perpendicular component Bperp. b Schematics of the Landau level (LL) splitting as a function of θ. The ν = 2 IQHF occurs when the two LLs with orbital index n = 0 and n =1 intersects, where the Zeeman splitting \({E}_{z} \sim {B}_{{{{{{\rm{perp}}}}}}}/{{\cos }}\theta\) and the cyclotron gap \({E}_{{{{{{\mathrm{c}}}}}}} \sim {B}_{{{{{{\rm{perp}}}}}}}\) are made equal by adjusting θ. An upward (downward) arrow is for spin-up (spin-down). c Four-terminal resistance \({R}_{{kl},{mn}}\) versus B at T = 1 K, \({I}_{{kl}}\) = 31.6 nA, and θ = 64°. The longitudinal resistance RSD,12 is given by the ratio of a voltage drop V12 between terminals 1 and 2 to current flow ISD between terminals S and D, and the nonlocal resistance R23,14 (R14,23) is calculated by V14/I23 (V23/I14). Note the different scales for longitudinal and nonlocal resistances. The shaded area represents the ν = 2 (B ~ 13 T, dashed line) IQHF region. d R23,14 versus B at different temperatures. The arrow indicates increasing temperature. e Localization length ξ versus B obtained from a fit (Supplementary Fig. 4) to the B dependence of R23,14 (red solid dots), R14,23 (blue solid dots), and RSD,12 (open dots) at different temperatures, respectively. f Current (I23) dependence of R23,14 versus B at T = 1 K.

The NRDNMR signals of the ν = 2 IQHF for \({R}_{23,14}^{\pm B}\) and \({R}_{14,23}^{\pm B}\) obtained at \({I}^{N}\) ~10 nA and the RDNMR signals of the ν = 2 IQHF for RSD,12 at \({I}^{N}\) ~ 500 nA are shown in Supplementary Figs. 8 and 9, respectively. Note that both RDNMR and NRDNMR signals are expected to be caused by the current-induced DNP14,27 responsible for the polarization of nuclei in the ν = 2/3 QHF of GaAs 2DEGs, but the underlying mechanism is still unclear15,16. A shift of the ν = 2 spike induced by the DNP with the degree of nuclear polarization of ~10% near the DW is assigned to cause the resistance change in the RDNMR measurement of the InSb 2DEG15. Relatively low nonlocal resistance (Supplementary Note 2) responds more sensitively to subtle changes in the sample, thus accounting for the DNP performed by very low current in the NRDNMR measurement and for increased detection sensitivity: the maximum signal amplitude for \({R}_{23,14}^{\pm B}\) and \({R}_{14,23}^{\pm B}\) is ten times larger than that for RSD,12. Furthermore, the NRDNMR signals for \({R}_{23,14}^{+B}\) and \({R}_{14,23}^{-B}\) are found to spread out over the ν = 2 IQHF region but are sparsely distributed for \({R}_{23,14}^{-B}\) and \({R}_{14,23}^{+B}\). This suggests that the DW structure in which the DNP occurs depends on the measurement configurations, which is further supported by measurements of nuclear spin-lattice relaxation time T1 and spin dephasing time T2 (see “Methods” section) as local probes of the low-frequency spin dynamics of the DW.

NRDNMR relaxation time measurements

We now discuss the T1 and T2 results of the ν = 2 IQHF for \({R}_{14,23}^{+B}\). Figure 2a shows that T1 near ν = 2 is very short (~8 s) and gradually increases until \(\left|2-\nu \right| \sim 0.023\) (dashed line). Then T1 increases sharply on the low (high)-ν side at 0.3 K (1 K) and decreases with decreasing ν (data on the high-ν side are not available). Details of the T-dependent T1 are presented in Fig. 2b in terms of the nuclear spin-lattice relaxation rate 1/T1. It is clear that the dependence of 1/T1 on T is linear at \(\left|2-\nu \right| < 0.023\) (called phase I, Fig. 2a) but following an Arrhenius-like behavior at \(\left|2-\nu \right| > 0.023\) (phase II), demonstrating that different relaxation mechanisms are responsible for the two phases. Fast 1/T1 at T → 0 in phase I is indicative of strong magnetic fluctuations coupling to the nuclei. Note that T1 in the ν = 2 IQHF obtained from the RDNMR measurement is independent of temperature as also observed in both the ν = 2/3 QHF and the two-subband QHF of the GaAs 2DEG28. All these results cannot be explained by a simple level crossing with single spin flips taking account of disorder and electron exchange interaction that accounts for a slow NSR dominated by the Korringa law: T1 T = const28. Furthermore, \({R}_{{{{{\mathrm{14,23}}}}}}^{+B}\) in phase I is found to be linear in T at large current for the NRDNMR measurement (inset, Fig. 2c), resulting in a linear dependence of 1/T1 on \({R}_{{{{{\mathrm{14,23}}}}}}^{+B}\) (Fig. 2c). Because the localization of quasiparticles in the DW is expected to account for less conducting DWs, an increase in 1/T1 with decreasing \({R}_{14,23}^{+B}\) suggests that quasiparticle localization rather than additional conducting states dominates the NSR in the DW. This is opposite to the Korringa relaxation where 1/T1 is proportional to the resistance of the QH state29. We point out that these findings are analogous to those of the ν ≈ 1 QH state in an extremely high-quality GaAs 2DEG as evidence for the formation of 2D skyrmion crystallization (SC)30, a Wigner crystal with lattice points occupied by charged skyrmionic spin textures31. The DW in the QHF is expected to have electron spins noncollinear to B in the excited states including the spin wave (SW) and the DW skyrmion6. The thermally-activated transport measurement of both the ν = 2/3 QHF and the two-subband QHF in the GaAs 2DEG suggests that low-lying excitations akin to the DW skyrmion may favor the NSR, but intricate composite-fermion interactions at ν = 2/3 and the subband degree of freedom complicate the interpretation28. Below we show that these excitations are necessarily related to the interpretation and discussion of our data.

Fig. 2: Temperature and magnetic field dependence of T1 and T2 in the NRDNMR measurement of the Ising quantum Hall ferromagnet (IQHF) at ν = 2.
figure 2

a \({R}_{{{{{\mathrm{14,23}}}}}}^{+B}\) (right y axis) as a function of B (or ν) swept upwards (solid line) and downwards (dash-dotted line) with a sweep rate of 1.7 mT s−1 at I14 = 3.16 μA, θ = 64°, and T = 1 K. The B oriented parallel (opposite) to the normal direction of the 2DEG is written as \(+\left(-\right)\)B. \(\triangle {R}_{{{{{{\rm{hys}}}}}}}\) measures the magnitude of the resistance hysteresis at a fixed B. The NRDNMR signal occurs around ν = 2 where the B (or ν) dependence of T1 (left y axis) is obtained. The dashed line marks the phase boundary (see main text). Note that error bars for all data in this study are the standard deviation of the mean. b 1/T1 versus T at different B. The solid line is a fit to the data at 12.82 T (13.18 T) using the Arrhenius law with an activation energy \({\triangle }_{{{{{{\rm{sw}}}}}}}\) of 1.36 (0.42) meV plus a T-independent constant of 0.02 s1. The dashed line is a linear fit, and its extrapolation to the x axis defines Tc at which 1/T1 is zero. Inset shows the ν dependence of Tc. c 1/T1 versus \({R}_{{{{{\mathrm{14,23}}}}}}^{+B}\) obtained by comparing the T dependence of \({R}_{{{{{\mathrm{14,23}}}}}}^{+B}\) at B = 13 T and I14 = 3.16 μA (inset) and that of 1/T1 (data at B = 13 T, b) point by point. The solid line is a guide to the eye. d \({\omega }_{{{{{{\rm{R}}}}}}}^{2}{T}_{2}\) versus B (or ν) at different temperatures. The dashed line marks the phase boundary that is the same as that in a. The Rabi frequency \({\omega }_{{{{{{\mathrm{R}}}}}}}=\gamma {B}_{{{{{{\rm{RF}}}}}}}/2\) (where \(\gamma\) = 9.36 MHz/T is the gyromagnetic ratio of 115In) is about 10 Hz (see “Methods” section). Inset shows \({\omega }_{{{{{{\rm{R}}}}}}}^{2}{T}_{2}\) as a function of BRF tuned by the RF output power. A good fit (\(\propto {B}_{{{{{{{\mathrm{RF}}}}}}}}^{2}\), solid line) to the data indicates a linear dependence of \({\omega }_{{{{{{\rm{R}}}}}}}\) on BRF.

The SOC in the IQHF treated as a perturbation due to its off-diagonal entries in the LL basis is expected to twist the planar (XY) component of \({{{{{\bf{n}}}}}}\) with a 360° rotation along the DW, which lifts the degeneracy of the DW ground state and forms the DW skyrmion excitation6. The energy of the DW skyrmion is determined by the SOC strength in the form of \({{{{{\rm{\alpha }}}}}}={\epsilon }_{{{{{{{\mathrm{so}}}}}}}}/\hslash {\omega }_{{{{{{\mathrm{c}}}}}}}\) [see ref. 6 and parameters given in Supplementary Table 1], in contrast to the case of QH skyrmions whose energy is determined by the effective g factor in the form of the Zeeman energy \({E}_{z}\). Note that the QH skyrmion cannot be available in the InSb 2DEG where a relatively large g factor of −39 to −8832 makes the ratio of the Zeeman and Coulomb energies 100 times larger than that required for its formation33. It is shown in Supplementary Table 1 that \({{{{{\rm{\alpha }}}}}}\) in the InSb 2DEG is only 3.8 times larger than that in the GaAs 2DEG while the difference in \({E}_{z}\) between these two 2DEGs is significant (over 100 times). The DW skyrmion in the IQHF is able to carry electric charges similar to the QH skyrmion13. The ν = 2 IQHF occurs at zero effective field as discussed later, where dissipative transport in 1D DW channels percolating through the sample is responsible for \({R}_{14,23}^{+B}\). We expect that single DW skyrmions with long-range Coulomb repulsions will form a crystalline state in the 1D DW channels with quasi-long-range order analogous to the 1D Wigner crystal in carbon nanotubes34,35, where a gapless XY spin wave mode with respect to the broken SO(2) symmetry would relax the nuclear spins more efficiently36,37 as shown by a short T1 at ν = 2 (Fig. 2a). Away from ν = 2 a decrease in the number of DW skyrmions due to a nonzero b* 21,38 (positive for ν > 2 and negative for ν < 2, Supplementary Fig. 10) makes the crystal harder to form. A crystal of the DW skyrmion with low density is expected to melt at low temperatures because it is easier to break their bonds, as indicated in the inset of Fig. 2b where the state of phase I is found to be thermodynamically most stable up to ~4 K at ν = 2. The NSR in less crystalline states is partially suppressed and therefore T1 becomes longer. When moving further into the region of \(\left|2-\nu \right| \; > \; 0.023\) (phase II), a sudden change in T1 signals a discontinuous transition to what we interpret as the SW-mediated NSR in neutral DWs with Q = 0. The Arrhenius-like behavior of 1/T1 in phase II (Fig. 2b) indicates that the NSR occurs in low-energy excited states separated from the ground state by an energy gap \({\triangle }_{{{{{{\rm{sw}}}}}}}\). It is predicted that the neutral DW is preferably supported by a nonzero \(b\)* 21, where symmetry breaking induced by the SO interaction leads to the opening of a gap in the otherwise gapless spectrum of SW excitations6,39. The number of neutral DWs increases with increasing \(\left|b\right|\), resulting in a reduced \({\triangle }_{{{{{{\rm{sw}}}}}}}\) by the SW–SW interaction40 (Supplementary Fig. 11) that accounts for a drop in T1 in phase II (Fig. 2a). Because of the difference in the NSR mode for phase I and phase II, the transformation between them is discontinuous that can be viewed as a first-order phase transition. This is responsible for a jump in T1 between the two phases at low T (Fig. 2a). As T increases, the jump is smeared out (Fig. 2a) due to a decrease in T1 in phase II caused by the Arrhenius behavior of thermal activation. Figure 2d shows the T and \(\nu\) dependence of T2 that is known to be dominated by dynamic fluctuating fields experienced by the nuclei in the quantum well associated with the dynamics of the 2DEG41. It is shown that T2 in phase I decreases with increasing T in contrast to the motional narrowing effect, suggesting that the dynamics of the DW skyrmion is in the frozen limit where a very small number of free DW skyrmions suppresses the motion effect42. The 1D lattice of localized DW skyrmions may be responsible for the frozen limit that is taken as evidence for the localization of skyrmions in the 2D lattice41. In addition, T2 drops off quickly with increasing \(\left|2-\nu \right|\) due to relatively strong nuclear dephasing in less crystalline states.

Figure 3a shows a contour plot of \(\triangle {R}_{{{{{{\rm{hys}}}}}}}/{R}_{{{{{\mathrm{14,23}}}}}}^{+B}\) as a function of B and θ. The position of maximum \(\triangle {R}_{{{{{{\rm{hys}}}}}}}/{R}_{{{{{\mathrm{14,23}}}}}}^{+B}\) is indicated by solid dots, where both 1/T1 and \(\triangle {R}_{14,23}/{R}_{14.23}^{{{{{{\rm{sat}}}}}}}\) are also at a maximum. The 1D DW channel percolating through the sample is expected to occur at these positions, in which more DW skyrmions are energetically favored. Based on many-body Hartree–Fock theory17, we have made a calculation of \({b}^{\ast }=0\) in the B-Bperp plane (Supplementary Fig. 10) at which the IQHF occurs. It is seen that there is indeed an association between the theoretical \({b}^{\ast }=0\) line and data points at these positions.

Fig. 3: T1 and T2 results of the quantum Hall ferromagnet at ν = 2 with zero effective field b*.
figure 3

a Contour plot of B and \(\theta\) dependence of \(\triangle {R}_{{{{{{\rm{hys}}}}}}}/{R}_{{{{{\mathrm{14,23}}}}}}^{+B}\) at a field-sweep rate of 1.7 mT s−1, I14 = 3.16 μA, and T = 1 K. Red dots indicate the position of maximum \(\triangle {R}_{{{{{{\rm{hys}}}}}}}/{R}_{{{{{\mathrm{14,23}}}}}}^{+B}\), 1/T1, and \(\triangle {R}_{14,23}/{R}_{14.23}^{{{{{{\rm{sat}}}}}}}\). The solid line is a guide to the eye and the dashed line is calculated for the B and \(\theta\) dependence of ν = 2. Temperature dependence of 1/T1 (b) and \({\omega }_{{{{{{\rm{R}}}}}}}^{2}{T}_{2}\) (c) versus ν where the data (symbols) are in one-to-one correspondence with the red dots in a by \(\nu ={n}_{{{{{{\rm{s}}}}}}}h/{eB}{{\cos }}\theta\) (where h is Planck’s constant and e the electron charge). It is clear that the data encircled by an oval in a exhibit a strong temperature dependence of T1 and T2.

However, the position of these data points is different from that of maximum \({R}_{{{{{\mathrm{14,23}}}}}}^{+B}\) (Supplementary Fig. 12) that is generally believed to be associated with \({b}^{\ast }=0\)21. Our results demonstrate that the resistance measurement cannot establish this association in a straightforward manner. The T dependence of 1/T1 and T2 at these positions in terms of \(\nu\) is plotted in Fig. 3b, c. Note that these results were obtained along the \({b}^{\ast }=0\) line rather than away from \({b}^{\ast }=0\) (Fig. 2). Although the Wigner crystal of the DW skyrmion occurs at \(\nu =2\) with \({b}^{\ast }=0\), it changes along the \({b}^{\ast }=0\) line as shown by a nonmonotonic filling dependence of 1/T1 and T2 with a peak around \(\nu =2\). This nonmonotonic behavior clearly rules out the possible role of exchange interactions at the Ising transition in 1/T1 and T2 because the exchange-correlation energy varies monotonously along the \({b}^{\ast }=0\) line20. Away from \(\nu =2\) backscattering due to a finite density of extended states increases as indicated by an increase in resistance (Supplementary Fig. 13), which effectively induces random potential fluctuations to localize the eigenstates of the DW skyrmion and therefore makes the crystal harder to form. Although the T and \(\nu\) dependence of 1/T1 and T2 for these less crystalline states are similar to that in Fig. 2, the \(\nu\) dependence is highly asymmetric because of the asymmetry of backscattering around \(\nu =2\) as indicated in Supplementary Fig. 13. Deep inside the strong backscattering region the DW skyrmions might have a clumped distribution between the potential wells behaving like a Wigner glass43, which breaks the quasi-long-range order of the 1D Wigner crystal. The NSR in this region is probably caused by fluctuations of low-lying electronic states with energies smaller than \({k}_{{{{{{\mathrm{B}}}}}}}T\) (where \({k}_{{{{{{\mathrm{B}}}}}}}\) is the Boltzmann constant), exhibiting the temperature-independent 1/T144. It is also seen that T2 in the glassy regime is very short and temperature-independent, indicating that more inhomogeneous local field fluctuations prevail over the motion effect for the nuclear dephasing. Note that the T1 and T2 properties of the glassy phase are the same as those of the \(\nu =2\) IQHF for other measurement configurations and also as those of the \(\nu =2\) IQHF obtained from the RDNMR measurement28, suggesting that the presence of 1D Wigner glass of the DW skyrmion is responsible for the NSR therein. Furthermore, the T1 property specified in the DW Wigner crystal and glass is qualitatively consistent with that predicted in the renormalized classical and quantum critical phases of the skyrmion lattice45,46, respectively, suggesting that the Wigner solid of the DW skyrmion is able to provide a rich test bed for experimental investigations and theoretical descriptions of quantum criticality.

Note that the width of the DW skyrmion is determined by the DW thickness of the order of magnetic length \({l}_{{{{{{\mathrm{B}}}}}}}=25.6{{{{{\rm{nm}}}}}}/\sqrt{B{\cos }\theta }\) (~11 nm at \(\nu =2\))15 and its length strongly depends on the SOC strength13. The formation of domain structure (several hundreds of nm in domain size15) in the InSb 2DEG with relatively low mobility does not allow the DW skyrmion to condense into the Wigner solid in 2D such as the 2D SC observed in a GaAs 2DEG with extremely high mobility30. As discussed above, the electrons in the 1D DW channels percolating through the sample move in a correlated fashion, which is expected to dominate the dissipative transport and thus the NSR in the \(\nu =2\) IQHF. It should be pointed out that this 1D skyrmionic Wigner crystal cannot be obtained by laterally confining the 2D SC because of the confinement-induced melting effect47. It is also impossible for the skyrmions in 1D nanowires48 and for the magnetic DW skyrmions11,12 to condense into a Wigner solid due to the absence of the long-range nature of their interactions. The 1D Wigner crystal consisting of topological solitons instead of electrons is thermally stable up to ~4 K that will be further increased by optimizing the size of the DW skyrmions and thus their interactions, providing a platform where the Wigner crystal can be formed at high temperatures. We note that the 1D skyrmionic Wigner crystal discovered in this study is only distributed over a very narrow range of filling factors (\(\left|2-\nu \right| < 0.023\)) in a certain measurement configuration with nonlocal resistance R23,14 but quite robust to condition changes. Crucial to the success of this discovery is the development of the NRDNMR technique that makes it possible to construct and characterize the DW with different configurations under equilibrium conditions. Our results highlight the uniqueness of the NRDNMR technique in the determination of intrinsic properties of quantum states and of their phase transition. Results of the present study strongly suggest that attention should be paid to previous transport and NMR measurements, where the role played by the Hall electric field has tended to be overlooked.


Nonlocal resistance measurement

The sample of 2DEGs in a 20-nm-wide InSb quantum well under study here was patterned into a Hall bar (100 μm length and 30 μm width) with Ti/Au as Ohmic contacts. A low-noise preamplifier (Stanford Research Systems, Model SR560) and an alternating current (AC) resistance bridge (Lake Shore Model 370) at 13.7 Hz were used for the direct current (DC) and AC measurements, respectively. It is shown in Supplementary Fig. 2 that the AC nonlocal resistance \({R}_{{{{{\mathrm{23,14}}}}}}\) of the ν = 2 IQHF and its DC counterpart \({R}_{{{{{\mathrm{23,14}}}}}}^{{{{{{\rm{DC}}}}}}+}\) (\({R}_{{{{{\mathrm{23,14}}}}}}^{{{{{{\rm{DC}}}}}}-}\)) for positive (negative) current satisfy \({R}_{{{{{\mathrm{23,14}}}}}}\approx \frac{1}{2}\left({R}_{{{{{\mathrm{23,14}}}}}}^{{{{{{\rm{DC}}}}}}+}{+R}_{{{{{\mathrm{23,14}}}}}}^{{{{{{\rm{DC}}}}}}-}\right)\), suggesting that the FSs are independent of the type of electric current. Because the signal-to-noise ratio of the DC NRDNMR measurement is relatively low, all other data shown in this manuscript were collected using the AC resistance bridge in a dilution refrigerator with in situ rotator stage. Low-temperature electron density of ns ~ 2.76 × 1015 m−2 and mobility of μ ~ 19.3 m2 Vs−1 of the InSb 2DEG were obtained from the fast Fourier transform (FFT) analysis of low-field Shubnikov-de-Haas (SdH) oscillations and from the value of RSD,12 at B = 0, respectively.

NRDNMR, T 1, and T 2 measurements

The details of the NRDNMR measurement are shown in Supplementary Fig. 7: a relatively large current is applied to polarize the nuclei in the ν = 2 IQHF region, as indicated by an exponential increase in the nonlocal resistance (for example, R14,23) on a time scale of hundreds of seconds. After R14,23 becomes saturated, a continuous-wave RF field at a resonance (Larmor) frequency of \({f}_{{{{{{{\mathrm{NMR}}}}}}}}=\gamma B\) is applied to irradiate the 2DEG (on resonance). The RF field is generated by a single coil with 2 turns (a cross-sectional area of 3 mm × 8 mm and a length of 4 mm) surrounding the sample connected to the RF generator using a 50 Ω coax cable. The change in R14,23 representing the depolarization of nuclei reaches a maximum ΔRnl at complete saturation (\({R}_{{{{{{\rm{nl}}}}}}}^{{{{{{\rm{sat}}}}}}}\)) with the rise time \({T}_{{{{{{\rm{r}}}}}}}={T}_{1}/\left(1+{\omega }_{{{{{{\rm{R}}}}}}}^{2}{T}_{1}{T}_{2}\right)\), where \({\omega }_{{{{{{\rm{R}}}}}}}\) is the Rabi frequency. The component of the RF field perpendicular to B (\({B}_{{{{{{\rm{RF}}}}}}}\)) is estimated to be ~2.2 μT for an RF output power of 0 dbm at θ = 64° that gives \({\omega }_{{{{{{\mathrm{R}}}}}}} \sim\)10 Hz. As the RF frequency is detuned from \({f}_{{{{{{\rm{NMR}}}}}}}\) (off resonance), R14,23 is decreased exponentially due to the repolarization of nuclei. The fall time of repolarization gives \({T}_{1}\) and thus T2 is obtained from Tr.

Data availability

The authors declare that data supporting the findings of this study are available within the paper and its Supplementary Information files.


  1. 1.

    Skyrme, T. H. R. A unified field theory of mesons and baryons. Nucl. Phys. 31, 556–569 (1962).

    MathSciNet  CAS  Article  Google Scholar 

  2. 2.

    Nagaosa, N. & Tokura, Y. Topological properties and dynamics of magnetic skyrmions. Nat. Nanotechnol. 8, 899–911 (2013).

    ADS  CAS  PubMed  Article  PubMed Central  Google Scholar 

  3. 3.

    Fert, A., Reyren, N. & Cros, V. Magnetic skyrmions: advances in physics and potential applications. Nat. Rev. Mater. 2, 17031 (2017).

    ADS  CAS  Article  Google Scholar 

  4. 4.

    Sondhi, S. L., Karlhede, A., Kivelson, S. A. & Rezayi, E. H. Skyrmions and the crossover from the integer to fractional quantum Hall effect at small Zeeman energies. Phys. Rev. B 47, 16419–16426 (1993).

    ADS  CAS  Article  Google Scholar 

  5. 5.

    Barrett, S. E., Dabbagh, G., Pfeiffer, L. N., West, K. W. & Tycko, R. Optically pumped NMR evidence for finite-size skyrmions in GaAs quantum wells near Landau level filling v = 1. Phys. Rev. Lett. 74, 5112–5115 (1995).

    ADS  CAS  PubMed  Article  Google Scholar 

  6. 6.

    Fal’ko, V. I. & Iordanskii, S. V. Topological defects and Goldstone excitations in domain walls between ferromagnetic quantum Hall liquids. Phys. Rev. Lett. 82, 402–405 (1999).

    ADS  Article  Google Scholar 

  7. 7.

    Nitta, M. Matryoshka Skyrmions. Nucl. Phys. B 872, 62–71 (2013).

    ADS  MathSciNet  CAS  MATH  Article  Google Scholar 

  8. 8.

    Jennings, P. & Sutcliffe, P. The dynamics of domain wall Skyrmions. J. Phys. A Math. Theor. 46, 465401 (2013).

    ADS  MathSciNet  MATH  Article  Google Scholar 

  9. 9.

    Danon, J., Balram, A. C., Sánchez, S. & Rudner, M. S. Charge and spin textures of Ising quantum Hall ferromagnet domain walls. Phys. Rev. B 100, 235406 (2019).

    ADS  CAS  Article  Google Scholar 

  10. 10.

    Cheng, R. et al. Magnetic domain wall skyrmions. Phys. Rev. B 99, 184412 (2019).

    ADS  CAS  Article  Google Scholar 

  11. 11.

    Li, M. et al. Magnetic domain wall structures in Pt/Co/Ni/Ir multi-layers. Preprint at (2020).

  12. 12.

    Nagase, T. et al. Observation of domain wall bimerons in chiral magnets. Nat. Commun. 12, 3490 (2021).

    ADS  CAS  PubMed  PubMed Central  Article  Google Scholar 

  13. 13.

    Brey, L. & Tejedor, C. Spins, charges, and currents at domain walls in a quantum Hall Ising ferromagnet. Phys. Rev. B 66, 041308(R) (2002).

    ADS  Article  CAS  Google Scholar 

  14. 14.

    Liu, H. W., Yang, K. F., Mishima, T. D., Santos, M. B. & Hirayama, Y. Dynamic nuclear polarization and nuclear magnetic resonance in the simplest pseudospin quantum Hall ferromagnet. Phys. Rev. B 82, 241304(R) (2010).

    ADS  Article  CAS  Google Scholar 

  15. 15.

    Yang, K. F. et al. Role of chiral quantum Hall edge states in nuclear spin polarization. Nat. Commun. 8, 15084 (2017).

    ADS  CAS  PubMed  PubMed Central  Article  Google Scholar 

  16. 16.

    Hirayama, Y. et al. Electron-spin/nuclear-spin interactions and NMR in semiconductors. Semicond. Sci. Technol. 24, 023001 (2009).

    ADS  Article  CAS  Google Scholar 

  17. 17.

    Jungwirth, T. & MacDonald, A. H. Pseudospin anisotropy classification of quantum Hall ferromagnets. Phys. Rev. B 63, 035305 (2000).

    ADS  Article  CAS  Google Scholar 

  18. 18.

    McEuen, P. L. et al. New resistivity for high-mobility quantum Hall conductors. Phys. Rev. Lett. 64, 2062–2065 (1990).

    ADS  CAS  PubMed  Article  Google Scholar 

  19. 19.

    Abanin, D. A. et al. Giant nonlocality near the Dirac point in graphene. Science 332, 328–330 (2011).

    ADS  CAS  PubMed  Article  Google Scholar 

  20. 20.

    De Poortere, E. P., Tutuc, E., Papadakis, S. J. & Shayegan, M. Resistance spikes at transitions between quantum Hall ferromagnets. Science 290, 1546–1549 (2000).

    ADS  PubMed  Article  Google Scholar 

  21. 21.

    Jungwirth, T. & MacDonald, A. H. Resistance spikes and domain wall loops in Ising quantum Hall ferromagnets. Phys. Rev. Lett. 87, 216801 (2001).

    ADS  CAS  PubMed  Article  Google Scholar 

  22. 22.

    Yang, K. F. et al. Resistively detected nuclear magnetic resonance via a single InSb two-dimensional electron gas at high temperature. Appl. Phys. Lett. 98, 142109 (2011).

    ADS  Article  CAS  Google Scholar 

  23. 23.

    Büttiker, M. Four-terminal phase-coherent conductance. Phys. Rev. Lett. 57, 1761–1764 (1986).

    ADS  PubMed  Article  Google Scholar 

  24. 24.

    Wang, J. K. & Goldman, V. J. Measurements and modeling of nonlocal resistance in the fractional quantum Hall effect. Phys. Rev. B 45, 13479–13487 (1992).

    ADS  CAS  Article  Google Scholar 

  25. 25.

    Main, P. C. et al. Resistance fluctuations in the quantum Hall regime. Phys. Rev. B 50, 4450–4455 (1994).

    ADS  CAS  Article  Google Scholar 

  26. 26.

    Trugman, S. A. Localization, percolation, and the quantum Hall effect. Phys. Rev. B 27, 7539–7546 (1983).

    ADS  Article  Google Scholar 

  27. 27.

    Korkusinski, M., Hawrylak, P., Liu, H. W. & Hirayama, Y. Manipulation of a nuclear spin by a magnetic domain wall in a quantum Hall ferromagnet. Sci. Rep. 7, 43553 (2017).

    ADS  CAS  PubMed  PubMed Central  Article  Google Scholar 

  28. 28.

    Yang, K. F. et al. Pump-probe nuclear spin relaxation study of the quantum Hall ferromagnet at filling factor ν = 2. New J. Phys. 21, 083004 (2019).

    ADS  CAS  Article  Google Scholar 

  29. 29.

    Desrat, W. et al. Resistively detected nuclear resonance in the quantum Hall regime: possible evidence for a Skyrme crystal. Phys. Rev. Lett. 88, 256807 (2002).

    ADS  CAS  PubMed  Article  Google Scholar 

  30. 30.

    Gervais, G. et al. Evidence for skyrmions crystallization from NMR relaxation experiments. Phys. Rev. Lett. 94, 196803 (2005).

    ADS  CAS  PubMed  Article  PubMed Central  Google Scholar 

  31. 31.

    Zhu, H. et al. Pinning-mode resonance of a Skyrme crystal near Landau-level filling factor ν = 1. Phys. Rev. Lett. 104, 226801 (2010).

    ADS  PubMed  Article  CAS  PubMed Central  Google Scholar 

  32. 32.

    Yang, K. F. et al. Nonlinear magnetic field dependence of spin polarization in high-density two-dimensional electron systems. New J. Phys. 13, 083010 (2011).

    ADS  Article  CAS  Google Scholar 

  33. 33.

    Fertig, H. A., Brey, L., Côté, R. & MacDonald, A. H. Charged spin-texture excitations and the Hartree-Fock approximation in the quantum Hall effect. Phys. Rev. B 50, 11018–11021 (1994).

    ADS  CAS  Article  Google Scholar 

  34. 34.

    Schulz, H. J. Wigner crystal in one dimension. Phys. Rev. Lett. 71, 1864–1867 (1993).

    ADS  CAS  PubMed  Article  PubMed Central  Google Scholar 

  35. 35.

    Shapir, I. et al. Imaging the electronics Wigner crystal in one dimension. Science 364, 870–875 (2019).

    ADS  CAS  PubMed  Article  PubMed Central  Google Scholar 

  36. 36.

    Côté, R. et al. Collective exciations, NMR, and phase transitions in Skyrme crystals. Phys. Rev. Lett. 78, 4825–4828 (1997).

    ADS  Article  Google Scholar 

  37. 37.

    Fertig, H. A. & Brey, L. Luttinger liquid at the edge of undoped graphene in a strong magnetic field. Phys. Rev. Lett. 97, 116805 (2006).

    ADS  CAS  PubMed  Article  Google Scholar 

  38. 38.

    Fal’ko, V. I. & Iordanskii, S. V. Spin-orbit coupling effect on quantum Hall ferromagnets with vanishing Zeeman energy. Phys. Rev. Lett. 84, 127–130 (2000).

    ADS  PubMed  Article  Google Scholar 

  39. 39.

    Mitra, A. & Girvin, S. M. Electron/nuclear spin domain walls in quantum Hall systems. Phys. Rev. B 67, 245311 (2003).

    ADS  Article  CAS  Google Scholar 

  40. 40.

    Kasner, M., Palacios, J. J. & MacDonald, A. H. Quasiparticle properties of quantum Hall ferromagnets. Phys. Rev. B 62, 2640–2658 (2000).

    ADS  CAS  Article  Google Scholar 

  41. 41.

    Khandelwal, P. et al. Spectroscopic evidence for the localization of skyrmions near ν = 1 as T→0. Phys. Rev. Lett. 86, 5353–5356 (2001).

    ADS  CAS  PubMed  Article  Google Scholar 

  42. 42.

    Villares Ferrer, A., Doretto, R. L. & Caldeira, A. O. NMR linewidth and skyrmion localization in quantum Hall ferromagnets. Phys. Rev. B 70, 045319 (2004).

    ADS  Article  CAS  Google Scholar 

  43. 43.

    Akhanjee, S. & Rudnick, J. Disorder induced transition into a one-dimensional Wigner glass. Phys. Rev. Lett. 99, 236403 (2007).

    ADS  PubMed  Article  CAS  Google Scholar 

  44. 44.

    Tycko, R., Barrett, S. E., Dabbagh, G., Pfeiffer, L. N. & West, K. W. Electronics states in gallium arsenide quantum wells probed by optically pumped NMR. Science 268, 1460–1463 (1995).

    ADS  CAS  PubMed  Article  Google Scholar 

  45. 45.

    Read, N. & Sachdev, S. Continuum quantum ferromagnets at finite temperature and the quantum Hall effect. Phys. Rev. Lett. 75, 3509–3512 (1995).

    ADS  CAS  PubMed  Article  Google Scholar 

  46. 46.

    Green, A. G. Quantum-critical dynamics of the Skyrmion lattice. Phys. Rev. B 61, R16299–R16302 (2000).

    ADS  CAS  Article  Google Scholar 

  47. 47.

    Kobayashi, T., Kumada, N., Ota, T., Sasaki, S. & Hirayama, Y. Low-frequency spin fluctuations in skyrmions confined by wires: measurement of local nuclear spin relaxation. Phys. Rev. Lett. 107, 126807 (2011).

    ADS  PubMed  Article  CAS  Google Scholar 

  48. 48.

    Lin, S. Z., Reichhardt, C., Batista, C. D. & Saxena, A. Particle model for skyrmions in metallic chiral magnets: dynamics, pinning, and creep. Phys. Rev. B 87, 214419 (2013).

    ADS  Article  CAS  Google Scholar 

Download references


We thank G. Yusa, K. Muraki, and N. Kumada for helpful discussions. This work was supported by the National Natural Science Foundation of China (11974132 and 11704144), the Jilin Natural Science Foundation (20180101286JC), the Fundamental Research Funds for the Central Universities, the JST-ERATO, the KAKENHI (15H05867 and 18H01811), and the CSRN (Y.H.) and the GP-Spin (K.N. and Y.H.) in Tohoku University.

Author information




H.L. designed and supervised the research. K.Y. performed the experiments and collected the data. K.N. and Y.H. fabricated InSb/AlInSb Hall bars. T.D.M. and M.B.S. grew InSb heterostructures. H.L. and K.Y. analyzed the data. H.L. wrote the manuscript.

Corresponding author

Correspondence to Hongwu Liu.

Ethics declarations

Competing interests

The authors declare no competing interests.

Additional information

Peer review information Nature Communications thanks Vesna Mitrovic and the other, anonymous, reviewer(s) for their contribution to the peer review of this work. Peer reviewer reports are available.

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Yang, K., Nagase, K., Hirayama, Y. et al. Wigner solids of domain wall skyrmions. Nat Commun 12, 6006 (2021).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.


Quick links

Nature Briefing

Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.

Get the most important science stories of the day, free in your inbox. Sign up for Nature Briefing