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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Generative Modeling for Machine Learning on the D-Wave

These are slides on Generative Modeling for Machine Learning on the D-Wave. The following topics are detailed: generative models; Boltzmann machines: a generative model; restricted Boltzmann machines; learning parameters: RBM training; practical ways to train RBM; D-Wave as a Boltzmann sampler; mapping RBM onto the D-Wave; Chimera restricted RBM; mapping binary RBM to Ising model; experiments; data; D-Wave effective temperature, parameters noise, etc.; experiments: contrastive divergence (CD) 1 step; after 50 steps of CD; after 100 steps of CD; D-Wave (experiments 1, 2, 3); D-Wave observations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Solving the 2D SUSY Gross-Neveu-Yukawa model with conformal truncation

We use Lightcone Conformal Truncation to analyze the RG flow of the two-dimensional supersymmetric Gross-Neveu-Yukawa theory, i.e. the theory of a real scalar superfield with a $\mathbb{Z}$ 2 -symmetric cubic superpotential, aka the 2d Wess-Zumino model. The theory depends on a single dimensionless coupling $\overline{g}$, and is expected to have a critical point at a tuned value ${\overline{g}}_{\ast }$ where it flows in the IR to the Tricritical Ising Model (TIM); the theory spontaneously breaks the $\mathbb{Z}$ 2 symmetry on one side of this phase transition, and breaks SUSY on the other side. We calculate the spectrum of energies as a function of $\overline{g}$ and see the gap close as the critical point is approached, and numerically read off the critical exponent ν in TIM. Beyond the critical point, the gap remains nearly zero, in agreement with the expectation of a massless Goldstino. We also study spectral functions of local operators on both sides of the phase transition and compare to analytic predictions where possible. In particular, we use the Zamolodchikov C -function to map the entire phase diagram of the theory. Crucial to this analysis is the fact that our truncation is able to preserve supersymmetry sufficiently to avoid any additional fine tuning.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1991

Ordering and phase transitions in random-field Ising systems

An exact analysis of the Ising model with infinite-range interactions in a random field and a local mean-field theory in three dimensions is carried out leading to a phase diagram with several coexistence surfaces and lines of critical points. The results show that the phase diagram depends crucially on whether the distribution of random fields is symmetric or not. Thus, Ising-like phase transitions in a porous medium (the asymmetric case) are in a different universality class from the conventional random-field model (symmetric case).

Maritan, Amos↗
DOE OSTI2015

Consistency tests of classical and quantum models for a quantum annealer

Lately the question of whether the D-Wave processors exhibit large-scale quantum behavior or can be described by a classical model has attracted significant interest. In this work we address this question by studying a 503 qubit D-Wave Two device in the "black box" model, i.e., by studying its input-output behavior. Our work generalizes an approach introduced in Boixo et al. [Nat. Commun. 4, 2067 (2013)], and uses groups of up to 20 qubits to realize a transverse Ising model evolution with a ground state degeneracy whose distribution acts as a sensitive probe that distinguishes classical and quantum models for the D-Wave device. Our findings rule out all classical models proposed to date for the device and provide evidence that an open system quantum dynamical description of the device that starts from a quantized energy level structure is well justified, even in the presence of relevant thermal excitations and a small value of the ratio of the single-qubit decoherence time to the annealing time.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
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 OSTI2024

Emergence of the isotropic Kitaev honeycomb lattice $α$–RuCl 3 and its magnetic properties

Here, we present a comprehensive investigation of the crystal and magnetic structures of the van der Waals antiferromagnet $α$–RuCl 3 using single crystal x-ray and neutron diffraction. The crystal structure at room temperature is a monoclinic (C2/m). However, with decreasing temperature, a remarkable first-order structural phase transition is observed, leading to the emergence of a rhombohedral ($R\bar3$) structure characterized by three-fold rotational symmetry forming an isotropic honeycomb lattice. On further cooling, a zigzag-type antiferromagnetic order develops below $T_{\textrm{N}} = 6\sim6.6$ K. The critical exponent of the magnetic order parameter was determined to be $\beta = 0.11(1)$, which is close to the two-dimensional Ising model. Additionally, the angular dependence of the magnetic critical field of the zigzag antiferromagnetic order for the polarized ferromagnetic phase reveals a six-fold rotational symmetry within the ab–plane. These findings reflect the symmetry associated with the Ising-like bond-dependent Kitaev spin interactions and underscore the universality of the Kitaev interaction-dominated antiferromagnetic system.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

