Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Ising spins”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Self-consistent mean-field quantum approximate optimization

We introduce a self-consistent mean-field quantum optimization algorithm that approximates the ground state of classical Ising Hamiltonians. The algorithm decomposes the problem into independent subproblems and treats the interactions between them in a mean-field manner. These interactions are captured by a common environment, constructed self-consistently through a variational quantum circuit, and which modifies the subproblems to account for mutual influence while maintaining computational independence. Consequently, subproblems can be solved individually, avoiding the computational cost of the full problem. We explore the properties of the generated environment and assess the algorithm's performance through extensive numerical simulations on Sherrington-Kirkpatrick spin glasses. Furthermore, we apply it experimentally to a weighted maximum clique problem applied to molecular docking. This framework enables the solution of problems that would otherwise exceed the qubit and gate counts of current quantum hardware.

Dupont, Maxime [Rigetti Computing] (ORCID:00000001↗
DOE OSTI2026

Self-consistent mean-field quantum approximate optimization

We introduce a self-consistent mean-field quantum optimization algorithm that approximates the ground state of classical Ising Hamiltonians. The algorithm decomposes the problem into independent subproblems and treats the interactions between them in a mean-field manner. These interactions are captured by a common environment, constructed self-consistently through a variational quantum circuit, and which modifies the subproblems to account for mutual influence while maintaining computational independence. Consequently, subproblems can be solved individually, avoiding the computational cost of the full problem. We explore the properties of the generated environment and assess the algorithm's performance through extensive numerical simulations on Sherrington-Kirkpatrick spin glasses. Furthermore, we apply it experimentally to a weighted maximum clique problem applied to molecular docking. This framework enables the solution of problems that would otherwise exceed the qubit and gate counts of current quantum hardware.

Dupont, Maxime [Rigetti Computing] (ORCID:00000001↗
DOE OSTI2026

Self-consistent mean-field quantum approximate optimization

We introduce a self-consistent mean-field quantum optimization algorithm that approximates the ground state of classical Ising Hamiltonians. The algorithm decomposes the problem into independent subproblems and treats the interactions between them in a mean-field manner. These interactions are captured by a common environment, constructed self-consistently through a variational quantum circuit, and which modifies the subproblems to account for mutual influence while maintaining computational independence. Consequently, subproblems can be solved individually, avoiding the computational cost of the full problem. We explore the properties of the generated environment and assess the algorithm's performance through extensive numerical simulations on Sherrington-Kirkpatrick spin glasses. Furthermore, we apply it experimentally to a weighted maximum clique problem applied to molecular docking. This framework enables the solution of problems that would otherwise exceed the qubit and gate counts of current quantum hardware.

Dupont, Maxime [Rigetti Computing] (ORCID:00000001↗
NASA NTRS1981

Ion composition and energy distribution during 10 magnetic storms

Data from the plasma composition experiment of ISEE 1 are used to investigate the relative quantities and energy characteristics of H(+), He(++), He(+), and O(+) ions in the near-equatorial magnetosphere during magnetic storm conditions. The ions in the study had energies between 0.1 and 17 keV/e and pitch angles between 45 and 135 deg. The data were obtained during 10 storms, for the most part at or immediately following the peak Dst, covering all major local time sectors and geocentric distances between 2 and 15 earth radii. The ion fluxes are averaged over the spacecraft spin angle and over time for periods ranging from about 20 min close to the earth to more than an hour in most distant regions. The inferred 'isotropic' number densities are characterized by a large to dominant fraction of terrestrial ions throughout the energy range covered. The data are found to be consistent with a terrestrial origin for all of the O(+), most of the He(+), and a large but varying fraction of the H(+), whereas the He(++) and part of the H(+) appear to be of solar wind origin.

Lennartsson, W.↗
NASA NTRS1982

Inductive electric field at the magnetopause

The electric field data for two crossings of the magnetopause by ISEE-1 on November 20, 1977, have been analyzed with high time resolution. In both cases the electric field has a negative dawn-dusk component in the boundary layer, so it must reverse somewhere within the current layer to the positive value outside. If there is a component parallel to the moving magnetopause current it is small, and by no means obvious. In the case of the exit crossing from the boundary layer to the magnetosheath the data show that the electric field vector is turning for about two seconds at roughly the satellite spin rate; this changing direction suggests that the electric field has a curl. Such a curl could be caused by a travelling localized perturbation of the magnetopause surface current associated with impulsive plasma transport through the magnetopause.

Heikkila, W. J.↗
DOE OSTI2024

Confinement and Kink Entanglement Asymmetry on a Quantum Ising Chain

In this work, we explore the interplay of confinement, string breaking and entanglement asymmetry on a 1D quantum Ising chain. We consider the evolution of an initial domain wall and show that, surprisingly, while the introduction of confinement through a longitudinal field typically suppresses entanglement, it can also serve to increase it beyond a bound set for free particles. Our model can be tuned to conserve the number of domain walls, which gives an opportunity to explore entanglement asymmetry associated with link variables. We study two approaches to deal with the non-locality of the link variables, either directly or following a Kramers-Wannier transformation that maps bond variables (kinks) to site variables (spins). We develop a numerical procedure for computing the asymmetry using tensor network methods and use it to demonstrate the different types of entanglement and entanglement asymmetry.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Accelerating Scientific Computing in the Post-Moore’s Era

Novel uses of graphical processing units for accelerated computation revolutionized the field of high-performance scientific computing by providing specialized workflows tailored to algorithmic requirements. As the era of Moore’s law draws to a close, many new non–von Neumann processors are emerging as potential computational accelerators, including those based on the principles of neuromorphic computing, tensor algebra, and quantum information. While development of these new processors is continuing to mature, the potential impact on accelerated computing is anticipated to be profound. We discuss how different processing models can advance computing in key scientific paradigms: machine learning and constraint satisfaction. Significantly, each of these new processor types utilizes a fundamentally different model of computation, and this raises questions about how to best use such processors in the design and implementation of applications. While many processors are being developed with a specific domain target, the ubiquity of spin-glass models and neural networks provides an avenue for multi-functional applications. Furthermore, this also hints at the infrastructure needed to integrate next-generation processing units into future high-performance computing systems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Entanglement features of random neural network quantum states

Restricted Boltzmann machines (RBMs) are a class of neural networks that have been successfully employed as a variational ansatz for quantum many-body wave functions. Here, we develop an analytic method to study quantum many-body spin states encoded by random RBMs with independent and identically distributed complex Gaussian weights. By mapping the computation of ensemble-averaged quantities to statistical mechanics models, we are able to investigate the parameter space of the RBM ensemble in the thermodynamic limit. We discover qualitatively distinct wave functions by varying RBM parameters, which correspond to distinct phases in the equivalent statistical mechanics model. Notably, there is a regime in which the typical RBM states have near-maximal entanglement entropy in the thermodynamic limit, similar to that of Haar-random states. However, these states generically exhibit nonergodic behavior in the Ising basis, and do not form quantum state designs, making them distinguishable from Haar-random states.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Lieb-Robinson Bounds with Exponential-in-Volume Tails

Lieb-Robinson bounds demonstrate the emergence of locality in many-body quantum systems. Intuitively, Lieb-Robinson bounds state that, with local or exponentially decaying interactions, the correlation that can be built up between two sites separated by distance 𝑟 after a time 𝑡 decays as exp (𝑣⁢𝑡 −𝑟), where 𝑣 is the emergent Lieb-Robinson velocity. In many problems, it is important to also capture how much of an operator grows to act on 𝑟 𝑑 sites in 𝑑 spatial dimensions. Perturbation theory and cluster expansion methods suggest that, at short times, these volume-filling operators are suppressed as exp (−𝑟 𝑑 ). We confirm this intuition, showing that, for 𝑟 >𝑣⁢𝑡, the volume-filling operator is suppressed by exp (−(𝑟−𝑣⁢𝑡) 𝑑 /(𝑣⁢𝑡) 𝑑−1 ). This closes a conceptual and practical gap between the cluster expansion and the Lieb-Robinson bound. We then present two very different applications of this new bound. Firstly, we obtain improved bounds on the classical computational resources necessary to simulate many-body dynamics with error tolerance 𝜀 for any finite time 𝑡: as 𝜀 becomes sufficiently small, only 𝜀 −O⁡(𝑡 𝑑−1 ) resources are needed. A protocol that likely saturates this bound is given. Secondly, we prove that disorder operators have volume-law suppression near the “solvable (Ising) point” in quantum phases with spontaneous symmetry breaking, which implies a new diagnostic for distinguishing many-body phases of quantum matter.

computational complexity↗
DOE OSTI2020

The matrix product approximation for the dynamic cavity method

Stochastic dynamics of classical degrees of freedom, defined on vertices of locally tree-like graphs, can be studied in the framework of the dynamic cavity method which is exact for tree graphs. Such models correspond for example to spin-glass systems, Boolean networks, neural networks, and other technical, biological, and social networks. The central objects in the cavity method are edge messages—conditional probabilities of two vertex variable trajectories. In this paper, we discuss a rather pedagogical derivation for the dynamic cavity method, give a detailed account of the novel matrix product edge message (MPEM) algorithm for the solution of the dynamic cavity equation as introduced, and present optimizations and extensions. Matrix product approximations of the edge messages are constructed recursively in an iteration over time. Computation costs and precision can be tuned by controlling the matrix dimensions of the MPEM in truncations. Without truncations, the dynamics is exact. Data for Glauber–Ising dynamics shows a linear growth of computation costs in time. In contrast to Monte Carlo simulations, the approach has a much better error scaling. Hence, it gives for example access to low probability events and decaying observables like temporal correlations. Here, we discuss optimized truncation schemes and an extension that allows to capture models which have a continuum time limit.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Superdiffusion resilience in Heisenberg chains with two-dimensional interactions on a quantum processor

Superdiffusive spin transport in the one-dimensional (1D) Heisenberg model is a key theoretical discovery in nonequilibrium quantum many-body physics. Although extensively studied in 1D systems, the breakdown and sustenance of superdiffusion in two-dimensional (2D) lattices with integrability-breaking terms, as found in real materials, remains an open question. To address this, we develop a toy model that extends the 1D Heisenberg model with a representative set of 2D interaction types and tunable strengths. Our model exhibits varying degrees of superdiffusion breakdown depending on the interaction type, spanning ballistic to diffusive regimes. We establish and justify a hierarchy of 2D interactions based on their resilience against superdiffusion breakdown: Heisenberg >𝑋⁢𝑋 > Ising. This precise control over the superdiffusive behavior also enables rigorous benchmarking of quantum hardware, and our simulations on IBM's Heron devices confirm the hardware's ability to accurately capture these many-body nonequilibrium phenomena. Overall, our results are relevant not only to simulating superdiffusion in real materials, such as the 1D Heisenberg compound KCuF3, which contains modest nonintegrable 2D terms, but also to extending superdiffusive behavior to larger 2D qubit lattices and other 2D materials.

Alagarsamy Manikandan, Keerthi Kumaran [ORNL]↗
DOE OSTI2025

Exact solution of the frustrated Potts model with next-nearest-neighbor interactions in one dimension via AI bootstrapping

The one-dimensional (1D) 𝐽 1 −𝐽 2 𝑞-state Potts model is solved exactly for arbitrary 𝑞 by analytically block-diagonalizing the original 𝑞 2 ×𝑞 2 transfer matrix into a simple 2 × 2 maximally symmetric subspace, based on using OpenAI's reasoning model o3-mini-high to exactly solve the 𝑞 = 3 case. Furthermore, by matching relevant subspaces, we map the Potts model onto a simpler effective 1D 𝑞-state Potts model, where 𝐽 2 acts as the nearest-neighbor interaction and 𝐽 1 as an effective magnetic field, nontrivially generalizing a 56-year-old theorem previously limited to the simplest case (𝑞 = 2, the Ising model). Our exact results provide insights to phenomena such as atomic or electronic order stacking in layered materials and the emergence of dome-shaped phases in complex phase diagrams. In conclusion, this work is anticipated to fuel both research in 1D frustrated magnets for recently discovered finite-temperature application potentials and the fast moving topic area of AI in science.

1-dimensional spin chains↗
DOE OSTI2025

Frustrated charge density wave and quasi-long-range bond-orientational order in the magnetic kagome FeGe

The intrinsic frustrated nature of a kagome lattice is amenable to the realization of exotic phases of matter, such as quantum spin liquids or spin ices, and the multiple-q charge density waves (CDW) in the kagome metals. Despite intense efforts to understand the mechanism driving the electronic modulations, its origin is still unknown and obscured by competing interactions and intertwined orders. Here, we identify a dimerization-driven 2D hexagonal charge-diffuse precursor in the antiferromagnetic kagome metal FeGe and demonstrate that the fraction of dimerized/undimerized states is the relevant order parameter of the multiple-q CDW of a continuous phase transition. The pretransitional charge fluctuations with propagation vector q = q M at T CDW < T < T*(125 K) are anisotropic, hence holding a quasi-long-range bond-orientational order. The broken translational symmetry emerges from the anisotropic diffuse precursor, akin to the Ising scenario of antiferromagnetic triangular lattices. The temperature and momentum dependence of the critical scattering show parallels to the stacked hexatic B-phases reported in liquid crystals and transient states of CDWs and highlight the key role of the topological defect-mediated melting of the CDW in FeGe.

Subires, D. [Donostia International Physics Center↗
DOE OSTI2021

Beyond magnons in Nd 2 ScNbO 7 : An Ising pyrochlore antiferromagnet with all-in–all-out order and random fields

We report the low-temperature magnetic properties of Nd 3+ pyrochlore Nd 2 ScNbO 7 . Susceptibility and magnetization show an easy-axis moment, and heat capacity reveals a phase transition to long-range order at T N =371(2) mK with a fully recovered ΔS=Rln(2), 53% of it recovered for T>T N . Elastic neutron scattering shows a long-range all-in all-out magnetic order with low-Q diffuse elastic scattering. Inelastic neutron scattering shows a low-energy flat band, indicating a magnetic Hamiltonian similar to Nd 2 Zr 2 O 7 . Nuclear hyperfine excitations measured by ultra-high-resolution neutron backscattering indicate a distribution of static electronic moments below T N , which may be due to B-site disorder influencing Nd crystal electric fields. Analysis of heat-capacity data shows an unexpected T-linear or T 3/2 term which is inconsistent with conventional magnon quasiparticles, but is consistent with fractionalized spinons or gapless local spin excitations. We use legacy data to show similar behavior in Nd 2 Zr 2 O 7 . Comparing local static moments also reveals a suppression of the nuclear Schottky anomaly in temperature, evidencing a fraction of Nd sites with nearly zero static moment, consistent with exchange-disorder-induced random singlet formation. Finally, taken together, these measurements suggest an unusual fluctuating magnetic ground state which mimics a spin liquid, but may not actually be one.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Single-ion properties of the transverse-field Ising model material CoNb 2 O 6

We report CoNb 2 O 6 is one of the few materials that is known to approximate the one-dimensional transverse-field Ising model (1D-TFIM) near its quantum critical point. It has been inferred that Co 2+ acts as a pseudospin 1/2 with anisotropic exchange interactions that are largely Ising-like, enabling the realization of the TFIM. However, the behavior of CoNb 2 O 6 is known to diverge from the ideal TFIM under transverse magnetic fields that are far from the quantum critical point, requiring the consideration of additional anisotropic, bond-dependent (Kitaev-like) terms in the microscopic pseudospin-1/2 Hamiltonian. These terms are expected to be controlled in part by single-ion physics, namely the wave function for the pseudospin-1/2 angular momentum doublet. Here, we present the results of both inelastic neutron scattering measurements and electron paramagnetic resonance spectroscopy on CoNb 2 O 6 , which elucidate the single-ion physics of Co 2+ in CoNb 2 O 6 for the first time. We find that the system is well-described by an intermediate spin-orbit coupled Hamiltonian, and the ground state is a well-isolated Kramers doublet with an anisotropic g tensor. We provide the approximate wave functions for this doublet, which we expect will be useful in theoretical investigations of the anisotropic exchange interactions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Tuning the Magnetic Properties of CrI 3 Using Ni Thin Film Deposition for Applications in Spintronic Devices

Chromium triiodide (CrI 3 ), a van der Waals magnet, has recently been shown to host Ising ferromagnetism down to the monolayer limit. It is a potentially important material in 2D magnet-based applications, such as magnetic sensors and spintronic devices. Prior studies have revealed the coexistence of two different types of interlayer magnetic coupling, the antiferromagnetic (AFM) coupling near the surface and ferromagnetic (FM) coupling in the deep bulk layers, in pristine CrI 3 crystals below the Curie temperature. In this study, we used Ni thin film deposition to tune the surface magnetic states in bulk CrI 3 . A nanometer thickness (4 nm thick) of Ni was deposited on the surface of CrI 3 using electron-beam evaporation to form a Ni/CrI 3 heterostructure. The deposited Ni thin layer forms nanoclusters that completely cover the CrI 3 surface. Magnetic states of CrI 3 before and after Ni deposition are probed by ultralow-frequency magneto-Raman spectroscopy. Instead of seeing three spin wave branches in Raman scattering below 2 T, as in pristine bulk CrI 3 , in Ni/CrI 3 , we observe only a single spin wave branch, which is softened compared to that in pristine CrI 3 and displays a Zeeman shift under an out-of-plane magnetic field up to 7 T. This observation reveals that the AFM layers on pristine CrI 3 crystal surfaces transit into FM layers, so the entire CrI 3 crystal is in the FM state after Ni deposition. First-principles calculations show that Ni atoms tend to diffuse into the CrI 3 lattice, and the FM interlayer coupling has a much lower energy than AFM coupling in the presence of Ni atoms. As a result, our studies show that the magnetic state of CrI 3 can be modified through deposition of a thin metal layer on the surface, offering a route for controlling the spin degree of freedom in van der Waals magnets.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Layer pseudospin magnetism in a transition metal dichalcogenide double-moiré system

Spontaneous order of layer pseudospins in two-dimensional bilayers is common in quantum Hall systems, where it is responsible for hysteretic responses to gate fields in states with Ising order and giant drag voltages in states with XY (spontaneous interlayer phase coherence) order. In this article we predict that layer pseudospin order will also occur in double-moiré strongly correlated two-dimensional electron systems. Here, we comment on similarities and differences in the competition between the two types of order in quantum Hall and double-moiré systems, and relate our findings to previous work on Falicov-Kimball models of electronic ferroelectrics.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Topological phase diagram of generalized Su-Schrieffer-Heeger models with interactions

We investigate interacting Su-Schrieffer-Heeger (SSH) chains with two- and three-site unit cells using density matrix renormalization group (DMRG) simulations. By selecting appropriate filling fractions and sweeping across interaction strength 𝐽 𝑧 and dimerization 𝛿, we map out their phase diagrams and identify transition lines via entanglement entropy and magnetization measurements. In the two-site model, we observe the emergence of an interaction-induced antiferromagnetic intermediate phase between the topologically trivial and nontrivial regimes, as well as a critical region at negative 𝐽 𝑧 with suppressed magnetization and finite-size scaling of entanglement entropy. In contrast, the three-site model lacks an intermediate phase and exhibits asymmetric edge localization and antiferromagnetic ordering in both positive and negative 𝐽 𝑧 regimes. We further examine the response of edge states to Ising perturbations. In the two-site model, zero-energy edge modes are topologically protected and remain robust up to a finite interaction strength. However, in the three-site model, where the edge states reside at finite energy, this protection breaks down. Despite this, the edge-localized nature of these states survives in the form of polarized modes whose spatial profiles reflect the noninteracting limit.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