One-dimensional Ising Model Research Articles

Non-Hermitian Ising model at finite temperature

As a very simple model, the Ising model plays an important role in statistical physics. In the paper, with the help of quantum Liouvillian statistical theory, we study the one-dimensional non-Hermitian Ising model at finite temperature and give its analytical solutions. We find that the non-Hermitian Ising model shows quite different properties from those of its Hermitian counterpart. For example, the ‘pseudo-phase transition’ is explored between the ‘topological’ phase and the ‘non-topological’ phase, at which the Liouvillian energy gap is closed rather than the usual energy gap. In particular, we point out that the one-dimensional non-Hermitian Ising model at finite temperature can be equivalent to an effective anisotropic XY model in the transverse field. This work will help people understand quantum statistical properties of non-Hermitian systems at finite temperatures.

Signatures of Quantum Criticality in the Complex Inverse Temperature Plane

Concepts of the complex partition functions and the Fisher zeros provide intrinsic statistical mechanisms for finite temperature and real time dynamical phase transitions. We extend the utility of these complexifications to quantum phase transitions. We exactly identify different Fisher zeros on lines or closed curves and elucidate their correspondence with domain-wall excitations or confined mesons for the one-dimensional transverse field Ising model. The crossover behavior of the Fisher zeros provides a fascinating picture for criticality near the quantum phase transition, where the excitation energy scales are quantitatively determined. We further confirm our results by tensor network calculations and demonstrate a clear signal of deconfined meson excitations from the disruption of the closed zero curves. Our results unambiguously show significant features of Fisher zeros for a quantum phase transition and open up a new route to explore quantum criticality.

Navigating the noise-depth tradeoff in adiabatic quantum circuits

Adiabatic quantum algorithms solve computational problems by slowly evolving a trivial state to the desired solution. On an ideal quantum computer, the solution quality improves monotonically with increasing circuit depth. By contrast, increasing the depth in current noisy computers introduces more noise and eventually deteriorates any computational advantage. What is the optimal circuit depth that provides the best solution? Here, we address this question by investigating an adiabatic circuit that interpolates between the paramagnetic and ferromagnetic ground states of the one-dimensional quantum Ising model. We characterize the quality of the final output by the density of defects $d$, as a function of the circuit depth $N$ and noise strength $\sigma$. We find that $d$ is well-described by the simple form $d_\mathrm{ideal}+d_\mathrm{noise}$, where the ideal case $d_\mathrm{ideal}\sim N^{-1/2}$ is controlled by the Kibble-Zurek mechanism, and the noise contribution scales as $d_\mathrm{noise}\sim N\sigma^2$. It follows that the optimal number of steps minimizing the number of defects goes as $\sim\sigma^{-4/3}$. We implement this algorithm on a noisy superconducting quantum processor and find that the dependence of the density of defects on the circuit depth follows the predicted non-monotonous behavior and agrees well with noisy simulations. Our work allows one to efficiently benchmark quantum devices and extract their effective noise strength $\sigma$.

Villain model with long-range couplings

The nearest-neighbor Villain, or periodic Gaussian, model is a useful tool to understand the physics of the topological defects of the two-dimensional nearest-neighbor XY model, as the two models share the same symmetries and are in the same universality class. The long-range counterpart of the two-dimensional XY has been recently shown to exhibit a non-trivial critical behavior, with a complex phase diagram including a range of values of the power-law exponent of the couplings decay, σ, in which there are a magnetized, a disordered and a critical phase [1]. Here we address the issue of whether the critical behavior of the two-dimensional XY model with long-range couplings can be described by the Villain counterpart of the model. After introducing a suitable generalization of the Villain model with long-range couplings, we derive a set of renormalization-group equations for the vortex-vortex potential, which differs from the one of the long-range XY model, signaling that the decoupling of spin-waves and topological defects is no longer justified in this regime. The main results are that for σ < 2 the two models no longer share the same universality class. Remarkably, within a large region of its the phase diagram, the Villain model is found to behave similarly to the one-dimensional Ising model with 1/r2 interactions.

Do additional target points speed up evolutionary algorithms?

