Pith. sign in

REVIEW 3 major objections 5 minor 26 references

Quantum computing of magnetic-skyrmion-like patterns in Heisenberg ferromagnets

T0 review · 3 major / 5 minor · reviewed 2026-08-07 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read A variational quantum eigensolver on a noiseless simulator maps the field-driven transition into skyrmion-like order in a two-dimensional Heisenberg ferromagnet and beats classical diagonalization beyond 17 sites.

desk verdict VQE study of DMI Heisenberg model with a plausible but unverified transition claim; missing exact-diagonalization check at the transition is the key gap. read the letter →

arxiv 2505.19808 v1 pith:NOENOKHN submitted 2025-05-26 quant-ph

Matej Komelj , Vinko Sršan , Kristina Žužek , Sašo Šturm This is my paper · ORCID
classification quant-ph PACS 74.20.pq74.25.Kc71.15.Mb
keywords variationalquantumeigensolverDzyaloshinskii-MoriyainteractionHeisenbergferromagnetskyrmionstopologicalchargehardware-efficientansatzzero-temperaturephasetransitionspin-1/2lattice
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

The paper sets out to show that the ground state of a two-dimensional spin-1/2 Heisenberg ferromagnet with Dzyaloshinskii-Moriya interaction (DMI) can be found with a variational quantum eigensolver (VQE) running on a noiseless simulator, and that this quantum route becomes faster than classical Lanczos diagonalization once the lattice exceeds 17 sites. The computed energy, magnetization, and topological charge all jump abruptly as the perpendicular magnetic field $B_z$ crosses a critical value, which the authors read as a zero-temperature transition into a vortex-like, skyrmion-like magnetic texture. The transition is controlled by the competition between exchange coupling and DMI, and it produces Bloch-type or Néel-type patterns depending on the direction of the DMI vector. If the claim holds, quantum skyrmion-like states are stable, field-switchable objects with a measurable magnetization jump, making them plausible building blocks for spintronics or information storage.

What carries the argument

The central object is the Hamiltonian $$H=J_\parallel\sum_{\langle i,j\rangle}(\$\sigma$^x_i\$\sigma$^x_j+\$\sigma$^y_i\$\sigma$^y_j)+J_\perp\sum_{\langle i,j\rangle}\$\sigma$^z_i\$\sigma$^z_j+\sum_{\langle i,j\rangle}\vec D_{ij}\cdot(\vec\sigma_i\times\vec\sigma_j)+B_z\sum_i\$\sigma$^z_i,$$ expressed in Pauli matrices and mapped directly onto qubits. The VQE approximates the ground state with a hardware-efficient ansatz, a shallow circuit made of two layers of single-qubit Euler rotations surrounding one entangler, with $6N$ variational parameters optimized against the energy. The topological charge is obtained from the Pauli expectation values at the lattice nodes by triangulating the lattice and summing the per-triangle solid-angle contributions, giving a discrete version of the winding number. The mechanism that carries the argument is the competition among ferromagnetic exchange, DMI, and the Zeeman field: DMI favors neighboring moments perpendicular to each other, the field favors out-of-plane alignment, and the transition appears as a discontinuity in the observables.

What would settle it

Compute the exact ground state for the same 16-site lattice at fields just below and above $B_z = 1.884$ (and $B_z = 1.222$) and check whether the energy, $m_x$, and topological charge jump. If the exact observables vary smoothly across those fields, the discontinuity seen by VQE is a variational failure rather than a physical transition.

Watch

Extended reading notes

Core claim

