Search NASA⌕ Search

SEARCH · Search NASA

Results for “Quantum spin models”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Sunny.jl: A Julia Package for Spin Dynamics

Sunny is a Julia package designed to serve the needs of the quantum magnetism community. It supports the specification of a very broad class of spin models and a diverse suite of numerical solvers. These include powerful methods for simulating spin dynamics both in and out of equilibrium. Uniquely, it features a broad generalization of classical and semiclassical approaches to SU(N) coherent states, which is useful for studying systems exhibiting strong spin-orbit coupling or local entanglement effects. Sunny also offers a well-developed framework for calculating the dynamical spin structure factor, enabling direct comparison with scattering experiments. Ease of use is a priority, with tools for symmetry-guided modeling and interactive visualization.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Spin Polarization from Circularly Polarized Light Induced Charge Transfer

Here, we show how a spin polarization can be generated through the photoinduced electron transfer of an achiral donor–acceptor complex following chiral light excitation. In particular, we illustrate the basic energetic and symmetry requirements for chirality induced spin selectivity where the chirality emerges from the electronic degrees of freedom following excitation with circularly polarized light. We study this effect in a simple model of a metalloporphyrin complex with an axial acceptor ligand using quantum mechanical rate theories and numerical simulations. We find that the spin polarization emerges due to the selective excitation of a ring current within the porphryin, breaking the degeneracy of the two degenerate spin states. The resultant spin polarization increases with the spin orbit coupling between the metal in the porphyrin and the axial ligand, and is transient, with a lifetime dependent on the rate of dephasing from the Jahn–Teller distortion mode. This proposed effect should be observable in spin-resolved photoemission spectroscopy.

charge transfer↗
DOE OSTI2025

Field-tunable BKT and quantum phase transitions in spin-$\frac{1}{2}$ triangular lattice antiferromagnet

Quantum magnetism is one of the most active fields for exploring exotic phases and phase transitions. The recently synthesized Na 2 BaCo(PO 4 ) 2 (NBCP) is an ideal material incarnation of the spin-$\frac{1}{2}$ easy-axis triangular lattice antiferromagnet (TLAF). Experimental evidence shows that NBCP hosts the spin supersolid state with a giant magnetocaloric effect. Theory further predicts that magnetic fields can drive NBCP through Berezinskii-Kosterlitz-Thouless (BKT) and other richer quantum phase transitions. However, detecting these transitions is challenging, as they onset at ultralow temperatures near 60 mK and require high magnetization sensitivity. Using a newly developed gradient force magnetometer in a dilution refrigerator, we mapped the magnetic susceptibility phase diagram down to 30 mK. Our results provide a more comprehensive and accurate understanding of BKT melting of spin supersolidity and several field-tunable quantum phase transitions, which establish NBCP as a model platform for frustrated magnetism and highlight potential applications of its giant magnetocaloric effects.

BKT transition↗
NASA NTRS2013

Room-Temperature Operation of Index-Coupled Distributed-Feedback 4.75 Micron Quantum Cascade Lasers Fabricated Without Epitaxial Regrowth

The Mars Science Laboratory (MSL) was spin-stabilized during its cruise to Mars. We discuss the effects of spin on the radiometric data and how the orbit determination team dealt with them. Additionally, we will discuss the unplanned benefits of detailed spin modeling including attitude estimation and spacecraft clock correlation.

quantum cascade lasers↗
DOE OSTI2024

Exactly solvable model of light-scattering errors in quantum simulations with metastable trapped-ion qubits

Here, we analytically solve a model for light scattering in Ising dynamics of metastable atomic qubits, generalizing the approach of Foss-Feig et al. to include leakage outside the qubit manifold. We analyze the influence of these fundamental errors in simulations of proposed experiments with metastable levels of 40 Ca + ions in a Penning trap. We find that “effective magnetic fields” generated by leaked qubits have significant impacts on spin-spin correlation functions for Greenberger-Horne-Zeilinger state preparation or for quantum simulations with strong coupling, while spin squeezing uses a much weaker coupling and is largely insensitive to the simulated leakage errors, even with a few hundred ions. Our theory and results are expected to be useful in modeling a variety of metastable qubit experiments in the future.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Determining the N -Representability of a Reduced Density Matrix via Unitary Evolution and Stochastic Sampling

