Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Ising 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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Coarsening dynamics of Ising-nematic order in a frustrated Heisenberg antiferromagnet

We study the phase ordering dynamics of the classical antiferromagnetic 𝐽 1 −𝐽 2 (nearest-neighbor and next-nearest-neighbor couplings) Heisenberg model on the square lattice in the strong frustration regime (𝐽 2 /𝐽 1 > 1/2). While thermal fluctuations preclude any long-range magnetic order at finite temperatures, the system exhibits a long-range spin-driven nematic phase at low temperatures. The transition into the nematic phase is further shown to belong to the two-dimensional Ising universality class based on the critical exponents near the phase transition. Our large-scale stochastic Landau-Lifshitz-Gilbert simulations find a two-stage phase ordering when the system is quenched from a high-temperature paramagnetic state into the nematic phase. In the early stage, collinear alignments of spins lead to a locally saturated Ising-nematic order. Once domains of well-defined Ising order are developed, the late-stage relaxation is dominated by curvature-driven domain coarsening, as described by the Allen-Cahn equation. The characteristic size of Ising-nematic domains scales as the square root of time, similar to the kinetic Ising model described by the time-dependent Ginzburg-Landau theory. Our results confirm that the late-stage ordering kinetics of the spin-driven nematic, which is a vestigial order of the frustrated Heisenberg model, belongs to the dynamical universality class of a nonconserved Ising order. Interestingly, the system shows no violation of the superuniversality hypothesis under weak bond disorder. The dynamic scaling invariance is preserved in the presence of weak bond disorder. Here, we also discuss possible applications of our results to materials for which vestigial Ising-nematic order is realized.

Antiferromagnets↗
DOE OSTI2022

Density Functional Theory Study of Controllable Optical Absorptions and Magneto-Optical Properties of Magnetic CrI 3 Nanoribbons: Implications for Compact 2D Magnetic Devices

A chromium triiodide (CrI 3 ) monolayer has an interesting ferromagnetic ground-state. In this work, we calculate band structures and magnetic moments of tensile-strained and bent zigzag CrI 3 nanoribbons with density functional theory. The edge iodine atoms form flat low-lying conduction bands and couple with chromium atoms ferromagnetically, while the non-edge iodine atoms weakly couple antiferromagnetically. Narrow CrI 3 nanoribbons have two locally stable magnetic moment orientations, namely out-of-plane and in-plane (along the nanoribbon periodic direction) configurations. This enables four magnetization states in CrI 3 nanoribbons, including two out-of-plane ones (up and down) and two in-plane ones (forward and backward along the nanoribbon periodical direction), increasing the operating controllability. Based on the one-dimensional Ising spin chain model, the spin correlation length of the narrow CrI 3 nanoribbon is estimated as about 10 Å at its estimated Curie temperature of 27 Kelvin, which is lower than the measured 45 Kelvin of the monolayer CrI 3 . The optical absorption and magneto-optical properties of CrI 3 nanoribbons are investigated with many-body perturbation GW-BSE (Bethe-Salpeter equation), including magnetic dichroism, Faraday and magneto-optical Kerr effects. Here, the low-energy dark excitons are mainly from transitions between electrons and holes with unlike spins and are non-Frenkel-like, while the bright excitons have mixed spin configurations. The intrinsic lifetime of excitons can be over one nanosecond, suitable for quantum information processes. Tensile strains and bending manifestly modulate the absorption spectra and magneto-optical properties of CrI 3 nanoribbons within a technologically important photon energy-range of ~1.0-2.0 eV. The CrI 3 nanoribbons can be used in 1D or 2D magnetic storage nanodevices, tunable magnetic optoelectronics, and spin-based quantum information controls.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Variational fast forwarding for quantum simulation beyond the coherence time

Abstract Trotterization-based, iterative approaches to quantum simulation (QS) are restricted to simulation times less than the coherence time of the quantum computer (QC), which limits their utility in the near term. Here, we present a hybrid quantum-classical algorithm, called variational fast forwarding (VFF), for decreasing the quantum circuit depth of QSs. VFF seeks an approximate diagonalization of a short-time simulation to enable longer-time simulations using a constant number of gates. Our error analysis provides two results: (1) the simulation error of VFF scales at worst linearly in the fast-forwarded simulation time, and (2) our cost function’s operational meaning as an upper bound on average-case simulation error provides a natural termination condition for VFF. We implement VFF for the Hubbard, Ising, and Heisenberg models on a simulator. In addition, we implement VFF on Rigetti’s QC to demonstrate simulation beyond the coherence time. Finally, we show how to estimate energy eigenvalues using VFF.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Single-ancilla ground state preparation via Lindbladians

