Open quantum systems have drawn attention over decades due to its applicability in the foundation of theoretical physics, e.g. statistical mechanics, quantum mechanics and condensed matter physics. The dynamics of open quantum systems has been described as separate entities from their surrounding environment that consists of a very large number of modes, somehow coupled to the mode of the system. Even though the exact solution of the dynamical behavior of the system is impossible to calculate, we obtain a tentative solution using the crucial Markov approximation. The input-output formalism of the quantum Langevin equation (QLE) has been considered as a useful tool which provides a semi-classical description of the dynamics of the system, whereas the master equation provides a complete picture of the dynamics of the system expressed in terms of density matrix. While studying the dynamics of nonlinear system/environment coupling using QLE, we see, for a small value of external field, that the steady state system field does not change much from the steady state field obtained in the absence of nonlinear dissipation. However, in a case where a stronger external field is applied, we see that the deviation becomes substantial from the solution of linear system. We also see that the nonlinear coupling introduces significant difference in the cavity fluctuation spectrum. The description, therefore, provides a potential explanation of parametric effects in terms of nonlinear dissipation phenomena associated with the nonlinear coupling. Even though the theories developed in the context of open quantum systems have proven to be powerful tools, they do not provide a satisfactory platform to be implemented on non-linear Hamiltonians. We often approximate it by linearizing over nonlinear steady state field amplitude, and therefore, the interesting effects are often overlooked. The limitation of the analytics provokes us to simulate open quantum dynamics numerically. The numerical method consists of transformation of the environmental degrees of freedom to a one-dimensional many-body chain, and the computational technique includes numerical diagonalization and renormalization process. The time-adaptive density matrix renormalisation group (t-DMRG) is known as one of the most powerful techniques for the simulation of strongly-correlated many-body quantum systems. In this thesis, along with the theoretical modeling, we implement DMRG numerical scheme for the simulation of canonical S/B model by mapping it to one-dimensional harmonic chain with nearest neighbor interactions, and use the method to investigate the dynamics of the free dissipative system. The thermalization of open quantum systems is also studied by generating minimally entangled typical thermal states (METTS) through imaginary time evolution, and real-time evolving an empty system in the presence of the thermal bath. Further, we simulate coherently driven free dissipative Kerr nonlinear system numerically using Euler's method by solving Heisenberg equation of motion and t-DMRG algorithm, and demonstrate how the numerical results are analogous to classical bistability. By comparing with analytics, we see that the DMRG numerics is analogous to the quantum-mechanical exact solution obtained by mapping the equation of motion of the density matrix of the system to a Fokker-Plank equation. The comparison between two different numerical techniques shows that the semi-classical Euler's method determines the dynamics of the system field of one among two coherent branches, whereas DMRG numerics gives the superposition of both of them. Hence, DMRG-determined time dynamics undergoes generating non-classical states. Our approach of dealing with nonlinearity represents an important contribution in the developments of technique to study the dynamical and steady-state behavior of open quantum systems, which is a fundamental aspect of quantum physics.
Open quantum systems have drawn attention over decades due to its applicability in the foundation of theoretical physics, e.g. statistical mechanics, quantum mechanics and condensed matter physics. The dynamics of open quantum systems has been described as separate entities from their surrounding environment that consists of a very large number of modes, somehow coupled to the mode of the system. Even though the exact solution of the dynamical behavior of the system is impossible to calculate, we obtain a tentative solution using the crucial Markov approximation. The input-output formalism of the quantum Langevin equation (QLE) has been considered as a useful tool which provides a semi-classical description of the dynamics of the system, whereas the master equation provides a complete picture of the dynamics of the system expressed in terms of density matrix. While studying the dynamics of nonlinear system/environment coupling using QLE, we see, for a small value of external field, that the steady state system field does not change much from the steady state field obtained in the absence of nonlinear dissipation. However, in a case where a stronger external field is applied, we see that the deviation becomes substantial from the solution of linear system. We also see that the nonlinear coupling introduces significant difference in the cavity fluctuation spectrum. The description, therefore, provides a potential explanation of parametric effects in terms of nonlinear dissipation phenomena associated with the nonlinear coupling. Even though the theories developed in the context of open quantum systems have proven to be powerful tools, they do not provide a satisfactory platform to be implemented on non-linear Hamiltonians. We often approximate it by linearizing over nonlinear steady state field amplitude, and therefore, the interesting effects are often overlooked. The limitation of the analytics provokes us to simulate open quantum dynamics numerically. The numerical method consists of transformation of the environmental degrees of freedom to a one-dimensional many-body chain, and the computational technique includes numerical diagonalization and renormalization process. The time-adaptive density matrix renormalisation group (t-DMRG) is known as one of the most powerful techniques for the simulation of strongly-correlated many-body quantum systems. In this thesis, along with the theoretical modeling, we implement DMRG numerical scheme for the simulation of canonical S/B model by mapping it to one-dimensional harmonic chain with nearest neighbor interactions, and use the method to investigate the dynamics of the free dissipative system. The thermalization of open quantum systems is also studied by generating minimally entangled typical thermal states (METTS) through imaginary time evolution, and real-time evolving an empty system in the presence of the thermal bath. Further, we simulate coherently driven free dissipative Kerr nonlinear system numerically using Euler's method by solving Heisenberg equation of motion and t-DMRG algorithm, and demonstrate how the numerical results are analogous to classical bistability. By comparing with analytics, we see that the DMRG numerics is analogous to the quantum-mechanical exact solution obtained by mapping the equation of motion of the density matrix of the system to a Fokker-Plank equation. The comparison between two different numerical techniques shows that the semi-classical Euler's method determines the dynamics of the system field of one among two coherent branches, whereas DMRG numerics gives the superposition of both of them. Hence, DMRG-determined time dynamics undergoes generating non-classical states. Our approach of dealing with nonlinearity represents an important contribution in the developments of technique to study the dynamical and steady-state behavior of open quantum systems, which is a fundamental aspect of quantum physics.
Remarkable experimental progress and theoretical advancements in the study of strongly corre- lated systems have opened up new doors towards a better understanding of interesting states of matter. Such states could possibly be used to develop novel technological applications based on the principles of quantum mechanics. In recent years, the investigation of both complex solid state and ultracold atom systems has lead to a number of fascinating discoveries with regards to their non-equilibrium properties. One route to probe non-equilibrium physics is to consider the evolution of open systems defined by the coupling to an external environment introducing dissipation. Engineering the environment according to predefined requirements can then be used to dynamically prepare states which are difficult to access otherwise. This thesis explores the dynamics of many-body open quantum systems aiming at under- standing the effects resulting from the interplay of kinetic motion, interactions and dissipation in different physical settings. First, we build upon recent developments concerning the coupling of ultracold atoms to optical resonators which are subject to photon losses. We show how the dissipation can trigger an attractor dynamics towards a topologically non-trivial steady state that is robust against perturbations. The numerically obtained full open system evolution en- ables us to not only characterize in detail the steady state but also to analyze the features of the approaching dynamics. Second, we study the propagation of two-time correlation functions in interacting quantum spin-1/2 chains in contact with an environment causing dephasing. Pro- vided with the numerically exact full time-evolution, we identify an algebraic scaling regime where time-translation invariance breaks down and uncover the emergence of aging dynam- ics. Finally, we study a system of interacting fermionic atoms in a one-dimensional lattice with particle loss at the central site. There we report on the dissipation-induced deceleration of the non-equilibrium dynamics known as the quantum Zeno effect. As in the majority of cases quantum many-body systems are not solvable analytically, another central aspect of this thesis deals with adapting numerical methods based on matrix product state algorithms to investigate open quantum systems dynamics.
The nonequilibrium effects of dissipation and drive play a key role in an immense variety of nanoscale and mesoscale quantum systems. To understand the behaviour of open quantum systems, we need accurate methods that capture the influence of the environment on the system, while managing the exponentially large Hilbert space required to describe the system. Tensor network algorithms offer an efficient way to approach this challenge. In this thesis, we develop and apply tensor network techniques to study the dynamics and steady states of various open quantum systems. The first part of the thesis focuses on the driven dissipative many body physics in coupled cavity arrays described by Born-Markov master equations. We extend transfer matrix product operator methods to Liouvillian dynamics, and utilize them to compute dynamical correlation functions and fluorescence spectrum of an infinite coupled cavity array in 1D. We also investigate thermalization, and observe the emergence of a quasi-thermal steady state with a negative effective temperature. In another study, we use infinite projected entangled pair state (iPEPS) methods to compute steady states of coupled cavity lattices in 2D. We find that a straightforward adaptation of iPEPS to Liouvillian dynamics is unstable, contradicting a recent publication in the field. The second part investigates more general systems involving strong couplings and structured environments that induce non-Markovian dynamics. We develop a powerful time-evolving matrix product operator (TEMPO) algorithm that builds on Feynman-Vernon influence functional formalism, and uses matrix product states (MPS) to represent the temporal non-Markovian correlations efficiently. We apply TEMPO to study the localization phase transition of the spin-boson model and the dynamics of two spatially separated two-level systems coupled to a common environment. Finally, we propose the Toblerone TEMPO algorithm, which extends TEMPO to many-body systems interacting with general bosonic environments.
The Bose-Hubbard model provides a foundational paradigm for exploring quantum phase transitions in many-body systems. This thesis presents a systematic investigation into this model and its extensions, designed to provide a clear path for understanding both the rich underlying physics and the sophisticated numerical methods required for its study. Our approach is sequential: we begin with the standard one-dimensional Bose-Hubbard model, introduce complexity via the extended model with long-range interactions, and conclude by analyzing the impact of dissipation in a cavity-coupled open quantum system.To navigate these regimes, we implement and apply a suite of powerful computational techniques, including Density Matrix Renormalization Group (DMRG), finite-temperature Quantum Monte Carlo (QMC), and Mean-Field Theory (MFT). A significant component of this work is the detailed exposition of the implementation and logic of these numerical methods. Through their comparative application, we precisely characterize the equilibrium phase diagrams and benchmark the capabilities of each method. Crucially, our analysis of the open-system dynamics reveals that dissipation replaces the equilibrium phase diagram with one governed by non-equilibrium steady states (NESS), where the system evolves to either a stable solution or a persistent limit cycle.As such, this thesis is intended to serve a dual purpose: it is both a rigorous physical study of quantum phase transitions in a condensed matter system and a selfcontained pedagogical guide for researchers beginning to explore the phenomenology of quantum phase transitions through numerical schemes.
Understanding the fate of many-body quantum physics in the presence of noise is of both practical and fundamental interest. Practically, all quantum systems (especially near-term experiments) are interacting with their environment, and so understanding the role that noise plays is crucial to making predictions of observations. Fundamentally, one can find novel physical phenomena in open systems that exists only out of equilibrium. This motivates trying to get a basic understanding of mixed state phases of matter. This dissertation will be in two parts, linked by the common theme of understanding the dynamics and steady states of open, many-body quantum systems. We will first ask about dissipative state preparation, or thinking about how one can use openness in a quantum system to engineer interesting quantum states. We will present a number of exact solutions for spins, free bosons, and free fermions that have unique properties in the steady state, be they highly entangled or topological in nature. Then, we will discuss their dynamics, and present fundamental limits on the time it takes to prepare highly entangled states under local Markovian channels, and also ways that one can use additional ingredients to evade these bounds. The second part will focus on a novel technique to simulate the dynamics of open spin-1/2 systems whose coupling to the environment will formally make them interacting. This technique relies on encoding the many-body density matrix into a stochastic covariance matrix of the Jordan-Wigner fermions which we show can be sampled from efficiently. We use this technique to uncover rich phenomena outlined in a number of many-body systems, as well as provide certain situations in which the stochastic observables can be analytically re-averaged to get closed form solutions for the entire interacting system.
In this thesis, we discuss aspects of the simulation, characterisation, and control of the dynamics of open quantum systems. The latter are systems governed by quantum mechanical laws and interacting with an environment, which is typically much larger than the system itself. According to the standard approach, one is usually able to perform a suitable average over the environmental degrees of freedom, resulting in an effective description of the main system dynamics, which accounts for the effects of the interaction with the surroundings. In standard scenarios, usually studied within the so-called Born-Markov approximation, the system-environment coupling is weak and such that we can perform a neat separation of timescales: the environmental dynamics is assumed to be fast compared to the typical evolution timescale of the system of interest. However, nowadays we are able to inspect physical scenarios where the usual Born-Markov approximation breaks down: the interaction between system and environment can be strong, likely leading to non-negligible memory (non-Markovian) effects. In this work, we discuss some instances in which we cannot work in the standard Born-Markov regime. We introduce and analyse some numerical and analytical techniques to characterise and simulate non-standard scenarios. We show that, even in the Markovian regime, the traditional formulation of thermodynamic irreversibility shows flaws and inconsistencies. We thus use phase-space methods to assess the role of initial correlations shared by the two parties of a bipartite harmonic system in the entropy production rate, including non-Markovian effects. Furthermore, we show that certain interactions enable to redraw the boundaries between the system and the environment in an effective manner, resulting in a new picture where the system degrees of freedom are augmented, while the residual environment is rearranged in such a way that the Born-Markov approximation is recovered. This analytical technique, known as reaction coordinate mapping approach, is employed in this work to show that, upon an accurate choice of the parameter regime, a spin-boson model can serve as a quantum analogue simulator of non-Markovian multiphoton Jaynes-Cummings models. These systematic studies shed light on the thermodynamic characterisation of open quantum systems, as well as on their numerical simulation.
Recent advances in quantum technologies and related experiments have created a need for highly accurate, versatile, and computationally efficient simulation techniques for the dynamics of open quantum systems. Long-lived correlation effects (non-Markovianity), system-environment hybridization, and the necessity for accuracy beyond the Born-Markov approximation form particular challenges. Approaches to meet these challenges have been introduced, originating from different fields, such as hierarchical equations of motion, Lindblad-pseudomode formulas, chain-mapping approaches, quantum Brownian motion master equations, stochastic unravelings, and refined quantum master equations. This diversity, while indicative of the field's relevance, has inadvertently led to a fragmentation that hinders cohesive advances and their effective cross-community application to current problems for complex systems. How are different approaches related to each other? What are their strengths and limitations? Here we give a systematic overview and concise discussion addressing these questions. We make use of a unified framework which very conveniently allows to link different schemes and, this way, may also catalyze further progress. In line with the state of the art, this framework is formulated not in a fully reduced space of the system but in an extended state space which in a minimal fashion includes effective reservoir modes. This in turn offers a comprehensive understanding of existing methods, elucidating their physical interpretations, interconnections, and applicability.