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.
We develop variational continuous unitary transformations (VCUTs), which integrate Wegner-Wilson flow equations with tensor network techniques to approximately diagonalize many-body localized (MBL) Hamiltonians. The diagonalizing unitary is represented as a matrix product operator whose bond dimension controls the accuracy. For the disordered Heisenberg chain, VCUTs accurately reproduce the full spectrum across the ergodic-to-MBL crossover at small system sizes and scale to L = 48 sites. Beyond eigenenergies, the method can track the spatial entanglement structure of the diagonalizing unitary U ( l ) at each flow step, enabling identification of local integrals of motion deep in the MBL phase.
Tensor network methods are playing a central role in multiple disciplines of modern quantum physics, such as condensed matter physics, quantum information and quantum gravity. Based on renormalization group (RG) ideas, they provide a large class of variational ansatz for quantum many-body systems. There are various tensor network structures which are suited for different systems based on their entanglement structures, such as matrix product states (MPS), projected entangled pair states (PEPS), multi-scale entanglement renormalization ansatz (MERA). Numerous techniques involving coarse-graining tensor networks have been developed, such as tensor renormalization group (TRG), tensor network renormalization (TNR) and its variants. \n \nFirst, we study the generalization of MERA to continuous systems, or cMERA, which is expected to become a powerful variational ansatz for the ground state of strongly interacting quantum field theories. We investigate, in the simpler context of Gaussian cMERA for free theories, the extent to which the cMERA state
ketPsicM with finite UV cut-off
Lambda can capture the spacetime symmetries of the ground state
ketPsi. For a free boson conformal field theory (CFT) in 1+1 dimensions as a concrete example, we build a quasi-local unitary transformation
V that maps
ketPsi into
ketPsicM and show two main results. (i) Any spacetime symmetry of the ground state
ketPsi is also mapped by
V into a spacetime symmetry of the cMERA
ketPsicM. However, while in the CFT the stress-energy tensor
Tmunu(x) (in terms of which all the spacetime symmetry generators are expressed) is local, the corresponding cMERA stress-energy tensor
TmunucM(x)=VTmunu(x)Vdagger is quasi-local. (ii) From the cMERA, we can extract quasi-local scaling operators
OalphacM(x) characterized by the exact same scaling dimensions
Deltaalpha, conformal spins
salpha, operator product expansion coefficients
Calphabetagamma, and central charge
c as the original CFT. We argue that these results should also apply to interacting theories. \n \nSecond we extend TNR to field theories in the continuum. A short-distance length scale
1/Lambda is introduced in the continuum partition function by smearing the fields. The resulting object is still defined in the continuum but has no fluctuations at distances shorter than
1/Lambda. An infinitesimal coarse-graining step is then generated by the combined action of a \\textit{rescaling} operator
mathfrakL and a \\textit{disentangling} operator
mathfrakK that implements a quasilocal field redefinition. As demonstrated for a free boson in two dimensions, continuous TNR exactly preserves translation and rotation symmetries and can generate a proper RG flow. Moreover, from a critical fixed point of this RG flow one can then extract the conformal data of the underlying conformal field theory. \n \nLeaving behind tensor network methods, we explore the universal structure that emerges from the eigenvectors of the reduced density matrix (RDM) for 1+1 dimensional critical lattice systems. Using the Ising model as an example, we demonstrate that the large-weight eigenvectors are related to the low-lying states of the underlying CFT with suitable boundary conditions by a conformal map, and they transform accordingly under the Virasoro generators. \n \nFinally, using the quantum renormalization group (QRG), we derive the bulk geometry that emerges in the holographic dual of the fermionic U(
N) vector model at a nonzero charge density. The obstruction that prohibits the metallic state from being smoothly deformable to the direct product state under the renormalization group flow gives rise to a horizon at a finite radial coordinate in the bulk. The region outside the horizon is described by the Lifshitz geometry with a higher-spin hair determined by microscopic details of the boundary theory. On the other hand, the interior of the horizon is not described by any Riemannian manifold, as it exhibits an algebraic non-locality. The non-local structure inside the horizon carries the information on the shape of the filled Fermi sea.
We introduce the concept of concatenated tensor networks to efficiently describe quantum states. We show that the corresponding concatenated tensor network states can efficiently describe time evolution and possess arbitrary block-wise entanglement and long-ranged correlations. We illustrate the approach for the enhancement of matrix product states, i.e. one-dimensional (1D) tensor networks, where we replace each of the matrices of the original matrix product state with another 1D tensor network. This procedure yields a 2D tensor network, which includes—already for tensor dimension 2—all states that can be prepared by circuits of polynomially many (possibly non-unitary) two-qubit quantum operations, as well as states resulting from time evolution with respect to Hamiltonians with short-ranged interactions. We investigate the possibility of efficiently extracting information from these states, which serves as the basic step in a variational optimization procedure. To this aim, we utilize the known exact and approximate methods for 2D tensor networks and demonstrate some improvements thereof, which are also applicable e.g. in the context of 2D projected entangled pair states. We generalize the approach to higher dimensional and tree tensor networks.
Directly evaluated enhanced perturbative continuous unitary transformations (deepCUTs) are used to calculate non-perturbatively extrapolated numerical data for the ground-state energy and the energy gap. The data coincide with numerical evaluations of the truncated perturbative series and provide robust extrapolations beyond the perturbative regime. We develop a general scheme to extract quantum-critical properties from the deepCUT data based on critical scaling and a strict correspondence between the truncation used for deepCUT and the length scale of correlations at the critical point. We apply our approach to transverse-field Ising models (TFIMs) as paradigmatic systems for quantum phase transitions of various universality classes depending on the lattice geometry and the choice of antiferromagnetic or ferromagnetic coupling. In particular, we focus on the quantum phase diagram of the bilayer antiferromagnetic TFIM on the triangular lattice with an Ising-type interlayer coupling. Without a field, the model is known to host a classically disordered ground state, and in the limit of decoupled layers it exhibits the 3d-XY ‘order by disorder’ transition of the corresponding single-layer model. Our starting point for the unknown parts of the phase diagram is a high-order perturbative calculation about the limit of isolated dimers where the model is in a gapped phase.
Abstract Tensor network methods as presented in our open source matrix product states software have opened up the possibility to study many-body quantum physics in one and quasi-one-dimensional systems in an easily accessible package similar to density functional theory codes but for strongly correlated dynamics. Here, we address methods which allow one to capture the full entanglement without truncation of the Hilbert space. Such methods are suitable for validation of and comparisons to tensor network algorithms, but especially useful in the case of new kinds of quantum states with high entanglement violating the truncation in tensor networks. Quantum cellular automata are one example for such a system, characterized by tunable complexity, entanglement, and a large spread over the Hilbert space. Beyond the evolution of pure states as a closed system, we adapt the techniques for open quantum systems simulated via the Lindblad master equation. We present three algorithms for solving closed-system many-body time evolution without truncation of the Hilbert space. Exact diagonalization methods have the advantage that they not only keep the full entanglement but also require no approximations to the propagator. Seeking the limits of a maximal number of qubits on a single core, we use Trotter decompositions or Krylov approximation to the exponential of the Hamiltonian. All three methods are also implemented for open systems represented via the Lindblad master equation built from local channels. We show their convergence parameters and focus on efficient schemes for their implementations including Abelian symmetries, e.g. symmetry used for number conservation in the Bose–Hubbard model or discrete symmetries in the quantum Ising model. We present the thermalization timescale in the long-range quantum Ising model as a key example of how exact diagonalization contributes to novel physics.