{"id":"87e14775-52c4-40bd-8924-53ee77783e81","arxiv_id":"2505.19808","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"A simulator-based quantum eigensolver reveals field-driven discontinuities in a small Heisenberg ferromagnet, which the authors read as hints of zero-temperature skyrmion-like magnetic patterns.","lead":"The authors used a variational quantum eigensolver on a simulator to study a small magnetic lattice with exchange and Dzyaloshinskii-Moriya interactions. They report sharp jumps in energy, magnetization, and topological charge as the magnetic field changes, which they interpret as hints of zero-temperature skyrmion-like magnetic patterns.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The claimed field-driven transitions rest on the VQE ansatz reaching the exact ground state; no overlap or convergence check is shown at the transition fields, so the discontinuities could be variational artifacts.","rationale":"The reader's weakest assumption correctly identifies the load-bearing point: the HEA must reach the true ground state at every field value, especially near the transition, and the paper gives no evidence for this. My reading of the manuscript confirms that the only accuracy statement is a vague reference to agreement with Lanczos, with no numerical comparison shown. The proposed exact-diagonalization test is the decisive check because N = 16 is small enough that the true ground state is readily computable, and it would immediately distinguish a physical level crossing from a variational jump. I therefore agree with the reader's CONDITIONAL verdict; no change is needed, but the condition should be the overlap/energy comparison at the transition fields. I considered the non-integer topological charge and the extrapolated efficiency comparison as alternative concerns, but both are explicitly hedged or peripheral to the main physical claim. The absence of a ground-state fidelity check is the single most load-bearing weakness.","tokens_in":6902,"tokens_out":5655,"duration_ms":65185,"concrete_test":"Run exact Lanczos diagonalization for the same N = 16 square lattice, open boundary conditions, and the same parameters as Figs. 3 and 5, scanning B_z densely across the claimed transitions (e.g., B_z = 1.86-1.91 in steps of 0.002 for D parallel, and B_z = 1.20-1.25 for D perpendicular). Compute the exact ground-state energy, m_x, and the triangulated Q, and for every field value compute the squared overlap |<ψ_VQE|ψ_exact>|^2 and the energy difference E_VQE - E_exact. If the exact results reproduce the discontinuities and the VQE overlaps stay above 0.99, the concern is resolved; if the exact energy is continuous or the overlap drops sharply near the claimed transition, the observed jumps are variational artifacts.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central physics conclusion is the discontinuous jump in E, m_x, and Q at B_z = 1.884 (Fig. 3) and B_z = 1.222 (Fig. 5), which is interpreted as a transition to skyrmion-like ordering. These results come solely from the hardware-efficient ansatz with 6N variational parameters and a single entangler. The paper says \"the only criterion is an agreement with the Lanczos-method results within the numerical accuracy,\" but it never reports overlaps, energy differences, or convergence data, especially at the transition fields. Since the model is finite and the DMI term breaks total-S_z conservation, an exact finite-system ground-state energy is generically continuous; a sharp VQE jump therefore needs independent confirmation before being assigned a physical meaning. For N = 16, exact Lanczos diagonalization is completely feasible, so the absence of that comparison is the weakest link: a variational collapse between two local minima of the HEA could produce exactly the kind of jump shown in Figs. 3 and 5 without any corresponding feature in the true ground state. The non-integer topological charges are acknowledged by the authors and are secondary; the efficiency claim in Fig. 2 is based on extrapolated fits from five lattice sizes, but even if that claim were wrong the skyrmion prediction would stand, whereas a variational artifact would invalidate it.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":7182,"tokens_out":3489,"duration_ms":39512,"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":[{"comment":"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.","section":"Section 3, Figs. 3 and 5"},{"comment":"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.","section":"Section 2, Fig. 2"},{"comment":"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.","section":"Section 4 and Fig. 7"}],"minor_comments":[{"comment":"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.","section":"Eq. (2) and following text"},{"comment":"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.","section":"Ansatz definition, Section 2"},{"comment":"The caption contains the stray word 'colorblack' between 'and' and 'J_perp', which should be removed.","section":"Fig. 5 caption"},{"comment":"There is a typo 'Zeemnan interaction' and an awkward 'Hence. a further optimization' with a period instead of a comma; both should be fixed.","section":"Section 2 text"},{"comment":"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.","section":"Eq. (1) and Section 2"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the scope of quantum magnetism and VQE methods, and the core idea is interesting. The decisive issue is verifiability: the N=16 system is small enough that exact Lanczos results are cheap, yet they are not reported at the transition fields. I would urge the editor to require those exact comparisons as a condition of acceptance, since they directly test the variational-artifact interpretation. The efficiency claim also needs a more careful statistical presentation."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe one thing to know: this is a proof-of-concept VQE-on-simulator study of a 16-site DMI Heisenberg model, and the advertised field-driven transition is plausible but not yet substantiated. The authors never show that their hardware-efficient ansatz actually reaches the true ground state at the transition fields, so the observed discontinuities could be variational artifacts. For N=16, exact Lanczos is trivial, so the omission is the paper's weakest point.\n\nThe new content is real: they run VQE on a noiseless simulator for this model, compare runtime scaling against Lanczos, and produce a field sweep with a clear jump in energy, magnetization, and topological charge at B_z ~ 1.884 (D parallel) and ~1.222 (D perpendicular). The magnetization patterns before and after the jump do resemble the start of Bloch- or Néel-type textures. They also honestly acknowledge the topological charges are non-integer and that a 4x4 lattice cannot host a proper skyrmion. That intellectual honesty deserves credit.\n\nThe soft spots, in order. First, the load-bearing one: no overlap, energy difference, or convergence data between VQE and the exact ground state is shown for any field, especially not at the transitions. With 6N parameters and a single entangler, the HEA can easily collapse into a local minimum and produce a spurious jump. Second, the VQE efficiency claim is based on runtime fits from only five lattice sizes (including triangular ones) and extrapolated to N>17; suggestive, but not a proof. Third, the final paragraph says the results 'predict' skyrmion-like structures, which oversells what a 16-site lattice with fractional Q can support. In a finite system the exact ground-state energy is generically continuous, so a sharp jump needs numerical confirmation before being interpreted as a phase transition.\n\nWho gets value: readers working on VQE applications to quantum magnetism and on the quantum-skyrmion question. The paper deserves serious peer review, but the referee should demand the missing exact-diagonalization checks at the transition fields, error bars, and a more cautious conclusion. With those additions, this would be a useful benchmark. As it stands, the transition claim is unverified.\n\nRecommendation: send to peer review; the gap is fixable and the topic is worth a second look.","headline":"VQE study of DMI Heisenberg model with a plausible but unverified transition claim; missing exact-diagonalization check at the transition is the key gap.","tokens_in":7745,"tokens_out":4646,"would_cite":false,"duration_ms":44113,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["74.20.pq","74.25.Kc","71.15.Mb"],"model":"deepseek-v4-flash","headline":"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.","keywords":["variational quantum eigensolver","Dzyaloshinskii-Moriya interaction","Heisenberg ferromagnet","quantum skyrmions","topological charge","hardware-efficient ansatz","zero-temperature phase transition","spin-1/2 lattice"],"falsifier":"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.","tokens_in":6675,"feed_emoji":"🧲","tokens_out":18145,"duration_ms":154644,"temperature":0.7,"pith_summary":"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.","feed_headline":"Quantum algorithm finds zero-temperature skyrmion-like spin textures","feed_subtitle":"On a simulator, the algorithm sees a sharp jump into skyrmion-like order and beats classical solvers past 17 sites.","key_machinery":"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.","core_discovery":"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$.","pith_inferences":["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."],"forward_implications":["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."],"supporting_citations":[{"why":"Introduces the variational quantum eigensolver, the algorithm the paper applies to the spin Hamiltonian.","marker":"[13]"},{"why":"Supplies the hardware-efficient ansatz circuit of Euler rotations around one entangler.","marker":"[18]"},{"why":"Provides the open-source VQE implementation and the noiseless simulator backend used for all calculations.","marker":"[17]"},{"why":"Predicts that quantum fluctuations stabilize skyrmion textures, a basis for expecting zero-temperature quantum skyrmions.","marker":"[7]"},{"why":"Computes quantum skyrmion lattices in Heisenberg ferromagnets and sets the parameter regime used here.","marker":"[8]"},{"why":"Defines the Dzyaloshinskii-Moriya interaction, the antisymmetric exchange term that generates skyrmion-like order.","marker":"[10, 11]"},{"why":"Gives the topological-charge definition used to identify skyrmion-like patterns.","marker":"[5]"},{"why":"Provides the lattice definition of a topological winding number that underlies the triangulated Q computation.","marker":"[15]"},{"why":"Gives the solid-angle formula for a triangle used to evaluate the topological charge on the discrete mesh.","marker":"[16]"},{"why":"Shows that a similar energy landscape yields long lifetimes for single skyrmionic bits, supporting the stability claim.","marker":"[25]"}],"fun_headline_variants":["VQE reveals quantum skyrmion-like transitions in Heisenberg model","Quantum solver beats classical past 17 sites, sees skyrmion jump","Skyrmion precursors at zero temperature from quantum eigensolver","Sharp magnetization jump hints at stable quantum skyrmions","VQE outruns Lanczos, spots magnetic transition in spin lattice"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["VQE reveals quantum skyrmion-like transitions in Heisenberg model","Quantum solver beats classical past 17 sites, sees skyrmion jump","Skyrmion precursors at zero temperature from quantum eigensolver","Sharp magnetization jump hints at stable quantum skyrmions","VQE outruns Lanczos, spots magnetic transition in spin lattice"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000218,"raw_usage":{"total_tokens":1458,"prompt_tokens":980,"completion_tokens":478,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":596,"completion_tokens_details":{"reasoning_tokens":389}},"tokens_in":596,"tokens_out":478,"duration_ms":4891,"temperature":1.0,"reasoning_tokens":389,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T14:05:50.949394+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[{"cited_title":"Rold´ an-Molina, M","cited_arxiv_id":null,"evidence_quote":"Predicts that quantum fluctuations stabilize skyrmion textures, a basis for expecting zero-temperature quantum skyrmions."},{"cited_title":"Haller, S","cited_arxiv_id":null,"evidence_quote":"Computes quantum skyrmion lattices in Heisenberg ferromagnets and sets the parameter regime used here."},{"cited_title":"del Ser, I","cited_arxiv_id":null,"evidence_quote":"Gives the topological-charge definition used to identify skyrmion-like patterns."},{"cited_title":"Van Oosterom and J","cited_arxiv_id":null,"evidence_quote":"Gives the solid-angle formula for a triangle used to evaluate the topological charge on the discrete mesh."},{"cited_title":"Hagemeister, N","cited_arxiv_id":null,"evidence_quote":"Shows that a similar energy landscape yields long lifetimes for single skyrmionic bits, supporting the stability claim."}],"review_version":1}