Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “quantum spin model”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Using Hyperoptimized Tensor Networks and First-Principles Electronic Structure to Simulate the Experimental Properties of the Giant {Mn 84 } Torus

The single-molecule magnet {Mn 84 } is a challenge to theory because of its high nuclearity. Here, we directly compute two experimentally accessible observables, the field-dependent magnetization up to 75 T and the temperature-dependent heat capacity, using parameter-free theory. In particular, we use first-principles calculations to derive short- and long-range exchange interactions and compute the exact partition function of the resulting classical Potts and Ising spin models for all 84 Mn S = 2 spins to obtain observables. The latter computation is made possible by using hyperoptimized tensor network contractions, a technique developed to simulate quantum supremacy circuits. We also synthesize the magnet and measure its heat capacity and magnetization, observing qualitative agreement between theory and experiment and identifying an unusual bump in the heat capacity and a plateau in the magnetization. Our work also identifies some limitations of current theoretical modeling in large magnets, such as sensitivity to small, long-range exchange couplings.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Improving Schrödinger Equation Implementations with Gray Code for Adiabatic Quantum Computers

We reformulate the continuous-space Schrödinger equation in terms of spin Hamiltonians. For the kinetic energy operator, the critical concept facilitating the reduction in model complexity is the idea of position encoding. A binary encoding of position produces a spin-1/2 Heisenberg-like model and yields exponential improvement in space complexity when compared to classical computing. Encoding with a binary reflected Gray code (BRGC), and a Hamming-distance-2 Gray code (H2GC) reduces the model complexity down to the 𝑋⁢𝑍 and transverse Ising model, respectively. For 𝐴 qubits BRGC yields 2 𝐴 positions and is reduced to its 2-local form with O⁡(𝐴) ancillary qubits. H2GC yields 2 𝐴/2+1 positions with O⁡(𝐴 2 ) three-local penalty terms. We also identify the bijective mapping between diagonal unitaries and the Walsh series, producing the mapping of any real potential to a series of 𝑘 -local Ising models through the fast Walsh transform. Finally, in a finite volume, we provide some numerical evidence to support the claim that the total time needed for adiabatic evolution is protected by the infrared cutoff of the system. As a result, initial state preparation from a free-field wave function to an interacting system is expected to exhibit polynomial time complexity with volume and constant scaling with respect to lattice discretization for all encodings. For H2GC, if the evolution starts with the transverse Hamiltonian due to hardware restrictions, then penalties are dynamically introduced such that the low-lying spectrum reproduces the energy levels of the Laplacian. The adiabatic evolution of the penalty Hamiltonian is therefore sensitive to the ultraviolet scale. It is expected to exhibit polynomial time complexity with lattice discretization, or exponential time complexity with respect to the number of qubits given a fixed volume.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Static magnetic order with strong quantum fluctuations in spin-1/2 honeycomb magnet Na 2 Co 2 TeO 6

Kitaev interactions, arising from the interplay of frustration and bond anisotropy, can lead to strong quantum fluctuations and, in an ideal case, to a quantum-spin-liquid state. However, in many nonideal materials, spurious non-Kitaev interactions typically promote a zigzag antiferromagnetic order in the d-orbital transition-metal compounds. Here, by combining neutron scattering with muon-spin rotation and relaxation techniques, we provide mechanism insights into the exotic properties of Na 2 Co 2 TeO 6 , a candidate material of the Kitaev model. Below T N , the zero-field muon-spin relaxation rate becomes almost constant (~0.45 μs –1 ). We attribute this temperature-independent relaxation rate to the strong quantum fluctuations, as well as to the frustrated Kitaev interactions. As the magnetic field increases, neutron scattering data indicate a broader spin-wave excitation at the K-point. Therefore, quantum fluctuations seem not only robust but are even enhanced by the applied magnetic field. Our findings provide valuable hints for understanding the onset of the quantum-spin-liquid state in Kitaev materials.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Quantum Phase Engineering of Two-Dimensional Post-Transition Metals by Substrates: Toward a Room-Temperature Quantum Anomalous Hall Insulator

