Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2011.11059 · Journal of Physics A Mathematical and Theoretical

Non-equilibrium dynamics of a dissipative two-site Hubbard model simulated on IBM quantum computers

Sabine Tornow, Wolfgang Gehrke, U. Helmbrecht

Published Apr 29, 2022• 15 citations
arXivHTMLPDF

Abstract

Abstract Many-body physics is one very well suited field for testing quantum algorithms and for finding working heuristics on present quantum computers. We have investigated the non-equilibrium dynamics of one- and two-electron systems, which are coupled to an environment that introduces decoherence and dissipation. In our approach, the electronic system is represented in the framework of a two-site Hubbard model while the environment is modelled by a spin bath. To simulate the non-equilibrium population probabilities of the different states on a quantum computer we have encoded the electronic states and environmental degrees of freedom into qubits and ancilla qubits (bath), respectively. The total evolution time was divided into short time intervals, during which the system evolves. After each of these time steps, the system interacts with ancilla qubits representing the bath in thermal equilibrium. We have specifically studied spin baths leading to both, unital and non-unital dynamics of the electronic system and have found that electron correlations clearly enhance the electron transfer rates in the latter case. For short time periods, the simulation on the quantum computer is found to be in very good agreement with the exact results if error mitigation methods are applied. Our method to simulate also non-unitary time-evolution on a quantum computer can be well extended to simulate electronic systems in correlated spin baths as well as in bosonic and fermionic baths.

DOI
10.1088/1751-8121/ac6bd0

Related papers

Non-Equilibrium Dynamics of a Dissipative Two-Site Hubbard Model Simulated on the IBM Quantum Computer

Sabine Tornow, Wolfgang Gehrke, U. Helmbrecht

Many-body physics is one very well suited field for testing quantum algorithms and for finding working heuristics on present quantum computers. We have investigated the non-equilibrium dynamics of one- and two-electron systems, which are coupled to an environment that introduces decoherence and dissipation. In our approach, the electronic system is represented in the framework of a two-site Hubbard model while the environment is modelled by a spin bath. To simulate the non-equilibrium population probabilities of the different states on a quantum computer we have encoded the electronic states and environmental degrees of freedom into qubits and ancilla qubits (bath), respectively. The total evolution time was divided into short time intervals, during which the system evolves. After each of these time steps, the system interacts with ancilla qubits representing the bath in thermal equilibrium. We have specifically studied spin baths leading to both, unital and non-unital dynamics of the electronic system and have found that electron correlations clearly enhance the electron transfer rates in the latter case. For short time periods, the simulation on the quantum computer is found to be in very good agreement with the exact results if error mitigation methods are applied. Our method to simulate also non-unitary time-evolution on a quantum computer can be well extended to simulate electronic systems in correlated spin baths as well as in bosonic and fermionic baths.

•Nov 22, 2020•3 citations

Dissipative two-electron transfer: A numerical renormalization group study

Sabine Tornow, R. Bulla, Frithjof B. Anders +1 more

We investigate nonequilibrium two-electron transfer in a model redox system represented by a two-site extended Hubbard model and embedded in a dissipative environment. The influence of the electron-electron interactions and the coupling to a dissipative bosonic bath on the electron transfer is studied in different temperature regimes. At high temperatures, Marcus transfer rates are evaluated, and at low temperatures, we calculate equilibrium and nonequilibrium population probabilities of the donor and acceptor with the nonperturbative numerical renormalization group approach. We obtain the nonequilibrium dynamics of the system prepared in an initial state of two electrons at the donor site and identify conditions under which the electron transfer involves one concerted two-electron step or two sequential single-electron steps. The rates of the sequential transfer depend nonmonotonically on the difference between the intersite and on-site Coulomb interaction, which become renormalized in the presence of the bosonic bath. If this difference is much larger than the hopping matrix element, the temperature as well as the reorganization energy, simultaneous transfer of both electrons between donor and acceptor can be observed.

arXiv:0803.4066•Physical Review B•Jul 18, 2008•24 citations

Driven-dissipative quantum mechanics on a lattice: Simulating a fermionic reservoir on a quantum computer

