---
title: Non-Hermitian Skyrmions in Complex Systems
url: https://www.emergentmind.com/topics/non-hermitian-skyrmion
type: topic
---

# Non-Hermitian Skyrmions in Complex Systems

Searching arXiv for recent papers on non-Hermitian skyrmions and closely related formulations.
In recent arXiv literature, the term **non-Hermitian skyrmion** denotes several distinct constructions in which a skyrmionic map or skyrmion number is defined for a non-Hermitian generator, open-wave system, driven–dissipative field, or complexified Skyrme theory. The phrase is used for a real-space magnetic skyrmion with balanced gain and loss or a PT-symmetric anisotropy, for momentum-space pseudo-spin textures in leaky photonic bands, for dynamical skyrmions on generalized-Brillouin-zone–time manifolds, for scalar-field textures in optical microcavities, and for complex BPS Skyrmion configurations with real energies. Across these settings, the common problem is not merely whether a topological charge exists, but which state space it is built from, whether it remains quantized, and how non-Hermitian singularities, biorthogonality, damping, or antilinear symmetry modify its protection [2606.27810], [2211.09884], [2101.07963], [2606.02265], [2605.07435], [2102.05781].

## 1. Real-space magnetic formulation

The most direct use of the term concerns a magnetic \(CP^1\) texture \(z\) with balanced gain/loss along the easy axis, introduced by the anti-Hermitian term \(i\gamma \sigma_z\). The canonical PT-symmetric deformation is
$$
i\,\partial_t z = (H_0 + i\gamma\,\sigma_z)\,z ,
$$
and, because this evolution does not preserve the norm of \(z\), the physical degree of freedom is the normalized spinor \(\hat z=z/|z|\). Writing the projector \(P=\hat z\hat z^\dagger\) and \(\Gamma=\gamma\sigma_z\), the trace-preserving projective evolution is
$$
\dot P = \{\Gamma,P\} - 2\,\mathrm{Tr}(\Gamma P)\,P .
$$
To probe the left–right structure spatially, the local non-Hermitian generator is taken pointwise as
$$
H(\mathbf{r}) = \mathbf{d}(\mathbf{r})\cdot\boldsymbol{\sigma} + i\gamma\,\sigma_z ,
$$
where \(\mathbf d(\mathbf r)\) encodes the Hermitian texture field and \(i\gamma\sigma_z\) implements PT-symmetric balanced gain/loss or anisotropy [2606.27810].

In the Hermitian limit \(\gamma=0\), the magnetization is a unit Bloch vector \(\mathbf m(\mathbf r)\in S^2\). A standard axisymmetric skyrmion is parameterized by
$$
\mathbf{m}(r,\phi) = \bigl(\sin\theta(r)\cos\phi,\;\sin\theta(r)\sin\phi,\;\cos\theta(r)\bigr),
$$
with \(\theta(0)=\pi\) and \(\theta(\infty)=0\). The equator ring is the locus \(r=r_0\) where \(\theta(r_0)=\pi/2\), equivalently \(m_z=0\). Its Hermitian topological charge is the Pontryagin index
$$
Q = \frac{1}{4\pi}\int d^2 r\;\mathbf{m}\cdot\bigl(\partial_x \mathbf{m}\times \partial_y \mathbf{m}\bigr)\in\mathbb{Z}.
$$
The non-Hermitian construction begins from this standard skyrmion but alters the state space in which the charge is defined.

## 2. Charge splitting, biorthogonality, and exceptional-point breakdown

