Semilocal Density Functional Theory Research Articles

Molecular orbital projectors in non-empirical jmDFT recover exact conditions in transition-metal chemistry.

Low-cost, non-empirical corrections to semi-local density functional theory are essential for accurately modeling transition-metal chemistry. Here, we demonstrate the judiciously modified density functional theory (jmDFT) approach with non-empirical U and J parameters obtained directly from frontier orbital energetics on a series of transition-metal complexes. We curate a set of nine representative Ti(III) and V(IV) d1 transition-metal complexes and evaluate their flat-plane errors along the fractional spin and charge lines. We demonstrate that while jmDFT improves upon both DFT+U and semi-local DFT with the standard atomic orbital projectors (AOPs), it does so inefficiently. We rationalize these inefficiencies by quantifying hybridization in the relevant frontier orbitals. To overcome these limitations, we introduce a procedure for computing a molecular orbital projector (MOP) basis for use with jmDFT. We demonstrate this single set of d1 MOPs to be suitable for nearly eliminating all energetic delocalization and static correlation errors. In all cases, MOP jmDFT outperforms AOP jmDFT, and it eliminates most flat-plane errors at non-empirical values. Unlike DFT+U or hybrid functionals, jmDFT nearly eliminates energetic delocalization and static correlation errors within a non-empirical framework.

Data-driven approach to parameterize SCAN+U for an accurate description of 3d transition metal oxide thermochemistry

Semilocal density-functional theory (DFT) methods exhibit significant errors for the phase diagrams of transition-metal oxides that are caused by an incorrect description of molecular oxygen and the large self-interaction error in materials with strongly localized electronic orbitals. Empirical and semiempirical corrections based on the $\mathrm{DFT}+U$ method can reduce these errors, but the parameterization and validation of the correction terms remains an on-going challenge. We develop a systematic methodology to determine the parameters and to statistically assess the results by considering interlinked thermochemical data across a set of transition metal compounds. We consider three interconnected levels of correction terms: (1) a constant oxygen binding correction, (2) Hubbard-$U$ correction, and (3) $\mathrm{DFT}/\mathrm{DFT}+U$ compatibility correction. The parameterization is expressed as a unified optimization problem. We demonstrate this approach for $3d$ transition metal oxides, considering a target set of binary and ternary oxides. With a total of 37 measured formation enthalpies taken from the literature, the dataset is augmented by the reaction energies of 1710 unique reactions that were derived from the formation energies by systematic enumeration. To ensure a balanced dataset across the available data, the reactions were grouped by their similarity using clustering and suitably weighted. The parameterization is validated using leave-one-out cross validation (CV), a standard technique for the validation of statistical models. We apply the methodology to the strongly constrained and appropriately normed (SCAN) density functional. Based on the CV score, the error of binary (ternary) oxide formation energies is reduced by 40% (75%) to 0.10 (0.03) eV/atom. A simplified correction scheme that does not involve $\mathrm{SCAN}/\mathrm{SCAN}+U$ compatibility terms still achieves an error reduction of 30% (25%). The method and tools demonstrated here can be applied to other classes of materials or to parameterize the corrections to optimize $\mathrm{DFT}+U$ performance for other target physical properties.

Exploring cesium-tellurium phase space via high-throughput calculations beyond semi-local density-functional theory.

Boosted by the relentless increase in available computational resources, high-throughput calculations based on first-principles methods have become a powerful tool to screen a huge range of materials. The backbone of these studies is well-structured and reproducible workflows efficiently returning the desired properties given chemical compositions and atomic arrangements as sole input. Herein, we present a new workflow designed to compute the stability and the electronic properties of crystalline materials from density-functional theory using the strongly constrained and appropriately normed approximation (SCAN) for the exchange-correlation potential. We show the performance of the developed tool exploring the binary Cs-Te phase space that hosts cesium telluride, a semiconducting material widely used as a photocathode in particle accelerators. Starting from a pool of structures retrieved from open computational material databases, we analyze formation energies as a function of the relative Cs content and for a few selected crystals, we investigate the band structures and density of states unraveling interconnections among the structure, stoichiometry, stability, and electronic properties. Our study contributes to the ongoing research on alkali-based photocathodes and demonstrates that high-throughput calculations based on state-of-the-art first-principles methods can complement experiments in the search for optimal materials for next-generation electron sources.