On a 16-site square lattice, the VQE ground state of the DMI-Heisenberg Hamiltonian exhibits a sharp transition as a function of the external field: at $B_z = 1.884$ (for $\vec D_{ij}\parallel \vec R_{ij}$ with $J_\perp = 0.5|J_\parallel|$) and at $B_z = 1.222$ (for $\vec D_{ij}\perp \vec R_{ij}$ with $J_\perp = 0.25|J_\parallel|$), the energy, the in-plane magnetization component $m_x$, and the topological charge $Q$ all jump. Below the transition the magnetization has a preferred in-plane direction; above it, the pattern becomes a rotating vortex resembling a Bloch-type skyrmion (helicity $\gamma=\pi/2$) or a Néel-type skyrmion (helicity $\gamma=0$). The topological charges are not integer because the 16-site lattice is too small to host a complete skyrmion, but a $5\times 5$ lattice already shows a Bloch-type skyrmion-like pattern, supporting the interpretation that these are precursors of genuine quantum skyrmions at $T=0$. The paper also reports that VQE scales roughly as $O(N^{1.1})$ while Lanczos scales as $O(N^{2.1})$, making VQE faster for $N>17$.

Load-bearing premise

The result stands or falls on whether the optimized variational circuit truly reaches the ground state at every field value, especially across the transition, since the paper checks only against a standard classical solver and does not report overlaps or convergence data at the transition points.

Editorial extensions

If this is right

  • For lattices with more than 17 sites, VQE on a noiseless simulator finds the ground state faster than serial Lanczos diagonalization, with measured scalings of about $O(N^{1.1})$ versus $O(N^{2.1})$.
  • The ground state of a DMI-Heisenberg ferromagnet in a perpendicular field undergoes a sharp zero-temperature transition into a skyrmion-like texture, so skyrmion-like order can exist without thermal fluctuations.
  • The transition is accompanied by a jump in magnetization large enough to be measured, which would allow experimental detection and possible use as a stable information carrier.
  • A $5\times 5$ lattice calculation shows a Bloch-type skyrmion-like pattern, indicating that the small-lattice non-integer topological charges are finite-size effects and that larger systems should host more complete skyrmionic textures.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • A decisive test of the ansatz's expressiveness would be to repeat the 16-site calculation with a different variational circuit or with exact diagonalization and check whether the discontinuity in energy and in-plane magnetization persists; this would separate a physical transition from a variational collapse.
  • If the transition is genuine, the same VQE pipeline could map the full phase diagram over exchange anisotropy, DMI strength, and external field for lattice sizes beyond classical reach, where Lanczos becomes impractical.
  • Because the simulator is noiseless, the crossover at 17 sites likely shifts on real quantum hardware once gate noise, measurement overhead, and optimization costs are included; the efficiency claim should be read as algorithm-scaling evidence rather than a hardware benchmark.
  • The non-integer topological charges suggest a finite-size scaling study toward larger lattices could reveal whether the charge approaches the integer skyrmion value in the thermodynamic limit, and whether the transition sharpens into a true phase boundary.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 5 minor

Summary. The paper applies the variational quantum eigensolver (VQE) on a noiseless simulator to a spin-1/2 anisotropic Heisenberg (XXZ) model with Dzyaloshinskii-Moriya interaction on small open-boundary square lattices. It reports runtime scaling fits for VQE versus the Lanczos method and claims that VQE becomes faster for more than 17 sites. For a 4x4 lattice, the paper finds discontinuities in the energy, in the magnetization components, and in the topological charge as a function of the perpendicular magnetic field, at B_z = 1.884 for D parallel to R and at B_z = 1.222 for D perpendicular to R. These jumps are interpreted as transitions to Bloch- and Neel-type skyrmion-like states, while the authors acknowledge that the computed topological charges are not integer and hence the patterns are not proper skyrmions. The final claim is that VQE is a promising tool for quantum magnetism and that the results call for experiments on skyrmion-like information carriers.

Significance. If the observed transitions are genuine ground-state features, the paper would support the existence of zero-temperature quantum skyrmion-like phases and demonstrate a useful VQE application in quantum magnetism. The study is commendably concrete in benchmarking VQE against Lanczos on the same Hamiltonian, and no parameter is fitted to the target result; the field-driven transition is read directly from computed observables. The open-source Tangelo implementation also aids reproducibility. However, the significance is conditional on variational convergence and on the representativeness of the 4x4 open-boundary lattice. The absence of exact comparisons at the transition fields and the acknowledged non-integer topological charges leave the central physics claim in need of stronger support before it can be regarded as established.

