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.

The heavy rare earth elements crystallize into hexagonally close packed (h.c.p.) structures and share a common outer electronic configuration, differing only in the number of 4f electrons they have. These chemically inert 4f electrons set up localized magnetic moments, which are coupled via an indirect exchange interaction involving the conduction electrons. This leads to the formation of a wide variety of magnetic structures, the periodicities of which are often incommensurate with the underlying crystal lattice. Such incommensurate ordering is associated with a 'webbed' topology of the momentum space surface separating the occupied and unoccupied electron states (the Fermi surface). The shape of this surface-and hence the magnetic structure-for the heavy rare earth elements is known to depend on the ratio of the interplanar spacing c and the interatomic, intraplanar spacing a of the h.c.p. lattice. A theoretical understanding of this problem is, however, far from complete. Here, using gadolinium as a prototype for all the heavy rare earth elements, we generate a unified magnetic phase diagram, which unequivocally links the magnetic structures of the heavy rare earths to their lattice parameters. In addition to verifying the importance of the c/a ratio, we find that the atomic unit cell volume plays a separate, distinct role in determining the magnetic properties: we show that the trend from ferromagnetism to incommensurate ordering as atomic number increases is connected to the concomitant decrease in unit cell volume. This volume decrease occurs because of the so-called lanthanide contraction, where the addition of electrons to the poorly shielding 4f orbitals leads to an increase in effective nuclear charge and, correspondingly, a decrease in ionic radii.

Extreme pressure strongly affects the superconducting properties of "simple" elemental metals, such as Li, K, and Al. Pressure induces superconductivity in Li (as high as 17 K) while suppressing it in Al. We report first-principles investigations of the superconducting properties of dense Li, K, and Al based on a recently proposed, parameter-free, method. Our results show an unprecedented agreement with experiments, assess the predictive power of the method over a wide range of densities and electron-phonon couplings, and provide predictions for K, where no experiments exist so far. More importantly, our results help uncover the physics of the different behaviors of Li and Al in terms of phonon softening and Fermi surface nesting in Li.

We propose a simplified version of self-interaction corrected local spin-density (SIC-LSD) approximation, based on multiple scattering theory, which implements self-interaction correction locally, within the KKR method. The multiple scattering aspect of this new SIC-LSD method allows for the description of crystal potentials which vary from site to site in a random fashion and the calculation of physical quantities averaged over ensembles of such potentials using the coherent potential approximation (CPA). This facilitates applications of the SIC to alloys and pseudoalloys which could describe disordered local moment systems, as well as intermediate valences. As a demonstration of the method, we study the well-known α-γ phase transition in Ce, where we also explain how SIC operates in terms of multiple scattering theory.

An ab initio study of magnetic exchange interactions in antiferromagnetic and strongly correlated 3d transition metal monoxides is presented. Their electronic structure is calculated using the local self-interaction correction approach, implemented within the Korringa-Kohn-Rostoker band structure method, which is based on multiple scattering theory. The Heisenberg exchange constants are evaluated with the magnetic force theorem. Based on these the corresponding Néel temperatures TN and spin wave dispersions are calculated. The Néel temperatures are obtained using mean field approximation, random phase approximation and Monte Carlo simulations. The pressure dependence of TN is investigated using exchange constants calculated for different lattice constants. All the calculated results are compared to experimental data.

Solid MgB(2) has rather interesting and technologically important properties, such as a very high superconducting transition temperature. Focusing on this compound, we report the first nontrivial application of a novel density-functional-type theory for superconductors, recently proposed by the authors. Without invoking any adjustable parameters, we obtain the transition temperature, the gaps, and the specific heat of MgB(2) in very good agreement with experiment. Moreover, our calculations show how the Coulomb interaction acts differently on sigma and pi states, thereby stabilizing the observed superconducting phase.