## Abstract

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.

## Introduction

The skyrmion, a topological soliton solution for the description of hadrons in 1960s^{1}, has attracted renewed interest in recent years due to its potential applications in information and communication technology demonstrated in magnetic materials^{2,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 skyrmion^{6,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 2DEG^{14,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 field^{17} 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 states^{18,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 *R*_{SD,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 IQHF^{20,21}. Note that, for a simple level crossing, the two approaching LL peaks will merge into one^{14}. The IQHF was also identified by the RDNMR measurement of the Knight shift at *ν* = 2^{15} and by the calculation based on the Hartree–Fock theory (see below). Large Zeeman and cyclotron splittings^{22} 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 channels^{6,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 (*R*_{LL}) of the LL peak outside the IQHF region. The FSs of two nonlocal configurations *R*_{23,14} and *R*_{14,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 relation^{23}. Furthermore, the response of all FSs to temperature (*T*) is monotonic and in the opposite direction to that of *R*_{LL} (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 spike^{22}, while the decrease in *R*_{LL} at high temperatures originates from the phonon-assisted inter-edge-bulk scattering^{24}. 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 IQHF^{15}, 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 state^{25}. 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., *R*_{23,14}, *R*_{14,23}, and *R*_{SD,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 *R*_{23,14} (*R*_{SD,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 considered^{26}. This conclusion is further supported by the data in Fig. 1f, where the FSs for *R*_{23,14} become a single peak as *I*_{23} increases up to 3.16 μA with *ε*_{H} ~ 10 V m^{−}^{1} 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 *R*_{LL} 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).

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 *R*_{SD,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 DNP^{14,27} responsible for the polarization of nuclei in the *ν* = 2/3 QHF of GaAs 2DEGs, but the underlying mechanism is still unclear^{15,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 2DEG^{15}. 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 *R*_{SD,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 *T*_{1} and spin dephasing time *T*_{2} (see “Methods” section) as local probes of the low-frequency spin dynamics of the DW.

### NRDNMR relaxation time measurements

We now discuss the *T*_{1} and *T*_{2} results of the *ν* = 2 IQHF for \({R}_{14,23}^{+B}\). Figure 2a shows that *T*_{1} near *ν* = 2 is very short (~8 s) and gradually increases until \(\left|2-\nu \right| \sim 0.023\) (dashed line). Then *T*_{1} 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 *T*_{1} are presented in Fig. 2b in terms of the nuclear spin-lattice relaxation rate 1/*T*_{1}. It is clear that the dependence of 1/*T*_{1} 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/*T*_{1} at *T* → 0 in phase I is indicative of strong magnetic fluctuations coupling to the nuclei. Note that *T*_{1} 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 2DEG^{28}. 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: *T*_{1} *T* = const^{28}. 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/*T*_{1} 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/*T*_{1} 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/*T*_{1} is proportional to the resistance of the QH state^{29}. 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 textures^{31}. 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 skyrmion^{6}. 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 interpretation^{28}. Below we show that these excitations are necessarily related to the interpretation and discussion of our data.

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 excitation^{6}. 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 −88^{32} makes the ratio of the Zeeman and Coulomb energies 100 times larger than that required for its formation^{33}. 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 skyrmion^{13}. 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 nanotubes^{34,35}, where a gapless *XY* spin wave mode with respect to the broken SO(2) symmetry would relax the nuclear spins more efficiently^{36,37} as shown by a short *T*_{1} 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 *T*_{1} becomes longer. When moving further into the region of \(\left|2-\nu \right| \; > \; 0.023\) (phase II), a sudden change in *T*_{1} 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/*T*_{1} 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 excitations^{6,39}. The number of neutral DWs increases with increasing \(\left|b\right|\), resulting in a reduced \({\triangle }_{{{{{{\rm{sw}}}}}}}\) by the SW–SW interaction^{40} (Supplementary Fig. 11) that accounts for a drop in *T*_{1} 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 *T*_{1} between the two phases at low *T* (Fig. 2a). As *T* increases, the jump is smeared out (Fig. 2a) due to a decrease in *T*_{1} in phase II caused by the Arrhenius behavior of thermal activation. Figure 2d shows the *T* and \(\nu\) dependence of *T*_{2} that is known to be dominated by dynamic fluctuating fields experienced by the nuclei in the quantum well associated with the dynamics of the 2DEG^{41}. It is shown that *T*_{2} 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 effect^{42}. 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 lattice^{41}. In addition, *T*_{2} 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/*T*_{1} 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 theory^{17}, we have made a calculation of \({b}^{\ast }=0\) in the *B*-*B*_{perp} 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.

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/*T*_{1} and *T*_{2} 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/*T*_{1} and *T*_{2} with a peak around \(\nu =2\). This nonmonotonic behavior clearly rules out the possible role of exchange interactions at the Ising transition in 1/*T*_{1} and *T*_{2} because the exchange-correlation energy varies monotonously along the \({b}^{\ast }=0\) line^{20}. 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/*T*_{1} and *T*_{2} 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 glass^{43}, 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/*T*_{1}^{44}. It is also seen that *T*_{2} 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 *T*_{1} and *T*_{2} 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 measurement^{28}, suggesting that the presence of 1D Wigner glass of the DW skyrmion is responsible for the NSR therein. Furthermore, the *T*_{1} 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 lattice^{45,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 strength^{13}. The formation of domain structure (several hundreds of nm in domain size^{15}) 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 mobility^{30}. 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 effect^{47}. It is also impossible for the skyrmions in 1D nanowires^{48} and for the magnetic DW skyrmions^{11,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 *R*_{23,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.

## Methods

### 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 *n*_{s} ~ 2.76 × 10^{15} m^{−2} and mobility of μ ~ 19.3 m^{2} 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 *R*_{SD,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, *R*_{14,23}) on a time scale of hundreds of seconds. After *R*_{14,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 *R*_{14,23} representing the depolarization of nuclei reaches a maximum Δ*R*_{nl} 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), *R*_{14,23} is decreased exponentially due to the repolarization of nuclei. The fall time of repolarization gives \({T}_{1}\) and thus *T*_{2} is obtained from *T*_{r}.

## Data availability

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

## References

- 1.
Skyrme, T. H. R. A unified field theory of mesons and baryons.

*Nucl. Phys.***31**, 556–569 (1962). - 2.
Nagaosa, N. & Tokura, Y. Topological properties and dynamics of magnetic skyrmions.

*Nat. Nanotechnol.***8**, 899–911 (2013). - 3.
Fert, A., Reyren, N. & Cros, V. Magnetic skyrmions: advances in physics and potential applications.

*Nat. Rev. Mater.***2**, 17031 (2017). - 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). - 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). - 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). - 7.
Nitta, M. Matryoshka Skyrmions.

*Nucl. Phys. B***872**, 62–71 (2013). - 8.
Jennings, P. & Sutcliffe, P. The dynamics of domain wall Skyrmions.

*J. Phys. A Math. Theor.***46**, 465401 (2013). - 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). - 10.
Cheng, R. et al. Magnetic domain wall skyrmions.

*Phys. Rev. B***99**, 184412 (2019). - 11.
Li, M. et al. Magnetic domain wall structures in Pt/Co/Ni/Ir multi-layers. Preprint at https://arXiv.org/abs/2004.07888 (2020).

- 12.
Nagase, T. et al. Observation of domain wall bimerons in chiral magnets.

*Nat. Commun.***12**, 3490 (2021). - 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). - 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). - 15.
Yang, K. F. et al. Role of chiral quantum Hall edge states in nuclear spin polarization.

*Nat. Commun.***8**, 15084 (2017). - 16.
Hirayama, Y. et al. Electron-spin/nuclear-spin interactions and NMR in semiconductors.

*Semicond. Sci. Technol.***24**, 023001 (2009). - 17.
Jungwirth, T. & MacDonald, A. H. Pseudospin anisotropy classification of quantum Hall ferromagnets.

*Phys. Rev. B***63**, 035305 (2000). - 18.
McEuen, P. L. et al. New resistivity for high-mobility quantum Hall conductors.

*Phys. Rev. Lett.***64**, 2062–2065 (1990). - 19.
Abanin, D. A. et al. Giant nonlocality near the Dirac point in graphene.

*Science***332**, 328–330 (2011). - 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). - 21.
Jungwirth, T. & MacDonald, A. H. Resistance spikes and domain wall loops in Ising quantum Hall ferromagnets.

*Phys. Rev. Lett.***87**, 216801 (2001). - 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). - 23.
Büttiker, M. Four-terminal phase-coherent conductance.

*Phys. Rev. Lett.***57**, 1761–1764 (1986). - 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). - 25.
Main, P. C. et al. Resistance fluctuations in the quantum Hall regime.

*Phys. Rev. B***50**, 4450–4455 (1994). - 26.
Trugman, S. A. Localization, percolation, and the quantum Hall effect.

*Phys. Rev. B***27**, 7539–7546 (1983). - 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). - 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). - 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). - 30.
Gervais, G. et al. Evidence for skyrmions crystallization from NMR relaxation experiments.

*Phys. Rev. Lett.***94**, 196803 (2005). - 31.
Zhu, H. et al. Pinning-mode resonance of a Skyrme crystal near Landau-level filling factor

*ν*= 1.*Phys. Rev. Lett.***104**, 226801 (2010). - 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). - 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). - 34.
Schulz, H. J. Wigner crystal in one dimension.

*Phys. Rev. Lett.***71**, 1864–1867 (1993). - 35.
Shapir, I. et al. Imaging the electronics Wigner crystal in one dimension.

*Science***364**, 870–875 (2019). - 36.
Côté, R. et al. Collective exciations, NMR, and phase transitions in Skyrme crystals.

*Phys. Rev. Lett.***78**, 4825–4828 (1997). - 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). - 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). - 39.
Mitra, A. & Girvin, S. M. Electron/nuclear spin domain walls in quantum Hall systems.

*Phys. Rev. B***67**, 245311 (2003). - 40.
Kasner, M., Palacios, J. J. & MacDonald, A. H. Quasiparticle properties of quantum Hall ferromagnets.

*Phys. Rev. B***62**, 2640–2658 (2000). - 41.
Khandelwal, P. et al. Spectroscopic evidence for the localization of skyrmions near

*ν*= 1 as*T*→0.*Phys. Rev. Lett.***86**, 5353–5356 (2001). - 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). - 43.
Akhanjee, S. & Rudnick, J. Disorder induced transition into a one-dimensional Wigner glass.

*Phys. Rev. Lett.***99**, 236403 (2007). - 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). - 45.
Read, N. & Sachdev, S. Continuum quantum ferromagnets at finite temperature and the quantum Hall effect.

*Phys. Rev. Lett.***75**, 3509–3512 (1995). - 46.
Green, A. G. Quantum-critical dynamics of the Skyrmion lattice.

*Phys. Rev. B***61**, R16299–R16302 (2000). - 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). - 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).

## Acknowledgements

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

### Affiliations

### Contributions

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

## 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 http://creativecommons.org/licenses/by/4.0/.

## About this article

### Cite this article

Yang, K., Nagase, K., Hirayama, Y. *et al.* Wigner solids of domain wall skyrmions.
*Nat Commun* **12, **6006 (2021). https://doi.org/10.1038/s41467-021-26306-8

Received:

Accepted:

Published:

DOI: https://doi.org/10.1038/s41467-021-26306-8

## Comments

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.