Variational and perturbative extensions of continuous unitary transformations for low-dimensional spin systems
Nils A. Drescher
Abstract
The emergence of complex behavior out of simple elementary building blocks is a typical property of condensed matter systems. While the individual constituents may be well-studied, strongly interacting many-particle systems can give rise to novel collective behavior that can not be understood by considering the individual properties alone. Effective models are a powerful tool to handle this complexity: Instead of aiming at a full description a complex many-particle system on the high-energy scale, an effective model covers only those degrees of freedom which are essential for the low-energy physics. While the high-energy degrees of freedom are integrated out, they influence the low-energy degrees of freedom due to a renormalization of the effective coupling constants. In particular, the concept of the quasi-particle allows for the description of a strongly-interacting many-particle excitations in terms of a simple objects with similar degrees of freedom as a dispersive elementary particle, which gives rise to phenomena like quasi-particle bound states or quasi-particle decay. This thesis aims at the application and methodological improvement of a systematic and well-controlled framework to derive effective low-energy models, denoted as Continuous Unitary Transformations (CUT). In this method, a Hamiltonian is mapped continuously to an effective basis, in which the low-energy degrees of freedom can be interpreted as quasi-particles. Mathematically, the transformation is parametrized by a flow parameter and is represented as a differential equation, denoted as flow equation, while the direction of the transformation is controlled by the generator scheme. As foundation of our studies, we use the sCUT method, in which the types of possible interaction terms are fixed by a truncation scheme. In a common project, we established a new flavor of CUT, called epCUT. It allows us to derive the effective Hamiltonian in a perturbation series. In contrast to the already existent pCUT method, epCUT can handle system with a non-equidistant spectrum as well. Moreover, we discovered the deepCUT method, the non-perturbative twin of epCUT. While it bears large similarities to sCUT, it replaces the cumbersome definition of an appropriate truncation scheme by the order of the flow equation system. Compared to sCUT, deepCUT offers a larger robustness while the ambiguities of the truncation scheme are eliminated. My contribution comprises the computational, parallelized construction of the flow equation system as well. In another project, we investigated the breakdown of the conventional, particle-sorting generator schemes in system with strongly overlapping quasi-particle sub-spaces and identified the emergence of a spectrum unbounded from below due to truncation errors as origin of the divergence. As a testing ground, we used a system of two coupled bosonic degrees of freedom. To overcome these problems, we developed a family of three novel variational generator schemes, denoted as scalar, vectorial and tensorial optimization. After this, we investigated the dimerized, two-dimensional S=1/2 Heisenberg model with an emphasis on the quantum phase transition between a magnetically disordered and a magnetically ordered phase using sCUT. Parametrizing the starting point by a variational parameter, we were able to describe the spontaneous breaking of the SU(2) spin symmetry by CUT for the first time. We were able to calculate the ground state energy and the dispersion of triplons, and magnons, respectively, in both phases. However, we were not able to provide quantitative results for the spontaneous magnetization. Furthermore, we found a way to derive an optimal starting point for the integration of the flow equations that is independent on the particularities of the model and does not require a parametrization by a variational parameter. As a further application, we analyzed the S=1 Heisenberg chain by means of deepCUT. To this end, we mapped the system to an extended S=1/2 Heisenberg ladder with the help of a variational parameter. Using the variational generator schemes developed before, we were able to decouple the low-energy spectrum and to calculate spectral densities. In this investigation, we found an indication for a S=0 bound state. Furthermore, we were able to find some signs of spontaneous quasi-particle decay, but a quantitatively correct description is pending.