We present a framework for the realization of dissipative evolutions of spin-boson models, including multiphoton exchange dynamics, as well as nonlinear transition rates. Our approach is based on the implementation of a generalized version of a dissipative linear quantum Rabi model. The latter comprises a linearly coupled spin-boson term, spin rotations, and standard dissipators. We provide numerical simulations of illustrative cases supporting the good performance of our method. Our work allows for the simulation of a large class of fundamentally different quantum models where the effect of distinct dissipative processes can be easily investigated.
We present a framework for the realization of dissipative evolutions of spin-boson models, including multiphoton exchange dynamics, as well as nonlinear transition rates. Our approach is based on the implementation of a generalized version of a dissipative linear quantum Rabi model. The latter comprises a linearly coupled spin-boson term, spin rotations, and standard dissipators. We provide numerical simulations of illustrative cases supporting the good performance of our method. Our work allows for the simulation of a large class of fundamentally different quantum models where the effect of distinct dissipative processes can be easily investigated.
The paradigmatic spin-boson model considers a spin degree of freedom interacting with an environment typically constituted by a continuum of bosonic modes. This ubiquitous model is of relevance in a number of physical systems where, in general, one has neither control over the bosonic modes, nor the ability to tune distinct interaction mechanisms. Despite this apparent lack of control, we present a suitable transformation that approximately maps the spin-boson dynamics into that of a tunable multiphoton Jaynes-Cummings model undergoing dissipation. Interestingly, the latter model describes the coherent interaction between a spin and a single bosonic mode via the simultaneous exchange of n bosons per spin excitation. Resorting to the so-called reaction coordinate method, we identify a relevant collective bosonic mode in the environment, which is then used to generate multiphoton interactions following the proposed theoretical framework. Moreover, we show that spin-boson models featuring structured environments can lead to non-Markovian multiphoton Jaynes-Cummings dynamics. We discuss the validity of the proposed method depending on the parameters and analyse its performance, which is supported by numerical simulations. In this manner, the spin-boson model serves as a good analogue quantum simulator for the inspection and realization of multiphoton Jaynes-Cummings models, as well as the interplay of non-Markovian effects and, thus, as a simulator of light-matter systems with tunable interaction mechanisms.
Abstract We establish an approximate equivalence between a generalised quantum Rabi model and its nth order counterparts, where spin-boson interactions are nonlinear as they comprise a simultaneous exchange of n bosonic excitations. Although there exists no unitary transformation between these models, we demonstrate their equivalence to a good approximation in a wide range of parameters. This shows that nonlinear spin-boson couplings, i.e., nth order quantum Rabi models, are accessible to quantum systems with only linear coupling between boson and spin modes by simply adding spin rotations and after an appropriate transformation. Furthermore, our result prompts novel approximate analytical solutions to the dynamics of the quantum Rabi model in the ultrastrong coupling regime improving previous approaches.
We establish an approximate equivalence between a generalised quantum Rabi model and its nth order counterparts, where spinboson interactions are nonlinear as they comprise a simultaneous exchange of n bosonic excitations. Although there exists no unitary transformation between these models, we demonstrate their equivalence to a good approximation in a wide range of parameters. This shows that nonlinear spin-boson couplings, i.e., nth order quantum Rabi models, are accessible to quantum systems with only linear coupling between boson and spin modes by simply adding spin rotations and after an appropriate transformation. Furthermore, our result prompts novel approximate analytical solutions to the dynamics of the quantum Rabi model in the ultrastrong coupling regime improving previous approaches.
Spin-boson models involving many interacting spins and bosons are ubiquitous in quantum simulation platforms. At the same time, characterizing the dynamics of these quantum systems represents a significant challenge. Here, we consider general spin-boson models where bosons are subject to Markovian dissipation (e.g., due to cavity loss). We present an exact hybrid quantum-classical stochastic approach where the solution of a classical stochastic equation—mimicking the bosonic modes—is input into a quantum stochastic equation for the spins. Furthermore, the spins are effectively decoupled for each stochastic realization, which nevertheless comes at the expense of sampling over unphysical states. In contrast with existing stochastic approaches based on the influence functional formalism, we place no restriction (factorizability or Gaussianity) on the initial state, or the spin-boson coupling (except that it be linear in the bosonic operator). Markovian dissipation, being at the heart of our approach, renders the stochastic equations Markovian even in the strong coupling regime. Furthermore, it ensures hermiticity (though not positivity) of the density matrix for each realization, thus improving the convergence of stochastic sampling. Interestingly, we find a condition on the classical simulability of the system based solely on the single atom cooperativity even in a many-body setting. We benchmark and showcase the utility of our approach in several examples, specifically in cases where a direct numerical computation is unfeasible.
Digital quantum computers have the potential to simulate complex quantum systems. The spin-boson model is one of such systems, used in disparate physical domains. Importantly, in a number of setups, the spin-boson model is open, i.e., the system is in contact with an external environment which can, for instance, cause the decay of the spin state. Here, we study how to simulate such open quantum dynamics in a digital quantum computer, for which we use an IBM hardware. We consider in particular how accurate different implementations of the evolution result as a function of the level of noise in the hardware and of the parameters of the open dynamics. For the regimes studied, we show that the key aspect is to simulate the unitary portion of the dynamics, while the dissipative part can lead to a more noise-resistant simulation. We consider both a single spin coupled to a harmonic oscillator, and also two spins coupled to the oscillator. In the latter case, we show that it is possible to simulate the emergence of correlations between the spins via the oscillator.
We extend a standard stochastic theory to study open quantum systems coupled to a generic quantum environment. We exemplify the general framework by studying a two-level quantum system coupled bilinearly to the three fundamental classes of non-interacting particles: bosons, fermions, and spins. In this unified stochastic approach, the generalized stochastic Liouville equation (SLE) formally captures the exact quantum dissipations when noise variables with appropriate statistics for different bath models are applied. Anharmonic effects of a non-Gaussian bath are precisely encoded in the bath multi-time correlation functions that noise variables have to satisfy. Starting from the SLE, we devise a family of generalized hierarchical equations by averaging out the noise variables and expand bath multi-time correlation functions in a complete basis of orthonormal functions. The general hierarchical equations constitute systems of linear equations that provide numerically exact simulations of quantum dynamics. For bosonic bath models, our general hierarchical equation of motion reduces exactly to an extended version of hierarchical equation of motion which allows efficient simulation for arbitrary spectral densities and temperature regimes. Similar efficiency and flexibility can be achieved for the fermionic bath models within our formalism. The spin bath models can be simulated with two complementary approaches in the present formalism. (I) They can be viewed as an example of non-Gaussian bath models and be directly handled with the general hierarchical equation approach given their multi-time correlation functions. (II) Alternatively, each bath spin can be first mapped onto a pair of fermions and be treated as fermionic environments within the present formalism.