We demonstrate that with a stepwise introduction of complexity to a model of an electron system embedded in a photonic cavity and a carefully controlled stepwise truncation of the ensuing many-body space it is possible to describe the time-dependent transport of electrons through the system with a non-Markovian generalized quantum master equation. We show how this approach retains effects of an external magnetic field and the geometry of an anisotropic electronic system. The Coulomb interaction between the electrons and the full electromagnetic coupling between the electrons and the photons are treated in a non-perturbative way using "exact numerical diagonalization". * vidar@hi.is †

We use a non-Markovian master equation to describe the transport of Coulomb interacting electrons through an electromagnetic cavity with one quantized photon mode. The central system is a finite parabolic quantum wire that is coupled weakly to external parabolic quasi-one-dimensional leads at $t=0$. With a stepwise introduction of complexity to the description of the system and a corresponding stepwise truncation of the ensuing many-body spaces we are able to describe the time-dependent transport of Coulomb-interacting electrons through a geometrically complex central system. We take into account the full electromagnetic interaction of electrons and cavity photons without resorting to the rotating wave approximation or reduction of the electron states to two levels. We observe that the number of initial cavity photons and their polarization can have important effects on the transport properties of the system. The quasiparticles formed in the central system have a lifetime limited by the coupling to the leads and radiation processes active on a much longer timescale.Comment: RevTeX (pdf-LaTeX) 11 pages with 12 jpg-figures include

We outline a rigorous method which can be used to solve the many-body Schrödinger equation for a Coulomb interacting electronic system in an external classical magnetic field as well as a quantized electromagnetic field. Effects of the geometry of the electronic system as well as the polarization of the quantized electromagnetic field are explicitly taken into account. We accomplish this by performing repeated truncations of many-body spaces in order to keep the size of the many particle basis on a manageable level. The electron-electron and electron-photon interactions are treated in a nonperturbative manner using "exact numerical diagonalization". Our results demonstrate that including the diamagnetic term in the photon-electron interaction Hamiltonian drastically improves numerical convergence. Additionally, convergence with respect to the number of photon states in the joint photon-electron Fock space basis is fast. However, the convergence with respect to the number of electronic states is slow and is the main bottleneck in calculations.

We investigate the coupling between a quantized electromagnetic field in a cavity resonator and a Coulomb interacting electronic system in a nanostructure in an external magnetic field. The effects caused by the geometry of the electronic system and the polarization of the electromagnetic field are explicitly taken into account. Our numerical results demonstrate that the twolevel system approximation and the Jaynes-Cummings model remain valid in the weak electron-photon coupling regime, while the quadratic vector potential in the diamagnetic part of the charge current leads to significant correction to the energy spectrum in the strong coupling regime. Furthermore, we find that coupling to a strong cavity photon mode polarizes the charge distribution of the system, requiring a large basis of single-electron eigenstates to be included in the model.

We derive a Markovian master equation for the single-electron density matrix, applicable to quantum cascade lasers (QCLs). The equation conserves the positivity of the density matrix, includes off-diagonal elements (coherences) as well as in-plane dynamics, and accounts for electron scattering with phonons and impurities. We use the model to simulate a terahertzfrequency QCL, and compare the results with both experiment and simulation via nonequilibrium Green's functions (NEGF). We obtain very good agreement with both experiment and NEGF when the QCL is biased for optimal lasing. For the considered device, we show that the magnitude of coherences can be a significant fraction of the diagonal matrix elements, which demonstrates their importance when describing THz QCLs. We show that the in-plane energy distribution can deviate far from a heated Maxwellian distribution, which suggests that the assumption of thermalized subbands in simplified density-matrix models is inadequate. We also show that the current density and subband occupations relax towards their steady-state values on very different time scales.