Search

The electronic structure of semiconducting 2D materials such as monolayer transition metal dichalcogenides (TMDs) are known to be tunable via environment and external fields, and van der Waals (vdW) heterostructures consisting of stacks of distinct types of 2D materials offer the possibility to further tune and optimize the electronic properties of 2D materials. In this work, we use density functional theory (DFT) calculations to calculate the structure and electronic properties of a vdW heterostructure of Janus monolayer WSSe with monolayer ZnO, both of which possess out of plane dipole moments. The effects of alignment, biaxial and uniaxial strain, orientation, and electric field on dipole moments and band edge energies of this heterostructure are calculated and examined. We find that the out of plane dipole moment of the ZnO monolayer is highly sensitive to strain, leading to the broad tunability of the heterostructure band edge energies over a range of experimentally-relevant strains. The use of strain-tunable 2D materials to control band offsets and alignment is a general strategy applicable to other vdW heterostructures, one that may be advantageous in the context of clean energy applications, including photocatalytic applications, and beyond.

Recent experiments on the antiferromagnetic intercalated transition metal dichalcogenide Fe1/3NbS2 have demonstrated reversible resistivity switching by application of orthogonal current pulses below its magnetic ordering temperature, making Fe1/3NbS2 promising for spintronics applications. Here, we perform density functional theory calculations with Hubbard U corrections of the magnetic order, electronic structure, and transport properties of crystalline Fe1/3NbS2, clarifying the origin of the different resistance states. The two experimentally proposed antiferromagnetic ground states, corresponding to in-plane stripe and zigzag ordering, are computed to be nearly degenerate. In-plane cross sections of the calculated Fermi surfaces are anisotropic for both magnetic orderings, with the degree of anisotropy sensitive to the Hubbard U value. The in-plane resistance, computed within the Kubo linear response formalism using a constant relaxation time approximation, is also anisotropic, supporting a hypothesis that the current-induced resistance changes are due to a repopulating of antiferromagnetic domains. Our calculations indicate that the transport anisotropy of Fe1/3NbS2 in the zigzag phase is reduced relative to stripe, consistent with the relative magnitudes of resistivity changes in experiment. Finally, our calculations reveal the likely directionality of the current-domain response, specifically, which domains are energetically stabilized for a given current direction.

Ab initio many-body perturbation theory within the GW approximation is a Green's function formalism widely used in the calculation of quasiparticle excitation energies of solids. In what has become an increasingly standard approach, Kohn–Sham eigenenergies, generated from a DFT calculation with a strategically-chosen exchange–correlation functional “starting point”, are used to construct G and W, and then perturbatively corrected by the resultant GW self-energy. In practice, there are several ways to construct the GW self-energy, and these can lead to variations in predicted quasiparticle energies. For example, for ZnO and TiO2, the GW fundamental gaps reported in the literature can vary by more than 1 eV depending on the GW code used. In this work, we calculate and analyze GW quasiparticle (QP) energies of these and other systems with three different GW codes: BERKELEYGW, ABINIT and YAMBO. Through a systematic analysis of the GW implementation of these three codes, we identify the primary origin of major discrepancies between codes reported in prior literature to be the different implementations the Coulomb divergence in the Fock exchange term and the frequency integration scheme of the GW self-energy. We then eliminate these discrepancies by using common numerical methods and algorithms, demonstrating that the same quasiparticle energies for a given material can be obtained with different codes, within numerical differences ascribable to the technical details of the underling implementations. This work will be important for users and developers in assessing the precision of future GW applications and methods.

Understanding the physical properties and controlling the generation of intrinsic and extrinsic defects is central to the technological adoption of 2D materials in devices. Here we identify a charged carbon-hydrogen complex at a chalcogen site (CHX) as a common, charged impurity in synthetically grown transition metal dichalcogenides (TMDs). This conclusion is drawn by comparing high resolution scanning probe microscopy measurements of nominally undoped and intentionally carbon doped TMD samples. While CH impurity densities in undoped CVD-grown WS2 and MOCVD-grown WSe2 can range anywhere from parts per million to parts per thousand, CH densities in the percentage levels were selectively generated by a post-synthetic methane plasma treatment. Our study indicates that methane plasma treatment is a selective and clean method for the controlled introduction of a charged carbon-hydrogen complex at a surface chalcogen site, a defect that is commonly present in synthetic TMDs.

