Numerous examples in condensed matter physics possess an effective Hamiltonian structure involving two-body anharmonic interactions between elementary excitations. The effective second-quantised Hamiltonian used to model these systems is essentially of the same form except for the particulars of the interaction parameters and the statistics of the particles which are interacting. The author investigates such systems in the limit of strongly anharmonic interactions which may possibly lead to symmetry breaking. He introduces quantum field operators and calculates their equations of motion. The interaction coefficients are expanded in a Taylor series about the critical mode to second order to illustrate the method, but he gives the expansion to infinite order for completeness. It is demonstrated that both in first and second order the equations of motion are reduced to a system of highly nonlinear coupled partial differential equations (PDE) for the field operator. He shows that the form of the second-order equation is exact due to renormalisation theory, in the case where the Hamiltonian coefficients represent operators in real space which are isotropic, and incorporates symmetries and topologies which would appear in an expansion to infinite order. This special situation arises in a very wide range of physical examples. These equations can be solved exactly using a particular ansatz in conjunction with recent mathematical developments in the field of non-linear analysis. In particular, the method of symmetry reduction for PDE provides a complete set of symmetry variables for the type of equations used in the study.
Numerous examples in condensed matter physics possess an effective Hamiltonian structure involving two-body anharmonic interactions between elementary excitations. The effective second-quantised Hamiltonian used to model these systems is essentially of the same form except for the particulars of the interaction parameters and the statistics of the particles which are interacting. The author investigates such systems in the limit of strongly anharmonic interactions which may possibly lead to symmetry breaking. He introduces quantum field operators and calculates their equations of motion. The interaction coefficients are expanded in a Taylor series about the critical mode to second order to illustrate the method, but he gives the expansion to infinite order for completeness. It is demonstrated that both in first and second order the equations of motion are reduced to a system of highly nonlinear coupled partial differential equations (PDE) for the field operator. He shows that the form of the second-order equation is exact due to renormalisation theory, in the case where the Hamiltonian coefficients represent operators in real space which are isotropic, and incorporates symmetries and topologies which would appear in an expansion to infinite order. This special situation arises in a very wide range of physical examples. These equations can be solved exactly using a particular ansatz in conjunction with recent mathematical developments in the field of non-linear analysis. In particular, the method of symmetry reduction for PDE provides a complete set of symmetry variables for the type of equations used in the study.
For pt.I see ibid. vol.22, p.4877 (1989). Previously (I) the author considered a second-quantised Hamiltonian which can be used to model a large number of important strongly interacting many-body systems. He has shown that the dynamics of such systems can be described exactly using quantum fields. The equations of motion derived in this way are highly nonlinear partial differential equations (PDE). He now uses a standard field theoretical approach where initially the fields are treated as classical functions. He studies the solutions of these equations which remarkably, may be obtained making full use of very recent mathematical discoveries (the symmetry reduction method). He finds exact solutions of the equations of motion for the classical field for each of the four main cases separately found in the preceding paper. These cases may correspond to different physical situations and the last, the second-order case, corresponds to the most general physical situation. Subsequently he outlines possible quantisation procedures for these fields. He then briefly discusses the type of boundary conditions which may be applied to both the classical field and the quantum fluctuations which arise.
Mean-field Hartree theory is a central tool for mapping interacting many-body dynamics into an effective nonlinear one-particle evolution. This approximation has been also employed when the Hamiltonian that governs the many-body dynamics is not Hermitian. Indeed, non-Hermitian Hamiltonians model particle gain/loss or the evolution of open quantum systems between consecutive quantum jumps. Furthermore, the validity of the Hartree approximation for generic non-Hermitian Hamiltonians lies at the basis of a quantum algorithm for nonlinear differential equations, where the nonlinear equation is mapped into a many-body Schr\"odinger equation, which may be non-Hermitian. In this application, the mean-field approximation can be seen as a linearization method, and at the end of the algorithm, the obtained linear system is solved using the Harrow-Hassidim-Lloyd (HHL) algorithm. The non-Hermitian mean-field Hartree approximation is far from being completely understood, and no rigorous results justify its applications. Instead, it is usually presented as a heuristic extension of the Hermitian case. In the first part of this thesis, we show that this approximation can fail. We consider the most general bosonic mean-field Hamiltonian with two-body interactions and we prove that, if the Hamiltonian is not Hermitian, the Hartree equation does not capture the correct second time derivative of the one-particle marginal state at t=0 in the limit of infinitely many particles. We construct a series of examples that we can solve analytically or numerically, which support the failure of the Hartree approximation and give further insight into the new phenomena that arise in non-Hermitian physics. Furthermore, we consider the mean-field evolution of open quantum systems and we show how to use the mean-field approximation of open quantum systems as a linearization method of nonlinear differential equations. In the second part of this thesis, we study the mean-field dynamics of a Hermitian system of N interacting bosons starting from an initially condensed state. For a broad class of mean-field Hamiltonians, including models with arbitrary bounded interactions and models with unbounded interaction potentials, we prove that the probability of having n particles outside the condensate decays exponentially in n for any finite evolution time. Our results strengthen previously known bounds that provide only polynomial control on the probability of having n excitations. Finally, we consider the HHL algorithm. Several circuit-level variants have been proposed to simplify its implementation, including constructions in which the clock register is factorized via Hadamard gates to reduce circuit depth and entanglement. We provide a rigorous analysis of the convergence of the output towards the ideal solution, and we show that this commonly deployed variant can fundamentally alter the asymptotic convergence of the algorithm. This phenomenon is not a finite-size artifact and compromises the correctness guarantees of the factorized-clock implementation. We provide an analytic characterization of the underlying mechanism and propose a minimal and experimentally feasible modification that suffices to restore convergence while preserving the practical advantages of the simplified circuit. Our results clarify an overlooked correctness issue in a widely used HHL implementation pattern and delineate the conditions under which simplified realizations faithfully reproduce the intended algorithmic behavior.
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.
The dynamics of interacting many-particle systems on quantum mechanical scale is a broad subject of research in modern physics. The understanding of the quantum correlations, or entanglement, between the particles in an isolated many-body quantum system may however be challenging, as the amount of interactions may be large. To be able to describe the system in a simpler way, mean field theory is often applied; the interactions of all individual particles are substituted by an averaged effect. First and second order mean field approximations are applied to the equations of motion of the amplitudes <ân> for each mode of the many-body system. In first order mean field, under the assumption that each mode has the same magnitude of the (constant) amplitude, the Kuramoto equation follows for the phases of the complex amplitudes. To study this first order approximation and to refine the mean field solution, the equations of motion have been extended to second order mean field. After that, numerical integration is used to solve the system of coupled differential equations for two modes. The solutions seem unstable, which is verified by defining quantum fluctuations beyond mean field. These are shown not to be negligible for a two mode system.It was discovered that the mean field equations of motion show divergent fluctuations for a system of two modes, which opens a whole set of questions on the validity of the first order mean field approximation and therefore the use of the Kuramoto model for quantum many-body dynamics. The extension to the second order mean field solution might be promising as this does include correlations between the operators.
This thesis investigates the rigorous derivation of non-linear effective equations describing the time evolution of Bose-Einstein condensation from the N-particle linear Schrödinger equation, focusing on the one-particle reduced density matrix. The first chapter outlines the mathematical framework of many-body quantum mechanics, including the Fock space formalism. The second chapter studies the mean-field regime, proving that the condensed state converges to the non-linear Hartree equation by analyzing fluctuation dynamics to quantify the convergence rate. The third chapter addresses the physically realistic but mathematically challenging Gross-Pitaevskii regime, characterized by rare but strong interactions. Here, Bogoliubov transformations are used to capture crucial microscopic short-scale correlations, allowing for a quantitative derivation of the cubic non-linear Schrödinger equation. While the preceding chapters focus exclusively on repulsive interactions, the fourth chapter explores the focusing case with attractive potentials, which can lead to finite-time blow-up. Operating in an intermediate scaling regime, the effective focusing cubic non-linear Schrödinger equation is rigorously derived. Since the negative energy cannot bound many-body fluctuations, convergence is proven by approximating the dynamics via a quadratic Hamiltonian, in the spirit of Bogoliubov theory, and avoiding Grönwall estimates through the analysis of a coupled system of one-body linear equations for the one-particle density matrices.
In this thesis we study three examples of interacting many-body systems undergoing a non equilibrium time evolution. Firstly we consider the time evolution in an integrable system: the sine-Gordon field theory in the repulsive regime. We will focus on the one point function of the semi-local vertex operator eiβφ(x)/2 on a specific class of initial states. By analytical means we show that the expectation value considered decays exponentially to zero at late times and we determine the decay time. The method employed is based on a form-factor expansion and uses the "Representative Eigenstate Approach" of Ref. [73] (a.k.a. "Quench Action"). In a second example we study the time evolution in models close to "special" integrable points characterised by hidden symmetries generating infinitely many local conservation laws that do not commute with one another, in addition to the infinite commuting family implied by integrability. We observe that both in the case where the perturbation breaks the integrability and when it breaks only the additional symmetries maintaining integrability, the local observables show a crossover behaviour from an initial to a final quasi stationary plateau. We investigate a weak coupling limit, identify a time window in which the effects of the perturbations become significant and solve the time evolution through a mean-field mapping. As an explicit example we study the XYZ spin-1/2 chain with additional perturbations that break integrability. Finally, we study the effects of integrability breaking perturbations on the non-equilibrium evolution of more general many-particle quantum systems, where the unperturbed integrable model is generic. We focus on a class of spinless fermion models with weak interactions. We employ equation of motion techniques that can be viewed as generalisations of quantum Boltzmann equations. We benchmark our method against time dependent density matrix renormalisation group computations and find it to be very accurate as long as interactions are weak. For small integrability breaking, we observe robust prethermalisation plateaux for local observables on all accessible time scales. Increasing the strength of the integrability breaking term induces a "drift" away from the prethermalisation plateaux towards thermal behaviour. We identify a time scale characterising this crossover.