Transverse-field Model Research Articles

Entanglement Dynamics of Random GUE Hamiltonians

In this work, we consider a model of a subsystem interacting with a reservoir and study dynamics of entanglement assuming that the overall time-evolution is governed by non-integrable Hamiltonians. We also compare to an ensemble of Integrable Hamiltonians. To do this, we make use of unitary invariant ensembles of random matrices with either Wigner-Dyson or Poissonian distributions of energy. Using the theory of Weingarten functions, we derive universal average time evolution of the reduced density matrix and the purity and compare these results with several physical Hamiltonians: randomized versions of the transverse field Ising and XXZ models, Spin Glass and, Central Spin and SYK model. The theory excels at describing the latter two. Along the way, we find general expressions for exponential n-point correlation functions in the gas of GUE eigenvalues.

Real- and Imaginary-Time Evolution with Compressed Quantum Circuits

The current generation of noisy intermediate scale quantum computers introduces new opportunities to study quantum many-body systems. In this paper, we show that quantum circuits can provide a dramatically more efficient representation than current classical numerics of the quantum states generated under non-equilibrium quantum dynamics. For quantum circuits, we perform both real- and imaginary-time evolution using an optimization algorithm that is feasible on near-term quantum computers. We benchmark the algorithms by finding the ground state and simulating a global quench of the transverse field Ising model with a longitudinal field on a classical computer. Furthermore, we implement (classically optimized) gates on a quantum processing unit and demonstrate that our algorithm effectively captures real time evolution.

Two-parameter counter-diabatic driving in quantum annealing

We introduce a two-parameter approximate counter-diabatic term into the Hamiltonian of the transverse-field Ising model for quantum annealing to accelerate convergence to the solution, generalizing an existing single-parameter approach. The protocol is equivalent to unconventional diabatic control of the longitudinal and transverse fields in the transverse-field Ising model and thus makes it more feasible for experimental realization than an introduction of new terms such as non-stoquastic catalysts toward the same goal of performance enhancement. We test the idea for the $p$-spin model with $p=3$, which has a first-order quantum phase transition, and show that our two-parameter approach leads to significantly larger ground-state fidelity and lower residual energy than those by traditional quantum annealing as well as by the single-parameter method. We also find a scaling advantage in terms of the time to solution as a function of the system size in a certain range of parameters as compared to the traditional methods.

Z2 Topological Order and Topological Protection of Majorana Fermion Qubits

The Kitaev chain model exhibits topological order that manifests as topological degeneracy, Majorana edge modes and Z2 topological invariant of the bulk spectrum. This model can be obtained from a transverse field Ising model(TFIM) using the Jordan–Wigner transformation. TFIM has neither topological degeneracy nor any edge modes. Topological degeneracy associated with topological order is central to topological quantum computation. In this paper, we explore topological protection of the ground state manifold in the case of Majorana fermion models which exhibit Z2 topological order. We show that there are at least two different ways to understand this topological protection of Majorana fermion qubits: one way is based on fermionic mode operators and the other is based on anti-commuting symmetry operators. We also show how these two different ways are related to each other. We provide a very general approach to understanding the topological protection of Majorana fermion qubits in the case of lattice Hamiltonians. We then show how in topological phases in Majorana fermion models gives rise to new braid group representations. So, we give a unifying and broad perspective of topological phases in Majorana fermion models based on anti-commuting symmetry operators and braid group representations of Majorana fermions as anyons.

Optimal Protocols in Quantum Annealing and Quantum Approximate Optimization Algorithm Problems.

Quantum annealing (QA) and the quantum approximate optimization algorithm (QAOA) are two special cases of the following control problem: apply a combination of two Hamiltonians to minimize the energy of a quantum state. Which is more effective has remained unclear. Here we analytically apply the framework of optimal control theory to show that generically, given a fixed amount of time, the optimal procedure has the pulsed (or "bang-bang") structure of QAOA at the beginning and end but can have a smooth annealing structure in between. This is in contrast to previous works which have suggested that bang-bang (i.e., QAOA) protocols are ideal. To support this theoretical work, we carry out simulations of various transverse field Ising models, demonstrating that bang-anneal-bang protocols are more common. The general features identified here provide guideposts for the nascent experimental implementations of quantum optimization algorithms.