Organic molecular crystals are expected to feature appreciable electron-phonon interactions that influence their electronic properties at zero and finite temperature. In this work, we report first-principles calculations and an analysis of the electron-phonon self-energy in naphthalene crystals. We compute the zero-point renormalization and temperature dependence of the fundamental band gap, and the resulting scattering lifetimes of electronic states near the valence- and conduction-band edges employing density functional theory. Further, our calculated phonon renormalization of the GW-corrected quasiparticle band structure predicts a fundamental band gap of 5 eV for naphthalene at room temperature, in good agreement with experiments. From our calculated phonon-induced electron lifetimes, we obtain the temperature-dependent mobilities of electrons and holes in good agreement with experimental measurements at room temperature. Finally, we show that an approximate energy self-consistent computational scheme for the electron-phonon self-energy leads to the prediction of strong satellite bands in the electronic band structure. We find that a single calculation of the self-energy can reproduce the self-consistent results of the band gap renormalization and electrical mobilities for naphthalene, provided that the on-the-mass-shell approximation is used, i.e., if the self-energy is evaluated at the bare eigenvalues.

Ferroelectric materials have technological applications in information storage and electronic devices. The ferroelectric polar phase can be controlled with external fields, chemical substitution and size-effects in bulk and ultrathin film form, providing a platform for future technologies and for exploratory research. In this work, we integrate spin-polarized density functional theory (DFT) calculations, crystal structure databases, symmetry tools, workflow software, and a custom analysis toolkit to build a library of known, previously-proposed, and newly-proposed ferroelectric materials. With our automated workflow, we screen over 67,000 candidate materials from the Materials Project database to generate a dataset of 255 ferroelectric candidates, and propose 126 new ferroelectric materials. We benchmark our results against experimental data and previous first-principles results. The data provided includes atomic structures, output files, and DFT values of band gaps, energies, and the spontaneous polarization for each ferroelectric candidate. We contribute our workflow and analysis code to the open-source python packages atomate and pymatgen so others can conduct analogous symmetry driven searches for ferroelectrics and related phenomena.

A detailed understanding of charged defects in two-dimensional semiconductors is needed for the development of ultrathin electronic devices. Here, we study negatively charged acceptor impurities in monolayer WS2 using a combination of scanning tunneling spectroscopy and large-scale atomistic electronic structure calculations. We observe several localized defect states of hydrogenic wave function character in the vicinity of the valence band edge. Some of these defect states are bound, while others are resonant. The resonant states result from the multivalley valence band structure of WS2, whereby localized states originating from the secondary valence band maximum at Γ hybridize with continuum states from the primary valence band maximum at K/K′. Resonant states have important consequences for electron transport as they can trap mobile carriers for several tens of picoseconds.

The importance of metal oxide photoanodes in solar fuels technology has garnered concerted efforts in photoanode discovery in recent decades, which complement parallel efforts in development of analytical techniques and optimization strategies using standard photoanodes such as TiO2, Fe2O3 and BiVO4. Theoretical guidance of highthroughput experiments has been particularly effective in dramatically increasing the portfolio of metal oxide photoanodes, motivating a new era of photoanode development where the characterization and optimization techniques developed on traditional materials are applied to nascent photoanodes that exhibit visible light photoresponse. The compendium of metal oxide photoanodes presented in the present work can also serve as the basis for further technique development, with a primary goal to establish workflows for discovery of materials that perform better against the critical criteria of operational stability, visible light photoresponse, and photovoltage suitable for tandem absorber architectures.

A series of structurally diverse alcoholamine- and alkoxyalkylamine-functionalized variants of the metal-organic framework Mg2 (dobpdc) are shown to adsorb CO2 selectively via cooperative chain-forming mechanisms. Solid-state NMR spectra and optimized structures obtained from van der Waals-corrected density functional theory calculations indicate that the adsorption profiles can be attributed to the formation of carbamic acid or ammonium carbamate chains that are stabilized by hydrogen bonding interactions within the framework pores. These findings significantly expand the scope of chemical functionalities that can be utilized to design cooperative CO2 adsorbents, providing further means of optimizing these powerful materials for energy-efficient CO2 separations.

