Self-interaction corrected SCAN functional for molecules and solids in the numeric atom-center orbital framework.

Semilocal density-functional approximations (DFAs), including the state-of-the-art SCAN functional, are plagued by the self-interaction error (SIE). While this error is explicitly defined only for one-electron systems, it has inspired the self-interaction correction method proposed by Perdew and Zunger (PZ-SIC), which has shown promise in mitigating the many-electron SIE. However, the PZ-SIC method is known for its significant numerical instability. In this study, we introduce a novel constraint that facilitates self-consistent localization of the SIC orbitals in the spirit of Edmiston-Ruedenberg orbitals [Rev. Mod. Phys. 35, 457 (1963)]. Our practical implementation within the all-electron numeric atom-centered orbitals code FHI-aims guarantees efficient and stable convergence of the self-consistent PZ-SIC equations for both molecules and solids. We further demonstrate that our PZ-SIC approach effectively mitigates the SIE in the meta-generalized gradient approximation SCAN functional, significantly improving the accuracy for ionization potentials, charge-transfer energies, and bandgaps for a diverse selection of molecules and solids. However, our PZ-SIC method does have its limitations. It cannot improve the already accurate SCAN results for properties such as cohesive energies, lattice constants, and bulk modulus in our test sets. This highlights the need for new-generation DFAs with more comprehensive applicability.

Evaluating the Role of Anharmonic Vibrations in Zeolite ß Materials.

The characterization of zeolitic materials is often facilitated by spectroscopic analysis of vibrations, which informs about the bonding character of the substrate and any adsorbents. Computational simulations aid the interpretation of the spectra but often ignore anharmonic effects that can affect the spectral characteristics significantly. Here, the impact of anharmonicity is demonstrated with a combination of dynamical and static simulations applied to the structures formed during the synthesis of Sn-BEA via solid-state incorporation (SSI): the initial siliceous BEA (Si-ß), aluminosilicate BEA (H-ß), dealuminated BEA (deAl-ß), and Sn-BEA (Sn-ß). Heteroatom and defect-containing BEA are shown to have strong anharmonic vibrational contributions, with atomic and elemental resolution highlighting particularly the prevalence for H atoms (H-ß, deAl-ß) as well as localization to heteroatoms at defect sites. We simulate the vibrational spectra of BEA accounting for anharmonic contributions and observe an improved agreement with experimental data compared to harmonic methods, particularly at wavenumbers below 1500 cm-1. The results demonstrate the importance of incorporating anharmonic effects in simulations of vibrational spectra, with consequences toward future characterization and application of zeolitic materials.

Anharmonicity in Thermal Insulators: An Analysis from First Principles.

The anharmonicity of atomic motion limits the thermal conductivity in crystalline solids. However, a microscopic understanding of the mechanisms active in strong thermal insulators is lacking. In this Letter, we classify 465 experimentally known materials with respect to their anharmonicity and perform fully anharmonic ab initio Green-Kubo calculations for 58 of them, finding 28 thermal insulators with κ<10 W/mK including 6 with ultralow κ≲1 W/mK. Our analysis reveals that the underlying strong anharmonic dynamics is driven by the exploration of metastable intrinsic defect geometries. This is at variance with the frequently applied perturbative approach, in which the dynamics is assumed to evolve around a single stable geometry.

A machine learning route between band mapping and band structure.

The electronic band structure and crystal structure are the two complementary identifiers of solid-state materials. Although convenient instruments and reconstruction algorithms have made large, empirical, crystal structure databases possible, extracting the quasiparticle dispersion (closely related to band structure) from photoemission band mapping data is currently limited by the available computational methods. To cope with the growing size and scale of photoemission data, here we develop a pipeline including probabilistic machine learning and the associated data processing, optimization and evaluation methods for band-structure reconstruction, leveraging theoretical calculations. The pipeline reconstructs all 14 valence bands of a semiconductor and shows excellent performance on benchmarks and other materials datasets. The reconstruction uncovers previously inaccessible momentum-space structural information on both global and local scales, while realizing a path towards integration with materials science databases. Our approach illustrates the potential of combining machine learning and domain knowledge for scalable feature extraction in multidimensional data.

Accessing the Anisotropic Nonthermal Phonon Populations in Black Phosphorus.

