The dissipative variant of the Ising model in a transverse field is one of the most important models in the analysis of open quantum many-body systems, due to its paradigmatic character for understanding driven-dissipative quantum phase transitions, as well as its relevance in modeling diverse experimental platforms in atomic physics and quantum simulation. Here, we present an exact solution for the steady state of the transverse-field Ising model in the limit of infinite-range interactions, with local dissipation and inhomogeneous transverse fields. Our solution holds despite the lack of any collective spin symmetry or even permutation symmetry. It allows us to investigate first- and second-order dissipative phase transitions, driven-dissipative criticality, and captures the emergence of a surprising "spin blockade" phenomenon. The ability of the solution to describe spatially varying local fields provides a new tool to study disordered open quantum systems in regimes that would be extremely difficult to treat with numerical methods.
We study quantum criticality in the infinite range transverse-field Ising model. We find subtle differences with respect to the well-known single-site mean-field theory, especially in terms of gap, entanglement and quantum criticality. The calculations are based on numerical diagonalization of Hamiltonians with up to a few thousand spins. This is made possible by the enhanced symmetries of the model, which divide the Hamiltonian into many block-diagonal sectors. The finite temperature phase diagram and the characteristic jump in heat capacity closely resemble the behavior in mean-field theory. However, unlike mean-field theory where excitations are always gapped, the excitation gap in the infinite range model goes to zero from both the paramagnetic side and from the ferromagnetic side on approach to the quantum critical point. Also, contrary to mean-field theory, at the quantum critical point the quantum Fisher information becomes large, implying long-range multipartite entanglement. We find that the main role of temperature is to shift statistical weights from one conserved sector to another. However, low energy excitations in each sector arise only near the quantum critical point implying that low energy quantum fluctuations can arise only in the vicinity of the quantum critical field where they can persist up to temperatures of order the exchange constant.
Quantum simulation of spin models can provide insight into problems that are difficult or impossible to study with classical computers. Trapped ions are an established platform for quantum simulation, but only systems with fewer than 20 ions have demonstrated quantum correlations. We studied quantum spin dynamics arising from an engineered, homogeneous Ising interaction in a two-dimensional array of (9)Be(+) ions in a Penning trap. We verified entanglement in spin-squeezed states of up to 219 ions, directly observing 4.0 ± 0.9 decibels of spectroscopic enhancement, and observed states with non-Gaussian statistics consistent with oversqueezed states. The good agreement with ab initio theory that includes interactions and decoherence lays the groundwork for simulations of the transverse-field Ising model with variable-range interactions, which are generally intractable with classical methods.
Ising models, and the physical systems described by them, play a central role in generating entangled states for use in quantum metrology and quantum information.In particular, ultracold atomic gases, trapped ion systems, and Rydberg atoms realize long-ranged Ising models, which even in the absence of a transverse field can give rise to highly non-classical dynamics and long-range quantum correlations.In the first part of this paper, we present a detailed theoretical framework for studying the dynamics of such systems driven (at time t = 0) into arbitrary unentangled non-equilibrium states, thus greatly extending and unifying the work of Foss-Feig et al (2013 Phys.Rev. A 87 042101).Specifically, we derive exact expressions for closed-time-path ordered correlation functions, and use these to study experimentally relevant observables, e.g.Bloch vector and spin-squeezing dynamics.In the second part, these correlation functions are then used to derive closed-form expressions for the dynamics of arbitrary spin-spin correlation functions in the presence 5
We study the
XY model with infinite-range interactions (Lipkin-Meshkov-Glick model) in the presence of dissipation from spontaneous decay. We show that independent and collective decay lead to qualitatively different phase transitions of the steady state, even though the phase boundary is the same. Independent decay leads to a second-order phase transition to a ferromagnet, while collective decay leads to a first-order transition to a time-dependent oscillatory phase. Then we show that the addition of a drive leads to infinite spin squeezing for collective decay in the thermodynamic limit. Our results can be experimentally seen in trapped-ion and cavity-QED experiments.
Phase transitions in dissipative quantum systems have been investigated using various analytical approaches, particularly in the mean-field (MF) limit. However, analytical results often depend on specific methodologies. For instance, the Keldysh formalism shows that the dissipative transverse Ising (DTI) model exhibits a discontinuous transition at the upper critical dimension,
dc=3; however, the fluctuationless MF approach shows that a continuous transition occurs in infinite dimensions (
d∞). Thus, analytical results must be confirmed by numerical verification. However, this may not be feasible owing to the exponential increase in computational complexity as
O(22N) with system size
N. However, if we assume that the spins are fully connected in infinite dimensions (
d∞), the spin indices are permutation-symmetric. This symmetry significantly reduces the complexity, to
O(N3). Using this symmetry, we simulate the DTI model on fully connected networks with large system sizes. Using finite-size scaling analysis, we show that the upper critical dimension of the DTI model is not
3 but a non-integer dimension
3.5; thus, the discontinuous transition that occurs at high dissipation strength is not an MF result, but is induced by fluctuation effects below
dc. We conclude that the DTI model with strong dissipation exhibits a continuous transition above
dc, as predicted using the fluctuationless MF approach, and we determine the critical exponents.