Two-dimensional (2D) materials have emerged as promising candidates for next-generation electronic and optoelectronic applications. Yet, only a few dozen 2D materials have been successfully synthesized or exfoliated. Here, we search for 2D materials that can be easily exfoliated from their parent compounds. Starting from 108,423 unique, experimentally known 3D compounds, we identify a subset of 5,619 compounds that appear layered according to robust geometric and bonding criteria. High-throughput calculations using van der Waals density functional theory, validated against experimental structural data and calculated random phase approximation binding energies, further allowed the identification of 1,825 compounds that are either easily or potentially exfoliable. In particular, the subset of 1,036 easily exfoliable cases provides novel structural prototypes and simple ternary compounds as well as a large portfolio of materials to search from for optimal properties. For a subset of 258 compounds, we explore vibrational, electronic, magnetic and topological properties, identifying 56 ferromagnetic and antiferromagnetic systems, including half-metals and half-semiconductors.

Wannier90 is an open-source computer program for calculating maximally-localised Wannier functions (MLWFs) from a set of Bloch states. It is interfaced to many widely used electronicstructure codes thanks to its independence from the basis sets representing these Bloch states. In the past few years the development of Wannier90 has transitioned to a community-driven model; this has resulted in a number of new developments that have been recently released in Wannier90 v3.0. In this article we describe these new functionalities, that include the implementation of new features for wannierisation and disentanglement (symmetry-adapted Wannier functions, selectivelylocalised Wannier functions, selected columns of the density matrix) and the ability to calculate new properties (shift currents and Berry-curvature dipole, and a new interface to many-body perturbation theory); performance improvements, including parallelisation of the core code; enhancements in functionality (support for spinor-valued Wannier functions, more accurate methods to interpolate quantities in the Brillouin zone); improved usability (improved plotting routines, integration with arXiv:1907.09788v1 [cond-mat.mtrl-sci]

yambo is an open source project aimed at studying excited state properties of condensed matter systems from first principles using many-body methods. As input, yambo requires ground state electronic structure data as computed by density functional theory codes such as Quantum ESPRESSO and Abinit. yambo's capabilities include the calculation of linear response quantities (both independentparticle and including electron-hole interactions), quasi-particle corrections based on the GW formalism, optical absorption, and other spectroscopic quantities. Here we describe recent developments ranging from the inclusion of important but oft-neglected physical effects such as electron-phonon interactions to the implementation of a real-time propagation scheme for simulating linear and nonlinear optical properties. Improvements to numerical algorithms and the user interface are outlined. Particular emphasis is given to the new and efficient parallel structure that makes it possible to exploit modern high performance computing architectures. Finally, we demonstrate the possibility to automate workflows by interfacing with the yambopy and AiiDA software tools. CONTENTS

Despite the enormous success and popularity of density-functional theory, systematic verification and validation studies are still limited in number and scope. Here, we propose a protocol to test publicly available pseudopotential libraries, based on several independent criteria including verification against all-electron equations of state and plane-wave convergence tests for phonon frequencies, band structure, cohesive energy and pressure. Adopting these criteria we obtain curated pseudopotential libraries (named SSSP or standard solid-state pseudopotential libraries), that we target for highthroughput materials screening ("SSSP efficiency") and high-precision materials modelling ("SSSP precision"). This latter scores highest among open-source pseudopotential libraries available in the ∆-factor test of equations of states of elemental solids.

Fundamental research and technological applications of topological insulators are hindered by the rarity of materials exhibiting a robust topologically nontrivial phase, especially in two dimensions. Here, by means of extensive first-principles calculations, we propose a novel quantum spin Hall insulator with a sizable band gap of ∼0.5  eV that is a monolayer of jacutingaite, a naturally occurring layered mineral first discovered in 2008 in Brazil and recently synthesized. This system realizes the paradigmatic Kane-Mele model for quantum spin Hall insulators in a potentially exfoliable two-dimensional monolayer, with helical edge states that are robust and that can be manipulated exploiting a unique strong interplay between spin-orbit coupling, crystal-symmetry breaking, and dielectric response.

Quantum spin Hall insulators are a class of two-dimensional materials with a finite electronic band gap in the bulk and gapless helical edge states. In the presence of time-reversal symmetry, Z2 topological order distinguishes the topological phase from the ordinary insulating one. Some of the phenomena that can be hosted in these materials, from one-dimensional low-dissipation electronic transport to spin filtering, could be very promising for many technological applications in the fields of electronics, spintronics and topological quantum computing. Nevertheless, the rarity of twodimensional materials that can exhibit non-trivial Z2 topological order at room temperature hinders development. Here, we screen a comprehensive database we recently created of 1825 monolayers that can be exfoliated from experimentally known compounds, to search for novel quantum spin Hall insulators. Using density-functional and many-body perturbation theory simulations, we identify 13 monolayers that are candidates for quantum spin Hall insulators, including high-performing materials such as AsCuLi2 and jacutingaite (Pt2HgSe3). We also identify monolayer Pd2HgSe3 as a novel Kane-Mele quantum spin Hall insulator, and compare it with jacutingaite. Such a handful of promising materials are mechanically stable and exhibit Z2 topological order, either unpertubed or driven by a small amount of strain. Such screening highlights a relative abundance of Z2 topological order of around 1%, and provides an optimal set of candidates for experimental efforts. arXiv:1908.08334v1 [cond-mat.mtrl-sci]

Maximally-localised Wannier functions (MLWFs) are routinely used to compute from first-principles advanced materials properties that require very dense Brillouin zone integration and to build accurate tight-binding models for scale-bridging simulations. At the same time, high-throughput (HT) computational materials design is an emergent field that promises to accelerate reliable and cost-effective design and optimisation of new materials with target properties. The use of MLWFs in HT workflows has been hampered by the fact that generating MLWFs automatically and robustly without any user intervention and for arbitrary materials is, in general, very challenging. We address this problem directly by proposing a procedure for automatically generating MLWFs for HT frameworks. Our approach is based on the selected columns of the density matrix method and we present the details of its implementation in an AiiDA workflow. We apply our approach to a dataset of 200 bulk crystalline materials that span a wide structural and chemical space. We assess the quality of our MLWFs in terms of the accuracy of the band-structure interpolation that they provide as compared to the band-structure obtained via full first-principles calculations. Finally, we provide a downloadable virtual machine that can be used to reproduce the results of this paper, including all first-principles and atomistic simulations as well as the computational workflows.

We report high-resolution angle resolved photoemission measurements on single crystals of Pt2HgSe3 grown by high-pressure synthesis. Our data reveal a gapped Dirac nodal line whose (001)-projection separates the surface Brillouin zone in topological and trivial areas. In the nontrivial k-space range we find surface states with multiple saddle-points in the dispersion resulting in two van Hove singularities in the surface density of states. Based on density functional theory calculations, we identify these surface states as signatures of a topological crystalline state which coexists with a weak topological phase.