Protective Coating Interfaces for Perovskite Solar Cell Materials: A First-Principles Study.

The protection of halide perovskites is important for the performance and stability of emergent perovskite-based optoelectronic technologies. In this work, we investigate the potential inorganic protective coating materials ZnO, SrZrO3, and ZrO2 for the CsPbI3 perovskite. The optimal interface registries are identified with Bayesian optimization. We then use semilocal density functional theory (DFT) to determine the atomic structure at the interfaces of each coating material with the clean CsI-terminated surface and three reconstructed surface models with added PbI2 and CsI complexes. For the final structures, we explore the level alignment at the interface with hybrid DFT calculations. Our analysis of the level alignment at the coating–substrate interfaces reveals no detrimental mid-gap states but rather substrate-dependent valence and conduction band offsets. While ZnO and SrZrO3 act as insulators on CsPbI3, ZrO2 might be suitable as an electron transport layer with the right interface engineering.

Improving the transferability of density functional theory predictions through correlation analysis: Structural and energetic properties of NiX alloys (X=C, Si, Ge, and Sn)

This work reports on the performance of density functional theory (DFT) for a series of single and binary systems, aiming for a quantitative description of $\mathrm{Ni}X$ ($X=\mathrm{C}$, Si, Ge, and Sn) alloys. Both semilocal GGA and a meta-GGA density functional, with and without dispersion corrections, are tested. We found in our study that no single functional simultaneously provides an accurate quantitative description of the investigated structural and energetic properties. However, the spread in computed DFT data could be rationalized in terms of the distribution of reduced density gradients and differences in the evolution of the exchange enhancement factors for different functionals. We demonstrate how to construct a regression model based on data from several density functionals that increases the predictivity of semilocal DFT. We foresee that the use of regression models (or extensions of it) can be valuable in the development of more accurate density functionals that in the future could provide a quantitative accuracy for complex multicomponent systems.

Eliminating Delocalization Error to Improve Heterogeneous Catalysis Predictions with Molecular DFT + U

Approximate semilocal density functional theory (DFT) is known to underestimate surface formation energies yet paradoxically overbind adsorbates on catalytic transition-metal oxide surfaces due to delocalization error. The low-cost DFT + U approach only improves surface formation energies for early transition-metal oxides or adsorption energies for late transition-metal oxides. In this work, we demonstrate that this inefficacy arises due to the conventional usage of metal-centered atomic orbitals as projectors within DFT + U. We analyze electron density rearrangement during surface formation and O atom adsorption on rutile transition-metal oxides to highlight that a standard DFT + U correction fails to tune properties when the corresponding density rearrangement is highly delocalized across both metal and oxygen sites. To improve both surface properties simultaneously while retaining the simplicity of a single-site DFT + U correction, we systematically construct multi-atom-centered molecular-orbital-like projectors for DFT + U. We demonstrate this molecular DFT + U approach for tuning adsorption energies and surface formation energies of minimal two-dimensional models of representative early (i.e., TiO2) and late (i.e., PtO2) transition-metal oxides. Molecular DFT + U simultaneously corrects adsorption energies and surface formation energies of multilayer models of rutile TiO2(110) and PtO2(110) to resolve the paradoxical description of surface stability and surface reactivity of semilocal DFT.

Elastic stability of Ga2O3 : Addressing the β to α phase transition from first principles