major comments (3)
  1. [Section 3, Figs. 3 and 5] The central claim of field-driven transitions rests entirely on VQE results at B_z = 1.884 and B_z = 1.222, but no comparison with exact Lanczos ground states is reported at or near these transition fields. Since the text states that for N = 16 the classical approach is still slightly faster, exact diagonalization is completely feasible. Reporting the energy difference between the VQE state and the exact ground state, or the squared overlap, at every field point would distinguish a genuine level crossing or avoided crossing from a variational collapse between two local minima of the hardware-efficient ansatz. Without such data, the discontinuities in Figs. 3 and 5 could be artifacts of the optimization failing to reach the true ground state on one side of the transition.
  2. [Section 2, Fig. 2] The efficiency claim that VQE is faster than Lanczos for more than 17 sites is based on scaling fits O(N^1.1) versus O(N^2.1) obtained from only five lattice sizes (7, 9, 16, 19, 25), with no error bars, no convergence tolerances for either method, no optimizer details, and no definition of the 'arbitrary units' for runtime. The inset in Fig. 2 places the crossover at N > 17, but the only data points beyond N = 16 are the VQE point at N = 25; there is no Lanczos point at N = 25 to support the fitted crossover. The abstract's statement that VQE 'turns out to be a more efficient approach' should be qualified accordingly, or the scaling analysis should be strengthened with more sizes and with error estimates.
  3. [Section 4 and Fig. 7] The inference from the 4x4 calculations to skyrmion-like phases is weakened by the authors' own admission that the patterns in Figs. 4(b) and 6(b) are not proper skyrmions because their topological charges are non-integer. The supporting 5x5 example in Fig. 7 is a single magnetization pattern at B_z = 1.5 and does not demonstrate a transition or a field-driven evolution on that lattice size. A systematic finite-size study, such as a field sweep on a 5x5 or 6x6 lattice with corresponding Lanczos checks for the smaller sizes, would be needed to ascertain whether the discontinuities persist and whether the topological charge approaches an integer value. In the absence of such data, the extrapolation to stable skyrmionic phases is qualitative.
minor comments (5)
  1. [Eq. (2) and following text] The text says the transition is confirmed by 'the third term in Eq. (2)', but the Zeeman term is the fourth term in the displayed Hamiltonian; please correct the reference.
  2. [Ansatz definition, Section 2] The sentence defining U(θ) reads 'expressed by the R_x(theta) and R_x(theta) gates' but then gives U = R_z(theta1) R_x(theta2) R_z(theta3); the first gate should presumably be R_z, not R_x.
  3. [Fig. 5 caption] The caption contains the stray word 'colorblack' between 'and' and 'J_perp', which should be removed.
  4. [Section 2 text] There is a typo 'Zeemnan interaction' and an awkward 'Hence. a further optimization' with a period instead of a comma; both should be fixed.
  5. [Eq. (1) and Section 2] The topological charge in Eq. (1) is defined for a continuous normalized vector field, while the lattice calculation uses Pauli expectation values that are not explicitly normalized before interpolation; a sentence clarifying the normalization convention would help readers reproduce the values in Figs. 3 and 5.

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity: the transition is read off from VQE-computed observables, and the VQE results are benchmarked against Lanczos on the same Hamiltonian.

full rationale

The paper's load-bearing steps are: (i) optimize a hardware-efficient ansatz to approximate the ground state of Hamiltonian (2) for fixed model parameters; (ii) compute energy, magnetization, and topological charge as functions of B_z; and (iii) interpret a discontinuity in those observables as evidence of a transition to a skyrmion-like state. None of these steps fits a parameter to the claimed result. The ansatz is a generic HEA from the cited literature, and the exchange, DMI, and field parameters are chosen independently rather than tuned to produce a jump. The paper explicitly states that the only validation criterion is agreement with Lanczos-method results within numerical accuracy; Lanczos diagonalization is an independent external check on the same Hamiltonian, so it does not make the argument circular. The efficiency claim is likewise an empirical runtime comparison, O(N^1.1) for VQE versus O(N^2.1) for Lanczos, not an input into the physical conclusion. The admitted non-integer topological charges and the absence of shown convergence data or overlap checks are legitimate correctness and robustness concerns, not circularity, because the target discontinuity is not baked into the ansatz, the optimizer, or any fitted parameter. There is no self-citation chain, no imported uniqueness argument, and no renaming of a known result as a new derivation.

