We develop a numerical procedure to efficiently model the nonequilibrium steady state of one-dimensional arrays of open quantum systems based on a matrix-product operator ansatz for the density matrix. The procedure searches for the null eigenvalue of the Liouvillian superoperator by sweeping along the system while carrying out a partial diagonalization of the single-site stationary problem. It bears full analogy to the density-matrix renormalization-group approach to the ground state of isolated systems, and its numerical complexity scales as a power law with the bond dimension. The method brings considerable advantage when compared to the integration of the time-dependent problem via Trotter decomposition, as it can address arbitrarily long-ranged couplings. Additionally, it ensures numerical stability in the case of weakly dissipative systems thanks to a slow tuning of the dissipation rates along the sweeps. We have tested the method on a driven-dissipative spin chain, under various assumptions for the Hamiltonian, drive, and dissipation parameters, and compared the results to those obtained both by Trotter dynamics and Monte Carlo wave function methods. Accurate and numerically stable convergence was always achieved when applying the method to systems with a gapped Liouvillian and a nondegenerate steady state.

We study the physics of quantum phase transitions from the perspective of non-equilibrium thermodynamics. For first order quantum phase transitions, we find that the average work done per quench in crossing the critical point is discontinuous. This leads us to introduce the quantum latent work in analogy with the classical latent heat of first order classical phase transitions. For second order quantum phase transitions the irreversible work is closely related to the fidelity susceptibility for weak sudden quenches of the system Hamiltonian. We demonstrate our ideas with numerical simulations of first, second, and infinite order phase transitions in various spin chain models.

We study high-temperature spin transport through an anisotropic spin-1 2 Heisenberg chain in which integrability is broken by a single impurity close to the center of the chain. For a finite impurity strength, the level spacing statistics of this model is known to be Wigner-Dyson. Our aim is to understand if this integrability breaking is manifested in the high-temperature spin transport. We focus first on the nonequilibrium steady state (NESS), where the chain is connected to spin baths that act as sources and sinks for spin excitations at the boundaries. Using a combination of open quantum system theory and matrix product operators techniques, we extract the transport properties by means of a finite-size scaling of the spin current in the NESS. Our results indicate that, despite of the formation of a partial domain wall in the steady state magnetization (and despite the Wigner-Dyson level spacing distribution of the model), transport remains ballistic. We contrast this behavior with the one produced by a staggered magnetic field in the XXZ chain, for which it is known that transport is diffusive. By performing a numerical computation of the real part of the spin conductivity, we show that our findings are consistent with linear response theory. We discuss subtleties associated with the apparent vanishing of the Drude in the presence of an impurity.

Optical transport represents a natural route towards fast communications, and it is currently used in large scale data transfer. The progressive miniaturization of devices for information processing calls for the microscopic tailoring of light transport and confinement at length scales appropriate for upcoming technologies. With this goal in mind, we present a theoretical analysis of a one-dimensional Fabry-Perot interferometer built with two highly saturable nonlinear mirrors: a pair of two-level systems. Our approach captures nonlinear and nonreciprocal effects of light transport that were not reported previously. Remarkably, we show that such an elementary device can operate as a microscopic integrated optical rectifier.

By using a fully quantum approach based on an input-output formulation of the stochastic Schrödinger equation, we show rectification of radiation fields in a one-dimensional waveguide doped with a pair of ideal two-level systems for three topical cases: classical driving, under the action of noise, and single photon pulsed excitation. We show that even under the constant action of unwanted noise the device still operates effectively as an optical isolator, which is of critical importance for noise resistance. Finally, harnessing stimulated emission allows for non-reciprocal behavior for single photon inputs, thus showing purely quantum rectification at the single-photon level. The latter is a considerable step towards the ultimate goal of devising an unconditional quantum rectifier for arbitrary quantum states.