A central result of the magnetic non-Hermitian skyrmion construction is that the Hermitian skyrmion charge splits into two inequivalent quantities. The first is the right-state Bloch vector \(n_a(\mathbf r)=\hat z^\dagger \sigma_a \hat z\in S^2\), with right-state charge
$$
Q_R = \frac{1}{4\pi}\int n\cdot(\partial_x n\times\partial_y n)\,d^2r\in\mathbb{Z}.
$$
The second is the biorthogonal Bloch vector
$$
m_a(\mathbf{r}) = \frac{\langle z_L|\sigma_a|z_R\rangle}{\langle z_L|z_R\rangle},
$$
which is generally complex and satisfies \(m\cdot m=1\), so that it lives on the complex quadric \(\{m\in\mathbb{C}^3\mid m\cdot m=1\}\cong T^*S^2\). The corresponding biorthogonal charge is
$$
Q_B = \frac{1}{4\pi}\int m\cdot(\partial_x m\times\partial_y m)\,d^2r .
$$
These two charges coincide in the Hermitian limit and part ways for \(\gamma\neq 0\): \(Q_R\) stays integer, while \(Q_B\) is complex and nonquantized once gain/loss is turned on [2606.27810].

The right-state sector remains homotopy protected because the PT flow reduces exactly to a Gilbert-type relaxation on the target sphere. In the paper’s formulation, the PT deformation viewed through the right state is indistinguishable from a relaxation of the texture toward the \(z\) axis, and the right Bloch vector remains a smooth map into \(S^2\). Under smooth evolution, \(Q_R\) therefore cannot change. The only stated way it could change is a true singularity such as \(z\to 0\), where the projective class becomes undefined.

The biorthogonal sector behaves differently because its target manifold is noncompact and because the local generator can reach an exceptional point. The squared eigenvalues of the local generator are
$$
\varepsilon^2 = \|\mathbf{d}\|^2 - \gamma^2 + 2i\gamma\,d_z ,
$$
and the exceptional-point condition is
$$
d_z = 0 \quad\text{and}\quad d_x^2 + d_y^2 = \gamma^2 .
$$
On a skyrmion, \(d_z=0\) occurs on the equator ring. For a hedgehog of constant modulus \(\|\mathbf d\|=h\), the exceptional-point criterion becomes \(\gamma=h\). At that ring the left and right eigenvectors coalesce, the phase rigidity
$$
r = \frac{|\langle z_L|z_R\rangle|}{\sqrt{\langle z_L|z_L\rangle\,\langle z_R|z_R\rangle}}
$$
tends to zero, and the biorthogonal Bloch field diverges. In the equatorial two-level caricature,
$$
H = h\,\sigma_x + i\gamma\,\sigma_z
$$
yields
$$
m = \left(\frac{-h}{\sqrt{h^2-\gamma^2}},\; 0,\; \frac{-i\gamma}{\sqrt{h^2-\gamma^2}}\right),
$$
so \(m\) diverges as \(\gamma\to h\). The paper identifies this loss of quantization as the real-space topological counterpart of the analyticity breakdown that a causal response function suffers at an exceptional point. In this sense, topological protection is no longer a single statement in the non-Hermitian magnetic problem: it splits into a right-state homotopy invariant and a biorthogonal quantity that becomes complex, nonquantized, and singular at the exceptional point.

## 3. Momentum-space and dynamical skyrmions

A second usage of the term appears in open photonic band structures. In a 1D optical lattice with a symmetry-protected bound state in the continuum at \(\Gamma\), the familiar polarization-vortex winding number of the dark branch is \(W=\mathrm{sgn}(U)\) and flips when a symmetry-preserving band inversion changes the sign of \(U\). To obtain a conserved invariant, the band identity is incorporated into a momentum-space pseudo-spin
$$
\epsilon(\mathbf q)=(E_x,E_y)/E_0,\qquad M(\mathbf q)=(\epsilon_x,\epsilon_y,U),\qquad n(\mathbf q)=M(\mathbf q)/\|M(\mathbf q)\| .
$$
The non-Hermitian skyrmion number is then
$$
Q_{\mathrm{sk}}=\frac{1}{4\pi}\int n(\mathbf q)\cdot\bigl[\partial_{q_x}n(\mathbf q)\times\partial_{q_y}n(\mathbf q)\bigr]\,dq_x\,dq_y
= \frac{pW}{2}=\frac{1}{2},
$$
with \(p=\mathrm{sgn}(U)\). The resulting texture is a half-skyrmion: the paper describes a transition from an antimeron to a meron-like texture through band inversion, while the half-charge skyrmion number remains conserved even though the BIC polarization-vortex winding does not [2211.09884].