The elastic and structural properties of $\ensuremath{\beta}\text{\ensuremath{-}}{\mathrm{Ga}}_{2}{\mathrm{O}}_{3}$ and $\ensuremath{\alpha}\text{\ensuremath{-}}{\mathrm{Ga}}_{2}{\mathrm{O}}_{3}$ are investigated from first principles. The full elastic tensors and elastic moduli of both phases at $0\phantom{\rule{0.16em}{0ex}}\mathrm{K}$ are computed in the framework of semilocal density-functional theory. We determine mechanical instabilities of $\ensuremath{\beta}\text{\ensuremath{-}}{\mathrm{Ga}}_{2}{\mathrm{O}}_{3}$ by evaluating the full stiffness tensor under load for a range of hydrostatic pressure values. While a phase transition from the $\ensuremath{\beta}$ to $\ensuremath{\alpha}$ phase is found to be energetically favored at $2.6\phantom{\rule{0.16em}{0ex}}\mathrm{G}\mathrm{Pa}$, we show that the $\ensuremath{\beta}$ phase is mechanically unstable only for much higher pressures $(>30\phantom{\rule{0.16em}{0ex}}\mathrm{G}\mathrm{Pa})$, which agrees well with experimental results. Our employed approach is based on the Born stability criterion, is independent of crystal symmetry, and thus can be readily applied to different materials.

Beyond-dipole van der Waals contributions within the many-body dispersion framework

By introducing a suitable range-separation of the Coulomb coupling in analogy to Ambrosetti et al (2014 J. Chem. Phys. 140 18A508), here we extend the many-body dispersion approach to include beyond-dipole van der Waals (vdW) interactions at a full many-body level, in combination with semi-local density functional theory. A reciprocal-space implementation is further introduced in order to efficiently treat periodic systems. Consistent reliability is found from molecular dimers to large supramolecular complexes and two-dimensional systems. The large weight of both many-body effects and multipolar terms illustrates how a correct description of vdW forces in large-scale systems requires full account of both contributions, beyond standard pairwise dipolar approaches.

Enabling Large-Scale Condensed-Phase Hybrid Density Functional Theory-Based Ab Initio Molecular Dynamics II: Extensions to the Isobaric-Isoenthalpic and Isobaric-Isothermal Ensembles.

In the previous paper of this series [Ko, H.-Y. et al. J. Chem. Theory Comput. 2020, 16, 3757-3785], we presented a theoretical and algorithmic framework based on a localized representation of the occupied space that exploits the inherent sparsity in the real-space evaluation of the exact exchange (EXX) interaction in finite-gap systems. This was accompanied by a detailed description of exx, a massively parallel hybrid message-passing interface MPI/OpenMP implementation of this approach in Quantum ESPRESSO (QE) that enables linear scaling hybrid density functional theory (DFT)-based ab initio molecular dynamics (AIMD) in the microcanonical/canonical (NVE/NVT) ensembles of condensed-phase systems containing 500-1000 atoms (in fixed orthorhombic cells) with a wall time cost comparable to semi-local DFT. In this work, we extend the current capabilities of exx to enable hybrid DFT-based AIMD simulations of large-scale condensed-phase systems with general and fluctuating cells in the isobaric-isoenthalpic/isobaric-isothermal (NpH/NpT) ensembles. The theoretical extensions to this approach include an analytical derivation of the EXX contribution to the stress tensor for systems in general simulation cells with a computational complexity that scales linearly with system size. The corresponding algorithmic extensions to exx include optimized routines that (i) handle both static and fluctuating simulation cells with non-orthogonal lattice symmetries, (ii) solve Poisson's equation in general/non-orthogonal cells via an automated selection of the auxiliary grid directions in the Natan-Kronik representation of the discrete Laplacian operator, and (iii) evaluate the EXX contribution to the stress tensor. Using this approach, we perform a case study on a variety of condensed-phase systems (including liquid water, a benzene molecular crystal polymorph, and semi-conducting crystalline silicon) and demonstrate that the EXX contributions to the energy and stress tensor simultaneously converge with an appropriate choice of exx parameters. This is followed by a critical assessment of the computational performance of the extended exx module across several different high-performance computing architectures via case studies on (i) the computational complexity due to lattice symmetry during NpT simulations of three different ice polymorphs (i.e., ice Ih, II, and III) and (ii) the strong/weak parallel scaling during large-scale NpT simulations of liquid water. We demonstrate that the robust and highly scalable implementation of this approach in the extended exx module is capable of evaluating the EXX contribution to the stress tensor with negligible cost (<1%) as well as all other EXX-related quantities needed during NpT simulations of liquid water (with a very tight 150 Ry planewave cutoff) in ≈5.2 s ((H2O)128) and ≈6.8 s ((H2O)256) per AIMD step. As such, the extended exx module presented in this work brings us another step closer to routinely performing hybrid DFT-based AIMD simulations of sufficient duration for large-scale condensed-phase systems across a wide range of thermodynamic conditions.

