Abstract We simulate coherent driven free dissipative Kerr nonlinear system numerically using Euler’s method by solving Heisenberg equation of motion and time evolving block decimation (TEBD) algorithm, and demonstrate how the numerical results are analogous to classical bistability. The comparison with analytics show that the TEBD numerics follow the quantum mechanical exact solution obtained by mapping the equation of motion of the density matrix of the system to a Fokker-Plank equation . Comparing between two different numerical techniques, we see that the semi-classical Euler’s method gives the dynamics of the system field of one among two coherent branches, whereas TEBD numerics generate the superposition of both of them. Therefore, the time dynamics determined by TEBD numerical method undergoes through a non-classical state which is also shown by determining second order correlation function.
Abstract We simulate coherent driven free dissipative Kerr nonlinear system numerically using Euler’s method by solving Heisenberg equation of motion and time evolving block decimation (TEBD) algorithm, and demonstrate how the numerical results are analogous to classical bistability. The comparison with analytics show that the TEBD numerics follow the quantum mechanical exact solution obtained by mapping the equation of motion of the density matrix of the system to a Fokker-Plank equation . Comparing between two different numerical techniques, we see that the semi-classical Euler’s method gives the dynamics of the system field of one among two coherent branches, whereas TEBD numerics generate the superposition of both of them. Therefore, the time dynamics determined by TEBD numerical method undergoes through a non-classical state which is also shown by determining second order correlation function.
Abstract We simulate coherent driven free dissipative Kerr nonlinear system numerically using time-evolving block decimation (TEBD) algorithm and time propagation on the Heisenberg equation of motion using Euler’s method to study how the numerical results are analogous to classical bistability . The system evolves through different trajectories to stabilize different branches for different external drives and initial conditions. The Wigner state reprentation confirms the system to suffer a residual effect of initial state throughout the non-classical dynamical evolution and the metastable states of the system . Furthermore, we also see the numerically simulated spectral density remains significantly different from analytical counterparts when initial states do not lie to the same branch of the final state.
We simulate a free dissipative and coherent-driven Kerr nonlinear system using a time-evolving block decimation (TEBD) algorithm to study the impact of the initial state on the exact quantum dynamics of the system. The superposition of two coherent branches results in non-classical time dynamics. The Wigner state representation confirms that the system ends up saturating to two different branches, through evolving different trajectories, resulting in de-Gaussification throughout evolution. Furthermore, we also see that the time evolution suffers the residual effect of the initial state.
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 generation of squeezed states of light in a two-mode Kerr nonlinear directional coupler (NLDC) was investigated using two different methods in quantum mechanics. First, the analytical method, a Heisenberg-picture-based method where the operators are evolving in time but the state vectors are time-independent. In this method, an analytical solution to the coupled Heisenberg equations of motion for the propagating modes was proposed based on the Baker–Hausdorff (BH) formula. Second, the phase space method, a Schrödinger-picture-based method in which the operators are constant and the density matrix evolves in time. In this method, the quantum mechanical master equation of the density matrix was converted to a corresponding classical Fokker–Planck (FP) equation in positive-P representation. Then, the FP equation was converted to a set of stochastic differential equations using Ito rules. The strengths and weaknesses of each method are discussed. Good agreement between both methods was achieved, especially at early evolution stages and lower values of linear coupling coefficient. On one hand, the analytical method seems insensitive to higher values of nonlinear coupling coefficients. Nevertheless, it demonstrated better numerical stability. On the other hand, the solution of the stochastic equations resulting from the phase space method is numerically expensive as it requires averaging over thousands of trajectories. Besides, numerically unstable trajectories appear with positive-P representation at higher values of nonlinearity.
Driven-dissipative systems have recently attracted great attention due to the existence of novel physical phenomena with no analog in the equilibrium case. The Keldysh path-integral theory is a powerful tool to investigate these systems. However, it has still been a challenge to study strong nonlinear effects implemented by recent experiments, since in this case the photon number is few and quantum fluctuations play a crucial role in the dynamics of the system. Here we develop an approach for deriving exact steady states of driven-dissipative systems by introducing the Keldysh partition function in the Fock-state basis and then mapping the standard saddle-point equations into Keldysh-Heisenberg equations. We take the strong Kerr nonlinear resonators with and without the nonlinear driving as two examples to illustrate our method. It is found that, in the absence of the nonlinear driving, the exact steady state obtained does not exhibit bistability and agrees well with the complex P-representation solution. While in the presence of the nonlinear driving, the multiphoton resonance effects are revealed and are consistent with the qualitative analysis. Our method provides an intuitive way to explore a variety of driven-dissipative systems especially with strong correlations.
For quantum computers to become useful tools to physicists, engineers and computational scientists, quantum algorithms for solving nonlinear differential equations need to be developed. Despite recent advances, the quest for a solver that can integrate nonlinear dynamical systems with a quantum advantage, while being realizable on available (or near-term) quantum hardware, is an open challenge. In this paper, we transform a nonlinear dynamical system into a linear system, which we integrate with quantum algorithms. Key to the method is the Fokker–Planck equation, which is a non-normal partial differential equation. Three approaches are proposed: (i) forward-Euler stepping by block encoding; (ii) Schrödingerization and (iii) forward-Euler stepping by a linear combination of unitaries. We emulate the integration of prototypical nonlinear systems with the proposed quantum solvers and compare the output with the benchmark solutions of classical integrators. We find that classical and quantum outputs are in good agreement. This paper opens opportunities for solving nonlinear differential equations with quantum algorithms.