Lorenzo Del Re, Brian Rost, Alexander F. Kemper +1 more

The driven-dissipative many-body problem remains one of the most challenging unsolved problems in quantum mechanics. The advent of quantum computers may provide a unique platform for efficiently simulating such driven-dissipative systems. But, there are many choices for how one can engineer the reservoir. One can simply employ ancilla qubits to act as a reservoir and then digitally simulate them via algorithmic cooling. A more attractive approach, which allows one to simulate an infinite reservoir, is to integrate out the bath degrees of freedom and describe the driven-dissipative system via a master equation, that can also be simulated on a quantum computer. In this work, we consider the particular case of noninteracting electrons on a lattice driven by an electric field and coupled to a fermionic thermostat. Then, we provide two different quantum circuits: the first one reconstructs the full dynamics of the system using Trotter steps, while the second one dissipatively prepares the final nonequilibrium steady state in a single step. We run both circuits on the IBM quantum experience. For circuit (i), we achieved up to five Trotter steps. When partial resets become available on quantum computers, we expect that the maximum simulation time can be significantly increased. The methods developed here suggest generalizations that can be applied to simulating interacting driven-dissipative systems.

arXiv:1912.08310•Physical review. B./Physical review. B•Sep 10, 2020•67 citations

Long-time error-mitigating simulation of open quantum systems on near term quantum computers

Brian Rost, Lorenzo Del Re, Nathan Earnest +3 more

We study an open quantum system simulation on quantum hardware, which demonstrates robustness to hardware errors even with deep circuits containing up to two thousand entangling gates. We simulate two systems of electrons coupled to an infinite thermal bath: 1) a system of dissipative free electrons in a driving electric field; and 2) the thermalization of two interacting electrons in a single orbital in a magnetic field—the Hubbard atom. These problems are solved using IBM quantum computers, showing no signs of decreasing fidelity at long times. Our results demonstrate that algorithms for simulating open quantum systems are able to far outperform similarly complex non-dissipative algorithms on noisy hardware. Our two examples show promise that the driven-dissipative quantum many-body problem can eventually be solved on quantum computers.

arXiv:2108.01183•npj Quantum Information•Jan 21, 2025•20 citations

Nonequilibrium many-body dynamics along a dissipative Hubbard chain: Symmetries and quantum Monte Carlo simulations

Lothar Mühlbacher, Joachim Ankerhold

The nonequilibrium dynamics of correlated charge transfer along a one-dimensional chain in the presence of a phonon environment is investigated within a dissipative Hubbard model. For this generalization of the ubiquitous spin-boson model the crucial role of symmetries is analyzed in detail and corresponding invariant subspaces are identified. It is shown that the time evolution typically occurs in each of the disjunct subspaces independently leading, e.g., asymptotically to a non-Boltzmann equilibrium state. Based on these findings explicit results are obtained for two interacting electrons by means of a substantially improved real-time quantum Monte Carlo approach. In the incoherent regime an appropriate mapping of the many-body dynamics onto an isomorphic single particle motion allows for an approximate description of the numerical data in terms of rate equations. These results may lead to control schemes of charge transport in tailored quantum systems as, e.g., molecular chains or quantum dot arrays.

arXiv:cond-mat/0605166•Physical Review B•Oct 12, 2006•3 citations

Lessons on electronic decoherence in molecules from exact modeling

Wenxiang Hu, Bing Gu, Ignacio Franco