The N-representability problem consists in determining whether, for a given p-body matrix, there exists at least one N-body density matrix from which the p-body matrix can be obtained by contraction, that is, if the given matrix is a p-body reduced density matrix (p-RDM). The knowledge of all necessary and sufficient conditions for a p-body matrix to be N-representable allows the constrained minimization of a many-body Hamiltonian expectation value with respect to the p-body density matrix and, thus, the determination of its exact ground state. However, the number of constraints that complete the N-representability conditions grows exponentially with system size, and hence, the procedure quickly becomes intractable for practical applications. This work introduces a hybrid quantum-stochastic algorithm to effectively replace the N-representability conditions. The algorithm consists of applying to an initial N-body density matrix a sequence of unitary evolution operators constructed from a stochastic process that successively approaches the reduced state of the density matrix on a p-body subsystem, represented by a p-RDM, to a target p-body matrix, potentially a p-RDM. The generators of the evolution operators follow the well-known adaptive derivative-assembled pseudo-Trotter method (ADAPT), while the stochastic component is implemented by using a simulated annealing process. The resulting algorithm is independent of any underlying Hamiltonian, and it can be used to decide whether a given p-body matrix is N-representable, establishing a criterion to determine its quality and correcting it. We apply the proposed hybrid ADAPT algorithm to alleged reduced density matrices from a quantum chemistry electronic Hamiltonian, from the reduced Bardeen–Cooper–Schrieffer model with constant pairing, and from the Heisenberg XXZ spin model. In all cases, the proposed method behaves as expected for 1-RDMs and 2-RDMs, evolving the initial matrices toward different targets.

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

Evidence for field induced quantum spin liquid behavior in a spin-1/2 honeycomb magnet

One of the most important issues in modern condensed matter physics is the realization of fractionalized excitations, such as the Majorana excitations in the Kitaev quantum spin liquid. To this aim, the 3d-based Kitaev material Na 2 Co 2 TeO 6 is a promising candidate whose magnetic phase diagram of B // a* contains a field-induced intermediate magnetically disordered phase within 7.5 T < |B| < 10 T. The experimental observations, including the restoration of the crystalline point group symmetry in the angle-dependent torque and the coexisting magnon excitations and spinon-continuum in the inelastic neutron scattering spectrum, provide strong evidence that this disordered phase is a field induced quantum spin liquid with partially polarized spins. Our variational Monte Carlo simulation with the effective K-J 1 -Γ-Γ'-J 3 model reproduces the experimental data and further supports this conclusion.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Synthesis and Characterization of Metastable Cobalt Honeycomb KCoAsO 4

The Kitaev model has served as a long-sought-after target in the realization of a quantum spin liquid that could host Majorana Fermions. Such non-Abelian anyons could revolutionize quantum computing if properly implemented to overcome decoherence. A 3d 7 electronic configuration, like Co 2+ , has been explored by theory and experimental work to design Kitaev materials. Here, in this study, we report the synthesis of a new cobaltate honeycomb material KCoAsO 4 . The compound is synthesized through a low-temperature solution route and crystallized in space group R¯3 with lattice parameters a = 5.0394(1) and c = 28.6790(1) as determined by neutron powder diffraction. The crystal structure follows motifs similar to those of the honeycomb compound BaCo2(AsO 4 ) 2 but presents differing magnetic behavior. Magnetization/heat capacity measurements on the powder show antiferromagnetic transition T N = 14 K. Two lower-temperature transitions are present in susceptibility at low field that resemble spin reorientations. Magnetization data as a function of field have curvature indicative of metamagnetic behavior below the magnetic ordering temperature, with the magnetic ordering suppressed upon application of a higher magnetic field. Computational studies suggest the presence of a weak nearest-neighbor Kitaev term, K 1 , consistent with related honeycomb cobaltates. Together, the data suggest that this material should present a new platform for developing Kitaev quantum spin liquids.

cobaltate↗
DOE OSTI2026

Quantifying the phase diagram and Hamiltonian of S = 1/2 kagome antiferromagnets: bridging theory and experiment