We design a quantum algorithm for ground state preparation in the early fault tolerant regime. As a Monte Carlo style quantum algorithm, our method features a Lindbladian where the target state is stationary. The construction of this Lindbladian is algorithmic and should not be seen as a specific approximation to some weakly coupled system-bath dynamics in nature. Our algorithm can be implemented using just one ancilla qubit and efficiently simulated on a quantum computer. It can prepare the ground state even when the initial state has zero overlap with the ground state, bypassing the most significant limitation of methods like quantum phase estimation. As a variant, we also propose a discrete-time algorithm, demonstrating even better efficiency and providing a near-optimal simulation cost depending on the desired evolution time and precision. Numerical simulations using Ising and Hubbard models demonstrate the efficacy and applicability of our method. Published by the American Physical Society 2024

Ding, Zhiyan (ORCID:000000018863403X)↗
DOE OSTI2023

Network construction of non-Abelian chiral spin liquids

We use a network of chiral junctions to construct a family of topological chiral spin liquids in two spatial dimensions. The chiral spin liquid phase harbors SU(2) _k k anyons, which stem from the underlying SU(2) _k k WZW models that describe the constituent spin chains of the network. The network exhibits quantized spin and thermal Hall conductances. We illustrate our construction by inspecting the topological properties of the SU(2) _2 2 model. We find that this model has emergent Ising anyons, with spinons acting as vortex excitations that bind Majorana zero modes. We also show that the ground state of this network is threefold degenerate on the torus, asserting its non-Abelian character. Our results shed new light on the stability of non-Abelian topological phases in artificial quantum materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Kibble–Zurek mechanism of Ising domains

The formation of topological defects after a symmetry-breaking phase transition is an overarching phenomenon that encodes the underlying dynamics. The Kibble-Zurek mechanism (KZM) describes these non-equilibrium dynamics of second-order phase transitions and predicts a power-law relationship between the cooling rates and the density of topological defects. It has been verified as a successful model in a wide variety of physical systems, including structure formation in the early Universe and condensed-matter materials. However, it is uncertain if the KZM mechanism is also valid for topologically trivial Ising domains, one of the most common and fundamental types of domain in condensed-matter systems. Here we show that the cooling rate dependence of Ising domain density follows the KZM power law in two different three-dimensional structural Ising domains: ferro-rotation domains in NiTiO 3 and polar domains in BiTeI. However, although the KZM slope of NiTiO 3 agrees with the prediction of the 3D Ising model, the KZM slope of BiTeI exceeds the theoretical limit, providing an example of steepening KZM slope with long-range dipolar interactions. Finally, our results demonstrate the validity of KZM for Ising domains and reveal an enhancement of the power-law exponent for transitions of non-topological quantities with long-range interactions. The Kibble-Zurek mechanism is shown to apply to structural Ising domains in three-dimensional materials. Long-range interactions modify the critical exponents away from theoretical predictions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Adaptive variational ground state preparation for spin-1 models on qubit-based architectures

Here, we apply the adaptive variational quantum imaginary time evolution (AVQITE) method to prepare ground states of one-dimensional spin S = 1 models. We compare different spin-to-qubit encodings (standard binary, Gray, unary, and multiplet) with regard to the performance and quantum resource cost of the algorithm. Using state-vector simulations, we study two well-known spin-1 models: the Blume-Capel model of transverse-field Ising spins with single-ion anisotropy, and the XXZ model with single-ion anisotropy. We consider system sizes of up to 20 qubits, which corresponds to spin-1 chains up to length 10. We determine the dependence of the number of CNOT gates in the AVQITE state preparation circuit on the encoding, the initial state, and the choice of operator pool in the adaptive method. Independent of the choice of encoding, we find that the CNOT gate count scales cubically with the number of spins for the Blume-Capel model and quartically for the anisotropic XXZ model. However, the multiplet and Gray encodings present smaller prefactors in the scaling relations. These results provide useful insights for the implementation of AVQITE on quantum hardware.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Unified theory of hastatic order and antiferromagnetism in URu 2 Si 2

The hidden-order phase of URu 2 Si 2 has eluded identification for over 35 years. A compelling proposal that explains the Ising heavy-fermion nature of the material is hastatic order: a symmetry-breaking heavy Fermi liquid arising from a spinorial microscopic hybridization. The original hastatic proposal cannot microscopically model the pressure-induced antiferromagnetic phase; and while it predicts a spinorial order parameter, it does not provide any detectable signatures of the spinorial nature. Here, we present a more realistic microscopic model of hastatic order in URu 2 Si 2 based on two conducting electron bands and explore its phase diagram in detail. Our model nontrivially preserves the Ising heavy-fermion physics of the original, while allowing us to tune between the antiferromagnet and hidden order using pressure analogs and magnetic field. Furthermore, our model is also consistent with recent phenomenological predictions of not one, but two, vector order parameters associated with the spinorial order that provide the first microscopic predictions for detecting the spinorial nature.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