Neutron scattering investigation of proposed Kosterlitz-Thouless transitions in the triangular-lattice Ising antiferromagnet TmMgGaO4

The transverse-field Ising model on the triangular lattice is expected to host an intermediate finite-temperature Kosterlitz-Thouless (KT) phase through a mapping of the spins on each triangular unit to a complex order parameter. TmMgGaO$_4$ is a candidate material to realize such physics due to the non-Kramers nature of Tm$^{3+}$ ion and the resulting two-singlet single-ion ground state. Using inelastic neutron scattering, we confirm this picture by determining the leading parameters of the low-energy effective Hamiltonian of TmMgGaO$_4$. Subsequently, we track the predicted KT phase and related transitions by inspecting the field and temperature dependence of the ac susceptibility. We further probe the spin correlations in both reciprocal space and real space via single crystal neutron diffraction and magnetic total scattering techniques, respectively. Magnetic pair distribution function analysis provides evidence for the formation of vortex-antivortex pairs that characterize the proposed KT phase around 5~K. Although structural disorder influences the field-induced behavior of TmMgGaO$_4$, the magnetism in zero field appears relatively free from these effects. These results position TmMgGaO$_4$ as a strong candidate for a solid-state realization of KT physics in a dense spin system.

Theory of Trotter Error with Commutator Scaling

The Lie-Trotter formula, together with its higher-order generalizations, provides a direct approach to decomposing the exponential of a sum of operators. Despite significant effort, the error scaling of such product formulas remains poorly understood. We develop a theory of Trotter error that overcomes the limitations of prior approaches based on truncating the Baker-Campbell-Hausdorff expansion. Our analysis directly exploits the commutativity of operator summands, producing tighter error bounds for both real- and imaginary-time evolutions. Whereas previous work achieves similar goals for systems with geometric locality or Lie-algebraic structure, our approach holds in general. We give a host of improved algorithms for digital quantum simulation and quantum Monte Carlo methods, including simulations of second-quantized plane-wave electronic structure, $k$-local Hamiltonians, rapidly decaying power-law interactions, clustered Hamiltonians, the transverse field Ising model, and quantum ferromagnets, nearly matching or even outperforming the best previous results. We obtain further speedups using the fact that product formulas can preserve the locality of the simulated system. Specifically, we show that local observables can be simulated with complexity independent of the system size for power-law interacting systems, which implies a Lieb-Robinson bound as a byproduct. Our analysis reproduces known tight bounds for first- and second-order formulas. Our higher-order bound overestimates the complexity of simulating a one-dimensional Heisenberg model with an even-odd ordering of terms by only a factor of $5$, and is close to tight for power-law interactions and other orderings of terms. This suggests that our theory can accurately characterize Trotter error in terms of both asymptotic scaling and constant prefactor.

Quantum phase transitions in long-range interacting hyperuniform spin chains in a transverse field

Hyperuniform states of matter are characterized by anomalous suppression of long-wavelength density fluctuations. While most of interesting cases of disordered hyperuniformity are provided by complex many-body systems like liquids or amorphous solids, classical spin chains with certain long-range interactions have been shown to demonstrate the same phenomenon. It is well-known that the transverse field Ising model shows a quantum phase transition (QPT) at zero temperature. Under the quantum effects of a transverse magnetic field, classical hyperuniform spin chains are expected to lose their hyperuniformity. High-precision simulations of these cases are complicated because of the presence of highly nontrivial long-range interactions. We perform extensive analysis of these systems using density matrix renormalization group to study the possibilities of phase transitions and the mechanism by which they lose hyperuniformity. We discover first-order QPTs in the hyperuniform spin chains. An interesting feature of the phase transitions in these disordered hyperuniform spin chains is that, depending on the parameter values, the presence of transverse magnetic field may remarkably lead to increase in the order of the ground state as measured by the "$\tau$ order metric," even if hyperuniformity is lost. Therefore, it would be possible to design materials to target specific novel quantum behaviors in the presence of a transverse magnetic field. Our numerical investigations suggest that these spin chains can show no more than two QPTs. We further analyze the long-range interacting spin chains via the Jordan-Wigner mapping, showing that under the pairwise interacting approximation and a mean-field treatment, there can be at most two QPTs. Based on these numerical and theoretical explorations, we conjecture that these spin chains can show a maximum of two QPTs at zero temperature.

