We report on the background, current status, and current lines of development of the octopus project. This program materializes the main equations of density-functional theory in the ground state, and of timedependent density-functional theory for dynamical effects. The focus is nowadays placed on the optical (i.e. electronic) linear response properties of nanostructures and biomolecules, and on the non-linear response to high-intensity fields of finite systems, with particular attention to the coupled ionic-electronic motion (i.e. photo-chemical processes). In addition, we are currently extending the code to the treatment of periodic systems (both to one-dimensional chains, two-dimensional slabs, or fully periodic solids), magnetic properties (ground state properties and excitations), and to the field of quantum-mechanical transport or "molecular electronics." In this communication, we concentrate on the development of the methodology: we review the essential numerical schemes used in the code, and report on the most recent implementations, with special attention to the introduction of adaptive coordinates, to the extension of our real-space technique to tackle periodic systems, and on large-scale parallelization. More information on the code, as well as the code itself, can be found at http://www.tddft.org/programs/octopus/.

Coherent ultrafast charge transfer in an organic photovoltaic blend

We present a reciprocal space analytical method to cut off the long range interactions in supercell calculations for systems that are infinite and periodic in one or two dimensions, generalizing previous work to treat finite systems. The proposed cutoffs are functions in Fourier space, that are used as a multiplicative factor to screen the bare Coulomb interaction. The functions are analytic everywhere except in a subdomain of the Fourier space that depends on the periodic dimensionality. We show that the divergences that lead to the nonanalytical behavior can be exactly canceled when both the ionic and the Hartree potential are properly screened. This technique is exact, fast, and very easy to implement in already existing supercell codes. To illustrate the performance of the scheme, we apply it to the case of the Coulomb interaction in systems with reduced periodicity ͑as one-dimensional chains and layers͒. For these test cases, we address the impact of the cutoff on different relevant quantities for ground and excited state properties, namely: the convergence of the ground state properties, the static polarizability of the system, the quasiparticle corrections in the GW scheme, and the binding energy of the excitonic states in the Bethe-Salpeter equation. The results are very promising and easy to implement in all available first-principles codes.

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).

The optical excitation of organic semiconductors not only generates charge-neutral electron-hole pairs (excitons), but also charge-separated polaron pairs with high yield. The microscopic mechanisms underlying this charge separation have been debated for many years. Here we use ultrafast two-dimensional electronic spectroscopy to study the dynamics of polaron pair formation in a prototypical polymer thin film on a sub-20-fs time scale. We observe multi-period peak oscillations persisting for up to about 1 ps as distinct signatures of vibronic quantum coherence at room temperature. The measured two-dimensional spectra show pronounced peak splittings revealing that the elementary optical excitations of this polymer are hybridized exciton-polaron-pairs, strongly coupled to a dominant underdamped vibrational mode. Coherent vibronic coupling induces ultrafast polaron pair formation, accelerates the charge separation dynamics and makes it insensitive to disorder. These findings open up new perspectives for tailoring light-to-current conversion in organic materials.