No Arabic abstract
We identify a fundamental challenge for non-perturbative linked cluster expansions (NLCEs) resulting from the reduced symmetry on graphs, most importantly the breaking of translational symmetry, when targeting the properties of excited states. A generalized notion of cluster additivity is introduced, which is used to formulate an optimized scheme of graph-based continuous unitary transformations (gCUTs) allowing to solve and to physically understand this fundamental challenge. Most importantly, it demands to go beyond the paradigm of using the exact eigenvectors on graphs.
We develop a numerical linked cluster expansion (NLCE) method that can be applied directly to inhomogeneous systems, for example Hamiltonians with disorder and dynamics initiated from inhomogeneous initial states. We demonstrate the method by calculating dynamics for single-spin expectations and spin correlations in two-dimensional spin models on a square lattice, starting from a checkerboard state. We show that NLCE can give moderate to dramatic improvement over an exact diagonalization of comparable computational cost, and that the advantage in computational resources grows exponentially as the size of the clusters included grows. Although the method applies to any type of NLCE, our explicit benchmarks use the rectangle expansion. Besides showing the capability to treat inhomogeneous systems, these benchmarks demonstrate the rectangle expansions utility out of equilibrium.
We develop strong-coupling series expansion methods to study two-particle spectra of quantum lattice models. At the heart of the method lies the calculation of an effective Hamiltonian in the two-particle subspace. We explicitly consider an orthogonality transformation to generate this block diagonalization, and find that maintaining orthogonality is crucial for systems where the ground state and the two-particle subspace are characterized by identical quantum numbers. We discuss the solution of the two-particle Schrodinger equation by using a finite lattice approach in coordinate space or by an integral equation in momentum space. These methods allow us to precisely determine the low-lying excitation spectra of the models at hand, including all two-particle bound/antibound states. Further, we discuss how to generate series expansions for the dispersions of the bound/antibound states. These allow us to employ series extrapolation techniques, whereby binding energies can be determined even when the expansion parameters are not small. We apply the method to the (1+1)D transverse Ising model and the two-leg spin-$case 1/2$ Heisenberg ladder. For the latter model, we also calculate the coherence lengths and determine the critical properties where bound states merge with the two-particle continuum.
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 Pade 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 generalize the family of approximate momentum average methods to formulate a numerically exact, convergent hierarchy of equations whose solution provides an efficient algorithm to compute the Greens function of a particle dressed by bosons suitable in the entire parameter regime. We use this approach to extract ground-state properties and spectral functions. Our approximation-free framework, dubbed the generalized Greens function cluster expansion (GGCE), allows access to exact numerical results in the extreme adiabatic limit, where many standard methods struggle or completely fail. We showcase the performance of the method, specializing three important models of charge-boson coupling in solids and molecular complexes: the molecular Holstein model, which describes coupling between charge density and local distortions, the Peierls model, which describes modulation of charge hopping due to intersite distortions, and a more complex Holstein+Peierls system with couplings to two different phonon modes, paradigmatic of charge-lattice interactions in organic crystals. The GGCE serves as an efficient approach that can be systematically extended to different physical scenarios, thus providing a tool to model the frequency dependence of dressed particles in realistic settings.
High order perturbation theory has seen an unexpected recent revival for controlled calculations of quantum many-body systems, even at strong coupling. We adapt integration methods using low-discrepancy sequences to this problem. They greatly outperform state-of-the-art diagrammatic Monte Carlo. In practical applications, we show speed-ups of several orders of magnitude with scaling as fast as $1/N$ in sample number $N$; parametrically faster than $1/sqrt{N}$ in Monte Carlo. We illustrate our technique with a solution of the Kondo ridge in quantum dots, where it allows large parameter sweeps.