High-throughput computational stacking reveals emergent properties in natural van der Waals bilayers.

Stacking of two-dimensional (2D) materials has emerged as a facile strategy for realising exotic quantum states of matter and engineering electronic properties. Yet, developments beyond the proof-of-principle level are impeded by the vast size of the configuration space defined by layer combinations and stacking orders. Here we employ a density functional theory (DFT) workflow to calculate interlayer binding energies of 8451 homobilayers created by stacking 1052 different monolayers in various configurations. Analysis of the stacking orders in 247 experimentally known van der Waals crystals is used to validate the workflow and determine the criteria for realisable bilayers. For the 2586 most stable bilayer systems, we calculate a range of electronic, magnetic, and vibrational properties, and explore general trends and anomalies. We identify an abundance of bistable bilayers with stacking order-dependent magnetic or electrical polarisation states making them candidates for slidetronics applications.

Oxygen Vacancies Nucleate Charged Domain Walls in Ferroelectrics.

We study the influence of oxygen vacancies on the formation of charged 180° domain walls in ferroelectric BaTiO_{3} using first principles calculations. We show that it is favorable for vacancies to assemble in crystallographic planes, and that such clustering is accompanied by the formation of a charged domain wall. The domain wall has negative bound charge, which compensates the nominal positive charge of the vacancies and leads to a vanishing density of free charge at the wall. This is in contrast to the positively charged domain walls, which are nearly completely compensated by free charge from the bulk. The results thus explain the experimentally observed difference in electronic conductivity of the two types of domain walls, as well as the generic prevalence of charged domain walls in ferroelectrics. Moreover, the explicit demonstration of vacancy driven domain wall formation implies that specific charged domain wall configurations may be realized by bottom-up design for use in domain wall based information processing.