We propose a generalization of the linked-cluster expansions to study driven-dissipative quantum lattice models, directly accessing the thermodynamic limit of the system. Our method leads to the evaluation of the desired extensive property onto small connected clusters of a given size and topology. We first test this approach on the isotropic spin-1/2 Hamiltonian in two dimensions, where each spin is coupled to an independent environment that induces incoherent spin flips. Then we apply it to the study of an anisotropic model displaying a dissipative phase transition from a magnetically ordered to a disordered phase. By means of a Pad\'e analysis on the series expansions for the average magnetization, we provide a viable route to locate the phase transition and to extrapolate the critical exponent for the magnetic susceptibility.
We propose a generalization of the linked-cluster expansions to study driven-dissipative quantum lattice models, directly accessing the thermodynamic limit of the system. Our method leads to the evaluation of the desired extensive property onto small connected clusters of a given size and topology. We first test this approach on the isotropic spin-1/2 Hamiltonian in two dimensions, where each spin is coupled to an independent environment that induces incoherent spin flips. Then we apply it to the study of an anisotropic model displaying a dissipative phase transition from a magnetically ordered to a disordered phase. By means of a Pad'e analysis on the series expansions for the average magnetization, we provide a viable route to locate the phase transition and to extrapolate the critical exponent for the magnetic susceptibility.
Abstract We generalize the numerical linked-cluster expansion (NLCE) method to study the dissipative quantum many-body system. We apply the NLCE to the triangular strip and two-dimensional lattice system. We investigate the dynamics and steady-state properties of the dissipative XYZ model where the coherent dynamics is governed by the anisotropic Heisenberg Hamiltonian while the nonunitary process is induced by the incoherent spin flips. By comparing with the quantum trajectory simulations, the NLCE results show good performance in capturing the dynamics of system with short-range correlations. For strong and long-range correlated system, the larger size clusters in the series should be included. The NLCE study for the magnetic susceptibility also signals the steady-state paramagnetic-ferromagnetic phase transition in the two-dimensional case.
Similar to various series expansions that are used to approximate mathematical func- tions, the linked-cluster expansion is an approximation method that allows us to approach the actual values of a very large physical system’s different physical quan- tities by systematically studying smaller systems embedded in this larger system. The main concept in linked-cluster expansion, weight, represents the additional con- tribution to a certain physical quantity by increasing the system size by one unit. These weights are used to eventually build up the result on a larger system. In our case, we focus on the partition function, a quantity that can be used to calculate several essential thermodynamic aspects of the system such as average energy, spe- cific heat, magnetic susceptibility, etc. We study these weights for the Ising and Potts models on a one dimensional lattice and Bethe lattice, with and without an external magnetic field. Previous studies have shown that in the one dimensional Ising model, adding more than two lattice sites to the system does not create any ad- ditional contribution to the partition function (i.e. their weights are all zero), giving us the result for an infinitely large system with just two lattice sites, a remarkable simplification. In our study, we prove that this property holds not only for the Ising model, but also its generalization, the Potts model, and show that it is a result of a spin flip symmetry inherent in the system. In order to test this, we break this symmetry using an externally applied magnetic field and show that in this case, for magnetic dipole energies comparable with the exchange energy between neighboring spins, this special property vanishes.
For arbitrary space dimension d, we investigate the quantum phase transitions of two paradigmatic spin models defined on a hypercubic lattice, the coupled-dimer Heisenberg model and the transverse-field Ising model. To this end, high-order linked-cluster expansions for the ground-state energy and the one-particle gap are performed up to order 9 about the decoupled-dimer and high-field limits, respectively. Extrapolations of the high-order series yield the location of the quantum phase transition and the correlation-length exponent $\ensuremath{\nu}$ as a function of space dimension d. The results are complemented by 1/d expansions to next-to-leading order of observables across the phase diagrams. Remarkably, our analysis of the extrapolated linked-cluster expansion allows to extract the coefficients of the full 1/d expansion for the phase-boundary location in both models exactly in leading order and quantitatively for subleading corrections.
We use numerical linked cluster expansions (NLCEs) to study correlated lattice systems in one and two dimensions.For fermions in the honeycomb lattice, we study its finite-temperature properties and short-range spin correlations using NLCEs and determinantal quantum Monte Carlo (DQMC) simulations.For the homogeneous system, we analyze a number of thermodynamic quantities, including the entropy, the specific heat, uniform and staggered spin susceptibilities, short-range spin correlations, and the double occupancy at and away from half filling.Employing a local density approximation (LDA), we examine the viability of adiabatic cooling by increasing the interaction strength for homogeneous as well as for trapped systems.We also use NLCEs to study thermodynamic properties of the two-dimensional spin-1/2 Ising, XY, and Heisenberg models with bimodal random-bond disorder on the square and honeycomb lattices.In all cases, the nearest-neighbor coupling between the spins takes values ±J with equal probability.We obtain the disorder averaged (over all disorder configurations) energy, entropy, specific heat, and uniform magnetic susceptibility in each case.These results are compared with the corresponding ones in the clean models.Analytic expressions are obtained for low orders in the expansion of these thermodynamic quantities in inverse temperature.For many-body localization in disordered isolated systems, we study quantum quenches in the thermodynamic limit.By a quantum quench it is meant that the initial state is stationary with respect iii to an initial Hamiltonian, which is suddenly changed to a new (time-independent) Hamiltonian.The latter then drives the (unitary) dynamics of the system.We are interested in the time average of observables after relaxation following the quench, which can be computed under the so-called diagonal ensemble.
The linked-cluster series expansion is developed fur use on models of quantum spin systems with single-ion potentials. The linked-cluster series is defined as a statistical mechanical perturbation series in which the mean-field Hamiltonian and the single-ion potential are treated exactly with the perturbation representing the effects of fluctuation correlations. All terms of the free energy up to a given finite order are calculated, and estimates of thermodynamic functions are made based on extrapolation methods including the ratio method and Pade approximants.
Imperfections in correlated materials can alter their ground state as well as finite-temperature properties in significant ways. Here, we develop a method based on numerical linked-cluster expansions for calculating exact finite-temperature properties of disordered lattice models directly in the thermodynamic limit. We show that a continuous distribution for disordered parameters can be achieved using a set of carefully chosen discrete modes in the distribution, which allows for the averaging of properties over all disorder realizations. We benchmark our results for thermodynamic properties of the square-lattice Ising and quantum Heisenberg models with bond disorder against Monte Carlo simulations and study them as the strength of disorder changes. We also apply the method to the disordered Heisenberg model on the frustrated checkerboard lattice, which is closely connected to Sr2Cu(Te0.5W0.5)O6. Our method can be used to study finite-temperature properties of other disordered quantum lattice models, including those for interacting lattice fermions.