Electronic decoherence processes in molecules and materials are usually thought and modeled via schemes for the system-bath evolution in which the bath is treated either implicitly or approximately. Here we present computations of the electronic decoherence dynamics of a model many-body molecular system described by the Su-Schrieffer-Heeger Hamiltonian with Hubbard electron-electron interactions using an exact method in which both electronic and nuclear degrees of freedom are taken into account explicitly and fully quantum mechanically. To represent the electron-nuclear Hamiltonian in matrix form and propagate the dynamics, the computations employ the Jordan-Wigner transformation for the fermionic creation/annihilation operators and the discrete variable representation for the nuclear operators. The simulations offer a standard for electronic decoherence that can be used to test approximations. They also provide a useful platform to answer fundamental questions about electronic decoherence that cannot be addressed through approximate or implicit schemes. Specifically, through simulations, we isolate basic mechanisms for electronic coherence loss and demonstrate that electronic decoherence is possible even for one-dimensional nuclear bath. Furthermore, we show that (i) decreasing the mass of the bath generally leads to faster electronic decoherence; (ii) electron-electron interactions strongly affect the electronic decoherence when the electron-nuclear dynamics is not pure-dephasing; (iii) classical bath models with initial conditions sampled from the Wigner distribution accurately capture the short-time electronic decoherence dynamics; (iv) model separable initial superpositions often used to understand decoherence after photoexcitation are only relevant in experiments that employ delta-like laser pulses to initiate the dynamics. These insights can be employed to interpret and properly model coherence phenomena in molecules.