Assumptions & free parameters 4 free parameters · 5 assumptions · 0 invented entities

The central claim depends on several unverified choices: the model Hamiltonian, the ansatz expressibility and optimization convergence, the reliability of the discrete topological charge, and the representativeness of 4x4 and 5x5 open lattices. These are not fitted to the target result but are assumptions about the validity of the numerical workflow. The runtime scaling exponents are fitted to the authors' own benchmarks.

free parameters (4)
  • DMI strength / exchange ratio |D|/|J_parallel| = 1
    Set to unity for simplicity in the paragraph after Eq. (2); the DMI strength is not varied.
  • easy-plane anisotropy J_perp / |J_parallel| = 0.5 (D parallel to R) and 0.25 (D perpendicular to R)
    Chosen by hand; the text calls the choice 'quite arbitrary' and changes it between the two DMI orientations to obtain a transition.
  • VQE runtime scaling exponent = 1.1
    Fitted to measured runtime versus lattice size in FIG. 2; used to claim efficiency for N>17.
  • Lanczos runtime scaling exponent = 2.1
    Fitted to measured runtime versus lattice size in FIG. 2 for comparison.
assumptions (5)
  • domain assumption The Hamiltonian in Eq. (2) is a faithful model of zero-temperature quantum skyrmion physics.
    The paper assumes DMI and easy-plane exchange are the relevant interactions, citing Refs. 7 and 8; this is standard but not derived.
  • standard math The VQE variational principle gives the ground-state energy.
    Uses the Rayleigh-Ritz bound that an expectation value is always above the ground energy, but convergence of the optimizer is not proven.
  • ad hoc to paper The hardware-efficient ansatz is expressive enough to represent the ground state at all field values.
    Not established; the paper only compares with Lanczos where stated, and no overlap or convergence data are shown.
  • domain assumption Piecewise-linear interpolation of magnetization from 16 or 25 lattice nodes gives a meaningful topological charge Q.
    The method is a simplified winding number, but for very few nodes the continuum definition may not be reliable; non-integer Q in the results supports this concern.
  • ad hoc to paper Open-boundary finite lattices of 4x4 and 5x5 are representative of thermodynamic-limit behavior.
    The authors themselves attribute non-integer Q to restricted geometry, so finite-size representativeness is assumed for the extrapolation.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Quantum computing of magnetic-skyrmion-like patterns in Heisenberg ferromagnets." pith.science (2026). https://pith.science/paper/NOENOKHN

@misc{pith2026250519808,
  author       = {Pith},
  title        = {Pith review of: Quantum computing of magnetic-skyrmion-like patterns in Heisenberg ferromagnets},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/NOENOKHN}},
  note         = {Machine review of arXiv:2505.19808}
}
read the original abstract

We diagonalize the quantum two-dimensional spin-1/2 Heisenberg model with Dzyaloshinskii-Moriya interaction (DMI) by applying the variational quantum eigensolver, running on a quantum-computer simulator, which turns out to be a more efficient approach than a classical direct diagonalization for systems with more than 17 sites. The calculated external-magnetic-field dependence of the total energy, of the magnetization, as well as of the topological charge exhibits a distinctive discontinuity which hints for the existence of zero-temperature magnetic skyrmions-like structures at the quantum level, controlled by the combination of the exchange-coupling and the DMI parameters. The potentially measurable jump in the magnetization upon changing the field indicates the investigated objects as stable enough for eventual applications in spintronics or even as information carriers.

Figures

Figures reproduced from arXiv: 2505.19808 by the authors.

