The CECAM electronic structure library and the modular software development paradigm.

First-principles electronic structure calculations are now accessible to a very large community of users across many disciplines, thanks to many successful software packages, some of which are described in this special issue. The traditional coding paradigm for such packages is monolithic, i.e.

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.

Entropic bonding of the type 1 pilus from experiment and simulation.

The type 1 pilus is a bacterial filament consisting of a long coiled proteic chain of subunits joined together by non-covalent bonding between complementing -strands. Its strength and structural stability are critical for its anchoring function in uropathogenic bacteria. The pulling and unravelling of the FimG subunit of the pilus was recently studied by atomic force microscopy experiments and steered molecular dynamics simulations (Alonso-Caballero 2018 .

The ONETEP linear-scaling density functional theory program.

We present an overview of the onetep program for linear-scaling density functional theory (DFT) calculations with large basis set (plane-wave) accuracy on parallel computers. The DFT energy is computed from the density matrix, which is constructed from spatially localized orbitals we call Non-orthogonal Generalized Wannier Functions (NGWFs), expressed in terms of periodic sinc (psinc) functions. During the calculation, both the density matrix and the NGWFs are optimized with localization constraints.

QuantumATK: an integrated platform of electronic and atomic-scale modelling tools.

QuantumATK is an integrated set of atomic-scale modelling tools developed since 2003 by professional software engineers in collaboration with academic researchers. While different aspects and individual modules of the platform have been previously presented, the purpose of this paper is to give a general overview of the platform. The QuantumATK simulation engines enable electronic-structure calculations using density functional theory or tight-binding model Hamiltonians, and also offers bonded or reactive empirical force fields in many different parametrizations.

SIESTA-SIPs: Massively parallel spectrum-slicing eigensolver for an ab initio molecular dynamics package.

Integration of Shift-and-Invert Parallel Spectral Transformation (SIPs) eigensolver (as implemented in the SLEPc library) into an ab initio molecular dynamics package, SIESTA, is described. The effectiveness of the code is demonstrated on applications to polyethylene chains, boron nitride sheets, and bulk water clusters. For problems with the same number of orbitals, the performance of the SLEPc eigensolver depends on the sparsity of the matrices involved, favoring reduced dimensional systems such as polyethylene or boron nitride sheets in comparison to bulk systems like water clusters.

Mechanical architecture and folding of E. coli type 1 pilus domains.

Uropathogenic Escherichia coli attach to tissues using pili type 1. Each pilus is composed by thousands of coiled FimA domains followed by the domains of the tip fibrillum, FimF-FimG-FimH. The domains are linked by non-covalent β-strands that must resist mechanical forces during attachment.

Continuous melting through a hexatic phase in confined bilayer water.

Liquid water is not only of obvious importance but also extremely intriguing, displaying many anomalies that still challenge our understanding of such an a priori simple system. The same is true when looking at nanoconfined water: The liquid between constituents in a cell is confined to such dimensions, and there is already evidence that such water can behave very differently from its bulk counterpart. A striking finding has been reported from computer simulations for two-dimensionally confined water: The liquid displays continuous or discontinuous melting depending on its density.

Enhanced Configurational Entropy in High-Density Nanoconfined Bilayer Ice.

A novel kind of crystal order in high-density nanoconfined bilayer ice is proposed from molecular dynamics and density-functional theory simulations. A first-order transition is observed between a low-temperature proton-ordered solid and a high-temperature proton-disordered solid. The latter is shown to possess crystalline order for the oxygen positions, arranged on a close-packed triangular lattice with AA stacking.

Structural and configurational properties of nanoconfined monolayer ice from first principles.

Understanding the structural tendencies of nanoconfined water is of great interest for nanoscience and biology, where nano/micro-sized objects may be separated by very few layers of water. Here we investigate the properties of ice confined to a quasi-2D monolayer by a featureless, chemically neutral potential, in order to characterize its intrinsic behaviour. We use density-functional theory simulations with a non-local van der Waals density functional.

Performance analysis of electronic structure codes on HPC systems: a case study of SIESTA.

We report on scaling and timing tests of the SIESTA electronic structure code for ab initio molecular dynamics simulations using density-functional theory. The tests are performed on six large-scale supercomputers belonging to the PRACE Tier-0 network with four different architectures: Cray XE6, IBM BlueGene/Q, BullX, and IBM iDataPlex. We employ a systematic strategy for simultaneously testing weak and strong scaling, and propose a measure which is independent of the range of number of cores on which the tests are performed to quantify strong scaling efficiency as a function of simulation size.

A first-principles study of As doping at a disordered Si-SiO2 interface.

Understanding the interaction between dopants and semiconductor-oxide interfaces is an increasingly important concern in the drive to further miniaturize modern transistors. To this end, using a combination of first-principles density-functional theory and a continuous random network Monte Carlo method, we investigate electrically active arsenic donors at the interface between silicon and its oxide. Using a realistic model of the disordered interface, we find that a small percentage (on the order of ∼10%) of the atomic sites in the first few monolayers on the silicon side of the interface are energetically favourable for segregation, and that this is controlled by the local bonding and local strain of the defect centre.

Room temperature compressibility and diffusivity of liquid water from first principles.

The isothermal compressibility of water is essential to understand its anomalous properties. We compute it by ab initio molecular dynamics simulations of 200 molecules at five densities, using two different van der Waals density functionals. While both functionals predict compressibilities within ~30% of experiment, only one of them accurately reproduces, within the uncertainty of the simulation, the density dependence of the self-diffusion coefficient in the anomalous region.

Optimal finite-range atomic basis sets for liquid water and ice.

Finite-range numerical atomic orbitals are the basis functions of choice for several first principles methods, due to their flexibility and scalability. Generating and testing such basis sets, however, remains a significant challenge for the end user. We discuss these issues and present a new scheme for generating improved polarization orbitals of finite range.