The coherent quantum dynamics of a single bosonic spin variable, subject to a constraint derived from the quantum spherical model of a ferromagnet, and coupled to an external heat bath, is studied through the Lindblad equation for the reduced density matrix. Closed systems of equations of motion for several quantum observables are derived and solved exactly. The relationship to the single-mode Dicke model from quantum optics is discussed. The analysis of the interplay of the quantum fluctuation and the dissipation and their influence on the relaxation of the time-dependent magnetisation leads to the distinction of qualitatively different regimes of weak and strong quantum couplings. Considering the model's behaviour in an external field as a simple mean-field approximation of the dynamics of a quantum spherical ferromagnet, the magnetic phase diagram appears to be re-entrant and presents a quantum analogue of well-established classical examples of fluctuation-induced order.
The coherent quantum dynamics of a single bosonic spin variable, subject to a constraint derived from the quantum spherical model of a ferromagnet, and coupled to an external heat bath, is studied through the Lindblad equation for the reduced density matrix. Closed systems of equations of motion for several quantum observables are derived and solved exactly. The relationship to the single-mode Dicke model from quantum optics is discussed. The analysis of the interplay of the quantum fluctuation and the dissipation and their influence on the relaxation of the time-dependent magnetisation leads to the distinction of qualitatively different regimes of weak and strong quantum couplings. Considering the model's behaviour in an external field as a simple mean-field approximation of the dynamics of a quantum spherical ferromagnet, the magnetic phase diagram appears to be re-entrant and presents a quantum analogue of well-established classical examples of fluctuation-induced order.
Abstract The purely relaxational non-equilibrium dynamics of the quantum spherical model as described through a Lindblad equation is analysed. It is shown that the phenomenological requirements of reproducing the exact quantum equilibrium state as stationary solution and the associated classical Langevin equation in the classical limit g → 0 fix the form of the Lindblad dissipators, up to an overall time-scale. In the semi-classical limit, the models’ behaviour becomes effectively the one of the classical analogue, with a dynamical exponent z = 2 indicating diffusive transport, and an effective temperature T eff , renormalised by the quantum coupling g . A different behaviour is found for a quantum quench, at zero temperature, deep into the ordered phase g ≪ g c ( d ) , for d > 1 dimensions. Only for d = 2 dimensions, a simple scaling behaviour holds true, with a dynamical exponent z = 1 indicating ballistic transport, while for dimensions d ≠ 2 , logarithmic corrections to scaling arise. The spin–spin correlator, the growing length scale and the time-dependent susceptibility show the existence of several logarithmically different length scales.
The spin-boson model, describing a two-level system coupled to a bosonic bath, is extensively studied as a paradigmatic dissipative quantum system, exhibiting rich dynamical behavior and even a localization transition in the strong coupling regime. Here, we additionally consider dephasing as a source of Markovian dissipation on top of the non-Markovian dynamics due to an Ohmic bath, and investigate the dynamics of the spin. We show that the characteristic frequency of the spin dynamics, while strongly renormalized by the bosonic bath, changes in a simple fashion (or does not change at all) with dephasing. To obtain these results, we develop an exact nonperturbative method known as the stochastic Schr\"odinger equation, mimicking the Ohmic bath via a stochastic magnetic field combined with the Lindblad quantum master equation due to dephasing, which allows us to numerically compute the dynamics. Furthermore, we derive weak-coupling analytic results utilizing the well-known noninteracting blip approximation. Our findings are relevant to quantum simulation of the spin-boson model in the regime of strong coupling in trapped ions and circuit QED architectures among others.
Abstract Using a dissipative quantum Rabi model, we study the dynamics of a slow qubit coupled to a fast quantum harmonic oscillator interacting with a bosonic bath from weak to strong and ultra-strong coupling regimes. Solving the quantum Heisenberg equations of motion, perturbative in the internal coupling between qubit and oscillator, we derive functional relationships directly linking the qubit coordinates in the Bloch sphere to oscillator observables. We then perform accurate time-dependent Matrix Product State simulations and compare our results both with the analytical solutions of the Heisenberg equations of motion, and with numerical solutions of a Lindblad master equation, perturbative in the external coupling between oscillator and environment. Indeed, we show that, up to the strong coupling regime, the qubit state accurately fulfils the derived functional relationships. We analyse in detail the case of a qubit starting with generic coordinates on the Bloch sphere of which we evaluate the three components of the Bloch vector through the averages of oscillator observables. Interestingly, a weak to intermediate oscillator coupling to the bath is able to simplify the Bloch vector evaluation since qubit-oscillator relationships are more immediate. Moreover, by monitoring the qubit fidelity with respect to free limit, we find the parameter regime where the combined effect of internal and external couplings is able to hinder the reliable evaluation of the qubit Bloch vector. Finally, in the ultra-strong coupling regime, non-Markovian effects become robust and the dynamics of qubit and oscillator are inextricably entangled making the qubit Bloch vector evaluation difficult. Graphical Abstract
We study the Hamiltonian dynamics of the spherical spin model with fully-connected two-body random interactions. In the statistical physics framework, the potential energy is of the so-called p = 2 kind, closely linked to the scalar field theory. Most importantly for our setting, the energy conserving dynamics are equivalent to the ones of the Neumann integrable model. We take initial conditions from the Boltzmann equilibrium measure at a temperature that can be above or below the static phase transition, typical of a disordered (paramagnetic) or of an ordered (disguised ferromagnetic) equilibrium phase. We subsequently evolve the configurations with Newton dynamics dictated by a different Hamiltonian, obtained from an instantaneous global rescaling of the elements in the interaction random matrix. In the limit of infinitely many degrees of freedom, , we identify three dynamical phases depending on the parameters that characterise the initial state and the final Hamiltonian. We next set the analysis of the system with finite number of degrees of freedom in terms of N non-linearly coupled modes. We argue that in the limit the modes decouple at long times. We evaluate the mode temperatures and we relate them to the frequency-dependent effective temperature measured with the fluctuation-dissipation relation in the frequency domain, similarly to what was recently proposed for quantum integrable cases. Finally, we analyse the N - 1 integrals of motion, notably, their scaling with N, and we use them to show that the system is out of equilibrium in all phases, even for parameters that show an apparent Gibbs-Boltzmann behaviour of the global observables. We elaborate on the role played by these constants of motion after the quench and we briefly discuss the possible description of the asymptotic dynamics in terms of a generalised Gibbs ensemble.
We present a general approach to derive Lindblad master equations for a subsystem whose dynamics is coupled to dissipative bosonic modes. The derivation relies on a Schrieffer-Wolff transformation which allows us to eliminate the bosonic degrees of freedom after self-consistently determining their state as a function of the coupled quantum system. We apply this formalism to the dissipative Dicke model and derive a Lindblad master equation for the atomic spins, which includes the coherent and dissipative interactions mediated by the bosonic mode. This master equation accurately predicts the Dicke phase transition and gives the correct steady state. In addition, we compare the dynamics using exact diagonalization and numerical integration of the master equation with the predictions of semiclassical trajectories. We finally test the performance of our formalism by studying the relaxation of a NOON state and show that the dynamics captures quantum metastability.
Abstract We study the Hamiltonian dynamics of the spherical spin model with fully-connected two-body random interactions. In the statistical physics framework, the potential energy is of the so-called p = 2 kind, closely linked to the scalar field theory. Most importantly for our setting, the energy conserving dynamics are equivalent to the ones of the Neumann integrable model. We take initial conditions from the Boltzmann equilibrium measure at a temperature that can be above or below the static phase transition, typical of a disordered (paramagnetic) or of an ordered (disguised ferromagnetic) equilibrium phase. We subsequently evolve the configurations with Newton dynamics dictated by a different Hamiltonian, obtained from an instantaneous global rescaling of the elements in the interaction random matrix. In the limit of infinitely many degrees of freedom, , we identify three dynamical phases depending on the parameters that characterise the initial state and the final Hamiltonian. We next set the analysis of the system with finite number of degrees of freedom in terms of N non-linearly coupled modes. We argue that in the limit the modes decouple at long times. We evaluate the mode temperatures and we relate them to the frequency-dependent effective temperature measured with the fluctuation-dissipation relation in the frequency domain, similarly to what was recently proposed for quantum integrable cases. Finally, we analyse the N − 1 integrals of motion, notably, their scaling with N , and we use them to show that the system is out of equilibrium in all phases, even for parameters that show an apparent Gibbs–Boltzmann behaviour of the global observables. We elaborate on the role played by these constants of motion after the quench and we briefly discuss the possible description of the asymptotic dynamics in terms of a generalised Gibbs ensemble.