Negative differential conductance induced by spin-charge separation

Spin-charge states of correlated electrons in a one-dimensional quantum dot attached to interacting leads are studied in the non-linear transport regime. With non-symmetric tunnel barriers, regions of negative differential conductance induced by spin-charge separation are found. They are due to a correlation-induced trapping of higher-spin states without magnetic field, and associated with a strong increase in the fluctuations of the electron spin.Comment: REVTEX, 4 pages including 3 figures; Accepted for publication on Physical Review Letter

Resonators coupled to voltage-biased Josephson junctions: From linear response to strongly driven nonlinear oscillations

Motivated by recent experiments, where a voltage biased Josephson junction is placed in series with a resonator, the classical dynamics of the circuit is studied in various domains of parameter space. This problem can be mapped onto the dissipative motion of a single degree of freedom in a nonlinear time-dependent potential, where in contrast to conventional settings the nonlinearity appears in the driving while the static potential is purely harmonic. For long times the system approaches steady states which are analyzed in the underdamped regime over the full range of driving parameters including the fundamental resonance as well as higher and sub-harmonics. Observables such as the dc-Josephson current and the radiated microwave power give direct information about the underlying dynamics covering phenomena as bifurcations, irregular motion, up- and down conversion. Due to their tunability, present and future set-ups provide versatile platforms to explore the changeover from linear response to strongly nonlinear behavior in driven dissipative systems under well defined conditions.Comment: 12 pages, 11 figure

Exact c-number Representation of Non-Markovian Quantum Dissipation