Halide perovskites constitute a chemically diverse class of crystals with great promise as photovoltaic absorber materials, featuring band gaps between about 1 and 3.5 eV depending on composition. Their diversity calls for a general computational approach to predicting their band gaps. However, such an approach is still lacking. Here, we use density functional theory (DFT) and ab initio many-body perturbation theory within the GW approximation to compute the quasiparticle or fundamental band gap of a set of ten representative halide perovskites: CH3NH3PbI3 (MAPbI3), MAPbBr3, CsSnBr3, (MA)2BiTlBr6, Cs2TlAgBr6, Cs2TlAgCl6, Cs2BiAgBr6, Cs2InAgCl6, Cs2SnBr6, and Cs2Au2I6. Comparing with recent measurements, we find that a standard generalized gradient exchange-correlation functional can significantly underestimate the experimental band gaps of these perovskites, particularly in cases with strong spin-orbit coupling (SOC) and highly dispersive band edges, to a degree that varies with composition. We show that these nonsystematic errors are inherited by one-shot G0W0 and eigenvalue self-consistent GW0 calculations, demonstrating that semilocal DFT starting points are insufficient for MAPbI3, MAPbBr3, CsSnBr3, (MA)2BiTlBr6, Cs2TlAgBr6, and Cs2TlAgCl6. On the other hand, we find that DFT with hybrid functionals leads to an improved starting point and GW0 results in better agreement with experiment for these perovskites. Our results suggest that GW0 with hybrid functional-based starting points are promising for predicting band gaps of systems with large SOC and dispersive bands in this technologically important class of semiconducting crystals.

Using first-principles calculations we examine the band structures of ferromagnetic hexagonal manganites YXO3 (X=V, Cr, Mn, Fe, and Co) in the nonpolar nonsymmorphic P63/mmc space group. For YVO3 and YCrO3 we find a band inversion near the Fermi energy that generates a nodal ring in the kz=0 mirror plane. We perform a more detailed analysis for these compounds and predict the existence of the topological "drumhead" surface states. Finally, we briefly discuss the low-symmetry polar phases (space group P63cm) of these systems, and show they can undergo a P63/mmc→P63cm transition by condensation of soft K3 and Γ2- phonons. Based on our findings, stabilizing these compounds in the hexagonal phase could offer a promising platform for studying the interplay of topology and multiferroicity, and the coexistence of real-space and reciprocal-space topological protection in the same phase.

With the aim of systematically comparing two popular approaches to density functional theory–all-electron calculations with local basis sets, and periodic calculations employing plane wave basis sets and norm-conserving pseudopotentials–we have computed complete-basis binding energies across the S22 set of intermolecular interactions, a dataset consisting of noncovalent interactions of small- and medium-sized molecules containing first- and second-row atoms, using the Troullier-Martins norm-conserving pseudopotentials with SPW92, a local spin-density approximation; and PBE, a generalised gradient approximation. We have found that it is challenging to reach the basis set limit with these periodic calculations; for the methods and systems examined, a minimum vacuum distance of 30 Å between a system and its nearest images is necessary–unless some form of dipole correction is employed–as is a kinetic energy cutoff of at least 80 Ry. The trends in convergence with respect to vacuum size and kinetic energy cutoff are largely independent of the level of density functional approximation employed. A sense of the impact of each hyperparameter on basis set error provides a foundation for ensuring quality calculations in future studies and allows us to quantify the basis set errors incurred in existing studies on similar systems.

While long-theorized, the direct observation of multiple highly dispersive C60 valence bands has eluded researchers for more than two decades due to a variety of intrinsic and extrinsic factors. Here we report a realization of multiple highly dispersive (330-520 meV) valence bands in pure thin film C60 on a novel substrate - the three-dimensional topological insulator Bi2Se3 - through the use of angle-resolved photoemission spectroscopy (ARPES) and first-principles calculations. The effects of this novel substrate reducing C60 rotational disorder are discussed. Our results provide important considerations for past and future band structure studies as well as the increasingly popular C60 electronic device applications, especially those making use of heterostructures.

We demonstrate the guiding principles behind simple two dimensional self-assembly of MOF nanoparticles (NPs) and oleic acid capped iron oxide (Fe3 O4 ) NCs into a uniform two-dimensional bi-layered superstructure. This self-assembly process can be controlled by the energy of ligand-ligand interactions between surface ligands on Fe3 O4 NCs and Zr6 O4 (OH)4 (fumarate)6 MOF NPs. Scanning transmission electron microscopy (TEM)/energy-dispersive X-ray spectroscopy and TEM tomography confirm the hierarchical co-assembly of Fe3 O4 NCs with MOF NPs as ligand energies are manipulated to promote facile diffusion of the smaller NCs. First-principles calculations and event-driven molecular dynamics simulations indicate that the observed patterns are dictated by combination of ligand-surface and ligand-ligand interactions. This study opens a new avenue for design and self-assembly of MOFs and NCs into high surface area assemblies, mimicking the structure of supported catalyst architectures, and provides a thorough fundamental understanding of the self-assembly process, which could be a guide for designing functional materials with desired structure.