Emergent order in the kagome Ising magnet Dy 3 Mg 2 Sb 3 O 14

The Ising model—in which degrees of freedom (spins) are binary valued (up/down)—is a cornerstone of statistical physics that shows rich behaviour when spins occupy a highly frustrated lattice such as kagome. Here we show that the layered Ising magnet Dy 3 Mg 2 Sb 3 O 14 hosts an emergent order predicted theoretically for individual kagome layers of in-plane Ising spins. Neutron-scattering and bulk thermomagnetic measurements reveal a phase transition at ~ 0.3 K from a disordered spin-ice-like regime to an emergent charge ordered state, in which emergent magnetic charge degrees of freedom exhibit three-dimensional order while spins remain partially disordered. Monte Carlo simulations show that an interplay of inter-layer interactions, spin canting and chemical disorder stabilizes this state. Our results establish Dy 3 Mg 2 Sb 3 O 14 as a tuneable system to study interacting emergent charges arising from kagome Ising frustration.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Adaptive variational quantum minimally entangled typical thermal states for finite temperature simulations

Scalable quantum algorithms for the simulation of quantum many-body systems in thermal equilibrium are important for predicting properties of quantum matter at finite temperatures. Here we describe and benchmark a quantum computing version of the minimally entangled typical thermal states (METTS) algorithm for which we adopt an adaptive variational approach to perform the required quantum imaginary time evolution. The algorithm, which we name AVQMETTS, dynamically generates compact and problem-specific quantum circuits, which are suitable for noisy intermediate-scale quantum (NISQ) hardware. We benchmark AVQMETTS on statevector simulators and perform thermal energy calculations of integrable and nonintegrable quantum spin models in one and two dimensions and demonstrate an approximately linear system-size scaling of the circuit complexity. We further map out the finite-temperature phase transition line of the two-dimensional transverse field Ising model. Finally, we study the impact of noise on AVQMETTS calculations using a phenomenological noise model.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Magnetic dynamics in NiTiO 3 honeycomb antiferromagnet using neutron scattering

The ilmenite NiTiO 3 consists of a buckled honeycomb lattice, with the Ni spins aligned ferromagnetically in-plane and antiferromagnetically out-of-plane. Using neutron spectroscopy, the magnetic structure and the dynamics were investigated as a function of temperature. Dispersive acoustic bands and nearly dispersionless optical bands at ≈ 3.7 meV are described by a highly anisotropy Heisenberg model with stronger antiferromagnetic (AFM) out-of-plane, weaker ferromagnetic (FM) in-plane interactions and an anisotropy gap of 0.95 meV. Furthermore, the order parameter yields a critical exponent between the Heisenberg and two-dimensional Ising models, consistent with highly anisotropic Heisenberg systems. The frustration parameter ≈ 2 supports a weakly frustrated system.

Antiferromagnets↗
DOE OSTI2021

Dynamics of a fractal set of first-order magnetic phase transitions in frustrated Lu 2 CoMnO 6

The axial next-nearest-neighbor Ising model predicts a fractal (infinite) set of phases with incommensurate wave vectors that are separated by first-order phase boundaries. This complexity results from a simple frustration condition between nearest- and next-nearest-neighbor interactions along a chain of Ising spins. Using x-ray photon correlation spectroscopy (XPCS), we investigate the surprising antiferromagnetic dynamics that emerge from such a complex phase diagram over a wide range of temperatures. In this work, we present XPCS measurements of the frustrated magnetic chain compound Lu 2 CoMnO 6 and Monte Carlo simulations. Incommensurate magnetic Bragg peaks slide towards commensurate “up-up-down-down” spin order with decreasing temperature and increasing time. Both simulation and experiment support a counterintuitive “upside-down” temperature dependence of the magnetic dynamics: at higher temperatures in the region of first-order phase boundaries, slower dynamics are observed where the speckle maintains its coherence. At the lowest temperatures, where part of the sample adopts commensurate order, the dynamics speed up and result in fast decoherence.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Experimental error mitigation using linear rescaling for variational quantum eigensolving with up to 20 qubits

