We propose an alternative concept for the dynamics of a quantum bath, the Chebyshev space, and a method based on this concept, the Chebyshev space method. The Chebyshev space is an abstract vector space that exactly represents the fermionic or bosonic bath degrees of freedom, without a discretization of the bath density of states. Relying on Chebyshev expansions, the Chebyshev space representation of a bath has very favorable properties with respect to extremely precise and efficient calculations of ground state properties, static and dynamical correlations, and time evolution for a great variety of quantum systems. The aim of the present work is to introduce the Chebyshev space in detail and to demonstrate the capabilities of the Chebyshev space method. Although the central idea is derived in full generality, the focus is on model systems coupled to fermionic baths. In particular, we address quantum impurity problems, such as an impurity in a host or a bosonic impurity with a static barrier, and the motion of a wave packet on a chain coupled to leads. For the bosonic impurity, the phase transition from a delocalized electron to a localized polaron in arbitrary dimension is detected. For the wave packet on a chain, we show how the Chebyshev space method implements different boundary conditions, including transparent boundary conditions replacing infinite leads. Furthermore, the self-consistent solution of the Holstein model in infinite dimension is calculated. With the examples, we demonstrate how highly accurate results for system energies, correlation and spectral functions, and time dependence of observables are obtained with modest computational effort.
We propose an alternative concept for the dynamics of a quantum bath, the Chebyshev space, and a method based on this concept, the Chebyshev space method. The Chebyshev space is an abstract vector space that exactly represents the fermionic or bosonic bath degrees of freedom, without a discretization of the bath density of states. Relying on Chebyshev expansions, the Chebyshev space representation of a bath has very favorable properties with respect to extremely precise and efficient calculations of ground state properties, static and dynamical correlations, and time evolution for a great variety of quantum systems. The aim of the present work is to introduce the Chebyshev space in detail and to demonstrate the capabilities of the Chebyshev space method. Although the central idea is derived in full generality, the focus is on model systems coupled to fermionic baths. In particular, we address quantum impurity problems, such as an impurity in a host or a bosonic impurity with a static barrier, and the motion of a wave packet on a chain coupled to leads. For the bosonic impurity, the phase transition from a delocalized electron to a localized polaron in arbitrary dimension is detected. For the wave packet on a chain, we show how the Chebyshev space method implements different boundary conditions, including transparent boundary conditions replacing infinite leads. Furthermore, the self-consistent solution of the Holstein model in infinite dimension is calculated. With the examples, we demonstrate how highly accurate results for system energies, correlation and spectral functions, and time dependence of observables are obtained with modest computational effort.
We present a detailed description of the recently proposed numerical renormalization group method for models of quantum impurities coupled to a bosonic bath. Specifically, the method is applied to the spin-boson model, both in the Ohmic and sub-Ohmic cases. We present various results for static as well as dynamic quantities and discuss details of the numerical implementation, e.g., the discretization of a bosonic bath with arbitrary continuous spectral density, the suitable choice of a finite basis in the bosonic Hilbert space, and questions of convergence with respect to truncation parameters. The method is shown to provide high-accuracy data over the whole range of model parameters and temperatures, which are in agreement with exact results and other numerical data from the literature.
Computer simulations of many-body quantum dynamics of indistinguishable particles is a challenging task for computational physics. In this paper we demonstrate that the method of coupled coherent states (CCS) developed previously for multidimensional quantum dynamics of distinguishable particles can be used to study indistinguishable bosons in the second-quantization formalism. To prove its validity, the technique termed here coupled coherent states for indistinguishable bosons (CCSB) is tested on two model problems. The first is a system-bath problem consisting of a tunneling mode coupled to a harmonic bath, previously studied by CCS and other methods in distinguishable representation in 20 dimensions. The harmonic bath is comprised of identical oscillators, and may be second quantized for use with CCSB, so that this problem may be thought of as a bosonic bath with an impurity. The cross-correlation function for the dynamics of the system and Fourier transform spectrum compare extremely well with a benchmark calculation, which none of the prior methods of studying the problem achieved. The second model problem involves 100 bosons in a shifted harmonic trap. Breathing oscillations in the one-body density are calculated and shown to compare favorably to a multiconfigurational time-dependent Hartree for bosons calculation, demonstrating the applicability of the method as a new formally exact way to study the quantum dynamics of Bose-Einstein condensates.
Theoretical studies of electron and energy transport in nanostructures usually involve nonequilibrium quantum impurity models, which feature a locally interacting quantum system (impurity) coupled to noninteracting bosonic and/or electronic baths each at its own equilibrium. In these models, the interplay of strong correlation, nonequilibrium effects and dissipation leads to rich and complex phenomena. As a nonequilibrium strongly correlated many-body problem, these models, except at some rare cases, have no exact analytical solution, and hence their treatment heavily relies upon numerical techniques. The goal of this thesis is to develop numerically accurate and relatively low-cost techniques for treatment of strongly correlated open nonequilibirum quantum systems. Such impurity solvers are in high demand both as standalone methods for simulation of response in nanoscale junctions and as part of divide-and-conquer schemes such as dynamical mean field theory (DMFT).For nonequilibirum open quantum systems, two characteristic energy scales governing the physics are intra-system interactions and strength of the system-baths couplings. Thus, systems strongly coupled to their baths with weak intra-system interactions or weakly coupled systems with strong local interactions can be treated within developed perturbation theory approaches. My research focuses on development of nonequilibirum dual techniques, which are capable to accurately treat situations where perturbative techniques are inapplicable due to absence of a small parameter in the system and which reduce to standard perturbative expansions in the limits of weak coupling or weak intra-system interaction. \tDual techniques rely on ability to solve exactly a simplified reference system. The latter includes exactly all intra-system interactions and accounts approximately for system-bath couplings. Dual techniques (described below) allow to account for deviation from the correct hybridization function in a organized and controlled manner of superperturbative expansion. Here, I use Markovian Lindbald quantum master equation (QME) as a solver for the reference system. In the thesis I discuss possibility of mapping between true non-Markov dynamics and simpler Markov Lindblad QME (mapping between physical and reference systems), numerical methods for the mapping and for solution of the resulting Lindblad QME, and application of the mapping in novel dual formulations: auxiliary-QME dual fermion and dual boson approaches. \tFirst I analytically prove that the continuous fermionic baths that the impurity is coupled to can be mapped onto a set of dissipative auxiliary modes under Lindblad QME. Based on this mapping, complicated non-Markovian dynamics of the open fermionic system can be treated by a much simpler Markovian Lindblad QME. However, nice mapping generally requires many auxiliary modes, and numerically only a small number of modes can be exactly diagonalized due to exponential increase of the Hilbert space. To perturbatively treat the difference between the true bath and the small-size auxiliary one, we introduce dual-fermion method. It can directly target the steady state and avoid long-time propagation. Besides current-voltage relation, it can also provide spectral functions, indicating itself as a potential impurity solver for DMFT. Furthermore, we generalize dual-fermion method to consider additionally the effects of bosonic environment; hence it is called dual-boson method. It can be used for studying electron-phonon interaction and light-matter interaction in nanojunctions. By comparing with numerically exact results we argue that both methods are very accurate and relatively cheap. \tFinally, matrix product state representation of the Lindblad QME enables us to treat a large number of auxiliary modes, which provide an exponentially improved mapping. Adding counting field to the auxiliary QME, we can compute the generating function of time-dependent full-counting statistics of electron transport, which yields not only electronic current but also higher-order cumulants (e.g., zero-frequency noise).
Solving quantum dynamics is an exponentially difficult problem. Thus, an exact numerical solution is inaccessible for any condensed matter system.\nA promising approach is to divide the system into a quantum subsystem containing degrees of freedom which are of greater interest or those which have more profound quantum character (e.g., have smaller mass) and a classical bath containing the rest of the system.\nImposing such a partition and treating the bath classically results in quantum-classical dynamics. The quantum-classical Liouville equation is a general equation in the Hilbert space of quantum degrees of freedom while it resides in the phase space of the classical degrees of freedom.\nAny numerical solution to this equation requires representation of the quantum subsystem in some basis. Solutions to this equation have been already proposed in the subsystem, adiabatic and force bases, each with its own cons and pros.\n\nIn this work, the quantum-classical equations of motion are cast in the subsystem basis and subsequently mapped to a number of fictitious harmonic oscillators.\nThe result is quantum-classical dynamics in the mapping basis which treats both quantum and classical degrees of freedom on the same footing, i.e., in phase space. Neglecting a portion of the back reaction of the quantum-subsystem to classical bath results in an expression for the time evolution of an operator (density matrix) equal to its Poisson bracket with the Hamiltonian.\nThis equation can be solved in terms of characteristics to provide a computationally tractable method for calculating quantum-classical dynamical properties. The expressions for expectation values and correlation functions in this formalism are derived.\nCalculations on spin-boson system, barrier crossing models---the so called Tully models---and the Fenna-Mathews-Olson pigments show very good agreement between the results of this method and numerical solutions to the Schrödinger equation.
In the first part of the Thesis we mostly concentrate on spectral properties of strongly correlated systems and on their equilibrium properties. This is accomplished by the general concept of imaginary-time dynamics which we apply to a number of different problems in which different strengths of this approach emerge. In Chapter 1 we introduce the formalism that allows for a connection between the quantum and the classical worlds. The connection is established by means of the imaginary-time quantum evolution which, under certain circumstances, is shown to be equivalent to a classical stochastic process. It is further shown that exact static and spectral properties of correlated systems can be obtained when this mapping is feasible. The relationship between the imaginary-time dynamics in different frameworks such as the path-integral and the perturbative one is also underlined. In Chapter 2 we present a specific implementation of the general ideas previously presented. In particular we introduced an extension to lattice systems of the Reptation Monte Carlo algorithm [30] which benefits of a sampling scheme based on directed updates. Specific improvements over the existing methodologies consist in the unbiased evaluation of the imaginary-time path integrals for bosons and a systematic scheme to improve over the Fixed-node approximation for fermions. Applications to the Hubbard and the Heisenberg models are presented. In Chapter 3 we demonstrate the application of the imaginary-time dynamics to the exact study of spectral properties. Subject of our attention is a highly anharmonic and correlated quantum crystal such as Helium 4 at zero temperature.[33] Concerning this system, we have obtained the first ab-initio complete phonon dispersion in good agreement with neutron spectroscopy experiments. Moreover, we have also studied the density excitations of solid helium in a region of wave-vectors in between the collective (phonon) and the single-particle regimes, where the presence of residual coherence in the dynamics shows analogies between the highly anharmonic crystal and the superfluid phase. In Chapter 4 we introduce a novel method, based on the imaginary-time dynamics, to obtain unbiased estimates of fermionic properties.[34] By means of this method and of a very accurate variational state, we provide strong evidence for the stability of a saturated ferromagnetic phase in the high-density regime of the two-dimensional infinite-U Hubbard model. By decreasing the electron density, we observe a discontinuous transition to a paramagnetic phase, accompanied by a divergence of the susceptibility on the paramagnetic side. This behavior, resulting from a high degeneracy among different spin sectors, is consistent with an infinite-order phase transition scenario. In Chapter 5 the use of imaginary-time dynamics in the context of finite-temperature response functions is highlighted. As an application, we study an intriguing quantum phase featuring both glassy order and Bose-Einstein condensation. [35] We introduce and validate a model for the role of geometrical frustration in the coexistence of off-diagonal long range order with an amorphous density profile. The exact characterization of the response of the system to an external density perturbation is what allows here to establish the existence of a spin-glass phase. The differences between such a phase and the otherwise insulating "Bose glasses" are further elucidated in the Chapter. In the second part of the Thesis we focus our attention on the dynamics of closed systems out of equilibrium. This is accomplished by both non-stochastic exact methods for the dynamics and the introduction of a novel time-dependent Variational Monte Carlo scheme. In Chapter 6 exact diagonalization schemes and renormalization-based methods for one-dimensional systems are introduced. We identify key phenomenological traits resulting from the many-body correlation in closed systems driven sufficiently away from equilibrium.[31] We provide evidences that the dynamics of interacting lattice bosons away from equilibrium can be trapped into extremely long-lived inhomogeneous metastable states. The slowing down of incoherent density excitations above a threshold energy, much reminiscent of a dynamical arrest on the verge of a glass transition, is identified as the key feature of this phenomenon. In Chapter 7 we present an extension to dynamical properties of the Variational Quantum Monte Carlo method.[32] This is accomplished by introducing a general class of time-dependent variational states which is based on the mapping of the many-body dynamics onto an instantaneous ground-state problem. The application of the method to the experimentally relevant quantum quenches of interacting bosons reveals the accuracy and the reliability of the introduced numerical scheme. We indeed obtain for the first time a consistent variational description of the approach to the equilibrium of local observables and underline the origin of the metastability and glassy behavior previously identified. In the very last part we draw our conclusions and show some possible paths for stimulating future research.
We develop a method to study quantum impurity models, small interacting quantum systems bilinearly coupled to an environment, in the presence of an additional Markovian quantum bath, with a generic nonlinear coupling to the impurity. We aim at computing the evolution operator of the reduced density matrix of the impurity, obtained after tracing out all the environmental degrees of freedom. First, we derive an exact real-time hybridization expansion for this quantity, which generalizes the result obtained in the absence of the additional Markovian dissipation and which could be amenable to stochastic sampling through diagrammatic Monte Carlo. Then, we obtain a Dyson equation for this quantity and we evaluate its self-energy with a resummation technique known as the noncrossing approximation. We apply this novel approach to a simple fermionic impurity coupled to a zero temperature fermionic bath and in the presence of Markovian pump, losses, and dephasing.