In this work, we propose a new strategy to engineer topological and magnetic properties of two-dimensional (2D) hexagonal lattices consisting of post-transition metals. Our first-principles calculations demonstrate that substrates serve as templates to form 2D lattices with high thermodynamic stability, where their topological properties as well as magnetic properties sensitively change as a function of lattice constants, i.e., the system undergoes a first-order phase transition from nonmagnetic to ferromagnetic state above a critical lattice constant. Consequently, substrates can be used to explore versatile magnetic, electronic, and quantum topological properties. We establish phase diagrams of versatile quantum phases in terms of exchange coupling and spin–orbit coupling effectively tuned by the lattice constants. We further reveal the first room-temperature quantum anomalous Hall (QAH) effect, i.e., Sn on 2√3 × 2√3 graphane is a QAH insulator with a large spin–orbit coupling gap of ~0.2 eV and a Curie temperature of ~380 K by using the 2D anisotropic Heisenberg model.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Variational quantum simulation of the critical Ising model with symmetry averaging

Here we investigate the use of deep multiscale entanglement renormalization ansatz (DMERA) circuits as a variational ansatz. We use the exactly solvable one-dimensional critical transverse-field Ising model as a test bed. Numerically exact simulation of the quantum circuit ansatz can in this case be carried out to hundreds of qubits by exploiting efficient classical algorithms for simulating matchgate circuits. We find that, for this system, the DMERA strongly outperforms a standard quantum approximate optimization algorithm (QAOA)–style ansatz, and that a major source of systematic error in correlation functions approximated using the DMERA is the breaking of the translational and Kramers-Wannier symmetries of the transverse-field Ising model. We are able to reduce this error by up to four orders of magnitude by symmetry averaging, without incurring additional cost in qubits or circuit depth. Here, we propose that this technique for mitigating systematic error could be applied to noisy intermediate-scale quantum (NISQ) simulations of physical systems with other symmetries.

1-dimensional spin chains↗
DOE OSTI2024

Out-of-plane magnetic phase diagram of the Kitaev quantum magnet Na 2 Co 2 TeO 6

We have investigated the magnetic properties and mapped out the phase diagram of the honeycomb magnet Na 2⁢ Co 2 ⁢TeO 6 with Co 3d 7 in out-of-plane magnetic fields. This material has previously been proposed to show nearest-neighbor Kitaev interactions between Co spins and maybe even Kitaev quantum spin liquid behavior in high fields. At 0.01 T, we observe a thermal phase transition at T N = 27K, transitioning from a paramagnetic state to a canonical ferrimagnetic state. Upon increasing the magnetic field, a spin-floplike phase transition occurred before saturation of J = 1/2 between 10 K and T N . Below 10 K, a peak-dip-peak structure emerges between 10 and 17 T in the magnetic susceptibility (d⁢M/⁢dH) before the magnetic saturation, reminiscent of magnetization plateau behavior. The measurement of the magnetocaloric effect also shows dip-peak-dip behavior in this field range. Our data can be explained by an XXZ model with a single ion anisotropy and possibly small Kitaev and Γ exchange interactions. Finally, we also determined the magnetization saturation field that helps constrain the energy scale of the exchange interactions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Auxiliary-state-facilitated phase synchronization phenomena in isolated spin systems

Extending classical synchronization to the quantum domain is of great interest both from the fundamental physics point of view and with a view toward quantum technology applications. This work characterizes phase synchronization of an effective spin-1 system, which is realized by coupling three quantum states with infinite lifetime to auxiliary excited states that have a finite lifetime. Integrating out the excited states, the effective spin-1 model features coherent and incoherent effective couplings. The following are our key findings. (i) Phase synchronization can be controlled by adjusting the phases of the couplings to the excited states. (ii) Unlike in the paradigmatic spin-1 system studied in the literature, where the dissipative couplings describe decay into the limit-cycle state, the effective spin-1 model investigated in this work is governed by a competition between dissipative decay into and out of the limit-cycle state, with the dissipative decay out of the limit-cycle state playing a critical role. (iii) We identify a parameter regime where phase synchronization of the effective spin-1 system is, in the absence of coherent effective couplings, governed entirely by the effective dissipators. The effective spin-1 model is benchmarked through comparisons with master-equation calculations for the full Hilbert space. Physical insights are gained through analytical perturbation theory calculations. In conclusion, our findings, which are expected to hold for a broad class of energy-level and coupling schemes, are examined using hyperfine states of 87 Rb as an example system.