We combine ultrafast electron diffuse scattering experiments and first-principles calculations of the coupled electron-phonon dynamics to provide a detailed momentum-resolved picture of lattice thermalization in black phosphorus. The measurements reveal the emergence of highly anisotropic nonthermal phonon populations persisting for several picoseconds after exciting the electrons with a light pulse. Ultrafast dynamics simulations based on the time-dependent Boltzmann formalism are supplemented by calculations of the structure factor, defining an approach to reproduce the experimental signatures of nonequilibrium structural dynamics. The combination of experiments and theory enables us to identify highly anisotropic electron-phonon scattering processes as the primary driving force of the nonequilibrium lattice dynamics in black phosphorus. Our approach paves the way toward unravelling and controlling microscopic energy flows in two-dimensional materials and van der Waals heterostructures, and may be extended to other nonequilibrium phenomena involving coupled electron-phonon dynamics such as superconductivity, phase transitions, or polaron physics.

Local Atomic Arrangements and Band Structure of Boron Carbide.

Boron carbide, the simple chemical combination of boron and carbon, is one of the best-known binary ceramic materials. Despite that, a coherent description of its crystal structure and physical properties resembles one of the most challenging problems in materials science. By combining ab initio computational studies, precise crystal structure determination from diffraction experiments, and state-of-the-art high-resolution transmission electron microscopy imaging, this concerted investigation reveals hitherto unknown local structure modifications together with the known structural alterations. The mixture of different local atomic arrangements within the real crystal structure reduces the electron deficiency of the pristine structure CBC+B12 , answering the question about electron precise character of boron carbide and introducing new electronic states within the band gap, which allow a better understanding of physical properties.

Formation of Vacancies in Si- and Ge-based Clathrates: Role of Electron Localization and Symmetry Breaking.

The formation of framework vacancies in Si- and Ge-based type-I clathrates is studied using density-functional theory as a function of filling the cages with K and Ba atoms. Our analysis reveals the relevance of structural disorder, geometric relaxation, and electronic saturation as well as vibrational and configurational entropy. In the Si clathrates, we find that vacancies are unstable, but very differently, in Ge clathrates, up to three vacancies per unit cell can be stabilized. This contrasting behavior is largely driven by the different energy gain on populating the electronic vacancy states, which originates from the different degree of localization of the valence orbitals of Si and Ge. This also actuates a qualitatively different atomic relaxation of the framework.

Ab Initio Green-Kubo Approach for the Thermal Conductivity of Solids.

We herein present a first-principles formulation of the Green-Kubo method that allows the accurate assessment of the phonon thermal conductivity of solid semiconductors and insulators in equilibrium ab initio molecular dynamics calculations. Using the virial for the nuclei, we propose a unique ab initio definition of the heat flux. Accurate size and time convergence are achieved within moderate computational effort by a robust, asymptotically exact extrapolation scheme. We demonstrate the capabilities of the technique by investigating the thermal conductivity of extreme high and low heat conducting materials, namely, Si (diamond structure) and tetragonal ZrO_{2}.

Evidence for photogenerated intermediate hole polarons in ZnO.

Despite their pronounced importance for oxide-based photochemistry, optoelectronics and photovoltaics, only fairly little is known about the polaron lifetimes and binding energies. Polarons represent a crucial intermediate step populated immediately after dissociation of the excitons formed in the primary photoabsorption process. Here we present a novel approach to studying photoexcited polarons in an important photoactive oxide, ZnO, using infrared (IR) reflection-absorption spectroscopy (IRRAS) with a time resolution of 100 ms. For well-defined (10-10) oriented ZnO single-crystal substrates, we observe intense IR absorption bands at around 200 meV exhibiting a pronounced temperature dependence. On the basis of first-principles-based electronic structure calculations, we assign these features to hole polarons of intermediate coupling strength.

Fingerprints for spin-selection rules in the interaction dynamics of O2 at Al(111).

We perform mixed quantum-classical molecular dynamics simulations based on first-principles potential-energy surfaces to demonstrate that the scattering of a beam of singlet O2 molecules at Al(111) will enable an unambiguous assessment of the role of spin-selection rules for the adsorption dynamics. At thermal energies we predict a sticking probability that is substantially less than unity, with the repelled molecules exhibiting characteristic kinetic, vibrational and rotational signatures arising from the nonadiabatic spin transition.