Assessing the G0W0Γ0(1) Approach: Beyond G0W0 with Hedin's Full Second-Order Self-Energy Contribution.

We present and benchmark a self-energy approach for quasiparticle energy calculations that goes beyond Hedin's GW approximation by adding the full second-order self-energy (FSOS-W) contribution. The FSOS-W diagram involves two screened Coulomb interaction (W) lines, and adding the FSOS-W to the GW self-energy can be interpreted as first-order vertex correction to GW (GWΓ(1)). Our FSOS-W implementation is based on the resolution-of-identity technique and exhibits better than O(N5) scaling with system size for small- to medium-sized molecules. We then present one-shot GWΓ(1) (G0W0Γ0(1)) benchmarks for the GW100 test set and a set of 24 acceptor molecules. For semilocal or hybrid density functional theory starting points, G0W0Γ0(1) systematically outperforms G0W0 for the first vertical ionization potentials and electron affinities of both test sets. Finally, we demonstrate that a static FSOS-W self-energy significantly underestimates the quasiparticle energies.

Machine-learning semilocal density functional theory for many-body lattice models at zero and finite temperature

We introduce a machine-learning density-functional-theory formalism for the spinless Hubbard model in one dimension at both zero and finite temperature. In the zero-temperature case this establishes a one-to-one relation between the site occupation and the total energy, which is then minimised at the ground-state occupation. In contrast, at finite temperature the same relation is defined between the Helmholtz free energy and the equilibrium site occupation. Most importantly, both functionals are semi-local, so that they are independent from the size of the system under investigation and can be constructed over exact data for small systems. These 'exact' functionals are numerically defined by neural networks. We also define additional neural networks for finite-temperature thermodynamical quantities, such as the entropy and heat capacity. These can be either a functional of the ground-state site occupation or of the finite-temperature equilibrium site occupation. In the first case their equilibrium value does not correspond to an extremal point of the functional, while it does in the second case. Our work gives us access to finite-temperature properties of many-body systems in the thermodynamic limit.

Efficient Band Structure Calculation of Two-Dimensional Materials from Semilocal Density Functionals.

The experimental and theoretical realization of two-dimensional (2D) materials is of utmost importance in semiconducting applications. Computational modeling of these systems with satisfactory accuracy and computational efficiency is only feasible with semilocal density functional theory methods. In the search for the most useful method in predicting the band gap of 2D materials, we assess the accuracy of recently developed semilocal exchange–correlation (XC) energy functionals and potentials. Though the explicit forms of exchange–correlation (XC) potentials are very effective against XC energy functionals for the band gap of bulk solids, their performance needs to be investigated for 2D materials. In particular, the LMBJ [J. Chem. Theory Comput.2020, 16, 265432097004] and GLLB-SC [Phys. Rev. B82, 2010, 115106] potentials are considered for their dominance in bulk band gap calculation. The performance of recently developed meta generalized gradient approximations, like TASK [Phys. Rev. Res.1, 2019, 033082] and MGGAC [Phys. Rev. B. 100, 2019, 155140], is also assessed. We find that the LMBJ potential constructed for 2D materials is not as successful as its parent functional, i.e., MBJ [Phys. Rev. Lett.102, 2009, 22640119658882] in bulk solids. Due to a contribution from the derivative discontinuity, the band gaps obtained with GLLB-SC are in a certain number of cases, albeit not systematically, larger than those obtained with other methods, which leads to better agreement with the quasi-particle band gap obtained from the GW method. The band gaps obtained with TASK and MGGAC can also be quite accurate.

Revisiting Understanding of Electrochemical CO2 Reduction on Cu(111): Competing Proton-Coupled Electron Transfer Reaction Mechanisms Revealed by Embedded Correlated Wavefunction Theory.