Most runtime analyses of randomised search heuristics focus on the expected number of function evaluations to find a unique global optimum. We ask a fundamental question: if additional search points are declared optimal, or declared as desirable target points, do these additional optima speed up evolutionary algorithms? More formally, we analyse the expected hitting time of a target set OPT∪S where S is a set of non-optimal search points and OPT is the set of optima and compare it to the expected hitting time of OPT.We show that the answer to our question depends on the number and placement of search points in S. For all black-box algorithms and all fitness functions with polynomial expected optimisation times we show that, if additional optima are placed randomly, even an exponential number of optima has a negligible effect on the expected optimisation time. Considering Hamming balls around all global optima gives an easier target for some algorithms and functions and can shift the phase transition with respect to offspring population sizes in the (1,λ) EA on OneMax. However, for the one-dimensional Ising model the time to reach Hamming balls of radius (1/2−ε)n around optima does not reduce the asymptotic expected optimisation time in the worst case. Finally, on functions where search trajectories typically join in a single search point, turning one search point into an optimum drastically reduces the expected optimisation time.

Simulating the Hamiltonian of dimer atomic spin model of one-dimensional optical lattice on quantum computers

The one-dimensional Ising model with its connections to several physical concepts plays a vital role in comprehension of several principles, phenomena and numerical methods. The Hamiltonian of a coupled one-dimensional dissipative spin system in the presence of magnetic field can be obtained from the Ising model. We simulate the above Hamiltonian by designing a quantum circuit with precise gate measurement and execute with the IBMQ experience platform through different [Formula: see text] states with controlled energy separation where we can check quantum synchronization in a dissipative lattice system. Our result shows the relation between various entangled states, the relation between the different energy separation ([Formula: see text]) with the spin–spin coupling ([Formula: see text]) in the lattice, along with fidelity calculations for several iterations of the model used. We also estimate the ground and first excited energy states of Ising-Hamiltonian using VQE algorithm and investigate the lowest energy values varying the number of layers of ansatz.

Quantum critical points, lines, and surfaces

In this paper we promote the idea of quantum critical lines ({\em inter alia} surfaces) as opposed to points. A quantum critical line obtains when criticality at zero temperature is extended over a continuum in a one-dimensional line. We base our ideas on a simple but exactly solved model introduced by one of the authors involving a one-dimensional quantum transverse field Ising model with added 3-spin interaction. While many of the ideas are quite general, there are other aspects that are not. In particular, a line of criticality with continuously varying exponents is not captured. However, the exact solvability of the model gives us considerable confidence in our results. Although the pure system is analytically exactly solved, the disorder case requires numerical analyses based on exact computation of the correlation function in the Pfaffian representation. The disorder case leads to dynamic structure factor as a function of frequency and wave vector. We expect that the model is experientally realizable and perhaps many other similar models will be found to explore quantum critical lines.

Finite-Size Scaling Theory at a Self-Dual Quantum Critical Point

The nondivergence of the generalized Grüneisen ratio (GR) at a quantum critical point (QCP) has been proposed to be a universal thermodynamic signature of self-duality. We study how the Kramers–Wannier-type self-duality manifests itself in the finite-size scaling behavior of thermodynamic quantities in the quantum critical regime. While the self-duality cannot be realized as a unitary transformation in the total Hilbert space for the Hamiltonian with the periodic boundary condition, it can be implemented in certain symmetry sectors with proper boundary conditions. Therefore, the GR and the transverse magnetization of the one-dimensional transverse-field Ising model exhibit different finite-size scaling behaviors in different sectors. This implies that the numerical diagnosis of self-dual QCP requires identification of the proper symmetry sectors.

Detecting Quantum Critical Points of Correlated Systems by Quantum Convolutional Neural Network Using Data from Variational Quantum Eigensolver

Machine learning has been applied to a wide variety of models, from classical statistical mechanics to quantum strongly correlated systems, for classifying phase transitions. The recently proposed quantum convolutional neural network (QCNN) provides a new framework for using quantum circuits instead of classical neural networks as the backbone of classification methods. We present the results from training the QCNN by the wavefunctions of the variational quantum eigensolver for the one-dimensional transverse field Ising model (TFIM). We demonstrate that the QCNN identifies wavefunctions corresponding to the paramagnetic and ferromagnetic phases of the TFIM with reasonable accuracy. The QCNN can be trained to predict the corresponding ‘phase’ of wavefunctions around the putative quantum critical point even though it is trained by wavefunctions far away. The paper provides a basis for exploiting the QCNN to identify the quantum critical point.

Calculation of the Free Energy of the Ising Model on a Cayley Tree via the Self-Similarity Method