Transient absorption (TA) spectroscopy is uniquely suited for understanding kinetic processes initiated by light over vast ranges of time. In combination with white light probes, the recorded differential absorption spectra can contain spectroscopic signatures characteristic of specific charge carrier population densities. However, disentangling the often-complex and convoluted spectra is made challenging without robust analysis methods relating the underlying physical mechanisms to the spectral components. In this work, we address the origin of the transient spectra of a model system of emerging solar energy harvesting materials using a monoclinic BiVO4 thin film. Using ground-state optical properties of the semiconductor, we find the main derivative-like spectral response to be related to shifting and broadening of oscillators, rather than specific carrier-related transitions. However, by using the Drude optical model of free carriers, we also identify the transient response related to free-hole density. Importantly, sample heating from the optical pumping, which begins at â10 ps and plateaus by â200 ps, dominates the overall spectral response at longer times. On the basis of a physical model of the spectral response, a kinetic model is developed that describes the pump power dependence of the free-hole density, as well as the temporal evolution of the spectral changes associated with shifting and broadening of oscillators. First-principles density functional theory calculations are used to rationalize experimental measurements. This comprehensive approach to analyzing and modeling the TA spectra offers a generalizable basis for understanding the complex pump-probe data, reveals thermal heating artifacts that are frequently erroneously assigned to long-lived photocarriers, and offers a path to eliminating ambiguity in analysis of photocarrier dynamics in solid-state systems.

Improving the efficiency of solar-powered oxygen evolution is both critical for development of solar fuels technologies and challenging due to the broad set of properties required of a solar fuels photoanode. Bismuth vanadate, in particular the monoclinic clinobisvanite phase, has received substantial attention and has exhibited the highest radiative efficiency among metal oxides with a band gap in the visible range. Efforts to further improve its photoelectrochemical performance have included alloying one or more metals onto the Bi and/or V sites, with progress on this frontier stymied by the difficulty in computational modelling of substitutional alloys and the high dimensionality of co-alloying composition spaces. Since substitutional alloying concurrently changes multiple materials properties, understanding the underlying cause for performance improvements is also challenging, motivating our application of combinatorial materials science techniques to map photoelectrochemical performance of 948 unique bismuth vanadate alloy compositions comprising 0 to 8% alloys of P, Ca, Mo, Eu, Gd, and W along with a variety of compositions from each pairwise combination of these elements. Upon identification of substantial improvements in the (Mo,Gd) co-alloying space, structural mapping was performed to reveal a remarkable correlation between performance enhancement and a lowered monoclinic distortion. First-principles density functional theory calculations indicate that the improvements are due to a lowered hole effective mass and hole polaron formation energy, and collectively, our results identify the monoclinic distortion as a critical parameter in the optimization and understanding of bismuth vanadate-based photoanodes.

We combine transport, magnetization, and torque magnetometry measurements to investigate the electronic structure of ZrTe5, a system that is thought to be near a topological phase transition. At fields beyond the quantum limit, we observe a magnetization reversal from paramagnetic to diamagnetic response, which is characteristic of a Dirac semimetal. However, on increasing temperature across a corresponding transport anomaly, all signatures of this Dirac-like nature are completely suppressed, providing the first thermodynamic evidence of a possible topological phase transition in this compound. ZrTe5 may thus provide a rare, experimentally accessible example in which such phase transitions can be studied directly.

We present the crystal structure, electronic structure, and transport properties of the material YbMnSb2, a candidate system for the investigation of Dirac physics in the presence of magnetic order. Our measurements reveal that this system is a low-carrier-density semimetal with a two-dimensional Fermi surface arising from a Dirac dispersion, consistent with the predictions of density-functional-theory calculations of the antiferromagnetic system. The low temperature resistivity is very large, suggesting that scattering in this system is highly efficient at dissipating momentum despite its Dirac-like nature.

We propose the use of three-dimensional Dirac materials as targets for direct detection of sub-MeV dark matter. Dirac materials are characterized by a linear dispersion for low-energy electronic excitations, with a small band gap of O(meV) if lattice symmetries are broken. Dark matter at the keV scale carrying kinetic energy as small as a few meV can scatter and excite an electron across the gap. Alternatively, bosonic dark matter as light as a few meV can be absorbed by the electrons in the target. We develop the formalism for dark matter scattering and absorption in Dirac materials and calculate the experimental reach of these target materials. We find that Dirac materials can play a crucial role in detecting dark matter in the keV to MeV mass range that scatters with electrons via a kinetically mixed dark photon, as the dark photon does not develop an in-medium effective mass. The same target materials provide excellent sensitivity to absorption of light bosonic dark matter in the meV to hundreds of meV mass range, superior to all other existing proposals when the dark matter is a kinetically mixed dark photon.