Copper (Cu) electrodes, as the most efficacious of CO2 reduction reaction (CO2RR) electrocatalysts, serve as prototypes for determining and validating reaction mechanisms associated with electrochemical CO2 reduction to hydrocarbons. As in situ electrochemical mechanism determination by experiments is still out of reach, such mechanistic analysis typically is conducted using density functional theory (DFT). The semilocal exchange-correlation (XC) approximations most often used to model such catalysis unfortunately engender a basic error: predicting the wrong adsorption site for CO (a key CO2RR intermediate) on the most ubiquitous facet of Cu, namely, Cu(111). This longstanding inconsistency casts lingering doubt on previous DFT predictions of the attendant CO2RR kinetics. Here, we apply embedded correlated wavefunction (ECW) theory, which corrects XC functional error, to study the CO2RR on Cu(111) via both surface hydride (*H) transfer and proton-coupled electron transfer (PCET). We predict that adsorbed CO (*CO) reduces almost equally to two intermediates, namely, hydroxymethylidyne (*COH) and formyl (*CHO) at -0.9 V vs the RHE. In contrast, semilocal DFT approximations predict a strong preference for *COH. With increasing applied potential, the dominance of *COH (formed via potential-independent surface *H transfer) diminishes, switching to the competitive formation of both *CHO and *COH (both formed via potential-dependent PCET). Our results also demonstrate the importance of including explicitly modeled solvent molecules in predicting electron-transfer barriers and reveal the pitfalls of overreliance on simple surface *H transfer models of reduction reactions.

Photoelectron spectra of early 3d-transition metal dioxide molecular anions from GW calculations.

Photoelectron spectra of early 3d-transition metal dioxide anions, ScO2 -, TiO2 -, VO2 -, CrO2 -, and MnO2 -, are calculated using semilocal and hybrid density functional theory (DFT) and many-body perturbation theory within the GW approximation using one-shot perturbative and eigenvalue self-consistent formalisms. Different levels of theory are compared with each other and with available photoelectron spectra. We show that one-shot GW with a PBE0 starting point (G0W0@PBE0) consistently provides very good agreement for all experimentally measured binding energies (within 0.1 eV-0.2 eV or less). We attribute this to the success of PBE0 in mitigating self-interaction error and providing good quasiparticle wave functions, which renders a first-order perturbative GW correction effective. One-shot GW calculations with a Perdew-Burke-Ernzerhof (PBE) starting point do poorly in predicting electron removal energies by underbinding orbitals with typical errors near 1.5 eV. A higher exact exchange amount of 50% in the DFT starting point of one-shot GW does not provide very good agreement with experiment by overbinding orbitals with typical errors near 0.5 eV. While not as accurate as G0W0@PBE0, the G-only eigenvalue self-consistent GW scheme with W fixed to the PBE level provides a reasonably predictive level of theory (typical errors near 0.3 eV) to describe photoelectron spectra of these 3d-transition metal dioxide anions. Adding eigenvalue self-consistency also in W, on the other hand, worsens the agreement with experiment overall. Our findings on the performance of various GW methods are discussed in the context of our previous studies on other transition metal oxide molecular systems.

R2SCAN-3c: A "Swiss army knife" composite electronic-structure method.

