Over the last few years, extraordinary advances in experimental and theoretical tools have allowed us to monitor and control matter at short time and atomic scales with a high degree of precision. An appealing and challenging route toward engineering materials with tailored properties is to find ways to design or selectively manipulate materials, especially at the quantum level. To this end, having a state-of-the-art ab initio computer simulation tool that enables a reliable and accurate simulation of light-induced changes in the physical and chemical properties of complex systems is of utmost importance. The first principles real-space-based Octopus project was born with that idea in mind, i.e., to provide a unique framework that allows us to describe non-equilibrium phenomena in molecular complexes, low dimensional materials, and extended systems by accounting for electronic, ionic, and photon quantum mechanical effects within a generalized time-dependent density functional theory. This article aims to present the new features that have been implemented over the last few years, including technical developments related to performance and massive parallelism. We also describe the major theoretical developments to address ultrafast light-driven processes, such as the new theoretical framework of quantum electrodynamics density-functional formalism for the description of novel light–matter hybrid states. Those advances, and others being released soon as part of the Octopus package, will allow the scientific community to simulate and characterize spatial and time-resolved spectroscopies, ultrafast phenomena in molecules and materials, and new emergent states of matter (quantum electrodynamical-materials).

We review the GPAW open-source Python package for electronic structure calculations. GPAW is based on the projector-augmented wave method and can solve the self-consistent density functional theory (DFT) equations using three different wave-function representations, namely real-space grids, plane waves, and numerical atomic orbitals. The three representations are complementary and mutually independent and can be connected by transformations via the real-space grid. This multi-basis feature renders GPAW highly versatile and unique among similar codes. By virtue of its modular structure, the GPAW code constitutes an ideal platform for the implementation of new features and methodologies. Moreover, it is well integrated with the Atomic Simulation Environment (ASE), providing a flexible and dynamic user interface. In addition to ground-state DFT calculations, GPAW supports many-body GW band structures, optical excitations from the Bethe–Salpeter Equation, variational calculations of excited states in molecules and solids via direct optimization, and real-time propagation of the Kohn–Sham equations within time-dependent DFT. A range of more advanced methods to describe magnetic excitations and non-collinear magnetism in solids are also now available. In addition, GPAW can calculate non-linear optical tensors of solids, charged crystal point defects, and much more. Recently, support for graphics processing unit (GPU) acceleration has been achieved with minor modifications to the GPAW code thanks to the CuPy library. We end the review with an outlook, describing some future plans for GPAW.

Experiments at the interface of quantum-optics and chemistry have revealed that strong coupling between light and matter can substantially modify chemical and physical properties of molecules and solids. While the theoretical description of such situations is usually based on non-relativistic quantum electrodynamics, which contains quadratic light-matter coupling terms, it is commonplace to disregard these terms and restrict to purely bilinear couplings. In this work we clarify the physical origin and the substantial impact of the most common quadratic terms, the diamagnetic and self-polarization terms, and highlight why neglecting them can lead to rather unphysical results. Specifically we demonstrate its relevance by showing that neglecting it leads to the loss of gauge invariance, basis-set dependence, disintegration (loss of bound states) of any system in the basis set-limit, unphysical radiation of the ground state and an artificial dependence on the static dipole. Besides providing important guidance for modeling strongly coupled light-matter systems, the presented results do also indicate under which conditions those effects might become accessible.

Strong light–matter interaction in cavity environments is emerging as a promising approach to control chemical reactions in a non-intrusive and efficient manner. The underlying mechanism that distinguishes between steering, accelerating, or decelerating a chemical reaction has, however, remained unclear, hampering progress in this frontier area of research. We leverage quantum-electrodynamical density-functional theory to unveil the microscopic mechanism behind the experimentally observed reduced reaction rate under cavity induced resonant vibrational strong light-matter coupling. We observe multiple resonances and obtain the thus far theoretically elusive but experimentally critical resonant feature for a single strongly coupled molecule undergoing the reaction. While we describe only a single mode and do not explicitly account for collective coupling or intermolecular interactions, the qualitative agreement with experimental measurements suggests that our conclusions can be largely abstracted towards the experimental realization. Specifically, we find that the cavity mode acts as mediator between different vibrational modes. In effect, vibrational energy localized in single bonds that are critical for the reaction is redistributed differently which ultimately inhibits the reaction.

Significance: Excitation and charge transfer are fundamental processes in nature, and controlling these processes is a major goal of quantum chemistry. While these processes are well understood for the usual free-space case, when the electromagnetic vacuum is changed due to, e.g., a cavity, these processes can be dramatically different. We consider these changes in transfer processes with real-space donor–acceptor models, where we put an emphasis on the impact of electron–electron correlations. We find results in line with recent experiments, where strong light–matter interaction leads to enhanced transfer reactions, even when in the corresponding free-space situation no transfer should be possible. We highlight that the processes depend crucially on the Coulomb and self-polarization interactions.