Research output: Contribution to journal › Article

**Exactly Thermalised Quantum Dynamics of the Spin-Boson Model coupled to a Dissipative Environment.** / Lane, Matt; Matos, Daniel; Ford, Ian J.; Kantorovich, Lev.

Research output: Contribution to journal › Article

Lane, M, Matos, D, Ford, IJ & Kantorovich, L 2020, 'Exactly Thermalised Quantum Dynamics of the Spin-Boson Model coupled to a Dissipative Environment', *Physical Review B*.

Lane, M., Matos, D., Ford, I. J., & Kantorovich, L. (Accepted/In press). Exactly Thermalised Quantum Dynamics of the Spin-Boson Model coupled to a Dissipative Environment. *Physical Review B*.

Lane M, Matos D, Ford IJ, Kantorovich L. Exactly Thermalised Quantum Dynamics of the Spin-Boson Model coupled to a Dissipative Environment. Physical Review B. 2020 Apr 8.

@article{fbaf7ac900614b3aad396f6e5fea3d7d,

title = "Exactly Thermalised Quantum Dynamics of the Spin-Boson Model coupled to a Dissipative Environment",

abstract = "We present an application of the Extended Stochastic Liouville-von Neumann equations (ESLN) method introduced earlier [PRB 95, 125124 (2017); PRB 97, 224310 (2018)] which describes the dynamics of an exactly thermalised open quantum system reduced density matrix coupled to a non- Markovian harmonic environment. Critically, the combined system of the open system fully coupled to its environment is thermalised at finite temperature using an imaginary time evolution procedure before the application of real time evolution. This initialises the combined system in the correct canonical equilibrium state rather than being initially decoupled. Here we apply our theory to the spin-boson Hamiltonian and develop a number of competing ESLN variants designed to reduce the numerical divergence of the trace of the open system density matrix. We find that a careful choice of the driving noises is essential for improving numerical stability. We also investigate the effect of applying higher order numerical schemes for solving stochastic differential equations, such as the Stratonovich-Heun scheme, and conclude that stochastic sampling dominates convergence with the improvement associated with the numerical scheme being less important for short times but required for late times. To verify the method and its numerical implementation, we first consider evolution under a fixed Hamiltonian and show that the system either remains in, or approaches, the correct canonical equilibrium state at long times. Additionally, evolution of the open system under non-equilibrium Landau-Zener (LZ) driving is considered and the asymptotic convergence to the LZ limit was observed for vanishing system-environment coupling and temperature. When coupling and temperature are non-zero, initially thermalising the combined system at a finite time in the past was found to be a better approximation of the true LZ initial state than starting in a pure state.",

author = "Matt Lane and Daniel Matos and Ford, {Ian J.} and Lev Kantorovich",

year = "2020",

month = apr,

day = "8",

language = "English",

journal = "Physical Review B (Condensed Matter and Materials Physics)",

issn = "1098-0121",

publisher = "American Physical Society",

}

TY - JOUR

T1 - Exactly Thermalised Quantum Dynamics of the Spin-Boson Model coupled to a Dissipative Environment

AU - Lane, Matt

AU - Matos, Daniel

AU - Ford, Ian J.

AU - Kantorovich, Lev

PY - 2020/4/8

Y1 - 2020/4/8

N2 - We present an application of the Extended Stochastic Liouville-von Neumann equations (ESLN) method introduced earlier [PRB 95, 125124 (2017); PRB 97, 224310 (2018)] which describes the dynamics of an exactly thermalised open quantum system reduced density matrix coupled to a non- Markovian harmonic environment. Critically, the combined system of the open system fully coupled to its environment is thermalised at finite temperature using an imaginary time evolution procedure before the application of real time evolution. This initialises the combined system in the correct canonical equilibrium state rather than being initially decoupled. Here we apply our theory to the spin-boson Hamiltonian and develop a number of competing ESLN variants designed to reduce the numerical divergence of the trace of the open system density matrix. We find that a careful choice of the driving noises is essential for improving numerical stability. We also investigate the effect of applying higher order numerical schemes for solving stochastic differential equations, such as the Stratonovich-Heun scheme, and conclude that stochastic sampling dominates convergence with the improvement associated with the numerical scheme being less important for short times but required for late times. To verify the method and its numerical implementation, we first consider evolution under a fixed Hamiltonian and show that the system either remains in, or approaches, the correct canonical equilibrium state at long times. Additionally, evolution of the open system under non-equilibrium Landau-Zener (LZ) driving is considered and the asymptotic convergence to the LZ limit was observed for vanishing system-environment coupling and temperature. When coupling and temperature are non-zero, initially thermalising the combined system at a finite time in the past was found to be a better approximation of the true LZ initial state than starting in a pure state.

AB - We present an application of the Extended Stochastic Liouville-von Neumann equations (ESLN) method introduced earlier [PRB 95, 125124 (2017); PRB 97, 224310 (2018)] which describes the dynamics of an exactly thermalised open quantum system reduced density matrix coupled to a non- Markovian harmonic environment. Critically, the combined system of the open system fully coupled to its environment is thermalised at finite temperature using an imaginary time evolution procedure before the application of real time evolution. This initialises the combined system in the correct canonical equilibrium state rather than being initially decoupled. Here we apply our theory to the spin-boson Hamiltonian and develop a number of competing ESLN variants designed to reduce the numerical divergence of the trace of the open system density matrix. We find that a careful choice of the driving noises is essential for improving numerical stability. We also investigate the effect of applying higher order numerical schemes for solving stochastic differential equations, such as the Stratonovich-Heun scheme, and conclude that stochastic sampling dominates convergence with the improvement associated with the numerical scheme being less important for short times but required for late times. To verify the method and its numerical implementation, we first consider evolution under a fixed Hamiltonian and show that the system either remains in, or approaches, the correct canonical equilibrium state at long times. Additionally, evolution of the open system under non-equilibrium Landau-Zener (LZ) driving is considered and the asymptotic convergence to the LZ limit was observed for vanishing system-environment coupling and temperature. When coupling and temperature are non-zero, initially thermalising the combined system at a finite time in the past was found to be a better approximation of the true LZ initial state than starting in a pure state.

M3 - Article

JO - Physical Review B (Condensed Matter and Materials Physics)

JF - Physical Review B (Condensed Matter and Materials Physics)

SN - 1098-0121

ER -

King's College London - Homepage

© 2018 King's College London | Strand | London WC2R 2LS | England | United Kingdom | Tel +44 (0)20 7836 5454