Nonlinear optics↗
DOE OSTI2024

Programmable Heisenberg interactions between Floquet qubits

Abstract The trade-off between robustness and tunability is a central challenge in the pursuit of quantum simulation and fault-tolerant quantum computation. In particular, quantum architectures are often designed to achieve high coherence at the expense of tunability. Many current qubit designs have fixed energy levels and consequently limited types of controllable interactions. Here by adiabatically transforming fixed-frequency superconducting circuits into modifiable Floquet qubits, we demonstrate an XXZ Heisenberg interaction with fully adjustable anisotropy. This interaction model can act as the primitive for an expressive set of quantum operations, but is also the basis for quantum simulations of spin systems. To illustrate the robustness and versatility of our Floquet protocol, we tailor the Heisenberg Hamiltonian and implement two-qubit iSWAP, CZ and SWAP gates with good estimated fidelities. In addition, we implement a Heisenberg interaction between higher energy levels and employ it to construct a three-qubit CCZ gate, also with a competitive fidelity. Our protocol applies to multiple fixed-frequency high-coherence platforms, providing a collection of interactions for high-performance quantum information processing. It also establishes the potential of the Floquet framework as a tool for exploring quantum electrodynamics and optimal control.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Variational Simulation of the Lipkin-Meshkov-Glick Model on a Neutral Atom Quantum Computer

We simulate the Lipkin-Meshkov-Glick model using the variational-quantum-eigensolver algorithm on a neutral atom quantum computer. We test the ground-state energy of spin systems with up to 15 spins. Two different encoding schemes are used: an individual spin encoding where each spin is represented by one qubit, and an efficient Gray code encoding scheme that only requires a number of qubits that scales with the logarithm of the number of spins. This more efficient encoding, together with zero-noise extrapolation techniques, is shown to improve the fidelity of the simulated energies with respect to exact solutions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Beyond-classical computation in quantum simulation

Quantum computers hold the promise of solving certain problems that lie beyond the reach of conventional computers. However, establishing this capability, especially for impactful and meaningful problems, remains a central challenge. Here, we show that superconducting quantum annealing processors can rapidly generate samples in close agreement with solutions of the Schrödinger equation. We demonstrate area-law scaling of entanglement in the model quench dynamics of two-, three-, and infinite-dimensional spin glasses, supporting the observed stretched-exponential scaling of effort for matrix-product-state approaches. We show that several leading approximate methods based on tensor networks and neural networks cannot achieve the same accuracy as the quantum annealer within a reasonable time frame. Thus, quantum annealers can answer questions of practical importance that may remain out of reach for classical computation.

