The density functional theory for superconductors developed in the preceding article is applied to the calculation of superconducting properties of several elemental metals. In particular, we present results for the transition temperature, for the gap at zero temperature, and for thermodynamic properties like the specific heat. We obtain an unprecedented agreement with experimental results. Superconductors both with strong and weak electron-phonon coupling are equally well described. This demonstrates that, as far as conventional superconductivity is concerned, the first-principles prediction of superconducting properties is feasible.

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

A novel approach to the description of superconductors in thermal equilibrium is developed within a formally exact density-functional framework. The theory is formulated in terms of three "densities": the ordinary electron density, the superconducting order parameter, and the diagonal of the nuclear N -body density matrix. The electron density and the order parameter are determined by Kohn-Sham equations that resemble the Bogoliubov-de Gennes equations. The nuclear density matrix follows from a Schrödinger equation with an effective N -body interaction. These equations are coupled to each other via exchange-correlation potentials which are universal functionals of the three densities. Approximations of these exchange-correlation functionals are derived using the diagrammatic techniques of many-body perturbation theory. The bare Coulomb repulsion between the electrons and the electron-phonon interaction enter this perturbative treatment on the same footing. In this way, a truly ab-initio description is achieved which does not contain any empirical parameters.