We introduce a response theory for open quantum systems within nonequilibrium steady states subject to a Hamiltonian perturbation. Working in the weak system-bath coupling regime, our results are derived within the Lindblad-Gorini-Kossakowski-Sudarshan formalism. We find that the response of the system to a small perturbation is not simply related to a correlation function within the system, unlike traditional linear response theory in closed systems or expectations from the fluctuation-dissipation theorem. In limiting cases, when the perturbation is small relative to the coupling to the surroundings or when it does not lead to a change of the eigenstructure of the system, a perturbative expansion exists where the response function is related to a sum of a system correlation functions and additional forces induced by the surroundings. Away from these limiting regimes, however, the secular approximation results in a singular response that cannot be captured within the traditional approach but can be described by reverting to a microscopic Hamiltonian description. These findings are illustrated by explicit calculations in coupled qubits and anharmonic oscillators in contact with bosonic baths at different temperatures.

We present a strategy for mapping the dynamics of a fermionic quantum system to a set of classical dynamical variables. The approach is based on imposing the correspondence relation between the commutator and the Poisson bracket, preserving Heisenberg's equation of motion for one-body operators. In order to accommodate the effect of two-body terms, we further impose quantization on the spin-dependent occupation numbers in the classical equations of motion, with a parameter that is determined self-consistently. Expectation values for observables are taken with respect to an initial quasiclassical distribution that respects the original quantization of the occupation numbers. The proposed classical map becomes complete under the evolution of quadratic Hamiltonians and is extended for all even order observables. We show that the map provides an accurate description of the dynamics for an interacting quantum impurity model in the coulomb blockade regime, at both low and high temperatures. The numerical results are aided by a novel importance sampling scheme that employs a reference system to reduce significantly the sampling effort required to converge the classical calculations.

In this work, we investigate the characteristics of the electric current in the so-called symmetric Anderson impurity model. We study the nonequilibrium model using two complementary approximate methods, the perturbative quantum master equation approach to the reduced density matrix and a self-consistent equation of motion approach to the nonequilibrium Green's function. We find that, at a particular symmetry point, an interacting Anderson impurity model recovers the same steady-state current as an equivalent noninteracting model, akin a two-band resonant level model. We show this in the Coulomb blockade regime for both high and low temperatures, where either the approximate master equation approach or the Green's function method provides accurate results for the current. We conclude that the steady-state current in the symmetric Anderson model at this regime does not encode characteristics of a many-body interacting system.