A third usage appears in non-Bloch quench dynamics with non-Hermitian skin effects. There the relevant manifold is not real space or momentum space alone, but the generalized Brillouin zone times time. Writing \(\beta=r e^{ik}\) on the generalized Brillouin zone, the dynamical pseudospin is defined biorthogonally with a metric operator \(\eta(k)\), and one constructs a real unit vector \(n(k,t)\in S^2\). The dynamic skyrmion charge is
$$
Q_{\mathrm{sk}}=\frac{1}{4\pi}\int d\kappa\int dt\; n(\kappa,t)\cdot\bigl[\partial_\kappa n(\kappa,t)\times\partial_t n(\kappa,t)\bigr],
$$
or, on submanifolds between fixed points \(k_m\) and \(k_n\),
$$
C_{mn}=\frac{1}{4\pi}\int_{k_m}^{k_n}dk\int_{0}^{T_k}dt\; [n(k,t)\times\partial_t n(k,t)]\cdot \partial_k n(k,t).
$$
When the generalized Brillouin zones of the pre- and post-quench Hamiltonians coincide and both spectra are real, fixed points occur at \(k=0,\pm\pi\). If the non-Bloch winding numbers differ, \(\nu_i\neq \nu_f\), the fixed points are of different types and dynamic skyrmions with \(C_{mn}=\pm 1\) appear; when \(\nu_i=\nu_f\), the dynamic Chern numbers vanish. Even when the generalized Brillouin zones do not coincide, global signatures persist through the time-averaged polarization \(\bar n_3\), so the skyrmion construction remains a detection scheme for non-Bloch topology in the presence of skin effects [2101.07963].

## 4. Scalar-field skyrmions in driven–dissipative microcavities

In optical microcavities the term is extended to a scalar complex field rather than a spin texture. The coherent exciton–polariton condensate is described by a single complex order parameter \(\psi(\mathbf r,t)\), and the skyrmionic texture is built from amplitude and phase gradients of that scalar field. For a stationary state \(\psi_s(\mathbf r)\), the construction introduces
$$
R=(\partial_x\psi_s,\partial_y\psi_s,\psi_s)^T,\qquad
\mathcal R=\mathrm{Re}\{e^{i\omega t}R\}=(\mathcal X,\mathcal Y,\mathcal Z)^T,
$$
and then the unit vector
$$
n(\mathbf r)=\mathcal R(\mathbf r)/|\mathcal R(\mathbf r)|.
$$
Over a single Skyrmion cell, the charge is
$$
Q=\frac{1}{4\pi}\iint n\cdot(\partial_x n\times \partial_y n)\,dx\,dy,
$$
with the cell boundary chosen by critical lines of \(\mathcal R\). In this formulation the Skyrmion charge is a bulk area integral mapping \(\mathbb R^2\to S^2\) and is explicitly distinct from vortex winding \(m=(1/2\pi)\oint \nabla\phi\cdot d\ell\) [2606.02265].

The driven–dissipative dynamics is governed by the open-dissipative Gross–Pitaevskii equation with an incoherent reservoir,
$$
i\hbar \partial_t \psi = \Big[ -\frac{\hbar^2}{2m} \nabla^2 + g_c |\psi|^2 + g_r n + \frac{i\hbar}{2}(R_c n - \gamma_c) \Big] \psi + E,
$$
$$
\partial_t n = P - (\gamma_r + R_c |\psi|^2) n.
$$
For nonresonant pumping, the continuity equation
$$
\partial_t I + \nabla\cdot j = (R_c n - \gamma_c) I
$$
shows that spatially inhomogeneous gain/loss requires a compensating current divergence in steady state. Near threshold, a Gaussian pump reduces the problem to a non-Hermitian Schrödinger equation with a complex harmonic frequency \(\omega_c\), so the threshold mode necessarily carries nonzero phase curvature through \(\mathrm{Im}\,\omega_c\neq 0\). That phase curvature produces radial sign reversals in \(\mathcal Z\) and nodal rings in \(\mathcal R_2=(\mathcal X,\mathcal Y)\), yielding an isolated Skyrmion of charge \(Q=1\) inside the first Skyrmion cell. The paper emphasizes that nonlinear interactions are not required for existence: Skyrmions persist even with \(g_c=g_r=0\), and the essential ingredients are localized gain and finite lifetime.