In this study, an interactive Ising model having the nearest and prolonged next-nearest neighbors defined on a Cayley tree is considered. Inspired by the results obtained for the one-dimensional Ising model, we will construct the partition function and then calculate the free energy of the Ising model having the prolonged next nearest and nearest neighbor interactions and external field on a two-order Cayley tree using the self-similarity of the semi-infinite Cayley tree. The phase transition problem for the Ising system is investigated under the given conditions.

Statistics of the number ofdefects after quantum annealing in a thermal environment.

We study the statistics of the kink number generated by quantum annealing ina one-dimensional transverse Ising model coupled to a bosonic thermal bath. Using the freezing ansatz for quantum annealing in the thermal environment, we show the relation between the ratio of the second to the first cumulant of the kink number distribution and the average kink density. The theoretical result is confirmed thoroughly by numerical simulation using the non-Markovian infinite time-evolving block decimation which we proposed recently. The simulation using D-Wave's quantum annealer is also discussed. This article is part of the theme issue 'Quantum annealing and computation: challenges and perspectives'.

Effect of Substrate on One-Dimensional Multiferroic Properties

We present a theoretical study of the substrate influence on the electrical and magnetic properties of a one-dimensional multiferroic. We used a one-dimensional axial next-nearest neighbor Ising model (1D ANNNI model). The effect of the substrate was modeled using the periodic Frenkel–Kontorova potential. It is shown that the periodic potential of the substrate reduces the polarization of the multiferroic at low temperatures. The substrate potential significantly affects the structural changes near the magnetic phase transition temperature.

Flocking in one dimension: Asters and reversals.

We study the one-dimensional active Ising model in which aligning particles undergo diffusion biased by the signs of their spins. The phase diagram obtained varying the density of particles, their hopping rate, and the temperature controlling the alignment shows a homogeneous disordered phase but no homogeneous ordered one, as well as two phases with localized dense structures. In the flocking phase, large ordered aggregates move ballistically and stochastically reverse their direction of motion. In what we termed the "aster" phase, dense immobile aggregates of opposite magnetization face each other, exchanging particles, without any net motion of the aggregates. Using a combination of numerical simulations and mean-field theory, we study the evolution of the shapes of the flocks, the statistics of their reversal times, and their coarsening dynamics. Solving exactly for the zero-temperature dynamics of an aster allows us to understand their coarsening, which shows extremal dynamics, while mean-field equationsaccount for their shape.

Exact expressions for the partition function of the one-dimensional Ising model in the fixed-M ensemble.

We obtain exact closed-form expressions for the partition function of the one-dimensional Ising model in the fixed-M ensemble, for three commonly used boundary conditions: periodic, antiperiodic, and Dirichlet. These expressions allow for the determination of fluctuation-induced forces in the canonical ensemble, which we term Helmholtz forces. The thermodynamic expressions and the calculations flowing from them should provide insights into the nature and behavior of fluctuation-induced forces in interesting and as-yet unexplored regimes.

Multipartite nonlocality in an Ising model with a tilted magnetic field at zero and finite temperatures

We use the multipartite nonlocality to detect the quantum phase transitions (QPTs) in the one-dimensional Ising model with a tilted magnetic field. For a given transverse magnetic field, we investigate the impact of the coupling strength, the longitudinal magnetic field, and the temperature on the multipartite nonlocality and its first derivative. We show that both the multipartite nonlocality and its first derivative of the ground-state can perfectly characterize the QPTs of the model. We further show that the multipartite nonlocality at finite temperature exhibits different behavior from that of the ground-state. We find that the thermal fluctuations can weaken the multipartite nonlocality. However, a high longitudinal magnetic field can destroy the weakening effect of the thermal fluctuations, particularly near the QPTs point. In particular, we find that the logarithm measure of multipartite nonlocality both at zero and finite temperature shows a linear behavior for any fixed longitudinal magnetic field.

Quantum criticality using a superconducting quantum processor