Hidden symmetries, the Bianchi classification and geodesics of the quantum geometric ground-state manifolds

We study the Killing vectors of the quantum ground-state manifold of a parameter-dependent Hamiltonian. We find that the manifold may have symmetries that are not visible at the level of the Hamiltonian and that different quantum phases of matter exhibit different symmetries. We propose a Bianchi-based classification of the various ground-state manifolds using the Lie algebra of the Killing vector fields. Moreover, we explain how to exploit these symmetries to find geodesics and explore their behaviour when crossing critical lines. We briefly discuss the relation between geodesics, energy fluctuations and adiabatic preparation protocols. Our primary example is the anisotropic transverse-field Ising model. We also analyze the Ising limit and find analytic solutions to the geodesic equations for both cases.

Lattice vibration as a knob on exotic quantum criticality

Control of quantum coherence in a many-body system is one of the key issues in modern condensed matter, and conventional wisdom is that lattice vibration is an innate source of decoherence. Much research has been conducted to eliminate lattice effects. Challenging this wisdom, we show that lattice vibration may not be a decoherence source but an impetus of a novel coherent quantum many-body state. We demonstrate the possibility by studying the transverse-field Ising model on a chain with renormalization group and density-matrix renormalization group methods and theoretically discover a stable $\mathcal{N}=1$ supersymmetric quantum criticality with central charge $c=3/2$. Thus, we propose an Ising spin chain with strong spin-lattice coupling as a candidate to observe supersymmetry. Generic precursor conditions of novel quantum criticality are obtained by generalizing the Larkin-Pikin criterion of thermal transitions. Our work provides the perspective that lattice vibration may be a knob for exotic quantum many-body states.

Multi-critical topological transition at quantum criticality

The investigation and characterization of topological quantum phase transition between gapless phases is one of the recent interest of research in topological states of matter. We consider transverse field Ising model with three spin interaction in one dimension and observe a topological transition between gapless phases on one of the critical lines of this model. We study the distinct nature of these gapless phases and show that they belong to different universality classes. The topological invariant number (winding number) characterize different topological phases for the different regime of parameter space. We observe the evidence of two multi-critical points, one is topologically trivial and the other one is topologically active. Topological quantum phase transition between the gapless phases on the critical line occurs through the non-trivial multi-critical point in the Lifshitz universality class. We calculate and analyze the behavior of Wannier state correlation function close to the multi-critical point and confirm the topological transition between gapless phases. We show the breakdown of Lorentz invariance at this multi-critical point through the energy dispersion analysis. We also show that the scaling theories and curvature function renormalization group can also be effectively used to understand the topological quantum phase transitions between gapless phases. The model Hamiltonian which we study is more applicable for the system with gapless excitations, where the conventional concept of topological quantum phase transition fails.

First- and second-order quantum phase transitions in the one-dimensional transverse-field Ising model with boundary fields

The one-dimensional transverse field Ising model with boundary fields is studied analytically and numerically. The phase diagram in the ordered state is obtained. We find that there exist two types of phase transitions. The longitudinal boundary fields applied to the left and right ends of the Ising chain are ${h}_{L}$ and ${h}_{R}$, respectively. For $|{h}_{L}|,|{h}_{R}|<\sqrt{1\ensuremath{-}g}$, where $g$ is the transverse field and ${h}_{L}{h}_{R}<0$, a first-order phase transition occurs with the changing ${h}_{L}$ or ${h}_{R}$. The energy gap and boundary magnetization are solved exactly. The analytical expressions of the finite-size scaling for the first-order phase transition are obtained. For $|{h}_{R}|>\sqrt{1\ensuremath{-}g}$ and ${h}_{L}{h}_{R}<0$, there exists a continuous phase transition with changing ${h}_{L}$ and vice versa. This transition is identified as a quantum wetting transition. The singularity of the boundary magnetization in this phase transition is explicitly shown. A simple computational procedure with high accuracy and efficiency is proposed to calculate the magnetization. The magnetization profile, correlation functions, and wetting layer thickness are studied numerically.

A real-space renormalization-group calculation for the quantum gauge theory on a square lattice