Solar-driven oxygen evolution is a critical technology for renewably synthesizing hydrogen- and carbon-containing fuels in solar fuel generators. New photoanode materials are needed to meet efficiency and stability requirements, motivating materials explorations for semiconductors with (i) band-gap energy in the visible spectrum and (ii) stable operation in aqueous electrolyte at the electrochemical potential needed to evolve oxygen from water. Motivated by the oxygen evolution competency of many Mn-based oxides, the existence of several Bi-containing ternary oxide photoanode materials, and the variety of known oxide materials combining these elements with Sm, we explore the Bi-Mn-Sm oxide system for new photoanodes. Through the use of a ferri/ferrocyanide redox couple in high-throughput screening, BiMn2O5 and its alloy with Sm are identified as photoanode materials with a near-ideal optical band gap of 1.8 eV. Using density functional theory-based calculations of the mullite Bi3+Mn3+Mn4+O5 phase, we identify electronic analogues to the well-known BiVO4 photoanode and demonstrate excellent Pourbaix stability above the oxygen evolution Nernstian potential from pH 4.5 to 15. Our suite of experimental and computational characterization indicates that BiMn2O5 is a complex oxide with the necessary optical and chemical properties to be an efficient, stable solar fuel photoanode.

The solar photoelectrochemical generation of hydrogen and carbon-containing fuels comprises a critical energy technology for establishing sustainable energy resources. The photoanode, which is responsible for solar-driven oxygen evolution, has persistently limited technology advancement due to the lack of materials that exhibit both the requisite electronic properties and operational stability. Efforts to extend the lifetime of solar fuel devices increasingly focus on mitigating corrosion in the highly oxidizing oxygen evolution environment, motivating our development of a photoanode discovery pipeline that combines electronic structure calculations, Pourbaix stability screening, and high-throughput experiments. By applying the pipeline to ternary metal oxides containing manganese, we identify a promising class of corrosion-resistant materials and discover five oxygen evolution photoanodes, including the first demonstration of photoelectrocatalysis with Mn-based ternary oxides and the introduction of alkaline earth manganates as promising photoanodes for establishing a durable solar fuels technology.

Half-metallic ferromagnetism (HMFM) occurs rarely in materials and yet offers great potential for spintronic devices. Recent experiments suggest a class of compounds with the ThCr2Si2 (122) structure (isostructural and containing elements common with Fe pnictide-based superconductors) can exhibit HMFM. Here, we use ab initio density functional theory calculations to understand the onset of half-metallicity in this family of materials and explain the appearance of ferromagnetism at a quantum critical point. We also predict new candidate materials with HMFM and potentially high Curie temperatures through A-site alloying.

Antiferroelectric PbZrO3 is being considered for a wide range of applications where the competition between centrosymmetric and noncentrosymmetric phases is important to the response. Here, we focus on the epitaxial growth of PbZrO3 thin films and understanding the chemistry-structure coupling in Pb1+δZrO3 (δ = 0, 0.1, 0.2). High-quality, single-phase Pb1+δZrO3 films are synthesized via pulsed-laser deposition. Although no significant lattice parameter change is observed in X-ray studies, electrical characterization reveals that while the PbZrO3 and Pb1.1ZrO3 heterostructures remain intrinsically antiferroelectric, the Pb1.2ZrO3 heterostructures exhibit a hysteresis loop indicative of ferroelectric response. Further X-ray scattering studies reveal strong quarter-order diffraction peaks in PbZrO3 and Pb1.1ZrO3 heterostructures indicative of antiferroelectricity, while no such peaks are observed for Pb1.2ZrO3 heterostructures. Density functional theory calculations suggest the large cation nonstoichiometry is accommodated by incorporation of antisite PbZr defects, which drive the Pb1.2ZrO3 heterostructures to a ferroelectric phase with R3c symmetry. In the end, stabilization of metastable phases in materials via chemical nonstoichiometry and defect engineering enables a novel route to manipulate the energy of the ground state of materials and the corresponding material properties.