Quantum computers have the potential to help solve a range of physics and chemistry problems, but noise in quantum hardware currently limits our ability to obtain accurate results from the execution of quantum-simulation algorithms. Various methods have been proposed to mitigate the impact of noise on variational algorithms, including several that model the noise as damping expectation values of observables. In this work, we benchmark various methods, including a new method proposed here. Here, we compare their performance in estimating the ground-state energies of several instances of the 1D mixed-field Ising model using the variational-quantum-eigensolver algorithm with up to 20 qubits on two of IBM's quantum computers. We find that several error-mitigation techniques allow us to recover energies to within 10% of the true values for circuits containing up to about 25 ansatz layers, where each layer consists of CNOT gates between all neighboring qubits and Y-rotations on all qubits.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Thermal tensor renormalization group simulations of square-lattice quantum spin models

In this work, we benchmark the well-controlled and numerically accurate exponential thermal tensor renormalization group (XTRG) in the simulation of interacting spin models in two dimensions. Finite temperature introduces a finite thermal correlation length ξ, such that for system sizes L >> ξ finite-size calculations actually simulate the thermodynamic limit. Here in this paper, we focus on the square lattice Heisenberg antiferromagnet (SLH) and quantum Ising models (QIM) on open and cylindrical geometries up to width W = 10 . We explore various one-dimensional mapping paths in the matrix product operator (MPO) representation, whose performance is clearly shown to be geometry dependent. We benchmark against quantum Monte Carlo (QMC) data, yet also the series-expansion thermal tensor network results. Thermal properties including the internal energy, specific heat, and spin structure factors, etc. are computed with high precision, obtaining excellent agreement with QMC results. XTRG also allows us to reach remarkably low temperatures. For SLH, we obtain an energy per site u$^*_g$ ≃ -0.6694 (4) and a spontaneous magnetization m$^*_S$ ≃ 0.30 ( 1 ) already consistent with the ground-state properties, which is obtained from extrapolated low-T thermal data on W ≤ 8 cylinders and W ≤ 10 open strips, respectively. We extract an exponential divergence versus T of the structure factor S (M), as well as the correlation length ξ, at the ordering wave vector M = ($π,π$) , which represents the renormalized classical behavior and can be observed over a narrow but appreciable temperature window, by analyzing the finite-size data by XTRG simulations. For the QIM with a finite-temperature phase transition, we employ several thermal quantities, including the specific heat, Binder ratio, as well as the MPO entanglement to determine the critical temperature $T_c$.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Universal topological quantum computation with strongly correlated Majorana edge modes

Abstract Majorana-based quantum gates are not complete for performing universal topological quantum computation while Fibonacci-based gates are difficult to be realized electronically and hardly coincide with the conventional quantum circuit models. In reference Hu and Kane (2018 Phys. Rev. Lett. 120 066801), it has been shown that a strongly correlated Majorana edge mode in a chiral topological superconductor can be decomposed into a Fibonacci anyon τ and a thermal operator anyon ɛ in the tricritical Ising model. The deconfinement of τ and ɛ via the interaction between the fermion modes yields the anyon collisions and gives the braiding of either τ or ɛ . With these braidings, the complete members of a set of universal gates, the Pauli gates, the Hadamard gate and extra phase gates for one-qubit as well as controlled-NOT (CNOT) gate for two-qubits, are topologically assembled. Encoding quantum information and reading out the computation results can be carried out through electric signals. With the sparse-dense mixed encodings, we set up the quantum circuit where the CNOT gate turns out to be a probabilistic gate and design the corresponding devices with thin films of the chiral topological superconductor. As an example of the universal topological quantum computing, we show the application to Shor’s integer factorization algorithm.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Extending a scaling equation of state to QCD

Whether quantum chromodynamics (QCD) exhibits a phase transition at finite temperature and density is an open question. It is important for hydrodynamic modeling of heavy ion collisions and neutron-star mergers. Lattice QCD simulations have definitively shown that the transition from hadrons to quarks and gluons is a crossover when the baryon chemical potential is zero or small. Here, we combine the parametric scaling equation of state, usually associated with the three-dimensional Ising model, with a background equation of state based on a smooth crossover from hadrons to quarks and gluons. Comparison to experimental data from the Beam Energy Scan II at the Relativistic Heavy Ion Collider or in heavy ion experiments at other accelerators may allow the critical exponents and amplitudes in the scaling equation of state to be determined for QCD if a critical point exists.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2019