The recently proposed r2SCAN meta-generalized-gradient approximation (mGGA) of Furness and co-workers is used to construct an efficient composite electronic-structure method termed r2SCAN-3c. To this end, the unaltered r2SCAN functional is combined with a tailor-made triple-ζ Gaussian atomic orbital basis set as well as with refitted D4 and geometrical counter-poise corrections for London-dispersion and basis set superposition error. The performance of the new method is evaluated for the GMTKN55 database covering large parts of chemical space with about 1500 data points, as well as additional benchmarks for non-covalent interactions, organometallic reactions, and lattice energies of organic molecules and ices, as well as for the adsorption on polar salt and non-polar coinage-metal surfaces. These comprehensive tests reveal a spectacular performance and robustness of r2SCAN-3c: It by far surpasses its predecessor B97-3c at only twice the cost and provides one of the best results of all semi-local density-functional theory (DFT)/QZ methods ever tested for the GMTKN55 database at one-tenth of the cost. Specifically, for reaction and conformational energies as well as non-covalent interactions, it outperforms prominent hybrid-DFT/QZ approaches at two to three orders of magnitude lower cost. Perhaps, the most relevant remaining issue of r2SCAN-3c is self-interaction error (SIE), owing to its mGGA nature. However, SIE is slightly reduced compared to other (m)GGAs, as is demonstrated in two examples. After all, this remarkably efficient and robust method is chosen as our new group default, replacing previous composite DFT and partially even expensive high-level methods in most standard applications for systems with up to several hundreds of atoms.

Random-Phase Approximation in Many-Body Noncovalent Systems: Methane in a Dodecahedral Water Cage.

The many-body expansion (MBE) of energies of molecular clusters or solids offers a way to detect and analyze errors of theoretical methods that could go unnoticed if only the total energy of the system was considered. In this regard, the interaction between the methane molecule and its enclosing dodecahedral water cage, CH4···(H2O)20, is a stringent test for approximate methods, including density functional theory (DFT) approximations. Hybrid and semilocal DFT approximations behave erratically for this system, with three- and four-body nonadditive terms having neither the correct sign nor magnitude. Here, we analyze to what extent these qualitative errors in different MBE contributions are conveyed to post-Kohn-Sham random-phase approximation (RPA), which uses approximate Kohn-Sham orbitals as its input. The results reveal a correlation between the quality of the DFT input states and the RPA results. Moreover, the renormalized singles energy (RSE) corrections play a crucial role in all orders of the many-body expansion. For dimers, RSE corrects the RPA underbinding for every tested Kohn-Sham model: generalized-gradient approximation (GGA), meta-GGA, (meta-)GGA hybrids, as well as the optimized effective potential at the correlated level. Remarkably, the inclusion of singles in RPA can also correct the wrong signs of three- and four-body nonadditive energies as well as mitigate the excessive higher-order contributions to the many-body expansion. The RPA errors are dominated by the contributions of compact clusters. As a workable method for large systems, we propose to replace those compact contributions with CCSD(T) energies and to sum up the remaining many-body contributions up to infinity with supermolecular or periodic RPA. As a demonstration of this approach, we show that for RPA(PBE0)+RSE it suffices to apply CCSD(T) to dimers and 30 compact, hydrogen-bonded trimers to get the methane-water cage interaction energy to within 1.6% of the reference value.

New insights into the 1D carbon chain through the RPA.

We investigated the electronic and structural properties of the infinite linear carbon chain (carbyne) using density functional theory (DFT) and the random phase approximation (RPA) to the correlation energy. The studies are performed in vacuo and for carbyne inside a carbon nano tube (CNT). In the vacuum, semi-local DFT and RPA predict bond length alternations of about 0.04 Å and 0.13 Å, respectively. The frequency of the highest optical mode at the Γ point is 1219 cm-1 and about 2000 cm-1 for DFT and the RPA. Agreement of the RPA to previous high level quantum chemistry and diffusion Monte-Carlo results is excellent. For the RPA we calculate the phonon-dispersion in the full Brillouine zone and find marked quantitative differences to DFT calculations not only at the Γ point but also throughout the entire Brillouine zone. To model carbyne inside a carbon nanotube, we considered a (10,0) CNT. Here the DFT calculations are even qualitatively sensitive to the k-points sampling. At the limes of a very dense k-points sampling, semi-local DFT predicts no bond length alternation (BLA), whereas in the RPA a sizeable BLA of 0.09 Å prevails. The reduced BLA leads to a significant red shift of the vibrational frequencies of about 350 cm-1, so that they are in good agreement with experimental estimates. Overall, the good agreement between the RPA and previously reported results from correlated wavefunction methods and experimental Raman data suggests that the RPA provides reliable results at moderate computational costs. It hence presents a useful addition to the repertoire of correlated wavefunction methods and its accuracy clearly prevails for low dimensional systems, where semi-local density functionals struggle to yield even qualitatively correct results.