We combine high-energy resolution fluorescence detection (HERFD) X-ray absorption spectroscopy (XAS) measurements with first-principles density functional theory (DFT) calculations to provide a molecular-scale understanding of local structure, and its role in defining optoelectronic properties, in CH3NH3Pb(I1-xBrx)3 perovskites. The spectra probe a ligand field splitting in the unoccupied d states of the material, which lie well above the conduction band minimum and display high sensitivity to halide identity, Pb-halide bond length, and Pb-halide octahedral tilting, especially for apical halide sites. The spectra are also sensitive to the organic cation. We find that the halides in these mixed compositions are randomly distributed, rather than having preferred octahedral sites, and that thermal tilting motions dominate over any preferred structural distortions as a function of halide composition. These findings demonstrate the utility of the combined HERFD XAS and DFT approach for determining structural details in these materials and connecting them to optoelectronic properties observed by other characterization methods.

The alignment of the frontier orbital energies of an adsorbed molecule with the substrate Fermi level at metal-organic interfaces is a fundamental observable of significant practical importance in nanoscience and beyond. Typical density functional theory calculations, especially those using local and semi-local functionals, often underestimate level alignment leading to inaccurate electronic structure and charge transport properties. In this work, we develop a new fully self-consistent predictive scheme to accurately compute level alignment at certain classes of complex heterogeneous molecule-metal interfaces based on optimally tuned range-separated hybrid functionals. Starting from a highly accurate description of the gas-phase electronic structure, our method by construction captures important nonlocal surface polarization effects via tuning of the long-range screened exchange in a range-separated hybrid in a non-empirical and system-specific manner. We implement this functional in a plane-wave code and apply it to several physisorbed and chemisorbed molecule-metal interface systems. Our results are in quantitative agreement with experiments, the both the level alignment and work function changes. Our approach constitutes a new practical scheme for accurate and efficient calculations of the electronic structure of molecule-metal interfaces.

A parameter-free version of the recently developed driven Liouville-von Neumann equation [T. Zelovich et al., J. Chem. Theory Comput. 10(8), 2927-2941 (2014)] for electronic transport calculations in molecular junctions is presented. The single driving rate, appearing as a fitting parameter in the original methodology, is replaced by a set of state-dependent broadening factors applied to the different single-particle lead levels. These broadening factors are extracted explicitly from the self-energy of the corresponding electronic reservoir and are fully transferable to any junction incorporating the same lead model. The performance of the method is demonstrated via tight-binding and extended Hückel calculations of simple junction models. Our analytic considerations and numerical results indicate that the developed methodology constitutes a rigorous framework for the design of "black-box" algorithms to simulate electron dynamics in open quantum systems out of equilibrium.

We summarize the MOLGW code that implements density-functional theory and many-body perturbation theory in a Gaussian basis set. The code is dedicated to the calculation of the many-body self-energy within the GW approximation and the solution of the Bethe–Salpeter equation. These two types of calculations allow the user to evaluate physical quantities that can be compared to spectroscopic experiments. Quasiparticle energies, obtained through the calculation of the GW self-energy, can be compared to photoemission or transport experiments, and neutral excitation energies and oscillator strengths, obtained via solution of the Bethe–Salpeter equation, are measurable by optical absorption. The implementation choices outlined here have aimed at the accuracy and robustness of calculated quantities with respect to measurements. Furthermore, the algorithms implemented in MOLGW allow users to consider molecules or clusters containing up to 100 atoms with rather accurate basis sets, and to choose whether or not to apply the resolution-of-the-identity approximation. Finally, we demonstrate the parallelization efficacy of the MOLGW code over several hundreds of processors. Program title: MOLGW Catalogue identifier: AFAW_v1_0 Program summary URL: http://cpc.cs.qub.ac.uk/summaries/AFAW_v1_0.html Program obtainable from: CPC Program Library, Queen's University, Belfast, N. Ireland Licensing provisions: GNU GPL v3 No. of lines in distributed program, including test data, etc.: 167871 No. of bytes in distributed program, including test data, etc.: 1309269 Distribution format: tar.gz Programming language: Fortran 2003 with a few C subroutines, Python scripts. Classification: 7.3, 16.6, 16.10. External routines: libint [2], libxc [3], SCALAPACK [4] (optional) Nature of problem: Prediction of the electronic structure of atoms, molecules, clusters with a particular interest in their spectroscopic features, such as quasiparticle energies and optical spectra. Solution method: Using the GW approximation to many-body perturbation theory, MOLGW calculates total-energies, quasiparticle energies, and optical excitations. Additional comments: Python3 is required to run the test suite provided. Running time: From 30 s to a few hours References: [1] http://www.gnu.org/copyleft/gpl.txt[2] https://github.com/evaleev/libint[3] http://www.tddft.org/programs/octopus/wiki/index.php/Libxc[4] http://www.netlib.org/scalapack