Memristive networks: From graph theory to statistical physics

In this work, we provide an introduction to a very specific toy model of memristive networks, for which an exact differential equation for the internal memory which contains the Kirchhoff laws is known. In particular, we highlight how the circuit topology enters the dynamics via an analysis of directed graph. We try to highlight in particular the connection between the asymptotic states of memristors and the Ising model, and the relation to the dynamics and statics of disordered systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Comparing three generations of D-Wave quantum annealers for minor embedded combinatorial optimization problems

Abstract Quantum annealing (QA) is a novel type of analog computation that aims to use quantum mechanical fluctuations to search for optimal solutions of Ising problems. QA in the transverse Ising model, implemented on D-Wave quantum processing units, are available as cloud computing resources. In this study we report concise benchmarks across three generations of D-Wave quantum annealers, consisting of four different devices, for the NP-hard discrete combinatorial optimization problems unweighted maximum clique and unweighted maximum cut on random graphs. The Ising, or equivalently quadratic unconstrained binary optimization, formulation of these problems do not require auxiliary variables for order reduction, and their overall structure and weights are not highly variable, which makes these problems simple test cases to understand the sampling capability of current D-Wave quantum annealers. All-to-all minor embeddings of size 52, with relatively uniform chain lengths, are used for a direct comparison across the Chimera, Pegasus, and Zephyr device topologies. A grid-search over annealing times and the minor embedding chain strengths is performed in order to determine the level of reasonable performance for each device and problem type. Experiment metrics that are reported are approximation ratios for non-broken chain samples, chain break proportions, and time-to-solution for the maximum clique problem instances. How fairly the quantum annealers sample optimal maximum cliques, for instances which contain multiple maximum cliques, is quantified using entropy of the measured ground state distributions. The newest generation of quantum annealing hardware, which has a Zephyr hardware connectivity, performed the best overall with respect to approximation ratios and chain break frequencies.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Polynomial time algorithms for estimating spectra of adiabatic Hamiltonians

Much research regarding quantum adiabatic optimization has focused on stoquastic Hamiltonians with Hamming-symmetric potentials, such as the well-studied “spike” example. Due to the large amount of symmetry in these potentials such problems are readily open to analysis both analytically and computationally. However, more realistic potentials do not have such a high degree of symmetry and may have many local minima. Here we present a somewhat more realistic class of problems consisting of many individually Hamming-symmetric potential wells. For two or three such wells we demonstrate that such a problem can be solved exactly in time polynomial in the number of qubits and wells. For greater than three wells, we present a tight-binding approach with which to efficiently analyze the performance of such Hamiltonians in an adiabatic computation. Here, we provide several basic examples designed to highlight the usefulness of this toy model and to give insight into using the tight-binding approach to examining it, including (1) an adiabatic unstructured search with a transverse field driver and a prior guess to the marked item and (2) a scheme for adiabatically simulating the ground states of small collections of strongly interacting spins, with an explicit demonstration for an Ising-model Hamiltonian.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Formulation of Multicomponent Lattice Gas Model Cluster Expansions Parameterized on Ab Initio Data: An Introduction to the Ab Initio Mean-Field Augmented Lattice Gas Modeling Code

A formalism is presented for the construction of multicomponent lattice gas models parameterized with ab initio (typically density functional theory) data. The leave-multiple-out and the leave-one-out cross validation score are showcased, and practical algorithms are developed and implemented in a new code called the Ab initio Mean-field Augmented Lattice Gas Modeling (AMALGM) code. Here, the functionality of these algorithms is demonstrated with a fully worked out example using the O/Fe(100) system. AMALGM and the formalism on which it is based are envisioned as a surface-oriented lattice gas alternative to other cluster expansion codes that are typically geared toward bulk systems whose lateral interactions between components are parameterized using an Ising model. While the formalism is created within the context of surfaces, it is equally applicable to bulk materials.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