Critical assessment of Co-Cu phase diagram from first-principles calculations

Using first-principles alloy theory, we perform a systematic study of the Co-Cu phase diagram. Calculations are carried out for ferromagnetic and paramagnetic ${\mathrm{Co}}_{1\ensuremath{-}x}{\mathrm{Cu}}_{x}$ solid solutions with face-centered-cubic (fcc) crystal structure. We find that the equilibrium volumes and magnetic states are crucial for a quantitative description of the thermodynamics of the Co-Cu system at temperatures up to 1400 K. In particular, the paramagnetic state of Cu-rich alloys with persisting local magnetic moments is shown to be responsible for the solubility of a small amount of Co in fcc Cu whereas the excess entropy in the ferromagnetic Co-rich region critically depends on the adopted lattice parameters. None of the common local or semilocal density functional theory approximations have the necessary accuracy for the lattice parameters when compared to the experimental data. The predicted ab initio Co-Cu phase diagram is in good agreement with the measurements and CALPHAD data, making it possible to gain a deep insight into the various contributions to the Gibbs free energy. The present study provides an atomic-level description of the thermodynamic quantities controlling the limited mutual solubility of Co and Cu and highlights the importance of high-temperature magnetism.

Termination Effects in Aluminosilicate and Aluminogermanate Imogolite Nanotubes: A Density Functional Theory Study

We investigate termination effects in aluminosilicate (AlSi) and aluminogermanate (AlGe) imogolite nanotubes (NTs) by means of semi-local and range-corrected hybrid Density Functional Theory (DFT) simulations. Following screening and identification of the smallest finite model capable of accommodating full relaxation of the NT terminations around an otherwise geometrically and electrostatically unperturbed core region, we quantify and discuss the effects of physical truncation on the structure, relative energy, electrostatics and electronic properties of differently terminated, finite-size models of the NTs. In addition to composition-dependent changes in the valence (VB) and conduction band (CB) edges and resultant band gap (BG), the DFT simulations uncover longitudinal band bending and separation in the finite AlSi and AlGe models. Depending on the given termination of the NTs, such longitudinal effects manifest in conjunction with the radial band separation typical of fully periodic AlSi and AlGe NTs. The strong composition dependence of the longitudinal and radial band bending in AlSi and AlGe NTs suggests different mechanisms for the generation, relaxation and separation of photo-generated holes in AlSi and AlGe NTs, inviting further research in the untapped potential of imogolite compositional and structural flexibility for photo-catalytic applications.

Photophysics and Electronic Structure of Lateral Graphene/MoS2 and Metal/MoS2 Junctions.

Integration of semiconducting transition metal dichalcogenides (TMDs) into functional optoelectronic circuitries requires an understanding of the charge transfer across the interface between the TMD and the contacting material. Here, we use spatially resolved photocurrent microscopy to demonstrate electronic uniformity at the epitaxial graphene/molybdenum disulfide (EG/MoS2) interface. A 10× larger photocurrent is extracted at the EG/MoS2 interface when compared to the metal (Ti/Au)/MoS2 interface. This is supported by semi-local density functional theory (DFT), which predicts the Schottky barrier at the EG/MoS2 interface to be ∼2× lower than that at Ti/MoS2. We provide a direct visualization of a 2D material Schottky barrier through combination of angle-resolved photoemission spectroscopy with spatial resolution selected to be ∼300 nm (nano-ARPES) and DFT calculations. A bending of ∼500 meV over a length scale of ∼2-3 μm in the valence band maximum of MoS2 is observed via nano-ARPES. We explicate a correlation between experimental demonstration and theoretical predictions of barriers at graphene/TMD interfaces. Spatially resolved photocurrent mapping allows for directly visualizing the uniformity of built-in electric fields at heterostructure interfaces, providing a guide for microscopic engineering of charge transport across heterointerfaces. This simple probe-based technique also speaks directly to the 2D synthesis community to elucidate electronic uniformity at domain boundaries alongside morphological uniformity over large areas.