Strongly Bound Excitons in Metal-Organic Framework MOF-5: A Many-Body Perturbation Theory Study

Mixed quantum-classical modeling of exciton-phonon scattering in solids: Application to optical linewidths of monolayer MoS2

We present a mixed quantum-classical framework for the microscopic and non-Markovian modeling of exciton-phonon scattering in solid-state materials and apply it to calculate the optical linewidths of monolayer MoS2. Within this framework, we combine reciprocal-space mixed quantum-classical dynamics with models for the quasiparticle band structure as well as the electron-hole and carrier-phonon interactions, parametrized against ab initio calculations, although noting that a direct interfacing with ab initio calculations is straightforward in principle. We introduce various parameters for truncating the Brillouin zone to select regions of interest. Variations of these parameters allow us to determine linewidths in the limit of asymptotic material sizes. The obtained asymptotic linewidths are found to agree favorably with experimental measurements across a range of temperatures. As such, our framework establishes itself as a promising route toward unraveling the non-Markovian and microscopic principles governing the nonadiabatic dynamics of solids.

Informative Training Data for Efficient Property Prediction in Metal-Organic Frameworks by Active Learning

In recent data-driven approaches to material discovery, scenarios where target quantities are expensive to compute and measure are often overlooked. In such cases, it becomes imperative to construct a training set that includes the most diverse, representative, and informative samples. Here, a novel regression tree-based active learning algorithm is employed for such a purpose. It is applied to predict the band gap and adsorption properties of metal-organic frameworks (MOFs), a novel class of materials that results from the virtually infinite combinations of their building units. Simpler and low dimensional descriptors, such as those based on stoichiometric and geometric properties, are used to compute the feature space for this model owing to their ability to better represent MOFs in the low data regime. The partitions given by a regression tree constructed on the labeled part of the data set are used to select new samples to be added to the training set, thereby limiting its size while maximizing the prediction quality. Tests on the QMOF, hMOF, and dMOF data sets reveal that our method constructs small training data sets to learn regression models that predict the target properties more efficiently than existing active learning approaches, and with lower variance. Specifically, our active learning approach is highly beneficial when labels are unevenly distributed in the descriptor space and when the label distribution is imbalanced, which is often the case for real world data. The regions defined by the tree help in revealing patterns in the data, thereby offering a unique tool to efficiently analyze complex structure-property relationships in materials and accelerate materials discovery.

Impact of Ligand Substitution and Metal Node Exchange in the Electronic Properties of Scandium Terephthalate Frameworks

The search for sustainable alternatives to established materials is a sensitive topic in materials science. Due to their unique structural and physical characteristics, the composition of metal-organic frameworks (MOFs) can be tuned by the exchange of metal nodes and the functionalization of organic ligands, giving rise to a large configurational space. Considering the case of scandium terephthalate MOFs and adopting an automatized computational framework based on density-functional theory, we explore the impact of metal substitution with the earth-abundant isoelectronic elements Al and Y, and ligand functionalization of varying electronegativity. We find that structural properties are strongly impacted by metal ion substitution and only moderately by ligand functionalization. In contrast, the energetic stability, the charge density distribution, and the electronic properties, including the size of the band gap, are primarily affected by the termination of the linker molecules. Functional groups such as OH and NH2 lead to particularly stable structures thanks to the formation of hydrogen bonds and affect the electronic structure of the MOFs by introducing midgap states.

Assessing the Role of the Kohn-Sham Density in the Calculation of the Low-Lying Bethe-Salpeter Excitation Energies

We adopt the GW many-body perturbation theory in conjunction with the Bethe-Salpeter equation (BSE) to compute 57 excitation energies of a set of 37 molecules. By using the PBEh global hybrid functional and a self-consistent scheme on the eigenvalues in GW, we show a strong dependence of the BSE energy on the starting Kohn-Sham (KS) density functional. This arises from both the quasiparticle energies and the spatial localization of the frozen KS orbitals employed to compute the BSE. In order to address the arbitrariness in the mean field choice, we adopt an orbital-tuning scheme where the amount of Fock exchange, α, is tuned to impose that the KS HOMO matches the GW quasiparticle eigenvalue, thus fulfilling the ionization potential theorem in DFT. The performance of the proposed scheme yields excellent results and it is similar to M06-2X and PBEh with α = 75%, consistent with tuned values of α ranging between 60% and 80%.

Connecting the dots for fundamental understanding of structure-photophysics-property relationships of COFs, MOFs, and perovskites using a Multiparticle Holstein Formalism

Photoactive organic and hybrid organic-inorganic materials such as conjugated polymers, covalent organic frameworks (COFs), metal-organic frameworks (MOFs), and layered perovskites, display intriguing photophysical signatures upon interaction with light. Elucidating structure-photophysics-property relationships across a broad range of functional materials is nontrivial and requires our fundamental understanding of the intricate interplay among excitons (electron-hole pair), polarons (charges), bipolarons, phonons (vibrations), inter-layer stacking interactions, and different forms of structural and conformational defects. In parallel with electronic structure modeling and data-driven science that are actively pursued to successfully accelerate materials discovery, an accurate, computationally inexpensive, and physically-motivated theoretical model, which consistently makes quantitative connections with conceptually complicated experimental observations, is equally important. Within this context, the first part of this perspective highlights a unified theoretical framework in which the electronic coupling as well as the local coupling between the electronic and nuclear degrees of freedom can be efficiently described for a broad range of quasiparticles with similarly structured Holstein-style vibronic Hamiltonians. The second part of this perspective discusses excitonic and polaronic photophysical signatures in polymers, COFs, MOFs, and perovskites, and attempts to bridge the gap between different research fields using a common theoretical construct - the Multiparticle Holstein Formalism. We envision that the synergistic integration of state-of-the-art computational approaches with the Multiparticle Holstein Formalism will help identify and establish new, transformative design strategies that will guide the synthesis and characterization of next-generation energy materials optimized for a broad range of optoelectronic, spintronic, and photonic applications.

Recommendation of Orbitals for G0W0 Calculations on Molecules and Crystals

The many-body GW approximation, especially the G₀W₀ method, has been widely used for condensed matter and molecules to calculate quasiparticle energies for ionization, electron attachment, and band gaps. Because G₀W₀ calculations are well-known to have a strong dependence on the orbitals, the goal of the present work is to provide guidance on the choice of density functional used to generate orbitals and to recommend a choice that gives the most broadly accurate results. We have systematically investigated the dependence of G₀W₀ calculations on the orbitals for 100 molecules and 8 crystals by considering orbitals obtained with a diverse set of Kohn-Sham (KS) and generalized KS (GKS) functionals (63 functionals plus Hartree-Fock). The percentage of Hartree-Fock exchange employed in density functionals has been found to have strong influence on the predicted molecular ionization energy and crystal fundamental band gaps (with optimum values between 40 and 56%), but to have less effect on predicting molecular electron affinities. The low cost of the Gaussian implementation, even with hybrid functionals in periodic calculations, the better performance of global hybrids as compared to range-separated hybrids of either than screened exchange or long-range-corrected type, and the relatively low cost of global-hybrid-functional periodic calculations using Gaussians means that one can employ global-hybrid functionals at a very reasonable cost and obtain more accurate band gaps of semiconductors than are obtained by the methods currently widely employed, namely local gradient approximations. We single out three global-hybrid functionals that give especially good results for both molecules (100 in the test set) and crystals (8 in the test set, for all of which our benchmark data are the proper band gap rather than an optical band gap uncorrected for exciton effects).