arXiv:1801.05846•The Journal of Chemical Physics•Apr 3, 2018•35 citations
Abstract Many-body physics is one very well suited field for testing quantum algorithms and for finding working heuristics on present quantum computers. We have investigated the non-equilibrium dynamics of one- and two-electron systems, which are coupled to an environment that introduces decoherence and dissipation. In our approach, the electronic system is represented in the framework of a two-site Hubbard model while the environment is modelled by a spin bath. To simulate the non-equilibrium population probabilities of the different states on a quantum computer we have encoded the electronic states and environmental degrees of freedom into qubits and ancilla qubits (bath), respectively. The total evolution time was divided into short time intervals, during which the system evolves. After each of these time steps, the system interacts with ancilla qubits representing the bath in thermal equilibrium. We have specifically studied spin baths leading to both, unital and non-unital dynamics of the electronic system and have found that electron correlations clearly enhance the electron transfer rates in the latter case. For short time periods, the simulation on the quantum computer is found to be in very good agreement with the exact results if error mitigation methods are applied. Our method to simulate also non-unitary time-evolution on a quantum computer can be well extended to simulate electronic systems in correlated spin baths as well as in bosonic and fermionic baths.
Many-body physics is one very well suited field for testing quantum algorithms and for finding working heuristics on present quantum computers. We have investigated the non-equilibrium dynamics of one- and two-electron systems, which are coupled to an environment that introduces decoherence and dissipation. In our approach, the electronic system is represented in the framework of a two-site Hubbard model while the environment is modelled by a spin bath. To simulate the non-equilibrium population probabilities of the different states on a quantum computer we have encoded the electronic states and environmental degrees of freedom into qubits and ancilla qubits (bath), respectively. The total evolution time was divided into short time intervals, during which the system evolves. After each of these time steps, the system interacts with ancilla qubits representing the bath in thermal equilibrium. We have specifically studied spin baths leading to both, unital and non-unital dynamics of the electronic system and have found that electron correlations clearly enhance the electron transfer rates in the latter case. For short time periods, the simulation on the quantum computer is found to be in very good agreement with the exact results if error mitigation methods are applied. Our method to simulate also non-unitary time-evolution on a quantum computer can be well extended to simulate electronic systems in correlated spin baths as well as in bosonic and fermionic baths.
We investigate nonequilibrium two-electron transfer in a model redox system represented by a two-site extended Hubbard model and embedded in a dissipative environment. The influence of the electron-electron interactions and the coupling to a dissipative bosonic bath on the electron transfer is studied in different temperature regimes. At high temperatures, Marcus transfer rates are evaluated, and at low temperatures, we calculate equilibrium and nonequilibrium population probabilities of the donor and acceptor with the nonperturbative numerical renormalization group approach. We obtain the nonequilibrium dynamics of the system prepared in an initial state of two electrons at the donor site and identify conditions under which the electron transfer involves one concerted two-electron step or two sequential single-electron steps. The rates of the sequential transfer depend nonmonotonically on the difference between the intersite and on-site Coulomb interaction, which become renormalized in the presence of the bosonic bath. If this difference is much larger than the hopping matrix element, the temperature as well as the reorganization energy, simultaneous transfer of both electrons between donor and acceptor can be observed.
The driven-dissipative many-body problem remains one of the most challenging unsolved problems in quantum mechanics. The advent of quantum computers may provide a unique platform for efficiently simulating such driven-dissipative systems. But, there are many choices for how one can engineer the reservoir. One can simply employ ancilla qubits to act as a reservoir and then digitally simulate them via algorithmic cooling. A more attractive approach, which allows one to simulate an infinite reservoir, is to integrate out the bath degrees of freedom and describe the driven-dissipative system via a master equation, that can also be simulated on a quantum computer. In this work, we consider the particular case of noninteracting electrons on a lattice driven by an electric field and coupled to a fermionic thermostat. Then, we provide two different quantum circuits: the first one reconstructs the full dynamics of the system using Trotter steps, while the second one dissipatively prepares the final nonequilibrium steady state in a single step. We run both circuits on the IBM quantum experience. For circuit (i), we achieved up to five Trotter steps. When partial resets become available on quantum computers, we expect that the maximum simulation time can be significantly increased. The methods developed here suggest generalizations that can be applied to simulating interacting driven-dissipative systems.
We study an open quantum system simulation on quantum hardware, which demonstrates robustness to hardware errors even with deep circuits containing up to two thousand entangling gates. We simulate two systems of electrons coupled to an infinite thermal bath: 1) a system of dissipative free electrons in a driving electric field; and 2) the thermalization of two interacting electrons in a single orbital in a magnetic field—the Hubbard atom. These problems are solved using IBM quantum computers, showing no signs of decreasing fidelity at long times. Our results demonstrate that algorithms for simulating open quantum systems are able to far outperform similarly complex non-dissipative algorithms on noisy hardware. Our two examples show promise that the driven-dissipative quantum many-body problem can eventually be solved on quantum computers.
The nonequilibrium dynamics of correlated charge transfer along a one-dimensional chain in the presence of a phonon environment is investigated within a dissipative Hubbard model. For this generalization of the ubiquitous spin-boson model the crucial role of symmetries is analyzed in detail and corresponding invariant subspaces are identified. It is shown that the time evolution typically occurs in each of the disjunct subspaces independently leading, e.g., asymptotically to a non-Boltzmann equilibrium state. Based on these findings explicit results are obtained for two interacting electrons by means of a substantially improved real-time quantum Monte Carlo approach. In the incoherent regime an appropriate mapping of the many-body dynamics onto an isomorphic single particle motion allows for an approximate description of the numerical data in terms of rate equations. These results may lead to control schemes of charge transport in tailored quantum systems as, e.g., molecular chains or quantum dot arrays.
Electronic decoherence processes in molecules and materials are usually thought and modeled via schemes for the system-bath evolution in which the bath is treated either implicitly or approximately. Here we present computations of the electronic decoherence dynamics of a model many-body molecular system described by the Su-Schrieffer-Heeger Hamiltonian with Hubbard electron-electron interactions using an exact method in which both electronic and nuclear degrees of freedom are taken into account explicitly and fully quantum mechanically. To represent the electron-nuclear Hamiltonian in matrix form and propagate the dynamics, the computations employ the Jordan-Wigner transformation for the fermionic creation/annihilation operators and the discrete variable representation for the nuclear operators. The simulations offer a standard for electronic decoherence that can be used to test approximations. They also provide a useful platform to answer fundamental questions about electronic decoherence that cannot be addressed through approximate or implicit schemes. Specifically, through simulations, we isolate basic mechanisms for electronic coherence loss and demonstrate that electronic decoherence is possible even for one-dimensional nuclear bath. Furthermore, we show that (i) decreasing the mass of the bath generally leads to faster electronic decoherence; (ii) electron-electron interactions strongly affect the electronic decoherence when the electron-nuclear dynamics is not pure-dephasing; (iii) classical bath models with initial conditions sampled from the Wigner distribution accurately capture the short-time electronic decoherence dynamics; (iv) model separable initial superpositions often used to understand decoherence after photoexcitation are only relevant in experiments that employ delta-like laser pulses to initiate the dynamics. These insights can be employed to interpret and properly model coherence phenomena in molecules.