Table of Contents
Fetching ...

On the Noisy Road to Open Quantum Dynamics: The Place of Stochastic Hamiltonians

Pietro De Checchi, Federico Gallina, Barbara Fresch, Giulio G. Giusteri

TL;DR

The paper provides a unified, self-contained analysis of stochastic Hamiltonians, the stochastic Schrödinger equation, and the stochastic Liouville equation for open quantum systems. It contrasts Itō and Stratonovich calculi to show how unitary, trajectory-level simulations and Lindblad-type master equations emerge, including exact and approximate links to Redfield and other dissipative models. By detailing continuous versus white-noise processes and presenting density-matrix as well as wavefunction formalisms, the work offers actionable recipes for classical trajectory integration and quantum-algorithm implementations. The results illuminate when stochastic Hamiltonians map to QMEs, how to implement random-unitary gates on quantum hardware, and how to leverage the SLE framework to incorporate dissipation and measurements. Overall, the paper furnishes a cohesive toolkit for modeling and simulating open quantum dynamics across classical and quantum computational platforms, with guidance on framework selection based on physical insight and computational context.

Abstract

Stochastic evolution underpins several approaches to the dynamics of open quantum systems, such as random modulation of Hamiltonian parameters, the stochastic Schrödinger equation (SSE), and the stochastic Liouville equation (SLE). In a stochastic formulation, the open-system problem is reduced from a coupled system-environment dynamics to an effective system-only description, with dissipative evolution recovered by ensemble averaging. In this work, we aim at a self-contained and accessible presentation of these approaches to further elaborate on their common roots in essential concepts of stochastic calculus and to delineate the conditions under which they are equivalent. We also discuss how different formulations naturally lead to different numerical time-integration schemes, better suited for either classical simulation platforms, based on finite-difference approximations, or quantum algorithms, that employ random unitary maps. Our analysis supplies a unified perspective and actionable recipes for classical and quantum implementations of stochastic evolution in the simulation of open quantum systems.

On the Noisy Road to Open Quantum Dynamics: The Place of Stochastic Hamiltonians

TL;DR

The paper provides a unified, self-contained analysis of stochastic Hamiltonians, the stochastic Schrödinger equation, and the stochastic Liouville equation for open quantum systems. It contrasts Itō and Stratonovich calculi to show how unitary, trajectory-level simulations and Lindblad-type master equations emerge, including exact and approximate links to Redfield and other dissipative models. By detailing continuous versus white-noise processes and presenting density-matrix as well as wavefunction formalisms, the work offers actionable recipes for classical trajectory integration and quantum-algorithm implementations. The results illuminate when stochastic Hamiltonians map to QMEs, how to implement random-unitary gates on quantum hardware, and how to leverage the SLE framework to incorporate dissipation and measurements. Overall, the paper furnishes a cohesive toolkit for modeling and simulating open quantum dynamics across classical and quantum computational platforms, with guidance on framework selection based on physical insight and computational context.

Abstract

Stochastic evolution underpins several approaches to the dynamics of open quantum systems, such as random modulation of Hamiltonian parameters, the stochastic Schrödinger equation (SSE), and the stochastic Liouville equation (SLE). In a stochastic formulation, the open-system problem is reduced from a coupled system-environment dynamics to an effective system-only description, with dissipative evolution recovered by ensemble averaging. In this work, we aim at a self-contained and accessible presentation of these approaches to further elaborate on their common roots in essential concepts of stochastic calculus and to delineate the conditions under which they are equivalent. We also discuss how different formulations naturally lead to different numerical time-integration schemes, better suited for either classical simulation platforms, based on finite-difference approximations, or quantum algorithms, that employ random unitary maps. Our analysis supplies a unified perspective and actionable recipes for classical and quantum implementations of stochastic evolution in the simulation of open quantum systems.

Paper Structure

This paper contains 14 sections, 55 equations, 3 figures.

Figures (3)

  • Figure 1: A schematic map of the two mathematical frameworks and relative workflows for the stochastic Hamiltonian applied to a wavefunction dynamics, based on the Itō (upper green branch) and Stratonovich (lower blue branch) interpretations. The Stratonovich pathway connects stochastic Hamiltonians to direct simulations of the density matrix dynamics through random unitary gates, while the Itō pathway provides the correct formulation for analytically obtaining quantum master equations. The random ordinary differential equation (RODE) case appears as a special limit when dealing with continuous stochastic processes as Hamiltonian fluctuations.
  • Figure 2: An example of a quantum circuit for the propagation of stochastic dynamics. Single qubits gates implement the unitary operators $R_\alpha(2 \theta) = \exp(-i \theta \sigma_\alpha)$, with $\alpha=\{Z,X\}$, where $\sigma_\alpha$ are Pauli matrices. While parameters $\theta_\varepsilon$ and $\theta_\Omega$ are deterministic, $\theta_\xi$ is intended as a stochastic term changing at each time step and each repetition of the circuit. The Magnus-Trotter block evolving the system is performed $n$ times to evolve the system up to the measurement.
  • Figure 3: A schematic map of the stochastic Liouville equation (SLE) framework, where averaging the commutator dynamics yields the exact evolution of the density matrix (upper branch), which can be solved numerically, while approximate models provide analytical solutions (lower branch).