Figure 1
Figure 1. FIG. 1. The DMI vector [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. FIG. 2 [PITH_FULL_IMAGE:figures/full_fig_p002_2.png] view at source ↗
Figure 4
Figure 4. FIG. 4. The magnetization pattern for [PITH_FULL_IMAGE:figures/full_fig_p003_4.png] view at source ↗
Figures from the paper (3 more)
Figure 5
Figure 5. Figure 5: FIG. 5 [PITH_FULL_IMAGE:figures/full_fig_p003_5.png]
Figure 6
Figure 6. Figure 6: FIG. 6. The magnetization pattern for [PITH_FULL_IMAGE:figures/full_fig_p003_6.png]
Figure 7
Figure 7. Figure 7: FIG. 7 [PITH_FULL_IMAGE:figures/full_fig_p004_7.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

26 extracted references · 14 canonical work pages

  1. [1]

    40 years of quantum computing, Nature Reviews Physics 4, 1 (2022)

    work page 2022

  2. [2]

    Alexeev, D

    Y. Alexeev, D. Bacon, K. R. Brown, R. Calderbank, L. D. Carr, F. T. Chong, B. DeMarco, D. Englund, E. Farhi, B. Fefferman, A. V. Gorshkov, A. Houck, J. Kim, S. Kim- mel, M. Lange, S. Lloyd, M. D. Lukin, D. Maslov, P. Maunz, C. Monroe, J. Preskill, M. Roetteler, M. J. Savage, and J. Thompson, Quantum computer systems for scientific discovery, PRX Quantum2,...

    2021

  3. [3]

    Psaroudaki and C

    C. Psaroudaki and C. Panagopoulos, Skyrmion qubits: A new class of quantum logic elements based on nanoscale magnetization, Phys. Rev. Lett.127, 067201 (2021)

    2021

  4. [4]

    J. Xia, X. Zhang, X. Liu, Y. Zhou, and M. Ezawa, Univer- sal quantum computation based on nanoscale skyrmion helicity qubits in frustrated magnets, Phys. Rev. Lett. 130, 106701 (2023)

    work page 2023

  5. [5]

    del Ser, I

    N. del Ser, I. El Achchi, and A. Rosch, Fractional topo- logical charges in two-dimensional magnets, Phys. Rev. B110, 094442 (2024)

    work page 2024

  6. [6]

    G. Yin, Y. Li, L. Kong, R. K. Lake, C. L. Chien, and J. Zang, Topological charge analysis of ultrafast single skyrmion creation, Phys. Rev. B93, 174403 (2016)

    work page 2016

  7. [7]

    Rold´ an-Molina, M

    A. Rold´ an-Molina, M. J. Santander, A. S. Nunez, and J. Fern´ andez-Rossier, Quantum fluctuations stabilize skyrmion textures, Phys. Rev. B92, 245436 (2015)

    work page 2015

  8. [8]

    Haller, S

    A. Haller, S. Groenendijk, A. Habibi, A. Michels, and T. L. Schmidt, Quantum skyrmion lattices in heisenberg ferromagnets, Phys. Rev. Res.4, 043113 (2022)

    work page 2022

Show all 26 references
  1. [9]

    Siegl, E

    P. Siegl, E. Y. Vedmedenko, M. Stier, M. Thorwart, and T. Posske, Controlled creation of quantum skyrmions, Phys. Rev. Res.4, 023111 (2022)

    2022

  2. [10]

    Dzyaloshinsky, A thermodynamic theory of “weak” fer- romagnetism of antiferromagnetics, Journal of Physics and Chemistry of Solids4, 241 (1958)

    I. Dzyaloshinsky, A thermodynamic theory of “weak” fer- romagnetism of antiferromagnetics, Journal of Physics and Chemistry of Solids4, 241 (1958)

    1958

  3. [11]

    Moriya, Anisotropic Superexchange Interaction and Weak Ferromagnetism, Physical Review120, 91 (1960)

    T. Moriya, Anisotropic Superexchange Interaction and Weak Ferromagnetism, Physical Review120, 91 (1960)

    1960

  4. [12]

    Lohani, C

    V. Lohani, C. Hickey, J. Masell, and A. Rosch, Quantum skyrmions in frustrated ferromagnets, Phys. Rev. X9, 041063 (2019)

    2019

  5. [13]

    Peruzzo, J

    A. Peruzzo, J. McClean, P. Shadbolt, M.-H. Yung, X.-Q. Zhou, P. J. Love, A. Aspuru-Guzik, and J. L. O’Brien, A variational eigenvalue solver on a photonic quantum processor, Nature Communications5, 4213 (2014)

    2014

  6. [14]

    Tilly, H

    J. Tilly, H. Chen, S. Cao, D. Picozzi, K. Setia, Y. Li, E. Grant, L. Wossnig, I. Rungger, G. H. Booth, and J. Tennyson, The Variational Quantum Eigensolver: A review of methods and best practices, Phys. Rep.986, 1 (2022), arXiv:2111.05176 [quant-ph]

    2022 arXiv

  7. [15]

    Berg and M

    B. Berg and M. L¨ uscher, Definition and statistical dis- tributions of a topological number in the lattice o(3)σ- model, Nuclear Physics B190, 412 (1981)

    1981

  8. [16]

    Van Oosterom and J

    A. Van Oosterom and J. Strackee, The solid angle of a plane triangle, IEEE Transactions on Biomedical Engi- neeringBME-30, 125 (1983)

    1983

  9. [17]

    Senicourt, J

    V. Senicourt, J. Brown, A. Fleury, R. Day, E. Lloyd, M. P. Coons, K. Bieniasz, L. Huntington, A. J. Garza, S. Matsuura, R. Plesch, T. Yamazaki, and A. Zarib- afiyan, Tangelo: An open-source python package for end-to-end chemistry workflows on quantum computers 10.48550/arXiv.2...

    2022 doi arXiv

  10. [18]

    Kandala, A

    A. Kandala, A. Mezzacapo, K. Temme, M. Takita, M. Brink, J. M. Chow, and J. M. Gambetta, Hardware- efficient variational quantum eigensolver for small molecules and quantum magnets, Nature549, 242 (2017)

    2017

  11. [19]

    Arpack software, https://github.com/opencollab/arpack- ng

  12. [20]

    R. B. Lehoucq, D. C. Sorensen, and C. Yang, ARPACK USERS GUIDE: Solution of Large Scale Eigenvalue Problems by Implicitly Restarted Arnoldi Methods(SIAM, Philadelphia, PA, 1998)

    1998

  13. [21]

    Tsuchimochi, M

    T. Tsuchimochi, M. Taii, T. Nishimaki, and S. L. Ten-no, Adaptive construction of shallower quantum circuits with quantum spin projection for fermionic systems, Phys. Rev. Res.4, 033100 (2022)

    2022

  14. [22]

    L. W. Bertels, H. R. Grimsley, S. E. Economou, E. Barnes, and N. J. Mayhall, Symmetry breaking slows convergence of the adapt variational quantum eigen- solver, Journal of Chemical Theory and Computation18, 6656 (2022)

    2022

  15. [23]

    R. A. Istomin and A. S. Moskvin, Overlap integral for quantum skyrmions, Journal of Experimental and Theo- retical Physics Letters71, 338 (2000)

    2000

  16. [24]

    S. A. D ´ ıaz and D. P. Arovas, Quantum nucleation of skyrmions in magnetic films by inhomogeneous fields, in Memorial Volume for Shoucheng Zhang, Chap. Chapter 2, pp. 19–33

  17. [25]

    Hagemeister, N

    J. Hagemeister, N. Romming, K. von Bergmann, E. Y. Vedmedenko, and R. Wiesendanger, Stability of single skyrmionic bits, Nature Communications6, 8455 (2015)

    2015

  18. [26]

    A. P. Petrovi´ c, C. Psaroudaki, P. Fischer, M. Garst, and C. Panagopoulos, Colloquium: Quantum proper- ties and functionalities of magnetic skyrmions (2024), arXiv:2410.11427 [cond-mat.mes-hall]

    2024 arXiv

Pith tools

Reviewed August 7, 2026 · model on record in the stance chip above.