We have studied the dynamics of a generalized toric code based on qudits at finite temperature by finding the master equation coupling the code's degrees of freedom to a thermal bath. We find that in the case of qutrits, new types of anyons and thermal processes appear that are forbidden for qubits. These include creation, annihilation and diffusion throughout the system code. It is possible to solve the master equation in a short-time regime and find expressions for the decay rates as a function of the dimension d of the qudits. While we provide an explicit proof that the system relaxes to the Gibbs state for arbitrary qudits, we also prove that above a certain crossover temperature the qutrits' initial decay rate is smaller than the original case for qubits. Surprisingly, this behavior only happens for qutrits and not for other qudits with d > 3.
We have studied the dynamics of a generalized toric code based on qudits at finite temperature by finding the master equation coupling the code's degrees of freedom to a thermal bath. We find that in the case of qutrits, new types of anyons and thermal processes appear that are forbidden for qubits. These include creation, annihilation and diffusion throughout the system code. It is possible to solve the master equation in a short-time regime and find expressions for the decay rates as a function of the dimension d of the qudits. While we provide an explicit proof that the system relaxes to the Gibbs state for arbitrary qudits, we also prove that above a certain crossover temperature the qutrits' initial decay rate is smaller than the original case for qubits. Surprisingly, this behavior only happens for qutrits and not for other qudits with d > 3.
We present an analysis of the relaxation dynamics of finite-size topological qubits in contact with a thermal bath. Using a continuous-time Monte Carlo method, we explicitly compute the low-temperature nonequilibrium dynamics of the toric code on finite lattices. In contrast to the size-independent bound predicted for the toric code in the thermodynamic limit, we identify a low-temperature regime on finite lattices below a size-dependent crossover temperature with nontrivial finite-size and temperature scaling of the relaxation time. We demonstrate how this nontrivial finite-size scaling is governed by the scaling of topologically nontrivial two-dimensional classical random walks. The transition out of this low-temperature regime defines a dynamical finite-size crossover temperature that scales inversely with the log of the system size, in agreement with a crossover temperature defined from equilibrium properties. We find that both the finite-size and finite-temperature scaling are stronger in the low-temperature regime than above the crossover temperature. Since this finite-temperature scaling competes with the scaling of the robustness to unitary perturbations, this analysis may elucidate the scaling of memory lifetimes of possible physical realizations of topological qubits.
We address the out-of-equilibrium dynamics arising from quantum-quench (QQ) protocols (instantaneous changes of the Hamiltonian parameters) in many-body systems within their quantum critical regime and in contact with (homogeneously coupled) thermal baths. We consider two classes of QQ protocols. In one of them, the thermal bath is used to prepare the initial finite-temperature Gibbs state; then, after quenching, the thermal bath is removed and the dynamics of the system is unitary. We also address a more complex QQ protocol where the thermal bath is not removed after quenching, thus the quantum evolution is also driven by the interaction with the bath, which may be described by appropriate master equations for the density matrix of the system, where a further relevant time scale, or inverse decay rate, characterizes the system-bath coupling. Under these QQ protocols, the critical system develops out-of-equilibrium scaling behaviors, which extend those for isolated critical systems, by introducing further scaling variables proportional to the temperature of the thermal bath and the decay rate of the system-bath interactions. These out-of-equilibrium scaling behaviors are checked by analyzing QQ protocols within fermionic Kitaev wires, or equivalently quantum Ising chains, supplemented with a particular modelization of thermal bath that guarantees the asymptotic thermalization within the Lindblad master equation for the dynamics of open systems.
Providing the microscopic behavior of a thermalization process has always been an intriguing issue. There are several models of thermalization, which often require interaction of the system under consideration with the microscopic constituents of the macroscopic heat bath. With the aim to simulate such a thermalization process, here we look at the thermalization of a two-level quantum system under the action of a Markovian master equation corresponding to memory-less action of a heat bath, kept at a certain temperature, using a single-qubit ancilla. A two-qubit interaction Hamiltonian ( H th , say) is then designed —with a single-qubit thermal state as the initial state of the ancilla— which gives rise to thermalization of the system qubit in the infinite time limit. Further, we study the general form of Hamiltonian, of which ours is a special case, and look for the conditions for thermalization to occur. We also derive a Lindblad-type non-Markovian master equation for the system dynamics under the general form of system-ancilla Hamiltonian.
We study the nonequilibrium dynamics of the extended toric-code model (both ordered and disordered) to probe the existence of dynamical quantum phase transitions (DQPTs). We show that in the case of the ordered toric-code model, the zeros of Loschmidt overlap (generalized partition function) occur at critical times when DQPTs occur, which is confirmed by the nonanalyticities in the dynamical counterpart of the free-energy density. Moreover, we show that DQPTs occur for any nonzero field strength if the initial state is the excited state of the toric-code model. In the disordered case, we show that it is imperative to study the behavior of the first time derivative of the dynamical free-energy density averaged over all the possible configurations to characterize the occurrence of DQPTs in the disordered toric-code model since the disorder parameter itself acts as a new artificial dimension. We also show that for the case where anyonic excitations are present in the initial state, the conditions for DQPTs to occur are the same as what happens in the absence of any excitation.
We propose and study a model of a quantum memory that features self-correcting properties and a lifetime growing arbitrarily with system size at nonzero temperature. This is achieved by locally coupling a two-dimensional $L\ifmmode\times\else\texttimes\fi{}L$ toric code to a three-dimensional (3D) bath of bosons hopping on a cubic lattice. When the stabilizer operators of the toric code are coupled to the displacement operator of the bosons, we solve the model exactly via a polaron transformation and show that the energy penalty to create anyons grows linearly with L. When the stabilizer operators of the toric code are coupled to the bosonic density operator, we use perturbation theory to show that the energy penalty for anyons scales with ln(L). For a given error model, these energy penalties lead to a lifetime of the stored quantum information growing, respectively, exponentially and polynomially with L. Furthermore, we show how to choose an appropriate coupling scheme in order to hinder the hopping of anyons (and not only their creation) with energy barriers that are of the same order as the anyon creation gaps. We argue that a toric code coupled to a 3D Heisenberg ferromagnet realizes our model in its low-energy sector. Finally, we discuss the delicate issue of the stability of topological order in the presence of perturbations. While we do not derive a rigorous proof of topological order, we present heuristic arguments suggesting that topological order remains intact when perturbative operators acting on the toric code spins are coupled to the bosonic environment.
We analyze quantum state preservation in open quantum systems using quantum error-correcting (QEC) codes explicitly embedded in microscopic system-bath models. Rather than assuming abstract quantum channels, we consider multi-qubit registers coupled to bosonic thermal environments, derive a second-order master equation for the reduced dynamics, and use it to benchmark the five-qubit, Steane, and toric codes under local and collective noise. We compute state fidelities as functions of system-bath coupling strength, bath temperatures, and the number of correction cycles. In the low-temperature regime, repeated error correction with the five-qubit code significantly suppresses decoherence and relaxation for weak-to-moderate couplings. In the high-temperature regime, thermal excitations reduce the effectiveness of all codes, although within the parameter range studied, the five-qubit code still yields the highest fidelities among the three codes. For two-qubit Werner states, we identify a critical evolution time associated with an early-time crossover, before which the overhead of QEC does not compensate for the noise-induced degradation; this critical time increases with entanglement, reflecting the greater fragility of strongly entangled states. Overall, our results provide a microscopic master-equation-based framework for benchmarking QEC performance in realistic open-system environments and for assessing code behavior in near-term noisy quantum architectures.