We revisit Fradkin and Raby’s real-space renormalization-group method to study the quantum gauge theory defined on links forming a two-dimensional square lattice. Following an old suggestion of theirs, a systematic perturbation expansion developed by Hirsch and Mazenko is used to improve the algorithm to second order in an intercell coupling, thereby incorporating the effects of discarded higher energy states. A careful derivation of gauge-invariant effective operators is presented in the Hamiltonian formalism. Renormalization group equations are analyzed near the nontrivial fixed point, reaffirming old work by Hirsch on the dual transverse field Ising model. In addition to recovering Hirsch’s previous findings, critical exponents for the scaling of the spatial correlation length and energy gap in the electric free (deconfined) phase are compared. Unfortunately, their agreement is poor. The leading singular behavior of the ground state energy density is examined near the critical point: we compute both a critical exponent and estimate a critical amplitude ratio.

Superuniversality from disorder at two-dimensional topological phase transitions

We investigate the effects of quenched randomness on topological quantum phase transitions in strongly interacting two-dimensional systems. We focus first on transitions driven by the condensation of a subset of fractionalized quasiparticles (`anyons') identified with `electric charge' excitations of a phase with intrinsic topological order. All other anyons have nontrivial mutual statistics with the condensed subset and hence become confined at the anyon condensation transition. Using a combination of microscopically exact duality transformations and asymptotically exact real-space renormalization group techniques applied to these two-dimensional disordered gauge theories, we argue that the resulting critical scaling behavior is `superuniversal' across a wide range of such condensation transitions, and is controlled by the same infinite-randomness fixed point as that of the 2D random transverse-field Ising model. We validate this claim using large-scale quantum Monte Carlo simulations that allow us to extract zero-temperature critical exponents and correlation functions in (2+1)D disordered interacting systems. We discuss generalizations of these results to a large class of ground-state and excited-state topological transitions in systems with intrinsic topological order as well as those where topological order is either protected or enriched by global symmetries. When the underlying topological order and the symmetry group are Abelian, our results provide prototypes for topological phase transitions between distinct many-body localized phases.

Quantum critical phase transition between two topologically ordered phases in the Ising toric code bilayer

We demonstrate that two toric code layers on the square lattice coupled by an Ising interaction display two distinct phases with intrinsic topological order. The second-order quantum phase transition between the weakly-coupled $\mathbb{Z}_2\times\mathbb{Z}_2$ and the strongly-coupled $\mathbb{Z}_2$ topological order can be described by the condensation of bosonic quasiparticles from both sides and belongs to the 3d Ising$^*$ universality class. This can be shown by an exact duality transformation to the transverse-field Ising model on the square lattice, which builds on the existence of an extensive number of local $\mathbb{Z}_2$ conserved parities. These conserved quantities correspond to the product of two adjacent star operators on different layers. Notably, we show that the low-energy effective model derived about the limit of large Ising coupling is given by an effective single-layer toric code in terms of the conserved quantities of the Ising toric code bilayer. The two topological phases are further characterized by the topological entanglement entropy which serves as a non-local order parameter.

Superradiance, bosonic Peierls distortion, and lattice gauge theory in a generalized Rabi-Hubbard chain

We investigate a one-dimensional Rabi-Hubbard type of model, arranged such that a qdot is sandwiched between every cavity. The role of the qdot is to transmit photons between neighboring cavities, while simultaneously acting as a photon non-linearity. We consider three-level qdots in the $\Lambda$ configuration, where the left and right leg couples exclusively to the left or right cavity. This non-commuting interaction leads to two highly entangled incompressible phases, separated by a second order quantum phase transition: the qdot degrees-of-freedom act as a dynamical lattice for the photons and a Peierls instability breaks a second $\mathbb{Z}_2$ symmetry which leads to a dimerization in entanglement and photon number. We also find a normal insulating phase and a superfluid phase that acts as a quantum many-body superradiant phase. In the superradiant phase, a $\mathbb{Z}_2$ symmetry is broken and the phase transition falls within the transverse field Ising model universality class. Finally, we show that a limit of the model can be interpreted as a $\mathbb{Z}_2$ lattice gauge theory.

Detecting Out-of-Time-Order Correlations via Quasiadiabatic Echoes as a Tool to Reveal Quantum Coherence in Equilibrium Quantum Phase Transitions.

We propose a new dynamical method to connect equilibrium quantum phase transitions and quantum coherence using out-of-time-order correlations (OTOCs). Adopting the iconic Lipkin-Meshkov-Glick and transverse-field Ising models as illustrative examples, we show that an abrupt change in coherence and entanglement of the ground state across a quantum phase transition is observable in the spectrum of multiple quantum coherence intensities, which are a special type of OTOC. We also develop a robust protocol to obtain the relevant OTOCs using quasi-adiabatic quenches through the ground state phase diagram. Our scheme allows for the detection of OTOCs without time reversal of coherent dynamics, making it applicable and important for a broad range of current experiments where time reversal cannot be achieved by inverting the sign of the underlying Hamiltonian.

Exploring Entanglement and Optimization within the Hamiltonian Variational Ansatz

Quantum variational algorithms are one of the most promising applications of near-term quantum computers; however, recent studies have demonstrated that unless the variational quantum circuits are configured in a problem-specific manner, optimization of such circuits will most likely fail. In this paper, we focus on a special family of quantum circuits called the Hamiltonian Variational Ansatz (HVA), which takes inspiration from the quantum approximation optimization algorithm and adiabatic quantum computation. Through the study of its entanglement spectrum and energy gradient statistics, we find that HVA exhibits favorable structural properties such as mild or entirely absent barren plateaus and a restricted state space that eases their optimization in comparison to the well-studied "hardware-efficient ansatz." We also numerically observe that the optimization landscape of HVA becomes almost trap free when the ansatz is over-parameterized. We observe a size-dependent "computational phase transition" as the number of layers in the HVA circuit is increased where the optimization crosses over from a hard to an easy region in terms of the quality of the approximations and speed of convergence to a good solution. In contrast with the analogous transitions observed in the learning of random unitaries which occur at a number of layers that grows exponentially with the number of qubits, our Variational Quantum Eigensolver experiments suggest that the threshold to achieve the over-parameterization phenomenon scales at most polynomially in the number of qubits for the transverse field Ising and XXZ models. Lastly, as a demonstration of its entangling power and effectiveness, we show that HVA can find accurate approximations to the ground states of a modified Haldane-Shastry Hamiltonian on a ring, which has long-range interactions and has a power-law entanglement scaling.

Probing eigenstate thermalization in quantum simulators via fluctuation-dissipation relations

The eigenstate thermalization hypothesis (ETH) offers a universal mechanism for the approach to equilibrium of closed quantum many-body systems. So far, however, experimental studies have focused on the relaxation dynamics of observables as described by the diagonal part of ETH, whose verification requires substantial numerical input. This leaves many of the general assumptions of ETH untested. Here, we propose a theory-independent route to probe the full ETH in quantum simulators by observing the emergence of fluctuation-dissipation relations, which directly probe the off-diagonal part of ETH. We discuss and propose protocols to independently measure fluctuations and dissipations as well as higher-order time ordered correlation functions. We first show how the emergence of fluctuation dissipation relations from a nonequilibrium initial state can be observed for the 2D Bose-Hubbard model in superconducting qubits or quantum gas microscopes. Then we focus on the long-range transverse field Ising model (LTFI), which can be realized with trapped ions. The LTFI exhibits rich thermalization phenomena: For strong transverse fields, we observe prethermalization to an effective magnetization-conserving Hamiltonian in the fluctuation dissipation relations. For weak transverse fields, confined excitations lead to non-thermal features resulting in a violation of the fluctuation-dissipation relations up to long times. Moreover, in an integrable region of the LTFI, thermalization to a generalized Gibbs ensemble occurs and the fluctuation-dissipation relations enable an experimental diagonalization of the Hamiltonian. Our work presents a theory-independent way to characterize thermalization in quantum simulators and paves the way to quantum simulate condensed matter pump-probe experiments.

Asymptotically Exact Theory for Nonlinear Spectroscopy of Random Quantum Magnets.

We study nonlinear response in quantum spin systems near infinite-randomness critical points. Nonlinear dynamical probes, such as two-dimensional (2D) coherent spectroscopy, can diagnose the nearly localized character of excitations in such systems. We present exact results for nonlinear response in the 1D random transverse-field Ising model, from which we extract information about critical behavior that is absent in linear response. Our analysis yields exact scaling forms for the distribution functions of relaxation times that result from realistic channels for dissipation in random magnets. We argue that our results capture the scaling of relaxation times and nonlinear response in generic random quantum magnets in any spatial dimension.