Critical analysis of the response function in low-dimensional materials.

The presence of sharp peaks in the real part of the static dielectric response function are usually accepted as indication of charge or spin instabilities in a material. However, there are misconceptions that Fermi surface (FS) nesting guarantees a peak in the response function like in one-dimensional systems, and, in addition, response function matrix elements between empty and occupied states are usually considered of secondary importance and typically set to unity like in the free electron gas case. In this work, we explicitly show, through model systems and real materials, within the framework of density functional theory, that predictions about the peaks in the response function, using FS nesting and constant matrix elements yields erroneous conclusions. We find that the inclusion of the matrix elements completely alters the structure of the response function. In all the cases studied other than the one-dimensional case we find that the inclusion of matrix elements washes out the structure found with constant matrix elements. Our conclusion is that it is imperative to calculate the full response function, with matrix elements, when making predictions about instabilities in novel materials.

Siesta: Recent developments and applications.

A review of the present status, recent enhancements, and applicability of the Siesta program is presented. Since its debut in the mid-1990s, Siesta's flexibility, efficiency, and free distribution have given advanced materials simulation capabilities to many groups worldwide. The core methodological scheme of Siesta combines finite-support pseudo-atomic orbitals as basis sets, norm-conserving pseudopotentials, and a real-space grid for the representation of charge density and potentials and the computation of their associated matrix elements. Here, we describe the more recent implementations on top of that core scheme, which include full spin-orbit interaction, non-repeated and multiple-contact ballistic electron transport, density functional theory (DFT)+U and hybrid functionals, time-dependent DFT, novel reduced-scaling solvers, density-functional perturbation theory, efficient van der Waals non-local density functionals, and enhanced molecular-dynamics options. In addition, a substantial effort has been made in enhancing interoperability and interfacing with other codes and utilities, such as wannier90 and the second-principles modeling it can be used for, an AiiDA plugin for workflow automatization, interface to Lua for steering Siesta runs, and various post-processing utilities. Siesta has also been engaged in the Electronic Structure Library effort from its inception, which has allowed the sharing of various low-level libraries, as well as data standards and support for them, particularly the PSeudopotential Markup Language definition and library for transferable pseudopotentials, and the interface to the ELectronic Structure Infrastructure library of solvers. Code sharing is made easier by the new open-source licensing model of the program. This review also presents examples of application of the capabilities of the code, as well as a view of on-going and future developments.

Band unfolding made simple.

We present a simple view on band unfolding of the energy bands obtained from supercell calculations. It relies on the relationship between the local density of states in reciprocal space and the fully unfolded band structure. This provides an intuitive and valid approach not only for periodic, but also for systems with no translational symmetry. By refolding into the primitive Brillouin zone of the pristine crystal we recover the conventional unfolded bands. We implement our algorithm in the Siesta package. As an application, we study a set of benchmark examples, ranging from simple defects on crystals to systems with increasing complexity and of current interest, as the effect of external pressure on rotated graphene bilayers.