NeuroSymbolic Approaches as a Vector for Assured Artificial Intelligence

The deployment of artificial intelligence systems in critical applications requires higher levels of assurance for safety, security, and interpretability. While neurosymbolic (NESY) approaches combining neural networks with symbolic reasoning offer potential advantages for assured AI, existing differentiable neurosymbolic frameworks face significant limitations including computational overhead and performance constraints. This report investigates the ISED (InferSampleEstimateDescend) framework as an alternative approach that enables neurosymbolic learning without requiring endtoend differentiability. We evaluate ISED’s utility for geointelligence applications by comparing neurosymbolic models against standard neural networks on aircraft classification tasks using the RarePlanes and MTARSI imagery datasets. Our results demonstrate that while ISEDbased models achieve slightly lower accuracy (89.7% vs 92.1% on RarePlanes; 91.1% vs 92.5% on MTARSI), they provide critical explainability capabilities that enable tracing incorrect predictions back to specific attribute misclassifications. We also present an automated pipeline that generates both attributeclass mappings and neurosymbolic model architectures from natural language descriptions, significantly reducing the manual effort required for NESY model deployment. These findings suggest that ISED offers a promising direction for developing assured AI systems where interpretability and reasoning transparency are prioritized alongside performance.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Dynamic scaling of order parameter fluctuations in model B

Brookhaven National Laboratory has performed a dedicated set of experiments, called the Beam Energy Scan (BES), with the goal of location a possible critical point in the phase diagram of nuclear matter, This point is analogous to the endpoint of the liquid-gas phase transition in water, where the liquid phase corresponds to a phase of hadrons, and the gas phase to a plasma of quarks and gluons. In order to describe the data generated by these experiments we have studied numerical simulations of the stochastic diffusion equation with a conserved charge. Here, we focus on the dynamics in the vicinity of a critical point in the Ising universality class. The model we consider is expected to describe the critical dynamics near a possible QCD critical point if the coupling of the order parameter to the momentum density of the fluid can be neglected. The simulations are performed on a spatial lattice, and the time evolution is performed using a Metropolis algorithm. We determine the dynamical critical exponent z ≃ 3.972(2), which agrees with predictions of the epsilon expansion. We also study non-equilibrium sweeps of the reduced temperature and observe approximate Kibble-Zurek scaling.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Surface magnetism in Fe 3 GeTe 2 van der Waals ferromagnet

The surface magnetization of Fe 3 GeTe 2 was examined by low-energy electron microscopy (LEEM) using an off-normal incidence electron beam. We found that the 180o domain walls are of Bloch type. Temperature-dependent LEEM measurements yield a surface magnetization with a surface critical exponent β1 = 0.79±0.02. This result is consistent with surface magnetism in the 3D semi-infinite Heisenberg (β1 = 0.84±0.01) or Ising (β1 = 0.78±0.02) models, which is distinctly different from the bulk exponent (β = 0.34±0.07). Further, the measurements reveal the power of LEEM with a tilted beam to determine magnetic domain structure in quantum materials, without need for the use of spin-polarized electrons. Single crystal diffraction measurements reveal inversion symmetry-breaking weak peaks and yield space group P-6m2. This Fe site defect-derived loss of inversion symmetry enables the formation of skyrmions in this Fe 3 GeTe 2 crystal.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Phase diagram of generalized XY model using the tensor renormalization group

We use the higher-order tensor renormalization group method to study the two-dimensional generalized XY model that admits integer and half-integer vortices. This model is the deformation of the classical XY model and has a rich phase structure consisting of nematic, ferromagnetic, and disordered phases and three transition lines belonging to the Berezinskii-Kosterlitz-Thouless and Ising class. We explore the model for a wide range of temperatures, T , and the deformation parameter, Δ , and compute specific heat along with integer and half-integer magnetic susceptibility, finding both Berezinskii-Kosterlitz-Thouless-like and Ising-like transitions and the region where they meet. Published by the American Physical Society 2024

2-dimensional systems↗
DOE OSTI2023

Spin dynamics of the generalized quantum spin compass chain

