We introduce the Computational 2D Materials Database (C2DB), which organises a variety of structural, thermodynamic, elastic, electronic, magnetic, and optical properties of around 1500 two-dimensional materials distributed over more than 30 different crystal structures. Material properties are systematically calculated by density functional theory and many-body perturbation theory (G 0 W 0 and the Bethe-Salpeter Equation for ∼250 materials) following a semi-automated workflow for maximal consistency and transparency. The C2DB is fully open and can be browsed online at c2db.fysik.dtu.dk or downloaded in its entirety. In this paper, we describe the workflow behind the database, present an overview of the properties and materials currently available, and explore trends and correlations in the data. Moreover, we identify a large number of new potentially synthesisable 2D materials with interesting properties targeting applications within spintronics, (opto-)electronics, and plasmonics. The C2DB offers a comprehensive and easily accessible overview of the rapidly expanding family of 2D materials and forms an ideal platform for computational modeling and design of new 2D materials and van der Waals heterostructures.

Magnetic order in two-dimensional (2D) materials is intimately coupled to magnetic anisotropy (MA) since the Mermin-Wagner theorem implies that rotational symmetry cannot be spontaneously broken at finite temperatures in 2D. Large MA thus comprises a key ingredient in the search for magnetic 2D materials that retains the magnetic order above room temperature. Magnetic interactions are typically modeled in terms of Heisenberg models and the temperature dependence on magnetic properties can be obtained with the Random Phase Approximation (RPA), which treats magnon interactions at the mean-field level. In the present work we show that large MA gives rise to strong magnon-magnon interactions that leads to a drastic failure of the RPA. We then demonstrate that classical Monte Carlo (MC) simulations correctly describe the critical temperatures in the large MA limit and agree with RPA when the MA becomes small. A fit of the MC results leads to a simple expression for the critical temperatures as a function of MA and exchange coupling constants, which significantly simplifies the theoretical search for new 2D magnetic materials with high critical temperatures. The expression is tested on a monolayer of CrI3, which were recently observed to exhibit ferromagnetic order below 45 K and we find excellent agreement with the experimental value. * tolsen@fysik.dtu.dk

The recent observation of ferromagnetic order in two-dimensional (2D) materials has initiated a booming interest in the subject of 2D magnetism. In contrast to bulk materials, 2D materials can only exhibit magnetic order in the presence of magnetic anisotropy. In the present work we have used the Computational 2D Materials Database (C2DB) to search for new ferromagnetic 2D materials using the spinwave gap as a simple descriptor that accounts for the role of magnetic anisotropy. In addition to known compounds we find 12 novel insulating materials that exhibit magnetic order at finite temperatures. For these we evaluate the critical temperatures from classical Monte Carlo simulations of a Heisenberg model with exchange and anisotropy parameters obtained from first principles. Starting from 150 stable ferromagnetic 2D materials we find five candidates that are predicted to have critical temperatures exceeding that of CrI3. We also study the effect of Hubbard corrections in the framework of DFT+U and find that the value of U can have a crucial influence on the prediction of magnetic properties. Our work provides new insight into 2D magnetism and identifies a new set of promising monolayers for experimental investigation.

We perform a computational screening for two-dimensional (2D) magnetic materials based on experimental bulk compounds present in the Inorganic Crystal Structure Database and Crystallography Open Database. A recently proposed geometric descriptor is used to extract materials that are exfoliable into 2D derivatives and we find 85 ferromagnetic and 61 antiferromagnetic materials for which we obtain magnetic exchange and anisotropy parameters using density functional theory. For the easy-axis ferromagnetic insulators we calculate the Curie temperature based on a fit to classical Monte Carlo simulations of anisotropic Heisenberg models. We find good agreement with the experimentally reported Curie temperatures of known 2D ferromagnets and identify 10 potentially exfoliable 2D ferromagnets that have not been reported previously. In addition, we find 18 easy-axis antiferromagnetic insulators with several compounds exhibiting very strong exchange coupling and magnetic anisotropy.

We have performed a computational screening of topological two-dimensional (2D) materials from the Computational 2D Materials Database (C2DB) employing density functional theory. A full ab initio scheme for calculating hybrid Wannier functions directly from the Kohn-Sham orbitals has been implemented and the method was used to extract Z2 indices, Chern numbers and Mirror Chern numbers of 3331 2D systems including both experimentally known and hypothetical 2D materials. We have found a total of 48 quantum spin Hall insulators, 7 quantum anomalous Hall insulators and 21 crystalline topological insulators. Roughly 75 % are predicted to be dynamically stable and one third was known prior to the screening. The most interesting of the novel topological insulators are investigated in more detail. We show that the calculated topological indices of the quantum anomalous Hall insulators are highly sensitive to the approximation used for the exchangecorrelation functional and reliable predictions of the topological properties of these materials thus require methods beyond density functional theory. We also performed GW calculations, which yield a gap of 0.65 eV for the quantum spin Hall insulator PdSe2 in the MoS2 crystal structure. This is significantly higher than any known 2D topological insulator and three times larger than the Kohn-Sham gap.