This scalar non-Hermitian skyrmion framework supports isolated Skyrmions, self-organized Skyrmion lattices, and Skyrmion bags in moiré lattices. Resonant interference can generate hexagonal lattices with \(Q=1\) per unit cell, whereas structured incoherent pumping can induce spontaneous phase synchronization and nonresonant Skyrmion bags on a \(1500\times 1500\,\mu\mathrm m^2\) domain. The same platform provides nonlinear all-optical control: an additional nonresonant ring pump changes the propagation phase through the reservoir-induced blueshift and switches the central Skyrmion number between distinct integers. The derivative field \(R_2\) also determines two force densities, \(F_\phi\propto I\nabla\phi=\mathrm{Im}\{\psi_s^*R_2\}\) and \(F_I\propto \nabla I=2\,\mathrm{Re}\{\psi_s^*R_2\}\), so isolated scalar Skyrmions generate inward intensity-gradient forces and outward phase-gradient forces.

## 5. Damped magnons and non-Hermitian topology on skyrmion strings

A different but related direction studies a skyrmion texture as the background that generates non-Hermitian topology for its excitations. In a three-dimensional chiral magnet hosting a skyrmion-string lattice, the static texture is a periodic array of skyrmion tubes extended along \(z\). The spin dynamics is described by a generalized Landau–Lifshitz–Gilbert equation with a real symmetric positive-semidefinite damping kernel \(\alpha_{r,r'}\), allowing both local Gilbert damping and nonlocal damping between different sites. Linearization about the stationary texture produces a bosonic Bogoliubov–de Gennes problem for the Nambu spinor \((\psi,\psi^*)^T\),
$$
i\hbar \partial_t
\begin{pmatrix}
\psi\\ \psi^*
\end{pmatrix}
=
H_\tau(\alpha)
\begin{pmatrix}
\psi\\ \psi^*
\end{pmatrix},
\qquad
H_\tau(\alpha)=(1+iA_\tau)^{-1}H_\tau^{(0)},
$$
with complex magnon frequencies under damping [2605.07435].

To leading order in small damping,
$$
E_{nk}(\alpha)=E_{nk}^{(0)}\,[1-i\,\alpha_{\mathrm{eff},n}(k)] + O(\alpha^2),\qquad
\alpha_{\mathrm{eff},n}(k)=\alpha_k\,\eta_{nk},
$$
where \(\eta_{nk}\ge 1\) quantifies particle–hole mixing and elliptical precession. The point-gap winding number is
$$
W_n(\omega_0)=\frac{1}{2\pi i}\int_{BZ} dk\; \partial_k \ln[\omega_n(k)-\omega_0].
$$
A nontrivial \(k\)-dependence of \(\alpha_{\mathrm{eff},n}(k)\) produces loops in the complex spectrum and enables the non-Hermitian skin effect. The paper stresses that this can occur even without nonlocal damping: with purely local damping \(\alpha_k\equiv \alpha_0\), the texture-induced \(k\)-dependence of \(\eta_{nk}\) alone can give \(W_n\neq 0\). With nearest-neighbor nonlocal damping along one direction,
$$
\alpha_k=\alpha_0+2\alpha_\parallel \cos(kc),
$$
and weak variation of \(\eta_{nk}\), the sign rule becomes
$$
W_n(E)=\mathrm{sgn}(k_0)
$$
for a stable particle band with minimum at \(k_0\neq 0,\pm \pi/c\). Under open boundary conditions, \(W=+1\) and \(W=-1\) correspond to top-boundary and bottom-boundary skin localization, respectively.

In the skyrmion-string lattice used in the demonstrations, the counterclockwise gyration mode shows strong \(k\)-dependence of \(\eta_{nk}\) and, with local damping only, yields a single prominent spectral loop with \(W_3(E)=+1\). With nonlocal damping, the localization side becomes band selective: at \(B/b=0.7\), low-energy bands with \(k_0\ge 0\) localize at the top boundary, whereas at \(B/b=0.2\) the seventh band has \(k_0<0\) and localizes at the bottom boundary. The resulting directional propagation after a magnetic-field pulse is a dynamical signature of point-gap topology. This framework sharply separates two notions of topology: the skyrmion number of the static texture and the spectral winding number of damped magnons. The paper also states explicitly that exceptional points are not required for this point-gap topology or for the non-Hermitian skin effect in the perturbative regime considered.

## 6. Complex BPS Skyrmions and conceptual distinctions

In Skyrme-type field theory, non-Hermitian skyrmions appear as complex-valued \(SU(2)\)-based or complexified field configurations whose energies remain real because of an antilinear CPT symmetry. The general setting uses
$$
U(x)=\exp\!\big[i\,\xi(x)\,\hat n(x)\cdot \sigma\big],
$$
with Skyrme-type Lagrangian
$$
\mathcal L=\mathcal L_0+\mathcal L_2+\mathcal L_4+\mathcal L_6+\mathcal L_{\mathrm{pot}},
$$
topological charge
$$
B=-\frac{1}{24\pi^2}\int d^3x\,\epsilon^{ijk}\,\mathrm{Tr}(L_iL_jL_k),
$$
and static energy \(E[U]=\int d^3x\,\mathcal H\). In the BPS sector, the Hamiltonian density can be written as \(\mathcal H=A^2+\tilde A^2\), so the Bogomolny equation \(A=\pm \tilde A\) reduces the problem to first-order radial equations. Non-Hermiticity is introduced by complex field shifts, complex potentials, or complex couplings, while energy reality follows when a modified antilinear CPT symmetry maps the Hamiltonian functional to its parity time-reversed complex conjugate and maps a solution either to itself or to a degenerate-energy partner [2102.05781].

The paper develops several concrete families. In the complex-boosted BPS model, \(\xi\mapsto \xi-i\kappa\) with \(\kappa=\mathrm{arctanh}\,\varepsilon\) yields a non-Hermitian Hamiltonian density that remains pseudo-Hermitian and is related by a Dyson map to a Hermitian counterpart. This sector contains complex compactons, step solutions, cusp solutions, shell solutions, and purely imaginary compactons, all with real energies. Other constructions include semi-kink solutions, massless shell solutions with \(E=0\), Bender–Boettcher-type potentials \(V_{\mathrm{BB}}(\xi)=(i\xi)^\varepsilon \sin^4\xi\) with real-energy sectors on discrete sets of \(\varepsilon\), and a solvable submodel with purely imaginary coupling \(\lambda=i\Lambda\) and real energy
$$
E^{(2)}=\pm 2\pi n\Lambda\,[\sin\bar c-\bar c].
$$
In this field-theoretic usage, the skyrmion is non-Hermitian because the fields, couplings, or potentials are complex, not because of gain/loss or an exceptional spectrum.

Taken together, these formulations show that **non-Hermitian skyrmion** is a genuinely plural concept. It does not always mean a real-space magnetic spin texture, and its “charge” is not universally the integer Pontryagin index of a Hermitian map. Depending on the setting, the relevant invariant may remain integer and homotopy protected, become complex and nonquantized, stay fixed as a half-charge meron across band inversion, appear as a dynamic Chern number on generalized-Brillouin-zone–time space, or arise from a scalar derivative field rather than a spin field. Likewise, exceptional points are central to the breakdown of the biorthogonal magnetic charge, but they are not the organizing mechanism in every non-Hermitian skyrmion construction. Current literature therefore treats the subject less as a single topological object than as a family of skyrmionic structures whose definition is inseparable from the non-Hermitian state space in which they are embedded [2606.27810], [2211.09884], [2101.07963], [2606.02265], [2605.07435].

Source: https://www.emergentmind.com/topics/non-hermitian-skyrmion