Spin-1/2 kagome antiferromagnets are leading candidates for realizing quantum spin liquid (QSL) ground states. While QSL ground states are predicted for the pure Heisenberg model, understanding the robustness of the QSL to additional interactions that may be present in real materials is a forefront question in the field. Here we employ large-scale density-matrix renormalization group simulations to investigate the effects of next-nearest neighbor exchange couplings J 2 and Dzyaloshinskii-Moriya interactions D, which are relevant to understanding the prototypical kagome materials herbertsmithite and Zn-barlowite. By utilizing clusters as large as XC12 and extrapolating the results to the thermodynamic limit, we precisely delineate the scope of the QSL phase, which remains robust across an expanded parameter range of J 2 and D. Direct comparison of the simulated static and dynamic spin structure factors with inelastic neutron scattering reveals the parameter space of the Hamiltonians for herbertsmithite and Zn-barlowite, and, importantly, provides compelling evidence that both materials exist within the QSL phase. These results establish a powerful convergence of theory and experiment in this most elusive state of matter.

Jiang, Shengtao [SLAC National Accelerator Laborat↗
DOE OSTI2024

Kondo enabled transmutation between spinons and superconducting vortices: Origin of magnetic memory in 4 H b − TaS 2

Recent experiments [Persky , ] demonstrate a magnetic memory effect in 4 H b − TaS 2 above its superconducting transition temperature, where Abriokosov vortices are spontaneously generated by lowering the temperature at zero magnetic field after field training the normal state. Motivated by the experiment, we propose the chiral quantum spin liquid (QSL) stabilized in the constituent layers of 4 H b − TaS 2 as a mechanism. We model 4 H b − TaS 2 as coupled layers of the chiral QSL and the superconductor. Through the Kondo coupling between the localized moments and conduction electrons, there is mutual transmutation between spinons and vortices during the thermal-cycling process, which yields the magnetic memory effect observed in experiments. We also propose a mechanism to stabilize chiral and nematic superconductivity in 4 H b − TaS 2 through the Kondo coupling of conduction electrons to the chiral QSL. Our picture suggests 4 H b − TaS 2 as an exciting platform to explore the interplay between QSL and superconductivity through the Kondo effect. Published by the American Physical Society 2024

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Krylov winding and emergent coherence in operator growth dynamics

The operator wavefunction provides a fine-grained description of quantum chaos and of the irreversible growth of simple operators into increasingly complex ones. Remarkably, at finite temperature this wavefunction can acquire a phase that increases linearly with the operator’s size, a phenomenon called . Although size winding occurs naturally in a holographic setting, the emergence of a coherent phase in a scrambled operator remains mysterious from the standpoint of a thermalizing quantum many-body system. Here, in this article, we elucidate this phenomenon by introducing the related concept of , whereby the operator wavefunction acquires a phase which winds linearly with the Krylov index. We show that Krylov winding is a generic feature of quantum chaotic systems and is a direct consequence of the universal operator growth bound hypothesis. It gives rise to size winding under two additional conditions: (i) a low-rank mapping between the Krylov and size bases, which ensures phase alignment among operators of the same size, and (ii) the saturation of the "chaos-operator growth" bound 𝜆 𝐿 ≤ 2⁢𝛼 (with 𝜆 𝐿 the Lyapunov exponent and 𝛼 the growth rate), which ensures a linear phase dependence on size. For systems which do not saturate this bound, with ℎ = 𝜆 𝐿 /2⁢𝛼 < 1, the winding with Pauli size ℓ becomes superliner, behaving as ℓ 1/ℎ . We illustrate these results with two classes of microscopic models: the Sachdev-Ye-Kitaev (SYK) model and its variants, and a disordered 𝑘-local spin model.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Unconventional bipartite entanglement in the quantum dimer magnet Yb2Be 2 SiO 7

The quantum dimer magnet, with antiferromagnetic intradimer and interdimer Heisenberg exchange between spin-1/2 moments, is known to host an (|$\uparrow$$\downarrow$$\rangle$ - |$\downarrow$$\uparrow$$\rangle$)/$\sqrt{2}$ singlet ground state when the intradimer exchange is dominant. Rare-earth-based quantum dimer systems with strong spin-orbit coupling offer the opportunity for tuning their magnetic properties by using magnetic anisotropy as a control knob. Here, we present bulk characterization and neutron scattering measurements of the quantum dimer magnet Yb 2 Be 2 SiO 7 . We find that the Yb3+ ions can be described by an effective spin-1/2 model at low temperatures and the system does not show signs of magnetic order down to 50 mK. The magnetization, heat capacity, and neutron spectroscopy data can be well-described by an isolated dimer model with highly anisotropic exchange that stabilizes a singlet ground state with a wavefunction (|$\uparrow$$\uparrow$$\rangle$ - |$\downarrow$$\downarrow$$\rangle$)/$\sqrt{2}$ or (|$\uparrow$$\uparrow$$\rangle$ + |$\downarrow$$\downarrow$$\rangle$)/$\sqrt{2}$. Our results show that strong spin-orbit coupling can induce unusual entangled states of matter in quantum dimer magnets.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Simulating plasma wave propagation on a superconducting quantum chip

Quantum computers may one day enable the efficient simulation of strongly coupled plasmas that lie beyond the reach of classical computation in regimes where quantum effects are important and the scale separation is large. Here, in this article, we take a first step toward efficient simulation of quantum plasmas by demonstrating linear plasma wave propagation on a superconducting quantum chip. Using high-fidelity and highly expressive device-native gates, combined with an error-mitigation technique, we simulate the scattering of laser pulses from inhomogeneous plasmas. Our approach is made feasible by the identification of a suitable local spin model whose excitations mimic plasma waves, and whose circuit implementation requires a lower gate count than other proposed approaches that would require a future fault-tolerant quantum computer. This work opens avenues to study more complicated phenomena that cannot be simulated efficiently on classical computers, such as nonlinear quantum dynamics when strongly coupled plasmas are driven out of equilibrium.

general physics↗
DOE OSTI2024

Quantum subspace expansion in the presence of hardware noise

Finding ground state energies on current quantum processing units (QPUs) using algorithms such as the variational quantum eigensolver (VQE) continues to pose challenges. Hardware noise severely affects both the expressivity and trainability of parameterized quantum circuits, limiting them to shallow depths in practice. Here, we demonstrate that both issues can be addressed by synergistically integrating VQE with a quantum subspace expansion, allowing for an optimal balance between quantum and classical computing capabilities and costs. We perform a systematic benchmark analysis of the iterative quantum-assisted eigensolver in the presence of hardware noise. We determine ground state energies of 1D and 2D mixed-field Ising spin models on noisy simulators and the IBM QPUs ibmq_quito (5 qubits) and ibmq_guadalupe (16 qubits). To maximize accuracy, we propose a suitable criterion to select the subspace basis vectors according to the trace of the noisy overlap matrix. Finally, we show how to systematically approach the exact solution by performing controlled quantum error mitigation based on probabilistic error reduction on the noisy backend fake_guadalupe.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

On the discretization error of the discrete generalized quantum master equation

The transfer tensor method (TTM) [Cerrillo and Cao, Phys. Rev. Lett. 112 , 110401 (2014)] can be considered a discrete-time formulation of the Nakajima–Zwanzig quantum master equation (NZ-QME) for modeling non-Markovian quantum dynamics. A recent paper [Makri, J. Chem. Theory Comput. 21 , 5037 (2025)] raised concerns regarding the consistency of the TTM discretization, particularly a spurious term at the initial time t = 0. Here, this work presents a detailed analysis of the discretization structure of the TTM, clarifying the origin of the initial-time correction and establishing a consistent relationship between the TTM discrete-time memory kernel K N and the continuous-time NZ-QME kernel $\mathscr{K}$( N Δ t ). This relationship is validated numerically using the spin-boson model, demonstrating convergence of reconstructed memory kernels and accurate dynamical evolution as Δ t → 0. While the TTM provides a consistent discretization, we note that alternative schemes are also viable, such as the midpoint derivative/midpoint integral scheme proposed in Makri’s work. The relative performance of various schemes for either computing accurate $\mathscr{K}$( N Δ t ) from exact dynamics or obtaining accurate dynamics from exact $\mathscr{K}$( N Δ t ) warrants further investigation.

Density-matrix↗
DOE OSTI2025

Magnon damping and mode softening in quantum double-exchange ferromagnets

We present a comprehensive analysis of the magnetic excitations and electronic properties of fully quantum double-exchange ferromagnets, i.e. systems where ferromagnetic (FM) ordering emerges from the competition between spin, charge, and orbital degrees of freedom, but without the canonical approximation of using classical localized spins. Specifically, we investigate spin excitations within the Kondo lattice-like model, as well as a two-orbital Hubbard Hamiltonian in proximity to the orbital-selective Mott phase. Computational analysis of the magnon dispersion, damping, and spectral weight within these models reveals unexpected phenomena, such as magnon mode softening and the anomalous decoherence of magnetic excitations as observed in earlier experimental efforts, but explained here without the use of the phononic degrees of freedom. We show that these effects are intrinsically linked to incoherent spectral features near the Fermi level, which arise due to the quantum nature of the local (on-site) triplets. This incoherent spectrum leads to a Stoner-like continuum on which spin excitations scatter, governing magnon lifetime and strongly influencing the dynamical spin structure factor. Our study explores the transition from coherent to incoherent magnon spectra by varying the electron density. Furthermore, we demonstrate that the magnitude of the localized spin mitigates decoherence by suppressing the incoherent spectral contributions near the Fermi level. We also discuss the effective J 1 –J 2 spin Hamiltonian, which can accurately describe the large doping region characterized by the magnon-mode softening. Finally, we show that this behavior is also present in multiorbital models with partially filled orbitals, namely, in systems without localized spin moments, provided that the model is in a strong coupling regime. Our results potentially have far-reaching implications for understanding FM ordering in various multi-band systems. These findings establish a previously unknown direct connection between the electronic correlations of those materials and spin excitations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Two-channel multi-impurity Kondo model: RKKY-induced criticality

We show that the Ruderman-Kittel-Kasua-Yoisida interaction between overscreened spins in two-channel Kondo impurity systems is a relevant perturbation when the number of impurities N is greater than 3 driving the system to a new quantum critical point with anomalous dimensions $\frac{1}{(𝑁+1)}$ for the spin operator and the Sommerfeld coefficient of the specific heat scales as 𝛾 ∼ 𝑇 −$\frac{3}{𝑁+1}$ . The critical point universal properties are relevant to many strong correlation problems, such as impurity placed in a Majorana metal and the multichannel Kondo lattice model of heavy fermion materials. In conclusion, we discuss relevance of our results for cluster DMFT studies of quantum criticality.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Biased degenerate ground-state sampling of small Ising models with converged quantum approximate optimization algorithm

The quantum alternating operator ansatz, a generalization of the quantum approximate optimization algorithm (QAOA), is a quantum algorithm used for approximately solving combinatorial optimization problems. QAOA typically uses the transverse field mixer as the driving Hamiltonian. One of the interesting properties of the transverse field driving Hamiltonian is that it results in nonuniform sampling of degenerate ground states of optimization problems. In this study, we numerically examine the fair sampling properties of the transverse field mixer QAOA, and Grover mixer QAOA (GM-QAOA), which provides theoretical guarantees of fair sampling of degenerate optimal solutions, up to a large enough p such that the mean expectation value converges to an optimal approximation ratio of 1. This comparison is performed with high-quality heuristically computed, but not necessarily optimal, QAOA angles, which give strictly monotonically improving solution quality as p increases. These angles are computed using the Julia based numerical simulation software JuliQAOA. Fair sampling of degenerate ground states is quantified using the Shannon entropy of the ground-state amplitudes distribution. The fair sampling properties are reported on several quantum signature Hamiltonians from previous quantum annealing fair sampling studies. Small random fully connected spin glasses are shown, which exhibit exponential suppression of some degenerate ground states with transverse field mixer QAOA. The transverse field mixer QAOA simulations show that some problem instances clearly saturate the Shannon entropy of 0 with a maximally biased distribution that occurs when the learning converges to an approximation ratio of 1 while other problem instances never deviate from a maximum Shannon entropy (uniform distribution) at any p step. Published by the American Physical Society 2025

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