We present a local control scheme to construct the external potential v that, for a given initial state, produces a prescribed time-dependent density in an interacting quantum many-body system. This numerical method is efficient and stable even for large and rapid density variations irrespective of the initial state and the interactions. It can at the same time be used to answer fundamental v-representability questions in density functional theory. In particular, in the absence of interactions, it allows us to construct the exact time-dependent Kohn-Sham potential for arbitrary initial states. We illustrate the method in a correlated one-dimensional two-electron system with different interactions, initial states and densities. For a Kohn-Sham system with a correlated initial state we demonstrate the interplay between memory and initial-state dependence as well as the failure of any adiabatic approximation.

We present a first-principles approach to electronic many-body systems strongly coupled to cavity modes in terms of matter–photon one-body reduced density matrices. The theory is fundamentally nonperturbative and thus captures not only the effects of correlated electronic systems but accounts also for strong interactions between matter and photon degrees of freedom. We do so by introducing a higher-dimensional auxiliary system that maps the coupled fermion-boson system to a dressed fermionic problem. This reformulation allows us to overcome many fundamental challenges of density-matrix theory in the context of coupled fermion-boson systems and we can employ conventional reduced density-matrix functional theory developed for purely fermionic systems. We provide results for one-dimensional model systems in real space and show that simple density-matrix approximations are accurate from the weak to the deep-strong coupling regime. This justifies the application of our method to systems that are too complex for exact calculations and we present first results, which show that the influence of the photon field depends sensitively on the details of the electronic structure.

When running time-dependent density functional theory (TDDFT) calculations for real-time simulations of non-equilibrium dynamics, the user has a choice of initial Kohn-Sham state, and typically a Slater determinant is used. We explore the impact of this choice on the exchange-correlation potential when the physical system begins in a 50 : 50 superposition of the ground and first-excited state of the system. We investigate the possibility of judiciously choosing a Kohn-Sham initial state that minimizes errors when adiabatic functionals are used. We find that if the Kohn-Sham state is chosen to have a configuration matching the one that dominates the interacting state, this can be achieved for a finite time duration for some but not all such choices. When the Kohn-Sham system does not begin in a Slater determinant, we further argue that the conventional splitting of the exchange-correlation potential into exchange and correlation parts has limited value, and instead propose a decomposition into a "single-particle" contribution that we denote v, and a remainder. The single-particle contribution can be readily computed as an explicit orbital-functional, reduces to exchange in the Slater determinant case, and offers an alternative to the adiabatic approximation as a starting point for TDDFT approximations.