King, Andrew D. [D-Wave Quantum Inc., Burnaby, BC ↗
DOE OSTI2026

Classical combinatorial optimization scaling for random Ising models on 2D heavy-hex graphs

Motivated by near term quantum computing hardware limitations, combinatorial optimization problems that can be addressed by current quantum algorithms and noisy hardware with little or no overhead are used to probe capabilities of quantum algorithms such as the quantum approximate optimization algorithm. In this study, a specific class of near term quantum computing hardware defined combinatorial optimization problems, Ising models on heavy-hex graphs both with and without geometrically local cubic terms, are examined for their classical computational hardness via empirical computation time scaling quantification. Specifically the time-to-solution (TTS) metric using the classical heuristic simulated annealing is measured for finding optimal variable assignments (ground states), as well as the time required for the optimization software Gurobi to find an optimal variable assignment. Because of the sparsity of these Ising models, the classical algorithms are able to find optimal solutions efficiently even for large instances (i.e. 100 000 spin variables). The Ising models both with and without geometrically local cubic terms exhibit average-case linear-time or weakly quadratic scaling when solved exactly using Gurobi, and the Ising models with no cubic terms show evidence of exponential-time TTS scaling when sampled using simulated annealing. These findings point to the necessity of developing and testing more complex, namely more densely connected, optimization problems in order for quantum computing to ever have a practical advantage over classical computing. Our results are another illustration that different classical algorithms can indeed have exponentially different running times, thus making the identification of the best practical classical technique important in any quantum computing vs. classical computing comparison.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Magnetic devil's staircaselike behavior in quasiperiodic qubit lattices

The devil's staircase (DS) phenomenon is a fractal response of magnetization to external fields, traditionally observed in periodic ferromagnetic systems, where the commensurability between spin arrangements, lattice parameters, and external magnetic fields governs abrupt changes in magnetization. Its occurrence in aperiodic, fractal-type systems has remained largely unexplored, despite their natural compatibility with such phenomena. Using a quantum annealing device, we uncover a wealth of abrupt magnetic transitions between spin manifolds driven by increasing external magnetic fields within a simple yet effective Ising-model framework. In contrast to periodic systems, where DS arises from long-range competing interactions, our findings reveal that short-range, purely antiferromagnetic couplings in aperiodic geometries produce equally rich ground-state magnetization patterns. Here, we demonstrate that while magnetic textures are determined by the lattice size, their formation remains remarkably robust and independent of scale, with commensurability emerging locally. Our results challenge the prevailing view that DS behavior is limited to periodic systems and establish quasiperiodic geometries as a natural host for this phenomenon.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Trotter Errors from Dynamical Structural Instabilities of Floquet Maps in Quantum Simulation

We study the behavior of errors in the quantum simulation of spin systems with long-range multibody interactions resulting from the Trotter-Suzuki decomposition of the time-evolution operator. We identify a regime where the Floquet operator underlying the Trotter decomposition undergoes sharp changes even for small variations in the simulation step size. This results in a time evolution operator that is very different from the dynamics generated by the targeted Hamiltonian, which leads to a proliferation of errors in the quantum simulation. These regions of sharp change in the Floquet operator, referred to as structural instability regions, appear typically at intermediate Trotter step sizes and in the weakly interacting regime, and are thus complementary to recently revealed quantum chaotic regimes of the Trotterized evolution [L. M. Sieberer et al. npj Quantum Inf. 5, 78 (2019); M. Heyl, P. Hauke, and P. Zoller, Sci. Adv. 5, eaau8342 (2019)]. We characterize these structural instability regimes in p-spin models, transverse-field Ising models with all-to-all p-body interactions, and analytically predict their occurrence based on unitary perturbation theory. We further show that the effective Hamiltonian associated with the Trotter decomposition of the unitary time-evolution operator, when the Trotter step size is chosen to be in the structural instability region, is very different from the target Hamiltonian, which explains the large errors that can occur in the simulation in the regions of instability. These results have implications for the reliability of near-term gate-based quantum simulators, and reveal an important interplay between errors and the physical properties of the system being simulated.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Plasmons in Z 2 topological insulators

Here, we study plasmonic excitations in the Kane-Mele model, a two-dimensional Z 2 topological insulator on the honeycomb lattice, using the random phase approximation (RPA). In the topologically nontrivial phase, the model has conducting edge states that traverse the bulk energy gap and display spin-momentum locking. Such a state of matter is called the quantum spin hall (QSH) phase, which is robust against time-reversal (TR) invariant perturbations. We find that in the QSH phase, gapless spin-polarized plasmons can be excited on the edges of the system. The propagation of these plasmons is chiral for each individual spin component and shows spin-momentum locking for both spin components on the same edge. Moreover, we study the effect of external magnetic fields on the gapless edge plasmons. Specifically, out-of-plane magnetic fields delocalize edge plasmons propagating in one direction without affecting the other one, while an in-plane magnetic field can be applied to selectively excite a specific spin-plasmon branch with proper doping or gating to the system. Our findings may have potential applications in novel plasmonic and spintronic devices. We also investigate plasmons in the Kane-Mele model on a finite-sized diamond-shaped nanoflake and observe low-energy plasmons circulating the boundary of the material.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Time-reversal symmetry-breaking charge order in a kagome superconductor

The kagome lattice 1 , which is the most prominent structural motif in quantum physics, benefits from inherent non-trivial geometry so that it can host diverse quantum phases, ranging from spin-liquid phases, to topological matter, to intertwined orders 2,3,4,5,6,7,8 and, most rarely, to unconventional superconductivity 6,9 . Recently, charge sensitive probes have indicated that the kagome superconductors AV 3 Sb 5 (A = K, Rb, Cs) 9,10,11 exhibit unconventional chiral charge order 12,13,14,15,16,17,18,19, which is analogous to the long-sought-after quantum order in the Haldane model 20 or Varma model 21 . However, direct evidence for the time-reversal symmetry breaking of the charge order remains elusive. Here we use muon spin relaxation to probe the kagome charge order and superconductivity in KV 3 Sb 5 . We observe a noticeable enhancement of the internal field width sensed by the muon ensemble, which takes place just below the charge ordering temperature and persists into the superconducting state. Notably, the muon spin relaxation rate below the charge ordering temperature is substantially enhanced by applying an external magnetic field. We further show the multigap nature of superconductivity in KV 3 Sb 5 and that the T c /λ$^{-2}_{ab}$ ratio (where T c is the superconducting transition temperature and λ ab is the magnetic penetration depth in the kagome plane) is comparable to those of unconventional high-temperature superconductors. Our results point to time-reversal symmetry-breaking charge order intertwining with unconventional superconductivity in the correlated kagome lattice.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Non-equilibrium critical scaling and universality in a quantum simulator

Universality and scaling laws are hallmarks of equilibrium phase transitions and critical phenomena. However, extending these concepts to non-equilibrium systems is an outstanding challenge. Despite recent progress in the study of dynamical phases, the universality classes and scaling laws for non-equilibrium phenomena are far less understood than those in equilibrium. In this work, using a trapped-ion quantum simulator with single-spin resolution, we investigate the non-equilibrium nature of critical fluctuations following a quantum quench to the critical point. We probe the scaling of spin fluctuations after a series of quenches to the critical Hamiltonian of a long-range Ising model. With systems of up to 50 spins, we show that the amplitude and timescale of the post-quench fluctuations scale with system size with distinct universal critical exponents, depending on the quench protocol. While a generic quench can lead to thermal critical behavior, we find that a second quench from one critical state to another (i.e. a double quench) results in a new universal non-equilibrium behavior, identified by a set of critical exponents distinct from their equilibrium counterparts. Our results demonstrate the ability of quantum simulators to explore universal scaling beyond equilibrium.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Spin liquid to spin glass crossover in the random quantum Heisenberg magnet

We study quantum SU($M$) spins with all-to-all and random Heisenberg exchange interactions of root-mean-square strength $J$. The $M \rightarrow \infty$ model has a spin liquid ground state with the spinons obeying the equations of the Sachdev-Ye-Kitaev (SYK) model. Numerical studies of the SU(2) model with $S=1/2$ spins show spin glass order in the ground state, but also display SYK spin liquid behavior in the intermediate frequency spin spectrum. We employ a $1/M$ expansion to describe the crossover from fractionalized fermionic spinons to a confining spin glass state with weak spin glass order $q_{EA}$. The SYK spin liquid behavior persists down to a frequency $\omega_\ast \sim J q_{EA}$, and for $\omega < \omega_\ast$, the spectral density is linear in $\omega$, thus quenching the extensive zero temperature entropy of the spin liquid. Furthermore, the linear $\omega$ spectrum is qualitatively similar to that obtained earlier using bosonic spinons for large $q_{EA}$. We argue that the extensive SYK spin liquid entropy is transformed as $T \rightarrow 0$ to an extensive complexity of the spin glass state.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