The reduced dynamics of a quantum system interacting with a linear heat bath finds an exact representation in terms of a stochastic Schr{\"o}dinger equation. All memory effects of the reservoir are transformed into noise correlations and mean-field friction. The classical limit of the resulting stochastic dynamics is shown to be a generalized Langevin equation, and conventional quantum state diffusion is recovered in the Born--Markov approximation. The non-Markovian exact dynamics, valid at arbitrary temperature and damping strength, is exemplified by an application to the dissipative two-state system.Comment: 4 pages, 2 figures. To be published in Phys. Rev. Let

Efficient low temperature simulations for fermionic reservoirs with the hierarchical equations of motion method: Application to the Anderson impurity model

The hierarchical equations of motion (HEOM) approach is an accurate method to simulate open system quantum dynamics, which allows for systematic convergence to numerically exact results. To represent the effects of the bath, the reservoir correlation functions are usually decomposed into the summation of multiple exponential terms in the HEOM method. Since the reservoir correlation functions become highly non-Markovian at low temperatures or when the bath has complex band structures, a present challenge is to obtain accurate exponential decompositions that allow efficient simulation with the HEOM. In this work, we employ the barycentric representation to approximate the Fermi function and hybridization functions in the frequency domain. The new method, by approximating these functions with optimized rational decomposition, greatly reduces the number of basis functions in decomposing the reservoir correlation functions, which further allows the HEOM method to be applied to ultra-low temperature and general bath structures. We demonstrate the efficiency, accuracy, and long-time stability of the new decomposition scheme by applying it to the Anderson impurity model (AIM) in the low-temperature regime with the Lorentzian and tight-binding hybridization functions.Comment: 39 pages, 12 figure

Taming quantum noise for efficient low temperature simulations of open quantum systems

The hierarchical equations of motion (HEOM), derived from the exact Feynman-Vernon path integral, is one of the most powerful numerical methods to simulate the dynamics of open quantum systems that are embedded in thermal environments. However, its applicability is restricted to specific forms of spectral reservoir distributions and relatively elevated temperatures. Here we solve this problem and introduce an effective treatment of quantum noise in frequency space by systematically clustering higher order Matsubara poles equivalent to an optimized rational decomposition. This leads to an elegant extension of the HEOM to arbitrary temperatures and very general reservoirs in combination with efficiency, high accuracy and long-time stability. Moreover, the technique can directly be implemented in alternative approaches such as Greens-function, stochastic, and pseudo-mode formulations. As one highly non-trivial application, for the sub-ohmic spin-boson model at vanishing temperature the Shiba relation is quantitatively verified which predicts the long-time decay of correlation functions.Comment: 11pages, 8 figures. v2: grant number correcte

Nonlinear response theory for lossy superconducting quantum circuits

We introduce a numerically exact and yet computationally feasible nonlinear response theory developed for lossy superconducting quantum circuits based on a framework of quantum dissipation in a minimally extended state space. Starting from the Feynman--Vernon path integral formalism for open quantum systems with the system degrees of freedom being the nonlinear elements of the circuit, we eliminate the temporally non-local influence functional of all linear elements by introducing auxiliary harmonic modes with complex-valued frequencies coupled to the non-linear degrees of freedom of the circuit. In our work, we propose a concept of time-averaged observables, inspired by experiment, and provide an explicit formula for producing their quasiprobability distribution. Furthermore, we systematically derive a weak-coupling approximation in the presence of a drive, and demonstrate the applicability of our formalism through a study on the dispersive readout of a superconducting qubit. The developed framework enables a comprehensive fully quantum-mechanical treatment of nonlinear quantum circuits coupled to their environment, without the limitations of typical approaches to weak dissipation, high temperature, and weak drive. Furthermore, we discuss the implications of our findings to the quantum measurement theory.Comment: 18 pages, 4 figures, 2 table

Brachistochrone of Entanglement for Spin Chains

We analytically investigate the role of entanglement in time-optimal state evolution as an appli- cation of the quantum brachistochrone, a general method for obtaining the optimal time-dependent Hamiltonian for reaching a target quantum state. As a model, we treat two qubits indirectly cou- pled through an intermediate qubit that is directly controllable, which represents a typical situation in quantum information processing. We find the time-optimal unitary evolution law and quantify residual entanglement by the two-tangle between the indirectly coupled qubits, for all possible sets of initial pure quantum states of a tripartite system. The integrals of the motion of the brachistochrone are determined by fixing the minimal time at which the residual entanglement is maximized. Entan- glement plays a role for W and GHZ initial quantum states, and for the bi-separable initial state in which the indirectly coupled qubits have a nonzero value of the 2-tangle.Comment: 9 pages, 4 figure

Scaling and universality in the anisotropic Kondo model and the dissipative two-state system

Scaling and universality in the Ohmic two-state system is investigated by exploiting the equivalence of this model to the anisotropic Kondo model. For the Ohmic two-state system, we find universal scaling functions for the specific heat, $C_{\alpha}(T)$, static susceptibility, $\chi_{\alpha}(T)$, and spin relaxation function $S_{\alpha}(\omega)$ depending on the reduced temperature $T/\Delta_{r}$ (frequency $\omega/\Delta_{r}$), with $\Delta_{r}$ the renormalized tunneling frequency, and uniquely specified by the dissipation strength $\alpha$ ($0<\alpha<1$). The scaling functions can be used to extract $\alpha$ and $\Delta_{r}$ in experimental realizations.Comment: 5 pages (LaTeX), 4 EPS figures. Minor changes, typos corrected, journal reference adde

Coherence correlations in the dissipative two-state system

We study the dynamical equilibrium correlation function of the polaron-dressed tunneling operator in the dissipative two-state system. Unlike the position operator, this coherence operator acts in the full system-plus-reservoir space. We calculate the relevant modified influence functional and present the exact formal expression for the coherence correlations in the form of a series in the number of tunneling events. For an Ohmic spectral density with the particular damping strength $K=1/2$, the series is summed in analytic form for all times and for arbitrary values of temperature and bias. Using a diagrammatic approach, we find the long-time dynamics in the regime $K<1$. In general, the coherence correlations decay algebraically as $t^{-2K}$ at T=0. This implies that the linear static susceptibility diverges for $K\le 1/2$ as $T\to 0$, whereas it stays finite for $K>1/2$ in this limit. The qualitative differences with respect to the asymptotic behavior of the position correlations are explained.Comment: 19 pages, 4 figures, to be published in Phys. Rev.