We provide direct evidence for the existence of isolated, one-dimensional charge density waves at mirror twin boundaries (MTBs) of single-layer semiconducting MoSe 2. Such MTBs have been previously observed by transmission electron microscopy and have been predicted to be metallic in MoSe 2 and MoS 2. Our low-temperature scanning tunnelling microscopy/spectroscopy measurements revealed a substantial bandgap of 100 meV opening at the Fermi energy in the otherwise metallic one-dimensional structures. We found a periodic modulation in the density of states along the MTB, with a wavelength of approximately three lattice constants. In addition to mapping the energy-dependent density of states, we determined the atomic structure and bonding of the MTB through simultaneous high-resolution non-contact atomic force microscopy. Density functional theory calculations based on the observed structure reproduced both the gap opening and the spatially resolved density of states.

Photovoltaic devices based on hybrid perovskite materials have exceeded 22% efficiency due to high charge-carrier mobilities and lifetimes. Properties such as photocurrent generation and open-circuit voltage are influenced by the microscopic structure and orientation of the perovskite crystals, but are difficult to quantify on the intra-grain length scale and are often treated as homogeneous within the active layer. Here, we map the local short-circuit photocurrent, open-circuit photovoltage, and dark drift current in state-of-the-art methylammonium lead iodide solar cells using photoconductive atomic force microscopy. We find, within individual grains, spatially correlated heterogeneity in short-circuit current and open-circuit voltage up to 0.6 V. These variations are related to different crystal facets and have a direct impact on the macroscopic power conversion efficiency. We attribute this heterogeneity to a facet-dependent density of trap states. These results imply that controlling crystal grain and facet orientation will enable a systematic optimization of polycrystalline and single-crystal devices for photovoltaic and lighting applications.

We apply a recently developed Wannier-Koopmans method (WKM), based on density functional theory (DFT), to calculate the electronic energy level alignment at an interface between a molecule and metal substrate. We consider two systems: benzenediamine on Au (111), and a bipyridine-Au molecular junction. The WKM calculated level alignment agrees well with the experimental measurements where available, as well as previous GW and DFT + Σ results. Our results suggest that the WKM is a general approach that can be used to correct DFT eigenvalue errors, not only in bulk semiconductors and isolated molecules, but also in hybrid interfaces.

We present accurate force fields developed from density functional theory (DFT) calculations with periodic boundary conditions for use in molecular simulations involving M2(dobdc) (M-MOF-74; dobdc4- = 2,5-dioxidobenzenedicarboxylate; M = Mg, Mn, Fe, Co, Ni, Zn) and frameworks of similar topology. In these systems, conventional force fields fail to accurately model gas adsorption due to the strongly binding open-metal sites. The DFT-derived force fields predict the adsorption of CO2, H2O, and CH4 inside these frameworks much more accurately than other common force fields. We show that these force fields can also be used for M2(dobpdc) (dobpdc4- = 4,4′-dioxidobiphenyl-3,3′-dicarboxylate), an extended version of MOF-74, and thus are a promising alternative to common force fields for studying materials similar to MOF-74 for carbon capture applications. Furthermore, it is anticipated that the approach can be applied to other metal-organic framework topologies to obtain force fields for different systems. We have used this force field to study the effect of contaminants such as H2O and N2 upon these materials' performance for the separation of CO2 from the emissions of natural gas reservoirs and coal-fired power plants. Specifically, mixture adsorption isotherms calculated with these DFT-derived force fields showed a significant reduction in the uptake of many gas components in the presence of even trace amounts of H2O vapor. The extent to which the various gases are affected by the concentration of H2O in the reservoir is quantitatively different for the different frameworks and is related to their heats of adsorption. Additionally, significant increases in CO2 selectivities over CH4 and N2 are observed as the temperature of the systems is lowered.