Here we calculate the dynamical spin structure factor of the generalized spin-1/2 compass spin chain using the density matrix renormalization group. The model, also known as the twisted Kitaev spin chain, was recently proposed to be relevant for the description of the spin chain compound CoNb 2 O 6 . It features bond-dependent interactions and interpolates between an Ising chain and a one-dimensional variant of Kitaev's honeycomb spin model. The structure factor, in turn, is found to interpolate from gapped and nondispersive in the Ising limit to gapless with nontrivial continua in the Kitaev limit. In particular, the component of the structure factor perpendicular to the Ising directions changes abruptly at the Kitaev point into a dispersionless continuum due to the emergence of an extensive ground-state degeneracy. We show this continuum is consistent with analytical Jordan-Wigner results. We also discuss implications for future inelastic scattering experiments and applications to materials, particularly CoNb 2 O 6 .

1-dimensional spin chains↗
DOE OSTI2024

Two-fermion negativity and confinement in the Schwinger model

We consider the fermionic (logarithmic) negativity between two fermionic modes in the Schwinger model. Recent results pointed out that fermionic systems can exhibit stronger entanglement than bosonic systems, exhibiting a negativity that decays only algebraically. The Schwinger model is described by fermionic excitations at short distances, while its asymptotic spectrum is the one of a bosonic theory. We show that the two-mode negativity detects this confining, fermion-to-boson transition, shifting from an algebraic decay to an exponential decay at distances of the order of the de Broglie wavelength of the first excited state. We derive analytical expressions in the massless Schwinger model and confront them with tensor network simulations. We also perform tensor network simulations in the massive model, which is not solvable analytically, and close to the Ising quantum critical point of the Schwinger model, where we show that the negativity behaves as its bosonic counterpart. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Quantum model learning agent: characterisation of quantum systems through machine learning

Accurate models of real quantum systems are important for investigating their behaviour, yet are difficult to distil empirically. Here, we report an algorithm—the quantum model learning agent (QMLA)—to reverse engineer Hamiltonian descriptions of a target system. We test the performance of QMLA on a number of simulated experiments, demonstrating several mechanisms for the design of candidate Hamiltonian models and simultaneously entertaining numerous hypotheses about the nature of the physical interactions governing the system under study. QMLA is shown to identify the true model in the majority of instances, when provided with limited a priori information, and control of the experimental setup. Our protocol can explore Ising, Heisenberg and Hubbard families of models in parallel, reliably identifying the family which best describes the system dynamics. We demonstrate QMLA operating on large model spaces by incorporating a genetic algorithm to formulate new hypothetical models. The selection of models whose features propagate to the next generation is based upon an objective function inspired by the Elo rating scheme, typically used to rate competitors in games such as chess and football. In all instances, our protocol finds models that exhibit F 1 score ≥ 0.88 when compared with the true model, and it precisely identifies the true model in 72% of cases, whilst exploring a space of over 250 000 potential models. By testing which interactions actually occur in the target system, QMLA is a viable tool for both the exploration of fundamental physics and the characterisation and calibration of quantum devices.

97 MATHEMATICS AND COMPUTING↗
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 OSTI2023

Emergent Halperin–Saslow mode, gauge glass and quenched disorders in quantum Ising magnet TmMgGaO 4

Quenched disorders could bring novel quantum excitations and models to certain quantum magnets. Motivated by recent experiments on the quantum Ising magnet TmMgGaO 4 , we explore the effects of the quenched disorder and the interlayer coupling in this triangular lattice Ising antiferromagnet. It is pointed out that the weak quenched (nonmagnetic) disorder would convert the emergent 2D Berezinskii–Kosterlitz–Thouless (BKT) phase and the critical region into a U(1) gauge glass. There will be an emergent Halperin–Saslow mode associated with this gauge glass. Using the Imry-Ma’s renormalization group result, we explain the fate of the finite-field [Formula: see text] symmetry breaking transition at the low temperatures. The ferromagnetic interlayer coupling would suppress the BKT phase and generate a tiny ferromagnetism. With quenched disorders, this interlayer coupling changes the 2D gauge glass into a 3D gauge glass, and the Halperin–Saslow mode persists. This work merely focuses on addressing a phase regime in terms of emergent U(1) gauge glass behaviors and hope to inspire future works and thoughts in weakly disordered frustrated magnets in general.

Physics↗
DOE OSTI2022

Universal quantum operation of spin-3/2 Blume-Capel chains

We propose a logical qubit based on the Blume-Capel model: A higher-spin generalization of the Ising chain which allows for an on-site anisotropy-preserving rotational invariance around the Ising axis. We show that such a spin-3/2 Blume-Capel model can also support localized Majorana zero modes at the ends of the chain. Inspired by known braiding protocols of these Majorana zero modes, upon appropriate manipulation of the system parameters, we demonstrate a set of universal gate operations which act on qubits encoded in the doubly degenerate ground states of the chain.

1-dimensional spin chains↗