Quantum criticality emerges from the collective behavior of many interacting quantum particles, often at the transition between different phases of matter. It is one of the cornerstones of condensed matter physics, which we access on noisy intermediate-scale (NISQ) quantum devices by leveraging a dynamically driven phenomenon. We probe the critical properties of the one-dimensional quantum Ising model on a programmable superconducting quantum chip via a Kibble-Zurek process, obtain scaling laws, and estimate critical exponents despite inherent sources of errors on the hardware. In addition, we investigate how the improvement of NISQ computers (more qubits, less noise) will consolidate the computation of those universal physical properties. A one-parameter noise model captures the effect of imperfections and reproduces the experimental data. Its systematic study reveals that the noise, analogously to temperature, induces a new length scale in the system. We introduce and successfully verify modified scaling laws, directly accounting for the noise without any prior knowledge. It makes data analyses for extracting physical properties transparent to noise. By understanding how imperfect quantum hardware modifies the genuine properties of quantum states of matter, we enhance the power of NISQ processors considerably for addressing quantum criticality and potentially other phenomena and algorithms.

Quantum Rényi entropy by optimal thermodynamic integration paths

Despite being a well-established operational approach to quantify entanglement, Rényi entropy calculations have been plagued by their computational complexity. We introduce here a theoretical framework based on an optimal thermodynamic integration scheme, where the Rényi entropy can be efficiently evaluated using regularizing paths. This approach avoids slowly convergent fluctuating contributions and leads to low-variance estimates. In this way, large system sizes and high levels of entanglement in model or first-principles Hamiltonians are within our reach. We demonstrate this approach in the one-dimensional quantum Ising model and perform an evaluation of entanglement entropy in the formic acid dimer, by discovering that its two shared protons are entangled even above room temperature.Received 29 December 2021Accepted 9 June 2022DOI:https://doi.org/10.1103/PhysRevResearch.4.L032002Published by the American Physical Society under the terms of the Creative Commons Attribution 4.0 International license. Further distribution of this work must maintain attribution to the author(s) and the published article's title, journal citation, and DOI.Published by the American Physical SocietyPhysics Subject Headings (PhySH)Research AreasChemical Physics & Physical ChemistryEntanglement entropyPath integralsQuantum phase transitionsTechniquesMonte Carlo methodsCondensed Matter, Materials & Applied PhysicsQuantum InformationStatistical Physics

Speeding up quantum adiabatic processes with a dynamical quantum geometric tensor

For adiabatic controls of quantum systems, the non-adiabatic transitions are reduced by increasing the operation time of processes. Perfect quantum adiabaticity usually requires the infinitely slow variation of control parameters. In this paper, we propose the dynamical quantum geometric tensor, as a metric in the control parameter space, to speed up quantum adiabatic processes and reach quantum adiabaticity in relatively short time. The optimal protocol to reach quantum adiabaticity is to vary the control parameter with a constant velocity along the geodesic path according to the metric. For the system initiated from the n-th eigenstate, the transition probability in the optimal protocol is bounded by P_{n}(t)\leq4\mathcal{L}_{n}^{2}/\tau^{2} with the operation time \tau and the quantum adiabatic length \mathcal{L}_{n} induced by the metric. Our optimization strategy is illustrated via two explicit models, the Landau-Zener model and the one-dimensional transverse Ising model.

Error measurements for a quantum annealer using the one-dimensional Ising model with twisted boundaries

A finite length ferromagnetic chain with opposite spin polarization imposed at its two ends is one of the simplest frustrated spin models. In the clean classical limit the domain wall inserted on account of the boundary conditions resides with equal probability on any one of the bonds, and the degeneracy is precisely equal to the number of bonds. If quantum mechanics is introduced via a transverse field, the domain wall will behave as a particle in a box, and prefer to be nearer the middle of the chain rather than the ends. A simple characteristic of a real quantum annealer is therefore which of these limits obtains in practice. Here we have used the ferromagnetic chain with antiparallel boundary spins to test a real flux qubit quantum annealer and discover that contrary to both expectations, the domain walls found are non-uniformly distributed on account of effective random longitudinal fields present notwithstanding tuning carried out to zero out such fields when the couplings between qubits are nominally zero. We present a simple derivation of the form of the distribution function for the domain walls, and show also how the effect we have discovered can be used to determine the strength of the effective random fields (noise) characterizing the annealer. The noise measured in this fashion is smaller than what is seen during the single-qubit tuning process, but nonetheless qualitatively affects the outcome of the simulation performed by the annealer.

Single-ion properties of the transverse-field Ising model materialCoNb2O6

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 pseudo-spin 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 pseudo-spin 1/2 Hamiltonian. These terms are expected to be controlled in part by single-ion physics, namely the wavefunction for the pseudo-spin 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 wavefunctions for this doublet, which we expect will be useful in theoretical investigations of the anisotropic exchange interactions.