Acene molecular crystals are of current interest in organic optoelectronics, both as active materials and for exploring and understanding new phenomena. Phonon scattering can be an important facilitator and dissipation mechanism in charge separation and carrier transport processes. Here, we carry out density functional theory (DFT) calculations of the structure and the full phonon dispersion of crystalline naphthalene, a well-characterized acene crystal for which detailed neutron-diffraction measurements, as well as infrared and Raman spectroscopy, are available. We evaluate the performance, relative to experiments, of DFT within the local density approximation (LDA); the generalized gradient approximation of Perdew, Burke, and Ernzerhof (PBE); and a recent van der Waals-corrected nonlocal correlation (vdW-DF-cx) functional. We find that the vdW-DF-cx functional accurately predicts lattice parameters of naphthalene within 1%. Intermolecular and intramolecular phonon frequencies across the Brillouin zone are reproduced within 7.8% and 1%, respectively. As expected, LDA (PBE) underestimates (overestimates) the lattice parameters and overestimates (underestimates) phonon frequencies, demonstrating their shortcomings for predictive calculations of weakly bound materials. If the unit cell is fixed to the experimental lattice parameters, PBE is shown to lead to improved phonon frequencies. Our study provides a detailed understanding of the phonon spectrum of naphthalene, and highlights the importance of including van der Waals dispersion interactions in predictive calculations of lattice parameters and phonon frequencies of molecular crystals and related organic materials.

Using density functional theory (DFT) with a van der Waals density functional, we calculate the adsorption energetics and geometry of benzenediamine (BDA) molecules on Au(111) surfaces. Our results demonstrate that the reported self-assembled linear chain structure of BDA, stabilized via hydrogen bonds between amine groups, is energetically favored over previously studied monomeric phases. Moreover, using a model, which includes nonlocal polarization effects from the substrate and the neighboring molecules and incorporates many-body perturbation theory calculations within the GW approximation, we obtain approximate self-energy corrections to the DFT highest occupied molecular orbital (HOMO) energy associated with BDA adsorbate phases. We find that, independent of coverage, the HOMO energy of the linear chain phase is lower relative to the Fermi energy than that of the monomer phase, and in good agreement with values measured with ultraviolet photoelectron spectroscopy and x-ray photoelectron spectroscopy.

Solar fuel generators entail a high degree of materials integration, and efficient photoelectrocatalysis of the constituent reactions hinges upon the establishment of highly functional interfaces. The recent application of high throughput experimentation to interface discovery for solar fuels photoanodes has revealed several surprising and promising mixed-metal oxide coatings for BiVO4. Using sputter deposition of composition and thickness gradients on a uniform BiVO4 film, we systematically explore photoanodic performance as a function of the composition and loading of Fe-Ce oxide coatings. This combinatorial materials integration study not only enhances the performance of this new class of materials but also identifies CeO2 as a critical ingredient that merits detailed study. A heteroepitaxial CeO2(001)/BiVO4(010) interface is identified in which Bi and V remain fully coordinated to O such that no surface states are formed. Ab initio calculations of the integrated materials and inspection of the electronic structure reveals mechanisms by which CeO2 facilitates charge transport while mitigating deleterious recombination. The results support the observations that addition of Ce to BiVO4 coatings greatly enhances photoelectrocatalytic activity, providing an important strategy for developing a scalable solar fuels technology.

Solar fuel generators entail a high degree of materials integration, and efficient photoelectrocatalysis of the constituent reactions hinges upon the establishment of highly functional interfaces. The recent application of high throughput experimentation to interface discovery for solar fuels photoanodes has revealed several surprising and promising mixed-metal oxide coatings for BiVO4. Using sputter deposition of composition and thickness gradients on a uniform BiVO4 film, we systematically explore photoanodic performance as a function of the composition and loading of Fe-Ce oxide coatings. This combinatorial materials integration study not only enhances the performance of this new class of materials but also identifies CeO2 as a critical ingredient that merits detailed study. A heteroepitaxial CeO2(001)/BiVO4(010) interface is identified in which Bi and V remain fully coordinated to O such that no surface states are formed. Ab initio calculations of the integrated materials and inspection of the electronic structure reveals mechanisms by which CeO2 facilitates charge transport while mitigating deleterious recombination. The results support the observations that addition of Ce to BiVO4 coatings greatly enhances photoelectrocatalytic activity, providing an important strategy for developing a scalable solar fuels technology.

Complex oxide β-Mn2V2O7 is identified as exhibiting near-optimal band energetics for solar fuel applications among known metal oxides. Experiments, corroborated by theory, indicate a bandgap near 1.8 eV. The calculations predict that β-Mn2V2O7 has well-aligned band edge energies for the hydrogen evolution reaction and oxygen evolution reaction. Photoelectrochemical measurements indicate appreciable photocurrent, corroborating the predictions.